Related papers: Controlling the Numerical Cerenkov Instability in …
Two-fluid plasma flow equations describe the flow of ions and electrons with different densities, velocities, and pressures. We consider the ideal plasma flow i.e. we ignore viscous, resistive, and collision effects. The resulting system of…
In this work, we propose staggered FDTD schemes based on the correction function method (CFM) to discretize Maxwell's equations with embedded perfect electric conductor (PEC) boundary conditions. The CFM uses a minimization procedure to…
The present study addresses the challenge of enhancing computational efficiency without compromising accuracy in numerical simulations of vacuum gas dynamics using the direct simulation Monte Carlo (DSMC) method. A technique termed "fixed…
The advent of high-power Hall thrusters and the increasing interest towards their use as a primary propulsion system for various missions have given a new boost to the efforts aiming at self-consistent predictive modeling of this thruster…
A novel 3-D higher-order finite-difference time-domain framework with complex frequency-shifted perfectly matched layer for the modeling of wave propagation in cold plasma is presented. Second- and fourth-order spatial approximations are…
Context: Calculating stellar pulsations requires a sufficient accuracy to match the quality of the observations. Many current pulsation codes apply a second order finite-difference scheme, combined with Richardson extrapolation to reach…
Noncollinear (NC) magnetism and spin-orbit coupling (SOC) are indispensable for predictive ab initio materials simulations with pronounced relativistic effects and magnetic frustration, yet they significantly increase the cost of…
The rotation profile of a magnetized plasma cylinder composed of multiple fluids is investigated analytically, expanding on previous results. The analytic steady-state solution is used as a benchmark for a time-dependent multiple-fluid…
Relativistic magnetic reconnection is a non-ideal plasma process that is a source of non-thermal particle acceleration in many high-energy astrophysical systems. Particle-in-cell (PIC) methods are commonly used for simulating reconnection…
We propose a high-order adaptive numerical solver for the semilinear elliptic boundary value problem modelling magnetic plasma equilibrium in axisymmetric confinement devices. In the fixed boundary case, the equation is posed on curved…
The simulation of fracture using continuum ductile damage models attains a pathological discretization dependence caused by strain localization, after loss of ellipticity of the problem, in regions whose size is connected to the spatial…
Fast, high-order accurate algorithms for electromagnetic scattering from axisymmetric objects are of great importance when modeling physical phenomena in optics, materials science (e.g. meta-materials), and many other fields of applied…
We present results of one-dimensional collisionless simulations of plasma turbulence and related anomalous electron current of the Electron Cyclotron Drift Instability (ECDI). Our highly resolved, long-term simulations of xenon plasma in…
In this work, we present a numerical method that remedies the instabilities of the conventional FDTD approach for solving Maxwell's equations in a space-time dependent magneto-electric medium with direct application to the simulation of the…
We present a multi-block finite-difference solver for massively parallel Direct Numerical Simulations (DNS) of incompressible flows. The algorithm combines the versatility of a multi-block solver with the method of eigenfunctions…
An axial-azimuthal two-dimensional Hall thruster discharge model was developed for analyzing gradient drift instability (GDI) and cross-field electron transport enhancement induced solely by the GDI. A hybrid particle-fluid model was used…
In collisionless and weakly collisional plasmas, the particle distribution function is a rich tapestry of the underlying physics. However, actually leveraging the particle distribution function to understand the dynamics of a weakly…
Entity is a new-generation, fully open-source particle-in-cell (PIC) code developed to overcome key limitations in astrophysical plasma modeling, particularly the extreme separation of scales and the performance challenges associated with…
A fundamental task in particle-in-cell (PIC) simulations of plasma physics is solving for charged particle motion in electromagnetic fields. This problem is especially challenging when the plasma is strongly magnetized due to numerical…
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…