Related papers: GEMPIC: Geometric ElectroMagnetic Particle-In-Cell…
We present a new class of structure-preserving semi-discrete continuous-discontinuous Galerkin (CG-DG) finite element schemes for linear and nonlinear hyperbolic systems of partial differential equations on unstructured simplex meshes that…
We introduce an efficient split finite element (FE) discretization of a y-independent (slice) model of the rotating shallow water equations. The study of this slice model provides insight towards developing schemes for the full 2D case.…
The most popular methods for self-consistent simulation of fields interacting with charged species is using finite difference time domain (FDTD) methods together with Newton's laws of motion to evolve locations and velocities of particles.…
We present a novel structure-preserving framework for solving the Vlasov-Poisson-Landau system of equations using a particle in cell (PIC) discretization combined with discrete gradient time integrators. The Vlasov-Poisson-Landau system is…
A new splitting is proposed for solving the Vlasov-Maxwell system. This splitting is based on a decomposition of the Hamiltonian of the Vlasov-Maxwell system and allows for the construction of arbitrary high order methods by composition…
The state of the art in electromagnetic Finite Element Particle-in-Cell (EM-FEMPIC) has advanced significantly in the last few years. These have included understanding function spaces that must be used to represent sources and fields…
Semi-discrete and fully discrete mixed finite element methods are considered for Maxwell-model-based problems of wave propagation in linear viscoelastic solid. This mixed finite element framework allows the use of a large class of existing…
In this paper we propose a novel approach to discretize linear port-Hamiltonian systems while preserving the underlying structure. We present a finite element exterior calculus formulation that is able to mimetically represent conservation…
This paper introduces an auto-stabilized weak Galerkin (WG) finite element method with a built-in stabilizer for Poisson equations. By utilizing bubble functions as a key analytical tool, our method extends to both convex and non-convex…
We extend the recently-developed explicit, energy-conserving particle-in-cell (PIC) scheme of [1] to the relativistic Vlasov-Maxwell system. As in the non-relativistic case, the method is built on an optimization problem that is…
Hamiltonian and action principle (HAP) formulations of plasma physics are reviewed for the purpose of explaining structure preserving numerical algorithms. Geometric structures associated with and emergent from HAP formulations are…
In the presence of an inhomogeneous oscillatory electric field, charged particles experience a net force, averaged over the oscillatory timescale, known as the ponderomotive force. We derive a one-dimensional Hamiltonian model which…
Explicit high-order non-canonical symplectic particle-in-cell algorithms for classical particle-field systems governed by the Vlasov-Maxwell equations are developed. The algorithm conserves a discrete non-canonical symplectic structure…
This paper considers weak Galerkin finite element approximations for a quasistatic Maxwell viscoelastic model. The spatial discretization uses piecewise polynomials of degree $k \ (k\geq 1)$ for the stress approximation, degree $k+1$ for…
The particle-in-cell (PIC) method has been widely used for plasma simulation, because of its noise-reduction capability and moderate computational cost. The immersed finite element (IFE) method is efficient for solving interface problems on…
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…
JAX-in-Cell is a fully electromagnetic, multispecies, and relativistic 1D3V Particle-in-Cell (PIC) framework implemented entirely in JAX. It provides a modern, Python-based alternative to traditional PIC frameworks. It leverages…
We study a system of Maxwell's equations that describes the time evolution of electromagnetic fields with an additional electric scalar variable to make the system amenable to a mixed finite element spatial discretization. We demonstrate…
In this work we propose a high-order structure-preserving discretization of the cold plasma model which describes the propagation of electromagnetic waves in magnetized plasmas. By utilizing B-Splines Finite Elements Exterior Calculus, we…
We present a geometric Particle-in-Cell (PIC) algorithm on two-dimensional (2D) unstructured meshes for studying electrostatic perturbations in magnetized plasmas. In this method, ions are treated as fully kinetic particles, and electrons…