Related papers: Exact local conservation of energy in fully implic…
In this paper, we extend the operator-split asymptotic-preserving, semi-Lagrangian algorithm for time dependent anisotropic heat transport equation proposed in [Chac\'on et al., JCP, 272, 719-746, 2014] to use a fully implicit time…
This paper introduces a novel formulation of the Particle-In-Cell (PIC) method for the modeling of relativistic plasmas, which leverages the ability of the Pseudo-Spectral Analytical Time-Domain solver (PSATD) to handle arbitrary time…
Orbital-free density functional theory (OF-DFT) runs at low computational cost that scales linearly with the number of simulated atoms, making it suitable for large-scale material simulations. It is generally considered that OF-DFT strictly…
We establish sharp energy decay rates for a large class of nonlinearly first-order damped systems, and we design discretization schemes that inherit of the same energy decay rates, uniformly with respect to the space and/or time…
In this paper, we present a linearly implicit energy-preserving scheme for the Camassa-Holm equation by using the multiple scalar auxiliary variables approach, which is first developed to construct efficient and robust energy stable schemes…
Energy loss of energetic ions in solid is crucial in many field, and accurate prediction of the ion stopping power is a long-time goal. Though great efforts have been made, it is still very difficult to find a universal prediction model to…
We introduce a new definition of discrete-time port-Hamiltonian systems (PHS), which results from structure-preserving discretization of explicit PHS in time. We discretize the underlying continuous-time Dirac structure with the collocation…
The robustness of dynamical systems against external perturbations is crucial in engineering; however, it is often overlooked for the lack of methods for rapidly computing it. This paper proposes a novel algorithm for estimating the…
We design a novel, exactly energy-conserving implicit non-symplectic integration method for an eight-dimensional Hamiltonian system with four degrees of freedom. In our algorithm, each partial derivative of the Hamiltonian with respect to…
A recent study by Gassner et al. [J. Sci. Comput. 90:79 (2022)] demonstrates that local energy stability--that is, ensuring the asymptotic numerical growth rate does not exceed the continuous growth rate--is crucial for achieving accurate…
It is shown analytically that the energy-conserving implicit nonsymplectic scheme of Bacchini, Ripperda, Chen and Sironi provides a first-order accuracy to numerical solutions of a six-dimensional conservative Hamiltonian system. Because of…
In this paper, we design, analyze, and numerically validate positive and energy-dissipating schemes for solving the time-dependent multi-dimensional system of Poisson-Nernst-Planck (PNP) equations, which has found much use in the modeling…
The integration of intermittent and volatile renewable energy resources requires increased flexibility in the operation of the electric grid. Storage, broadly speaking, provides the flexibility of shifting energy over time; network, on the…
In this paper we study arbitrarily high-order energy-conserving methods for simulating the dynamics of a charged particle. They are derived and studied within the framework of Line Integral Methods (LIMs), previously used for defining…
We present a fully iterative adaptive algorithm for the numerical minimization of strongly convex energy functionals in Hilbert spaces. The proposed approach, which we first present in abstract form, generates a hierarchical sequence of…
Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These…
Particle-in-Cell (PIC) simulation is the most important numerical tool in plasma physics. However, its long-term accuracy has not been established. To overcome this difficulty, we developed a canonical symplectic PIC method for the…
Energy extraction from quantum batteries by means of completely positive trace-preserving (CPTP) maps leads to the concept of CPTP-local passive states, which identify bipartite states from which no energy can be squeezed out by applying…
We demonstrate that the concept of a conservation law can be naturally extended from deterministic to probabilistic cellular automata (PCA) rules. The local function for conservative PCA must satisfy conditions analogous to conservation…
Particle-In-Cell (PIC) methods are frequently used for kinetic, high-fidelity simulations of plasmas. Implicit formulations of PIC algorithms feature strong conservation properties, up to numerical round-off errors, and are not subject to…