Collisionless magnetic reconnection with arbitrary guide-field

Size: px
Start display at page:

Download "Collisionless magnetic reconnection with arbitrary guide-field"

Transcription

1 Collisionless magnetic reconnection with arbitrary guide-field Richard Fitzpatrick Center for Magnetic Reconnection Studies Institute for Fusion Studies Department of Physics University of Texas at Austin Austin, TX 7871 USA Franco Porcelli Burning Plasma Research Group Dipartimento di Energetica Politecnico di Torino 1019 Torino Italy A new set of reduced equations governing two-dimensional, two-fluid, collisionless magnetic reconnection with arbitrary guide-field is derived. These equations represent a significant advance in magnetic reconnection theory, since the existing reduced equations used to investigate collisionless reconnection are only valid in the large guide-field limit. The new equations are used to calculate the linear growth-rate of a strongly unstable, spontaneously reconnecting, plasma instability, as well as the general linear dispersion relation for such an instability. rfitzp@farside.ph.utexas.edu porcelli@polito.it 1

2 I. INTRODUCTION Magnetic reconnection is a fundamental physical phenomenon which occurs in a wide variety of laboratory and space plasmas: e.g., magnetic fusion experiments, 1 the solar corona, and the Earth s magneto-tail. 3 The reconnection process gives rise to a change in magnetic field-line topology with an accompanying release of magnetic energy. Conventional collisional magnetohydrodynamical (MHD) theory is capable of accounting for magnetic reconnection, but generally predicts reconnection rates which are many orders of magnitude smaller than those observed in high-temperature plasmas. 4 On the other hand, more sophisticated plasma models that neglect collisions (which are irrelevant in high temperature plasmas), and treat electrons and ions as separate fluids, yield much faster reconnection rates which are fairly consistent with observations. 5,6 Magnetic reconnection is generally thought of as a basically two-dimensional (-D) phenomenon in which oppositely directed magnetic field-lines spontaneously merge together. However, in many cases of interest, there is also an essentially uniform magnetic field directed perpendicular to the merging field-lines. This latter field is usually termed the guide-field. Unfortunately, the addition of a guide-field significantly complicates the equations governing -D, two-fluid, collisionless magnetic reconnection. In particular, the guide-field couples the compressional Alfvén (CA) wave into these equations, despite the fact that this wave does not play a particularly significant role in magnetic reconnection. Furthermore, the CA wave becomes increasingly problematic (especially in numerical simulations) as the strength of the guide-field increases. Now, it is common practice in conventional MHD theory to reduce the equations governing MHD phenomena by removing the CA wave. 7 Indeed, this approach has been very successful over the years. In this paper, we present a new set of reduced equations governing -D, two-fluid, collisionless magnetic reconnection which is valid for both large and small guide-fields (measured with respect to the reconnecting field). Our new equations represent a significant advance in magnetic reconnection theory, since the existing reduced equations used to investigate collisionless reconnection are only valid

3 in the large guide-field limit. 8,9 In addition, we calculate the linear growth-rate of a strongly unstable, spontaneously reconnecting, plasma instability from our new equations, as well as the general linear dispersion relation for such an instability. II. DERIVATION OF REDUCED EQUATIONS A. Basic equations Standard right-handed Cartesian coordinates (x, y, z) are adopted. In the following, we employ a normalization scheme such that all lengths are normalized to a convenient length-scale a, all magnetic fields to a convenient scale field-strength B 0, and all times to τ A = a/v A, where V A = B 0 / µ 0 n m i. Here, n is the electron/ion number density, which is assumed to be constant. We also assume that the ions are singly charged. Furthermore, m i is the ion mass. Finally, the electrons are assumed to be hot, and the ions cold. The normalized electron equation of motion takes the form 10 [ ] Ve ɛ d i t + (V e )V e = d i P (E + V e B), (1) where E is the electric field, B the magnetic field, V e the electron velocity, P the electron pressure (which is assumed to be scalar), d i = (c/ω p i )/a (where ω p i = n e /ɛ 0 m i ) the normalized collisionless ion skin-depth, and ɛ = m e /m i the electron-ion mass ratio. The normalized sum of the electron and ion equations of motion yields (neglecting electron inertia) V t + (V )V = P + J B, () where V = V e +d i J is the ion velocity, and J = B. Also, B = 0, and E = B/ t. Our system of equations is closed by a simple energy equation: where Γ = 5/3 is the ratio of specific heats. ( ) t + V P = Γ V P, (3) 3

4 B. -D equations Suppose that / z 0. Without loss of generality we can write B = ψ ẑ + B z ẑ, (4) and E z = ψ/ t. The z-component of Eq. (1) yields ( ) ( ) ɛ d i t + V e V e z = t + V e ψ. (5) The z-component of the curl of (1) reduces to ( ) ( ) ɛ d i t + V e + V e ω e z = t + V e B z [V e z, ψ] + V e B z, (6) where and ω e = V e, (7) [A, B] A B ẑ. (8) To lowest order in ɛ, Eqs. (5) and (6) can be rearranged to give ( ) t + V ψ e = d i [ψ, B z ], (9) ( ) t + V B e z = [V z, ψ] + d i [ ψ, ψ] V B e z, (10) respectively. Here, where d e = ɛ d i is the collisionless electron skin-depth. ψ e = ψ d e ψ, (11) B e z = B z d e B z, (1) Equation () can be written ( ) t + V V = P B z B z ψ ψ, (13) ( ) t + V V z = [B z, ψ], (14) where V = (V x, V y, 0). 4

5 C. Reduction scheme Let P = P (0) + B (0) p 1 + p, (15) B z = B (0) + b z, (16) where P (0) and B (0) are uniform constants, p 1, p, and b z are all O(1), and P (0) B (0) 1. (17) The purpose of the above ordering is to make the compressional Alfvén (CA) wave propagate much faster than any other wave in the system. When this occurs, the CA wave effectively decouples from these other waves, and can therefore be eliminated from our system of equations. We expect the rapidly propagating CA wave to maintain approximate force balance within the plasma at all times: i.e., P + B z B z + ψ ψ 0. (18) Hence, it follows from (17) that p 1 b z. (19) We also expect the CA wave to make the plasma flow almost incompressible, so that V φ ẑ. (0) The curl of Eq. (13) yields where U = ψ. Now, to lowest order, Eqs. (3) and (10) become U t = [φ, U] + [ ψ, ψ], (1) 5

