Speaker
Description
See the full abstract here:
http://ocs.ciemat.es/EPS2019ABS/pdf/P1.1065.pdf
We construct fluid plasma models that combine physical consistency and numerical stability with a simple and flexible implementation [1]. This is achieved by exploiting the anti-symmetric property of the plasma flow operator, ∫dVφ(∇·v+v·∇)ψ=−∫dVψ(∇·v+v·∇)φ, which results in exact conservation. Since this expression has an equivalent discrete analog, computer implementations of fluid models using our representation automatically inherit the conservation properties of the continuous system. Additionally, it can be shown that the plasma velocity generates infinitesimal rotations of the fluid element, implying local time reversibility, e.g. existence of a discrete inverse. The anti-symmetric equations are written using generalized moments √n,√nv,√p, related to conserved quantities n,K= (1/2)mv2,U= (3/2)p, n, nv, p, related to conserved quantities n, K=
(1/2)mv2,
U
=
(3/2) p,
that ensure mass and energy positivity. Altogether, these properties lead to an intrinsically stable code implementation, provided that a few restrictions on the choice of numerical scheme are satisfied. We demonstrate this by applying our methodology to the Braginskii two-fluid model, the magnetohydrodynamics (MHD) equations, and also more computationally efficient drift-ordered models. The conservation properties are verified using representative cases such as single seeded blob dynamics and the Orzsag-Tang vortex, obtaining high-fidelity simulation results with negligible dissipation.
This material is based upon work supported by the U.S. Department of Energy, Office of Science, Office of Fusion Energy Sciences, Theory Program, under Award no. DE-FG02-95ER54309.
References
[1] F.D. Halpern and R.E. Waltz, Phys. Plasmas 25, 060703 (2018); F.D. Halpern and R.E. Waltz, A family of consistent and numerically stable fluid plasma models (to be submitted).