PaperlayerView exact revision on arXiv

Source-linked AI summary

Large-Scale Dynamos Driven by Shear-Flow-Induced Jets

B. Tripathi, A. E. Fraser, P. W. Terry, E. G. Zweibel, M. J. Pueschel, R. Fan

arXiv:2608.12530v1astro-ph.SRastro-ph.HEphysics.flu-dynphysics.plasm-phphysics.space-ph

TL;DR

Classical dynamo theory leaves some parameters unjustified from first principles, motivating a mechanistic account of large-scale field generation in shear turbulence. The paper combines analytic modeling with 3D incompressible magnetohydrodynamic simulations and finds that robust 3D jets drive a quasi-periodic large-scale dynamo through the mean-vorticity effect.

  • Problem

    The paper addresses how the jet-driven dynamo differs fundamentally from classical and other dynamos.

  • Method

    The authors combine analytic dynamo theory with 3D incompressible magnetohydrodynamic simulations of externally forced shear turbulence.

  • Results

    Robust large-scale 3D jets drive the dynamo through Yoshizawa’s mean-vorticity effect, producing polarity-flipping large-scale magnetic fields.

  • Takeaways & Limitations

    The jet-driven dynamo applies to shear-driven laboratory and astrophysical plasmas, including systems resembling binary neutron star merger interfaces.

Abstract

from arXiv · show

At every scale they occupy, magnetic fields affect various phenomena, including star formation, cosmic ray transport, charged particle acceleration, space weather, transport in planetary atmospheres, and laboratory plasmas. These fields are often generated and sustained by turbulent flows in a process called the dynamo. In 1955, E. N. Parker parameterized the effects of small-scale turbulence to propose a mean-field dynamo theory. The widely used theory reproduces observed large-scale fields but suffers from difficulty in tuning parameters as they are not justified from first principles: Studies of turbulent flows show tangled magnetic fields, which are folded and fragmented into small-scale structures due to shear-flow straining. Here, considering a shear flow that is unstable and driven, we develop analytic theory and perform three-dimensional (3D), advanced computer simulations of turbulence with up to 4096 x 4096 x 8192 grid points, showing ab initio generation of quasi-periodic, large-scale magnetic fields. The generation occurs via the mean-vorticity effect---an additional mean-field dynamo process postulated in 1990. Crucial to this dynamo is the prior generation of large-scale 3D jets, robustly produced as topologically protected and exact nonlinear solutions of the magnetohydrodynamic equations. The jet-driven dynamo applies to shear-driven laboratory and astrophysical systems. These include binary neutron star mergers, where the reported dynamo likely operates on microsecond timescales to produce in milliseconds some of the strongest magnetic fields in the Universe, providing signals for multimessenger astronomy.

Methods · I. Model and computational considerations.

The study uses extensive three-dimensional incompressible magnetohydrodynamic simulations to investigate dynamo behavior across varied domain and magnetic-field conditions. The model evolves velocity and magnetic fields under constrained equations, including vector-potential and solenoidal conditions.

  • I. Model and computational considerations.: 90 million CPU hours supported more than 90 KH dynamo simulations and analysis of 0.25 petabytes of turbulence data.The data came from evolving velocity u and magnetic field b according to the 3D incompressible MHD equations.
  • I. Model and computational considerations.: The simulations evolve velocity u and magnetic field b using the 3D incompressible MHD equations.The computational campaign generated the turbulence dataset analyzed in the study.
  • I. Model and computational considerations.: The formulation includes kinematic viscosity ν, electric resistivity η, fluid pressure P, scalar function Ψ, and external acceleration f.The fluid density satisfies μ0ρ = 1, with μ0 denoting vacuum magnetic permeability.
  • I. Model and computational considerations.: Both non-periodic and periodic domains are considered, with turbulence simulated with and without an external mean magnetic field.These alternatives define the domain and imposed-field conditions explored in the computational study.

II. Background profiles and non-dimensionalization.

The study uses single or double shear-layer backgrounds in non-periodic or periodic domains, with scales normalized by U0 and a. It tests dynamo behavior under varied initial magnetic fields and random perturbations using two pseudospectral solvers.

  • Background profiles: Single shear layers define non-periodic backgrounds, while double shear layers define periodic backgrounds in specified L-sized computational volumes.The periodic setup uses z1=L/4, z2=3L/4, and L=10𝜋.
  • Non-dimensionalization: Speed and length scales are measured with U0 and a, and the Reynolds numbers are defined as Re=U0a/𝜈 and Rm=U0a/𝜂.These normalizations are used throughout the work.
  • Initial magnetic field: Dynamo saturation is largely insensitive to MA variations below the Kelvin–Helmholtz threshold and unaffected by variations in θ.The imposed field is B0=(𝐞"xcosθ+𝐞"y sinθ)/MA, with θ measuring its angle relative to the mean flow.
  • Initialization and solvers: Broadband, phase-randomized, divergence-free perturbations seed the instability, and the jet-driven dynamo emerges from general random perturbations.Equations are solved with Dedalus in non-periodic and periodic domains and GHOST in periodic domains.

