Related papers: An exponential integrator for the drift-kinetic mo…
We propose, analyze, and implement a quadrature method for evaluating integrals of the form $\int_0^2 f(s)\exp(zs)\, {\rm d}s$, where $z$ is a complex number with a possibly large negative real part. The integrand may exhibit exponential…
The first part of this paper introduces sufficient conditions to determine conservation laws of diffusion equations of arbitrary fractional order in time. Numerical methods that satisfy a discrete analogue of these conditions have…
This paper deals with the construction and analysis of two integrators for (semi-linear) second-order partial differential-algebraic equations of semi-explicit type. More precisely, we consider an implicit-explicit Crank-Nicolson scheme as…
In this article, we present a new second order finite difference discrete scheme for fractal mobile/immobile transport model based on equivalent transformative Caputo formulation. The new transformative formulation takes the singular kernel…
We present a novel hierarchical formulation of the fourth-order forward symplectic integrator and its numerical implementation in the GPU-accelerated direct-summation N-body code FROST. The new integrator is especially suitable for…
This paper presents a hybrid numerical method for linear collisional kinetic equations with diffusive scaling. The aim of the method is to reduce the computational cost of kinetic equations by taking advantage of the lower dimensionality of…
We investigate a high-order, fully explicit, asymptotic-preserving scheme for a kinetic equation with linear relaxation, both in the hydrodynamic and diffusive scalings in which a hyperbolic, resp. parabolic, limiting equation exists. The…
Taking advantage of the flexibility of the variational method with coordinate transformations, we derive a self-consistent set of equations of motion from a discretized Lagrangian to study kinetic plasmas using a Fourier decomposed…
We suggest a numerical integration procedure for solving the equations of motion of certain classical spin systems which preserves the underlying symplectic structure of the phase space. Such symplectic integrators have been successfully…
This paper proposes a new non-oscillatory {\em energy-splitting} conservative algorithm for computing multi-fluid flows in the Eulerian framework. In comparison with existing multi-fluid algorithms in literatures, it is shown that the mass…
We show how to increase the order of one-dimensional discrete gradient numerical integrator without losing its advantages, such as exceptional stability, exact conservation of the energy integral and exact preservation of the trajectories…
This letter studies symmetric and symplectic exponential integrators when applied to numerically computing nonlinear Hamiltonian systems. We first establish the symmetry and symplecticity conditions of exponential integrators and then show…
This paper presents an extension of the hybrid scheme proposed by Wang et al. (J. Comput. Phys. 229 (2010) 169-180) for numerical simulation of compressible isotropic turbulence to flows with higher turbulent Mach numbers. The scheme still…
An optimally efficient explicit numerical scheme for solving fluid dynamics equations, or any other parabolic or hyperbolic system of partial differential equations, should allow local regions to advance in time with their own, locally…
In this paper, we consider the task of efficiently computing the numerical solution of evolutionary complex Ginzburg--Landau equations on Cartesian product domains with homogeneous Dirichlet/Neumann or periodic boundary conditions. To this…
Allowing for space- and time-dependence of mass in Klein--Gordon equations resolves the problem of negative probability density and of violation of Lorenz covariance of interaction in quantum mechanics. Moreover it extends their…
We point out that two of Milne's fourth-order integrators are well-suited to bit-reversible simulations. The fourth-order method improves on the accuracy of Levesque and Verlet's algorithm and simplifies the definition of the velocity $v$…
Dynamic systems have a fundamental relevance in the description of physical phenomena. The search for more accurate and faster numerical integration methods for the resolution of such systems is, therefore, an important topic of research.…
In this paper, we derive a variational integrator for certain highly oscillatory problems in mechanics. To do this, we take a new approach to the splitting of fast and slow potential forces: rather than splitting these forces at the level…
A new exponentially fitted version of the Discrete Variational Derivative method for the efficient solution of oscillatory complex Hamiltonian Partial Differential Equations is proposed. When applied to the nonlinear Schroedinger equation,…