Typeset in LATEX 2ε July 26, 2017 A Mathematical Introduction to Magnetohydrodynamics
Omar Maj
Max Planck Institute for Plasma Physics, D-85748 Garching, Germany. e-mail: [email protected]
1 Contents
Preamble 3
1 Basic elements of fluid dynamics 4 1.1 Kinematics of fluids...... 4 1.2 Lagrangian trajectories and flow of a vector field...... 5 1.3 Deformation tensor and vorticity...... 14 1.4 Advective derivative and Reynolds transport theorem...... 17 1.5 Dynamics of fluids...... 20 1.6 Relation to kinetic theory and closure...... 24 1.7 Incompressible flows ...... 32 1.8 Equations of state, isentropic flows and vorticity...... 34 1.9 Effects of Euler-type nonlinearities...... 35
2 Basic elements of classical electrodynamics 39 2.1 Maxwell’s equations...... 39 2.2 Lorentz force and motion of an electrically charged particle...... 52 2.3 Basic mathematical results for electrodynamics...... 57
3 From multi-fluid models to magnetohydrodynamics 69 3.1 A model for multiple electrically charged fluids...... 69 3.2 Quasi-neutral limit...... 74 3.3 From multi-fluid to a single-fluid model...... 82 3.4 The Ohm’s law for an electron-ion plasma...... 86 3.5 The equations of magnetohydrodynamics...... 90
4 Conservation laws in magnetohydrodynamics 95 4.1 Global conservation laws in resistive MHD...... 95 4.2 Global conservation laws in ideal MHD...... 98 4.3 Frozen-in law...... 100 4.4 Flux conservation...... 105 4.5 Topology of the magnetic field...... 109 4.6 Analogy with the vorticity of isentropic flows...... 117
5 Basic processes in magnetohydrodynamics 119 5.1 Linear MHD waves...... 119 5.2 Nonlinear shear Alfv´enwaves...... 128 5.3 Magnetic field diffusion...... 129
A Energy conservation in extended MHD models 137
B Magnetic vector potential in MHD 141
References 145
2 Preamble
Magnetohydrodynamics is the theory of electrically conducting, neutral fluids in the low-frequency regime. The basic equations of magnetohydrodynamics (MHD) have been proposed by Hannes Alfv´en[1, 2], who realized the importance of the electric currents carried by a plasma and the magnetic field they generate. Alfv´encombined the equations of fluid dynamics with Faraday’s and Amp`ere’s laws of electrodynamics, thus obtaining a novel mathematical theory, which helped understanding space plasmas in Earth and planetary magnetospheres, as well as the physics of the Sun, solar wind, and stellar atmospheres. In fusion research, MHD is crucial to the understanding of plasma equilibria and their stability. Liquid metals and electrolytes, like salt water, can also be modeled by MHD equations. Besides the important physical applications, MHD equations exhibit a re- markably beautiful mathematical structure, with connections to geometry and topology that allows us to understand some of the dynamics of magnetic fields in plasmas in terms of topological ideas [3, 4, 5, 6, 7, 8, 9]. As a dynamical system MHD is an example of infinite-dimensional Hamiltonian system [10, 11]. The scope of this lecture is, in this regard, extremely limited. The goal is to introduce MHD equations in a reasonably self-contained way and to discuss some of their most important features. The style of this lectures is quite similar to the mathematical introduction to fluid dynamics by Chorin and Marsden [12] as the title of this note suggests. Particularly, we shall attempt to introduce the physics modeling in a mathematically precise albeit not always rigorous way. The physics literature on the subject is vast. As a reference for further read- ing, the book by Biskamp [13] provides a clear and comprehensive exposition, while the lectures by Schnack [14] offer a more gradual learning curve. For MHD of the solar atmosphere one can refer to Priest [15] as well as to Aschwanden’s book on the solar corona [16]. For an introduction to magnetohydrodynamics with emphasis on equilibria and stability of fusion plasmas one can refer to the books by Freidberg [17] and Zohm [18], while Goedbloed, Poedts and Keppens address applications to both astrophysical and fusion plasmas [19, 20]. A nice introduction to MHD with a broader perspective which includes applications to metals can be found in Davidson’s book [21]. On the mathematical side, MHD has received considerable attention from ap- plied mathematicians. Its rich mathematical structure has become a paradigm for the application of geometry and topology [9, 22] as well as for structure preserving discretization [23, 24, 25, 26]. As a system of partial differential equations, well-posedness of the Cauchy problem for MHD equations subject to appropriate boundary conditions have been studied first by Duvaut and Lions [27, in French], Sermange and Temam [28], Secchi [29] and more recently, by Chen and coworkers [30] and by Fefferman and coworkers [31, and references therein].
3 1 Basic elements of fluid dynamics
The basic understanding of fluid dynamics is an essential prerequisite to the study of MHD. We shall start by recalling the basic elements thereof, following Chorin and Marsden [12], cf. also Marsden and Hughes [32]. First, we define the physical quantities that describe the dynamical state of a fluid (kinematics) and continue with the equations of motion (dynamics).
1.1 Kinematics of fluids. Fluid dynamics is built on the basis of the con- tinuum hypothesis: A fluid is a distribution of matter occupying a certain region of the continuous three-dimensional space. The considered region of space is a domain (i.e., an open and connected subset) Ω ⊆ R3. We neglect the fact that any fluid is ultimately made of atoms and molecules as we are interested in studying its collective motion on a much larger spatial scale. With the continuum hypothesis, one needs to quantify how matter is dis- tributed in Ω at any time t in a certain time interval I ⊆ R. Thus, the first physical quantity of interest is the mass density, which is a non-negative time- dependent scalar field ρ : I × Ω → R, such that
ρ(t, x) ≥ 0 (1.1) gives the mass per unit of volume at time t ∈ I at the spatial location x ∈ Ω. By definition, the amount of mass contained in an arbitrary volume W ⊆ Ω (to be referred to as a control volume) is given by Z (mass in W at time t) = ρ(t, x)dx, W which implies that the mass density must be at least locally integrable. Physically we think of an infinitesimal volume of fluid centered around a point x ∈ Ω. The volume of this infinitesimal region is mathematically repre- sented by the Lebesgue measure dx in Ω and the mass is represented by the measure ρ(t, x)dx. Such infinitesimal portions of fluid are referred to as fluid elements. In plasma physics the mass density is often replaced by an equivalent non- negative scalar field referred to as the particle number density or simply number density which is defined in terms of the mass density by
n(t, x) = ρ(t, x)/m.
Here the fluid is regarded as a collection of particles that all have the same mass m; thus, the number of particles contained in a control volume W is given by 1 Z Z (number of particles in W at time t) = ρ(t, x)dx = n(t, x)dx. m W W According to this definition, the number of particles does not need to be an integer, due to the continuum hypothesis. Next we need to describe the motion of the fluid. We introduce a velocity field defined as a time-dependent vector field, u : I × Ω → R3 such that,
u(t, x) ∈ R3 (1.2)
4 gives the velocity of the fluid element at the point x ∈ Ω and time t ∈ I. The vector u(t, x) is referred to as the fluid velocity. The mass times the velocity of the fluid element, namely, ρ(t, x)u(t, x)dx, gives the linear momentum of the fluid element, hence Z (momentum in W at time t) = ρ(t, x)u(t, x)dx. W In addition to the quantities ρ(t, x) and u(t, x), that can be regarded as the counterpart in fluid dynamics of mass and velocity of a particle in mechan- ics, we need to specify another scalar field for the internal energy of the fluid element. Differently from a point-mass particle, a fluid element is a thermody- namical system that can undergo expansions and compressions, thus absorbing and releasing energy. The thermodynamic status of a fluid element is specified by the internal energy density U : I × Ω → R≥0, where R≥0 denotes the set of non-negative real numbers. If each fluid element is regarded as an ideal gas composed by n(t, x)dx particles, we can equivalently express the internal energy density in terms of a new variable. Specifically the laws of thermodynamics for a perfect gas allow us to write the internal energy of a fluid element as 3 U(t, x)dx = n(t, x)k T (t, x)dx, 2 B where kB is the Boltzmann constant and T : I × Ω → R is a strictly positive scalar field, such that T (t, x) > 0 (1.3) represents the local temperature of the fluid element at time t ∈ I and position x ∈ Ω. Therefore, the total energy carried by a fluid element is the sum of the kinetic energy associated to its motion plus the internal energy associated to its thermodynamics, namely, Z 1 2 3 (energy in W at time t) = ρ(t, x)u(t, x) + n(t, x)kBT (t, x) dx. W 2 2 When the law of a perfect gas does not apply we can still define T as above, but now it has the meaning of an effective temperature, which is just a measure of the internal energy and does not imply thermodynamical equilibrium. Summarizing, we shall describe the dynamical state of a fluid by the triple of functions ρ, u, T , where • the mass density ρ is a non-negative scalar field, • the fluid velocity u is a vector field, and • the temperature T is a positive scalar field. The equations of fluid dynamics are a system of partial differential equations governing the time evolution of (ρ, u, T ).
1.2 Lagrangian trajectories and flow of a vector field. Under appro- priate hypotheses, we can associate to any velocity field u : I × Ω → R3 a one-parameter family of maps Ft :Ω → Ω, where time t, varying in a possibly smaller interval Iε ⊆ I, is the parameter. Such a family of maps is referred to
5 as the flow of the vector field. It gives an equivalent description of the motion of the fluid, i.e., the vector field u and its flow Ft contain the same information on the fluid motion. In this section we shall define the flow and prove some of its basic properties. Although it is often overlooked in the physics literature, the flow is a key concept in the mathematical theory of fluid dynamics and thus of MHD. Let us start from a given velocity field u : I × Ω → R3. The associated flow Ft is constructed from the solution of the Cauchy problem dx(t) = u t, x(t), x(t ) = x . (1.4) dt 0 0 Physically, the solution t 7→ x(t) represents the trajectory of a fluid element as it moves with the fluid velocity from the initial position x0 ∈ Ω at time t = t0 ∈ I. Such curves are referred to as Lagrangian trajectories. Basic results from the theory of ordinary differential equations (ODE) guar- antee the existence and uniqueness of the solution of the Cauchy problem (1.4) at least for a short time. A compact account of results on ODEs can be found, for instance, in Marsden et al. [33] as part of the theory of vector fields. The first chapter of both H¨ormander’s[34] and Tao’s [35] lectures on nonlinear par- tial differential equations gives a very nice and compact overview of the theory. Specifically we have the following standard result that we recall without proof. Let us fix constants τ, ρ ≥ 0 such that the interval Iτ = [t0 − τ, t0 + τ] and 3 the ball Bρ(x0) = {x ∈ R | |x−x0| ≤ ρ} are contained in I and Ω, respectively, and let u be continuous and V = sup |u(t, x)| be the maximum velocity in the restricted domain Iτ × Bρ(x0). Theorem 1.1 (Local existence and uniqueness for ODEs). If u is continuous on Iτ ×Bρ(x0), where we have |u(t, x)| ≤ V , and satisfies the Lipschitz condition
|u(t, x) − u(t, y)| ≤ L|x − y|, with constant L ≥ 0 uniformly in time, then for any positive ε ≤ min{τ, ρ/V } 1 there exists a solution x ∈ C ([t0−ε, t0+ε]) of the Cauchy problem (1.4) and any 1 other solution x˜ ∈ C ([t0 −ε,˜ t0 +ε ˜]) must satisfy x(t) =x ˜(t) on the intersection of the domains. We can see that the upper limit of the domain of definition is determined by the minimum time ρ/V needed to traverse the ball Bρ(x0). The interval [t0 − ε, t0 + ε] is referred to as the lifespan of the solution. This has a physical significance: The lifespan of the solution is determined by how fast the trajectory can travel up to the boundary of the considered ball. In general the maximum lifespan depends on the initial condition x0. For instance, if the initial condition is very close to the boundary of Ω, ρ and thus ε can be rather small. We can at most refine a bit this result and make the lifespan of the solution uniform for all initial conditions in a small neighborhood of x0. This can be established as a corollary of the basic existence result applied to a smaller ball centered on x0: For all initial conditions y0 ∈ Bρ/2(x0), theorem 1.1 with x0 and ρ replaced by y0 and ρ/2, respectively, gives a solution of the Cauchy problem with initial condition x(t0) = y0; then such a solution is contained in Bρ(x0) and the lifespan is ≤ min{τ, ρ/(2V )} for all y0 ∈ Bρ/2(x0). Hence,
6 Corollary 1.2. Let u, t0, x0, ρ, τ, and V be as in theorem 1.1. Then there exists a neighborhood U ⊂ Bρ(x0) and 0 < ε ≤ min{τ, ρ/(2V )} such that for every 1 y0 ∈ U the Cauchy problem (1.4) has a solution x ∈ C ([t0 − ε, t0 + ε]). We shall however work under the assumption that the lifespan of Lagrangian trajectories is uniform on the whole domain Ω, i.e., we assume that there is an ε > 0 depending only on t0, such that for every initial condition x0 ∈ Ω there is a Lagrangian trajectory x : Iε → Ω with Iε = [t0 − ε, t0 + ε]). For a generic ordinary differential equation, this is a very strong assumption. For our problem, however, this is not so strong because, in practice, it just means that the domain Ω and the boundary conditions for the vector field u have been chosen properly, in the sense that “Ω contains the fluid”. We shall fix the initial time to be t0 = 0 and let Iε = (−ε, ε) the interval of existence of the Lagrangian trajectories. At this point we are ready to define the flow of the velocity field.
Definition 1.1 (Flow). For every t ∈ Iε the map Ft :Ω → Ω is defined by
x0 7→ x(t) = Ft(x0), where x(t) is the Lagrangian trajectory passing through x0 at the time t = 0. In addition this defines a map F : Iε × Ω → Ω given by F (t, x) = Ft(x). The uniqueness of the solution of the Cauchy problem for Lagrangian tra- jectories is essential in the definition of Ft. In fact, for Ft to be unambiguously defined we need that Ft(x) 6= Ft(y) implies x 6= y; it is not admissible that the same point is mapped into two different points. The reader can check that definition 1.1 is well posed in this sense because of the uniqueness of Lagrangian trajectories. From a physical point of view the flow describes the displacement of the fluid as time advances, i.e., given a control volume W ⊆ Ω, then Ft(W ) ⊆ Ω is the volume occupied by the fluid initially in W after it has evolved for a time t. In summary, we have constructed a one-parameter family of maps Ft which can be used in two different ways, namely,
• t 7→ Ft(x) is the Lagrangian trajectory passing through x at time t = 0;
• x 7→ Ft(x) is the displacement of the point x after a time t. As a consequence of the definition, the flow satisfies, cf. equation (1.4)
d F (x ) = u t, F (x ), dt t 0 t 0 (1.5) F0(x0) = x0, where the initial point x0 ∈ Ω is regarded as a parameter, so that we write a total derivative instead of a partial derivative and consider this an ordinary differential equation rather than a partial differential equation. We shall now establish a few key properties of the flow Ft that essentially descend from equation (1.5) and conclude this section with two examples.
Proposition 1.3 (Semi-group property). For every t, s ∈ Iε such that t+s ∈ Iε we have Ft+s = Ft ◦ Fs = Fs ◦ Ft.
7 Proof. For every x0 ∈ Ω, let us consider the Lagrangian trajectory x(t) corre- sponding to the initial condition x(0) = x0. By definition, Ft+s(x0) = x(t + s). 0 0 0 0 0 0 Let us define the functions t 7→ cs(t ) = x(t + s) and t 7→ ct(t ) = x(t + t). We observe that cs solves the Cauchy problem dc s = u(t0, c ), c (0) = x(s) = F (x ), dt0 s s s 0
0 hence, cs(t ) = Ft0 Fs(x0) = Ft0 ◦ Fs(x0). Then, Ft+s(x0) = x(t + s) = cs(t) = Ft ◦ Fs(x0). The second identity is obtained by applying the same argument to 0 the curve ct(t ). This property is usually referred to as the semi-group property of the flow and it has an immediate consequence.
Corollary 1.4. For every t ∈ Iε, Ft :Ω → Ω is invertible and the inverse is −1 given by Ft = F−t.
Proof. We have F0(x) = x for all x ∈ Ω and proposition 1.3 with s = −t gives x = F0(x) = Ft ◦ F−t(x) = F−t ◦ Ft(x) which shows that F−t is both a left- and a right-inverse of Ft.
Another basic result from ODE theory implies that Ft :Ω → Ω is Lipschitz continuous in Ω for all t ∈ Iε.
Proposition 1.5. If u is Lipschitz continuous uniformly in time and Ft :Ω → Ω is defined on Ω for all t ∈ Iε, then
L|t| Ft(x) − Ft(y) ≤ |x − y|e , for all x, y ∈ Ω and t ∈ Iε. Here L is the Lipschitz constant of u. Proof. Let us first consider the half-interval t ≥ 0. For every x, y ∈ Ω fixed, let h(t) = Ft(x) − Ft(y) and, by the Lipschitz condition for u,
dh(t) = u t, Ft(x) − u t, Ft(y) ≤ L|h(t)|. dt We actually need to control the derivative of the norm, rather then the norm of the derivative. With this aim we can estimate
1 d 2 dh(t) dh(t) 2 h(t) = h(t) · ≤ h(t) · ≤ Lh(t) , 2 dt dt dt and for t ≥ 0, d h d i h(t)2e−2Lt = h(t)2 − 2Lh(t)2 e−2Lt ≤ 0, dt dt hence h(t)2e−2Lt ≤ h(0)2 which is equivalent to the claim for t ≥ 0. As for the other half-interval t ≤ 0, let s = −t ≥ 0 and h(s) = F−s(x) − F−s(y) and we notice that dh(s) = −u − s, F (x) + u − s, F (y), ds −s −s and repeat the argument, integrating in the variable s.
8 We now know that Ft is a continuous transformation of Ω into itself for all t ∈ Iε. We also shall need to understand when Ft is differentiable and in those cases have a convenient way to compute its Jacobian matrix and determinant, namely,