6 b e z = [φ, b z ] + [V z, ψ] + d i [ ψ, ψ] V B (0), t () p 1 t = [φ, p 1] Γ P (0) B (0) V, (3) where b e z = b z d e b z. We can eliminate V between the above two equations, making use of (19), to give Z e t where Z = b z /c β, Z e = Z c β d e Z, c β = = [φ, Z] + c β [V z, ψ] + d β [ ψ, ψ], (4) β/(1 + β), d β = c β d i, and β = Γ P (0) /[B (0) ]. D. Reduced equations where Our new set of reduced collisionless equations takes the form ψ e = [φ, ψ] + d β [ψ, Z], t (5) Z e = [φ, Z] + c β [V z, ψ] + d β [ ψ, ψ], t (6) U t = [φ, U] + [ ψ, ψ], (7) V z t = [φ, V z] + c β [Z, ψ], (8) φ = U, (9) ψ e = ψ d e ψ, (30) Z e = Z c β d e Z, (31) and d β = c β d i, (3) β c β = 1 + β, (33) β = Γ P (0) [B (0) ]. (34) 6

7 Note that the parameter β measures the strength of the guide-field. The large guide-field limit corresponds to β 1, whereas β 1 corresponds to the small guide-field limit. Note also that d i only occurs in our reduced equations in the combination d β = [β/(1+β)] 1/ d i. In the small guide-field limit, d β d i. However, in the large guide-field limit d β β d i ρ s. Eqs. (5) (8) are generalized versions of the reduced equations presented in Ref. 11 which take electron inertia into account. E. Zero guide-field limit In the zero guide-field limit, B (0) 0, which corresponds to β, we adopt the modified ordering P = P (0) + p 1, (35) B z = b z, (36) where P (0) 1 is a uniform constant, and p 1, b z are O(1). As before, the purpose of this ordering scheme is to make the compressional Alfvén wave propagate much faster than any other wave in the system. From Eq. (3), the above ordering scheme implies that V = 0. In this limit, the remaining two-fluid equations readily yield Eqs. (5) (8), with c β = 1 and d β = d i. We conclude that Eqs. (5) (8) hold in the zero, small, and large guide-field limits, as long as the compressional Alfvén wave propagates more rapidly than any other wave in the system. III. LARGE GROWTH-RATE A. Introduction Let us calculate the linear growth-rate of a reconnecting instability from Eqs. (5) (8) in the limit in which the tearing stability index 1 is very large. 7

8 Let where denotes a perturbed quantity. ψ(x, y, t) = 1 x + ψ(x) e i (k y+γ t), (37) Z(x, y, t) = Z(x) e i (k y+γ t), (38) φ(x, y, t) = φ(x) e i (k y+γ t), (39) V z (x, y, t) = Ṽz(x) e i (k y+γ t), (40) Throughout most of the plasma, the perturbation is governed by the equations of ideal- MHD, which yield (d /dx k ) ψ = 0. Ideal-MHD breaks down in a thin layer [of thickness O(d β ) or O(d e )] centered on x = 0. 1 Both d e and d β are assumed to be much less than unity. The layer solution must be asymptotically matched to the ideal-mhd solution via the tearing stability index, = [d ln ψ/dx] For x 1, Eqs. (5) (8) yield the following linearized layer equations: ( g 1 de d ) ψ = i x dx i d β x Z, (41) ( g 1 cβ d e d dx + c ) β g x Z = i d β x d ψ dx, (4) where g = γ/k. Let g d φ dx = i x d ψ dx, (43) φ(p) = + φ(x) e i p x dx, (44) etc. The above Fourier transform should be understood in the generalized sense, since φ(x) is not square-integrable. The Fourier transformed layer equations take the form g ( g (1 + de p ) ψ = d( φ d β Z), (45) dp 1 + cβ de p c β d ) d(p Z ψ) = d g dp β, (46) dp g p d(p φ ψ) =, (47) dp 8

9 respectively. B. Inner layer In the limit p c β /g, d 1, Eqs. (45) (47) reduce to β ( ) d r dw Q (1 + c dr 1 + r β dr r ) W = 0, (48) where r = d e p, W = p φ/(1 + c β d e p ), and Q = g/d β. In the limit r 0, the asymptotic solution of the above equation takes the form W = W 0 [ r 1 ν + α r ν + O(r 1 ν ) ], (49) where ν = (1/) ( Q 1), and α is determined by the condition that the solution to (48) be well-behaved as r. Suppose that Q 1, which implies that 0 < ν 1. For r < 1, Eq. (48) yields ( ) d r dw 0. (50) dr 1 + r dr Integrating, and comparing with (49) (recalling that ν 0), we obtain W = W 0 [ (r 1 r) + α + O(r ) ] (51) for r < 1. For r 1, Eq. (48) reduces to d W dr Q (1 + c β r ) W 0, (5) which is a parabolic cylinder equation. 13 The solution which is well-behaved as r is W U(Q/ c β, Q c β r), (53) where U(a, x) is a standard parabolic cylinder function. 13 For 1 r (Q c β ) 1/, we obtain 13 9

10 Γ(1/4 + Q/4 c W = W 1 β ) Q c β Γ(3/4 + Q/4 c β ) r + O(r ), (54) where Γ(x) is a standard gamma function. 13 Finally, comparing with Eq. (51), we find α = Q 1 G(Q/c β ), (55) where G(x) = x Γ(1/4 + x/4) Γ(3/4 + x/4). (56) Note that G(x) 1 as x. C. Outer layer Suppose that x d e. In this limit, the real-space layer equations (41) (43) reduce to [ d x d ( g dβ d ) φ (g + x ) d φ ] = 0. (57) dx dx g + cβ x dx dx Now, we expect the solution to the above equation to match to the ideal-mhd solution at (relatively) large x. Thus, φ φ [ ( )] 1 1 x + O x (58) at large x, where is the standard tearing stability index. 1 Let Y = g (d φ/dx)/ φ 0. It follows that at large x. Thus, integrating Eq. (57), we obtain Y g 1 x (59) x d ( g d ) β dy (g + x ) Y = g dx g + cβ x dx. (60) In the large limit (i.e., 1/ 0), the above equation reduces to s d ( ) 1 dy Q (1 + s ) Y = 0, (61) ds 1 + c β s ds 10

