Related papers: Orbital advection with magnetohydrodynamics and ve…
Magnetic helicity fluxes are investigated in a family of gauges in which the contribution from ideal magnetohydrodynamics takes the form of a purely advective flux. Numerical simulations of magnetohydrodynamic turbulence in this advective…
A conservative scheme has been formulated and verified for gyrokinetic particle simulations of electromagnetic waves and instabilities in magnetized plasmas. An electron continuity equation derived from drift kinetic equation is used to…
We present the theory, algorithms and implementation of a parallel finite-volume algorithm for the solution of the incompressible magnetohydrodynamic (MHD) equations using unstructured grids that are applicable for a wide variety of…
This paper presents applications of weighted meshless scheme for conservation laws to the Euler equations and the equations of ideal magnetohydrodynamics. The divergence constraint of the latter is maintained to the truncation error by a…
A new implementation for the time evolution of the magnetic vector potential is obtained for smoothed particle magnetohydrodynamics by considering the induction equation in integral form. Galilean invariance is achieved through proper gauge…
Orbital magnetic susceptibility involves rich physics such as interband effects despite of its conceptual simplicity. In order to appreciate the rich physics related to the orbital magnetic susceptibility, it is essential to derive a…
High order algorithms have emerged in numerical astrophysics as a promising avenue to reduce truncation error (proportional to a power of the linear resolution $\Delta x$) with only a moderate increase to computational expense. Significant…
We present a new numerical scheme for solving the advection equation and its application to Vlasov simulations. The scheme treats not only point values of a profile but also its zeroth to second order piecewise moments as dependent…
We propose and analyze two convection quasi-robust and pressure robust finite element methods for a fully nonlinear time-dependent magnetohydrodynamics problem. Both methods employ the $H_{\rm div}$ conforming BDM element coupled with an…
We present a rigorous derivation of the orbital magnetization for interacting electrons within the self-consistent Hartree-Fock approximation. Our method also allows us to derive formulas for the orbital magnetic susceptibility. The results…
Accretion disks are three-dimensional, turbulent, often self-gravitating, magnetohydrodynamic flows, which can be modeled in detail with numerical simulations. In this paper, we present a new algorithm that is based on a spectral…
Numerical methods for solving the ideal magnetohydrodynamic (MHD) equations in more than one space dimension must either confront the challenge of controlling errors in the discrete divergence of the magnetic field, or else be faced with…
This paper proposes a virtual element method (VEM) combined with a second-order implicit-explicit scheme based on the scalar auxiliary variable (SAV) method for the incompressible magnetohydrodynamics (MHD) equations. We employ the BDF2…
We present a robust computational framework for advective-diffusive-reactive systems that satisfies maximum principles, the non-negative constraint, and element-wise species balance property. The proposed methodology is valid on general…
We develop a new numerical scheme for ideal magnetohydrodynamic (MHD) simulations, which is robust against one- and multi-dimensional shocks, and is accurate for low Mach number flows and discontinuities. The scheme belongs to a family of…
Long simulation times in climate sciences typically require coarse grids due to computational constraints. Nonetheless, unresolved subscale information significantly influences the prognostic variables and can not be neglected for reliable…
We present a novel approach to solve the advection-diffusion equation under arbitrary transporting fields using a quantum-inspired 'Schrodingerisation' technique for Hamiltonian simulation. Although numerous methods exist for solving…
We present a constrained transport (CT) algorithm for solving the 3D ideal magnetohydrodynamic (MHD) equations on a moving mesh, which maintains the divergence-free condition on the magnetic field to machine-precision. Our CT scheme uses an…
Recent three-dimensional radiative hydrodynamics simulations of protoplanetary disks report disparate disk behaviors, and these differences involve the importance of convection to disk cooling, the dependence of disk cooling on metallicity,…
Grid-based magnetohydrodynamic (MHD) simulations have proven invaluable for the study of astrophysical accretion disks. However, the fact that angular momentum transport in disks is mediated by MHD turbulence (with structure down to very…