Related papers: Drift-preserving numerical integrators for stochas…
In this work, we aim at constructing numerical schemes, that are as efficient as possible in terms of cost and conservation of invariants, for the Vlasov--Fokker--Planck system coupled with Poisson or Amp\`ere equation. Splitting methods…
We propose an exponential integrator for the drift-kinetic equation in cylindrical geometry. This approach removes the CFL condition from the linear part of the system (which is often the most stringent requirement in practice) and treats…
We study the existence of log-canonical Poisson structures that are preserved by difference equations of special form. We also study the inverse problem, given a log-canonical Poisson structure to find a difference equation preserving this…
Stochastic contact Hamiltonian systems are a class of important mathematical models, which can describe the dissipative properties with odd dimensions in the stochastic environment. In this article, we investigate the numerical dynamics of…
Methods for discretizing port-Hamiltonian systems are of interest both for simulation and control purposes. Despite the large literature on mixed finite elements, no rigorous analysis of the connections between mixed elements and…
Geometric integrators of the Schr\"{o}dinger equation conserve exactly many invariants of the exact solution. Among these integrators, the split-operator algorithm is explicit and easy to implement, but, unfortunately, is restricted to…
In this paper, we develop two energy-preserving splitting methods for solving three-dimensional stochastic Maxwell equations driven by multiplicative noise. We use operator splitting methods to decouple stochastic Maxwell equations into…
A conforming finite element scheme with mixed explicit-implicit time discretization for quasi-incompressible Navier-Stokes-Maxwell-Stefan systems in a bounded domain with periodic boundary conditions is presented. The system consists of the…
There is a growing interest in the conservation of invariants when numerically solving a system of ordinary differential equations. Methods that exactly preserve these quantities in time are known as geometric integrators. In this paper we…
Many important physical systems can be described as the evolution of a Hamiltonian system, which has the important property of being conservative, that is, energy is conserved throughout the evolution. Physics Informed Neural Networks and…
We present a novel methodology for constructing arbitrarily high-order structure-preserving methods tailored for damped Hamiltonian systems. This method combines the idea of exponential integrator and energy-preserving collocation methods,…
The aim of this paper is to prove the Kolmogorov theorem of persistence of Diophantine flows for nearly-integrable Poisson systems associated to a real analytic Hamiltonian with aperiodic time dependence, provided that the perturbation is…
In this paper, we define arbitrarily high-order energy-conserving methods for Hamiltonian systems with quadratic holonomic constraints. The derivation of the methods is made within the so-called line integral framework. Numerical tests to…
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…
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…
A new splitting is proposed for solving the Vlasov-Maxwell system. This splitting is based on a decomposition of the Hamiltonian of the Vlasov-Maxwell system and allows for the construction of arbitrary high order methods by composition…
While symplectic integration methods based on operator splitting are well established in many branches of science, high order methods for Hamiltonian systems that split in more than two parts have not been studied in great detail. Here, we…
We propose a new explicit pseudo-energy and momentum conserving scheme for the time integration of Hamiltonian systems. The scheme, which is formally second-order accurate, is based on two key ideas: the integration during the time-steps of…
Variational integrators for Lagrangian dynamical systems provide a systematic way to derive geometric numerical methods. These methods preserve a discrete multisymplectic form as well as momenta associated to symmetries of the Lagrangian…
This paper aims to construct structure-preserving numerical schemes for multi-dimensional space fractional Klein-Gordon-Schr\"{o}dinger equation, which are based on the newly developed partitioned averaged vector field methods. First, we…