Related papers: Nystr\"om Type Exponential Integrators for Strongl…
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…
We describe an initial implementation of an electrostatic Particle-in-Cell (ES-PIC) module with adaptive Cartesian mesh in our Unified Flow Solver framework. Challenges of PIC method with cell-based adaptive mesh refinement (AMR) are…
Abridged. Kinetic plasma simulations are nowadays commonly used to study a wealth of non-linear behaviours and properties in laboratory and space plasmas. In particular, in high-energy physics and astrophysics, the plasma usually evolves in…
In this paper, we formulate and analyse exponential integrations when applied to nonlinear Schr\"{o}dinger equations in a normal or highly oscillatory regime. A kind of exponential integrators with energy preservation, optimal convergence…
We propose a class of Particle-In-Cell (PIC) methods for the Vlasov-Poisson system with a strong and inhomogeneous external magnetic field with fixed direction, where we focus on the motion of particles in the plane orthogonal to the…
The numerical integration of stiff equations is a challenging problem that needs to be approached by specialized numerical methods. Exponential integrators form a popular class of such methods since they are provably robust to stiffness and…
We present a new multi-fluid, multi-temperature plasma solver with adaptive Cartesian mesh (ACM) based on a full-Newton (non-linear, implicit) scheme for collisional low-temperature plasma. The particle transport is described using the…
This paper discusses a novel fully implicit formulation for a 1D electrostatic particle-in-cell (PIC) plasma simulation approach. Unlike earlier implicit electrostatic PIC approaches (which are based on a linearized Vlasov-Poisson…
Numerical simulations of relativistic plasmas have become more feasible, popular, and crucial for various astrophysical sources with the availability of computational resources. The necessity for high-accuracy particle dynamics is…
Two-fluid (electron-positron) plasma modelling has shown that inductive acceleration can convert Poynting flux directly into bulk kinetic energy in the relativistic flows driven by rotating magnetized neutron stars and black holes. Here, we…
Particle in Cell (PIC) simulations have become a vital tool for the investigation of kinetic processes in plasma physics. Many of the systems investigated with PIC simulations contain particles with relativistic velocities. The correct…
Molecular Dynamics method is based on solution of Newtonian differential equations of motion. A new very accurate and efficient time-reversible explicit integrator was derived on the basis of second order Tailor expansion of force. There is…
It is possible to excite various linear and non-linear low-frequency modes in dusty plasma which is an admixture of electrons, ions, gas atoms, and negatively charged solid particles. The experimental as well as theoretical study of these…
We present and analyze fully discrete Nystr\"om methods for the solution of three classes of well conditioned boundary integral equations for the solution of two dimensional scattering problems by homogeneous dielectric scatterers.…
Implicit particle-in-cell codes offer advantages over their explicit counterparts in that they suffer weaker stability constraints on the need to resolve the higher frequency modes of the system. This feature may prove particularly valuable…
Plasma supports collective modes and particle-wave interactions that leads to complex behavior in inertial fusion energy applications. While plasma can sometimes be modeled as a charged fluid, a kinetic description is useful towards the…
Modeling of kinetic plasmas using electromagnetic particle in cell methods (EM-PIC) is a problem that is well worn, in that methods developed have been used extensively both understanding physics and exploiting them for device design.…
In many problems involving particle accelerators and relativistic plasmas, the accurate modeling of relativistic particle motion is essential for accurate physical predictions. Here, we extend a charge-conserving finite element time-domain…
This paper deals with the numerical integration of Hamiltonian systems in which a stiff anharmonic potential causes highly oscillatory solution behavior with solution-dependent frequencies. The impulse method, which uses micro- and…
We present a novel numerical method for solving the elliptic partial differential equation problem for the electrostatic potential with piecewise constant conductivity. We employ an integral equation approach for which we derive a system of…