Related papers: Time-symmetry, symplecticity and stability of Eule…
As is well known, energy is generally deemed as one of the most important physical invariants in many conservative problems and hence it is of remarkable interest to consider numerical methods which are able to preserve it. In this paper,…
In this paper we propose and investigate a general approach to constructing local energy-preserving algorithms which can be of arbitrarily high order in time for solving Hamiltonian PDEs. This approach is based on the temporal…
This letter investigates the Lie point symmetries and conserved quantities of the Lagrangian systems on time scales, which unify the Lie symmetries of the two cases for the continuous and the discrete Lagrangian systems. By defining the…
Simulation of many-particle system evolution by molecular dynamics takes to decrease integration step to provide numerical scheme stability on the sufficiently large time interval. It leads to a significant increase of the volume of…
A new family of methods involving complex coefficients for the numerical integration of differential equations is presented and analyzed. They are constructed as linear combinations of symmetric-conjugate compositions obtained from a basic…
Noether and Lie symmetry analyses based on point transformations that depend on time and spatial coordinates will be reviewed for a general class of time-dependent Hamiltonian systems. The resulting symmetries are expressed in the form of…
The usual explicit finite-difference method of solving partial differential equations is limited in stability because it approximates the exact amplification factor by power-series. By adapting the same exponential-splitting method of…
The purpose of this paper is to describe geometrically discrete Lagrangian and Hamiltonian Mechanics on Lie groupoids. From a variational principle we derive the discrete Euler-Lagrange equations and we introduce a symplectic 2-section,…
We propose a linearly implicit structure-preserving numerical method for semilinear Hamiltonian systems with polynomial nonlinearities, combining Kahan's method and exponential integrator. This approach efficiently balances computational…
Time-reversible symplectic methods, which are precisely compatible with Liouville's phase-volume-conservation theorem, are often recommended for computational simulations of Hamiltonian mechanics. Lack of energy drift is an apparent…
On this paper, we have proposed an approach to observe the time-centered difference scheme for dissipative mechanical systems from a Hamiltonian perspective and to introduce the idea of symplectic algorithm to dissipative systems. The…
In recent decades, there have been many attempts to construct symplectic integrators with variable time steps, with rather disappointing results. In this paper we identify the causes for this lack of performance, and find that they fall…
Including the effect of thermal fluctuations in traditional computational fluid dynamics requires developing numerical techniques for solving the stochastic partial differential equations of fluctuating hydrodynamics. These Langevin…
A linear evolving surface partial differential equation is first discretized in space by an arbitrary Lagrangian Eulerian (ALE) evolving surface finite element method, and then in time either by a Runge-Kutta method, or by a backward…
An Euler-type framework with equidistant step sizes is proposed for a class of time-changed stochastic differential equations.We establish the strong convergence rate of the standard Euler--Maruyama method under the global Lipschitz…
By combining a standard symmetric, symplectic integrator with a new step size controller, we provide an integration scheme that is symmetric, reversible and conserves the values of the constants of motion. This new scheme is appropriate for…
Continuous Galerkin Petrov time discretization scheme is tested on some Hamiltonian systems including simple harmonic oscillator, Kepler's problem with different eccentricities and molecular dynamics problem. In particular, we implement the…
It has been argued that it is incompatible to maintain unitary time-evolution for time-dependent non-Hermitian Hamiltonians when the metric operator is explicitly time-dependent. We demonstrate here that the time-dependent Dyson equation…
We consider symplectic time integrators in numerical General Relativity and discuss both free and constrained evolution schemes. For free evolution of ADM-like equations we propose the use of the Stoermer-Verlet method, a standard…
Dynamic simulation of elastic bodies is a longstanding task in engineering and computer graphics. In graphics, numerical integrators like implicit Euler and BDF2 are preferred due to their stability at large time steps, but they tend to…