Lecture IV: Magnetized plasma stability

Size: px
Start display at page:

Download "Lecture IV: Magnetized plasma stability"

Transcription

1 Outline Lecture IV: Magnetized plasma stability Wojciech Fundamenski UKAEA Fusion, Culham Science Centre, Abingdon and Plasma Physics Group, Imperial College, London 22nd May 2009

2 Outline Part I: 1 Hydrodynamic waves and instabilities 2 Ideal MHD waves in a uniform plasma in a stratified plasma Ideal in a confined plasma Ideal in a general screw pinch Flute-reduced MHD 3 Non-homogenous shear Alfvén waves Current driven ideal MHD instabilities: kink modes Pressure driven ideal MHD instabilities: ballooning modes Resistive MHD instabilities: tearing modes

3 Outline Part II: Drift-fluid and kinetic waves and instabilities 4 Drift waves and instabilities 5 Kinetic waves and instabilities

4 Part I

5 Something to consider... Let the nature of a fluid be assumed to be such that of its parts, which lie evenly and are continous, that which is under lesser pressure is driven along by that under greater pressure. Archimedes (c. 500 BC) Dynamical equilibrium does not guarantee dynamical stability. In general, an equilibrium is said to be stable if the system remains bounded, i.e., confined to its neighbourhood, after being subjected to a small perturbation. There are many different classifications of stability, e.g., linear stability implies small perturbations, while non-linear stability allows for perturbations of arbitrary size.

6 Aspects of dynamical stability The former is equivalent to spectral stability, which occurs when all eigenvalues of the linearized dynamical operator have real parts which are positive or zero. The most general formulation of linear stability is the energy principle, which states that an equilibrium point is stable when it represents a minimum of the potential energy of the system. Most dynamical systems have both stable and unstable equilibria, e.g., the lowest and highest position of a pendulum. Having identified an equilibrium point, one should next investigate its stability properties and, if the point proves unstable, the physical mechanism responsible for the instability.

7 Hydrodynamic stability of a stratified fluid By way of introduction to plasma instabilities, let us consider the stability properties of a stratified neutral fluid in the presence of a gravitational field g = gê g. Note that MHD differs from hydrodynamics only by the appearance of the Lorentz force term J B in the momentum conservation equation and by an additional evolution equation for B, including the resistivity η. Moreover, the presence of B provides a preferred direction. When B 0, or β 1, isotropy is restored and MHD reduces to hydrodynamics.

8 Hydrodynamic (neutral fluid) equations The fluid is assumed to conserve mass, momentum and energy, d t ρ + ρ V = 0, (1) ρd t V + p + π = ρg, (2) d t p + γp V + (γ 1)[π : V + q] = 0, (3) where d t = t + V is the convective time derivative, γ = c p /c v is the ratio of heat capacities at constant pressure and volume, π is the viscous stress and q is the heat flow. The momentum equation, albeit with a specific formulation of the viscosity tensor, is known as the Navier-Stokes equation.

9 Collisional closure of hydrodynamic equations In Lecture V we will see that π and q can be related to the velocity and temperature gradients, respectively, as π = µw, q = κ T, T = p/n = pm/ρ, (4) W = 2{ V} 2 3 ( V)I, W ij = xj V i + xi V j 2 3 δ ij V,(5) W = 2 V ( V), { V} = 1 2 ( V + V ), (6) W : V = W : { V} = 2{ V} : { V} 2 3 ( V)2, (7) where µ is the dynamic viscosity, κ is the heat conductivity and W is the (traceless) rate-of-strain tensor ; it is also customary to introduce the kinematic viscosity, ν = µ/ρ, and heat diffusivity, χ = κ/n = κm/ρ. Note that since W = W, W : V = W : V and { V} = { V}, then W : V = 1 2 (W + W ) : V = 1 2 (W : V + W : V ) = W : { V}.

10 Hydrostatic equilibrium and the thin-layer approximation The equilibrium conditions of hydrostatics are obtained by setting d t = 0 and V = 0 in (1)-(3), which yields p = ρg, κ T = const. (8) If the flow is non-dissipative (ν = 0 = χ), then (1)-(3) reduce to d t ρ + ρ V = 0 = d t T + γt V, (9) d t V + ρ 1 p = g, (10) where (9) yields the polytropic relation, p ρ γ, or T ρ γ 1. Combined with (8), it leads to the so called byrotropic equilibrium ρ(z) = ρ(0)(1 (1 γ 1 )z/l g ) 1/(γ 1), where z is the vertical distance and L g = p(0)/ρ(0)g is a characteristic length scale. When γ = 1 one finds ρ(z) = ρ(0) exp( z/l g ). A version of (1)-(3) for a layer of thickness L L g is known as the thin layer, thin atmosphere, or Boussinesq, approximation.

11 Acoustic (sound) waves Together with (10) this relation allows the propagation of compressible (acoustic, sound) waves with a phase velocity equal to the sound speed, V S, ω 2 = (kv S ) 2 v ph ω k = V S = γp γt ρ = m. (11) Acoustic waves are longitudinal disturbances of the form ρ 1, p 1, V 1 k exp(ik x iωt), where k is the wave vector and ω is the wave frequency. Since their group velocity, v gr k ω = V S is independent of k, acoustic waves are non-dispersive a fact which can be ultimately traced to the absence of an intrinsic length scale in hydrodynamics. Acoustic waves, which are the only waves possible in a uniform neutral fluid, are the primary means of establishing and restoring the hydrostatic equilibrium (8).

12 Incompressible flows ( V = 0) and the vorticity eq n If the flow is incompressible ( V = 0), which pertains when the characteristic flow velocity is small compared to the sound speed, i.e., when the Mach number, M S = V /V S 1, then (1) and (3) reduce to d t ρ = 0, d t T = (γ 1)[2νm{ V} : { V} + χ 2 T ]. (12) Moreover, the Navier-Stokes equation (2) can then be replaced by, d t Ω = Ω V + ν 2 Ω, Ω V, (13) where Ω is the fluid vorticity. The vorticity equation (13) is obtained by taking the curl of (2), which eliminates the pressure gradient. Here we used the fact that V V = 1 2 (V V) V V = 1 2 V 2 V Ω and that (V Ω) = V Ω + Ω V for incompressible flows.

13 Kelvin s vorticity theorem It is worth noting that V may be computed a posteriori from 2 V = Ω. Since Ω is just the vorticity, (13) is known as the vorticity equation in the fluid frame of reference; it may also be written in the laboratory frame as t Ω = (V Ω) + ν 2 Ω. (14) The above equations imply that, for incompressible flows, vorticity can only be created/destroyed by vortex stretching Ω V and viscous dissipation ν 2 Ω. In the absence of viscosity, the total circulation of the flow must remain constant and the vorticity is effectively frozen into the flow, d t Ω ds = d t V ds = d t V dl = 0, (15) S a result known as Kelvin s vorticity theorem. S C

14 Incompressible and non-dissipative flows As a consequence of (15), an initially irrotational flow (Ω = 0) must remain so at all times. The velocity can then be written as the gradient of some potential (V = φ) and the incompressibility condition, V = 0, may be recast in the form of Laplace s equation, 2 φ = 0. Finally, if the flow is both incompressible and non-dissipative, then d t ρ = 0 = d t T, d t V + ρ 1 p = g, d t Ω = Ω V. (16)

15 MHD version of the vorticity eq n There is a clear resemblance between the vorticity, Ω = V, the current density, µ 0 J = B and the magnetic field, B = A. This similarity is based on a general version of Kelvin s theorem, which states that d t S Q ds = 0 for any vector field satisfying t Q = (V Q). Hence, the conservation of magnetic flux in ideal MHD, and of circulation in ideal HD, stems from same vector relation. For future reference we also write down the MHD version of (13), d t Ω = Ω V + (B J J B)/ρ + ν 2 Ω, (17) which shows that the Lorentz force, J B, provides a source of vorticity.

16 Rayleigh-Taylor (R-T) instability Consider two fluids having different densities and zero relative velocity. This system has two possible hydrostatic equilibria: the heavy fluid is 1 above (ρ a > ρ b ), 2 below (ρ a < ρ b ), the lighter fluid. Of these, only (2) is stable, since it clearly represents the lowest energy state of the system. In contrast, the gravitational potential energy stored in (1) can be converted into kinetic energy, driving the heavier fluid down, and the lighter up, in some complicated 3D motion. The physical mechanism leading to interchange of the lighter and heavier fluids is the bouyancy force, g ρ, experienced by elements of the lighter fluid, and the process is known as an interchange, or Rayleigh-Taylor (R-T), instability.

17 Rayleigh-Taylor (R-T) instability Source: MIT

18 Rayleigh-Taylor (R-T) instability Source: DOE SCIDAC

19 Brunt-Vaisala frequency It also applies to a stratified fluid in which the density increases gradually with height (ê g ρ < 0), e.g., in a fluid heated from below, provided the fluid expands with increasing temperature, which is the case for nearly all fluids, e.g., for an ideal gas, in which p ρt, one finds ln ρ/ ln T = 1 and α ρ = 1/T > 0. The instability occurs when the ê g T is super-adiabatic, i.e., when the specific entropy, s = (c p /γ) ln(pρ γ ), decreases with height, N 2 g ( ln T (1 1/γ) ln p) = g ln(pρ γ )/γ > 0, (18) where N is the Brunt-Vaisala frequency of internal surface waves, see below.

20 Rayleigh-Benard (R-B) convection Combined with equilibrium (8) and ideal gas law, this gives the Schwarzschild instability condition, ê g T > (1 1/γ)mg. (19) The above variant of the R-T instability, in which ρ is caused by T, is known after Rayleigh and Benard (R-B). In the presence of constant heating, (3) requires constant energy flow from the bottom to top of the system. When (19) is satisfied, the system becomes R-B unstable, leading to upward flow of warm, light fluid, and downward flow of cold, heavy fluid. The resulting thermally driven flow is known as R-B, or thermal, convection.

