Related papers: A high-order Boris integrator
This article considers non-relativistic charged particle dynamics in both static and non-static electromagnetic fields, which are governed by nonseparable, possibly time-dependent Hamiltonians. For the first time, explicit symplectic…
We apply high-order mixed finite element discretization techniques and their associated preconditioned iterative solvers to the Variable Eddington Factor (VEF) equations in two spatial dimensions. The mixed finite element VEF…
A new family of methods involving complex coefficients for the numerical integration of differential equations is presented and analyzed. They are constructed as linear combinations of symmetric-conjugate compositions obtained from a basic…
We are studying the efficient solution of the system of linear equation stemming from the mass conserving mixed stress (MCS) method discretization of the Stokes equations. To that end we perform static condensation to arrive at a system for…
A wide range of implicit time integration methods, including multi-step, implicit Runge-Kutta, and Galerkin finite-time element schemes, is evaluated in the context of chaotic dynamical systems. The schemes are applied to solve the Lorenz…
We present sixth- and eighth-order Hermite integrators for astrophysical $N$-body simulations, which use the derivatives of accelerations up to second order ({\it snap}) and third order ({\it crackle}). These schemes do not require previous…
The computation of the radiative transfer equation is expensive mainly due to two stiff terms: the transport term and the collision operator. The stiffness in the former comes from the fact that particles (such as photons) travels at the…
This work deals with the numerical approximation of plasmas which are confined by the effect of a fast oscillating magnetic field (see \cite{Bostan2012}) in the Vlasov model. The presence of this magnetic field induces oscillations (in…
We present a novel methodology for constructing arbitrarily high-order structure-preserving methods tailored for damped Hamiltonian systems. This method combines the idea of exponential integrator and energy-preserving collocation methods,…
A new approach is developed to integrate numerically the equations of motion for systems of interacting rigid polyatomic molecules. With the aid of a leapfrog framework, we directly involve principal angular velocities into the integration,…
We present an adaptive arbitrary-order accurate time-stepping numerical scheme for the flow of vesicles suspended in Stokesian fluids. Our scheme can be summarized as an approximate implicit spectral deferred correction (SDC) method.…
The Morse-Ingard equations of thermoacoustics are a system of coupled time-harmonic equations for the temperature and pressure of an excited gas. They form a critical aspect of modeling trace gas sensors. In this paper, we analyze a…
We present SuperSCS: a fast and accurate method for solving large-scale convex conic problems. SuperSCS combines the SuperMann algorithmic framework with the Douglas-Rachford splitting which is applied on the homogeneous self-dual embedding…
Coupled systems of free flow and porous media arise in a variety of technical and environmental applications. For laminar flow regimes, such systems are described by the Stokes equations in the free-flow region and Darcy's law in the porous…
We develop a preconditioned nonlinear conjugate-gradient solver for ground states of binary dipolar Bose-Einstein condensates within the extended Gross-Pitaevskii equation including Lee-Huang-Yang corrections. The optimization is carried…
We change a previous time-stepping algorithm for solving a multi-scale Vlasov-Poisson system within a Particle-In-Cell method, in order to do accurate long time simulations. As an exponential integrator, the new scheme allows to use large…
By numerical integration, we study the relaxation dynamics of degenerate harmonic oscillator modes dispersively coupled to particle positions. Depending on whether the effective inertial potential induced by the oscillators keep the…
The main result in this paper is a provably entropy stable shock capturing approach for the high order entropy stable DGSEM based on a hybrid blending with a subcell low order variant. Since it is possible to rewrite a high order SBP…
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…
In atomic, molecular, and nuclear physics, the method of complex coordinate rotation is a widely used theoretical tool for studying resonant states. Here, we propose a novel implementation of this method based on the gradient optimization…