Related papers: Structure-Preserving Integration for Magnetic Gaus…
Numerical methods that preserve geometric invariants of the system, such as energy, momentum or the symplectic form, are called geometric integrators. In this paper we present a method to construct symplectic-momentum integrators for…
In this paper, we propose a mass- and modified energy-conservative relaxation Crank-Nicolson finite element method for the Schr\"{o}dinger-Poisson equation. Utilizing only a single auxiliary variable, we simultaneously reformulate the…
We construct explicit integrators of arbitrary even orders of accuracy for massive point vortex dynamics in binary mixture of Bose--Einstein condensates proposed by Richaud et al. The integrators are symplectic and preserve the angular…
In this paper, we develop a structure-preserving discretization of the Lagrangian framework for electromagnetism, combining techniques from variational integrators and discrete differential forms. This leads to a general family of…
In this paper we present energy-conserving, mixed discontinuous Galerkin (DG) and continuous Galerkin (CG) schemes for the solution of a broad class of physical systems described by Hamiltonian evolution equations. These systems often arise…
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…
In this paper, we explore the chaotic signatures of the geodesic dynamics for particles moving in the slowly rotating Hartle-Thorne spacetime; an approximate solution of vacuum Einstein field equations describing the exterior of a massive,…
We revisit the Scovel-Weinstein framework (Scovel & Weinstein, CPAM 1994) for reducing the Vlasov-Poisson system while preserving its Hamiltonian structure. Standard particle-in-cell (PIC) algorithms approximate the distribution function by…
Variational time integrators are derived in the context of discrete mechanical systems. In this area, the governing equations for the motion of the mechanical system are built following two steps: (a) Postulating a discrete action; (b)…
We develop a geometric framework for the numerical integration of mechanical systems evolving on manifolds. After briefly reviewing classical numerical methods and highlighting their limitations and shortcomings in non-flat (non-Euclidean)…
We perform a numerical analysis of a class of randomly perturbed {H}amiltonian systems and {P}oisson systems. For the considered additive noise perturbation of such systems, we show the long time behavior of the energy and quadratic…
Simulation and analysis of multidimensional dynamics of a quantum non-Hmeritian system is a challenging problem. Gaussian wavepacket dynamics has proven to be an intuitive semiclassical approach to approximately solving the dynamics of…
We derive an exact analytical solution to the time-dependent Schr\"odinger equation for transmission of a Gaussian wave packet through an arbitrary potential of finite range. We consider the situation where the initial Gaussian wave packet…
We present a novel structure-preserving framework for solving the Vlasov-Poisson-Landau system of equations using a particle in cell (PIC) discretization combined with discrete gradient time integrators. The Vlasov-Poisson-Landau system is…
We present a construction of the Anisotropic Gaussian Semi-Classical Schr\"{o}dinger Propagator, emblematic of a class of Fourier Integral Operators of quadratic phase kernels related to the Schr\"{o}dinger equation. We deduce a set of…
Stokes-Dirac structures are infinite-dimensional Dirac structures defined in terms of differential forms on a smooth manifold with boundary. These Dirac structures lay down a geometric framework for the formulation of Hamiltonian systems…
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…
Several aspects of the connection between conserved integrals (invariants) and symmetries are illustrated within a hybrid Lagrangian-Hamiltonian framework for dynamical systems. Three examples are considered: a nonlinear oscillator with…
An effective operational approach to quantum mechanics is to focus on the evolution of wave-packets, for which the wave-function can be seen in the semi-classical regime as representing a classical motion dressed with extra degrees of…
We develop a structure-preserving numerical discretization for the electrostatic Euler-Poisson equations with a constant magnetic field. The scheme preserves positivity of the density, positivity of the internal energy and a minimum…