HR: 1340h
AN: SM13A-1118    [Abstracts]
TI: Plasma Equilibrium With Entropy Input and Applications
AU: Birn, J
EM: jbirn@lanl.gov
AF: Los Alamos National Laboratory, MS D466, P.O. Box 1663, Los Alamos, NM 87545, United States
AU: * Zaharia, S G
EM: szaharia@lanl.gov
AF: Los Alamos National Laboratory, MS D466, P.O. Box 1663, Los Alamos, NM 87545, United States
AB: During slow, quasi-static evolution in the collisionless magnetotail the entropy (S) as a function of magnetic flux (A) is a much better conserved quantity than the plasma pressure function P(A). The entropy profile S(A) has been found to be quasi-conserved even in scenarios involving topology change through magnetic reconnection, as long as the dissipation is very localized. Even though entropy conservation is a better constraint for space plasmas, most available plasma equilibrium solvers take the pressure profile P(A) as an input. This is mostly due to the difficulty of computing equilibria with S(A) input: as opposed to prescribing the current or pressure profiles for plasma equilibrium (which leads to the conventional Grad-Shafranov elliptic partial differential equation), prescribing the entropy results in a nonstandard differential equation whose solution by usual iterative techniques is not guaranteed. To solve this nonstandard equation we have developed a numerical code based on the 'alternating dimensions method' (ADM) suggested by Grad et al., [1975]; the code alternates between solving the 2D Grad-Shafranov equation to obtain the field configuration (flux function A) from the pressure profile P(A) and a 1D ordinary differential equation that uses the entropy profile to derive the pressure function P(A) from a flux surface average. We apply this method to 2 cases relevant to the magnetotail, both involving topological change in the transition from the initial to the final state: first, the code is applied to a forced magnetic reconnection problem (the 'Newton Challenge'), and the resulting equilibria obtained with various boundary deformations are compared with snapshots from dynamic MHD simulations. The second application starts from a more realistic magnetotail geometry that includes magnetic field line flaring. Assuming entropy conservation, we compute possible final equilibria with reconnected topology, and contrast them to equilibria obtained with unchanged P(A) profile. We also discuss how much magnetic energy is converted into thermal energy through topological change in both scenarios.
DE: 2744 Magnetotail
DE: 2764 Plasma sheet
DE: 7827 Kinetic and MHD theory
DE: 7833 Mathematical and numerical techniques (0500, 3200)
DE: 7835 Magnetic reconnection (2723, 7526)
SC: SPA-Magnetospheric Physics [SM]
MN: 2007 Fall Meeting