11 where s = x/g. Let V = (1 + c β s ) 1 dy/ds. The above equation can be transformed into ( ) d s dv Q (1 + c ds 1 + s β ds s ) V = 0. (6) Note that this outer layer equation has exactly the same form as the inner layer equation, (48). The small s asymptotic solution of Eq. (6) is V = V 0 [ s 1 ν + α 1 s ν + O(s 1 ν ) ], (63) where α 1 is determined by the condition that V not blow-up as s. However, since Eq. (6) has exactly the same form as Eq. (48), we immediately deduce that α 1 = α = Q 1 G(Q/c β ). (64) Fourier transforming (63), we obtain the asymptotic solution ( ) g 1+ νr W = W [α Γ(1 + ν) cos(ν π/) r 1 ν Γ( ν) sin(ν π/) ν + O(r )] 1 ν de (65) in Fourier space. Comparison with (49) yields α = π in the limit ν 0. It follows from (55) that g d e (66) g 3 = π d e d β G (g/c β d β ). (67) Our derivation of this expression is valid provided d e d β Q 1, (68) and d β d e. These constraints all reduce to the single constraint β m e /m i. D. Discussion We have derived the linear dispersion relation (67) from our reduced equations in the large limit. This dispersion relation yields 14 11

12 ( ) 1/3 g = de 1/3 d /3 β (69) π for m e /m i β (m e /m i ) 1/4, and g = 1 Γ(1/4) d e d i (70) π Γ(3/4) for (m e /m i ) 1/4 β. Note that our derivation of (67) is valid for all β values significantly above m e /m i. The dispersion relation (67) was first obtained by Mirnov, et al. 15 However, the derivation presented in Ref. 15 is only valid for β in the range (m e /m i ) β 1. The improved derivation presented in this paper removes the upper constraint on β, and, thereby, confirms that (67) is valid in both the large and small guide-field limits. IV. SMALL GROWTH-RATE Let us now calculate the linear growth-rate of a reconnecting instability from Eqs. (5) (8) in the limit in which the tearing stability index is relatively small. In the small limit, we expect the tearing layer to be localized on the d e length-scale, and to be governed solely by electron physics. Taking Eqs. (45) and (46), and neglecting the ion terms (i.e., the terms involving φ and d Z/dp ), we obtain ( d r d Z ) Q (1 + c dr 1 + r β dr r ) Z = 0, (71) and ψ = d β d Z 1 + de p dp, (7) where r = d e p and Q = g/d β. The solution to (71) which is well-behaved as r has the small r expansion Z = Z 0 [r 1 + α + O(r)], (73) assuming that Q 1. Here, α is given by Eq. (55), since Eq. (71) has exactly the same form as Eq. (48). Inverse Fourier transforming, we obtain 1

13 Z = Z [ π 0 d e α + 1 ( )] 1 x + O, (74) x and ψ x Z = Z [ ( )] π x 1 0 d e α O. (75) x It follows by matching to the standard ideal-mhd solution, that = π/d e α, or ψ = ψ 0 [ x O ( 1 x ) ], (76) = π d e d β g G(g/c β d β ). (77) Hence, we conclude that 14 g = d β π d e (78) for β ( d e ), and g = ( ) Γ(1/4) d i ( d e ) (79) π Γ(3/4) for β ( d e ). The dispersion relation (77) was first obtained by Mirnov, et al. 15 V. GENERAL DISPERSION RELATION Let us now derive a general dispersion relation from Eqs. (5) (8) which holds for large, small, and intermediate values of. For finite, Eqs. (60) and (6) yield s d ( ) 1 dy Q (1 + s ) Y = ˆλ ds 1 + cβ H Q, (80) s ds ( ) d s dv Q (1 + c ds 1 + s β s ) V = ˆλ H Q d ( ) 1, (81) ds ds 1 + s 13

14 where s = x/g, λ H = π/, and ˆλ H = (λ H /d β )(π/). Note, from Eq. (58), that Y (s) is subject to the constraint 0 Y (s) ds = 1. (8) Our basic expansion parameter is H = G(Q/c β )/Q [see Eq. (56)]. Note that H 1 provided Q 1. It can easily be verified from the following analysis that V O(Q ) and ˆλ H < O(Q/H). It follows that in Eqs. (80) and (81) the inhomogenous term on the right-hand side can be neglected for s H 1, whereas the second term on the left-hand side can be neglected for s H. In the limit s H, Eq. (81) reduces to which can be solved to give ( ) d s dv ˆλ ds 1 + s H Q d ( ) 1, (83) ds ds 1 + s V a 0 ˆλ H Q s + a 1 + a 0 s, (84) where a 0 (Q, c β ) and a 1 (Q, c β ) are constants of integration. Note that the first two terms on the right-hand side of the above expression are negligible for s H 1, and vice versa. In the limit s H 1, the above expression can be asymptotically matched to the inner layer solution described in Sect. III B. In real space, this solution yields (for Q 1) It therefore follows that V V 0 [ π ] d e G(Q/c β ) 1 d β Q s O(s). (85) a 0 ˆλ H Q a 1 = π d e G(Q/c β ). (86) d β Q In the limit s H 1, Eq. (81) reduces to d V ds Q (1 + c β s ) V 0. (87) This is a parabolic cylinder equation. 13 The solution which is well-behaved as s has the expansion 14