III. Dynamos in non-periodic and periodic domains.

The dynamo operates in both non-periodic and periodic domains because large-scale jets, strongest near the shear layer, enable it independently of boundary type. In non-periodic simulations, perfectly conducting no-slip walls confine the plasma while allowing a prescribed shear flow.

  • Non-periodic domains: Non-periodic domains use perfectly conducting, no-slip walls that confine the plasma and impose ux=±U0 at z=± L/2.The boundaries satisfy jx=jy=bz=0 and are represented by uz=uy=Ax=Ay=Ψ=0.
  • Non-periodic domains: With an initial uniform horizontal mean field B0, the mean field is isolated from the remaining fields before applying the magnetic boundary conditions.This setup is used for non-periodic domains with the specified perfectly conducting walls.
  • Domain dependence: The dynamo does not depend on domain type because large-scale jets enable it, with both jets and shear layer far from boundaries.The jets are strongest near the shear layer, which is identified as a minimal and sufficient ingredient for exciting the dynamo.

IV. Benchmarks and major code optimizations.

The study benchmarks Dedalus and GHOST simulations using different parallelization strategies, time-stepping schemes, and spectral bases. GHOST is further optimized with an exact integrating-factor treatment of visco-resistive terms and a mean-flow-maintaining forcing function.

  • Dedalus simulations: Dedalus uses MPI parallelization in periodic x and y directions, pseudospectral τ-methods, mixed RK443/SBDF2 steppers, and 3/2 dealiasing.Non-periodic domains use Fourier–Fourier–Chebyshev bases with harmonics ranging from 128^2×512 to 1,024^3.
  • GHOST simulations: GHOST parallelizes with MPI in z and OpenMP in y, binds threads to MPI tasks, and solves nonlinearities with second-order explicit Runge–Kutta.OpenMP threads are placed next to each other within NUMA nodes, and GHOST is benchmarked against Dedalus using identical initial conditions.
  • Major code optimizations: An exact integrating-factor technique solves GHOST’s visco-resistive terms, addressing instability risks from unresolved small-scale fluctuations.The implementation produces exact solutions for the visco-resistive terms.
  • Major code optimizations: A forcing function was implemented to maintain the mean flow during GHOST simulations.The forcing supplements the new visco-resistive-term integrator as a computational optimization.

V. External mean-flow forcing.

The simulations externally force only the mean flow while allowing magnetic fields to evolve freely, and produce essentially the same results across forcing timescales except for very large 𝜏f.

  • External mean-flow forcing: External forcing acts only on the mean flow ux(kx=0, ky=0), while magnetic fields evolve freely to isolate Kelvin–Helmholtz instability-driven fluctuations.The forcing is f=fx(z)𝐞"x in Eq. (S1a).
  • External mean-flow forcing: F0 removes the initial viscous relaxation of the mean shear flow, ensuring the initial state is a true MHD equilibrium.The equilibrium condition is F0+ 𝜈∇2⟨ux(t=0)⟩x,y=0.
  • External mean-flow forcing: Varying 𝜏f from 0 to infinity produces essentially the same results except when 𝜏f is very large.The tested limits are frozen mean flow at 𝜏f=0 and decaying mean flow at 𝜏f=infinity.
  • External mean-flow forcing: For very large 𝜏f, Kelvin–Helmholtz energy extraction flattens the mean-shear gradient faster than external forcing replenishes mean-flow energy.The instability extracts energy from the mean shear flow faster than the forcing injects it.

VI. Analytic quasilinear EMF model.

