Related papers: Nystr\"om Type Exponential Integrators for Strongl…
In this paper, we implement exponential integrators, specifically Integrating Factor (IF) and Exponential Time Differencing (ETD) methods, using pseudo-spectral techniques to solve phase-field equations within a Python framework. These…
Until recently, electromagnetic finite element PIC (EM-FEMPIC) methods that demonstrated charge conservation used explicit field solvers. It is only recently, that a series of papers developed the mathematics necessary for charge…
A model of an electron-beam-plasma system is introduced to model the electrical breakdown physics of low-pressure nitrogen irradiated by an intense pulsed electron beam. The rapidly rising beam current induces an electric field which drives…
The accurate prediction of occurrence and strength of kinetic instabilities in plasmas remains a significant challenge in nuclear fusion research. To accurately capture the plasmas dynamics one is required to solve the Vlasov equation for…
In recent work, we initiated a research program aimed at the systematic investigation of quantum superintegrable systems describing the interaction of two non-relativistic spin-$1/2$ particles in three-dimensional Euclidean space. In that…
In this paper, exponential energy-preserving methods are formulated and analysed for solving charged-particle dynamics in a strong and constant magnetic field. The resulting method can exactly preserve the energy of the dynamics. Moreover,…
Global particle-in-cell (PIC) simulations of pulsar magnetospheres are performed with a volume, surface and pair production-based plasma injection schemes to systematically investigate the transition between electrosphere and force-free…
In view of the existing limitations of sequential computing, parallelization has emerged as an alternative in order to improve the speedup of numerical simulations. In the framework of evolutionary problems, space-time parallel methods…
We present a MATLAB toolbox for five different classes of exponential integrators for solving (mildly) stiff ordinary differential equations or time-dependent partial differential equations. For the efficiency of such exponential…
Many problems in astrophysics involve relativistic outflows. The plasma dynamics in these scenarios is critical to determine the conditions for the self-consistent evolution of the fields and particle acceleration. Advances in computer…
We report on the observation of rapid particle acceleration in numerical simulations of relativistic jet-plasma interactions and discuss the underlying mechanisms. The dynamics of a charge-neutral, narrow, electron-positron jet propagating…
We have developed efficient numerical algorithms for solving 3D steady-state Poisson-Nernst-Planck (PNP) equations with excess chemical potentials described by the classical density functional theory (cDFT). The coupled PNP equations are…
We introduce an extension of the particle-in-cell (PIC) method that captures the Landau collisional effects in the Vlasov-Maxwell-Landau equations. The method arises from a regularisation of the variational formulation of the Landau…
In this paper we present a model based on dynamics of the electrons in the plasma using a simplified Boltzmann equation coupled with a Poisson equation. The motivation arose to simulate active plasma resonance spectroscopy which is used for…
Particle acceleration in collisionless plasma systems is a central question in astroplasma and astroparticle physics. The structure of the acceleration regions, electron-ion energy equilibration, preacceleration of particles at shocks to…
We present results from a fully kinetic particle-in-cell (PIC) simulation of direct current magnetron sputtering (dcMS) in a 2D cylindrically symmetric geometry. The particle-in-cell model assumes an electrostatic approximation and includes…
Kinetic Particle In Cell (PIC) methods can extend greatly their range of applicability if implicit time differencing and spatial adaption are used to address the wide range of time and length scales typical of plasmas. For implicit…
In this paper, we propose a $\mu$-mode integrator for computing the solution of stiff evolution equations. The integrator is based on a $d$-dimensional splitting approach and uses exact (usually precomputed) one-dimensional matrix…
Current physics models used to interpret experimental measurements of particle beams require either simplifying assumptions to be made in order to ensure analytical tractability, or black box optimization methods to perform model based…
Symplectic integrators are widely used for long-term integration of conservative astrophysical problems due to their ability to preserve the constants of motion; however, they cannot in general be applied in the presence of nonconservative…