Related papers: An exponential integrator for the drift-kinetic mo…
Exponential time differencing methods is a power tool for high-performance numerical simulation of computationally challenging problems in condensed matter physics, fluid dynamics, chemical and biological physics, where mathematical models…
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…
In this paper we consider an approach to improve the performance of exponential Runge--Kutta integrators and Lawson schemes} in cases where the solution of a related, but usually much simpler, problem can be computed efficiently. While for…
We introduce a class of exponential Runge-Kutta integration methods for kinetic equations. The methods are based on a decomposition of the collision operator into an equilibrium and a non equilibrium part and are exact for relaxation…
The anisotropic diffusion equation is imperative in understanding cosmic ray diffusion across the Galaxy, the heliosphere, and its interplay with the ambient magnetic field. This diffusion term contributes to the highly stiff nature of the…
The paper deals with numerical discretizations of separable nonlinear Hamiltonian systems with additive noise. For such problems, the expected value of the total energy, along the exact solution, drifts linearly with time. We present and…
To help reveal the complete picture of linear kinetic drift modes, four independent numerical approaches, based on integral equation, Euler initial value simulation, Euler matrix eigenvalue solution and Lagrangian particle simulation,…
In this paper we present a new ultra efficient numerical method for solving kinetic equations. In this preliminary work, we present the scheme in the case of the BGK relaxation operator. The scheme, being based on a splitting technique…
This paper presents a class of Crank-Nicolson (CN) type schemes enhanced by radial basis function (RBF) interpolation for the time integration of linear parabolic partial differential equations (PDEs). The resulting RBF-CN schemes preserve…
In this work, we aim at constructing numerical schemes, that are as efficient as possible in terms of cost and conservation of invariants, for the Vlasov--Fokker--Planck system coupled with Poisson or Amp\`ere equation. Splitting methods…
We present a collection of well-conditioned integral equation methods for the solution of electrostatic, acoustic or electromagnetic scattering problems involving anisotropic, inhomogeneous media. In the electromagnetic case, our approach…
In this paper we want to propose practical numerical methods to solve a class of initial-boundary problem of time-space fractional convection-diffusion equations (TSFCDEs). To start with, an implicit difference method based on two-sided…
The simulation of multi-body systems with frictional contacts is a fundamental tool for many fields, such as robotics, computer graphics, and mechanics. Hard frictional contacts are particularly troublesome to simulate because they make the…
A new method is proposed for integrating the equations of motion of an elastic filament. In the standard finite-difference and finite-element formulations the continuum equations of motion are discretized in space and time, but it is then…
We extend the asymptotic preserving and energy conserving time integrator for charged-particle motion developed in [Ricketson & Chac\'on, JCP, 2020] to include finite Larmor-radius (FLR) effects in the presence of electric-field…
In the last few decades, numerical simulation for nonlinear oscillators has received a great deal of attention, and many researchers have been concerned with the design and analysis of numerical methods for solving oscillatory problems. In…
We study the strong approximation of stochastic differential equations with discontinuous drift coefficients and (possibly) degenerate diffusion coefficients. To account for the discontinuity of the drift coefficient we construct an…
In this paper, we present a novel spectral renormalization exponential integrator method for solving gradient flow problems. Our method is specifically designed to simultaneously satisfy discrete analogues of the energy dissipation laws and…
This paper deals with a splitting method applied to a conservation law model of manufacturing system incorporating yield loss. A splitting scheme has been proposed. The yield loss term is treated by solving implicitly an ordinary…
We propose a family of integrators, Flow-Composed Implicit Runge-Kutta (FCIRK) methods, for perturbations of nonlinear ordinary differential equations, consisting of the composition of flows of the unperturbed part alternated with one step…