Related papers: Leapfrog methods for relativistic charged-particle…
In this paper, we rigorously analyze the energy, momentum and magnetic moment behaviours of two splitting methods for solving charged-particle dynamics. The near-conservations of these invariants are given for the system under constant…
A general Hamiltonian theory for the adiabatic motion of relativistic charged particles confined by slowly-varying background electromagnetic fields is presented based on a unified Lie-transform perturbation analysis in extended phase space…
A novel structure-preserving numerical method to solve random hyperbolic systems of conservation laws is presented. The method uses a concept of generalized, measure-valued solutions to random conservation laws. This yields a linear partial…
A new approach to the electroweak properties of two-particle composite systems is developed. The approach is based on the use of the instant form of relativistic Hamiltonian dynamics. The main novel feature of this approach is the new…
The Lorentz transformations are represented by Einstein velocity addition on the ball of relativistically admissible velocities. This representation is by projective maps. The Lie algebra of this representation defines the relativistic…
The goal of this paper is to develop energy-preserving variational integrators for time-dependent mechanical systems with forcing. We first present the Lagrange-d'Alembert principle in the extended Lagrangian mechanics framework and derive…
The relativistic fluid is a highly successful model used to describe the dynamics of many-particle, relativistic systems. It takes as input basic physics from microscopic scales and yields as output predictions of bulk, macroscopic motion.…
Numerical algorithms based on variational and symplectic integrators exhibit special features that make them promising candidates for application to general relativity and other constrained Hamiltonian systems. This paper lays part of the…
We propose a metriplectic reformulation of Lagrangian variational formulations for non-equilibrium thermodynamics. We prove that solutions to these constrained variational principles can be generated by the sum of a classic Poisson bracket…
Modern N-body techniques for planetary dynamics are generally based on symplectic algorithms specially adapted to the Kepler problem. These methods have proven very useful in studying planet formation, but typically require the timestep for…
We present a thermodynamically consistent theoretical framework for lyotropic liquid crystals (LCs) based on the GENERIC (General Equation for the Non-Equilibrium Reversible-Irreversible Coupling) formalism. This formalism ensures…
We develop a complete relativistic theory to describe the dynamics of electronic angular momentum including both spin (S) and orbital (L) contributions in magnetic systems. We start with the relativistic Dirac-Kohn-Sham Hamiltonian under…
This study proposes and analyses a novel higher-order, structure preserving discretization method for inviscid barotropic flows from a Lagrangian perspective. The method is built on a multisymplectic variational principle discretized over a…
We study the action and the dynamics of a relativistic particle, uncharged or charged, in multiscale spacetimes. Invariance under reparametrizations and Poincar\'e symmetries uniquely determine the action and the line element to be the…
In this article, a brief description of Discrete Mechanics and Variational Integrators which preserve the symplectic structure of the flow will be provided and a Newton-Raphson algorithm that can be used to solve implicit equations on the…
In this paper, for the first time we propose two linear, decoupled, energy-stable numerical schemes for multi-component two-phase compressible flow with a realistic equation of state (e.g. Peng-Robinson equation of state). The methods are…
We present a new method for solving the relativistic Vlasov--Maxwell system of equations, applicable to a wide range of extreme high-energy-density astrophysical and laboratory environments. The method directly discretizes the kinetic…
A very simple explicit integrator for the rotational motion of rigid linear molecules is presented which can preserve the rigidity of the molecules without requiring any constraint force. The integrator is time-reversible and symplectic,…
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…
The primary challenge in solving kinetic equations, such as the Vlasov equation, is the high-dimensional phase space. In this context, dynamical low-rank approximations have emerged as a promising way to reduce the high computational cost…