Related papers: Functionally-fitted energy-preserving integrators …
We discuss structure-preserving numerical discretizations for repulsive and attractive Euler-Poisson equations that find applications in fluid-plasma and self-gravitation modeling. The scheme is fully discrete and structure preserving in…
Variational integrators are a special kind of geometric discretisation methods applicable to any system of differential equations that obeys a Lagrangian formulation. In this thesis, variational integrators are developed for several…
The Volterra lattice equations are completely integrable and possess bi-Hamiltonian structure. They are integrated using partitioned Lobatto IIIA-B methods which preserve the Poisson structure. Modified equations are derived for the…
Suppression of rectification at metal--Mott-insulator interfaces, which is previously shown by numerical solutions to the time-dependent Schr\"odinger equation and experiments on real devices, is reinvestigated theoretically by…
We present a new class of exponential integrators for ordinary differential equations. They are locally exact, i.e., they preserve the linearization of the original system at every point. Their construction consists in modifying existing…
Hamiltonian integration methods for the Vlasov-Maxwell equations are developed by a Hamiltonian splitting technique. The Hamiltonian functional is split into five parts, i.e., the electrical energy, the magnetic energy, and the kinetic…
In this paper structure-preserving time-integrators for rigid body-type mechanical systems are derived from a discrete Hamilton-Pontryagin variational principle. From this principle one can derive a novel class of variational partitioned…
We give a wide class of Lie-Poisson systems for which explicit, Lie-Poisson integrators, preserving all Casimirs, can be constructed. The integrators are extremely simple. Examples are the rigid body, a moment truncation, and a new, fast…
The dynamical motion of mechanical systems possesses underlying geometric structures, and preserving these structures in numerical integration improves the qualitative accuracy and reduces the long-time error of the simulation. For a single…
We analyze a new framework for expressing finite element methods on arbitrarily many intersecting meshes: multimesh finite element methods. The multimesh finite element method, first presented in [40], enables the use of separate meshes to…
We investigate discretization strategies for a recently introduced class of energy-based models. The model class encompasses classical port-Hamiltonian systems, generalized gradient flows, and certain systems with algebraic constraints. Our…
We build a new mathematical model of shape optimization for maximizing ionic concentration governed by the multi-physical coupling steady-state Poisson-Nernst-Planck system. Shape sensitivity analysis is performed to obtain the Eulerian…
A mechanism for the reduction of dynamic energy dissipation in the computing circuit is described. The resonant circuit with controlled switches conserves the stored energy by recovering upto 90% of energy that would be otherwise lost…
We present a new linearly implicit exponential integrator that preserves the polynomial first integrals or Lyapunov functions for the conservative and dissipative stiff equations, respectively. The method is tested by both oscillated…
We study finite-difference approximations of both Poisson and Poisson-Boltzmann (PB) electrostatic energy functionals for periodic structures constrained by Gauss' law and a class of local algorithms for minimizing the finite-difference…
We introduce a family of fourth order two-step methods that preserve the energy function of canonical polynomial Hamiltonian systems. Each method in the family may be viewed as a correction of a linear two-step method, where the correction…
We show how the integrators used for the molecular dynamics step of the Hybrid Monte Carlo algorithm can be further improved. These integrators not only approximately conserve some Hamiltonian $H$ but conserve exactly a nearby shadow…
First-order energy dissipative schemes in time are available in literature for the Poisson-Nernst-Planck (PNP) equations, but second-order ones are still in lack. This work proposes novel second-order discretization in time and finite…
We present two novel classes of fully discrete energy-preserving algorithms for the sine-Gordon equation subject to Neumann boundary conditions. The cosine pseudo-spectral method is first used to develop structure-preserving spatial…
A novel class of high-order linearly implicit energy-preserving integrating factor Runge-Kutta methods are proposed for the nonlinear Schr\"odinger equation. Based on the idea of the scalar auxiliary variable approach, the original equation…