Vector potential-based MHD solver for non-periodic flows using Fourier continuation expansions
arXiv:2107.07077 · doi:10.1016/j.cpc.2022.108304
Abstract
A high-order method to evolve in time electromagnetic and velocity fields in conducting fluids with non-periodic boundaries is presented. The method has a small overhead compared with fast FFT-based pseudospectral methods in periodic domains. It uses the magnetic vector potential formulation for accurately enforcing the null divergence of the magnetic field, and allowing for different boundary conditions including perfectly conducting walls or vacuum surroundings, two cases relevant for many astrophysical, geophysical, and industrial flows. A spectral Fourier continuation method is used to accurately represent all fields and their spatial derivatives, allowing also for efficient solution of Poisson equations with different boundaries. A study of conducting flows at different Reynolds and Hartmann numbers, and with different boundary conditions, is presented to study convergence of the method and the accuracy of the solenoidal and boundary conditions.
References in corpus (4)
- Global-Scale Turbulent Convection and Magnetic Dynamo Action in the Solar Envelope
- A simulation of convective dynamo in the solar convective envelope: maintenance of the solar-like differential rotation and emerging flux
- Fourier continuation method for incompressible fluids with boundaries
- Adaptive mesh refinement with spectral accuracy for magnetohydrodynamics in two space dimensions