Self-consistent mean field forces in two-fluid models of turbulent plasmas

1. Self-consistent mean field forces in two-fluid models of turbulent plasmas C. C. Hegna University of Wisconsin Madison, WI Hall Dynamo Get-together PPPL via UW-Madison June 11, 2004

2. Theses • The properties of turbulent plasmas are described using the two-fluid equations. • Global constraints are derived for the fluctuation induced mean field forces that act on the ion and electron fluids. • Relationship between relaxation of parallel momentum flows and parallel currents C. C. Hegna, “Self-consistent mean-field forces in turbulent plasmas: current and momentum relaxation,” Physics of Plasmas5, 2257 (1998); 3480 (1998). --- RFP physics was largely the motivation

3. Outline • Brief review of mean field resistive MHD theory relevant to magnetized plasmas - applications to RFPs • Two-fluid theory • Constraints on the fluctuation induced mean-field forces • Heuristic derivations of local forms for the mean-field forces • A simple quasilinear theory Subsequent work - Mirnov, et al ‘03 - is much more complete • Relation to relaxation Steinhauer, Ishida, ‘98-’03; Mahajan and co-workers ‘01

4. In resistive MHD dynamo theory, a mean field force is identified • Fluctuations affects mean field dynamics in resistive MHD through a dynamo electric field Write all quantities as mean field and fluctuations The bracket <> notation denotes either an ensemble average or an average over the “small” spatial scales or “fast” time scales of the fluctuations Mean field Ohm’s Law

5. Global conservation laws have motivated local forms for the mean field force of resistive MHD • In resistive MHD, fluctuations do not dissipate helicity, but do dissipate energy. (Boozer, J. Plasma Physics, 1986; Bhattacharjee and Hameiri, PRL 1986; Phys. Fluids 1987; Strauss, Phys. Fluids 1985). These condititions are used to motivate a “local” form for the mean-field force in toroidal confinement devices --- fluctuations generate an additional electron viscosity or hyper-resitivity, not a “dynamo.” K2 is a profile dependent positive function satisfying boundary conditions. Consistent with the Taylor state, F gets large, ---> J||/B = constant

6. Two-fluid equations can be written in a concise form • The exact two-fluid momentum balance equations • These equations can be written more concisely with the identification of the canonical momentum. Momentum balance equations Plasma flow for each species

7. A pressure equation is also used for each species • The pressure evolution equations • Q = collision energy transfer and Ohmic heating, last term represents viscous heating. • In general, compressibility is allowed. This modifies the usual definition of the mean field force and allows for anomalous particle transport. • In what follows, the effects of heat flux q are simplified. A weakness in the theory and a potential new area of investigation.

8. Fluctuations induce mean field forces on both the ion and electron species. • For simplicity, we consider a cylindrical plasmas with all the usual boundary conditions. Quantities are split into equilibrium and fluctuating quantities, Nonlinearities produce fluctuation induced mean field forces (actually forces per unit charge) Note, the first term contains both the MHD and Hall dynamo terms. For the electrons, ve=u - J/ne + O(me/mi).

9. Three global properties of the mean-field forces can be shown • Mean field momentum balance equations • Three global constraints to be shown The last condition can also be written using F||M = F||i -F||e, F||O = (miF|||e+meF||i)/(mi+me)

10. A number of assumptions are used to prove the three global constraints • Simplifying assumptions used in the constraint derivations: • Fluctuation amplitudes are small compared to the mean magnetic field, typically valid in all MFE devices • The equilibrium quantities evolve on a slow diffusive time scale Viscosities and radial mean flow are ordered with resistivity. Parallel heat flux is ordered small to be consistent with the neglect of heat flux, (again, probably a weak point)

11. A number of assumptions are used to prove the three global constraints • Assumptions (continued) • The viscous force is dissipative for both species • All other equilibrium flows are ordered small - probably not a crucial assumption,may be generalized to equilibrium with flow • Ion and electron skin depths are small. With b ~ 1, rs ~ ds While velocity and magnetic field fluctuations are small, gradients of fluctuating quantities may be large, in general

