Related papers: A Charge Conserving Exponential Predictor Correcto…
A finite element discretization using a method of lines approached is proposed for approximately solving the Poisson-Nernst-Planck (PNP) equations. This discretization scheme enforces positivity of the computed solutions, corresponding to…
The so-called matrix-element method (MEM) has long been used successfully as a classification tool in particle physics searches. In the presence of invisible final state particles, the traditional MEM typically assigns probabilities to an…
A fundamental task in particle-in-cell (PIC) simulations of plasma physics is solving for charged particle motion in electromagnetic fields. This problem is especially challenging when the plasma is strongly magnetized due to numerical…
We present and analyze a linearized finite element method (FEM) for the dynamical incompressible magnetohydrodynamics (MHD) equations. The finite element approximation is based on mixed conforming elements, where Taylor--Hood type elements…
Peer methods are a comprehensive class of time integrators offering numerous degrees of freedom in their coefficient matrices that can be used to ensure advantageous properties, e.g. A-stability or super-convergence. In this paper, we show…
Recently a number of theoretical studies of the uniform electron gas (UEG) at finite temperature have appeared that are of relevance for dense plasmas, warm dense matter and laser excited solids and thermodynamic density functional theory…
In simulations of partial differential equations using particle-in-cell (PIC) methods, it is often advantageous to resample the particle distribution function to increase simulation accuracy, reduce compute cost, and/or avoid numerical…
Structure-preserving algorithms for solving conservative PDEs with added linear dissipation are generalized to systems with time-dependent damping/driving terms. This study is motivated by several PDE models of physical phenomena, such as…
We introduce the energy-stepping Monte Carlo (ESMC) method, a Markov chain Monte Carlo (MCMC) algorithm based on the conventional dynamical interpretation of the proposal stage but employing an energy-stepping integrator. The…
We develop a structure-preserving solution framework for the optimal control of the time-dependent Maxwell's equations. Building on a well-posedness theory for a weak form of the forward problem, we first analyze a forward solver that…
Local well-posedness for a nonlinear parabolic-hyperbolic coupled system modelling Micro-Electro-Mechanical System (MEMS) is studied. The particular device considered is a simple capacitor with two closely separated plates, one of which has…
In this paper, we construct a semi-implicit finite difference method for the time dependent Poisson-Nernst-Planck system. Although the Poisson-Nernst-Planck system is a nonlinear system, the numerical method presented in this paper only…
Energy-conserving Hermite methods for solving Maxwell's equations in dielectric and dispersive media are described and analyzed. In three space dimensions methods of order $2m$ to $2m+2$ require $(m+1)^3$ degrees-of-freedom per node for…
The rigorous stability analysis of high-order implicit-explicit multistep (IEMS) methods for nonlinear parabolic equations by using discrete energy arguments is a long standing open issue due to their non-A-stable property. A novel…
In this paper, a novel sixth order energy-conserved method is proposed for solving the three-dimensional time-domain Maxwell's equations. The new scheme preserves five discrete energy conservation laws, three momentum conservation laws,…
This article offers a new approach for analysing the dynamic behaviour of distributions of charged particles in an electromagnetic field. After discussing the limitations inherent in the Lorentz-Dirac equation for a single point particle a…
Many entropy-conservative and entropy-stable (summarized as entropy-preserving) methods for hyperbolic conservation laws rely on Tadmor's theory for two-point entropy-preserving numerical fluxes and its higher-order extension via flux…
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…
This paper is devoted to the numerical symplectic approximation of the charged-particle dynamics (CPD) with arbitrary electromagnetic fields. By utilizing continuous-stage methods and exponential integrators, a general class of symplectic…
Although persistent excitation is often acknowledged as a sufficient condition to exponentially converge in the field of adaptive parameter estimation, it must be noted that in practical applications this may be unguaranteed. Recently, more…