11,280 research outputs found
Simulating Supersonic Turbulence in Magnetized Molecular Clouds
We present results of large-scale three-dimensional simulations of weakly
magnetized supersonic turbulence at grid resolutions up to 1024^3 cells. Our
numerical experiments are carried out with the Piecewise Parabolic Method on a
Local Stencil and assume an isothermal equation of state. The turbulence is
driven by a large-scale isotropic solenoidal force in a periodic computational
domain and fully develops in a few flow crossing times. We then evolve the flow
for a number of flow crossing times and analyze various statistical properties
of the saturated turbulent state. We show that the energy transfer rate in the
inertial range of scales is surprisingly close to a constant, indicating that
Kolmogorov's phenomenology for incompressible turbulence can be extended to
magnetized supersonic flows. We also discuss numerical dissipation effects and
convergence of different turbulence diagnostics as grid resolution refines from
256^3 to 1024^3 cells.Comment: 10 pages, 3 figures, to appear in the proceedings of the DOE/SciDAC
2009 conferenc
Full sphere hydrodynamic and dynamo benchmarks
Convection in planetary cores can generate fluid flow and magnetic fields, and a number of sophisticated codes exist to simulate the dynamic behaviour of such systems. We report on the first community activity to compare numerical results of computer codes designed to calculate fluid flow within a whole sphere. The flows are incompressible and rapidly rotating and the forcing of the flow is either due to thermal convection or due to moving boundaries. All problems defined have solutions that allow easy comparison, since they are either steady, slowly drifting or perfectly periodic. The first two benchmarks are defined based on uniform internal heating within the sphere under the Boussinesq approximation with boundary conditions that are uniform in temperature and stress-free for the flow. Benchmark 1 is purely hydrodynamic, and has a drifting solution. Benchmark 2 is a magnetohydrodynamic benchmark that can generate oscillatory, purely periodic, flows and magnetic fields. In contrast, Benchmark 3 is a hydrodynamic rotating bubble benchmark using no slip boundary conditions that has a stationary solution. Results from a variety of types of code are reported, including codes that are fully spectral (based on spherical harmonic expansions in angular coordinates and polynomial expansions in radius), mixed spectral and finite difference, finite volume, finite element and also a mixed Fourier–finite element code. There is good agreement between codes. It is found that in Benchmarks 1 and 2, the approximation of a whole sphere problem by a domain that is a spherical shell (a sphere possessing an inner core) does not represent an adequate approximation to the system, since the results differ from whole sphere results
Nanoscale imaging of equilibrium quantum Hall edge currents and of the magnetic monopole response in graphene
The recently predicted topological magnetoelectric effect and the response to
an electric charge that mimics an induced mirror magnetic monopole are
fundamental attributes of topological states of matter with broken time
reversal symmetry. Using a SQUID-on-tip, acting simultaneously as a tunable
scanning electric charge and as ultrasensitive nanoscale magnetometer, we
induce and directly image the microscopic currents generating the magnetic
monopole response in a graphene quantum Hall electron system. We find a rich
and complex nonlinear behavior governed by coexistence of topological and
nontopological equilibrium currents that is not captured by the monopole
models. Furthermore, by utilizing a tuning fork that induces nanoscale
vibrations of the SQUID-on-tip, we directly image the equilibrium currents of
individual quantum Hall edge states for the first time. We reveal that the edge
states that are commonly assumed to carry only a chiral downstream current, in
fact carry a pair of counterpropagating currents, in which the topological
downstream current in the incompressible region is always counterbalanced by
heretofore unobserved nontopological upstream current flowing in the adjacent
compressible region. The intricate patterns of the counterpropagating
equilibrium-state orbital currents provide new insights into the microscopic
origins of the topological and nontopological charge and energy flow in quantum
Hall systems
Adaptive mesh refinement with spectral accuracy for magnetohydrodynamics in two space dimensions
We examine the effect of accuracy of high-order spectral element methods,
with or without adaptive mesh refinement (AMR), in the context of a classical
configuration of magnetic reconnection in two space dimensions, the so-called
Orszag-Tang vortex made up of a magnetic X-point centered on a stagnation point
of the velocity. A recently developed spectral-element adaptive refinement
incompressible magnetohydrodynamic (MHD) code is applied to simulate this
problem. The MHD solver is explicit, and uses the Elsasser formulation on
high-order elements. It automatically takes advantage of the adaptive grid
mechanics that have been described elsewhere in the fluid context [Rosenberg,
Fournier, Fischer, Pouquet, J. Comp. Phys. 215, 59-80 (2006)]; the code allows
both statically refined and dynamically refined grids. Tests of the algorithm
using analytic solutions are described, and comparisons of the Orszag-Tang
solutions with pseudo-spectral computations are performed. We demonstrate for
moderate Reynolds numbers that the algorithms using both static and refined
grids reproduce the pseudo--spectral solutions quite well. We show that
low-order truncation--even with a comparable number of global degrees of
freedom--fails to correctly model some strong (sup--norm) quantities in this
problem, even though it satisfies adequately the weak (integrated) balance
diagnostics.Comment: 19 pages, 10 figures, 1 table. Submitted to New Journal of Physic
- …