21 Rayleigh-Benard (R-B) convection Sources: U.Goettingen (left), Van Dyke (right)

22 Kelvin-Helmholtz (K-H) instability If the light and heavy fluids have a finite relative velocity along their horizontal interface, ê g (V a V b ) 0, then the kinetic energy of this flow can amplify vertical perturbations. In the non-linear phase of this process, the surface ripples evolve into complicated vortices, destroy the smooth interface and effectively mix the two fluids. This mechanism, known as the Kelvin-Helmholtz (K-H) instability, applies equally well to a sheared flow, e.g., when the horizontal velocity changes with height, ê g V 0.

23 Kelvin-Helmholtz (K-H) instability Sources: UCAR

24 R-T and K-H stability criteria The linear stability criteria for both R-T and K-H instabilities in ideal, incompressible fluids, can be derived by linearising (16) around the hydrostatic equilibrium (8), and assuming small, vertical perturbations of the two-fluid interface, ξ exp(ikx iωt), where x is the distance and k the wave number along the surface. Since the initial flows are irrotational (Ω = 0) everywhere, except at the interface, which is thus a thin vorticity sheet, the perturbed velocity may be written as V 1 = φ 1. Denoting quantities above and below the interface by subscripts a and b, respectively, and the vertical distance away from the surface as z, we choose φ 1a and φ 1b so that V 1a and V 1b decrease exponentially with z : φ 1a exp(ikx + kz iωt) and φ 1b exp(ikx kz iωt).

25 Surface gravity waves Inserting these quantities into (16), introducing the effect of surface tension, T S by adding the term (T S /ρ) 2 x ξ on the left hand side of (10), and matching the conditions at the interface, yields the phase velocity of surface gravity waves (Choudhuri, 1998), ω k = ρ [ ] bv b + ρ a V a g ρb ρ a ± + k2 T S ρ aρ b (V b V a ) 2 ρ b + ρ a k ρ b + ρ a g(ρ b + ρ a ) (ρ b + ρ a ) 2. (20) Since all perturbed quantities are assumed to vary with time as exp( iωt), the real part of ω represents the frequency of sinusoidal oscillations, while the imaginary part is the rate of exponential growth or decay, i.e., when ω = Re(ω) 0, the wave propagates along the surface with velocity (20), when Im(ω) < 0, it is quickly damped, and when Im(ω) > 0 it quickly grows without bound. The last case is just the instability we seek.

26 R-T instability criterion The condition for R-T instability follows from (20) with V a = V b, ρ b + k 2 T S /g < ρ a, (21) which agrees with the intuitive result (ρ b < ρ a ) for T S = 0 and shows that surface tension has a strong restoring effect on small scale (large k) perturbations. To appreciate this effect, imagine that an R-T unstable system is contained in a vertical cylinder. As the diameter d of the cylinder is reduced, k min = 2π/d increases until, for a sufficiently small d, k 2 min > g(ρ a ρ b )/T S is achieved and the R-T instability is suppressed. Indeed, since the volume to surface ratio scales linearly with d, one would expect surface tension to dominate at sufficiently small scales.

27 K-H instability criterion Similarly, we obtain the condition for the K-H instability, [ ρb ρ a g k + kt S g ] < ρ aρ b (V b V a ) 2, (22) (ρ b + ρ a ) which indicates that an R-T stable system (ρ b > ρ a ) becomes K-H unstable for sufficiently large velocity discontinuity V a V b. Although short scales (large k) are more unstable, they are also more effectively stabilized by surface tension, so that the K-H instability is most likely to occur at some intermediate k, determined from (22).

28 Internal gravity waves In a stratified fluid the surface gravity waves are replaced by internal gravity waves, and the stabilising effect of surface tension by that of viscosity and heat diffusivity. The dispersion relation for internal gravity waves may be derived by linearising (1)-(3) about the equilibrium point (8). In the absence of velocity shear, this yields (k g N/k) 2 = ω 2 + iωk 2 (ν + χ) νχk 4, (23) where k g = k ê g and N is the Brunt-Vaisala frequency (18), which in the incompressible limit (γ ) reduces to N 2 = g ln ρ = α ρ g T = α ρ g z T. The system becomes (R-T) unstable when ω 2 < 0 or (k g N/k) 2 < k 4 (ν 2 + χ 2 2νχ). (24)

29 in incompressible, ideal fluids Since the right hand side is always positive, the diffusive terms have a stabilising effect, predominantly at small scales (large k), allowing an incompressible fluid to support a finite temperature gradient, z T, and a compressible fluid to support a super-adiabatic gradient, in excess of (19). In an incompressible ( V = 0, γ ), non-dissipative (ν = 0 = χ), fluid, (23) and (24) become (ωk/k g N) 2 = 0 and ê g T > 0, or z T < 0, respectively. The relationship between waves and instabilities, as outlined above, is common to both fluid and plasma dynamics. Magnetized plasma equilibria are thus subject to a range of instabilities, most of which can be associated with corresponding plasma waves.

30 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Three types of MHD instabilities The underlying reason for fluid, and specifically magnetohydrodynamic, instabilities has been succinctly summarized by Biskamp (2003), A fluid, in particular a plasma, becomes unstable when the gradient of velocity, pressure, or magnetic field exceeds a certain threshold, which occurs, roughly speaking, when the convective transport of momentum, heat, or magnetic flux is more efficient than the corresponding diffusive transport by viscosity, thermal conduction, or resistivity. There are hence three types of instabilities, which play a fundamental role in macroscopic plasma dynamics: the Kelvin-Helmholtz instability driven by a velocity shear; the Rayleigh-Taylor instability caused by the buoyancy force in a stratified system; and current-driven MHD instabilities in a magnetized plasma.

31 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Linearized MHD equations Before investigating MHD instabilities, all of which are related to plasma non-uniformity, it is worthwhile to review the basic MHD waves in a uniform magnetized plasma. For this purpose we subject an MHD equilibrium to small perturbations of the form exp(ik x iωt) (denoting equilibrium and perturbed quantities by subscripts 0 and 1, respectively), which is equivalent to Fourier transforming the MHD equations in space and time and neglecting products of perturbed quantities, such as V 1 B 1, ωρ 1 = ρ 0 k V 1, E 1 + V 1 B 0 = 0, (25) ωρ 0 V 1 = kp 1 + B 0 J 1 /i ik B 1 = µ 0 J 1, (26) ωp 1 = γp 0 k V 1, k E 1 ωb 1 = 0. (27)

32 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Perturbed velocity eq n Note: Fourier transform reduces the differential operators to algebraic ones: t iω and ik. Since equilibrium is assumed to be stationary, V 0 = 0, then d t = t + V iω + iv 1 k. In the above, we also neglected the source terms and replaced the ideal gas ratio of specific heats 5 3 by some general value γ. Note that B = 0 implies k B 1 = 0, which also follows from (27). After simplification, (25)-(27) reduce to a single equation for the perturbed velocity, V 1, [ ] ω 2 B0 ρ 0 V 1 = (k B 0 ) + γp 0 k k V 1 k B0 (k V 1 ) B 0. µ 0 µ 0 (28)

33 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation MHD dispersion relation in a uniform plasma This velocity can have both longitudinal (k V 1 ) and transverse (k V 1 ) components, thus allowing both compressional and shear deformations. We next choose a Cartesian coordinate system (x, y, z) in which the the z-axis is aligned with the magnetic field (ê z = b 0 = B 0 /B). Without loss of generality, we can assume the wave to propagate in the in the y z plane, so that k = k ê z + k ê y. The vector equation (28) then yields three algebraic equations for the velocity components V 1x, V 1y, V 1z. The eigenmodes of this system may be found by requiring the determinant of the matrix to vanish, which yields the following dispersion relation, [ ω 2 (k V A ) 2] [ ω 4 (ωkv max ) 2 + (kk V S V A ) 2] = 0, (29)

34 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Roots of the MHD dispersion relation where V S is the sound speed (11), V A is the Alfvén speed, Vmax 2 = VS 2 + V A 2 and k2 = k 2 + k2. The three roots of this equation, namely, ωa 2 = k2 V A 2, ω2 ± = 1 2 k2 Vmax (1 2 ± ) 1 α 2, (30) α 2 = 4(k /k) 2 (V A V S /V 2 max) 2, (31) represent the three possible MHD waves in a uniform plasma. Since the discriminant in (30) is positive definite, all three waves are stable. This follows because k /k 1 and (a 2 + b 2 ) 2 /a 2 b 2 = 2 + a 2 /b 2 + b 2 /a 2 4, so that α 2 1.

35 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Shear Alfvén wave The first root (ω A ) corresponds to the shear Alfvén wave which travels along the field lines (k b) at a phase velocity equal to the Alfvén speed, ω A /k = ±V A The Alfvén waves are non-dispersive, with same group velocity v gr = k ω A = ±V A for all k. It consists of purely transverse, and hence incompressible (V 1 k = 0), plasma motions and magnetic field displacements (B 1 V 1 ). The origin of the shear Alfvén wave is the combination of ion inertia and field line tension, T B, as discussed in Lecture III.

36 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Magnetosonic waves The second (ω + ) and third (ω ) roots of (30) represent the fast and slow magnetosonic waves, which as the name suggests are a combination of Alfvén and sound waves. They can propagate both along and across B 0, k b 0 ω 2 + = k 2 V 2 max, ω 2 = 0, (32) k b < ω2 +/k 2 V 2 max < 1, 0 < ω 2 /k 2 V 2 max < 1 2.(33) The fast wave includes in-phase oscillations of plasma pressure and magnetic stress, p + T B, with T B being the Maxwell stress tensor, while the slow wave includes out-of-phase oscillations; the perturbations to p and T B are compressional, while those to T B are tensional, i.e., involve transverse shearing.

