Related papers: A New Time-Dependent Finite Difference Method for …
A full-F, isothermal, electromagnetic, gyro-fluid model is used to simulate plasma turbulence in a COMPASS-sized, diverted tokamak. A parameter scan covering three orders of magnitude of plasma resistivity and two values for the ion to…
This paper discusses temporally continuous and discrete forms of the speed-limited particle-in-cell (SLPIC) method first treated by Werner et al. [Phys. Plasmas 25, 123512 (2018)]. The dispersion relation for a 1D1V electrostatic plasma…
We explore a novel simulation route for Plasma Wakefield Acceleration (PWFA) by using the computational method known as the Lattice Boltzmann Method (LBM). LBM is based on a discretization of the continuum kinetic theory while assuring the…
In this paper, we present an efficient algorithm for the long time behavior of plasma simulations. We will focus on 4D drift-kinetic model, where the plasma's motion occurs in the plane perpendicular to the magnetic field and can be…
In this paper we apply a finite difference lattice Boltzmann model to study the phase separation in a two-dimensional liquid-vapor system. Spurious numerical effects in macroscopic equations are discussed and an appropriate numerical scheme…
We have developed a one-dimensional code to solve ultra-relativistic hydrodynamic problems, using the Glimm method for an accurate treatment of shocks and contact discontinuities. The implementation of the Glimm method is based on an exact…
This paper presents novel methods to approximate the nonlinear AC optimal power flow (OPF) into tractable linear/quadratic programming (LP/QP) based OPF problems that can be used for power system planning and operation. We derive a linear…
Diffusive transport processes in magnetized plasmas are highly anisotropic, with fast parallel transport along the magnetic field lines sometimes faster than perpendicular transport by orders of magnitude. This constitutes a major challenge…
Numerical resolution of high-dimensional nonlinear PDEs remains a huge challenge due to the curse of dimensionality. Starting from the weak formulation of the Lawson-Euler scheme, this paper proposes a stochastic particle method (SPM) by…
Monte-Carlo computations for highly relativistic parallel shock particle acceleration are presented for upstream flow gamma factors, $\Gamma=(1-V_{1}^{2}/c^{2})^{-0.5}$ with values between 5 and $10^{3}$. The results show that the spectral…
This study addresses the challenge of simulating realistic particle systems by proposing a novel particle decomposition scheme that improves the parallel performance of surface resolved particle simulations. Realistic particle systems often…
Shocks waves are a ubiquitous feature of many astrophysical plasma systems, and an important process for energy dissipation and transfer. The physics of these shock waves are frequently treated/modeled as a collisional, fluid MHD…
Suspensions with fiber-like particles in the low Reynolds number regime are modeled by two different approaches that both use a Lagrangian representation of individual particles. The first method is the well-established formulation based on…
We adapt and modify the eigenfunction method of computing the power-law spectrum of particles accelerated at a relativistic shock front via the first-order Fermi process (Kirk, J.G., Schneider, P., Astrophysical Journal 315, 425 (1987)) to…
This paper introduces a new model for highly accurate distribution voltage solutions, coined as a parameterized linear power flow model. The proffered model is grounded on a physical model of linear power flow equations, and uses…
In this paper we present an investigation of numerical Monte Carlo simulations of the diffusive shock acceleration in the test particle limit. Very high gamma flow astrophysical plasmas, have been used, from $\gamma_{up}$ $\sim50$ up to…
Relativistic collisionless shocks are associated with efficient particle acceleration when propagating into weakly magnetized homogeneous media; as the magnetization increases, particle acceleration becomes suppressed. We demonstrate that…
We study the rate of convergence of an explicit and an implicit-explicit finite difference scheme for linear stochastic integro-differential equations of parabolic type arising in non-linear filtering of jump-diffusion processes. We show…
In this paper, we derive a new drift-kinetic method for estimating gradients in the plasma properties through a velocity space distribution at a single point. The gradients are intrinsically related to agyrotropic features of the…
Cosmic rays are charged particles that are accelerated to relativistic speeds by astrophysical shocks. Numerical models have been successful in confirming the acceleration process for (quasi-)parallel shocks, which have the magnetic field…