540 likes | 713 Views
Non-helical MHD at 1024 3. Haugen, Brandenburg, & Dobler (2003, ApJ). Inverse cascade of magnetic helicity. argument due to Frisch et al. (197 5 ). and. Initial components fully helical:. and. k is forced to the left. Magnetic helicity. Maxwell eqns. Vector potential.
E N D
Non-helical MHD at 10243 Haugen, Brandenburg, & Dobler (2003, ApJ)
Inverse cascade of magnetic helicity argument due to Frisch et al. (1975) and Initial components fully helical: and k is forced to the left
Magnetic helicity Maxwell eqns Vector potential Uncurled induction eqn
Slow saturation Brandenburg (2001, ApJ 550, 824)
Periodic box, no shear: resistively limited saturation Brandenburg & Subramanian Phys. Rep. (2005, 417, 1-209) Significant field already after kinematic growth phase followed by slow resistive adjustment Blackman & Brandenburg (2002, ApJ 579, 397)
Magnetic helicity conservation Steady state, closed box Early times
Slow-down explained by magnetic helicity conservation molecular value!!
With hyperdiffusivity Brandenburg & Sarson (2002, PRL) for ordinary hyperdiffusion
Evidence from different simulations:strong fields only with helicity flux 3-D simulations, no mean-field modeling Forced turbulence in domain with solar-like shear Brandenburg (2005, ApJ 625, 539) Convective dynamo in a box with shear and rotation Käpylä, Korpi, Brandenburg (2008, A&A 491, 353) Only weak field if boxis closed
Nonlinear stage: consistent with … Brandenburg (2005, ApJ)
Best if W contours ^ to surface Example: convection with shear need small-scale helical exhaust out of the domain, not back in on the other side Magnetic Buoyancy? Tobias et al. (2008, ApJ) Käpylä et al. (2008, A&A)
To prove the point: convection with vertical shear and open b.c.s Magnetic helicity flux Käpylä, Korpi, Brandenburg (2008, A&A) Käpylä, Korpi, & myself (2008, A&A 491, 353) Effects of b.c.s only in nonlinear regime
Implications of tau approximation • MTA does not a priori break down at large Rm. (Strong fluctuations of b are possible!) • Extra time derivative of emf • hyperbolic eqn, oscillatory behavior possible! • t is not correlation time, but relaxation time with
Connection with a effect: writhe with internal twist as by-product a effect produces helical field W clockwise tilt (right handed) left handed internal twist both for thermal/magnetic buoyancy
… the same thing mathematically Two-scale assumption Production of large scale helicity comes at the price of producing also small scale magnetic helicity
Revised nonlinear dynamo theory(originally due to Kleeorin & Ruzmaikin 1982) Two-scale assumption Dynamical quenching Kleeorin & Ruzmaikin (1982) ( selective decay) Steady limit algebraic quenching:
General formula with magnetic helicity flux Rm also in the numerator
Mean field theory is predictive • Open domain with shear • Helicity is driven out of domain (Vishniac & Cho) • Mean flow contours perpendicular to surface! • Excitation conditions • Dependence on angular velocity • Dependence on b.c.: symmetric vs antisymmetric
Calculate full aij and hij tensors Original equation (uncurled) Mean-field equation fluctuations Response to arbitrary mean fields
Test fields Example:
Validation: Roberts flow SOCA SOCA result Brandenburg, Rädler, Schrinner (2009, A&A) normalize
Kinematic a and ht independent of Rm (2…200) Sur et al. (2008, MNRAS)
Scale-dependence: nonlocality cf talk by Alexander Nepomnyashchy
Time-dependent case Hubbard & Brandenburg (2009, ApJ)
From linear to nonlinear Brandenburg et al. (2008, ApJ) Use vector potential Mean and fluctuating U enter separately
Nonlinear aij and hij tensors Consistency check: consider steady state to avoid da/dt terms Expect: l=0 (within error bars) consistency check!
ht(Rm) dependence for B~Beq • l is small consistency • a1 and a2 tend to cancel • to decrease a • h2 is small
Application to passive vector eqn cf. Cattaneo & Tobias (2009) Verified by test-field method Tilgner & Brandenburg (2008)
Is the field in the Sun fibril? Käpylä et al (2008) with rotation without rotation
Takes many turnover times Rm=121, By, 512^3 LS dynamo not always excited