37 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Low and high beta MHD waves Recall that T B > 0 is indicative of magnetic pressure, which counteracts compressional deformations, while T B < 0 represents magnetic tension, which counteracts bending and/or twisting deformations. Hence, the MHD momentum equation can be written as ρd t V + (p + T B ) = 0, see Lecture III. In the limits of small and large plasma pressure, compared to the magnetic pressure, i.e., when β (V S /V A ) 2 1 or 1, the fast wave reduces to the Alfvén and acoustic waves, respectively, while the slow wave vanishes entirely, β 1 ω 2 + k 2 V 2 A, ω2 0, (34) β 1 ω 2 + k 2 V 2 S, ω2 0. (35) The low β perpendicular wave (32), involves purely longitudinal deformations and is known as the compressional Alfvén wave.

38 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation R-T instability in a magnetized fluid We next consider the effect of finite B 0 on the stability of a plasma in a gravitiation field, i.e., on the MHD version of the R-T and K-H instabilities. In ideal MHD (ν = χ = η = 0), (23) becomes (k g N/k) 2 + k 2 V 2 A = ω2, k g ê g (k ê g ) (36) which shows that magnetic field has a stabilising effect on the R-T instability provided that k = k b does not vanish. The effect is largest for k = k, when vertical motions (V 1 ê g ) involve bending deformations of the magnetic field lines.

39 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Flute-like perturbations most unstable Such deformations excite shear Alfvén waves, which transport the free energy released by the bouyancy force away from the point of disturbance. For k g = k = k, the phase velocity becomes ω 2 /k 2 = (N/k) 2 + VA 2 and the waves are R-T unstable when N 2 > (kv A ) 2. As a result, the R-T instability in a magnetized plasma proceeds by internal gravity waves perpendicular to both g and b, i.e., by B-field aligned (k g b, k 0, k = k ), or flute-like, perturbations for which the stabilising effect of magnetic tension T B, which plays a role analogous to the surface tension T S, is minimized. The name refers to the flutes in a grecian column which are resembled by field aligned perturbations in a cylindrical plasma column.

40 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation K-H instability in a magnetized fluid Similarly, neglecting surface tension and assuming ρ a = ρ b for simplicity, one finds the dispersion relation for the K-H instability across a vorticity sheet in a non-dissipative magnetized plasma, ω k = 1 2 (V b + V a ) ± which is the ideal MHD version of (20). V 2 A 1 4 (V b V a ) 2, (37) Once again magnetic field has a stabilising effect, with the K-H instability suppressed entirely when V A > V a V b /2, i.e., when the kinetic energy of the flow discontinuity can be transported away from the disturbance by the shear Alfvén wave.

41 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Linearized MHD force balance We are now ready to consider the stability properties of a magnetically confined plasma equilibrium, i.e., a non-uniform, toroidal plasma supporting a radial pressure gradient. We begin with the MHD momentum equation, recall Lecture II, replacing p by the divergence of the pressure tensor, p = p + π, f ρd t V + π = J B p, d t = t + V. (38) Here we denoted the left hand side, which represents the plasma inertia and/or the net force acting on a fluid element, by f; note that f = 0 in equilibrium, i.e., in the absence of flow, V = 0 = π. Let us consider f as an operator acting on a small displacement away from equilibrium ξ(x), f(ξ) = J 1 B 0 + J 0 B 1 p 1. (39)

42 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Linearized MHD force operator, f Here we linearized f, neglecting products of small quantities; the perturbed pressure, p 1, current, J 1, and field, B 1, are found by linearising the MHD fluid equations and integrating over time (for p 1 and B 1 only), µ 0 J 1 = B 1, B 1 = (ξ B 0 ), p 1 = γp 0 ξ ξ p 0. (40) Inserting (39) into (38) gives the explicit form of the linear force operator, f(ξ) = µ 1 0 ( B 0) B 1 +µ 1 0 ( B 1) B 0 + (γp 0 ξ+ξ p 0 ). (41)

43 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Properties of f This operator has several important properties: 1 it is self-adjoint, η f(ξ)dx = for any two displacements ξ and η; ξ f(η)dx, (42) 2 the squares of its eigenvalues ω 2 are purely real for all discrete normal modes ξ(x, t) = ξ ω (x) exp( iωt), such that its eigenvalues ω are either real or purely imaginary, 3 as a result, all discrete normal modes are orthogonal, (ωn 2 ωm) 2 ρξm ξ n dx = 0; (43)

44 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Properties of f 4 f(ξ) permits both discrete and continuum eigenvalues, although the latter are only allowed in the stable domain, i.e., for ω 2 0, and represent wave continua; 5 consequently, the transition from unstable to stable domains must occur via the marginally stable state, defined by ω = 0 = f(ξ ω ). The above properties suggest that normal mode analysis is well suited for investigating MHD spectral stability of confined plasmas, and indeed, it is the most common technique used for this purpose.

45 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation MHD energy principle I An alternative approach relies on variational techniques and the energy principle, which must be satisfied for MHD linear stability, namely δ 2 W [ξ(x)] > 0, (44) where δ 2 W is the second variation of the potential (free) energy; the first variation δw vanishes in equilibrium. This variation is equal to the work done in displacing the plasma element against the force f, given by (39), δ 2 W [ξ] = 1 ξ f(ξ)dx. (45) 2 The relation between the two approaches is precisely the same as that between the Schrödinger and Heisenberg pictures of quantum mechanics

46 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation MHD energy principle II Inserting (41) into (45) and simplifying yields the following expression for δ 2 W [ξ], 1 2 p γp 0 ξ 2 + (ξ p 0 )( ξ ) J 0 (B 1 ξ ) + B2 1 µ 0 dx s (p 1 + B 0 B1 µ 0 )ξ ds + 1 B1 2 2 v µ 0 dx, (46) where the subscripts p, v and s denote integrals over the plasma volume, δ 2 W p, the vacuum envelope, which is assumed to surround the plasma, δ 2 W v, and the plasma-vacuum surface, δ 2 W s, respectively. The surface integral, which follows from application of Gauss s theorem, vanishes for purely tangential displacements, e.g., for a perfectly conducting boundary.

47 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Internal and external MHD instabilities This suggests a natural division of MHD instabilities into internal, or fixed boundary, modes, which do not deform the plasma-vacuum surface, and external, or free boundary, modes, which include such deformations. The fixed boundary form of the MHD energy principle may be obtained by neglecting the surface and vacuum terms, separating ξ, B 1 and J 1 into parallel and perpendicular components, and effecting several cancelations involving ξ,

48 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation MHD energy principle III δ 2 W p [ξ] = 1 2 p γp 0 ξ 2 2(ξ p 0 )(ξ κ) J 0 (B 1 ξ ) + B2 1 µ 0 + B2 0 µ 0 ξ + 2ξ κ 2 dx > 0. (47) Let us consider the energy associated with each of the terms above: 1 fluid compression acoustic waves, 2 flux tube interchange R-T instability, 3 flux tube twisting kink instability; the name originates in the kinking (buckling or twisting into a spiral) of a current filament, 4 magnetic tension, or flux tube bending shear Alfvén waves, 5 magnetic compression compressional Alfvén waves.

49 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation MHD energy principle IV Only the kink and interchange terms above can lower the plasma free energy, thus destabilizing the equilibrium. These two terms are known collectively as the Newcomb terms. The interchange term typically changes sign as one travels poloidally around the flux surface, i.e., it is positive (destabilizing) on the outer side of the torus and negative (stabilizing) on the inner side. However, global stability is determined by the volume, i.e., flux surface, average of all the terms in (47). In contrast, the compressional and tensile terms, which perturb the total stress tensor, p + T B, always increase this energy thus providing a stabilising effect.

50 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation General plasma vorticity eq n As a result, the most stringent stability criteria are imposed by those perturbations which satisfy (p 1 + T B1 ) = 0. In other words, an equilibrium is MHD stable, if an only if, it is stable against incompressible displacements. Let us now return to (38) to derive the primary dynamical relation for plasma fluid instabilities. Taking a parallel projection of the curl of (38) yields the general plasma vorticity equation, b ( f 2κ f) = b [B 2 (J /B) + 2κ p ], (48) where f now denotes the fluid inertia, i.e., the left side of (38).

51 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Shear Alfvén law This relation, also known as the shear Alfvén law, is of fundamental importance in the theory of confined plasma stability. It is worth noting that (J /B) and p in (48) originate in the gradients of parallel and perpendicular currents. To illustrate this fact, (48) may be recast term by term into a charge conservation equation, J = 0, with J = J + J + J p, or J p = J J, (49) expressing the total current continuity (quasi-neutrality) as the sum of polarisation, parallel return and diamagnetic currents, J p = (ρ/b)b ( t + V E )V E, J = b p/b. (50)

52 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Implications of quasi-neutrality Here it is useful to recall the discussion of J = 0 in Lecture III: In steady state (d t = 0), we recover the equilibrium relation, (J + J ) = 0, whereas transiently the net compression of J + J gives rise to the polarisation current, and hence to perpendicular plasma motion. Depending on the ordering of v E, this motion can either be fast (v E /v ti 1), as in the MHD ordering assumed here, or slow (v E /v ti δ i ), as in the drift ordering. While global equilibrium can be maintained by a parallel return current, polarisation currents are always generated on small scales. They are the origin of so-called micro-instabilities (e.g., drift waves, interchange modes, etc.) which are responsible for plasma turbulence.

