Related papers: An \mathcal{O}(N) Maxwell solver with improved num…
We present a numerical model and a set of conservative algorithms for Non-Maxwellian plasma kinetics with inelastic collisions. These algorithms self-consistently solve for the time evolution of an isotropic electron energy distribution…
A numerical method for solving elliptic PDEs with variable coefficients on two-dimensional domains is presented. The method is based on high-order composite spectral approximations and is designed for problems with smooth solutions. The…
The self-consistent modeling of vacuum polarization due to virtual electron-positron fluctuations is of relevance for many near term experiments associated with high intensity radiation sources and represents a milestone in describing…
In this paper we present a newly developed GRPIC code framework called FPIC, providing a detailed description of the Maxwell-equations solver, of the particle ``pushers'', and of the other algorithms that are needed in this approach. We…
A method to compute guaranteed lower bounds to the eigenvalues of the Maxwell system in two or three space dimensions is proposed as a generalization of the method of Liu and Oishi [SIAM J. Numer. Anal., 51, 2013] for the Laplace operator.…
A new method for numerical solving of boundary problem for ordinary differential equations with slowly varying coefficients which is aimed at better representation of solutions in the regions of their rapid oscillations or exponential…
We introduce a new notion of "matrix potential" to nonlinear optical systems. In terms of a matrix potential $g$, we present a gauge field theoretic formulation of the Maxwell-Bloch equation that provides a semiclassical description of the…
Scientific codes are an indispensable link between theory and experiment; in (astro-)plasma physics, such numerical tools are one window into the universe's most extreme flows of energy. The discretization of Maxwell's equations - needed to…
In this paper, using the Newman-Penrose formalism, we find the Maxwell equations in NUT space and after separation into angular and radial components solve them analytically. All the angular equations are solved in terms of Jaccobi…
Smoothing short-wavelength charge density variations can stabilize explicit electrostatic particle-in-cell (PIC) plasma simulations against grid heating and cold beam instabilities, which cause unphysical heating when the Debye length is…
This article presents an $O(N\log N)$ method for numerical solution of Maxwell's equations for dielectric scatterers using a 3D boundary integral equation (BIE) method. The underlying BIE method used is based on a hybrid…
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…
Numerical modeling of electromagnetic waves is an important tool for understanding the interaction of light and matter, and lies at the core of computational electromagnetics. Traditional approaches to injecting and evolving electromagnetic…
The Particle-in-Cell (PIC) simulation has been a widely used method for studying plasma physics. However, fully three-dimensional PIC simulations always require huge computational resources. For problems with near azimuthal symmetry, recent…
We argue that the Maxwellian approximation can essentially underestimate the rates of some nuclear reactions in hot plasma under conditions very close to thermal equilibrium. This phenomenon is demonstrated explicitly on the example of…
We introduce a Galilean electromagnetic particle-in-cell (GEM-PIC) algorithm, which transforms the full set of Maxwell equations and the Vlasov equation into the boosted coordinates. This approach preserves the electromagnetic structure of…
In this paper, we propose a weak Galerkin (WG) finite element method for the Maxwell eigenvalue problem. By restricting subspaces, we transform the mixed form of Maxwell eigenvalue problem into simple elliptic equation. Then we give the WG…
The Vlasov--Maxwell equations are used for the kinetic description of magnetized plasmas. As they are posed in an up to 3+3 dimensional phase space, solving this problem is extremely expensive from a computational point of view. In this…
The mathematical description of laboratory fusion plasmas produced in Tokamaks is still challenging. Complete models for electrons and ions, as Vlasov-Maxwell systems, are computationally too expensive because they take into account all…
We introduce a quasi-static particle-in-cell (PIC) code -- WAND-PIC -- which does not suffer from some of the common limitations of many quasi-static PICs, such as the need for a predictor-corrector method in solving electromagnetic fields.…