A spectral element discretization for quasi‐static magnetohydrodynamic flows
Mattias Brynjell‐Rahkola
DOI: https://doi.org/10.1002/fld.5321
IF: 1.847
2024-07-20
International Journal for Numerical Methods in Fluids
Abstract:The staggered Legendre spectral element method is extended to the quasi‐static magnetohydrodynamic (MHD) setting, which relative to hydrodynamics involves an additional Poisson‐like problem for the electric potential. The proposed discretization is suitable for large‐scale simulations involving three‐dimensional geometries of arbitrary complexity that feature turbulent and/or transitional dynamics. The numerical solutions are shown to converge exponentially with polynomial order. The classical staggered PN‐PN−2 spectral element method (SEM) is revisited and extended to quasi‐static magnetohydrodynamic (MHD) flows. In this realm, which is valid in the limit of vanishing magnetic Reynolds number, the evaluation of the Lorentz force in the momentum equation requires the electric current density, governed by Ohm's law and a charge conservation condition derived from Ampère's law, to be determined. Once discretized with the SEM, this translates into solving one additional problem for the electric potential involving the so‐called consistent Poisson operator. The method is well suited for fully three‐dimensional flows in complex geometries. Changes in resolution requirements aside, consideration of the electromagnetic quantities is estimated to increase the computational cost associated with MHD by about 40% relative to hydrodynamics. The accuracy and the capabilities of the scheme is demonstrated on a set of common flows from the MHD literature. Exponential convergence with polynomial order is confirmed for the electric current density.
mechanics,physics, fluids & plasmas,computer science, interdisciplinary applications,mathematics