53 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Bending/twisting and interchange perturbations Due to the ever-present interchange forcing, micro-instabilities exhibit a ballooning character, i.e., are most active in regions of unfavourable magnetic curvature. In short, the two terms on the right hand side of (48), b ( f 2κ f) = b [B 2 (J /B) + 2κ p ], whose energies are given by the third and fourth and second terms in (47), represent vorticity sources due to gradients of parallel current and scalar pressure (or diamagnetic current) and give rise to instabilities involving flux tube bending/twisting and interchange, respectively.

54 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Pressure gradient driven instabilities I Let us first consider the pressure gradient term in (48). This term, which can also be written as 2b κ p, vanishes when κ p = 0, i.e., when κ p, and is largest when (b κ) p. It vanishes for the equilibrium pressure gradient p 0 at the inboard and outboard mid-planes of the torus, since there κ p 0 R, but generates a vorticity source of 2κ n k p 1 due to diamagnetic (poloidal) pressure variation. Here we assume an eikonal of S = k x, and decompose k into diamagnetic and radial components, ξ exp(ik x), k = k b + k ê = k b + k ê + k ê. (51)

55 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Pressure gradient driven instabilities II Substituting for p 1 with p 1 ξ p 0, (40), we see that such poloidal variation is caused by flute-like radial displacements, generating a vorticity source of 2k ξ κ n p 0. We thus expect the condition for interchange instability to be related to the sign of the product κ n p 0 or, more generally, of κ p 0. This is confirmed by the second term in (47), which indicates that ξ lowers the potential energy of the plasma, δ 2 W p [ξ], and hence drives the interchange instability, when (ξ κ)(ξ p 0 ) > 0.

56 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Pressure gradient driven instabilities III This corresponds to the outboard side of the torus, where κ p 0 > 0, which defines the region of unfavourable, or bad, magnetic curvature. Similarly, κ is said to be favourable, or good, when κ p 0 < 0, which occurs on the inboard side of the torus. Based on this preliminary analysis, we may conclude that the magnetic curvature κ plays the role of negative gravity, g, and subjects the plasma fluid to a buoyancy force proportional to κ p 0.

57 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Pressure gradient driven instabilities IV By expressing the grad-b and curvature GC drifts in gravitational form, recall Lecture II, this negative gravity is seen to represent the centripetal acceleration of a charged particle (fluid element) moving along a curved magnetic field line (flux tube). The buoyancy force is easily understood in the particle picture: following a flute-like perturbation ξ exp(ik θ) of the flux surface r = const, the curvature GC drift, v κ e s κ b e s ê, which has opposite sign for ions and electrons, leads to a poloidal electric field, and thus to a radial (outward) E B drift.

58 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Pressure gradient driven instabilities V It is worth stressing that the plasma remains quasi-neutral at all times. In effect, the radial plasma motion is required to prevent the separation (polarisation) of charge that would otherwise ensue due to an unchecked curvature drift i.e., v κ (J ) leads to V E (J p ) necessary to ensure that (50) is satisfied. Hence, the frequently cited explanation of the radial plasma motion as being due to charge polarisation is not, strictly speaking, correct.

59 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Current driven instabilities Let us next examine the current gradient term (48), which, when linearized, can be decomposed into two parts, B 2 b (J /B) B 0 B 1 (J 0 /B 0 ) + B 0 B 0 (J 1 /B 0 ), (52) corresponding to the Newcomb terms in (47). Here, the first part, which contains gradients of the equilibrium parallel current, corresponds to term (iii) in (47), associated with twisting deformations and the kink instability. Similarly, the second part, which contains gradients of the perturbed parallel current, corresponds to term (iv) in (47), associated with bending deformations (note that µ 0 J 1 = B 1 is largest when b 1 b 0 ) and the familiar stabilising effect of magnetic tension.

60 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Magnetic differential eq n Hence, the displacement ξ(x) (or more accurately, the normal mode amplitude ξ ω (x), as discussed earlier) is most likely to grow when it exhibits a flute-like character, b 0 = 0 0, b 0 = B 0 /B 0, (53) so that field line bending is minimized. The inhomogenous form of (53), namely f = A/B, is known as the magnetic differential equation and has a solution if, and only if, it satisfies the Newcomb solubility conditions, f = A/B (A/B)dl = 0 = (A/B θ ) mn. (54)

61 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Spectral MHD stability: Fourier expansion To analyse the spectral MHD stability of the plasma, we express ξ(x) in flux co-ordinates (r f, θ f, ζ f ), and Fourier expand ξ(r f, θ f, ζ f ) in terms of the poloidal and toroidal harmonics with mode numbers m and n, and Fourier components ξ mn (r f ), which we assume to be flux surface labels, ξ(r f, θ f, ζ f ) = m,n ξ mn (r f ) exp(imθ f inζ f ). (55) Combining (55) and (53) yields a criterion for flute-like perturbations, B θ ( θf +q ζf )ξ 0, ξ mn (r f ) 0, m/n q, (56) which states that ξ mn (r f ) must vanish on all ergodic flux surfaces,

62 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Resonant flux surfaces most unstable so that ξ becomes localized to resonant flux surfaces on which ξ mn (r f ) are field line labels. Since k = (m qn) θ f, (56) may also be written as k (m qn)ψ p/ gb k k, k 0 (57) which is clearly satisfied near resonant surfaces, q m/n. Note that k k implies that the frequency of the magnetosonic waves, ω + (32), greatly exceeds that of the shear Alfvén waves, ω A (30), such that the former equilibrate on the Alfvén time scales considered, and can thus be neglected in MHD stability analysis. This equilibration establishes a quasi-static perpendicular force balance which leads to a purely electrostatic perpendicular electric field, E = ϕ; this relation is employed in both flute reduced-mhd and in DHD turbulence, see Lecture VI.

63 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation MHD instabilities in a general screw pinch I* Although the above equations contain all the essential elements of MHD stability, they are too complicated for analytic treatment in all but a few simple systems. Below we examine two such simple systems, the general screw pinch and the cylindrical tokamak, both of which neglect the effect of toroidal curvature, so that κ = κ r (r)ê = b 2 Pê /r is purely poloidal, and hence a flux surface label. Likewise, both systems exhibit θ and z (or ζ) symmetry, so that the displacement ξ is easily Fourier transformed in θ and z, and requires only a single radial function, ξ(x) = ξ mn (r) exp(imθ inζ) = ξ mn (r) exp(ik P rθ ik T z).

64 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation MHD instabilities in a general screw pinch II* As a result of these symmetries, the Fourier components ξ mn (r) in (55) become decoupled, i.e., evolve separately, allowing the energy variation (46) and the full eigenmode equation for ξ mn (r) to be evaluated explicitly In short, the model retains the essential elements of both kink and interchange modes while dispensing with most of the algebraic complexities related to toroidal geometry. It is thus very useful for introducing both instabilities.

65 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation MHD instabilities in a general screw pinch III* Writing the displacement ξ(x) as ξ b + ξ ê + ξ ê, where b = b P ê θ + b T ê z, ê = b T ê θ b P ê z, (58) ξ = ξ P b P + ξ T b T, ξ = ξ P b T ξ T b P, (59) the plasma energy variation may be calculated as (Freidberg, 1991), δ 2 [ W p [ξ] a k 2π 2 = f ξ 2 R 0 /µ + gξ2 dr + k ] 0 k 2 Ba 2 ξ a 2, = d r, (60) 0 where k and k are the parallel and diamagnetic components of k, a

66 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation MHD instabilities in a general screw pinch IV* k = k T b T + k P b P, k = k T b P k P b T, (61) k = k T b T k P b P, k 2 = k 2 T + k2 P, (62) k T and k P are its axial (toroidal) and poloidal components, k T = n/r 0, k P = m/r, (63) and f (r) and g(r) are functions of radial position given by ( ) k B 2 ( f (r) = r = rk2 T B2 T k k 2 1 m ) 2 = B2 P qn rk 2 (qn m)2,(64) ( ) g(r) = 2 k T k µ 0p 0 + k 2 r k2 T k f (r) k 2 k r 2. (65)

67 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation MHD instabilities in a general screw pinch V* The vacuum energy variation is found as [ δ 2 W v [ξ] Λw a 2 k 2 ] 2π 2 = Ba 2 ξ a 2 (66) R 0 /µ 0 m where Λ w is a geometrical function of the minor radius, a, and the conducting wall radius, b, involving modified Bessel functions taking the argument k T a and k T b. For k T b 1, this function may be estimated as Λ w 1 + (a/b)2 m 1 (a/b) 2 m, k T b 1, a < b, (67) which for a/b 1 reduces to 1 + 2(a/b) 2 m 1, while for a/b 1 becomes singular as [(1 a/b) m ] 1 ; in this limit, a finite ξ a leads to an infinite increase of δ 2 W v and is thus prohibited, i.e., ξ (a) 0. a

68 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation MHD instabilities in a general screw pinch VI* Imposing the condition that ξ (r) minimise δ 2 W p (60), where we neglect the vacuum energy δ 2 W v (66) for simplicity, allows us to derive the corresponding eigenmode equation for ξ (r). Setting the variation of (60) with respect to ξ (r) to zero, yields an Euler-Lagrange equation, [f (r)ξ ] g(r)ξ = 0, (68) with f (r) and g(r) defined above. This MHD dispersion relation is the cylindrical version of (29), and describes all MHD waves and instabilities in a marginally stable screw pinch equilibrium.