The analytic quasilinear model derives mean-field evolution from interactions between externally imposed perturbations, mean flow, and mean magnetic field. It incorporates coherent-straining times and shows that pure Alfvénic states suppress α while maximizing β and ϒ.

  • Perturbation-driven mean field: An imposed perturbation with wavenumber k=(0,k_y) generates the dominant mean magnetic field and drives self-consistent zonal velocity and magnetic fluctuations.The initial perturbations are u′=(0,0,u′_z) and b′=(0,0,b′_z), while interactions generate additional components such as u′_x and b′_x.
  • Mean-field evolution: The averaged field B_x evolves through resistive diffusion and the vertical derivative of the electromotive force, (∂_t−η∂_z^2)B_x=−∂_zℰ_y.This equation defines how the mean electromotive force enters the quasilinear evolution of the averaged magnetic field.
  • Quasilinear closure: The model expresses generated zonal jet and magnetic-field fluctuations using coherent-straining times τ_ZF and τ_ZM, which can be replaced by nonlinearly modified times.This replacement extends the calculation beyond the small-amplitude limit to include nonlinear turbulent interactions.
  • Transport coefficients: For generally equal straining times τ_ZF=τ_ZM, pure Alfvénic states u∥±b make α zero while β and ϒ become maximally non-zero.The non-kinematic terms associated with τ_ZF arise from the time evolution of velocity fluctuations.

VII. Zonal jets in the KH dynamo.

The Kelvin–Helmholtz instability nonlinearly generates an x-invariant vertical motion that the mean flow stretches into energetically dominant, large-scale zonal jets. Unlike x-dependent fluctuations, this three-dimensional motion retains its vertical scale instead of fragmenting.

  • VII. Zonal jets in the KH dynamo.: Energetically dominant, large-scale jets arise when the mean flow stretches the KH-generated motion uz(kx=0,ky).The KH instability exists only at kx≠0, while vertical mean flow and vertical mean field are absent.
  • VII. Zonal jets in the KH dynamo.: The x-invariant fluctuation uz(kx=0,ky) avoids the vertical-scale fragmentation that shears uz(kx≠0,ky) into multiple fragments.Rapid distortion theory shows that the mean flow maintains the z-spatial scale of the three-dimensional x-invariant fluctuation.
  • VII. Zonal jets in the KH dynamo.: The mean flow preserves uz(kx=0,ky)’s vertical scale and produces strong jets ux(kx=0,ky).By contrast, uz(kx≠0,ky) is progressively shortened in vertical scale on the (x,z)-plane before breaking into multiple fragments.

VIII. Analytic semi-local dynamo characterization.

The semi-local analysis shows that shear-layer jets are strongest at z=0, generate a mean magnetic field reversed across the layer, and drive polarity reversals over time. It identifies Yoshizawa’s generic ϒ-effect as operating specifically through 3D jet formation.

  • Jet structure: Jets are strongest at the shear-layer interface z=0 and decrease symmetrically with increasing |z|.This spatial profile is shown in Figs. 1c and 1d.
  • Mean-field generation: The jets ux(kx=0,ky) generate the mean magnetic field through field-line stretching.The analysis connects the jet structure to the field-line stretching term in the mean-field equation.
  • Mean-field reversal: The generated mean field is reversed across z=0, with its polarity flipping over time because of a phase difference between flow components.The polarity reversal is associated with the mean flow being reversed across the shear layer.
  • ϒ-effect: The Yoshizawa’s postulated generic ϒ-effect operates specifically through the formation of 3D jets.The conclusion follows from the preceding semi-local and energy-transfer analyses and is also supported by Fig. 3.

IX. Spectral energy-flux computations.

The analysis uses cylindrical wavenumber shells aligned with the z axis because Kelvin–Helmholtz turbulence is strongly anisotropic and localized around the shear layer. Energy-flux computations show kinetic and magnetic energies cascading from large to small physical scales.

  • Spectral energy-flux computations: Because Kelvin–Helmholtz turbulence is strongly anisotropic and localized around the shear layer, the analysis uses cylindrical shells in the (kx,ky)-plane with axis along z.This replaces spherical wavenumber shells for post-processing.
  • Spectral energy-flux computations: The u-to-u energy flux Πuu(k0) measures velocity-energy transfer through wavenumber k0 using the full system velocity.The velocity field is decomposed using Fourier amplitudes in the (kx, ky) wavenumbers.
  • Spectral energy-flux computations: The b-to-b energy flux Πbb(k0) is defined analogously for the full magnetic field and transfer through wavenumber k0.The magnetic-field decomposition parallels the velocity-field definition.
  • Spectral energy-flux computations: Computed fluxes show kinetic and magnetic energies cascading from large to small physical scales.The cascades are reported in Extended Data Fig. 2.

X. Relevance of the ϒ-dynamo to laboratory plasmas. · XI. Relevance of the ϒ-dynamo to astrophysical plasmas and BNS mergers.

