Lecture 2 Basic MHD Equations MHD =

Maxwell’s Equation: v µ 1 ∂E v v v ρ : charge density ∇ × B = j + 2 c ∂t v  v v ε j : ∇ ⋅ B = 0  = 8.854×10−12 Fm−1 magnetic permittivity  v µ 0 v v ∂B π ∇ × E = − = 4 ×10−7 Hm−1 magnetic permeability  ∂t 0 ∗ 1  v ρv c = = 3×108 m / s speed of light ∇ ⋅ E = µ  ε 0ε 0 v v B H = µ If v << c non-relativistic v v D = εE Electric displaceme nt v v ∇× B = vj µ kT 1 λ = ( ) 2 << l in plasma d 4πne2 scale v v OHM' s law : vj = σ (E + vv × B) v E in the frame of reference moving with σ electric conductivity ohm m-1

A generalization of Ohm’s law e− p + 3-fluid model for , and neutral atom n ne n p

v v v ∇pe me ∂j v vv vv 1 1 v v v nee(E0 + ) = ( + ∇ ⋅ (vj + jv ) + + )Bj + j × B nee c ∂t Ω ei Ω en f 2Ωτ v v v v + in []∇P × B − ()vj × B × B B e v v v E = E + vv × B 0 τ eB τ Ω = electron gyration frequency me σ σ 2 v v v v f Ω in v v v If ignore electron inertia & E0 = j + j ×σB − ( j × B)× B nee τneeB σ pressure gradient: 2 −1 nee me = τ −1 −1 electron conductivity at B = 0 σ ei +τ en v v j // B: 0 E0 = j σ if v v v vj ⊥ B: E = vσj + 3 vj × B σ 3 0 n e σ e = cowling conductivity 3 τ 1+ f 2 BΩ in nee σ v v v B× E vj =σ E + 0 Solve above equation 1 0 2 B σ σ B = 3 = 3 1 B 2 1 2 nee σ1+ ( 3 σ) σ nee σ If gas is fully ionized n =0 σ n e 2 a = σ = Ω = = e ei 1 2 2 2 ei 1 3 1+ Ω ei me v τ σ Compare to vj = σE τ σ There is an additional term normal to v and v , it is calledσ Hall current B E σ τ Induction Equation  v v Ohm’s Law:  ∇× B = vj v v  v ∂B v v v j  v j v v = −∇×(−v × B + ) ⇐ E µ= − v × B ∂t σ  v ∂B v v σ = −∇× E η ∂t v v v v v v v v v v v v v v = ∇×(v × B) − ∇×η( ∇× B) ⇐µσ∇×(∇× B) = ∇(∇ ⋅ B) − ∇ ⋅∇B v ∂B v 2 v 1 = ∇×(vv × B) +η∇ B magnetic diffusivity ∂t = µ v v v ∇× B v v j j = E = vv × B + σ r v For condition of the sun : E ≈ vv × B

3 A typical active region, B=100G, v=10 m/s, E0=v0 B0=10 V/m If no motion, 1 jµ= = 8×10−4 Am−2 l =107 m σ l =103 ohmm−1 −7 −1 E0 = j0 /σ = 8×10 Vm Electric conductivity Fully ionized collision- dominated plasma n e2 σ = τe ei me 3/ 2 6 T Coulomb logarithm lnΛ is Collision time τ ei = 0.266×10 ne ln Λ between 5 and 20

T 3/ 2 σ =1.53×10−2 ohmm−1 ln Λ Diffusion constant η = 5.2×107 (ln Λ)T −3/ 2 m2s −1

typical value in solar chromosphere or coronaη: = 109T −3/ 2 nk T for turbulent plasma, anomalous collision time :τ * = w −1 B pe w w : plasma frequency σ pe n e2τ * anomalous conductivity : * = e me Condition for this is electron conduction speed >thermal speed Plasma equations ρ ∂ continuity equation : + ∇ρ⋅( vv) = 0 ∂t ρ Dvv ∂vv v v equation of motion : = ρ( + v ⋅∇vv) = −∇P + vj × B + F Dt ∂t D : derivative of following the motion Dt v vj × B : Lorentz force v v v F = Fg + Fv viscosity : coefficient viscosity v 2 v Fv = ∇ v ρν ν T 5/ 2 ρν = 2.21×10−16 kgm−1s −1 ln Λ