15 ] Q V V 0 [1 G(Q/c β ) s + O(s ) for s H. Asymptotic matching to Eq. (84) yields Note that for s H 1 we can write where U(a, x) is a parabolic cylinder function. 13 (88) a 0 = Q a 1 G(Q/c β ). (89) U(Q/ c β, Q c β s) V a 1, (90) U(Q/ c β, 0) Dividing Eq. (80) by 1 + s and then integrating from s = 0 to, we obtain 0 s dv 1 + s ds ds = Q ˆλ H Q π, (91) where use has been made of Eq. (8). The dominant contribution to the integral in the above expression comes from the region s 1. Furthermore, ˆλ H < O(Q/H), so the second term on the right-hand side is negligible. Hence, it follows from Eq. (90) that a 1 U(Q/ c β, s= Q c β s) Q, (9) U(Q/ c β, 0) or s=0 Equations (86), (89), and (93) can be combined to give Q G(Q/c β ) ˆλ H Q = π Hence, our general dispersion relation takes the form λ H = π a 1 Q. (93) Q d β G(Q/c β ) d e d e G(Q/c β ), (94) d β Q G(Q/c β ), (95) Q where λ H = π/. This result is valid provided Q 1, which implies that λ H d β, and β m e /m i. The above dispersion relation is consistent with the dispersion relations (67) and (77) which we previously derived in the large and small limits, respectively. Equation (95) was first obtained by Mirnov, et al. 15 using a heuristic argument. We, on the other hand, have derived this result in a rigorous fashion. In the limit β 1, the above expression reduces to a dispersion relation previously obtained by Porcelli

16 VI. SUMMARY We have derived a novel set of reduced equations, (5) (8), governing -D, two-fluid, collisionless magnetic reconnection with (essentially) arbitrary guide-field. These equations represent a major advance in reconnection theory, since the reduced equations previously used to investigate collisionless reconnection 8,9 are only valid in the large guide-field limit. Using our equations, we have calculated the linear growth-rate of a highly unstable, spontaneously reconnecting, plasma instability. Our expression for the growth-rate, (67), is valid in both the large and small guide-field limits. This expression was previously obtained by Mirnov, et al. 15 However, the derivation given in Ref. 15 is only valid in the large guidefield limit. Interestingly enough, expression (67) indicates that the guide-field ceases to have any effect on the linear growth-rate once the plasma β calculated with the guide-field significantly exceeds the rather small critical value (m e /m i ) 1/ (for hydrogen). We have also derived the general linear dispersion relation, (95), of a spontaneously reconnecting plasma instability. This result was previously obtained Mirnov, et al. 15 using a heuristic argument. In future work, we intend to investigate the nonlinear properties of our new equations. Acknowledgments This research was funded by the U.S. Department of Energy under contract DE-FG05-96ER as well as via Cooperative Agreement No. DE-FC0-01ER-5465 under the auspices of the program for Scientific Discovery through Advanced Computing (SIDAC). FP was supported in part by the Euratom Human Capital and Mobility Scheme. 1 F.L. Waelbroeck, Phys. Plasmas 1, 37 (1989). K. Shibata, Adv. Space Res. 17, 9 (1996). 16

17 3 V.M. Vasyliunas, Rev. Geophy. Space Phys. 13, 303 (1975). 4 D. Biskamp, Phys. Fluids 9, 150 (1986). 5 B. Coppi, Phys. Rev. Lett. 11, 6 (1964). 6 M. Ottaviani, and F. Porcelli, Phys. Rev. Lett. 71, 380 (1993). 7 H.R. Strauss, Phys. Fluids 0, 1357 (1977). 8 R.D. Hazeltine, M. Kotschenreuther, and P.G. Morrison, Phys. Fluids 8, 466 (1985). 9 D. Grasso, F. Pegoraro, F. Porcelli, and F. Califano, Plasma Phys. Control. Fusion 41, 1497 (1999). 10 S.I. Braginskii, in Reviews of plasma physics (Consultants Bureau, New York NY, 1965), Vol. 1, p R. Fitzpatrick, Scaling of forced magnetic reconnection in the Hall-magnetohydrodynamical Taylor problem with arbitrary guide-field, submitted to Physics of Plasmas (004). 1 H.P. Furth, J. Killeen, and M.N. Rosenbluth, Phys. Fluids 6, 459 (1963). 13 M. Abramowitz, and I.A. Stegun, Handbook of mathematical functions (Dover, New York NY, 1965). 14 F. Porcelli, Phys. Rev. Lett. 66, 45 (1991). 15 V.V. Mirnov, C.C. Hegna, and S.C. Prager, Two fluid tearing instability in force free magnetic configuration, submitted to Physics of Plasmas (004). 17

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

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

Two-fluid magnetic island dynamics in slab geometry

Two-fluid magnetic island dynamics in slab geometry Two-fluid magnetic island dynamics in slab geometry Richard Fitzpatrick and François L. Waelbroeck Institute for Fusion Studies Department of Physics University of Texas at Austin Austin, TX 78712 A novel

More information

Gyrokinetic Simulations of Tearing Instability

Gyrokinetic Simulations of Tearing Instability Gyrokinetic Simulations of Tearing Instability July 6, 2009 R. NUMATA A,, W. Dorland A, N. F. Loureiro B, B. N. Rogers C, A. A. Schekochihin D, T. Tatsuno A rnumata@umd.edu A) Center for Multiscale Plasma

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

Supersonic drift-tearing magnetic islands in tokamak plasmas

Supersonic drift-tearing magnetic islands in tokamak plasmas Supersonic drift-tearing magnetic islands in tokamak plasmas R. Fitzpatrick, and F.L. Waelbroeck Institute for Fusion Studies Department of Physics University of Texas at Austin Austin, TX 78712 A two-fluid

More information

Effect of local E B flow shear on the stability of magnetic islands in tokamak plasmas

Effect of local E B flow shear on the stability of magnetic islands in tokamak plasmas Effect of local E B flow shear on the stability of magnetic islands in tokamak plasmas R. Fitzpatrick and F. L. Waelbroeck Citation: Physics of Plasmas (1994-present) 16, 052502 (2009); doi: 10.1063/1.3126964

More information

Resistive MHD, reconnection and resistive tearing modes

Resistive MHD, reconnection and resistive tearing modes DRAFT 1 Resistive MHD, reconnection and resistive tearing modes Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford OX1 3NP, UK (This version is of 6 May 18 1. Introduction

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

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

Issues in Neoclassical Tearing Mode Theory

Issues in Neoclassical Tearing Mode Theory Issues in Neoclassical Tearing Mode Theory Richard Fitzpatrick Institute for Fusion Studies University of Texas at Austin Austin, TX Tearing Mode Stability in Tokamaks According to standard (single-fluid)

More information

Fast Secondary Reconnection and the Sawtooth Crash

Fast Secondary Reconnection and the Sawtooth Crash Fast Secondary Reconnection and the Sawtooth Crash Maurizio Ottaviani 1, Daniele Del Sarto 2 1 CEA-IRFM, Saint-Paul-lez-Durance (France) 2 Université de Lorraine, Institut Jean Lamour UMR-CNRS 7198, Nancy

More information

Measuring from electron temperature fluctuations in the Tokamak Fusion Test Reactor

Measuring from electron temperature fluctuations in the Tokamak Fusion Test Reactor PHYSICS OF PLASMAS VOLUME 5, NUMBER FEBRUARY 1998 Measuring from electron temperature fluctuations in the Tokamak Fusion Test Reactor C. Ren, a) J. D. Callen, T. A. Gianakon, and C. C. Hegna University

More information

High-m Multiple Tearing Modes in Tokamaks: MHD Turbulence Generation, Interaction with the Internal Kink and Sheared Flows

High-m Multiple Tearing Modes in Tokamaks: MHD Turbulence Generation, Interaction with the Internal Kink and Sheared Flows TH/P3-3 High-m Multiple Tearing Modes in Tokamaks: MHD Turbulence Generation, Interaction with the Internal Kink and Sheared Flows A. Bierwage 1), S. Benkadda 2), M. Wakatani 1), S. Hamaguchi 3), Q. Yu

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

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

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

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