The ϒ-dynamo shares key features with the Madison Dynamo Experiment and is suggested to operate in laboratory and astrophysical plasmas, including BNS mergers and the Sun. In BNS mergers, Kelvin–Helmholtz instability can generate strong large-scale fields in milliseconds from infinitesimal magnetic fluctuations.

  • X. Relevance of the ϒ-dynamo to laboratory plasmas.: The ϒ-dynamo and Madison Dynamo Experiment both exhibit large-scale vortical flows and jets, despite differences such as spherical versus Cartesian geometry.Both systems’ large-scale flows are externally maintained, by rotating impellers in MDE and mean-flow forcing in the simulations.
  • X. Relevance of the ϒ-dynamo to laboratory plasmas.: Both systems show identical EMF-component alignment, with an angle of ~90° between mean magnetic field and EMF, indicating a virtually zero α-effect.The MDE’s considerably large radial vorticity is consistent with the dominant radial EMF.
  • X. Relevance of the ϒ-dynamo to laboratory plasmas.: These similarities suggest that jet-driven dynamos occur in laboratory systems and in nature.The comparison specifically relates the reported ϒ-dynamo to the Madison Dynamo Experiment.
  • XI. Relevance of the ϒ-dynamo to astrophysical plasmas and BNS mergers.: The shear-flow-driven dynamo exists in a wide variety of astrophysical systems, with BNS mergers and solar magnetism presented as specific cases.The section addresses the relevance of the reported mechanism to astrophysical plasmas.
  • XI. Relevance of the ϒ-dynamo to astrophysical plasmas and BNS mergers.: BNS merger simulations contain a shear layer with half-width a≈10–15 m and flow velocity around U0≲0.1c, where curvature effects are very weak.The thin shear layer permits Newtonian Kelvin–Helmholtz simulations to model the local flow.
  • XI. Relevance of the ϒ-dynamo to astrophysical plasmas and BNS mergers.: The mechanism addresses BNS merger simulations’ use of dipolar initial fields of ~1013–1015 G, despite observations indicating initial total fields of <108–1010 G.The passage identifies both initial-field amplitude and topology as unsupported by astrophysical observations.
  • XI. Relevance of the ϒ-dynamo to astrophysical plasmas and BNS mergers.: Even with zero initial large-scale magnetic field, Kelvin–Helmholtz instability generates strong large-scale fields in milliseconds from infinitesimal magnetic fluctuations.The resulting fields are advected outward across the remnant surface and can affect electromagnetic emissions detected by LIGO-Virgo-KAGRA.
  • XI. Relevance of the ϒ-dynamo to astrophysical plasmas and BNS mergers.: For the Sun, the proposed model replaces the fine-tuned α-effect with a ϒ-effect from azimuthal vorticity in meridional circulation, contributing to the azimuthal EMF and generating poloidal fields.The mechanism is expressed as 𝓔∝ϒ∇´U.

XII. How does the jet-driven ϒ-dynamo differ from the classical and other dynamos?

The jet-driven ϒ-dynamo differs from classical dynamos by sourcing the mean magnetic field through turbulent cross-helicity aligned with a large-scale shear flow. Unlike the Ω-effect’s direct field straining, its self-consistent coupling to magnetic fluctuations produces exponential growth without requiring passive field components.

  • ϒ-dynamo mechanism: The generalized mean EMF includes ϒ∇´U, where ϒ measures alignment between turbulent velocity and magnetic-field fluctuations, rather than relying only on α and β terms.The mean field Bx evolves along the direction of the large-scale flow Ux.
  • ϒ-dynamo mechanism: The ϒ-dynamo equation couples only to the mean flow Ux, turbulent diffusion β, and ϒ, so passive components By or Bz need not be evolved.This contrasts with the α2-dynamo, whose loop couples Bx to another large-scale magnetic-field component.
  • Self-consistent growth: Because ϒ depends on magnetic-field fluctuations as well as turbulent flow, treating it as constant incorrectly predicts linear growth of Bx.The self-consistent evolution of ϒ measures turbulent cross-helicity and corrects this assumption.
  • Self-consistent growth: The coupled evolution equations for Bx and ϒ produce exponential solutions, with inhomogeneous mean flow requiring Fourier-space calculations.The detailed calculations are deferred to a forthcoming publication.
  • ϒ versus Ω: The ϒ-effect arises from averaged turbulent interactions in mean-field theory, whereas the Ω-effect directly represents large-scale flow straining the large-scale magnetic field.The Ω-effect is characterized by B·∇U and contributes zero in this work.
  • ϒ versus classical dynamos: The simulated turbulent flow and fields are non-helical and close to pure Alfvénic states, making α zero and explaining why the ϒ-effect dominates the mean EMF.In these states, β and ϒ are similar and are not affected by the non-kinematic effect in the same manner as α.
Loading 2608.12530v1…