Related papers: Explicit high-order noncanonical symplectic algori…
Explicit high-order non-canonical symplectic particle-in-cell algorithms for classical particle-field systems governed by the Vlasov-Maxwell equations are developed. The algorithm conserves a discrete non-canonical symplectic structure…
Symplectic schemes are powerful methods for numerically integrating Hamiltonian systems, and their long-term accuracy and fidelity have been proved both theoretically and numerically. However direct applications of standard symplectic…
Given a fluid equation with reduced Lagrangian $l$ which is a functional of velocity $\MM{u}$ and advected density $D$ given in Eulerian coordinates, we give a general method for semidiscretising the equations to give a canonical…
We study the non-canonical symplectic structure, or K-symplectic structure inherited by the charged particle dynamics. Based on the splitting technique, we construct non-canonical symplectic methods which is explicit and stable for the…
We present an adaptation of the so-called structural method \cite{CMM23} for Hamiltonian systems, and redesign the method for this specific context, which involves two coupled differential systems. Structural schemes decompose the problem…
Studying single-particle dynamics over many periods of oscillations is a well-understood problem solved using symplectic integration. Such integration schemes derive their update sequence from an approximate Hamiltonian, guaranteeing that…
We derive the dynamics of several rigid bodies of arbitrary shape in a 2-dimensional inviscid and incompressible fluid, whose vorticity field is given by point vortices. We adopt the idea of Vankerschaver et al. (2009) to derive the…
Symplectic integrators for Hamiltonian systems have been quite successful for studying few-body dynamical systems. These integrators are frequently derived using a formalism built on symplectic maps. There have been recent efforts to extend…
Dynamics of a charged particle in the canonical coordinates is a Hamiltonian system, and the well-known symplectic algorithm has been regarded as the de facto method for numerical integration of Hamiltonian systems due to its long-term…
We study the geometric particle-in-cell methods for an electrostatic hybrid plasma model. In this model, ions are described by the fully kinetic equations, electron density is determined by the Boltzmann relation, and space-charge effects…
Hamilton's equations are fundamental for modeling complex physical systems, where preserving key properties such as energy and momentum is crucial for reliable long-term simulations. Geometric integrators are widely used for this purpose,…
Symplectic integration methods based on operator splitting are well established in many branches of science. For Hamiltonian systems which split in more than two parts, symplectic methods of higher order have been studied in detail only for…
Relativistic dynamics of a charged particle in time-dependent electromagnetic fields has theoretical significance and a wide range of applications. It is often multi-scale and requires accurate long-term numerical simulations using…
We present a brief tutorial on the nuts and bolts computation of a multisymplectic particle-in-cell algorithm using the discretized Lagrangian approach. This approach, originated by Marsden, Shadwick, and others, brings the benefits of…
A numerically stable method to solve the discretized Boltzmann-Enskog equation describing the behavior of non ideal fluids under inhomogeneous conditions is presented. The algorithm employed uses a Lagrangian finite-difference scheme for…
Hamiltonian systems are known to conserve the Hamiltonian function, which describes the energy evolution over time. Obtaining a numerical spatio-temporal scheme that accurately preserves the discretized Hamiltonian function is often a…
In this article, we propose high-order finite-difference entropy stable schemes for the two-fluid relativistic plasma flow equations. This is achieved by exploiting the structure of the equations, which consists of three independent flux…
We introduce novel high order well-balanced finite volume methods for the full compressible Euler system with gravity source term. They require no a priori knowledge of the hydrostatic solution which is to be well-balanced and are not…
We present the symplectic algorithm in the Lagrangian formalism for the Hamiltonian systems by virtue of the noncommutative differential calculus with respect to the discrete time and the Euler--Lagrange cohomological concepts. We also show…
The barotropic ideal fluid with step and delta-function discontinuities coupled to Einstein's gravity is studied. The discontinuities represent star surfaces and thin shells; only non-intersecting discontinuity hypersurfaces are considered.…