Computational Methods in Plasma Physics

Computational Methods in Plasma Physics Computational Methods in Plasma Physics Richard Fitzpatrick Institute for Fusion Studies University of Texas at Austin Purpose of Talk Describe use of numerical methods to solve simple problem in plasma

More information

EXACT STEADY STATE RECONNECTION SOLUTIONS IN WEAKLY COLLISIONAL PLASMAS

EXACT STEADY STATE RECONNECTION SOLUTIONS IN WEAKLY COLLISIONAL PLASMAS The Astrophysical Journal, 617:1353 1360, 2004 December 20 # 2004. The American Astronomical Society. All rights reserved. Printed in U.S.A. EXACT STEADY STATE RECONNECTION SOLUTIONS IN WEAKLY COLLISIONAL

More information

An Explosive Scaling Law for Nonlinear Magnetic Reconnection and Its Insensitivity to Microscopic Scales )

An Explosive Scaling Law for Nonlinear Magnetic Reconnection and Its Insensitivity to Microscopic Scales ) An Explosive Scaling Law for Nonlinear Magnetic Reconnection and Its Insensitivity to Microscopic Scales ) Makoto IROTA Institute of Fluid Science, Tohoku University, Sendai, Miyagi 980-8577, Japan (Received

More information

Self-Organization of Plasmas with Flows

Self-Organization of Plasmas with Flows Self-Organization of Plasmas with Flows ICNSP 2003/ 9/10 Graduate School of Frontier Sciences,, National Institute for Fusion Science R. NUMATA, Z. YOSHIDA, T. HAYASHI ICNSP 2003/ 9/10 p.1/14 Abstract

More information

INTERACTION OF DRIFT WAVE TURBULENCE AND MAGNETIC ISLANDS

INTERACTION OF DRIFT WAVE TURBULENCE AND MAGNETIC ISLANDS INTERACTION OF DRIFT WAVE TURBULENCE AND MAGNETIC ISLANDS A. Ishizawa and N. Nakajima National Institute for Fusion Science F. L. Waelbroeck, R. Fitzpatrick, W. Horton Institute for Fusion Studies, University

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

Gyrokinetic simulations of collisionless magnetic reconnection

Gyrokinetic simulations of collisionless magnetic reconnection PHYSICS OF PLASMAS 14, 09110 007 Gyrokinetic simulations of collisionless magnetic reconnection B N Rogers a and S Kobayashi Department of Physics and Astronomy, Dartmouth College, Hanover, New Hampshire

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

Reconnection and the Formation of Magnetic Islands in MHD Models

Reconnection and the Formation of Magnetic Islands in MHD Models Reconnection and the Formation of Magnetic Islands in MHD Models N. F. Loureiro, D. A. Uzdensky, A. A. Schekochihin, R. Samtaney and S. C. Cowley Yosemite 2010 Reconnection Workshop Introduction (I) In

More information

Models for Global Plasma Dynamics

Models for Global Plasma Dynamics Models for Global Plasma Dynamics F.L. Waelbroeck Institute for Fusion Studies, The University of Texas at Austin International ITER Summer School June 2010 Outline 1 Models for Long-Wavelength Plasma

More information

A Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1

A Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1 A Comparison between the Two-fluid Plasma Model and Hall-MHD for Captured Physics and Computational Effort 1 B. Srinivasan 2, U. Shumlak Aerospace and Energetics Research Program University of Washington,

More information

A numerical study of forced magnetic reconnection in the viscous Taylor problem

A numerical study of forced magnetic reconnection in the viscous Taylor problem A numerical study of forced magnetic reconnection in the viscous Taylor problem Richard Fitzpatrick Citation: Physics of Plasmas (1994-present) 10, 2304 (2003); doi: 10.1063/1.1574516 View online: http://dx.doi.org/10.1063/1.1574516

More information

Magnetic Reconnection: Recent Developments and Future Challenges

Magnetic Reconnection: Recent Developments and Future Challenges Magnetic Reconnection: Recent Developments and Future Challenges A. Bhattacharjee Center for Integrated Computation and Analysis of Reconnection and Turbulence (CICART) Space Science Center, University

More information

Two-fluid magnetic island dynamics in slab geometry. II. Islands interacting with resistive walls or resonant magnetic perturbations

Two-fluid magnetic island dynamics in slab geometry. II. Islands interacting with resistive walls or resonant magnetic perturbations To-fluid magnetic isl dynamics in slab geometry. II. Isls interacting ith resistive alls or resonant magnetic perturbations Richard Fitzpatrick François L. Waelbroeck Citation: Physics of Plasmas (199-present)

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

Finite-time singularity formation at a magnetic neutral line in Hall magnetohydrodynamics

Finite-time singularity formation at a magnetic neutral line in Hall magnetohydrodynamics Finite-time singularity formation at a magnetic neutral line in Hall magnetohydrodynamics Yuri E. Litvinenko, Liam C. McMahon Department of Mathematics, University of Waikato, P. B. 3105, Hamilton, New

More information

Domain-decomposed Fully Coupled Implicit Methods for a Magnetohydrodynamics Problem

Domain-decomposed Fully Coupled Implicit Methods for a Magnetohydrodynamics Problem Domain-decomposed Fully Coupled Implicit Methods for a Magnetohydrodynamics Problem Serguei Ovtchinniov 1, Florin Dobrian 2, Xiao-Chuan Cai 3, David Keyes 4 1 University of Colorado at Boulder, Department

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

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

Anisotropic viscous dissipation in compressible magnetic X-points

Anisotropic viscous dissipation in compressible magnetic X-points Astronomy & Astrophysics manuscript no. icvis8 c ESO 28 June 8, 28 Anisotropic viscous dissipation in compressible magnetic X-points I. J. D. Craig Department of Mathematics, University of Waikato, P.

More information

Scaling of Magnetic Reconnection in Collisional and Kinetic Regimes

Scaling of Magnetic Reconnection in Collisional and Kinetic Regimes Scaling of Magnetic Reconnection in Collisional and Kinetic Regimes William Daughton Los Alamos National Laboratory Collaborators: Vadim Roytershteyn, Brian Albright H. Karimabadi, Lin Yin & Kevin Bowers

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

Plasma heating in stellarators at the fundamental ion cyclotron frequency

Plasma heating in stellarators at the fundamental ion cyclotron frequency PHYSICS OF PLASMAS VOLUME 7, NUMBER FEBRUARY 000 Plasma heating in stellarators at the fundamental ion cyclotron frequency V. A. Svidzinski and D. G. Swanson Department of Physics, Auburn University, Auburn,

More information

arxiv: v1 [physics.plasm-ph] 15 Mar 2011

arxiv: v1 [physics.plasm-ph] 15 Mar 2011 Gyromap for a two-dimensional Hamiltonian fluid model derived from raginskii s closure for magnetized plasmas arxiv:1103900v1 physicsplasm-ph 15 Mar 011 O Izacard, 1 C Chandre, 1 E Tassi, 1 and G Ciraolo

More information

Finite-Orbit-Width Effect and the Radial Electric Field in Neoclassical Transport Phenomena

Finite-Orbit-Width Effect and the Radial Electric Field in Neoclassical Transport Phenomena 1 TH/P2-18 Finite-Orbit-Width Effect and the Radial Electric Field in Neoclassical Transport Phenomena S. Satake 1), M. Okamoto 1), N. Nakajima 1), H. Sugama 1), M. Yokoyama 1), and C. D. Beidler 2) 1)

