Related papers: How to regularize a symplectic-energy-momentum int…
Many applications in computational physics that use numerical integrators based on splitting and composition can benefit from the development of optimized algorithms and from choosing the best ordering of terms. The cost in programming and…
We construct a symplectic integrator for non-separable Hamiltonian systems combining an extended phase space approach of Pihajoki and the symmetric projection method. The resulting method is semiexplicit in the sense that the main time…
Hamiltonian systems are differential equations which describe systems in classical mechanics, plasma physics, and sampling problems. They exhibit many structural properties, such as a lack of attractors and the presence of conservation…
Explicit and semi-explicit geometric integration schemes for dissipative perturbations of Hamiltonian systems are analyzed. The dissipation is characterized by a small parameter $\epsilon$, and the schemes under study preserve the…
Measurement data is often sampled irregularly i.e. not on equidistant time grids. This is also true for Hamiltonian systems. However, existing machine learning methods, which learn symplectic integrators, such as SympNets [20] and…
We present a simple algorithm to switch between $N$-body time integrators in a reversible way. We apply it to planetary systems undergoing arbitrarily close encounters and highly eccentric orbits, but the potential applications are broader.…
A new approach is developed to integrate numerically the equations of motion for systems of interacting rigid polyatomic molecules. With the aid of a leapfrog framework, we directly involve principal angular velocities into the integration,…
In this work, a kernel-based surrogate for integrating Hamiltonian dynamics that is symplectic by construction and tailored to large prediction horizons is proposed. The method learns a scalar potential whose gradient enters a…
Two documents have recently been posted on the arXiv describing a numerical integration algorithm: "symplectic orbit/spin tracking code for all-electric storage rings" [1] and some computational results therefrom [2]. This note comments…
Symplectic integrators are the preferred method of solving conservative $N$-body problems in cosmological, stellar cluster, and planetary system simulations because of their superior error properties and ability to compute orbital…
This paper presents transient numerical simulations of hydraulic systems in engineering applications using the spectral element method (SEM). Along with a detailed description of the underlying numerical method, it is shown that the SEM…
Method of moment estimators exhibit appealing statistical properties, such as asymptotic unbiasedness, for nonconvex problems. However, they typically require a large number of samples and are extremely sensitive to model misspecification.…
We develop a new numerical technique for approximating solutions of the Navier-Stokes equations on moving domains. The method aims at simulating an incompressible fluid past an object whose motion is assigned a priori using a level-set…
We study dynamics of an wheeled inverted pendulum under a proportional-integral-derivative controller on horizontal, inclined and soft surfaces. An oscillatory area and conditions of the stability for the control are shown on the phase…
Airborne transient electromagnetic (TEM) is a cost-effective method to image the distribution of electrical conductivity in the ground. We consider layered earth inversion to interpret large data sets of hundreds of kilometre. Different…
In order to solve continuous-time optimal control problems, direct methods transcribe the infinite-dimensional problem to a nonlinear program (NLP) using numerical integration methods. In cases where the integration error can be manipulated…
We continue the presentation of the pragmatic mode-sum regularization (PMR) method for computing the renormalized stress-energy tensor (RSET). We show in detail how to employ the $t$-splitting variant of the method, which was first…
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…
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…
Symplectic N-body integrators are widely used to study problems in celestial mechanics. The most popular algorithms are of 2nd and 4th order, requiring 2 and 6 substeps per timestep, respectively. The number of substeps increases rapidly…