Related papers: Symplectic Integrator Mercury: Bug Report
We present a structure preserving discretization of the fundamental spacetime geometric structures of fluid mechanics in the Lagrangian description in 2D and 3D. Based on this, multisymplectic variational integrators are developed for…
We derive, and discuss the properties of, a symplectic map for the dynamics of bodies on nearly parabolic orbits. The orbits are perturbed by a planet on a circular, coplanar orbit interior to the pericenter of the parabolic orbit. The map…
An extended Hamiltonian approach to conduct isothermal-isobaric molecular dynamics simulations with full cell flexibility is presented. The components of the metric tensor are used as the fictitious degrees of freedom for the cell, thus…
Integrability in string/field theories is known to emerge when considering dynamics in the moduli space of physical theories. This implies that one has to look at the dynamics with respect to unusual time variables like coupling constants…
We consider the simplest classical integrable model corresponding to a non-hyperelliptic spectral curve. We show that a certain complicated integral occurs when computing the average of observables in this model. This integral does not…
We present a hybrid mimetic finite-difference and virtual element formulation for coupled single-phase poromechanics on unstructured meshes. The key advantage of the scheme is that it is convergent on complex meshes containing highly…
An efficient geometric integrator is proposed for solving the perturbed Kepler motion. This method is stable and accurate over long integration time, which makes it appropriate for treating problems in astrophysics, like solar system…
The simulation of rare events is one of the key problems in atomistic simulations. Towards its solution a plethora of methods have been proposed. Here we combine two such methods metadynamics and inte-grated tempering sampling. In…
We present a numerically exact Inchworm Monte Carlo method for equilibrium multiorbital quantum impurity problems with general interactions and hybridizations. We show that the method, originally developed to overcome the dynamical sign…
For studies of the long-term evolution of small Solar System objects, it is fundamental to add the Yarkovsky and Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effects in the dynamical model. Still, implementations of these effects in…
There are two fundamental problems studied by the theory of hamiltonian integrable systems: integration of equations of motion, and construction of action-angle variables. The third problem, however, should be added to the list: separation…
State-of-the-art robotics simulators operate in discrete time. This requires users to choose a time step, which is both critical and challenging: large steps can produce non-physical artifacts, while small steps force the simulation to run…
Calculating the long term solution of ordinary differential equations, such as those of the $N$-body problem, is central to understanding a wide range of dynamics in astrophysics, from galaxy formation to planetary chaos. Because generally…
The paper deals with an integrodifferential operator which models numerous phenomena in superconductivity, in biology and in viscoelasticity. Initialboundary value problems with Neumann, Dirichlet and mixed boundary conditions are analyzed.…
We consider a three-dimensional Fourier integral in which the exponent in the exponential factor is the product of some phase function and a large parameter. The asymptotics of this integral is sought when the large parameter tends to…
The long-term dynamics of many dynamical systems evolve on an attracting, invariant "slow manifold" that can be parameterized by a few observable variables. Yet a simulation using the full model of the problem requires initial values for…
We develop the equations of motion for full body models that describe the dynamics of rigid bodies, acting under their mutual gravity. The equations are derived using a variational approach where variations are defined on the Lie group of…
Symplectic integrators evolve dynamical systems according to modified Hamiltonians whose error terms are also well-defined Hamiltonians. The error of the algorithm is the sum of each error Hamiltonian's perturbation on the exact solution.…
We study the optimal design of numerical integrators for dissipative systems, for which there exists an underlying thermodynamic structure known as GENERIC (general equation for the nonequilibrium reversible-irreversible coupling). We…
Variational integrators are well-suited for simulation of mechanical systems because they preserve mechanical quantities about a system such as momentum, or its change if external forcing is involved, and holonomic constraints. While they…