More information

Francesco Califano. Physics Department, University of Pisa. The role of the magnetic field in the interaction of the solar wind with a magnetosphere

Francesco Califano. Physics Department, University of Pisa. The role of the magnetic field in the interaction of the solar wind with a magnetosphere Francesco Califano Physics Department, University of Pisa The role of the magnetic field in the interaction of the solar wind with a magnetosphere Collaboration with M. Faganello & F. Pegoraro Vien na,

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

Collisionless nonideal ballooning modes

Collisionless nonideal ballooning modes PHYSICS OF PLASMAS VOLUME 6, NUMBER 1 JANUARY 1999 Collisionless nonideal ballooning modes Robert G. Kleva and Parvez N. Guzdar Institute for Plasma Research, University of Maryland, College Park, Maryland

More information

Two-fluid magnetic island dynamics in slab geometry. II. Islands interacting with resistive walls or resonant magnetic perturbations

Two-fluid magnetic island dynamics in slab geometry. II. Islands interacting with resistive walls or resonant magnetic perturbations PHYSICS OF PLASMAS 1, 0308 005 To-fluid magnetic isl dynamics in slab geometry. II. Isls interacting ith resistive alls or resonant magnetic perturbations Richard Fitzpatrick François L. Waelbroeck Institute

More information

MHD turbulence in the solar corona and solar wind

MHD turbulence in the solar corona and solar wind MHD turbulence in the solar corona and solar wind Pablo Dmitruk Departamento de Física, FCEN, Universidad de Buenos Aires Motivations The role of MHD turbulence in several phenomena in space and solar

More information

Magnetic reconnection in high-lundquist-number plasmas. N. F. Loureiro Instituto de Plasmas e Fusão Nuclear, IST, Lisbon, Portugal

Magnetic reconnection in high-lundquist-number plasmas. N. F. Loureiro Instituto de Plasmas e Fusão Nuclear, IST, Lisbon, Portugal Magnetic reconnection in high-lundquist-number plasmas N. F. Loureiro Instituto de Plasmas e Fusão Nuclear, IST, Lisbon, Portugal Collaborators: R. Samtaney, A. A. Schekochihin, D. A. Uzdensky 53 rd APS

More information

Answers to Problem Set Number MIT (Fall 2005).

Answers to Problem Set Number MIT (Fall 2005). Answers to Problem Set Number 6. 18.305 MIT (Fall 2005). D. Margetis and R. Rosales (MIT, Math. Dept., Cambridge, MA 02139). December 12, 2005. Course TA: Nikos Savva, MIT, Dept. of Mathematics, Cambridge,

More information

Two-fluid theory of collisionless magnetic reconnection

Two-fluid theory of collisionless magnetic reconnection Two-fluid theory of collisionless magnetic reconnection D. Biskamp and E. Schwarz Max-Planck-Institut für Plasmaphysik, 85748 Garching, Germany J. F. Drake Institute for Plasma Research, University of

More information

Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models

Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models 0-0 Applying Asymptotic Approximations to the Full Two-Fluid Plasma System to Study Reduced Fluid Models B. Srinivasan, U. Shumlak Aerospace and Energetics Research Program, University of Washington, Seattle,

More information

Exact solutions for magnetic annihilation in curvilinear geometry

Exact solutions for magnetic annihilation in curvilinear geometry Exact solutions for magnetic annihilation in curvilinear geometry E. Tassi b,, V.S. Titov and G. Hornig Theoretische Physik IV, Ruhr-Universität Bochum, 44780 Bochum, Germany b Theoretische Physik IV,

More information

Bifurcated states of a rotating tokamak plasma in the presence of a static error-field

Bifurcated states of a rotating tokamak plasma in the presence of a static error-field Bifurcated states of a rotating tokamak plasma in the presence of a static error-field Citation: Physics of Plasmas (1994-present) 5, 3325 (1998); doi: 10.1063/1.873000 View online: http://dx.doi.org/10.1063/1.873000

More information

Exponential Growth of Nonlinear Ballooning Instability. Abstract

Exponential Growth of Nonlinear Ballooning Instability. Abstract Exponential Growth of Nonlinear Ballooning Instability P. Zhu, C. C. Hegna, and C. R. Sovinec Center for Plasma Theory and Computation University of Wisconsin-Madison Madison, WI 53706, USA Abstract Recent

More information

ブラックホール磁気圏での 磁気リコネクションの数値計算 熊本大学 小出眞路 RKKコンピュー 森野了悟 ターサービス(株) BHmag2012,名古屋大学,

ブラックホール磁気圏での 磁気リコネクションの数値計算 熊本大学 小出眞路 RKKコンピュー 森野了悟 ターサービス(株) BHmag2012,名古屋大学, RKK ( ) BHmag2012,, 2012.2.29 Outline Motivation and basis: Magnetic reconnection around astrophysical black holes Standard equations of resistive GRMHD Test calculations of resistive GRMHD A simulation

More information

