Resistive MHD, reconnection and resistive tearing modes
|
|
- Veronica Pierce
- 5 years ago
- Views:
Transcription
1 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 Introduction In these notes we show how MagnetoHydroDynamics (MHD can be deduced from Braginskii equations. We then study how keeping some of the small Braginskii terms, such as resistivity, can be important. Resistivity and other Braginskii (or two-fluid effects break the frozen-in law that ties plasma and magnetic field lines together, allowing magnetic reconnection (breaks in magnetic field lines to happen. This new effect can lead to instabilities that cannot be described using MHD.. MagnetoHydroDynamics (MHD In MHD (Schekochihin 15 the plasma is modeled as a single fluid characterized by its density ρ(r, t, its average velocity u(r, t, its pressure p(r, t and its current density J(r, t. For a Braginskii plasma composed of one ion species of charge e and mass m i and electrons with charge e and mass m e, the plasma density is the plasma average velocity is the plasma total pressure is and the plasma current density is ρ = n i m i + n e m e = n e (m i + m e n e m i, (.1 u = n im i u i + n e m e u e n i m i + n e m e u i, (. p = n i T i + n e T e = n e (T i + T e, (.3 J = en e (u i u e. (.4 To obtain the MHD equations for ρ, u, p and J, we start from the Braginskii equations, derived using the orderings ρ i L λ ee L λ ei L λ ii L me 1, (.5 m i and u i u e u i u e v ti u i u e ρ i L v ti, (.6 T i T e T i T e, (.7 among other assumptions. To find the MHD equations, we will modify these assumptions
2 Felix I. Parra and perform a subsidiary expansion in ρ i L λ ee L λ ei L λ ii L me m i 1. (.8 In this limit, the parallel collisional friction force and the collisional temperature equilibration term, of order F ei, n e m e ν ei u i u e me m i L u i u e p e λ ei L v ti (.9 and W ie = 3n em e ν ei me L T e T i p e v ti (T e T i m i m i λ ei T e L, (.1 must be of the same order as the other terms in their respective equations, the parallel momentum equations and the energy equations of both ions and electrons. Thus, we assume that F ei, p e /L, giving mi λ ei u i u e v ti u i u e m e L v ti u i u e ρ i L v ti, (.11 and that W ie p e v ti /L, leading to T i T e T i T e mi m e λ ei L T e. (.1 These assumptions are different from Braginskii s more general ordering in (.6 and (.7. Equations (.11 and (.1 imply that and u e u i = u (.13 T i T e = T. (.14 The plasma is sufficiently collisional that both electrons and ions move at the same average velocity and have the same temperature, and can then be treated as a single fluid. We proceed to derive the MHD equations from Braginskii equations. Multiplying Braginskii s ion continuity equation by m i, we obtain the MHD continuity equation, ρ + (ρu =. (.15 t The equation for u is Braginskii s total momentum equation. Note that under the assumptions in (.5, we showed that the ion viscosity was small in λ ii /L 1, giving ( u ρ t + u u = p + J B, (.16 We also use the electron momentum equation. We have shown that under the assumptions in (.5, the electromagnetic forces dominate, giving en e (E + u e B =. Using (.13, this equation finally becomes E + u B =. (.17
3 Resistive MHD, reconnection and resistive tearing modes 3 Summing the electron and ion thermal energy equations and using (.13 and (.14, we find ( T 3n e t + u T = n e T u, (.18 where we have neglected the terms q i and Π i : u i because they are small in λ ii /L 1, and the terms q e and F ei (u i u e because they are small in m i /m e (λ ee /L 1. Using (.15 and p = n e (T i + T e n e T, equation (.18 can be rewritten as the usual MHD energy equation ( ( p t + u =. (.19 ρ 5/3 In addition to the plasma model in (.15, (.16, (.17 and (.19, we need Maxwell s equations for the evolution of E and B: the induction equation and Ampere s law without the displacement current, E = B t, (. B = µ J. (.1 3. Resistive MHD and reconnection We have neglected several terms to derive the MHD equations. Some of them are in fact important even though they are small. We will focus on the electron momentum equation. The electron momentum equation with next order terms is en e (E + u e B p e +.51n e m e ν ei (u i u e ˆb.71n eˆbˆb Te, (3.1 where we have neglected n e m e ( u e + u e u e, Π e and F ei, because they are much smaller than the electron pressure gradient. We rewrite this equation performing the following manipulations. According to (.11, u i u e u i u e, leading to (u i u e ˆb u i u e. Hence, we can write where is the plasma resistivity..51n e m e ν ei (u i u e ˆb.51n e m e ν ei (u i u e = en e ηj, (3. We can rewrite ˆbˆb T e as Using these results, equation (3.1 becomes where the velocity u B is η =.51m eν ei e n e (3.3 ˆbˆb T e = T e (ˆb T e ˆb. (3.4 E + u B B = p e en e + ηj.71 e T e, (3.5 u B = u e }{{} v ti.71 eb ˆb T e u. (3.6 }{{} ρ i L v ti
4 4 Felix I. Parra There are two differences between the usual MHD electron momentum equation, E + u B =, and equation (3.5: the difference between the plasma velocity u and the new velocity u B, and the terms in the right side of (3.5. Both of these differences are small, but some of the terms in the right side of (3.5 are important because they change a fundamental property of MHD: the frozen-in law. To see what happens to the frozen-in law, we obtain E from (3.5 and substitute it into the induction equation (. to find B t = (u B B n e p e en (ηj. (3.7 e The usual MHD frozen-in law is B/ t = (u B, and it implies that the magnetic field lines move with the plasma velocity u. Since plasma infinitesimal volume elements do not split into two (u would have to be discontinuous, a magnetic field line that moves with the plasma cannot be broken. Conversely, equation (3.7 implies that the magnetic field lines move with the velocity u B u, and in addition to this motion, the magnetic field changes in time due to two new terms: n e p e /en e and (ηj. In these notes, we focus on the term (ηj due to resistivity. The term n e p e /en e is in general as important as the resistivity term, but it can be made zero if we take the solution p/ρ 5/3 = constant to equation (.19, implying that p e n 5/3 e and hence n e p e =. Then, in the particular case p/ρ 5/3 = constant, and using (.1 to obtain J, equation (3.7 becomes B t = (u B B ( η B. (3.8 µ Using [(η/µ B] = (η/µ ( B+µ 1 B η [(η/µ B] and B =, we rewrite (3.8 as ( B η ignore t = (u B B + B 1µ B µ η. (3.9 Thus, if we ignore η to avoid unnecessary complications, the magnetic field moves with velocity u B, which is uninteresting because u B is almost the plasma velocity u, and it diffuses with a diffusion coefficient D η = η/µ. The fact that the magnetic field diffuses implies that magnetic field lines can be broken. The resistivity is not the only term that can break magnetic field lines, and there is a large area of research that studies the different physical effects that can lead to magnetic reconnection. From here on, we only consider resistivity, and instead of (3.5, we consider the simpler equation E + u B = ηj (3.1 that contains the necessary terms to obtain (3.9. To simplify the problem even further, we consider η a constant, i.e. we consider a plasma in which the temperature variation is small. We characterize the strength of the resistive term using the characteristic time of the magnetic field diffusion τ η = L /D η = µ L /η, where L is the length of the system. We measure the importance of the resistivity in the plasma by comparing this characteristic time with the Alfven time, τ A = L/v A, where v A = B/ ρµ is the Alfven speed. This comparison gives the Lundquist number, S = τ η = µ v A L τ A η = BL µ η ρ. (3.11
5 Resistive MHD, reconnection and resistive tearing modes (a (b J z /J z.5 J z /J z B y /B y1 1 B y /B y x/l x/l Figure 1. Sketches of B y(x (solid and J z(x (dashed normalized by their values at x, B y = B y(x, and at x =, J z = J z(x =, respectively. Figure (a represents a general current sheet, and figure (b is the specific example in equations (4. and (4.5. Since we are assuming that resistivity is a small effect, we are considering the case S Resistive tearing modes To study the effect of the magnetic diffusion, we consider a current sheet, which is a magnetic field configuration that is stable to ideal MHD modes, but unstable to small but finite resistivity. The unstable mode reconnects the magnetic field, converting the magnetic energy into plasma energy. In this section, we first describe the equilibrium that we will consider, and we then study its stability by linearizing the equations Current sheet equilibrium We consider a simple sheared magnetic field B = B z ẑ + B y (xŷ. (4.1 The function B y (x is of the form shown in figure 1. To simplify the calculation, we assume in these notes that B y (x is antisymmetric with respect to x =, i.e. B y (x = B y ( x, and hence B y = at x =. We also assume that B y tends to a constant for x ±. A possible example, shown in figure 1(b, is { By x/l for x L B y (x = B y for x > L. (4. Assuming that B y (x is antisymmetric with respect to x = is particularly restrictive and not necessary for the derivation, but this assumption simplifies the calculation around x =. The sheared magnetic field in (4.1 is similar to the magnetic field used in the seminal calculation by Furth, Killeen and Rosenbluth (FKR that we follow in these notes (Furth et al Using Ampere s law, given in (.1, we find the current density corresponding to the magnetic field in (4.1, J = 1 µ B = J z (xẑ, (4.3 where J z (x = 1 µ B y(x. (4.4
6 6 Felix I. Parra From here on, the prime indicates differentiation with respect to x. The dependence of J on x is sketched in figure 1. The localization of the current is the reason why this type of equilibria is called a current sheet. For the magnetic field in (4., we find { By /µ J z (x = L for x L for x > L. (4.5 This current profile is sketched in figure 1(b. The current must be carried by the plasma. For the plasma equilibrium, we assume that u =. Then, equation (.16 gives p = J B = (B y B y/µ ˆx, and we can integrate in x to obtain the equilibrium pressure p(x. Equation (3.1 gives E = E z (xẑ, with E z (x = ηj z (x. This electric field is externally imposed to drive the current in the current sheet. Note that the plasma density ρ(x is undetermined. We assume that it only depends on x. The equilibrium that we have just described is stable in the MHD model because the sheared magnetic field makes moving the magnetic field lines impossible (they cannot break and cross each other. Including resistivity into the analysis will change the stability properties. 4.. Perturbation equations Linearizing equations (.16, (3.1, (. and (.1, we find ρ u 1 t = p 1 + J 1 B + J B 1, (4.6 E 1 + u 1 B = ηj 1, (4.7 and E 1 = B 1 t (4.8 B 1 = µ J 1. (4.9 We will not need the continuity equation or the energy equation. Due to the symmetry of the equilibrium (it only depends on x and B z is a constant, we can choose a form of the perturbation that simplifies the derivation. We will assume that B 1 (x, y, t, J 1 (x, y, t, E 1 (x, y, t, u 1 (x, y, t and p 1 (x, y, t depend on x and y, but not on z. Moreover, we will assume that B 1 ẑ =. Since the perturbed magnetic field must satisfy B 1 =, we choose B 1 = ẑ ψ, (4.1 where ψ(x, y, t is the flux function. With this form of B 1 and (4.9, we obtain the perturbed current density J 1 = 1 B 1 = ψ ẑ. (4.11 µ µ The form of B 1 in (4.1 gives the form of the electric field E 1 and the velocity u 1. Equation (4.8 gives E 1 = ψ ẑ φ, (4.1 t where φ(x, y, t is the electrostatic potential. Using (4.1, the components of the vector equation (4.7 perpendicular to ẑ give u 1 = 1 B z ẑ φ, (4.13
7 Resistive MHD, reconnection and resistive tearing modes 7 where we have assumed u 1 ẑ = to simplify the derivation. We need equations for ψ and φ. The ẑ-component of (4.7 gives ψ t B y B z φ y = η µ ψ. (4.14 To find an equation for φ, we use (4.6. To eliminate the term p 1, we take a curl of this equation, obtaining the vorticity equation, t [ (ρu 1] = B J 1 J 1 B + B 1 J J B 1. (4.15 Using (4.1, (4.11 and (4.13, equation (4.15 gives ( µ ψ B z t [ (ρ φ] = B y y To simplify the equations even further, we take the forms ψ(x, y, t = ψ(x exp(γt + ik y y + c.c., B y ψ y. (4.16 φ(x, y, t = φ(x exp(γt + ik y y + c.c., (4.17 where γ is the growth rate of the instability and k y is the mode number in the direction of symmetry y. With this form of the solution, equations (4.14 and (4.16 become γ }{{} ψ ik y B y φ = η ( ψ k B }{{ z µ ψ y (4.18 }} {{} and ψ γτ A τ A φ v A ψ φ ψ γτ A v A ψ L ψ τ A S 1 ψ τ A µ γ [ (ρ ik y B z B φ ρk φ ] y = ψ k ψ y B y ψ. (4.19 y B y }{{}}{{} The order of magnitude estimates in these equations are obtained by assuming that k y L 1 and d/dx 1/L. The boundary conditions for equations (4.18 and (4.19 are that ψ, φ for x ±. For Lundquist number S 1, the mode that satisfies these boundary conditions and equations (4.18 and (4.19 is composed of two very distinct regions: an outer region in which the resistivity is negligible, and an inner region in which the resistivity is comparable to other terms. We solve the equations in these two regions and then match them to find the shape of the mode and the growth rate γ Outer region We will obtain that ψ L S 1 γτ A 1. (4. Then, the right side of equation (4.18 is small and can be neglected, leading to φ o = iγb z k y B y ψo, (4.1 where the subindex o indicates that this is only valid in the outer region. According to this result, the left side of (4.19 is of order (γτ A ψ/l ψ/l, and hence, equation
8 8 Felix I. Parra k y L small o x/l k y L large o x/l Figure. The functions ψ o(x and φ o(x in the outer region. For ψ o(x, we give typical solutions for k yl small (solid and k yl large (dashed. (4.19 becomes ψ o ky ψ o B y ψo =. (4. B y Note however, that the left side of (4.19 is small only away from x = because at x =, B y =. Thus, equation (4. is only valid for ], [ and ], [. We solve the equation with the boundary conditions ψ o for x ±, and ψ o ( = ψ o ( + (4.3 at x =. It is not obvious that ψ o should be continuous at x =, but we will see that this condition is consistent with the solution in the inner region. For x ±, B y /B y, and equation (4. becomes ψ o k y ψ o =. (4.4 Then, the solutions are ψ o exp(±k y x for x ±, and we must choose the solution that goes like exp(k y x for x <, and the solution that goes like exp( k y x for x > (here we have assumed without loss of generality that k y >. This choice fixes the derivative at x = and x = +. Thus, in general the solution is of the form shown in figure, with a discontinuous derivate at x =. Because we have assumed that B y (x is antisymmetric, the solutions are symmetric in x, giving ψ o (x ψ o ( (1 + x. (4.5 The parameter = ψ o( + ψ o(. (4.6 ψ o ( is the standard measure of the size of the discontinuity in ψ o. The function φ o (x, given by (4.1, diverges at x = as shown in figure due to B y (x B y(x for x close to zero. To find an example of outer region solution, we use the magnetic field in (4.. For this magnetic field, equation (4. becomes ψ o k y ψ o 1 L [δ(x + L δ(x L] ψ o =. (4.7
9 Resistive MHD, reconnection and resistive tearing modes L k y L Figure 3. The quantity as a function of k yl. We have used the particular case in (4.9. The solution to this equation is ψ o (x ψ o ( = cosh(k y x + 1 k yl k y L tanh(k y L k y L (1 k y L tanh(k y L sinh(k y x k y L k y L cosh(k y L (1 k y L sinh(k y L exp( k y( x L As expected, this solution has a discontinuous derivative at x =. We find for x L for x > L. (4.8 = k y[1 k y L k y L tanh(k y L] k y L (1 k y L tanh(k y L. (4.9 The quantity ranges from + at k y L 1 to for k y L 1, as can be seen in figure 3. In fact, the quantity ranges from to + when k y L decreases from to for any profile of B y (x that is of the form shown in figure 1 (see the different profiles of ψ o (x for different values of k y L in figure. We show that this is the case in Appendix A. The magnetic field is unknown in a region around x =. The discontinuity in ψ o indicates that near x =, the assumptions that we have used to obtain equations (4.1 and (4. fail. For this reason, we need to solve the equations in an inner region in which the resistivity becomes important Inner region In the outer region solution, the first derivative ψ o is discontinuous. This discontinuity is indicative of a fast rate of change of ψ i in the inner region, that is, ψ i is large (see figure 4. Here the subindex i refers to quantities in the inner region. If the size of the inner region is of order δ L, ψ i ψ δl. (4.3 The function φ i, given in (4.1 for the outer region, must be of the form shown in figure 4. Thus, we expect its second derivative to satisfy φ i φ δ. (4.31
10 1 Felix I. Parra i 3 i 1 / x/ 5 5 x/ Figure 4. The functions ψ(x and φ(x in the inner region for >. Using the estimates (4.3 and (4.31, and employing B y (x B y(x, B y (x B y (x and ρ(x ρ(, equations (4.18 and (4.19 become and γ }{{} ψ i ψ γτ A τ A ik y B y(x φi B }{{ z } φ δ ψ v A ψ L τ A = η µ ψ i } {{} S 1 L δ ψ τ A small by δ/l 1 η µ k ψ y i (4.3 } {{} S 1 ψ τ A µ γρ( φ ik y B z B y(x i small by (δ/l 1 µ γρ( ( ik y B z B y(x k φ y i = ψ i }{{} small by δ/l 1 ky + B y ( ψ B i. (4.33 y( }{{}}{{}}{{} φ γτ A v A ψ ( L δ 3 ψ L φ L ψ γτ A v A ψ δ L By balancing all the terms in equations (4.3 and (4.33, we find the order of magnitude of γτ A, δ/l and φ/v A ψ, γτ A S 3/5, L δ ψ L φ v A ψ S 1/5, ψ L δ L S /5. (4.34 Note that the result for γ is consistent with our assumption (4., needed to solve the outer region. According to the estimates in (4.34, δ S /5 L, and hence, since ψ i ψ/l, the total change in ψ i across the inner region is ψ i δ ψ i δ L ψ S /5 ψ ψ. (4.35 Thus, to lowest order in S /5 1, ψi (x ψ o (, where ψ o ( is determined by the outer region. The fact that ψ(x does not change much across the inner region justifies the boundary condition (4.3 for the outer region. The approximation ψ i (x ψ o ( is known as constant flux approximation. Using ψ i (x ψ o (, equation (4.3 becomes Using equation (4.33 to obtain γ ψ o ( ik yb y(x φi = η B z µ ψ i. (4.36 ψ i as a function of φ i, we can rewrite equation (4.36
11 Resistive MHD, reconnection and resistive tearing modes 11 Im( 1/4 Re( 1/4 Choice of 1/4 as Figure 5. Sketch showing our choice of γ 1/4. γ ψ o ( ik yb y(x φi = iηγρ( φ B z k y B z B y(x i. (4.37 For large x, the solution for φ i must match with the outer region solution. The outer region solution, given in (4.1, is φ o iγb z k y B y(x ψ o ( (4.38 for x around zero. Thus, the boundary condition for equation (4.37 is that φ i iγb z k y B y(x ψ o ( (4.39 for x ±. Once φ i (x is obtained from equation (4.37 with boundary condition (4.39, we can find ψ i (x by using (4.33, ψ i = iµ γρ( φ k y B z B y(x i. (4.4 We can integrate this equation across the inner region to find the size of the jump in ψ i, ψ i( ψ i( = ψ i dx = iµ γρ( k y B z B y( φ i x dx. (4.41 This jump in ψ i must be equal to the discontinuity observed in the outer region solution, characterized by the parameter. Then, iµ γρ( = k y B z B y( ψ o ( This condition is an equation for the growth rate γ Final solution To solve the equation for φ i in (4.37, we use the normalized variables φ i x dx. (4.4 U = φ i Φ, X = x δ, (4.43 where the quantities Φ and δ need to be determined. With this normalization, equation (4.37 becomes γ ψ o ( ik yb y(δφ XU = iηγρ(φ 1 d U B z k y B z B y(δ 3 X dx. (4.44
12 1 Felix I. Parra Imposing we obtain the normalization factors and the equation for U(X, γ ψ o ( = ik yb y(δφ B z = iηγρ(φ k y B z B y(δ 3, (4.45 δ = (ηρ(1/4 γ 1/4 (k y B y( 1/, Φ = iγ 3/4 B z (k y B y( 1/ (ηρ( 1/4 ψ o (, (4.46 d U = X(1 + XU. (4.47 dx With the results in (4.46, the boundary condition in (4.39 becomes U(X 1 X (4.48 for X ±. In (4.46, there are four possible choices for γ 1/4 : γ 1/4 exp(iθ γ /4 + ikπ/, with k =, 1,, 3. Here θ γ is the argument of γ, defined to be between π and π. To be consistent with the solution for U(X that we will obtain in (4.49, we choose the root of γ 1/4 = γ 1/4 exp(iθ γ /4 with positive real part and the imaginary part with smallest absolute value, that is, the root with argument between π/4 and π/4 (see figure 5. The solution to equation (4.47 with boundary conditions (4.48 is the Rutherford- Furth solution U(X = X 1 exp ( µx (1 µ 1/4 dµ. (4.49 (see Appendix B for a proof. This solution satisfies the boundary condition if X has positive real part, Re(X >, because then we can change variable of integration to s = µx / to find that for X large U(X = 1 X X / ( exp ( s 1 small when exp( s 1 1/4 4s X 4 ds 1 X / exp ( s ds X [ = 1 small for Re(X > 1 exp ( ] X X 1 X. (4.5 We note that X = x /δ, and that the argument of δ is the argument of γ 1/4. Thus, to ensure that Re(X >, we choose the root γ 1/4 = γ 1/4 exp(iθ γ /4 as in figure 5. Note that the root γ 1/4 exp(iθ γ /4 would also give Re(X >, but using this root would be equivalent to a simple change of sign of x, and it would give the same final φ i (x. With (4.46 and (4.49, the condition in (4.4 becomes = µ (ρ( 1/4 γ 5/4 (k y B y( 1/ η 3/4 1 d U X dx dx = µ ρ 1/4 γ 5/4 [Γ(3/4]. (4.51 (k y B y( 1/ η 3/4 where Γ(ν = t ν 1 exp( t dt is the Euler Gamma function (in Appendix C we show how to take the integral of the Rutherford-Furth solution. According to (4.51, is proportional to (γ 1/4 5. From our choice γ 1/4 = γ 1/4 exp(iθ γ /4 (see figure 5, the argument of must be between 5π/4 and 5π/4. Thus, positive corresponds to
13 Resistive MHD, reconnection and resistive tearing modes u 1 y/l 1 B + B 1 1 u Figure 6. Sketch of the perturbed magnetic field B + B 1 and perturbed velocity u 1 in the inner region. Note the different scale of x and y. exp(i, and negative to exp(±iπ. Solving for γ, we obtain ( 4/5 1 (k y B γ = y( /5 η 3/5 ( 4/5. (4.5 [Γ(3/4] µ 4/5 (ρ( 1/5 From this formula, we see that the current sheet is unstable for > ( k y L small since γ exp(i, and it is stable for < ( k y L large since γ exp(±4iπ/5. The perturbed magnetic field B + B 1 and the perturbed velocity u 1, defined in (4.1 and (4.13, are sketched in figure 6 for the unstable case. For k y L small, the current sheet breaks into a series of individual current channels. The drive for the instability is the fact that currents of the same sign attract. Resistivity is necessary because without it, the magnetic field cannot reconnect and form islands. In the case k y L large, in which <, the instability is quenched by the tension of the lines that prevents reconnection from happening by making it hard to bend magnetic field lines. x/ REFERENCES Furth, H.P., Killeen, J. & Rosenbluth, M.N Finite-resistivity instabilities of a sheet pinch. Phys. Fluids 6, 459. Schekochihin, A.A. 15 Lecture notes on Kinetic Theory and Magnetohydrodynamics of Plasmas.
14 14 Felix I. Parra Appendix A. Limits of ψ o for k y L 1 and k y L 1 In this Appendix, we show that scales as 1/k y for k y L 1, and as k y for k y L 1, hence proving that ranges from to as k y decreases from to.. For k y L 1, we solve equation (4. by trying solutions of the form ψ o = Ψ exp( k y x. Then, equation (4. becomes [ B y(x ( ] Ψ = x B y (x x k yb y (xψ. (A 1 The function Ψ must be continuous at x = and it has to become constant for x ± (recall that ψ o exp( k y x for large x. Expanding Ψ in k y L 1, we find Ψ(x = Ψ (x + Ψ 1 (x +... (A }{{} k ylψ Ψ To lowest order, equation (A 1 becomes [By(Ψ /B y ] =. One possible solution to this equation that satisfies the boundary conditions is Ψ (x = B y ( x. It will not become obvious why we choose this symmetric form until we solve the next order of equation (A 1, given by [ ( ] By(x Ψ1 = x B y (x x k yb y (xψ = k y B y (xb y(x. (A 3 We can integrate this equation imposing first that Ψ 1 for x ± to find ( ( Ψ1 B y = k y B y (x By(x 1, (A 4 where B y = B y (x is the value of B y at infinity. We integrate equation (A 4 again imposing that Ψ 1 for x ± to obtain ( B y Ψ 1 (x = k y B y ( x By(s 1 ds. (A 5 x (we could have chosen that Ψ 1 tend to a constant different from zero, but this constant can be absorbed into Ψ. The solution (A 5 is continuous at x = because we chose the lowest order symmetric function Ψ = B y ( x. Adding Ψ and Ψ 1, we find that ψ o is given by [ ( ] B y ψ o (x = B y ( x exp( k y x 1 + k y By(s 1 ds + O(kyL. (A 6 x The solution in (A 6 gives = [B y(] /k y By + for k y. For k y L 1, equation (4. becomes ψ o ky ψ o =. The solution to this equation that goes to zero at x ± and is continuous at x = is ψ o = exp( k y x. (A 7 This solution leads to = k y. We will see shortly that the mode is unstable when >, and stable when <. Thus, the mode will be unstable for small k y, and stable for large k y.
15 Resistive MHD, reconnection and resistive tearing modes 15 Appendix B. Rutherford-Furth solution In this Appendix we check by direct substitution that the function U(X in (4.49 is a solution to equation (4.47. By differentiating with respect to X twice, we find d U dx = 3X 1 1 exp ( µx µ(1 µ 1/4 dµ X3 exp ( µx µ (1 µ 1/4 dµ. (B 1 The first integral can be rewritten in a convenient form by integrating by parts in µ, 3X 1 [ ] µ=1 exp ( µx µ(1 µ 1/4 dµ = X exp ( µx (1 µ 3/4 µ= 1 X3 exp ( µx (1 µ 3/4 dµ 1 = X X3 exp ( µx (1 µ (1 µ 1/4 dµ. (B Using this result in (B 1, we find d U 1 dx = X X3 exp ( µx (1 µ 1/4 dµ = X(1 XU, (B 3 proving that the function U(X in (4.49 is a solution to equation (4.47. Appendix C. Integral of the Rutherford-Furth solution to obtain In this Appendix, we take the integral X 1 (d U/dX dx. Using (B 1, we find 1 d U 1 ( X dx dx = dx dµ exp ( µx 3 µx µ(1 µ 1/4. (C 1 We can take the integral in X first. Using the usual Gaussian integrals, we obtain 1 d U X dx dx = π 1 Changing to the new variable t = µ, the integral becomes 1 d U π 1 X dx dx = t 1/4 (1 t 1/4 dt = µ 1/ (1 µ 1/4 dµ. (C where B(ν, µ = 1 tν 1 (1 t µ 1 dt is the Euler Beta function. Since we finally obtain π B(3/4, 3/4, (C 3 B(ν, µ = Γ(νΓ(µ, (C 4 Γ(ν + µ 1 d U π X dx dx = [Γ(3/4] = [Γ(3/4]. (C 5 Γ(3/
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 information35. 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 informationSMR/ 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 informationFundamentals 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 informationGyrokinetic 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 informationMacroscopic 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 informationWave propagation in an inhomogeneous plasma
DRAFT Wave propagation in an inhomogeneous plasma Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford OX NP, UK This version is of 7 February 208. Introduction In
More informationMAGNETOHYDRODYNAMICS
Chapter 6 MAGNETOHYDRODYNAMICS 6.1 Introduction Magnetohydrodynamics is a branch of plasma physics dealing with dc or low frequency effects in fully ionized magnetized plasma. In this chapter we will study
More informationSW103: 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 informationCollisionless magnetic reconnection with arbitrary guide-field
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
More informationHeating 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 informationIntroduction 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 informationIssues 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 informationCurrent-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 informationFluid 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 information20. 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 informationRadiative & Magnetohydrodynamic Shocks
Chapter 4 Radiative & Magnetohydrodynamic Shocks I have been dealing, so far, with non-radiative shocks. Since, as we have seen, a shock raises the density and temperature of the gas, it is quite likely,
More informationFormation 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 informationTwo-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! e x2 erfi(x)!!=!! 1 2 i! e x2!erf(ix)!,
Solution to Problem 1 (a) Since ρ is constant and uniform, the continuity equation reduces to!iu!!=!! "u x "x!+! "u y "y!!=!!!.! which is satisfied when the expression for u x and u y are substituted into
More informationSimple 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 informationTwo 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 informationMHD RELATED TO 2-FLUID THEORY, KINETIC THEORY AND MAGANETIC RECONNECTION
MHD RELATED TO 2-FLUID THEORY, KINETIC THEORY AND MAGANETIC RECONNECTION Marty Goldman University of Colorado Spring 2017 Physics 5150 Issues 2 How is MHD related to 2-fluid theory Level of MHD depends
More informationKinetic, Fluid & MHD Theories
Lecture 2 Kinetic, Fluid & MHD Theories The Vlasov equations are introduced as a starting point for both kinetic theory and fluid theory in a plasma. The equations of fluid theory are derived by taking
More informationWaves 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 informationRecapitulation: 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 informationMHD 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 informationFokker-Planck collision operator
DRAFT 1 Fokker-Planck collision operator Felix I. Parra Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford OX1 3NP, UK (This version is of 16 April 18) 1. Introduction In these
More informationAPPENDIX Z. USEFUL FORMULAS 1. Appendix Z. Useful Formulas. DRAFT 13:41 June 30, 2006 c J.D Callen, Fundamentals of Plasma Physics
APPENDIX Z. USEFUL FORMULAS 1 Appendix Z Useful Formulas APPENDIX Z. USEFUL FORMULAS 2 Key Vector Relations A B = B A, A B = B A, A A = 0, A B C) = A B) C A B C) = B A C) C A B), bac-cab rule A B) C D)
More informationApplying 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 information14. Energy transport.
Phys780: Plasma Physics Lecture 14. Energy transport. 1 14. Energy transport. Chapman-Enskog theory. ([8], p.51-75) We derive macroscopic properties of plasma by calculating moments of the kinetic equation
More informationThe Particle in a Box
Page 324 Lecture 17: Relation of Particle in a Box Eigenstates to Position and Momentum Eigenstates General Considerations on Bound States and Quantization Continuity Equation for Probability Date Given:
More informationMAE210C: Fluid Mechanics III Spring Quarter sgls/mae210c 2013/ Solution II
MAE210C: Fluid Mechanics III Spring Quarter 2013 http://web.eng.ucsd.edu/ sgls/mae210c 2013/ Solution II D 4.1 The equations are exactly the same as before, with the difference that the pressure in the
More informationThe 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 informationSolar 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 informationStabilization 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 informationSupersonic 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 information12. MHD Approximation.
Phys780: Plasma Physics Lecture 12. MHD approximation. 1 12. MHD Approximation. ([3], p. 169-183) The kinetic equation for the distribution function f( v, r, t) provides the most complete and universal
More informationPlasma spectroscopy when there is magnetic reconnection associated with Rayleigh-Taylor instability in the Caltech spheromak jet experiment
Plasma spectroscopy when there is magnetic reconnection associated with Rayleigh-Taylor instability in the Caltech spheromak jet experiment KB Chai Korea Atomic Energy Research Institute/Caltech Paul M.
More information9. Boundary layers. Flow around an arbitrarily-shaped bluff body. Inner flow (strong viscous effects produce vorticity) BL separates
9. Boundary layers Flow around an arbitrarily-shaped bluff body Inner flow (strong viscous effects produce vorticity) BL separates Wake region (vorticity, small viscosity) Boundary layer (BL) Outer flow
More informationV (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 informationUNIVERSITY OF EAST ANGLIA
UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2007 2008 FLUIDS DYNAMICS WITH ADVANCED TOPICS Time allowed: 3 hours Attempt question ONE and FOUR other questions. Candidates must
More informationReduced 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 informationMagnetic reconnection in coronal plasmas
UW, 28 May, 2010 p.1/17 Magnetic reconnection in coronal plasmas I.J.D Craig Department of Mathematics University of Waikato Hamilton New Zealand UW, 28 May, 2010 p.2/17 Why reconnection? Reconnection
More informationChapter 5. Instabilities. 5.1 Maco-Instabilities - Ideal MHD Rayleigh-Taylor Instability
Chapter 5 Instabilities 5.1 Maco-Instabilities - Ideal MHD 5.1.1 Rayleigh-Taylor Instability Before we consider the so-caled normal mode approach let us first look at the result using th evariational approach:
More informationTransition 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 information3 Generation and diffusion of vorticity
Version date: March 22, 21 1 3 Generation and diffusion of vorticity 3.1 The vorticity equation We start from Navier Stokes: u t + u u = 1 ρ p + ν 2 u 1) where we have not included a term describing a
More informationDispersive 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 informationMagnetohydrodynamic 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 informationChapter 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 informationINTERACTION 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 informationMomentum transport from magnetic reconnection in laboratory an. plasmas. Fatima Ebrahimi
Momentum transport from magnetic reconnection in laboratory and astrophysical plasmas Space Science Center - University of New Hampshire collaborators : V. Mirnov, S. Prager, D. Schnack, C. Sovinec Center
More informationMeasuring 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 informationKinetic theory of ions in the magnetic presheath
Kinetic theory of ions in the magnetic presheath Alessandro Geraldini 1,2, Felix I. Parra 1,2, Fulvio Militello 2 1. Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford, Oxford
More informationWeek 2 Notes, Math 865, Tanveer
Week 2 Notes, Math 865, Tanveer 1. Incompressible constant density equations in different forms Recall we derived the Navier-Stokes equation for incompressible constant density, i.e. homogeneous flows:
More informationMagnetic 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 informationHeat 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 informationA 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 informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis
More informationReconnection 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 informationMicrotearing Simulations in the Madison Symmetric Torus
Microtearing Simulations in the Madison Symmetric Torus D. Carmody, P.W. Terry, M.J. Pueschel - University of Wisconsin - Madison dcarmody@wisc.edu APS DPP 22 Overview PPCD discharges in MST have lower
More information2/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 informationChapter 5 MAGNETIZED PLASMAS. 5.1 Introduction. 5.2 Diamagnetic current
Chapter 5 MAGNETIZED PLASMAS 5.1 Introduction We are now in a position to study the behaviour of plasma in a magnetic field. In the first instance we will re-examine particle diffusion and mobility with
More informationComparison 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(a) What are the probabilities associated with finding the different allowed values of the z-component of the spin after time T?
1. Quantum Mechanics (Fall 2002) A Stern-Gerlach apparatus is adjusted so that the z-component of the spin of an electron (spin-1/2) transmitted through it is /2. A uniform magnetic field in the x-direction
More informationブラックホール磁気圏での 磁気リコネクションの数値計算 熊本大学 小出眞路 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 informationAmplification 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 informationPROBLEM SET. Heliophysics Summer School. July, 2013
PROBLEM SET Heliophysics Summer School July, 2013 Problem Set for Shocks and Particle Acceleration There is probably only time to attempt one or two of these questions. In the tutorial session discussion
More informationEffect of magnetic reconnection on CT penetration into magnetized plasmas
Earth Planets Space, 53, 547 55, 200 Effect of magnetic reconnection on CT penetration into magnetized plasmas Yoshio Suzuki, Takaya Hayashi 2, and Yasuaki Kishimoto Naka Fusion esearch Establishment,
More informationIdeal 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 informationMagnetic Reconnection in Laboratory, Astrophysical, and Space Plasmas
Magnetic Reconnection in Laboratory, Astrophysical, and Space Plasmas Nick Murphy Harvard-Smithsonian Center for Astrophysics namurphy@cfa.harvard.edu http://www.cfa.harvard.edu/ namurphy/ November 18,
More informationCreation 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 informationThe RFP: Plasma Confinement with a Reversed Twist
The RFP: Plasma Confinement with a Reversed Twist JOHN SARFF Department of Physics University of Wisconsin-Madison Invited Tutorial 1997 Meeting APS DPP Pittsburgh Nov. 19, 1997 A tutorial on the Reversed
More informationSawteeth in Tokamaks and their relation to other Two-Fluid Reconnection Phenomena
Sawteeth in Tokamaks and their relation to other Two-Fluid Reconnection Phenomena S. C. Jardin 1, N. Ferraro 2, J. Chen 1, et al 1 Princeton Plasma Physics Laboratory 2 General Atomics Supported by the
More informationQuasi-separatrix layers and 3D reconnection diagnostics for linetied
Quasi-separatrix layers and 3D reconnection diagnostics for linetied tearing modes John M. Finn, LANL A. Steve Richardson, NRL CMSO, Oct 2011 Operated by Los Alamos National Security, LLC for the U.S.
More informationMAGNETIC 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 information6.1. Linearized Wave Equations in a Uniform Isotropic MHD Plasma. = 0 into Ohm s law yields E 0
Chapter 6. Linear Waves in the MHD Plasma 85 Chapter 6. Linear Waves in the MHD Plasma Topics or concepts to learn in Chapter 6:. Linearize the MHD equations. The eigen-mode solutions of the MHD waves
More informationEuler Equations: local existence
Euler Equations: local existence Mat 529, Lesson 2. 1 Active scalars formulation We start with a lemma. Lemma 1. Assume that w is a magnetization variable, i.e. t w + u w + ( u) w = 0. If u = Pw then u
More informationPLASMA 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 informationControl of linear modes in cylindrical resistive MHD with a resistive wall, plasma rotation, and complex gain
Control of linear modes in cylindrical resistive MHD with a resistive wall, plasma rotation, and complex gain Dylan Brennan 1 and John Finn 2 contributions from Andrew Cole 3 1 Princeton University / PPPL
More informationAccurate representation of velocity space using truncated Hermite expansions.
Accurate representation of velocity space using truncated Hermite expansions. Joseph Parker Oxford Centre for Collaborative Applied Mathematics Mathematical Institute, University of Oxford Wolfgang Pauli
More informationTransport by convection. Coupling convection-diffusion
Transport by convection. Coupling convection-diffusion 24 mars 2017 1 When can we neglect diffusion? When the Peclet number is not very small we cannot ignore the convection term in the transport equation.
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
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ρ c (2.1) = 0 (2.3) B = 0. (2.4) E + B
Chapter 2 Basic Plasma Properties 2.1 First Principles 2.1.1 Maxwell s Equations In general magnetic and electric fields are determined by Maxwell s equations, corresponding boundary conditions and the
More informationMHD 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 Turbulence, magnetic reconnection, particle acceleration Understand the mechanisms
More information7 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 informationRandom Walk on the Surface of the Sun
Random Walk on the Surface of the Sun Chung-Sang Ng Geophysical Institute, University of Alaska Fairbanks UAF Physics Journal Club September 10, 2010 Collaborators/Acknowledgements Amitava Bhattacharjee,
More informationLaminar Flow. Chapter ZERO PRESSURE GRADIENT
Chapter 2 Laminar Flow 2.1 ZERO PRESSRE GRADIENT Problem 2.1.1 Consider a uniform flow of velocity over a flat plate of length L of a fluid of kinematic viscosity ν. Assume that the fluid is incompressible
More information13. 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 informationBeyond Ideal MHD. Nick Murphy. Harvard-Smithsonian Center for Astrophysics. Astronomy 253: Plasma Astrophysics. February 8, 2016
Beyond Ideal MHD Nick Murphy Harvard-Smithsonian Center for Astrophysics Astronomy 253: Plasma Astrophysics February 8, 2016 These lecture notes are largely based on Plasma Physics for Astrophysics by
More informationNonlinear MRI. Jeremy Goodman. CMPD/CMSO Winter School 7-12 January 2008, UCLA
Nonlinear MRI Jeremy Goodman CMPD/CMSO Winter School 7-12 January 2008, UCLA Questions for the nonlinear regime How does the linear instability saturate? i.e., What nonlinear mode-mode interactions brake
More informationMagnetohydrodynamic 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[variable] = units (or dimension) of variable.
Dimensional Analysis Zoe Wyatt wyatt.zoe@gmail.com with help from Emanuel Malek Understanding units usually makes physics much easier to understand. It also gives a good method of checking if an answer
More informationMagnetic 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 informationBefore you begin read these instructions carefully.
MATHEMATICAL TRIPOS Part IB Thursday, 6 June, 2013 9:00 am to 12:00 pm PAPER 3 Before you begin read these instructions carefully. Each question in Section II carries twice the number of marks of each
More informationMHD 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 informationQuantum Physics III (8.06) Spring 2016 Assignment 3
Quantum Physics III (8.6) Spring 6 Assignment 3 Readings Griffiths Chapter 9 on time-dependent perturbation theory Shankar Chapter 8 Cohen-Tannoudji, Chapter XIII. Problem Set 3. Semi-classical approximation
More informationFinite Volume for Fusion Simulations
Finite Volume for Fusion Simulations Elise Estibals, Hervé Guillard, Afeintou Sangam To cite this version: Elise Estibals, Hervé Guillard, Afeintou Sangam. Finite Volume for Fusion Simulations. Jorek Meeting
More informationSOLAR 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 informationUNIVERSITY of LIMERICK
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4607 SEMESTER: Autumn 2012-13 MODULE TITLE: Introduction to Fluids DURATION OF
More information