Stabilization of numerical interchange in spectral-element magnetohydrodynamics
- Univ. of Wisconsin, Madison, WI (United States). Department of Engineering Physics
In this study, auxiliary numerical projections of the divergence of flow velocity and vorticity parallel to magnetic field are developed and tested for the purpose of suppressing unphysical interchange instability in magnetohydrodynamic simulations. The numerical instability arises with equal-order C0 finite- and spectral-element expansions of the flow velocity, magnetic field, and pressure and is sensitive to behavior at the limit of resolution. The auxiliary projections are motivated by physical field-line bending, and coercive responses to the projections are added to the flow-velocity equation. Their incomplete expansions are limited to the highest-order orthogonal polynomial in at least one coordinate of the spectral elements. Cylindrical eigenmode computations show that the projections induce convergence from the stable side with first-order ideal-MHD equations during h-refinement and p-refinement. Hyperbolic and parabolic projections and responses are compared, together with different methods for avoiding magnetic divergence error. Lastly, the projections are also shown to be effective in linear and nonlinear time-dependent computations with the NIMROD code [C. R. Sovinec, et al., J. Comput. Phys. 195 (2004) 355-386], provided that the projections introduce numerical dissipation.
- Research Organization:
- Univ. of Wisconsin, Madison, WI (United States)
- Sponsoring Organization:
- USDOE Office of Science (SC), Fusion Energy Sciences (FES)
- Grant/Contract Number:
- FC02-08ER54975; AC02-05CH11231
- OSTI ID:
- 1436975
- Alternate ID(s):
- OSTI ID: 1324850
- Journal Information:
- Journal of Computational Physics, Vol. 319, Issue C; ISSN 0021-9991
- Publisher:
- ElsevierCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Web of Science
Similar Records
Pressure-driven tearing and thermal transport in finite-beta reversed field pinch computations
Coupling extended magnetohydrodynamic fluid codes with radiofrequency ray tracing codes for fusion modeling