69 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation MHD instabilities in a general screw pinch VII* We next write down the corresponding full eigenmode problem for ξ(x), f[ξ(x)] + ρω 2 ξ(x) = 0 (69) where f[ξ(x)] is given by (41). For a general screw pinch this problem takes the form of a second order ODE for the radial displacement ξ, [F(r)(rξ ) ] G(r)rξ = 0, (70) where F(r) and G(r) are functions of frequency, wave number and radius, F(r) = ρv max 2 (ω 2 ωa 2 )(ω2 ωh 2 ) r (ω 2 ω+ 2 )(ω2 ω 2 ) (71)

70 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation MHD instabilities in a general screw pinch VIII* G(r) = ρ(ω2 ωa 2 ) ( 4k 2 + T V 2 ) A B2 P (ω 2 ωs 2) r µ 0 r 3 (ω 2 ω+ 2 )(ω2 ω 2 [ ) B 2 ( ) + P µ 0 r 2 2kT GB P V 2 max (ω 2 ωh 2 ) ] µ 0 r 2 (ω 2 ω+ 2 )(ω2 ω 2 ). (72) Here, the Alfvénic and magneto-sonic frequencies, ω A and ω ±, defined in (30), are supplemented by the sonic and hybrid frequencies, ω S and ω H, ω 2 S = k2 V 2 S, ω2 H = k2 V 2 H, V 2 H = V 2 A V 2 S /V 2 max, (73) where V A and V S are the Alfven and sound speeds, see Lecture III, and V 2 max = V 2 A + V 2 S ; finally, G = k PB T k T B P.

71 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation MHD instabilities in a general screw pinch IX* The eigenmode equation (70) is a more general version of the marginally stable dispersion relation (68), and describes all possible in a general screw pinch. Indeed, it reduces to the latter result when the marginal stability condition (ω = 0) is imposed on (71)-(72). The agreement between the two results confirms our earlier assertion that a marginally stable eigenmode of the MHD force operator f[ξ(x)] (41) must minimise the free energy of the plasma; this is a direct consequence of the self-adjointness of f, see (42).

72 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Flute-reduced MHD I* The above results indicate that toroidal plasmas are most unstable to flute-like perturbations ξ exp(ik x) which satisfy (53). These are in general radially localized to resonant flux surfaces, although in large aspect ratio (small ɛ) tokamaks, where all displacements satisfy (53), b 0 k /k r/r 0 q = ɛ/q O(ɛ) 1. (74) they also apply on ergodic flux surfaces. To show this, we first note that in cylindarical geometry, the poloidal and axial components of k are k P = m/r and k T = n/r 0, while parallel and diamagnetic components become k R 0 = m/q n, k r m, (75) so that the ratio k /k O(ɛ/q) is small everywhere when ɛ 1 and q O(1).

73 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Flute-reduced MHD II* Here we assume that the perturbation is localized to the flux surface r = const and thus exclude the radial component of k. Hence k = b (k b) reduces to the projection in the diamagnetic direction k = k ê. More generally, k should be replaced by k = k. Of course, k 1 is not true for all perturbations when ɛ 1, although it is always satisfied close to resonant flux surfaces, q m/n. Hence, it is useful to derive a version of the MHD (and more generally of DHD) model in which flute-like displacements are assumed from the outset.

74 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Flute-reduced MHD III* Such a model, known a flute-reduced MHD may be derived by imposing the flute-ordering (74) on the MHD equations. This is accomplished by applying a multi-scale expansion in the small parameter ɛ qk /k 1, to the MHD, or more generally to DHD, equations and retaining only leading order terms. The resulting equations are a useful reference point for the study of stability and dynamics of toroidal plasmas. Below we outline the derivation, broadly following the original account of Strauss (1976), yet making explicit the multiple-scale expansion.

75 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Flute-reduced MHD IV* Multi-scale expansion involves separation of slow, O(1), and fast, O(ɛ 1 ), spatial variations into (independent) co-ordinates X and x, and slow, O(ɛ 2 ), and fast, O(1), temporal variations into co-ordinates T and t, x ɛ 1 x + X, t t + ɛ 2 T. (76) Hence, the net gradient can be expressed as the sum of slow, 0 = X, and fast, 1 = x = + O(ɛ), gradients, and the net time derivative as the sum of slow, T, and fast, t, time derivatives, ɛ 1 + 0, t t + ɛ 2 T. (77) Similarly, any quantity A is expanded in powers of ɛ into an equilibrium part A 0 and the flute-perturbed parts A 1, A 2,...,

76 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Flute-reduced MHD V* satisfying (74), as follows A = A 0 (X, T ) + ɛ n A n (x, X, t, T ). (78) n=1 Hence the fast variables x and t measure the variation of the flute-perturbed quantities, and the slow variables X and T that of equilibrium quantities. First, we assume the existence of an MHD equilibrium for B 0, J 0 and p 0, J 0 B 0 = 0 p 0, 0 B 0 = 0, µ 0 J 0 = 0 B 0. (79) We next note that in a tokamak, parallel projection of Ohm s law expresses the balance between inductive EMF and collisional resistivity, so that the equilibrium electric field is purely parallel,

77 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Flute-reduced MHD VI* t A 0 = E A 0 = E 0 = η J 0, E 0 = ϕ 1 t A 1 = 0, (80) and the first order potentials must vanish. Flute-reduction then allows us to express the electromagnetic (vector ) field in terms of two scalar fields ψ 1 and ϕ 2. Due to the flute-ordering, (74), the perturbed electric field is mainly transverse and electrostatic E 1 E 1 ϕ 2. (81) The perturbed magnetic field, which follows from B = 0 and (74), implies B 1 0 to leading order. Hence, B 1 may be written as B 1 = A 2 = (A 2 +b 0 A 2 ) = B 1 b 0 + ψ 1 b 0, (82)

78 in a uniform plasma in a stratified plasma in a magnetically confined plasma in a general screw pinch in the flute-reduced approximation Flute-reduced MHD VII* where ψ 1 = A 2 denotes the parallel vector potential. As a result, the parallel gradient, (B/B 0 ), can be simplified as = B B 1 = b 0 0 B 1 0 {ψ 1, } b 0 0 { ψ 1, }, (83) where breve denotes division by B 0 and {A, } is the Poisson bracket, Ă = A/B 0, {A, } b 0 A = b 0 A. (84) The first term in the parallel gradient represents the variation along the equilibrium field, while the second term represents variation due to field line bending and perpendicular gradients. By nature of (74), only transverse deflection B 1 = ψ 1 b 0 contributes to the latter.

Dispersive Media, Lecture 7 - Thomas Johnson 1. Waves in plasmas. T. Johnson

Dispersive Media, Lecture 7 - Thomas Johnson 1. Waves in plasmas. T. Johnson 2017-02-14 Dispersive Media, Lecture 7 - Thomas Johnson 1 Waves in plasmas T. Johnson Introduction to plasmas as a coupled system Magneto-Hydro Dynamics, MHD Plasmas without magnetic fields Cold plasmas

More information

2/8/16 Dispersive Media, Lecture 5 - Thomas Johnson 1. Waves in plasmas. T. Johnson

2/8/16 Dispersive Media, Lecture 5 - Thomas Johnson 1. Waves in plasmas. T. Johnson 2/8/16 Dispersive Media, Lecture 5 - Thomas Johnson 1 Waves in plasmas T. Johnson Introduction to plasma physics Magneto-Hydro Dynamics, MHD Plasmas without magnetic fields Cold plasmas Transverse waves

More information

Macroscopic plasma description

Macroscopic plasma description Macroscopic plasma description Macroscopic plasma theories are fluid theories at different levels single fluid (magnetohydrodynamics MHD) two-fluid (multifluid, separate equations for electron and ion

More information

Recapitulation: Questions on Chaps. 1 and 2 #A

Recapitulation: Questions on Chaps. 1 and 2 #A Recapitulation: Questions on Chaps. 1 and 2 #A Chapter 1. Introduction What is the importance of plasma physics? How are plasmas confined in the laboratory and in nature? Why are plasmas important in astrophysics?

More information

Fluid equations, magnetohydrodynamics

Fluid equations, magnetohydrodynamics Fluid equations, magnetohydrodynamics Multi-fluid theory Equation of state Single-fluid theory Generalised Ohm s law Magnetic tension and plasma beta Stationarity and equilibria Validity of magnetohydrodynamics

More information

Magnetohydrodynamic waves in a plasma

Magnetohydrodynamic waves in a plasma Department of Physics Seminar 1b Magnetohydrodynamic waves in a plasma Author: Janez Kokalj Advisor: prof. dr. Tomaž Gyergyek Petelinje, April 2016 Abstract Plasma can sustain different wave phenomena.

More information

PHYS 432 Physics of Fluids: Instabilities

PHYS 432 Physics of Fluids: Instabilities PHYS 432 Physics of Fluids: Instabilities 1. Internal gravity waves Background state being perturbed: A stratified fluid in hydrostatic balance. It can be constant density like the ocean or compressible

More information

Lecture # 3. Introduction to Kink Modes the Kruskal- Shafranov Limit.

Lecture # 3. Introduction to Kink Modes the Kruskal- Shafranov Limit. Lecture # 3. Introduction to Kink Modes the Kruskal- Shafranov Limit. Steve Cowley UCLA. This lecture is meant to introduce the simplest ideas about kink modes. It would take many lectures to develop the

More information

Prototype Instabilities

Prototype Instabilities Prototype Instabilities David Randall Introduction Broadly speaking, a growing atmospheric disturbance can draw its kinetic energy from two possible sources: the kinetic and available potential energies

More information

Simple examples of MHD equilibria

Simple examples of MHD equilibria Department of Physics Seminar. grade: Nuclear engineering Simple examples of MHD equilibria Author: Ingrid Vavtar Mentor: prof. ddr. Tomaž Gyergyek Ljubljana, 017 Summary: In this seminar paper I will

More information

MHD Linear Stability Analysis Using a Full Wave Code

MHD Linear Stability Analysis Using a Full Wave Code US-Japan JIFT Workshop on Progress of Extended MHD Models NIFS, Toki,Japan 2007/03/27 MHD Linear Stability Analysis Using a Full Wave Code T. Akutsu and A. Fukuyama Department of Nuclear Engineering, Kyoto

More information

The Virial Theorem, MHD Equilibria, and Force-Free Fields

The Virial Theorem, MHD Equilibria, and Force-Free Fields The Virial Theorem, MHD Equilibria, and Force-Free Fields Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 10 12, 2014 These lecture notes are largely

More information

Current-driven instabilities

Current-driven instabilities Current-driven instabilities Ben Dudson Department of Physics, University of York, Heslington, York YO10 5DD, UK 21 st February 2014 Ben Dudson Magnetic Confinement Fusion (1 of 23) Previously In the last

More information

Heating and current drive: Radio Frequency

Heating and current drive: Radio Frequency Heating and current drive: Radio Frequency Dr Ben Dudson Department of Physics, University of York Heslington, York YO10 5DD, UK 13 th February 2012 Dr Ben Dudson Magnetic Confinement Fusion (1 of 26)

More information

Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas )

Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas ) Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas ) Yasutomo ISHII and Andrei SMOLYAKOV 1) Japan Atomic Energy Agency, Ibaraki 311-0102, Japan 1) University

