Related papers: A Charge Conserving Exponential Predictor Correcto…
We present a class of non-standard numerical schemes which are modifications of the discrete gradient method. They preserve the energy integral exactly (up to the round-off error). The considered class contains locally exact discrete…
We examine the stability of an Einstein-Maxwell perfect fluid configuration with a privileged direction of symmetry by means of a $1+1+2$-tetrad formalism. We use this formalism to cast, in a quasi linear symmetric hyperbolic form the…
In this work, we develop and rigorously analyze a new class of particle methods for the magnetized Vlasov--Poisson--Fokker--Planck system. The proposed approach addresses two fundamental challenges: (1) the curse of dimensionality, which we…
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…
Achieving large-scale kinetic modelling is a crucial task for the development and optimization of modern plasma devices. With the trend of decreasing pressure in applications such as plasma etching, kinetic simulations are necessary to…
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…
Simulations of the dynamics generated by partial differential equations (PDEs) provide approximate, numerical solutions to initial value problems. Such simulations are ubiquitous in scientific computing, but the correctness of the results…
Present day electromagnetic field calculations have limitations that are due to techniques employing edge-based discretization methods. While these vector finite element methods solve the issues of tangential continuity of fields and the…
This paper introduces a numerical scheme for time harmonic Maxwell's equations by using weak Galerkin (WG) finite element methods. The WG finite element method is based on two operators: discrete weak curl and discrete weak gradient, with…
A nonlinear Helmholtz equation (NLH) with high wave number and Sommerfeld radiation condition is approximated by the perfectly matched layer (PML) technique and then discretized by the linear finite element method (FEM).…
The leapfrog integrator is widely used because of its excellent stability in molecular dynamics simulation. This is recognized as being due to the existence of a discrete variational structure of the equations. We introduce a modified…
We present an efficient and accurate energy-conserving implicit particle-in-cell~(PIC) algorithm for the electrostatic Vlasov system, with particular emphasis on its high robustness for simulating complex plasma systems with multiple…
An asymptotic preserving and energy stable scheme for the Euler-Poisson system under the quasineutral scaling is designed and analysed. Correction terms are introduced in the convective fluxes and the electrostatic potential, which lead to…
Two-fluid plasma flow equations describe the flow of ions and electrons with different densities, velocities, and pressures. We consider the ideal plasma flow i.e. we ignore viscous, resistive, and collision effects. The resulting system of…
As is well known, energy is generally deemed as one of the most important physical invariants in many conservative problems and hence it is of remarkable interest to consider numerical methods which are able to preserve it. In this paper,…
We construct a finite element method (FEM) for the infinity Laplacian. Solutions of this problem may be singular, which has prompted us to conduct an a posteriori analysis of the method deriving residual based estimators to drive an…
The propagation of electromagnetic waves in general media is modeled by the time-dependent Maxwell's partial differential equations (PDEs), coupled with constitutive laws that describe the response of the media. In this work, we focus on…
The rapid deposition of energy by Edge Localised Modes (ELMs) onto plasma facing components, is a potentially serious issue for large Tokamaks such as ITER and DEMO. The trigger for ELMs is believed to be the ideal Magnetohydrodynamic…
In this paper, we propose an efficient exponential integrator finite element method for solving a class of semilinear parabolic equations in rectangular domains. The proposed method first performs the spatial discretization of the model…
Energy methods for constructing time-stepping algorithms are of increased interest in application to nonlinear problems, since numerical stability can be inferred from the conservation of the system energy. Alternatively, symplectic…