Related papers: Exponential Integrators for Resistive Magnetohydro…
This paper is concerned with geometric exponential energy-preserving integrators for solving charged-particle dynamics in a magnetic field from normal to strong regimes. We firstly formulate the scheme of the methods for the system in a…
This paper introduces Exp-ParaDiag, a novel time-parallel method that combines the strength of exponential integrators into the ParaDiag framework. We develop and analyze Exp-ParaDiag based on first and second order accurate exponential…
In this paper, we investigate the use of higher-order exponential Rosenbrock time integration methods on the shallow water equations on the sphere. This stiff, nonlinear model provides a testing ground for accurate and stable time…
We present a MATLAB toolbox for five different classes of exponential integrators for solving (mildly) stiff ordinary differential equations or time-dependent partial differential equations. For the efficiency of such exponential…
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…
Exponential integrators are time stepping schemes which exactly solve the linear part of a semilinear ODE system. This class of schemes requires the approxima- tion of a matrix exponential in every step, and one successful modern method is…
Among the family of fourth-order time integration schemes, the two-stage Gauss--Legendre method, which is an implicit Runge--Kutta method based on collocation, is the only superconvergent. The computational cost of this implicit scheme for…
We propose and analyse a numerical integrator that computes a low-rank approximation to large time-dependent matrices that are either given explicitly via their increments or are the unknown solution to a matrix differential equation.…
In this paper, we develop a new adaptive hyperbolic-cross-space mapped Jacobi (AHMJ) method for solving multidimensional spatiotemporal integrodifferential equations in unbounded domains. By devising adaptive techniques for sparse mapped…
We propose a computationally efficient and systematically convergent approach for elastodynamics simulations. We recast the second-order dynamical equation of elastodynamics into an equivalent first-order system of coupled equations, so as…
In this work, we describe various improved implementations of the mapping approach to surface hopping (MASH) for simulating nonadiabatic dynamics. These include time-reversible and piecewise-continuous integrators, which is only formally…
The electron inertia term and the off-diagonal electron pressure terms are well-known for the frozen-in condition breakdown in collisionless magnetic reconnection, which are naturally kinetic and difficult to be employed in…
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…
Seismic imaging is a major challenge in geophysics with broad applications. It involves solving wave propagation equations with absorbing boundary conditions (ABC) multiple times. This drives the need for accurate and efficient numerical…
In certain astrophysical systems the commonly employed ideal magnetohydrodynamics (MHD) approximation breaks down. Here, we introduce novel explicit and implicit numerical schemes of ohmic resistivity terms in the moving-mesh code AREPO. We…
We present a novel parallel implementation for large-scale three-dimensional electromagnetic inversion based on a Gauss-Newton framework combined with a rational near-best approximation of the matrix exponential for transient simulations.…
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…
We propose a high order adaptive-rank implicit integrators for stiff time-dependent PDEs, leveraging extended Krylov subspaces to efficiently and adaptively populate low-rank solution bases. This allows for the accurate representation of…
A variational integrator for ideal magnetohydrodynamics is derived by applying a discrete action principle to a formal Lagrangian. Discrete exterior calculus is used for the discretisation of the field variables in order to preserve their…