1 Lorentz force in MHD and its relation to the Maxwell stress tensor One of the main simplifications of the MHD approximation is the fact that the fact that the electro- magnetic force on the plasma (Lorentz force) can be expressed entirely in terms of the magnetic field, with no reference to the electric field E or the current density J. In its ground form the Lorentz force density is J × B f = ; (1.1) L c where Faraday's law for slow phenomena allows the elimination of the current density J through c J = (r × B) : (1.2) 4π In addition, the magnetic field is divergence-free: r · B = 0 : (1.3) In what follows, the following vector/tensor identities may prove useful: r (A · B) = A × (r × B) + B × (r × A) + (A · r) B + (B · r) A (1.4) r · (A ⊗ B) = (r · A) B + (A · r) B : A and B are two arbitrary vectors with components Ai and Bi respectively , and A ⊗ B is the dyadic 'direct product" tensor of the vectors A and B, with tensor components (A ⊗ B)ij = AiBj. 1 a. Show, by whatever means you see fit (e.g. vector manipulation, writing things out in terms of components ···) and using the above equations, that the Lorentz force density can be written as f L = −r · M ; (1.5) with B2 B ⊗ B M ≡ I − (1.6) 8π 4π the Maxwell stress tensor. The tensor I ≡ diag(1 ; 1 ; 1) is the 3 × 3 unit tensor. For those who forgot: if we put x1 = x, x2 = y and x3 = z the cartesian components of the vector V ≡ r · M are: @Mji Vi = ; (1.7) @xj where summation over repeated indices is implied (i.e. the Einstein summation convention). b. Consider a magnetic field that has a fixed direction, chosen along the x-axis, but varies in strength in the direction perpendicular to the field (i.e. in the z-direction): B(x) = B(z) x^ : (1.8) Calculate the Maxwell stress tensor for this case, and show explicitly that the resulting Lorentz force is perpendicular to B, as the ground form (1.1) obviously demands. 2 2 A simple two-dimensional MHD equilibrium (Exam question in 2014!) Consider a static (V = 0) and time-independent (@=@t = 0) MHD equilibrium with gravity. The gravitational acceleration is g = −rΦ, with Φ(x) the gravitational potential. The magnetic field B(x) lies in the x − z plane. The plasma has a density ρ and a pressure P . We will assume that the plasma has a constant temperature T0, so that pressure and density are related through RT P (x) = 0 ρ(x) ; (2.1) µ with R the universal gas constant and µ = m=mH the mass of the atoms in units of the hydrogen mass. The force balance for this simple equilibrium reads: J × B rP = + ρ g : (2.2) c a. Assume that gravity acts along the negative z-axis, g = −g z^ (2.3) with g > 0. Also assume that locally the magnetic field is inclined at an angle θ(s) with respect to the vertical direction, see Figure 1. Here s measures the distance along a given magnetic field line. Give the component of the force-balance equation (2.2) along the magnetic field, and show that it can be written as dP = −ρg (2.4) dz regardless the local direction of the magnetic field. b. Show that the density and pressure along a field line vary as ρ(s) = ρ0 exp(−z(s)=H) ;P (s) = P0 exp(−z(s)=H) (2.5) and determine the constant scale height H. Here s is the length measured along the magnetic field. 3 z g q(s) s z(s) B x Figure 1: The geometry as it applies to question 2: a magnetic field B = (Bx ; 0 ;Bz) in the x − z-plane is locally inclined at an angle θ(s) with respect to the vertical (z-)direction. The quantity s parametrizes the length along a given magnetic field line. Gravity acts in the negative z-direction: g = −g z^. c. The magnetic field is divergence-free: r · B = 0. Show that the magnetic field that is generated by a scalar function A(x ; z) through the relations @A @A B (x ; z) = ;B (x ; z) = − (2.6) x @z z @x is always divergence-free. The function A(x ; z) is known as the flux function. d. Show that such a magnetic field has a current density J equal to J(x) = J(x ; z) y^, where c @2 @2 J(x ; z) = (∆ A) with ∆ ≡ + : (2.7) 4π 2 2 @x2 @z2 4 e. Now show the following: { The flux function A(x ; y) is constant along the magnetic field; { Any solution of the force balance equation must have rP + ρ rΦ = F (A) rA; (2.8) where F (A) satisfies ∆2A + 4π F (A) = 0 : (2.9) 5 3 A buoyant slender flux tube in an isothermal atmosphere Consider a thin strand of magnetic field lines (thin flux tube or flux rope) with a length much larger than its width. This quasi-cylindrical flux rope will serve as a simple model of a Solar Prominence, see Figure 2, a magnetically confined strand of gas suspended above the Solar surface, or of a magnetic loop connecting two small Sunspots (see Figure). The magnetic field inside the rope is B = B `^: (3.1) with `^ a unit vector tangent to the rope's axis. If the rope traces its trajectory in the x − z plane we can always write `^= (cos θ ; 0 ; sin θ) (3.2) with θ the angle between `^ and the x-axis. We take the cross section to be circular, with radius a and cross sectional area πa2. The rope hangs in an atmosphere with a constant temperature T . This atmosphere is stratified in the vertical direction due to gravity with a gravitational acceleration equal to g = −g z^. As a result, the atmospheric pressure Pe and density ρe varies with height z as Pe(z) = Pe0 exp(−z=H) ; ρe(z) = ρe0 exp(−z=H) : (3.3) Here RT H = (3.4) µg is the pressure scale height. There is no magnetic field in the atmosphere that surrounds the flux rope. As always, pressure and density are related by the ideal gas law: ρRT P = nk T = : (3.5) b µ In this assignment you will calculate the trajectory of this flux rope. The trajectory is curved in general, but can be shown to lie in a plane, which we take to be the x − z plane. You are also asked to calculate how quantities vary along the length of the rope under the assumption that the local radius of the quasi-cylindrical rope is small and satisfies a H. 6 Figure 2: A large flux rope above the Solar surface, as seen by the Solar Dynamics Observatory on March 30, 2010. This rope shows some twist, an effect that is neglected in this calculation. 7 Figure 3: The geometry of a slender rope of magnetic flux that extends from the Solar photosphere into the corona. The rope runs from foot point to foot point, and reaches a maximum height above the photosphere equal to zmax. The inclination angle of the rope trajectory (defined by the tangent vector `^) with respect to the horizontal is θ(z). The unit vector n^ points towards the center of curvature, and is perpendicu- lar to the rope's axis. The length along the rope is measured by the parameter `, with ` = 0 at the left foot point. The shape of the rope is determined by the balance between the buoyancy force f buoy due to the density difference between the rope and its surroundings, and the magnetic tension force f tens due to the magnetic field that runs through the rope. 8 a. If the flux rope is sufficiently thin, the temperature of the material in the rope (internal temper- ature Ti) and that of the atmosphere ( external temperature T ) will quickly equilibrate. At the edge of the flux rope the total pressure (magnetic + gas pressure) must be the same on both sides. Show that this implies that the mass density inside the rope satisfies the following relation: µB2 ρ = ρ − : (3.6) i e 8πRT You may neglect the variation of the atmospheric pressure around the axis of the tube as a H. Along the magnetic field (in the direction of the rope's axis) the Lorentz force (J × B)=c vanishes. Therefore the force balance in this direction only involves the gradient of the internal gas pressure Pi and gravity: @P `^· r P ≡ i = −ρ g = −ρ g sin θ : (3.7) i @` i k i Here θ is the angle between the tangent to the tube and the horizontal direction (the x-direction), ^ gk = g · ` is the component of the gravitational acceleration along the rope's axis and ` is the length measured along the rope trajectory. b. Use the geometry of the problem, the ideal gas law and the isothermal assumption Ti = T = constant to show that this equation is solved by choosing the density inside the flux rope to have the same exponential profile as quantities outside the flux rope: for a section of rope at height z one has ρi(z) = ρi0 exp(−z=H) (3.8) (Hint: Use that every position on the rope's axes can be uniquely parametrized by `, but also by z: there is a one-to-one relationship between the two!) 9 c. The total magnetic flux through the rope, 2 Φm = πa B (3.9) is a conserved quantity that can not vary along the length of the (curved) tube.