Related papers: How to regularize a symplectic-energy-momentum int…
We describe a method for symplectic tracking of charged particles through static electric and magnetic fields. The method can be applied to cases where the fields have a dependence on longitudinal as well as transverse position, and where…
Seismic traveltime tomography using transmission data is widely used to image the Earth's interior from global to local scales. In seismic imaging, it is used to obtain velocity models for subsequent depth-migration or full-waveform…
The shearing sheet is a model dynamical system that is used to study the small-scale dynamics of astrophysical disks. Numerical simulations of particle trajectories in the shearing sheet usually employ the leapfrog integrator, but this…
Direct N-body simulations and symplectic integrators are effective tools to study the long-term evolution of planetary systems. The Wisdom-Holman (WH) integrator in particular has been used extensively in planetary dynamics as it allows for…
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…
Simulations of the dynamics generated by partial differential equations (PDEs) provide approximate, numerical solutions to initial value problems. Such simulations are ubiquitous in scientific computing, but the correctness of the results…
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…
Long-term stability studies of nonlinear Hamiltonian systems require symplectic integration algorithms which are both fast and accurate. In this paper, we study a symplectic integration method wherein the symplectic map representing the…
The characterization of an interior permanent magnet synchronous machine (IPMSM) requires numerical analysis of the nonlinear magnetic motor model in different load conditions. To obtain the case-specific best machine behavior, a strategy…
In this work, we investigate the regularized solutions and their finite element solutions to the inverse source problems governed by partial differential equations, and establish the stochastic convergence and optimal finite element…
We combine the ideas of atomistic/continuum energy blending and ghost force correction to obtain an energy-based atomistic/continuum coupling scheme which has, for a range of benchmark problems, the same convergence rates as optimal…
Reset control is known to be able to outperform PID and the like linear controllers. However, in motion control systems, quantization can cause severe performance degradation. This paper shows the application of time regularization to…
For the simulation of rectilinearly moving conductors across a magnetic field, the Galer-kin finite element method (GFEM) is generally employed. The inherent instability of GFEM is very often addressed by employing Streamline…
This article considers non-relativistic charged particle dynamics in both static and non-static electromagnetic fields, which are governed by nonseparable, possibly time-dependent Hamiltonians. For the first time, explicit symplectic…
The motion of a simple pendulum in a uniform gravitational field can be described by the solution of a second-order differential equation, nonlinear differential equation. In practice we solve this equation using the small angle…
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…
The point-splitting renormalization method offers a prescription to calculate finite expectation values of quadratic operators constructed from quantum fields in a general curved spacetime. It has been recently shown by Levi and Ori that…
This paper investigates ensemble Kalman inversion (EKI) for variational inverse problems with convex, potentially non-smooth regularization. While deterministic EKI and its Tikhonov-regularized variants have primarily been analyzed for…
There is a growing interest in the conservation of invariants when numerically solving a system of ordinary differential equations. Methods that exactly preserve these quantities in time are known as geometric integrators. In this paper we…
Symplectic integrators are widely used for the study of planetary dynamics and other $N$-body problems. In a study of the outer Solar system, we demonstrate that individual symplectic integrations can yield biased errors in the semi-major…