Related papers: A high-order Boris integrator
High-resolution simulations of particle-based kinetic plasma models typically require a high number of particles and thus often become computationally intractable. This is exacerbated in multi-query simulations, where the problem depends on…
In this study, we consider a numerical implementation of the nonlinear Rosenbluth-Trubnikov collision operator for particle simulations in plasma physics in the framework of the finite element method (FEM). The relevant particle evolution…
The transport of charged particles or photons in a scattering medium can be modelled with a Boltzmann equation. The mathematical treatment for scattering in such scenarios is often simplified if evaluated in a frame where the scattering…
In this study, we explore multiple higher-order pole solutions in spinor Bose--Einstein condensates. These solutions are associated with different pairs of higher-order poles of the transmission coefficient in the inverse scattering…
The efficient computation of the overdamped, random motion of micron and nanometre scale particles in a viscous fluid requires novel methods to obtain the hydrodynamic interactions, random displacements and Brownian drift at minimal cost.…
We develop a robust and efficient iterative method for hyper-elastodynamics based on a novel continuum formulation recently developed. The numerical scheme is constructed based on the variational multiscale formulation and the…
We revisit the Hierarchical Poincar\'e-Steklov (HPS) method in a preconditioned iterative setting for variable-coefficient Helmholtz problems with impedance boundary conditions. HPS is commonly presented as a direct solver based on nested…
We propose a computationally lean, two-stage approach that reliably predicts self-assembly behavior of complex charged molecules on a metallic surfaces under electrochemical conditions. Stage one uses ab initio simulations to provide…
In this work we propose and analyze a new Hybrid High-Order method for the Brinkman problem for fluids with power-law viscosity. The proposed method supports general meshes and arbitrary approximation orders and is robust in all regimes,…
Several different approaches are proposed for solving fully implicit discretizations of a simplified Boltzmann-Poisson system with a linear relaxation-type collision kernel. This system models the evolution of free electrons in…
This article deals with the numerical integration in time of nonlinear Schr\"odinger equations. The main application is the numerical simulation of rotating Bose-Einstein condensates. The authors perform a change of unknown so that the…
We propose a very sensitive method for measuring the electric dipole moment of the electron using polar molecules embedded in a cryogenic solid matrix of inert-gas atoms. The polar molecules can be oriented in the $\hat{\rm{z}}$ direction…
We show that applying any deterministic B-series method of order $p_d$ with a random step size to single integrand SDEs gives a numerical method converging in the mean-square and weak sense with order $\lfloor p_d/2\rfloor$.As an…
In this paper we generalize and improve a recently developed domain decomposition preconditioner for the iterative solution of discretized Helmholtz equations. We introduce an improved method for transmission at the internal boundaries…
In this paper, we apply the optimized Schwarz method to the two dimensional nonlinear Schr{\"o}dinger equation and extend this method to the simulation of Bose-Einstein condensates (Gross-Pitaevskii equation). We propose an extended version…
We construct a high-order adaptive time stepping scheme for vesicle suspensions with viscosity contrast. The high-order accuracy is achieved using a spectral deferred correction (SDC) method, and adaptivity is achieved by estimating the…
We study dynamic self-organisation and order-disorder transitions in a two-dimensional system of self-propelled particles. Our model is a variation of the Vicsek model, where particles align the motion to their neighbours but repel each…
Next-generation high-power lasers that can be focused to intensities exceeding 10^23 W/cm^2 are enabling new physics and applications. The physics of how these lasers interact with matter is highly nonlinear, relativistic, and can involve…
The paper introduces the sweeping preconditioner, which is highly efficient for iterative solutions of the variable coefficient Helmholtz equation including very high frequency problems. The first central idea of this novel approach is to…
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…