Related papers: A semi-implicit Hall-MHD solver using whistler wav…
We propose an optimally performant fully implicit algorithm for the Hall magnetohydrodynamics (HMHD) equations based on multigrid-preconditioned Jacobian-free Newton-Krylov methods. HMHD is a challenging system to solve numerically because…
There is no standard numerical implementation of the Hall effect, which is one of the non-ideal magnetohydrodynamic (MHD) effects. Numerical instability arises when a simple implementation is used, in which the Hall electric field is added…
When modeling astrophysical fluid flows, it is often appropriate to discard the canonical magnetohydrodynamic approximation thereby freeing the magnetic field to diffuse with respect to the bulk velocity field. As a consequence, however,…
In this paper we present a novel pressure-based semi-implicit finite volume solver for the equations of compressible ideal, viscous and resistive magnetohydrodynamics (MHD). The new method is conservative for mass, momentum and total energy…
We consider the application of the WaveHoltz iteration to time-harmonic elastic wave equations with energy conserving boundary conditions. The original WaveHoltz iteration for acoustic Helmholtz problems is a fixed-point iteration that…
It is well established that the magnetically structured solar atmosphere supports the propagation of MHD waves along various kind of jets including also the solar wind. It is well-known as well that under some conditions, namely high enough…
We present a deep learning-based iterative approach to solve the discrete heterogeneous Helmholtz equation for high wavenumbers. Combining classical iterative multigrid solvers and convolutional neural networks (CNNs) via preconditioning,…
In this paper we combine a flexible covariant formulation of the shallow water equations with the semi-implicit numerical scheme developed over the years by Casulli and collaborators. After adopting an orthogonal, but non-orthonormal,…
The choice of numerical integrator in approximating solutions to dynamic partial differential equations depends on the smallest time-scale of the problem at hand. Large-scale deformations in elastic solids contain both shear waves and bulk…
The magnetohydrodynamics (MHD) equations are generally known to be difficult to solve numerically, due to their highly nonlinear structure and the strong coupling between the electromagnetic and hydrodynamic variables, especially for high…
In this work we introduce a novel semi-implicit structure-preserving finite-volume/finite-difference scheme for the viscous and resistive equations of magnetohydrodynamics (MHD) based on an appropriate 3-split of the governing PDE system,…
Hall-MHD is a mixed hyperbolic-parabolic partial differential equation that describes the dynamics of an ideal two fluid plasma with massless electrons. We study the only shock wave family that exists in this system (the other…
In subsurface multiphase flow simulations, poor nonlinear solver performance is a significant runtime sink. The system of fully implicit mass balance equations is highly nonlinear and often difficult to solve for the nonlinear solver,…
A numerical method for the quasi-neutral two-fluid (QNTF) plasma model is described. The basic equations are ion and electron fluid equations and the Maxwell equations without displacement current. The neglect of displacement current is…
We propose here a novel stabilization strategy for the PMCHWT equation that cures its frequency and conductivity related instabilities and is obtained by leveraging quasi-Helmholtz projectors. The resulting formulation is well-conditioned…
We present a preconditioning method for the multi-dimensional Helmholtz equation with smoothly varying coefficient. The method is based on a frame of functions, that approximately separates components associated with different singular…
We propose and analyze a new method for the unsteady incompressible magnetohydrodynamics equations on convex domains with hybrid approximations of both vector-valued and scalar-valued fields. The proposed method is convection-semirobust,…
We propose a second-order accurate semi-implicit and well-balanced finite volume scheme for the equations of ideal magnetohydrodynamics (MHD) including gravitational source terms. The scheme treats all terms associated with the acoustic…
A robust multilevel preconditioner based on the hybridizable discontinuous Galerkin method for the Helmholtz equation with high wave number is presented in this paper. There are two keys in our algorithm, one is how to choose a suitable…
New implicit and implicit-explicit time-stepping methods for the wave equation in second-order form are described with application to two and three-dimensional problems discretized on overset grids. The implicit schemes are single step,…