More information

1 Energy dissipation in astrophysical plasmas

1 Energy dissipation in astrophysical plasmas 1 1 Energy dissipation in astrophysical plasmas The following presentation should give a summary of possible mechanisms, that can give rise to temperatures in astrophysical plasmas. It will be classified

More information

Contents. I Introduction 1. Preface. xiii

Contents. I Introduction 1. Preface. xiii Contents Preface xiii I Introduction 1 1 Continuous matter 3 1.1 Molecules................................ 4 1.2 The continuum approximation.................... 6 1.3 Newtonian mechanics.........................

More information

CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION

CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,

More information

Jet Stability: A computational survey

Jet Stability: A computational survey Jet Stability Galway 2008-1 Jet Stability: A computational survey Rony Keppens Centre for Plasma-Astrophysics, K.U.Leuven (Belgium) & FOM-Institute for Plasma Physics Rijnhuizen & Astronomical Institute,

More information

Prof. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit

Prof. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit Prof. dr. A. Achterberg, Astronomical Dept., IMAPP, Radboud Universiteit Rough breakdown of MHD shocks Jump conditions: flux in = flux out mass flux: ρv n magnetic flux: B n Normal momentum flux: ρv n

More information

SOLAR MHD Lecture 2 Plan

SOLAR MHD Lecture 2 Plan SOLAR MHD Lecture Plan Magnetostatic Equilibrium ü Structure of Magnetic Flux Tubes ü Force-free fields Waves in a homogenous magnetized medium ü Linearized wave equation ü Alfvén wave ü Magnetoacoustic

More information

Magnetohydrodynamic Waves

Magnetohydrodynamic Waves Magnetohydrodynamic Waves Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 17, 2016 These slides are largely based off of 4.5 and 4.8 of The Physics of

More information

Solar Physics & Space Plasma Research Center (SP 2 RC) MHD Waves

Solar Physics & Space Plasma Research Center (SP 2 RC) MHD Waves MHD Waves Robertus vfs Robertus@sheffield.ac.uk SP RC, School of Mathematics & Statistics, The (UK) What are MHD waves? How do we communicate in MHD? MHD is kind! MHD waves are propagating perturbations

More information

196 7 atmospheric oscillations:

196 7 atmospheric oscillations: 196 7 atmospheric oscillations: 7.4 INTERNAL GRAVITY (BUOYANCY) WAVES We now consider the nature of gravity wave propagation in the atmosphere. Atmospheric gravity waves can only exist when the atmosphere

More information

9 Fluid Instabilities

9 Fluid Instabilities 9. Stability of a shear flow In many situations, gaseous flows can be subject to fluid instabilities in which small perturbations can rapidly flow, thereby tapping a source of free energy. An important

More information

The Euler Equation of Gas-Dynamics

The Euler Equation of Gas-Dynamics The Euler Equation of Gas-Dynamics A. Mignone October 24, 217 In this lecture we study some properties of the Euler equations of gasdynamics, + (u) = ( ) u + u u + p = a p + u p + γp u = where, p and u

More information

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t) IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common

More information

Chapter 1. Governing Equations of GFD. 1.1 Mass continuity

Chapter 1. Governing Equations of GFD. 1.1 Mass continuity Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for

More information

Ideal MHD Equilibria

Ideal MHD Equilibria CapSel Equil - 01 Ideal MHD Equilibria keppens@rijnh.nl steady state ( t = 0) smoothly varying solutions to MHD equations solutions without discontinuities conservative or non-conservative formulation

More information

MAGNETOHYDRODYNAMIC EQUILIBRIUM AND STABILITY OF PLASMA

MAGNETOHYDRODYNAMIC EQUILIBRIUM AND STABILITY OF PLASMA University of Ljubljana Faculty of Mathematics and Physics Seminar 1 b -1st year, II. cycle MAGNETOHYDRODYNAMIC EQUILIBRIUM AND STABILITY OF PLASMA Author: Lino alamon Advisor: prof. dr. Tomaº Gyergyek

More information

Magnetohydrodynamics Stability of a Compressible Fluid Layer Below a Vacuum Medium

Magnetohydrodynamics Stability of a Compressible Fluid Layer Below a Vacuum Medium Mechanics and Mechanical Engineering Vol. 12, No. 3 (2008) 267 274 c Technical University of Lodz Magnetohydrodynamics Stability of a Compressible Fluid Layer Below a Vacuum Medium Emad E. Elmahdy Mathematics

More information

TURBULENT TRANSPORT THEORY

TURBULENT TRANSPORT THEORY ASDEX Upgrade Max-Planck-Institut für Plasmaphysik TURBULENT TRANSPORT THEORY C. Angioni GYRO, J. Candy and R.E. Waltz, GA The problem of Transport Transport is the physics subject which studies the physical

More information

MHD WAVES AND GLOBAL ALFVÉN EIGENMODES

MHD WAVES AND GLOBAL ALFVÉN EIGENMODES MHD WVES ND GLOBL LFVÉN EIGENMODES S.E. Sharapov Euratom/CCFE Fusion ssociation, Culham Science Centre, bingdon, Oxfordshire OX14 3DB, UK S.E.Sharapov, Lecture 3, ustralian National University, Canberra,

More information

CHAPTER 5. RUDIMENTS OF HYDRODYNAMIC INSTABILITY

CHAPTER 5. RUDIMENTS OF HYDRODYNAMIC INSTABILITY 1 Lecture Notes on Fluid Dynamics (1.63J/.1J) by Chiang C. Mei, 00 CHAPTER 5. RUDIMENTS OF HYDRODYNAMIC INSTABILITY References: Drazin: Introduction to Hydrodynamic Stability Chandrasekar: Hydrodynamic

More information

Theoretical Foundation of 3D Alfvén Resonances: Time Dependent Solutions

Theoretical Foundation of 3D Alfvén Resonances: Time Dependent Solutions Theoretical Foundation of 3D Alfvén Resonances: Time Dependent Solutions Tom Elsden 1 Andrew Wright 1 1 Dept Maths & Stats, University of St Andrews DAMTP Seminar - 8th May 2017 Outline Introduction Coordinates

More information

Plasma waves in the fluid picture I

Plasma waves in the fluid picture I Plasma waves in the fluid picture I Langmuir oscillations and waves Ion-acoustic waves Debye length Ordinary electromagnetic waves General wave equation General dispersion equation Dielectric response

More information

SMR/ Summer College on Plasma Physics. 30 July - 24 August, Introduction to Magnetic Island Theory.

SMR/ Summer College on Plasma Physics. 30 July - 24 August, Introduction to Magnetic Island Theory. SMR/1856-1 2007 Summer College on Plasma Physics 30 July - 24 August, 2007 Introduction to Magnetic Island Theory. R. Fitzpatrick Inst. for Fusion Studies University of Texas at Austin USA Introduction

More information

Chapter 1. Introduction to Nonlinear Space Plasma Physics

Chapter 1. Introduction to Nonlinear Space Plasma Physics Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large

More information

28. THE GRAVITATIONAL INTERCHANGE MODE, OR g-mode. We now consider the case where the magneto-fluid system is subject to a gravitational force F g

28. THE GRAVITATIONAL INTERCHANGE MODE, OR g-mode. We now consider the case where the magneto-fluid system is subject to a gravitational force F g 8. THE GRAVITATIONAL INTERCHANGE MODE, OR g-mode We now consider the case where the magneto-fluid system is subject to a gravitational force F g =!g, where g is a constant gravitational acceleration vector.

More information

Equations of linear stellar oscillations

Equations of linear stellar oscillations Chapter 4 Equations of linear stellar oscillations In the present chapter the equations governing small oscillations around a spherical equilibrium state are derived. The general equations were presented

More information

Waves in plasma. Denis Gialis

Waves in plasma. Denis Gialis Waves in plasma Denis Gialis This is a short introduction on waves in a non-relativistic plasma. We will consider a plasma of electrons and protons which is fully ionized, nonrelativistic and homogeneous.

More information

Toroidal flow stablization of disruptive high tokamaks

Toroidal flow stablization of disruptive high tokamaks PHYSICS OF PLASMAS VOLUME 9, NUMBER 6 JUNE 2002 Robert G. Kleva and Parvez N. Guzdar Institute for Plasma Research, University of Maryland, College Park, Maryland 20742-3511 Received 4 February 2002; accepted

