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 (10)
- Astrophysical magnetic fields and nonlinear dynamo theory
- Direct numerical simulation of turbulent channel flow up to
- 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
- An Intercomparison Between Divergence-Cleaning and Staggered Mesh Formulations for Numerical Magnetohydrodynamics
- Direct simulations of helical Hall-MHD turbulence and dynamo action
- Fourier continuation method for incompressible fluids with boundaries
- Adaptive mesh refinement with spectral accuracy for magnetohydrodynamics in two space dimensions
- Instability in electromagnetically driven flows Part II
- Magnetic structure, dipole reversals, and 1/f noise in resistive MHD spherical dynamos