Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userMach Number Dependence of Turbulent Magnetic Field Amplification: Solenoidal versus Compressive Flows
203
0
0.0
(
0
)
Added by
Christoph Federrath
Publication date
2011
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We study the growth rate and saturation level of the turbulent dynamo in magnetohydrodynamical simulations of turbulence, driven with solenoidal (divergence-free) or compressive (curl-free) forcing. For models with Mach numbers ranging from 0.02 to 20, we find significantly different magnetic field geometries, amplification rates, and saturation levels, decreasing strongly at the transition from subsonic to supersonic flows, due to the development of shocks. Both extreme types of turbulent forcing drive the dynamo, but solenoidal forcing is more efficient, because it produces more vorticity.
rate research
Read More
The excitation and further sustenance of large-scale magnetic fields in rotating astrophysical systems, including planets, stars and galaxies, is generally thought to involve a fluid magnetic dynamo effect driven by helical magnetohydrodynamic turbulence. While this scenario is appealing on general grounds, it however currently remains largely unconstrained, notably because a fundamental understanding of the nonlinear asymptotic behaviour of large-scale fluid magnetism in the astrophysically-relevant but treacherous regime of large magnetic Reynolds number $Rm$ is still lacking. We explore this problem using local high-resolution simulations of turbulent magnetohydrodynamics driven by an inhomogeneous helical forcing generating a sinusoidal profile of kinetic helicity, mimicking the hemispheric distribution of kinetic helicity in rotating turbulent fluid bodies. We identify a transition at large $Rm$ to an asymptotic nonlinear state, followed up to $Rmsimeq 3times 10^3$, characterized by an asymptotically small resistive dissipation of magnetic helicity, by its efficient spatial redistribution across the equator through turbulent fluxes driven by the hemispheric distribution of kinetic helicity, and by the presence in the tangled dynamical magnetic field of plasmoids typical of reconnection at large $Rm$.
The nonlinear and nonlocal coupling of vorticity and strain-rate constitutes a major hindrance in understanding the self-amplification of velocity gradients in turbulent fluid flows. Utilizing highly-resolved direct numerical simulations of isotropic turbulence in periodic domains of up to $12288^3$ grid points, and Taylor-scale Reynolds number $R_lambda$ in the range $140-1300$, we investigate this non-locality by decomposing the strain-rate tensor into local and non-local contributions obtained through Biot-Savart integration of vorticity in a sphere of radius $R$. We find that vorticity is predominantly amplified by the non-local strain coming beyond a characteristic scale size, which varies as a simple power-law of vorticity magnitude. The underlying dynamics preferentially align vorticity with the most extensive eigenvector of non-local strain. The remaining local strain aligns vorticity with the intermediate eigenvector and does not contribute significantly to amplification; instead it surprisingly attenuates intense vorticity, leading to breakdown of the observed power-law and ultimately also the scale-invariance of vorticity amplification, with important implications for prevailing intermittency theories.
With the aim of efficiently simulating three-dimensional multiphase turbulent flows with a phase-field method, we propose a new discretization scheme for the biharmonic term (the 4th-order derivative term) of the Cahn-Hilliard equation. This novel scheme can significantly reduce the computational cost while retaining the same accuracy as the original procedure. Our phase-field method is built on top of a direct numerical simulation solver, named AFiD (www.afid.eu) and open-sourced by our research group. It relies on a pencil distributed parallel strategy and a FFT-based Poisson solver. To deal with large density ratios between the two phases, a pressure split method [1] has been applied to the Poisson solver. To further reduce computational costs, we implement a multiple-resolution algorithm which decouples the discretizations for the Navier-Stokes equations and the scalar equation: while a stretched wall-resolving grid is used for the Navier-Stokes equations, for the Cahn-Hilliard equation we use a fine uniform mesh. The present method shows excellent computational performance for large-scale computation: on meshes up to 8 billion nodes and 3072 CPU cores, a multiphase flow needs only slightly less than 1.5 times the CPU time of the single-phase flow solver on the same grid. The present method is validated by comparing the results to previous studies for the cases of drop deformation in shear flow, including the convergence test with mesh refinement, and breakup of a rising buoyant bubble with density ratio up to 1000. Finally, we simulate the breakup of a big drop and the coalescence of O(10^3) drops in turbulent Rayleigh-Benard convection at a Rayleigh number of $10^8$, observing good agreement with theoretical results.
Phoresis, the drift of particles induced by scalar gradients in a flow, can result in an effective compressibility, bringing together or repelling particles from each other. Here, we ask whether this effect can affect the transport of particles in a turbulent flow. To this end, we study how the dispersion of a cloud of phoretic particles is modified when injected in the flow, together with a blob of scalar, whose effect is to transiently bring particles together, or push them away from the center of the blob. The resulting phoretic effect can be quantified by a single dimensionless number. Phenomenological considerations lead to simple predictions for the mean separation between particles, which are consistent with results of direct numerical simulations. Using the numerical results presented here, as well as those from previous studies, we discuss quantitatively the experimental consequences of this work and the possible impact of such phoretic mechanisms in natural systems.
Collisionless shocks are ubiquitous in the Universe and often associated with strong magnetic field. Here we use large-scale particle-in-cell simulations of non-relativistic perpendicular shocks in the high-Mach-number regime to study the amplification of magnetic field within shocks. The magnetic field is amplified at the shock transition due to the ion-ion two-stream Weibel instability. The normalized magnetic-field strength strongly correlates with the Alfvenic Mach number. Mock spacecraft measurements derived from PIC simulations are fully consistent with those taken in-situ at Saturns bow shock by the Cassini spacecraft.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