13. REDUCED MHD. Since the magnetic field is almost uniform and uni-directional, the field has one almost uniform component ( B z

13. REDUCED MHD. Since the magnetic field is almost uniform and uni-directional, the field has one almost uniform component ( B z 13. REDUCED MHD One often encounters situations in which the magnetic field is strong and almost unidirectional. Since a constant field does not produce a current density, these fields are sometimes said

More information

Influence of Generalized (r, q) Distribution Function on Electrostatic Waves

Influence of Generalized (r, q) Distribution Function on Electrostatic Waves Commun. Theor. Phys. (Beijing, China) 45 (2006) pp. 550 554 c International Academic Publishers Vol. 45, No. 3, March 15, 2006 Influence of Generalized (r, q) istribution Function on Electrostatic Waves

More information

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

Plasma Astrophysics Chapter 1: Basic Concepts of Plasma. Yosuke Mizuno Institute of Astronomy National Tsing-Hua University

Plasma Astrophysics Chapter 1: Basic Concepts of Plasma. Yosuke Mizuno Institute of Astronomy National Tsing-Hua University Plasma Astrophysics Chapter 1: Basic Concepts of Plasma Yosuke Mizuno Institute of Astronomy National Tsing-Hua University What is a Plasma? A plasma is a quasi-neutral gas consisting of positive and negative

More information

SW103: Lecture 2. Magnetohydrodynamics and MHD models

SW103: Lecture 2. Magnetohydrodynamics and MHD models SW103: Lecture 2 Magnetohydrodynamics and MHD models Scale sizes in the Solar Terrestrial System: or why we use MagnetoHydroDynamics Sun-Earth distance = 1 Astronomical Unit (AU) 200 R Sun 20,000 R E 1

More information

Kelvin-Helmholtz and Rayleigh-Taylor instabilities in partially ionised prominences

Kelvin-Helmholtz and Rayleigh-Taylor instabilities in partially ionised prominences Highlights of Spanish Astrophysics VII, Proceedings of the X Scientific Meeting of the Spanish Astronomical Society held on July 9-13, 2012, in Valencia, Spain. J. C. Guirado, L.M. Lara, V. Quilis, and

More information

Analysis of Jeans Instability of Partially-Ionized. Molecular Cloud under Influence of Radiative. Effect and Electron Inertia

Analysis of Jeans Instability of Partially-Ionized. Molecular Cloud under Influence of Radiative. Effect and Electron Inertia Adv. Studies Theor. Phys., Vol. 5, 2011, no. 16, 755-764 Analysis of Jeans Instability of Partially-Ionized Molecular Cloud under Influence of Radiative Effect and Electron Inertia B. K. Dangarh Department

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

RESISTIVE BALLOONING MODES AND THE SECOND REGION OF STABILITY

RESISTIVE BALLOONING MODES AND THE SECOND REGION OF STABILITY Plasma Physics and Controlled Fusion, Vol. 29, No. 6, pp. 719 to 121, 1987 Printed in Great Britain 0741-3335/87$3.00+.OO 1OP Publishing Ltd. and Pergamon Journals Ltd. RESISTIVE BALLOONING MODES AND THE

More information

Solar Flare. A solar flare is a sudden brightening of solar atmosphere (photosphere, chromosphere and corona)

Solar Flare. A solar flare is a sudden brightening of solar atmosphere (photosphere, chromosphere and corona) Solar Flares Solar Flare A solar flare is a sudden brightening of solar atmosphere (photosphere, chromosphere and corona) Flares release 1027-1032 ergs energy in tens of minutes. (Note: one H-bomb: 10

More information

ELECTROSTATIC ION-CYCLOTRON WAVES DRIVEN BY PARALLEL VELOCITY SHEAR

ELECTROSTATIC ION-CYCLOTRON WAVES DRIVEN BY PARALLEL VELOCITY SHEAR 1 ELECTROSTATIC ION-CYCLOTRON WAVES DRIVEN BY PARALLEL VELOCITY SHEAR R. L. Merlino Department of Physics and Astronomy University of Iowa Iowa City, IA 52242 December 21, 2001 ABSTRACT Using a fluid treatment,

More information

Ion plateau transport near the tokamak magnetic axis. Abstract

Ion plateau transport near the tokamak magnetic axis. Abstract Ion plateau transport near the tokamak magnetic axis K.C. Shaing and R.D. Hazeltine Institute for Fusion Studies, The University of Texas at Austin, Austin, Texas 78712 (December 2, 1997) Abstract Conventional

More information

Plasma instabilities. Dr Ben Dudson, University of York 1 / 37

Plasma instabilities. Dr Ben Dudson, University of York 1 / 37 Plasma instabilities Dr Ben Dudson, University of York 1 / 37 Previously... Plasma configurations and equilibrium Linear machines, and Stellarators Ideal MHD and the Grad-Shafranov equation Collisional

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

Large scale flows and coherent structure phenomena in flute turbulence

Large scale flows and coherent structure phenomena in flute turbulence Large scale flows and coherent structure phenomena in flute turbulence I. Sandberg 1, Zh. N. Andrushcheno, V. P. Pavleno 1 National Technical University of Athens, Association Euratom Hellenic Republic,

More information

Asymmetry and thermal effects due to parallel motion of electrons in collisionless magnetic reconnection

Asymmetry and thermal effects due to parallel motion of electrons in collisionless magnetic reconnection Asymmetry and thermal effects due to parallel motion of electrons in collisionless magnetic reconnection Motohiko Tanaka a) Research Laboratory of Electronics and Department of Physics, Massachusetts Institute

More information

Cold plasma waves. Waves in non-magnetized plasma Cold plasma dispersion equation Cold plasma wave modes

Cold plasma waves. Waves in non-magnetized plasma Cold plasma dispersion equation Cold plasma wave modes Cold plasma waves Waves in non-magnetized plasma Cold plasma dispersion equation Cold plasma wave modes EM wave propagation through and interaction with plasmas belong to central issues of plasma physics.

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

Fast magnetic reconnection via jets and current micro-sheets

Fast magnetic reconnection via jets and current micro-sheets Fast magnetic reconnection via jets and current micro-sheets P. G. Watson Center for Magnetic Reconnection Studies, Institute for Fusion Studies, Department of Physics, University of Texas at Austin, Austin,

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

arxiv: v1 [physics.plasm-ph] 23 Jan 2019

arxiv: v1 [physics.plasm-ph] 23 Jan 2019 arxiv:9.777v [physics.plasm-ph] 3 Jan 9 Gridless particle technique for the Vlasov-Poisson system in problems with high degree of symmetry E. Boella G. Coppa A. D Angola 3 and B. Peiretti Paradisi 4 Center

More information

Heat Transport in a Stochastic Magnetic Field. John Sarff Physics Dept, UW-Madison

Heat Transport in a Stochastic Magnetic Field. John Sarff Physics Dept, UW-Madison Heat Transport in a Stochastic Magnetic Field John Sarff Physics Dept, UW-Madison CMPD & CMSO Winter School UCLA Jan 5-10, 2009 Magnetic perturbations can destroy the nested-surface topology desired for

More information

Control of Neo-classical tearing mode (NTM) in advanced scenarios

Control of Neo-classical tearing mode (NTM) in advanced scenarios FIRST CHENGDU THEORY FESTIVAL Control of Neo-classical tearing mode (NTM) in advanced scenarios Zheng-Xiong Wang Dalian University of Technology (DLUT) Dalian, China Chengdu, China, 28 Aug, 2018 Outline

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

A simple model of the resistive wall mode in tokamaks

A simple model of the resistive wall mode in tokamaks A simple model of the resistive wall mode in tokamaks Richard Fitzpatrick Institute for Fusion Studies, Department of Physics, University of Texas at Austin, Austin TX 78712 (February 18, 2003) A simple

More information

Week 5-6: Lectures The Charged Scalar Field

Week 5-6: Lectures The Charged Scalar Field Notes for Phys. 610, 2011. These summaries are meant to be informal, and are subject to revision, elaboration and correction. They will be based on material covered in class, but may differ from it by

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

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

Modelling of the penetration process of externally applied helical magnetic perturbation of the DED on the TEXTOR tokamak

Modelling of the penetration process of externally applied helical magnetic perturbation of the DED on the TEXTOR tokamak INSTITUTE OF PHYSICS PUBLISHING Plasma Phys. Control. Fusion 8 (6) 69 8 PLASMA PHYSICS AND CONTROLLED FUSION doi:.88/7-/8// Modelling of the penetration process of externally applied helical magnetic perturbation

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

INTERACTION OF AN EXTERNAL ROTATING MAGNETIC FIELD WITH THE PLASMA TEARING MODE SURROUNDED BY A RESISTIVE WALL

INTERACTION OF AN EXTERNAL ROTATING MAGNETIC FIELD WITH THE PLASMA TEARING MODE SURROUNDED BY A RESISTIVE WALL INTERACTION OF AN EXTERNAL ROTATING MAGNETIC FIELD WITH THE PLASMA TEARING MODE SURROUNDED BY A RESISTIVE WALL S.C. GUO* and M.S. CHU GENERAL ATOMICS *Permanent Address: Consorzio RFX, Padova, Italy **The

More information

Plasma Flow in MST: Effects of Edge Biasing and Momentum Transport from Nonlinear Magnetic Torques

Plasma Flow in MST: Effects of Edge Biasing and Momentum Transport from Nonlinear Magnetic Torques Plasma Flow in MST: Effects of Edge Biasing and Momentum Transport from Nonlinear Magnetic Torques J.S. Sarff, A.F. Almagri, J.K. Anderson, B.E. Chapman, D. Craig, C-S. Chiang, N.A. Crocker, D.J. Den Hartog,

More information

Stabilization of a low-frequency instability inadipoleplasma

Stabilization of a low-frequency instability inadipoleplasma J. Plasma Physics: page 1 of 8. c 2008 Cambridge University Press doi:10.1017/s0022377808007071 1 Stabilization of a low-frequency instability inadipoleplasma D.T. GARNIER 1,A.C.BOXER 2, J.L. ELLSWORTH

More information

Introduction to Magnetohydrodynamics (MHD)

Introduction to Magnetohydrodynamics (MHD) Introduction to Magnetohydrodynamics (MHD) Tony Arber University of Warwick 4th SOLARNET Summer School on Solar MHD and Reconnection Aim Derivation of MHD equations from conservation laws Quasi-neutrality

More information

Theory and Simulation of Neoclassical Transport Processes, with Local Trapping

Theory and Simulation of Neoclassical Transport Processes, with Local Trapping Theory and Simulation of Neoclassical Transport Processes, with Local Trapping Daniel H. E. Dubin Department of Physics, University of California at San Diego, La Jolla, CA USA 92093-0319 Abstract. Neoclassical

More information

Plasma Physics for Astrophysics

Plasma Physics for Astrophysics - ' ' * ' Plasma Physics for Astrophysics RUSSELL M. KULSRUD PRINCETON UNIVERSITY E;RESS '. ' PRINCETON AND OXFORD,, ', V. List of Figures Foreword by John N. Bahcall Preface Chapter 1. Introduction 1

More information

Magnetospheric Physics - Final Exam - Solutions 05/07/2008

Magnetospheric Physics - Final Exam - Solutions 05/07/2008 Magnetospheric Physics - Final Exam - Solutions 5/7/8. Dipole magnetic field a Assume the magnetic field of the Earth to be dipolar. Consider a flux tube with a small quadratic cross section in the equatorial

More information

Cut-on, cut-off transition of sound in slowly varying flow ducts

Cut-on, cut-off transition of sound in slowly varying flow ducts Cut-on, cut-off transition of sound in slowly varying flow ducts Sjoerd W. Rienstra 19-3-1 Department of Mathematics and Computing Science, Eindhoven University of Technology, The Netherlands, s.w.rienstra@tue.nl

More information

Relativistic magnetohydrodynamics. Abstract

Relativistic magnetohydrodynamics. Abstract Relativistic magnetohydrodynamics R. D. Hazeltine and S. M. Mahajan Institute for Fusion Studies, The University of Texas, Austin, Texas 78712 (October 19, 2000) Abstract The lowest-order description of

More information

Weibel instability and filamentary structures of a relativistic electron beam in plasma

Weibel instability and filamentary structures of a relativistic electron beam in plasma Mitglied der Helmholtz-Gemeinschaft 7 th Direct Drive and Fast Ignition Workshop Prague, 3-6 May, 009 Weibel instability and filamentary structures of a relativistic electron beam in plasma Anupam Karmakar

More information

Magnetic field reconnection is said to involve an ion diffusion region surrounding an

Magnetic field reconnection is said to involve an ion diffusion region surrounding an The magnetic field reconnection site and dissipation region by P.L. Pritchett 1 and F.S. Mozer 2 1. Department of Physics and Astronomy, UCLA, Los Angeles, CA 90095-1547 2. Space Sciences Laboratory, University

More information

The role of polarization current in magnetic island evolution

The role of polarization current in magnetic island evolution The role of polarization current in magnetic island evolution J. W. Connor, F. L. Waelbroeck, and H. R. Wilson Citation: Phys. Plasmas 8, 835 (001); doi: 10.1063/1.137006 View online: http://dx.doi.org/10.1063/1.137006

More information

0.2. CONSERVATION LAW FOR FLUID 9

0.2. CONSERVATION LAW FOR FLUID 9 0.2. CONSERVATION LAW FOR FLUID 9 Consider x-component of Eq. (26), we have D(ρu) + ρu( v) dv t = ρg x dv t S pi ds, (27) where ρg x is the x-component of the bodily force, and the surface integral is

More information