Related papers: Exponential Integrators for Resistive Magnetohydro…
Recently, an extended version of magnetohydrodynamics that incorporates electron inertia, dubbed inertial magnetohydrodynamics, has been proposed. This model features a noncanonical Hamiltonian formulation with a number of conserved…
This paper investigates the equations of motion for a relativistic charged particle in a general magnetic field. By reformulating the dynamics in four-dimensional spacetime and separating the linear and nonlinear parts, we construct an…
We construct a relativistic resistive magneto-hydrodynamic (RRMHD) numerical simulation code for high-energy heavy-ion collisions. We split the system of differential equations into two parts, a non-stiff and a stiff part. For the non-stiff…
We present a finite-volume, genuinely 4th-order accurate numerical method for solving the equations of resistive relativistic magnetohydrodynamics (Res-RMHD) in Cartesian coordinates. In our formulation, the magnetic field is evolved in…
In the framework of a mixed finite element method, a structure-preserving formulation for incompressible magnetohydrodynamic (MHD) equations with general boundary conditions is proposed. A leapfrog-type temporal scheme fully decouples the…
We present a new method for evolving the equations of magnetohydrodynamics (both Newtonian and relativistic) that is capable of maintaining a divergence-free magnetic field ($\nabla \cdot \mathbf{B} = 0$) on adaptively refined, conformally…
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…
We introduce exponential numerical integration methods for stiff stochastic dynamical systems of the form $d\mathbf{z}_t = L(t)\mathbf{z}_tdt + \mathbf{f}(t)dt + Q(t)d\mathbf{W}_t$. We consider the setting of time-varying operators $L(t),…
Exponential integrators are a well-known class of time integration methods that have been the subject of many studies and developments in the past two decades. Surprisingly, there have been limited efforts to analyze their stability and…
We present a multiscale integrator for Hamiltonian systems with slowly varying quadratic stiff potentials that uses coarse timesteps (analogous to what the impulse method uses for constant quadratic stiff potentials). This method is based…
The paper describes an explicit multi-dimensional numerical scheme for Special Relativistic Two-Fluid Magnetohydrodynamics of electron-positron plasma and a suit of test problems. The scheme utilizes Cartesian grid and the third order WENO…
We present an extension to the special relativistic, ideal magnetohydrodynamics (MHD) equations, designed to capture effects due to resistivity. The extension takes the simple form of an additional source term which, when implemented…
In this paper, we mainly focus on the rigorous convergence analysis of two fully decoupled, unconditionally energy-stable methods for the diffuse interface two-phase magnetohydrodynamics (MHD) model. The two methods consist of the…
We propose two new classes of time integrators for stiff DEs: the implicit-explicit exponential (IMEXP) and the hybrid exponential methods. In contrast to the existing exponential schemes, the new methods offer significant computational…
Simulation of geothermal systems is challenging due to coupled physical processes in highly heterogeneous media. Combining the exponential Rosenbrock--Euler and Rosenbrock-type methods with control-volume (two-point flux approximation)…
The paper describes a new upwind conservative numerical scheme for special relativistic resistive magnetohydrodynamics with scalar resistivity. The magnetic field is kept approximately divergence free and the divergence of the electric…
We present an efficient algorithm for maximum likelihood estimation (MLE) of exponential family models, with a general parametrization of the energy function that includes neural networks. We exploit the primal-dual view of the MLE with a…
The objective of this work is the introduction and investigation of favourable time integration methods for the Gross--Pitaevskii equation with rotation term. Employing a reformulation in rotating Lagrangian coordinates, the equation takes…
In this paper we propose a novel way to integrate time-evolving partial differential equations that contain nonlinear advection and stiff linear operators, combining exponential integration techniques and semi-Lagrangian methods. The…
Magnetic quadrupoles are essential components of particle accelerators like the Large Hadron Collider. In order to study numerically the stability of the particle beam crossing a quadrupole, a large number of particle revolutions in the…