Related papers: Numerical solution of perturbed Kepler problem usi…
Prior to the recent development of symplectic integrators, the time-stepping operator $\e^{h(A+B)}$ was routinely decomposed into a sum of products of $\e^{h A}$ and $\e^{hB}$ in the study of hyperbolic partial differential equations. In…
We investigate perturbations in the Kepler problem. We offer an overview of the dynamical system using Newtonian, Lagrangian and Hamiltonian Mechanics to build a foundation for analyzing perturbations. We consider the effects of a…
Planar supersymmetric quantum mechanical systems with separable spectral problem in curvilinear coordinates are analyzed in full generality. We explicitly construct the supersymmetric extension of the Euler/Pauli Hamiltonian describing the…
We present an algorithm for the rapid numerical integration of smooth, time-periodic differential equations with small nonlinearity, particularly suited to problems with small dissipation. The emphasis is on speed without compromising…
The Allen-Cahn equation is solved numerically by operator splitting Fourier spectral methods. The basic idea of the operator splitting method is to decompose the original problem into sub-equations and compose the approximate solution of…
In this paper, an efficient parallel splitting method is proposed for the optimal control problem with parabolic equation constraints. The linear finite element is used to approximate the state variable and the control variable in spatial…
An unsteady problem is considered for a space-fractional equation in a bounded domain. A first-order evolutionary equation involves the square root of an elliptic operator of second order. Finite element approximation in space is employed.…
The port-Hamiltonian approach presents an energy-based modeling of dynamical systems with energy-conservative and energy-dissipative parts as well as an interconnection over the so-called ports. In this paper, we apply an operator splitting…
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…
We present two types of meta-algorithm that can greatly improve the accuracy of existing algorithms for integrating the equations of motion of dynamical systems. The first meta-algorithm takes an integrator that is time-symmetric only for…
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,…
The practical application of quantum technologies to chemical problems faces significant challenges, particularly in the treatment of realistic basis sets and the accurate inclusion of electron correlation effects. A direct approach to…
It is argued that, for motion in a central force field, polar reciprocals of trajectories are an elegant alternative to hodographs. The principal advantage of polar reciprocals is that the transformation from a trajectory to its polar…
We apply the method of spectral sequences to study classical problems in analysis. We illustrate the method by finding polynomial vector fields that commute with a given polynomial vector field and finding integrals of polynomial…
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.…
This paper introduces a time-domain combined field integral equation for electromagnetic scattering by a perfect electric conductor. The new equation is obtained by leveraging the quasi-Helmholtz projectors, which separate both the unknown…
The split common fixed-point problem is an inverse problem that consists in finding an element in a fixed-point set such that its image under a bounded linear operator belongs to another fixed-point set. Recently Censor and Segal proposed…
The hodograph of the Kepler-Coulomb problem, that is, the path traced by its velocity vector, is shown to be a circle and then it is used to investigate other properties of the motion. We obtain the configuration space orbits of the problem…
In many physical problems it is not possible to find an exact solution. However, when some parameter in the problem is small, one can obtain an approximate solution by expanding in this parameter. This is the basis of perturbative methods,…
A novel approach is developed for analytic orbit propagation based on asymptotic expansion with respect to a small perturbative acceleration. The method improves upon existing first order asymptotic expansions by leveraging on linear…