Related papers: A high-order Boris integrator
This paper proposes an implicit family of sub-step integration algorithms grounded in the explicit singly diagonally implicit Runge-Kutta (ESDIRK) method. The proposed methods achieve third-order consistency per sub-step and thus the…
When integrating the radiative transfer equation for polarized light, the necessity of high-order numerical methods is well known. In fact, well-performing high-order formal solvers enable higher accuracy and the use of coarser spatial…
A system of reduced equations is proposed for the electron motion in the strongly-radiation dominated regime for an arbitrary electromagnetic field configuration. The developed approach is used to analyze various scenarios of an electron…
We present a new molecular-dynamics algorithm for integrating the equations of motion for a system of particles interacting with mixed continuous/impulsive forces. This method, which we call Impulsive Verlet, is constructed using operator…
The boundary integral method is an efficient approach for solving time-harmonic obstacle scattering problems by a bounded scatterer. This paper presents the directional preconditioner for the iterative solution of linear systems of the…
The Convected Scheme (CS) is a `forward-trajectory' semi-Lagrangian method for solution of transport equations, which has been most often applied to the kinetic description of plasmas and rarefied neutral gases. In its simplest form, the CS…
We investigate vortex states in Bose-Einstein condensates under the combined action of the spin-orbit coupling (SOC), gradient magnetic field, and harmonic-oscillator trapping potential. The linear version of the system is solved exactly.…
In this paper, a well-conditioned collocation method is constructed for solving general $p$-th order linear differential equations with various types of boundary conditions. Based on a suitable Birkhoff interpolation, we obtain a new set of…
We provide an algorithm for evolving general spin-$s$ Gross-Pitaevskii / non-linear Schr\"odinger systems carrying a variety of interactions, where the $2s+1$ components of the `spinor' field represent the different spin-multiplicity…
We present a family of multistep integrators based on the Adams-Bashforth methods. These schemes can be constructed for arbitrary convergence order with arbitrary step size variation. The step size can differ between different subdomains of…
We consider the simulation of electromagnetic scattering by single and multiple isotropic homogeneous dielectric particles using boundary integral equations. Galerkin discretizations of the classical Poggio-Miller-Chang-Harrington-Wu-Tsai…
We propose an algorithm for low Mach number reacting flows subjected to electric field that includes the chemical production and transport of charged species. This work is an extension of a multi-implicit spectral deferred correction…
Conservative symmetric second-order one-step integrators are derived using the Discrete Multiplier Method for a family of vortex-blob models approximating the incompressible Euler's equations on the plane. Conservative properties and second…
Spectral deferred correction (SDC) methods are an attractive approach to iteratively computing collocation solutions to an ODE by performing so-called sweeps with a low-order time stepping method. SDC allows to easily construct high order…
The solvation model proposed by Fattebert and Gygi [Journal of Computational Chemistry 23, 662 (2002)] and Scherlis et al. [Journal of Chemical Physics 124, 074103 (2006)] is reformulated, overcoming some of the numerical limitations…
In this study, we develop an extended implicit moment method, namely, a coupled high-order low-order (HOLO) method and apply it to the electromagnetic Vlasov-Darwin model. The high-order (HO) system evolves particles in a manner that…
It is demonstrated that the inherent structure of Boris pusher may result in the accumulation of errors in numerical integration when estimating the works of electric fields on particles. This, in turn, leads to an incorrect estimations of…
The combined source integral equation (CSIE) for the electric field on the surface of a perfect electrically conducting scatterer can be discretized very accurately with lowest-order Rao-Wilton-Glisson basis and testing functions if the…
Electrokinetic phenomena in nanopore sensors and microfluidic devices require accurate simulation of coupled fluid-electrostatic interactions in geometrically complex domains with irregular boundaries and adaptive mesh refinement. We…
Network interactions between dynamical units are often subject to time delay. We develop a phase reduction method for delay-coupled oscillator networks. The method is based on rewriting the delay-differential equation as an ordinary…