Related papers: Exponential Integrators for Resistive Magnetohydro…
We propose a linearly implicit structure-preserving numerical method for semilinear Hamiltonian systems with polynomial nonlinearities, combining Kahan's method and exponential integrator. This approach efficiently balances computational…
A new code and methodology are introduced for solving the general relativistic magnetohydrodynamic (GRMHD) equations in fixed background spacetimes using time-explicit, finite-volume discretization. The code has options for solving the…
Maxwell equations describe the propagation of electromagnetic waves and are therefore fundamental to understanding many problems encountered in the study of antennas and electromagnetics. The aim of this paper is to propose and analyse an…
A version of extended magnetohydrodynamics (MHD) that incorporates electron inertia is obtained by constructing an action principle. Unlike MHD which freezes in magnetic flux, the present theory freezes in an alternative flux related to the…
We describe a method for incorporating ambipolar diffusion in the strong coupling approximation into a multidimensional magnetohydrodynamics code based on the total variation diminishing scheme. Contributions from ambipolar diffusion terms…
We propose several continuous data assimilation (downscaling) algorithms based on feedback control for the 2D magnetohydrodynamic (MHD) equations. We show that for sufficiently large choices of the control parameter and resolution and…
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 paper, we propose an efficient exponential integrator finite element method for solving a class of semilinear parabolic equations in rectangular domains. The proposed method first performs the spatial discretization of the model…
A Finite Volume (FV) scheme is developed for solving the extended magnetohydrodynamic (XMHD) equations, yielding accurate results in the ideal, resistive, and Hall MHD limits. This is accomplished by first re-writing the XMHD equations such…
An algorithm for a family of self-starting high-order implicit time integration schemes with controllable numerical dissipation is proposed for both linear and nonlinear transient problems. This work builds on the previous works of the…
We construct a family of solutions $(u,B)$ of the incompressible magnetohydrodynamic (MHD) system, the $L^\infty$ norm of which blows up instantaneously at the critical rate. The solutions remain smooth except at the blowup time. An inverse…
Three numerical algorithms are proposed to solve the time-dependent elastodynamic equations in elastic solids. All algorithms are based on approximating the solution of the equations, which can be written as a matrix exponential. By…
Perfectly matched layers are a very efficient and accurate way to absorb waves in media. We present a stable convolutional unsplit perfectly matched formulation designed for the linearized stratified Euler equations. However, the technique…
In this paper, we present exact divergence-free spectral method for solving the incompressible and resistive magneto-hydrodynamic (MHD) equations in two and three dimensions, as well as the efficient solution algorithm and unconditionally…
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,…
In this contribution we propose an exponential update algorithm for magnetic moments appearing in the framework of micromagnetics and the Landau-Lifshitz-Gilbert (LLG) equation. This algorithm can be interpreted as the geometric integration…
We present an algorithm for simulating the equations of ideal magnetohydrodynamics and other systems of differential equations on an unstructured set of points represented by sample particles. The particles move with the fluid, so the time…
Based on the algorithm Informed Importance Tempering (IIT) proposed by Li et al. (2023) we propose an algorithm that uses an adaptive bounded balancing function. We argue why implementing parallel tempering where each replica uses a…
The Smectic-A (SmA) phase is modeled by a modified Landau-de Gennes (mLdG) model proposed by Xia et al. [Phys. Rev. Lett., 126 (2021), 177801], in which a tensor order parameter Q for the orientational order is coupled with a real scalar…
Immersed boundary methods simplify mesh generation by embedding the domain of interest into an extended domain that is easy to mesh, introducing the challenge of dealing with cells that intersect the domain boundary. Combined with explicit…