More information

2 The incompressible Kelvin-Helmholtz instability

2 The incompressible Kelvin-Helmholtz instability Hydrodynamic Instabilities References Chandrasekhar: Hydrodynamic and Hydromagnetic Instabilities Landau & Lifshitz: Fluid Mechanics Shu: Gas Dynamics 1 Introduction Instabilities are an important aspect

More information

The Physics of Fluids and Plasmas

The Physics of Fluids and Plasmas The Physics of Fluids and Plasmas An Introduction for Astrophysicists ARNAB RAI CHOUDHURI CAMBRIDGE UNIVERSITY PRESS Preface Acknowledgements xiii xvii Introduction 1 1. 3 1.1 Fluids and plasmas in the

More information

Amplification of magnetic fields in core collapse

Amplification of magnetic fields in core collapse Amplification of magnetic fields in core collapse Miguel Àngel Aloy Torás, Pablo Cerdá-Durán, Thomas Janka, Ewald Müller, Martin Obergaulinger, Tomasz Rembiasz Universitat de València; Max-Planck-Institut

More information

Fundamentals of Magnetic Island Theory in Tokamaks

Fundamentals of Magnetic Island Theory in Tokamaks Fundamentals of Magnetic Island Theory in Tokamaks Richard Fitzpatrick Institute for Fusion Studies University of Texas at Austin Austin, TX, USA Talk available at http://farside.ph.utexas.edu/talks/talks.html

More information

Performance limits. Ben Dudson. 24 th February Department of Physics, University of York, Heslington, York YO10 5DD, UK

Performance limits. Ben Dudson. 24 th February Department of Physics, University of York, Heslington, York YO10 5DD, UK Performance limits Ben Dudson Department of Physics, University of York, Heslington, York YO10 5DD, UK 24 th February 2014 Ben Dudson Magnetic Confinement Fusion (1 of 24) Previously... In the last few

More information

Reduced MHD. Nick Murphy. Harvard-Smithsonian Center for Astrophysics. Astronomy 253: Plasma Astrophysics. February 19, 2014

Reduced MHD. Nick Murphy. Harvard-Smithsonian Center for Astrophysics. Astronomy 253: Plasma Astrophysics. February 19, 2014 Reduced MHD Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 19, 2014 These lecture notes are largely based on Lectures in Magnetohydrodynamics by Dalton

More information

(a) (b) (c) (d) (e) (f) r (minor radius) time. time. Soft X-ray. T_e contours (ECE) r (minor radius) time time

(a) (b) (c) (d) (e) (f) r (minor radius) time. time. Soft X-ray. T_e contours (ECE) r (minor radius) time time Studies of Spherical Tori, Stellarators and Anisotropic Pressure with M3D 1 L.E. Sugiyama 1), W. Park 2), H.R. Strauss 3), S.R. Hudson 2), D. Stutman 4), X-Z. Tang 2) 1) Massachusetts Institute of Technology,

More information

AST242 LECTURE NOTES PART 5

AST242 LECTURE NOTES PART 5 AST242 LECTURE NOTES PART 5 Contents 1. Waves and instabilities 1 1.1. Sound waves compressive waves in 1D 1 2. Jeans Instability 5 3. Stratified Fluid Flows Waves or Instabilities on a Fluid Boundary

More information

Stabilization of sawteeth in tokamaks with toroidal flows

Stabilization of sawteeth in tokamaks with toroidal flows PHYSICS OF PLASMAS VOLUME 9, NUMBER 7 JULY 2002 Stabilization of sawteeth in tokamaks with toroidal flows Robert G. Kleva and Parvez N. Guzdar Institute for Plasma Research, University of Maryland, College

More information

MAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT

MAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT MAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT ABSTRACT A. G. Tarditi and J. V. Shebalin Advanced Space Propulsion Laboratory NASA Johnson Space Center Houston, TX

More information

Comparison of Kinetic and Extended MHD Models for the Ion Temperature Gradient Instability in Slab Geometry

Comparison of Kinetic and Extended MHD Models for the Ion Temperature Gradient Instability in Slab Geometry Comparison of Kinetic and Extended MHD Models for the Ion Temperature Gradient Instability in Slab Geometry D. D. Schnack University of Wisconsin Madison Jianhua Cheng, S. E. Parker University of Colorado

More information

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit

More information

Fundamentals of Fluid Dynamics: Elementary Viscous Flow

Fundamentals of Fluid Dynamics: Elementary Viscous Flow Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research

More information

Theory of Ship Waves (Wave-Body Interaction Theory) Quiz No. 2, April 25, 2018

Theory of Ship Waves (Wave-Body Interaction Theory) Quiz No. 2, April 25, 2018 Quiz No. 2, April 25, 2018 (1) viscous effects (2) shear stress (3) normal pressure (4) pursue (5) bear in mind (6) be denoted by (7) variation (8) adjacent surfaces (9) be composed of (10) integrand (11)

More information

Stability of a plasma confined in a dipole field

Stability of a plasma confined in a dipole field PHYSICS OF PLASMAS VOLUME 5, NUMBER 10 OCTOBER 1998 Stability of a plasma confined in a dipole field Plasma Fusion Center, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 Received

More information

Conservation Laws in Ideal MHD

Conservation Laws in Ideal MHD Conservation Laws in Ideal MHD Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 3, 2016 These lecture notes are largely based on Plasma Physics for Astrophysics

More information

Euler equation and Navier-Stokes equation

Euler equation and Navier-Stokes equation Euler equation and Navier-Stokes equation WeiHan Hsiao a a Department of Physics, The University of Chicago E-mail: weihanhsiao@uchicago.edu ABSTRACT: This is the note prepared for the Kadanoff center

More information

Quick Recapitulation of Fluid Mechanics

Quick Recapitulation of Fluid Mechanics Quick Recapitulation of Fluid Mechanics Amey Joshi 07-Feb-018 1 Equations of ideal fluids onsider a volume element of a fluid of density ρ. If there are no sources or sinks in, the mass in it will change

More information

School and Conference on Analytical and Computational Astrophysics November, Angular momentum transport in accretion disks

School and Conference on Analytical and Computational Astrophysics November, Angular momentum transport in accretion disks 2292-13 School and Conference on Analytical and Computational Astrophysics 14-25 November, 2011 Angular momentum transport in accretion disks Gianluigi Bodo Osservatorio Astronomico, Torino Italy Angular

More information

CHAPTER 19. Fluid Instabilities. In this Chapter we discuss the following instabilities:

CHAPTER 19. Fluid Instabilities. In this Chapter we discuss the following instabilities: CHAPTER 19 Fluid Instabilities In this Chapter we discuss the following instabilities: convective instability (Schwarzschild criterion) interface instabilities (Rayleight Taylor & Kelvin-Helmholtz) gravitational

More information

Shear instabilities. Chapter Energetics of shear instabilities

Shear instabilities. Chapter Energetics of shear instabilities Chapter 7 Shear instabilities In this final Chapter, we continue our study of the stability of fluid flows by looking at another very common source of instability, shear. By definition, shear occurs whenever

More information

Transition From Single Fluid To Pure Electron MHD Regime Of Tearing Instability

Transition From Single Fluid To Pure Electron MHD Regime Of Tearing Instability Transition From Single Fluid To Pure Electron MHD Regime Of Tearing Instability V.V.Mirnov, C.C.Hegna, S.C.Prager APS DPP Meeting, October 27-31, 2003, Albuquerque NM Abstract In the most general case,

More information

Two Fluid Dynamo and Edge-Resonant m=0 Tearing Instability in Reversed Field Pinch

Two Fluid Dynamo and Edge-Resonant m=0 Tearing Instability in Reversed Field Pinch 1 Two Fluid Dynamo and Edge-Resonant m= Tearing Instability in Reversed Field Pinch V.V. Mirnov 1), C.C.Hegna 1), S.C. Prager 1), C.R.Sovinec 1), and H.Tian 1) 1) The University of Wisconsin-Madison, Madison,

More information

Computational Astrophysics

Computational Astrophysics 16 th Chris Engelbrecht Summer School, January 2005 3: 1 Computational Astrophysics Lecture 3: Magnetic fields Paul Ricker University of Illinois at Urbana-Champaign National Center for Supercomputing

More information

PHYS 643 Week 4: Compressible fluids Sound waves and shocks

PHYS 643 Week 4: Compressible fluids Sound waves and shocks PHYS 643 Week 4: Compressible fluids Sound waves and shocks Sound waves Compressions in a gas propagate as sound waves. The simplest case to consider is a gas at uniform density and at rest. Small perturbations

More information

20. Alfven waves. ([3], p ; [1], p ; Chen, Sec.4.18, p ) We have considered two types of waves in plasma:

20. Alfven waves. ([3], p ; [1], p ; Chen, Sec.4.18, p ) We have considered two types of waves in plasma: Phys780: Plasma Physics Lecture 20. Alfven Waves. 1 20. Alfven waves ([3], p.233-239; [1], p.202-237; Chen, Sec.4.18, p.136-144) We have considered two types of waves in plasma: 1. electrostatic Langmuir

More information

Ideal Magnetohydrodynamics (MHD)

Ideal Magnetohydrodynamics (MHD) Ideal Magnetohydrodynamics (MHD) Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 1, 2016 These lecture notes are largely based on Lectures in Magnetohydrodynamics

More information

Receiver. Johana Brokešová Charles University in Prague

Receiver. Johana Brokešová Charles University in Prague Propagation of seismic waves - theoretical background Receiver Johana Brokešová Charles University in Prague Seismic waves = waves in elastic continuum a model of the medium through which the waves propagate

More information

