Related papers: High Order Phase-fitted Discrete Lagrangian Integr…
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…
Network interactions between dynamical units are often subject to time delay. We develop a phase reduction method for delay-coupled oscillator networks. The method is based on rewriting the delay-differential equation as an ordinary…
Variational integrators for Lagrangian dynamical systems provide a systematic way to derive geometric numerical methods. These methods preserve a discrete multisymplectic form as well as momenta associated to symmetries of the Lagrangian…
On the basis of the previous work by Tang \& Zhang (Appl. Math. Comput. 323, 2018, p. 204--219), in this paper we present a more effective way to construct high-order symplectic integrators for solving second order Hamiltonian equations.…
We develop an algebraic approach for finding the eigenfunctions of a large class of few and many-body Hamiltonians, in one and higher dimensions, having linear spectra. The method presented enables one to exactly map these interacting…
Correlated many-body problems ubiquitously appear in various fields of physics such as condensed matter physics, nuclear physics, and statistical physics. However, due to the interplay of the large number of degrees of freedom, it is…
We analyze a fully discrete finite element numerical scheme for the Cahn-Hilliard-Stokes-Darcy system that models two-phase flows in coupled free flow and porous media. To avoid a well-known difficulty associated with the coupling between…
This paper addresses the problem of efficiently computing higher-order variational integrators in simulation and trajectory optimization of mechanical systems as those often found in robotic applications. We develop $O(n)$ algorithms to…
Strongly interacting electrons in solids are generically described by Hubbardtype models, and the impact of solar light can be modeled by an additional time-dependence. This yields a finite dimensional system of ordinary differential…
An explicit high order semi-Lagrangian method is developed for solving Lagrangian transport equations in Eulerian-Lagrangian formulations. To ensure a semi-Lagrangian approximation that is consistent with an explicit Eulerian, discontinuous…
This paper describes a fourth-order integration algorithm for the gravitational N-body problem based on discrete Lagrangian mechanics. When used with shared timesteps, the algorithm is momentum conserving and symplectic. We generalize the…
In the present work, a kind of trigonometric collocation methods based on Lagrange basis polynomials is developed for effectively solving multi-frequency oscillatory second-order differential equations…
In this paper, we consider the task of efficiently computing the numerical solution of evolutionary complex Ginzburg--Landau equations on Cartesian product domains with homogeneous Dirichlet/Neumann or periodic boundary conditions. To this…
We consider the numerical integration of non-autonomous separable parabolic equations using high order splitting methods with complex coefficients (methods with real coefficients of order greater than two necessarily have negative…
In this paper, we study the Lagrangian functions for a class of second-order differential systems arising from physics. For such systems, we present necessary and sufficient conditions for the existence of Lagrangian functions. Based on the…
Magnetic quadrupoles are essential components of particle accelerators like the Large Hadron Collider. In order to study numerically the stability of the particle beam crossing a quadrupole, a large number of particle revolutions in the…
We design a novel, exactly energy-conserving implicit non-symplectic integration method for an eight-dimensional Hamiltonian system with four degrees of freedom. In our algorithm, each partial derivative of the Hamiltonian with respect to…
Many combinatorial optimization problems can be mapped to finding the ground states of the corresponding Ising Hamiltonians. The physical systems that can solve optimization problems in this way, namely Ising machines, have been attracting…
We introduce a finite element method for numerical upscaling of second order elliptic equations with highly heterogeneous coefficients. The method is based on a mixed formulation of the problem and the concepts of the domain decomposition…
We show that when time-reversible symplectic algorithms are used to solve periodic motions, the energy error after one period is generally two orders higher than that of the algorithm. By use of correctable algorithms, we show that the…