Related papers: Higher Order Charge Conserving Electromagnetic Fin…
The high-frequency Helmholtz equation on the entire space is truncated into a bounded domain using the perfectly matched layer (PML) technique and subsequently, discretized by the higher-order finite element method (FEM) and the continuous…
MFEM is an open-source, lightweight, flexible and scalable C++ library for modular finite element methods that features arbitrary high-order finite element meshes and spaces, support for a wide variety of discretization approaches and…
We investigate the properties of the high-order discontinuous Galerkin spectral element method (DGSEM) with implicit backward-Euler time stepping for the approximation of hyperbolic linear scalar conservation equation in multiple space…
The EHP and the MCAP provide new rigorous weak variational formalism for a broad range of initial boundary value problems in mathematical physics and mechanics. Both approaches utilize the mixed formulation and lead to the development of…
The MFEM (Modular Finite Element Methods) library is a high-performance C++ library for finite element discretizations. MFEM supports numerous types of finite element methods and is the discretization engine powering many computational…
In recent years, the immersed finite element methods (IFEM) introduced in \cite{Li2003}, \cite{Li2004} to solve elliptic problems having an interface in the domain due to the discontinuity of coefficients are getting more attentions of…
We describe a charge-conserving scatter-gather algorithm for particle-in-cell simulations on unstructured grids. Charge conservation is obtained from first principles, i.e., without the need for any post-processing or correction steps. This…
The use of explicit particle-in-cell (PIC) method for relativistic plasma simulations is restricted by numerical heating and instabilities that may significantly constrain the choice of time and space steps. To partially eliminate these…
In this work, we present an efficient numerical implementation of the finite element method for modal analysis that leverages various symmetry operations, including spatial symmetry in point groups and space-time symmetry in…
A hyperbolic integro-differential equation is considered, as a model problem, where the convolution kernel is assumed to be either smooth or no worse than weakly singular. Well-posedness of the problem is studied in the context of semigroup…
Numerical homogenization for mechanical multiscale modeling by means of the finite element method (FEM) is an elegant way of obtaining structure-property relations, if the behavior of the constituents of the lower scale is well understood.…
The scalar and vector Laplacians are basic operators in physics and engineering. In applications, they show up frequently perturbed by lower-order terms. The effect of such perturbations on mixed finite element methods in the scalar case is…
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…
We consider the application of finite element exterior calculus (FEEC) methods to a class of canonical Hamiltonian PDE systems involving differential forms. Solutions to these systems satisfy a local multisymplectic conservation law, which…
In this paper, we develop a framework to construct energy-preserving methods for multi-components Hamiltonian systems, combining the exponential integrator and the partitioned averaged vector field method. This leads to numerical schemes…
Fast and accurate resolution of electromagnetic problems via the \ac{BEM} is oftentimes challenged by conditioning issues occurring in three distinct regimes: (i) when the frequency decreases and the discretization density remains constant,…
In this work, we present a novel family of explicit relaxation Particle-in-Cell (ER-PIC) methods for the Vlasov-Poisson equation with a strong magnetic field. These schemes achieve exact energy conservation by combining a splitting…
In recent years, several gauge-symmetric particle-in-cell (PIC) methods have been developed whose simulations of particles and electromagnetic fields exactly conserve charge. While it is rightly observed that these methods' gauge symmetry…
Constrained problems are frequently encountered in classical and quantum optimization. Particle conservation, in particular, is commonly imposed when studying energy spectra of chemical and solid state systems. Though particle…
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…