Related papers: Controlling the Numerical Cerenkov Instability in …
A new fluid model for runaway electron simulation based on fluid description is introduced and implemented in the magnetohydrodynamics code M3D-C1, which includes self-consistent interactions between plasma and runaway electrons. The model…
Plasmonic nanostructures significantly improve the performance of photoconductive devices (PCDs) in generating terahertz radiation. However, they are geometrically intricate and result in complicated electromagnetic (EM) field and carrier…
This work aims to construct an efficient and highly accurate numerical method to address the time singularity at $t=0$ involved in a class of time-fractional parabolic integro-partial differential equations in one and two dimensions. The…
Modern particle-in-cell (PIC) codes have become an integral tool in plasma astrophysics. As most plasma phenomena grow from initially small instabilities, it is important to ensure PIC codes can suppress noise and ensure that any growing…
In wave propagation problems, finite difference methods implemented on staggered grids are commonly used to avoid checkerboard patterns and to improve accuracy in the approximation of short-wavelength components of the solutions. In this…
The time evolution of a collisionless plasma is modeled by the Vlasov-Maxwell system which couples the Vlasov equation (the transport equation) with the Maxwell equations of electrodynamics. We only consider a 'two-dimensional' version of…
The numerical solution of high dimensional Vlasov equation is usually performed by particle-in-cell (PIC) methods. However, due to the well-known numerical noise, it is challenging to use PIC methods to get a precise description of the…
In this paper, we present three-dimensional (3D) Particle-In-Cell (PIC) simulations to study the stability of 2D Bernstein-Greene-Kruskal (BGK) modes in a magnetized plasma with a finite background magnetic field. The simulations were…
Differentiable programming allows for derivatives of functions implemented via computer code to be calculated automatically. These derivatives are calculated using automatic differentiation (AD). This thesis explores two applications of…
A finite-difference time-domain (FDTD) modelling of finite-size zero thickness space-time modulated Huygens' metasurfaces based on Generalized Sheet Transition Conditions (GSTCs), is proposed and numerically demonstrated. A typical…
In this work, several key scaling laws of the quasi-static magnetic compression of field reversed configuration (FRC) plasma [Spencer, Tuszewski, and Linford, 1983] are amended from a series of 2D FRC MHD equilibriums numerically obtained…
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…
The properties of linear instabilities in the Large Plasma Device [W. Gekelman et al., Rev. Sci. Inst., 62, 2875 (1991)] are studied both through analytic calculations and solving numerically a system of linearized collisional plasma fluid…
The first order by time partial differential equations are used as models in applications such as fluid flow, heat transfer, solid deformation, electromagnetic waves, and others. In this paper we propose the new numerical method to solve a…
Particle-in-Cell (PIC) simulation is the most important numerical tool in plasma physics. However, its long-term accuracy has not been established. To overcome this difficulty, we developed a canonical symplectic PIC method for the…
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…
We present a non-staggered method for the Maxwell equations in adaptively refined grids. The code is based on finite volume central scheme that preserves in a discrete form both divergence-free property of magnetic field and the Gauss law.…
This paper develops and analyzes an interior penalty discontinuous Galerkin (IPDG) method using piecewise linear polynomials for the indefinite time harmonic Maxwell equations with the impedance boundary condition in the three dimensional…
Coherent multimode instabilities are responsible for several phenomena of recent interest in semiconductor lasers, such as the generation of frequency combs and ultrashort pulses. These techonologies have proven disruptive in optical…
A class of nonstandard pseudospectral time domain (PSTD) schemes for solving time-dependent hyperbolic and parabolic partial differential equations (PDEs) is introduced. These schemes use the Fourier collocation spectral method to compute…