1/27/2010. With this method, all filed variables are separated into. from the basic state: Assumptions 1: : the basic state variables must

1/27/2010. With this method, all filed variables are separated into. from the basic state: Assumptions 1: : the basic state variables must Lecture 5: Waves in Atmosphere Perturbation Method With this method, all filed variables are separated into two parts: (a) a basic state part and (b) a deviation from the basic state: Perturbation Method

More information

PLASMA ASTROPHYSICS. ElisaBete M. de Gouveia Dal Pino IAG-USP. NOTES: (references therein)

PLASMA ASTROPHYSICS. ElisaBete M. de Gouveia Dal Pino IAG-USP. NOTES:  (references therein) PLASMA ASTROPHYSICS ElisaBete M. de Gouveia Dal Pino IAG-USP NOTES:http://www.astro.iag.usp.br/~dalpino (references therein) ICTP-SAIFR, October 7-18, 2013 Contents What is plasma? Why plasmas in astrophysics?

More information

Figure 1. adiabatically. The change in pressure experienced by the parcel is. dp = -ρ o gξ

Figure 1. adiabatically. The change in pressure experienced by the parcel is. dp = -ρ o gξ 6. Internal waves Consider a continuously stratified fluid with ρ o (z) the vertical density profile. z p' ξ p ρ ρ ο (z) Figure 1. Figure by MIT OpenCourseWare. At a point P raise a parcel of water by

More information

INTERNAL GRAVITY WAVES

INTERNAL GRAVITY WAVES INTERNAL GRAVITY WAVES B. R. Sutherland Departments of Physics and of Earth&Atmospheric Sciences University of Alberta Contents Preface List of Tables vii xi 1 Stratified Fluids and Waves 1 1.1 Introduction

More information

Entropy generation and transport

Entropy generation and transport Chapter 7 Entropy generation and transport 7.1 Convective form of the Gibbs equation In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?

More information

Lecture 3: The Navier-Stokes Equations: Topological aspects

Lecture 3: The Navier-Stokes Equations: Topological aspects Lecture 3: The Navier-Stokes Equations: Topological aspects September 9, 2015 1 Goal Topology is the branch of math wich studies shape-changing objects; objects which can transform one into another without

More information

Sound Generation from Vortex Sheet Instability

Sound Generation from Vortex Sheet Instability Sound Generation from Vortex Sheet Instability Hongbin Ju Department of Mathematics Florida State University, Tallahassee, FL.3306 www.aeroacoustics.info Please send comments to: hju@math.fsu.edu When

More information

CHAPTER 8 ENTROPY GENERATION AND TRANSPORT

CHAPTER 8 ENTROPY GENERATION AND TRANSPORT CHAPTER 8 ENTROPY GENERATION AND TRANSPORT 8.1 CONVECTIVE FORM OF THE GIBBS EQUATION In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?

More information

Soft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies

Soft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed

More information

Goals of this Chapter

Goals of this Chapter Waves in the Atmosphere and Oceans Restoring Force Conservation of potential temperature in the presence of positive static stability internal gravity waves Conservation of potential vorticity in the presence

More information

Chapter 5. Sound Waves and Vortices. 5.1 Sound waves

Chapter 5. Sound Waves and Vortices. 5.1 Sound waves Chapter 5 Sound Waves and Vortices In this chapter we explore a set of characteristic solutions to the uid equations with the goal of familiarizing the reader with typical behaviors in uid dynamics. Sound

More information

Problem Set Number 01, MIT (Winter-Spring 2018)

Problem Set Number 01, MIT (Winter-Spring 2018) Problem Set Number 01, 18.377 MIT (Winter-Spring 2018) Rodolfo R. Rosales (MIT, Math. Dept., room 2-337, Cambridge, MA 02139) February 28, 2018 Due Thursday, March 8, 2018. Turn it in (by 3PM) at the Math.

More information

35. RESISTIVE INSTABILITIES: CLOSING REMARKS

35. RESISTIVE INSTABILITIES: CLOSING REMARKS 35. RESISTIVE INSTABILITIES: CLOSING REMARKS In Section 34 we presented a heuristic discussion of the tearing mode. The tearing mode is centered about x = 0, where F x that B 0 ( ) = k! B = 0. Note that

More information

Bursty Transport in Tokamaks with Internal Transport Barriers

Bursty Transport in Tokamaks with Internal Transport Barriers Bursty Transport in Tokamaks with Internal Transport Barriers S. Benkadda 1), O. Agullo 1), P. Beyer 1), N. Bian 1), P. H. Diamond 3), C. Figarella 1), X. Garbet 2), P. Ghendrih 2), V. Grandgirard 1),

More information

Generalized Solovev equilibrium with sheared flow of arbitrary direction and stability consideration

Generalized Solovev equilibrium with sheared flow of arbitrary direction and stability consideration Generalized Solovev equilibrium with sheared flow of arbitrary direction and stability consideration D.A. Kaltsas and G.N. Throumoulopoulos Department of Physics, University of Ioannina, GR 451 10 Ioannina,

More information

Creation and destruction of magnetic fields

Creation and destruction of magnetic fields HAO/NCAR July 30 2007 Magnetic fields in the Universe Earth Magnetic field present for 3.5 10 9 years, much longer than Ohmic decay time ( 10 4 years) Strong variability on shorter time scales (10 3 years)

More information

Space Plasma Physics Thomas Wiegelmann, 2012

Space Plasma Physics Thomas Wiegelmann, 2012 Space Plasma Physics Thomas Wiegelmann, 2012 1. Basic Plasma Physics concepts 2. Overview about solar system plasmas Plasma Models 3. Single particle motion, Test particle model 4. Statistic description

More information

Chapter 5. The Differential Forms of the Fundamental Laws

Chapter 5. The Differential Forms of the Fundamental Laws Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations

More information

7 The Navier-Stokes Equations

7 The Navier-Stokes Equations 18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and

More information

Fluid Dynamics Exercises and questions for the course

Fluid Dynamics Exercises and questions for the course Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r

More information

The Linear Theory of Tearing Modes in periodic, cyindrical plasmas. Cary Forest University of Wisconsin

The Linear Theory of Tearing Modes in periodic, cyindrical plasmas. Cary Forest University of Wisconsin The Linear Theory of Tearing Modes in periodic, cyindrical plasmas Cary Forest University of Wisconsin 1 Resistive MHD E + v B = ηj (no energy principle) Role of resistivity No frozen flux, B can tear

More information

3.5 Vorticity Equation

3.5 Vorticity Equation .0 - Marine Hydrodynamics, Spring 005 Lecture 9.0 - Marine Hydrodynamics Lecture 9 Lecture 9 is structured as follows: In paragraph 3.5 we return to the full Navier-Stokes equations (unsteady, viscous

More information

General introduction to Hydrodynamic Instabilities

General introduction to Hydrodynamic Instabilities KTH ROYAL INSTITUTE OF TECHNOLOGY General introduction to Hydrodynamic Instabilities L. Brandt & J.-Ch. Loiseau KTH Mechanics, November 2015 Luca Brandt Professor at KTH Mechanics Email: luca@mech.kth.se

More information

Impact of diverted geometry on turbulence and transport barrier formation in 3D global simulations of tokamak edge plasma

Impact of diverted geometry on turbulence and transport barrier formation in 3D global simulations of tokamak edge plasma 1 Impact of diverted geometry on turbulence and transport barrier formation in 3D global simulations of tokamak edge plasma D. Galassi, P. Tamain, H. Bufferand, C. Baudoin, G. Ciraolo, N. Fedorczak, Ph.

More information

Numerical methods for the Navier- Stokes equations

Numerical methods for the Navier- Stokes equations Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:

More information

Vortex motion. Wasilij Barsukow, July 1, 2016

Vortex motion. Wasilij Barsukow, July 1, 2016 The concept of vorticity We call Vortex motion Wasilij Barsukow, mail@sturzhang.de July, 206 ω = v vorticity. It is a measure of the swirlyness of the flow, but is also present in shear flows where the

More information

NIMROD simulations of dynamo experiments in cylindrical and spherical geometries. Dalton Schnack, Ivan Khalzov, Fatima Ebrahimi, Cary Forest,

NIMROD simulations of dynamo experiments in cylindrical and spherical geometries. Dalton Schnack, Ivan Khalzov, Fatima Ebrahimi, Cary Forest, NIMROD simulations of dynamo experiments in cylindrical and spherical geometries Dalton Schnack, Ivan Khalzov, Fatima Ebrahimi, Cary Forest, 1 Introduction Two experiments, Madison Plasma Couette Experiment

More information

Chapter 4. Gravity Waves in Shear. 4.1 Non-rotating shear flow

Chapter 4. Gravity Waves in Shear. 4.1 Non-rotating shear flow Chapter 4 Gravity Waves in Shear 4.1 Non-rotating shear flow We now study the special case of gravity waves in a non-rotating, sheared environment. Rotation introduces additional complexities in the already

More information

Gravity Waves. Lecture 5: Waves in Atmosphere. Waves in the Atmosphere and Oceans. Internal Gravity (Buoyancy) Waves 2/9/2017

Gravity Waves. Lecture 5: Waves in Atmosphere. Waves in the Atmosphere and Oceans. Internal Gravity (Buoyancy) Waves 2/9/2017 Lecture 5: Waves in Atmosphere Perturbation Method Properties of Wave Shallow Water Model Gravity Waves Rossby Waves Waves in the Atmosphere and Oceans Restoring Force Conservation of potential temperature

More information

Waves and characteristics: Overview 5-1

Waves and characteristics: Overview 5-1 Waves and characteristics: Overview 5-1 Chapter 5: Waves and characteristics Overview Physics and accounting: use example of sound waves to illustrate method of linearization and counting of variables

More information