Related papers: Exact Charge-Conserving Scatter-Gather Algorithm f…
Until recently, electromagnetic finite element PIC (EM-FEMPIC) methods that demonstrated charge conservation used explicit field solvers. It is only recently, that a series of papers developed the mathematics necessary for charge…
Conventional explicit electromagnetic particle-in-cell (PIC) algorithms do not conserve discrete energy exactly. Time-centered fully implicit PIC algorithms can conserve discrete energy exactly, but may introduce large dispersion errors in…
The accuracy of Lagrangian point-particle models for simulation of particle-laden flows may degrade when the particle and fluid momentum equations are two-way coupled. In these cases the fluid velocity at the location of the particle, which…
In this paper we propose a new high order accurate space-time DG finite element scheme for the solution of the linear elastic wave equations in first order velocity-stress formulation in two and three-space dimensions on staggered…
In this paper, an implicit nonsymplectic exact energy-preserving integrator is specifically designed for a ten-dimensional phase-space conservative Hamiltonian system with five degrees of freedom. It is based on a suitable…
The approach presented in the second installment of this series is extended to multidimensional systems of conservation laws that are approximated via a Discontinuous Galerkin method on unstructured (triangular) grids. Special attention is…
This paper investigates the spectral structure, numerical dispersion, and observability of fully discrete approximations of the one-dimensional wave equation by $P^k$ (local) discontinuous Galerkin methods. Characterizing the coupled…
We address the generation of initial conditions (ICs) for GRAMSES, a code for nonlinear general relativistic (GR) $N$-body cosmological simulations recently introduced in Ref. [1]. GRAMSES adopts a constant mean curvature slicing with a…
Mathematical descriptions of flow phenomena usually come in the form of partial differential equations. The differential operators used in these equations may have properties such as symmetry, skew-symmetry, positive or negative…
Kinetic simulations of relativistic gases and plasmas are critical for understanding diverse astrophysical and terrestrial systems, but the accurate construction of the relativistic Maxwellian, the Maxwell-J\"uttner (MJ) distribution, on a…
This paper studies high-order accurate entropy stable nodal discontinuous Galerkin (DG) schemes for the ideal special relativistic magnetohydrodynamics (RMHD). It is built on the modified RMHD equations with a particular source term, which…
The Energy Conserving semi-implicit method (ECsim), presented by Lapenta in 2017, is a Particle in Cell (PIC) algorithm for the simulation of plasmas. Energy conservation is achieved within a semi-implicit formulation that does not require…
We extend the positivity-preserving method of Zhang & Shu (2010, JCP, 229, 3091-3120) to simulate the advection of neutral particles in phase space using curvilinear coordinates. The ability to utilize these coordinates is important for…
In this paper, we develop a high order structure-preserving local discontinuous Galerkin (DG) scheme for the compressible self-gravitating Euler equations, which pose great challenges due to the presence of time-dependent gravitational…
Depending on their sizes, dust grains store more or less charges, catalyse more or less chemical reactions, intercept more or less photons and stick more or less efficiently to form embryos of planets. Hence the need for an accurate…
This paper constitutes our initial effort in developing sparse grid discontinuous Galerkin (DG) methods for high-dimensional partial differential equations (PDEs). Over the past few decades, DG methods have gained popularity in many…
This work presents a novel interpolation-free mesh adaptation technique for the Euler equations within the arbitrary Lagrangian Eulerian framework. For the spatial discretization, we consider a residual distribution scheme, which provides a…
We propose an unsupervised machine-learning checkpoint-restart (CR) lossy algorithm for particle-in-cell (PIC) algorithms using Gaussian mixtures (GM). The algorithm features a particle compression stage and a particle reconstruction stage,…
We consider the discretization and subsequent model reduction of a system of partial differential-algebraic equations describing the propagation of pressure waves in a pipeline network. Important properties like conservation of mass,…
To cope with fast-fluctuating distributed energy resources (DERs) and uncontrolled loads, this paper formulates a time-varying optimization problem for distribution grids with DERs and develops a novel non-iterative algorithm to track the…