Related papers: A Charge Conserving Exponential Predictor Correcto…
In this paper, we introduce and discuss an exactly energy-conserving Particle-in-Cell method for arbitrary curvilinear coordinates. The flexibility provided by curvilinear coordinates enables the study of plasmas in complex-shaped domains…
We propose a time-adaptive predictor/multi-corrector method to solve hyperbolic partial differential equations, based on the generalized-$\alpha$ scheme that provides user-control on the numerical dissipation and second-order accuracy in…
The Poisson-Nernst-Planck (PNP) equations are one of the most effective model for describing electrostatic interactions and diffusion processes in ion solution systems, and have been widely used in the numerical simulations of biological…
The conformal formulation of the Einstein constraint equations is first reviewed, and we then consider the design, analysis, and implementation of adaptive multilevel finite element-type numerical methods for the resulting coupled nonlinear…
This paper presents a novel 3-D full electromagnetic particle-in-cell (PIC) code called JefiPIC, which uses Jefimenko's equations as the electromagnetic (EM) field solver through a full-space integration method. Leveraging the power of…
The trapezoidal rule, which is a special case of the Newmark family of algorithms, is one of the most widely used methods for transient hyperbolic problems. In this work, we show that this rule conserves linear and angular momenta and…
Ultracold-atom simulations of the Hubbard model provide insights into the character of charge and spin correlations in and out of equilibrium. The corresponding numerical simulations, on the other hand, remain a significant challenge. We…
Energy conserving particle-in-cell schemes are constructed for a class of reduced relativistic Vlasov--Maxwell equations of laser-plasma interaction. Discrete Poisson equation is also satisfied by the numerical solution. Specifically,…
This paper inspects more closely the problem of the momentum and energy of a bound (non-radiating) electromagnetic (EM) field. It has been shown that for an isolated system of non-relativistic mechanically free charged particles a…
For an arbitrary parameter $p\in [1,+\infty]$, we consider the problem of exponential stabilization in the spatial $L^{p}$-norm, and $W^{1,p}$-norm, respectively, for a class of anti-stable linear parabolic PDEs with space-time-varying…
This paper studies economic model predictive Control (EMPC) schemes, where the stage cost depends only on control inputs. Such problems arise in applications like water distribution networks and differ from standard EMPC since multiple…
In this study, we present a finite element solver for a thermodynamically consistent electrolyte model that accurately captures multicomponent ionic transport by incorporating key physical phenomena such as steric effects, solvation, and…
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 this work, a composite economic model predictive control (CEMPC) is proposed for the optimal operation of a stand-alone integrated energy system (IES). Time-scale multiplicity exists in IESs dynamics is taken into account and addressed…
In a previous paper, we developed a new particle-in-cell method for the Vlasov-Maxwell system in which the electromagnetic fields and the equations of motion for the particles were cast in terms of scalar and vector potentials through a…
This paper concerns the inclusion of Newton's method into an adaptive finite element method (FEM) for the solution of nonlinear partial differential equations (PDEs). It features an adaptive choice of the damping parameter in the Newton…
We consider the preconditioned conjugate gradient method (PCG) with optimal preconditioner in the frame of the boundary element method (BEM) for elliptic first-kind integral equations. Our adaptive algorithm steers the termination of PCG as…
We prove sharp wavenumber-explicit error bounds for first- or second-family-N\'ed\'elec-element (a.k.a. edge-element) conforming discretisations, of arbitrary (fixed) order, of the variable-coefficient time-harmonic Maxwell equations posed…
We use partial differential equations (PDEs) to describe physical systems. In general, these equations include evolution and constraint equations. One method used to find solutions to these equations is the Free-evolution approach, which…
We explore quantum corrections of electrically charged black holes subject to vacuum polarization effects of fermion fields in QED. Solving this problem exactly is challenging so we restrict to perturbative corrections that one can obtain…