Related papers: Symplectic fourth-order maps for the collisional N…
A series of robust and optimal mixed methods based on two mixed formulations of the fourth-order elliptic singular perturbation problem are developed in this paper. First, a mixed method based on a second-order system is proposed without…
In the framework of multi-body dynamics, successive encounters with a third body, even if well outside of its sphere of influence, can noticeably alter the trajectory of a spacecraft. Examples of these effects have already been exploited by…
We suggest a numerical integration procedure for solving the equations of motion of certain classical spin systems which preserves the underlying symplectic structure of the phase space. Such symplectic integrators have been successfully…
In this paper we study the performance of a symplectic numerical integrator based on the splitting method. This method is applied to a subtle problem i.e. higher order resonance of the elastic pendulum. In order to numerically study the…
We present a novel hierarchical formulation of the fourth-order forward symplectic integrator and its numerical implementation in the GPU-accelerated direct-summation N-body code FROST. The new integrator is especially suitable for…
In this work, we introduce a quadratically convergent and dynamically consistent integrator specifically designed for the replicator dynamics. The proposed scheme combines a two-stage rational approximation with a normalization step to…
Multibody simulation with frictional contact has been a challenging subject of research for the past thirty years. Rigid-body assumptions are commonly used to approximate the physics of contact, and together with Coulomb friction, lead to…
The dynamic equation of mass point in rotating coordinates is governed by Coriolis and centrifugal force, besides a corotating potential relative to frame. Such a system is no longer a canonical Hamiltonian system so that the construction…
We propose a new primal-dual splitting method for solving composite inclusions involving Lipschitzian, and parallel-sum-type monotone operators. Our approach extends the framework in \cite{Siopt4} to a more general class of monotone…
Modelling the cosmic large-scale structure can be done through numerical N-body simulations or by using perturbation theory. Here, we present an N-body approach that effectively implements a multi-step forward model based on Lagrangian…
A fixed time-step variational integrator cannot preserve momentum, energy, and symplectic form simultaneously for nonintegrable systems. This barrier can be overcome by treating time as a discrete dynamic variable and deriving adaptive…
In this paper, we consider the coupled N/TH problem, in which the termination criterion for the neutronics iteration adopts an adaptive tolerance with respect to the fuel temperature residual at each Picard iteration. We refer to this…
We derive a new criterion for estimating characteristic dynamical timescales in N-body simulations. The criterion uses the second, third, and fourth derivatives of particle positions: acceleration, jerk, and snap. It can be used for…
The paper considers a linear system of Boltzmann transport equations modelling the evolution of three species of particles, photons, electrons and positrons. The system is coupled because of the collision term (an integral operator). The…
We present several methods, which utilize symplectic integration techniques based on two and three part operator splitting, for numerically solving the equations of motion of the disordered, discrete nonlinear Schr\"odinger (DDNLS)…
We apply theta body relaxations to the $K_i$-cover problem and show polynomial time solvability for certain classes of graphs. In particular, we give an effective relaxation where all $K_i$-$p$-hole facets are valid, and study its relation…
In this article we develop a numerical scheme to deal with interfaces between touching numerical grids when solving Schr\"o{}dinger equation. In order to pass the information among grids we use the values of the fields only at the contact…
$N$-body simulations study the dynamics of $N$ particles under the influence of mutual long-distant forces such as gravity. In practice, $N$-body codes will violate Newton's third law if they use either an approximate Poisson solver or…
Efficient simulation of many-body quantum systems is central to advances in physics, chemistry, and quantum computing, with a key question being whether the simulation cost scales polynomially with the system size. In this work, we analyze…
We present a framework for the simulation of rigid and deformable bodies in the presence of contact and friction. Our method is based on a non-smooth Newton iteration that solves the underlying nonlinear complementarity problems (NCPs)…