1 / 52

A Few Issues in MHD Turbulence

A Few Issues in MHD Turbulence. Alex Alexakis * , Bill Matthaeus % , Pablo Mininni^ , Duane Rosenberg and Annick Pouquet. NCAR / CISL / TNT * Observatoire de Nice % Bartol, U. Delaware ^ Universidad de Buenos Aires. Warwick, September 17, 2007 pouquet@ucar.edu.

verne
Download Presentation

A Few Issues in MHD Turbulence

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. A Few Issues in MHD Turbulence Alex Alexakis*, Bill Matthaeus%, Pablo Mininni^, Duane Rosenberg and Annick Pouquet NCAR / CISL / TNT * Observatoire de Nice % Bartol, U. Delaware ^Universidad de Buenos Aires Warwick, September 17, 2007pouquet@ucar.edu

  2. A Few Issues in MHD Turbulence:The strong and the weak

  3. * Introduction • Some examples of MHD turbulence in astro, geophysics & the lab. • Equations, invariants, exact laws, phenomenologies • * Direct Numerical Simulations, mostly for the decay case • Is there any measurable difference with fluid turbulence? • What are the features of an MHD flow, both spatially and spectrally? • Discussion • Accessing high Reynolds numbers for better scaling laws • Can modeling of MHD flows help understand their properties? • An example : The generation of magnetic fields at low PM • Can adaptive mesh refinement help unravel characteristic features of MHD? • * Conclusion

  4. Observations of galactic magnetic fields (after Brandenburg & Subramanian, 2005)

  5. Cyclical reversal of the solar magnetic field over the last 130 years • Prediction of the next solar cycle, because of long-term memory in the system (Dikpati, 2007)

  6. HINODE / Sunrise November 2006

  7. Surface (1 bar) radial magnetic fields for Jupiter, Saturne & Earth, Uranus& Neptune(16-degree truncation, Sabine Stanley, 2006) Axially dipolar Quadrupole~dipole

  8. Reversal of the Earth’s magnetic field over the last 2Myrs (Valet, Nature, 2005) Brunhes Jamarillo Matuyama Olduvai Temporal assymmetry of reversal process

  9. Experimental dynamo within a constrained flow: Riga Gailitis et al., PRL 84 (2000) Karlsruhe, with a Roberts flow (see e.g.Magnetohydrodynamics 38, 2002, special issue)

  10. W R H=2R W Taylor-Green turbulent flow at Cadarache Bourgoin et al PoF 14 (‘02), 16 (‘04)… Numerical dynamo at a magnetic Prandtl number PM=/=1(Nore et al., PoP, 4, 1997) and PM ~ 0.01(Ponty et al., PRL, 2005). In liquid sodium, PM ~ 10-6 : does it matter? Experimental dynamo in 2007

  11. Small-scale ITER (Cadarache)

  12. The MHD equations [1] • P is the pressure, j = ∇ × B is the current, F is an external force, ν is the viscosity, ηthe resistivity, v the velocity and B the induction (in Alfvén velocity units); incompressibility is assumed, and .B = 0. ______ Lorentz force

  13. The MHD invariants (=  =0) * Energy: ET=1/2< v2 + B2 > (direct cascade to small scales, including in 2D) * Cross helicity: HC= < v.B > (direct cascade) And: * 3D: Magnetic helicity: HM=< A.B > with B=  x A (Woltjer, mid ‘50s) * 2D: EA= < A2 > (+)[A: magnetic potential] Both HM and EAundergo an inverse cascade (evidence: statistical mechanics, closure models and numerical simulations)

  14. The Elsässer variablesz± = v ± b t z+ + z- .z+ = - P + +∆ z+ + -∆ z- + F± with  ± = [ ±]/2 and F± = FV ± FM Invariants: E± = < z ±2 >/2 = < v2 + B2 ± 2 v.B >/2

  15. Parameters in MHD • RV = Urms L0 / ν >> 1 • Magnetic Reynolds number RM = Urms L0 / η * Magnetic Prandtl number: PM = RM / RV= ν / η PM ishigh in the interstellar medium. PM is low in the solar convection zone, in the liquid core of the Earth, in liquid metals and in laboratory experiments. • Energy ratio EM/EV • Uniform magnetic field B0 • Amount of magnetic helicity HM

  16. Two scaling laws in MHD [1]à la Yaglom (1949), and Antonia et al. (1997) F(r) = F(x+r) - F(x) : structure function for field F ; longitudinal component FL(r) = F . r/ |r| <  z-/+L i  zi±2 > = - [ 4/d ] ±r in dimension d, with z± = v ± b,± = - dtE± = T± C, omitting dissipation and forcing and assuming stationarity [Politano and AP, 1998, Geophys. Res. Lett. 25 (see also Phys. Rev. E 57)] Note: <z+ z+ z- > ~ < ( v +  B)2 ( v -B) > ~ + Alack of equipartition between kinetic and magnetic energy is observed in many instances (e.g., in the solar wind)

  17. Two scaling laws in MHD [2] • In terms of velocity and magnetic field, these two scaling laws become: • < vLvi2 >+ <vLbi2 > -2 < bLvibi > =-(4/d) Tr • - <bLbi2> - < bLvi2 >+2 < vLvibi > =-(4/d) cr • with T = - dtET andc = - dtHc • Strong V regime, or strong B regime or Alfvénic (v~B) regime (Ting et al 1986, …); • in the latter case, v-B correlations play a dynamical role • (Politano+AP,GRL 25, 1998; see also Boldyrev, 2006) • Where such laws apply, the energy inputTcan be measured,e.g.in the solar wind

  18. Phenomenologies for MHD turbulence • If MHD is like fluids Kolmogorov spectrum EK41(k) ~ k-5/3 Or • Slowing-down of energy transfer to small scales because of Alfvén waves propagation along a (quasi)-uniform field B0: EIK(k) ~ (TB0)1/2k-3/2 (Iroshnikov - Kraichnan (IK), mid ‘60s) transfer ~ NL * [NL/A] , or 3-wave interactions but still with isotropy. Eddy turn-over time NL~ l/ul and wave (Alfvén) time A ~ l/B0 And • Weak turbulence theory for MHD (Galtier et al PoP 2000): anisotropy develops and the exact spectrum is: EWT(k) = Cwk-2 f(k//) Note: WT is IK -compatible when isotropy (k// ~k ) is assumed: NL~ l/ul andA~l// /B0 Or k-5/3(Goldreich Sridhar, APJ 1995) ? Ork-3/2(Nakayama, ApJ 1999; Boldyrev, PRL 2006) ?

  19. Spectra of three-dimensional MHD turbulence • EK41(k) ~ k-5/3 as observed in the Solar Wind (SW) and in DNS (2D & 3D) Jokipii, mid 70s, Matthaeus et al, mid 80s, … • EIK(k) ~ k-3/2 as observed in SW, in DNS (2D & 3D), and in closure models Müller & Grappin 2005; Podesta et al. 2007; Mason et al. 2007; Yoshida 2007 • EWT(k) ~ k-2 as may have been observed in the Jovian magnetosphere, and recently in a DNS, Mininni & AP arXiv:0707.3620v1, astro-ph(see later) • Is there a lack of universality in MHD turbulence (Boldyrev; Schekochihin)? If so, what are the parameters that govern the (plausible) classes of universality? The presence of a strong guiding uniform magnetic field? * Can one have different spectra at different scales in a flow? * Is there it a lack of resolving power (instruments, computers)? * Is an energy spectrum the wrong way to analyze / understand MHD?

  20. Recent results using direct numerical simulations and models of MHD

  21. Numerical set-up • Periodic boundary conditions, pseudo-spectral code, de-aliased with the 2/3 rule • From 643 to 15363 gridpoints • No imposed uniform magnetic field • Either decay runs or forcing • ABC flow (Beltrami) + random noise at small scale, with V and B in equipartition (EV=EM), or Orszag-Tang (X-point configuration), or Taylor-Green flow

  22. Energy transfer and non-local interactions in Fourier • space • Temporal evolution of maximum of current and • vorticity • Energy dissipation and scaling laws • Piling, folding & rolling-up of current and vorticity sheets • Alignment of v and B fields in small-scale structures • Energy spectra and anisotropy • Intermittency • Beyond Moore’s law

  23. Energy Transfer Let uK(x) be the velocity field with wave numbers in the range K < |k| < K+1 Q K K+1 Q+1 ^ Isotropy ^ Sharp filters Fourier space

  24. Rate of energy transfer in MHD10243 runs, either T-G or ABC forcing(Alexakis, Mininni & AP; Phys. Rev. E 72, 0463-01 and 0463-02, 2005) Advection terms R~ 800 All scales contribute to energy transfer through the Lorentz force This plateau seems to be absent in decay runs (Debliquy et al., PoP 12, 2005)

  25. Rate of energy transfer Tub(Q,K) from u to b for different K shells K= 10 K= 20 K= 30 The magnetic field at a given scale receives energy in equal amounts from the velocity field from all larger scales (but more from the forcing scale). The nonlocal transfer represents ~25% of the total energy transfer at this Reynolds number (the scaling of this ratio with Reynolds number is not yet determined).

  26. Jmax for the Orszag-Tangflow Mininni et al., PRL 97, 244503 (2006)(res.:1283 to 10243) • Growth of the current maximum: linear phase, followed by an algebraic tnregime with n ~ 3

  27. Jmax for a random flow, resolutions up to 15363 grid points (RV from 690 to 10100) • Growth of the current maximum: linear phase followed by an algebraic tn regime with n ~ 3 OT flow

  28. AMR in incompressible 2D - MHD turbulence at R~1000 • Adaptive Mesh Refinement using spectral elements of different orders • Accuracy matters when looking at Max norms, here the current, although there are no noticeable differences on the L2 norms (for both the energy and enstrophy) Rosenberg et al., New J. Phys. 2007

  29. Energy dissipation rate in MHD for several RVOT- vortex Low Rv High Rv Orszag-Tang simulations at different Reynolds numbers (X factor of 10) • Is the energy dissipation rate Tconstant in MHD at large Reynolds number (Mininni + AP, arXiv:0707.3620v1, astro-ph), as presumably it is in 2D-MHD in the reconnection phase? There is evidence of constant  in the hydro case(Kaneda et al., 2003)

  30. Scaling with Reynolds number of times for max. current & dissipation Time Tmax(1) at which globalmaximum of dissipation is reached in (ABC+ random) flow and Time Tmax(2) at which the current reaches its firstmaximum Both scale as Rv0.08

  31. Scaling with Reynolds number of dissipation

  32. MHD decay simulation @ NCAR on 15363gridpoints Visualization freeware: VAPORhttp://www.cisl.ucar.edu/hss/dasg/software/vapor Zoom on individual current structures: folding and rolling-up Mininni et al.,PRL 97, 244503 (2006) Magnetic field lines in brown At small scale, long correlation length along the local mean magnetic field (k// ~ 0)

  33. Recent observations (and computations as well) of Kelvin-Helmoltz roll-up of current sheets Hasegawa et al., Nature (2004); Phan et al., Nature (2006), …

  34. Current and vorticity are strongly correlatedin the rolled-up sheet Current J2 Vorticity 2 15363 run, early time

  35. V and B are aligned in the rolled-up sheet, but not equal (B2 ~2V2): Alfvén vortices?(Petviashvili & Pokhotolov, 1992. Solar Wind: Alexandrova et al., JGR 2006) Current J2 cos(V, B) 15363 run, early time

  36. Hc : Velocity - magnetic field correlation PdFs of cos(v,B): • Flow with weak normalized total cross helicity Hc • Flow with strong Hc Matthaeus et al. arXiv.org/abs/ 0708.0801

  37. Velocity - magnetic field correlation [3] • Local map in 2D of v & B alignment: |cos (v,B)| > 0.7 (black/white) (otherwise, grey regions). even though the global normalized correlation coefficient is ~ 10-4 Quenching of nonlinear terms in MHD (Meneguzzi et al., J. Comp. Phys. 123, 32, 1996) similar to the Beltramisation (v // ) of fluids

  38. Vorticity =xv & Relative helicity intensity h=cos(v,) • Local v- alignment (Beltramization). Tsinober & Levich, Phys. Lett. (1983); Moffatt, J. Fluid Mech. (1985); Farge, Pellegrino, & Schneider, PRL (2001), Holm & Kerr PRL (2002). --> no mirror symmetry, together with weak nonlinearities in the small scales Blue h> 0.95 Red h<-0.95

  39. Strong relative magnetic helicity (~  1): change of topology across sheet Current J2 cos(A, B) , with B=xA 15363 run, early time

  40. Current at peak of dissipation Zoom Global view

  41. MHD scaling at peak of dissipation [1] Energy spectra compensated by k3/2 Solid: ET Dash: EM Dot: EV Insert: energy flux Dash-dot: k5/3-compensated LintM ~ 3.1 M ~ 0.4

  42. MHD scaling at peak of dissipation [2] • Anomalous exponents of structure functions for Elsässer variables, with isotropy assumed (similar results for v and B) Note 4 ~ 0.98  0.01, i.e. far from fluids and with more intermittency

  43. MHD decay run at peak of dissipation [3] Isotropy ratio RS = S2b / S2b// Isotropy obtains in the large scales, whereas anisotropy develops at smaller scales RS is proportional to the so-called Shebalin angles LintM ~ 3.1, M ~ 0.4

  44. (Mininni + AP, arXiv:0707.3620v1, astro-ph) LintM ~ 3.1 , M ~ 0.4 Structure functions at peak of dissipation [4] L1/2 -compensation of S2b  and //structure fns. Flat atlarge scales, with equipartition of the  and // components, hence E(k) ~ k-3/2 , isotropicIK (1965) spectrum Solid: perpendicular Dash: parallel Insert: l 2/3-compensated

  45. (Mininni + AP, arXiv:0707.3620v1, astro-ph) LintM ~ 3.1 , M ~ 0.4 Structure functions at peak of dissipation [5] Structure function S2 , with 3 ranges: L2 (regular) at small scale L at intermediate scale, as for weak turbulence: Ek~k-2. Is it the signature of weak MHD turbulence? L1/2 at largest scales (Ek~ k-3/2) Insert: anisotropy ratio Solid: perpendicular Dash: parallel

  46. Solid: parallel, Dash: perpendicular, LintM ~ 3.1, M ~ 0.4 Solid: perpendicular Dash: parallel

  47. Kolmogorov-compensated EnergySpectra: k5/3 E(k) Navier-Stokes, ABC forcing Small Kolmogorov k-5/3 law (flat part of the spectrum) K41 scaling increases in range, as the Reynolds number increases • Bottleneck at dissipation scale Solid: 20483, Rv= 104, R~1200 Dash: 10243, Rv=4000 Kolmogorov k-5/3 law Linear resolution: X 2 Cost: X 16 Mininni et al., submitted

  48. Evidence of weak MHD turbulence with a k-2 spectrum • in Reduced MHD computations (Dmitruk et al., PoP 10, 2003) • in the Jovian magnetosphere (Saur et al., Astron. Astrophys. 386, 2002)

  49. Extreme events at high Rv unraveled by high-resolution runs (grid from 483 to 15363 points)

More Related