12. The first two conditions can be shown from the generalized helicity evolution equations • Two separate ways to generate the evolution of the mean generalized helicity • The first from the total momentum balance • The last from the mean momentum balance • Subtracting the average of the first equation from the last equation.

13. With the assumed orderings and appropriate boundary conditions, the first two conditions are derived • The previously derived condition • All the terms on the right hand side are smaller than O(h) • C1 is the fluctuation induced generalized helicity flux. • In the me = 0 limit, the electron condition corresponds to the same as that derived for resistive MHD. In two-fluid theory, there two constraints, one for each fluid.

14. Energy balance relations are used to prove the third condition • Total energy conservation • Construct mean magnetic energy evolution from • Subtract this from O(h) average of the top equation • C2 is the leading-order energy flux caused by the fluctuations • Third term denotes anomalous cross-field transport --- similar bits show up in resistive MHD - Hameiri and Bhattacharjee, ‘87.

15. By accounting for the cross-field diffusion in our definition of F, the fluctuations are shown to dissipate energy • One can redefine the mean field force to account for turbulence induced cross field heat and particle transport • This redefinition doesn’t affect The first two conditions • The final condition is derived • This condition can also be written using FM = Fi - Fe, FO = (miFe + meFi)/(me + mi) • Fluctuations dissipate energy. Energy delivered to electrons and ions via Ohmic and viscous heating.

16. The global constraints imply a particular choice for the local form of the mean field forces • In analogy with the hyper-resistivity form implied by the resistive MHD global constraints, local forms are implied by the two-fluid global constraints. • The first two conditions imply where xs vanishes on the boundary • The third condition can then be written • A solution that guarantees the above is with conditions on the coefficients fst.

17. The implied local forms suggest a coupling between current and flow evolution • The parallel components of the turbulent mean field force • Can rewrite these equations as Those that appear in Ohm’s law And the total momentum balance Coefficients are spatially dependent functions that vanish on the boundary and satisfy ke2 > 0, ki2 > 0, (Le + Li)2 < ke2ki2/4 ke2 > 0, ki2 > 0, (Le + Li)2 < ke2ki2/4

18. The local forms are also derivable from quasi-linear tearing mode theory • Simplified, incompressible, low-b, me = 0, ps = 0 limit • Mean field forces are derived from quasilinear expression for the tearing mode resonant on some surface. < > = average over layer width • In the Ohm’s law, the resistive MHD and Hall dynamos • In the total momentum balance, the Reynolds and Maxwell stresses (to within a factor of ne).

19. A simple model is used in the linear layer analysis • A four field model • Equilibrium near the rational surface

20. The mean field forces derived from quasilinear theory exhibit the same structure as that implied from the global constraints • After linear theory • In the resistive MHD limit, Me2 is the largest coefficient - hyper-resistivity dominates and FM is small. • In the two-fluid limit, the dominant drive comes in the combination dVe/dx=dVo/dx - dJo/dx (ne)-1 • Highly simplified theory --- A much more complete job calculation of two-fluid tearing mode growth rates has been calculated by Mirnov, et al.

21. Relationship to relaxation theory • Constraints imply a relaxation theory via the minimization of W - leKe - liKi where Ks = ÷dV As.Bs • To lowest order in c/(wpia), the minimizing solutions yield A state discussed by many - Sudan, ‘79; Finn and Antonsen, ‘83; Avinash and Taylor ‘91; Steinhauer and Ishida, ‘98; Mahajan et al ‘01.

22. Summary • Using a modest number of assumptions, three global constraints are derived for turbulence induced mean field forces in two-fluid models of plasmas. • These constrains imply functional forms for the parallel mean-field forces in the Ohm’s law and the total momentum balance equations suggesting the fluctuations relax the plasma to states with field aligned current and bulk plasma momentum. • Applications to flow profile evolution during discrete dynamo events on MST?