P = nkT n = np + ne Perfect gas law Energyρ equation: ζ • Thermal conduction : or Ds ρT = − v v Dt q = −K∇T ρ v De p D K thermal conduction tensor = − ρ = −ζ Dt Dt • Radiation : l = ∇⋅q = −K ∇2T e: internal energy, s: entropy r r r Kr heat diffusivity υ optical thick, this is valid; Ideal gas: ρ optical thin, lr = nenH Q(T) e = c υT • ρε Heating : K ρc N + 2 c = c + γ = = H = + H r + H w m c N υ nuclear + viscous disipatation + waves N: # of degrees of freedom • Energetics : j 2 v v Different terms for ζ E ⋅ j = + vv ⋅ j × B σ A few dimension less parameters γ

l0V0 Reynolds number : Re = η

l0V0 Magnetic Reynolds number : Rm = coupling between flow and magnetic fields

Solar atmosphere: Rm >> 1 coupling is strong

V0 rp0 Machβ number : M = Cs = speed of sound wave Cs ρ 0 V Alfven mach number : M = 0 β V = 0 Alfven speed β µ A V A ( )1/ 2 A µρ 0

2 p0 Gas plasma : β = 2 pressure B0 Magnetic <<1 low magnetic dominant - corona β V0 inertial possby number : R0 = l0Ω Coriolis

Rm γ Viscoucity Prandtl number : pm = = Re η diffusion Example. Sunspot 7 3 -6 -1 l0 = 10 m v0 = 10 m/s Ω = 10 s 4 20 -3 3 T0 = 10 k n0 = 10 m β 0 = 10 G

7 4 Rm = 3×10 , Cs = 2×10 m/s (M = 0.05) 5 12 VA = 3×10 m/s Re =10 R0 =100 More about Induction Equation v ∂B v 2 = ∇×(vv × B) +η∇ B ∂t l V v 0 0 1st term can be dropped, ∂B 2 •if Rm = << 1 =η∇ B Ædiffusion equation η ∂t 3 τ η 2 2 2 diffusion time scale l0 −8 l0 T d = =1.9×10 η lnΛ 6 15 -3 7 14 T =10 K =10 m l0 =10 m τ d =10 S

For solar flare τ d = 100ρ to 1000S l0 = 100 to 1000 m ρ v •if Rm >>1 frozen-in condition, B Moves with flow -- generally true in solar condition v v ρ D B B D ( ) = ( ⋅∇)vv is similar to fluid motion = −ρ∇ ⋅v prove: ∂t Dt 2 ηt1 2 ηt2

A diffusing current sheet: (a): the variation with time of the strength; (b): a sketch of the magnetic field lines at three times µ

Lorentz Force v 2 v 1 v v v v B B F = vj × B = (∇× B)× B = (B ⋅∇) µ − ∇( ) L 2µ magnetic tention magnetic pressure

Example (1) uniform field v D× B v vj = = 0 F = 0µ µ L 2 v v B e v v B e2x (2) B = B e× yv j = 0 × zv j × B = − 0 xv 0 µ pressure term v zˆ (3) B = −yxˆ + yˆ vj = µ v xˆ (vj × B) = − tension term y=0 µ

(4) X-type neutral point

2 2 B0 = yxˆ + xyˆ (solution y − x = const)

T balances with P, but unstable. (a) (b) The magnetic pressure P and tension T forces due to: (a) A uniform field; (b) a unidirectional field whose strength increases along the x-axis. The magnetic field lines near an X- type neutral point in equilibrium with no current

The resultant magnetic force (R) due to a symmetrically curved field (Equation 2.59) Flux tubes and current sheets µ Flux tubes F = ∫ B ⋅ds =constant B2 P + = constant (dark sunspots) 2 current sheet : µ B varies a lot in a very short spatical scale L B current density j = L η if v = 0, current sheet diffuses at a speed of /l, η : magnetic diffusivity enhanced plasma pressure in the center of the sheet expels material

away at a speed of VA ,reconnection A section of a magnetic

(a): a current sheet in the yz plane across which the magnetic field rotates from B1 to B2; (b): a section across a neutral current sheet in the centre of which the magnetic field vanishes and the plasma pressure is P0; ©: the reconnection of magnetic field lines by their passage through a current sheet. The (a) central sheet bifurcates into two (b) (c) pairs of slow shocks. Homework

1. Using reasonable parameters in solar photosphere, chromosphere, corona and in sunspots, calculate the following parameters: Alfven speed, sound speed, and plasma ß.

2. Refer to Equation 2.3, assume a reasonable diffusion constant, and magnetic field strength, compute magnetic topology as a function of time using the induction equation (ignore the velocity term).