High-Order Symplectic and Symmetric Composition Methods for Multi-Frequency and Multi-Dimensional Oscillatory Hamiltonian Systems

The multi-frequency and multi-dimensional adapted Runge-Kutta-Nyström (ARKN) integrators, and multi-frequency and multi-dimensional extended Runge-Kutta-Nyström(ERKN) integrators have been developed to efficiently solve multi-frequency oscillatory Hamiltonian systems. The aim of this paper is to analyze and derive high-order symplectic and symmetric composition methods based on the ARKN integrators and ERKN integrators. We first consider the symplecticity conditions for the multi-frequency and multi-dimensional ARKN integrators. We then analyze the symplecticity of the adjoint integrators of the multi-frequency and multi-dimensional symplectic ARKN integrators and ERKN integrators, respectively. On the basis of the theoretical analysis and by using the idea of composition methods, we derive and propose four new high-order symplectic and symmetric methods for the multi-frequency oscillatory Hamiltonian systems. The numerical results accompanied in this paper quantitatively show the advantage and efficiency of the proposed high-order symplectic and symmetric methods.     

A long-term numerical energy-preserving analysis of symmetric and/or symplectic extended RKN integrators for efficiently solving highly oscillatory Hamiltonian systems

This paper presents a long-term analysis of one-stage extended Runge–Kutta–Nyström (ERKN) integrators for highly oscillatory Hamiltonian systems. We study the long-time numerical energy conservation not only for symmetric integrators but also for symplectic integrators. In the analysis, we neither assume symplecticity for symmetric methods, nor assume symmetry for symplectic methods. It turns out that these both types of integrators have a near conservation of the total and oscillatory energy over a long term. To prove the result for explicit integrators, a relationship between ERKN integrators and trigonometric integrators is established. For the long-term analysis of implicit integrators, the above approach does not work anymore and we use the technology of modulated Fourier expansion. By taking some adaptations of this technology for implicit methods, we derive the modulated Fourier expansion and show the near energy conservation.

     

New open modified Newton Cotes type formulae as multilayer symplectic integrators

New modified open Newton Cotes integrators are introduced in this paper. For the new proposed integrators the connection between these new algorithms, differential methods and symplectic integrators is studied. Much research has been done on one step symplectic integrators and several of them have obtained based on symplectic geometry. However, the research on multistep symplectic integrators is very poor. Zhu et al. [1] studied the well known open Newton Cotes differential methods and they presented them as multilayer symplectic integrators. Chiou and Wu [2] studied the development of multistep symplectic integrators based on the open Newton Cotes integration methods. In this paper we introduce a new open modified numerical method of Newton Cotes type and we present it as symplectic multilayer structure. The new obtained symplectic schemes are applied for the solution of Hamilton’s equations of motion which are linear in position and momentum. An important remark is that the Hamiltonian energy of the system remains almost constant as integration proceeds. We have applied also efficiently the new proposed method to a nonlinear orbital problem and an almost periodic orbital problem.     

Variational integrators of higher order for constrained dynamical systems

The variational integrators presented in [5] are applied to systems with holonomic constraints, yielding constrained higher order variational integrators that are an extension of the constrained Galerkin methods in [4]. The construction of the integrators bases on a discrete version of Hamilton's principle. The inheritance of qualitative properties associated to the solution of the dynamical system to the discrete solution is analysed. Furthermore, the convergence order of the integrators and the computational efficiency is investigated numerically. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Symmetric and symplectic exponential integrators for nonlinear Hamiltonian systems

This letter studies symmetric and symplectic exponential integrators when applied to numerically computing nonlinear Hamiltonian systems. We first establish the symmetry and symplecticity conditions of exponential integrators and then show that these conditions are extensions of the symmetry and symplecticity conditions of Runge–Kutta methods. Based on these conditions, some symmetric and symplectic exponential integrators up to order four are derived. Two numerical experiments are carried out and the results demonstrate the remarkable numerical behaviour of the new exponential integrators in comparison with some symmetric and symplectic Runge–Kutta methods in the literature.     

Volume preserving RK methods for linear systems

1. IntroductionIn recent yearss there has been a great interest in constructing numerical integrationschemes for ODEs in such a way that some qualitative geometrical properties of the solutionof the ODEs are exactly preserved. R.th[ll and Feng Kang[2'31 has proposed symplectic algorithms for Hamiltollian systems, and since then st ruct ure s- preserving me t ho ds fordynamical systems have been systematically developed[4--7]. The symplectic algorithms forHamiltonian systems, the volume-pre…     

基于辛格式的深度哈密尔顿神经网络

HNN是一类基于物理先验学习哈密尔顿系统的神经网络.本文通过误差分析解释使用不同积分器作为超参数对HNN的影响.如果我们把网络目标定义为在任意训练集上损失为零的映射,那么传统的积分器无法保证HNN存在网络目标.我们引进反修正方程,并严格证明基于辛格式的HNN具有网络目标,且它与原哈密尔顿量之差依赖于数值格式的精度.数值实验表明,由辛HNN得到的哈密尔顿系统的相流不能精确保持原哈密尔顿量,但保持网络目标;网络目标在训练集、测试集上的损失远小于原哈密尔顿量的损失;在预测问题上辛HNN较非辛HNN具备更强大的泛化能力和更高的精度.因此,辛格式对于HNN是至关重要的.     

Second-order accurate projective integrators for multiscale problems

We introduce new projective versions of second-order accurate Runge–Kutta and Adams–Bashforth methods, and demonstrate their use as outer integrators in solving stiff differential systems. An important outcome is that the new outer integrators, when combined with an inner telescopic projective integrator, can result in fully explicit methods with adaptive outer step size selection and solution accuracy comparable to those obtained by implicit integrators. If the stiff differential equations are not directly available, our formulations and stability analysis are general enough to allow the combined outer–inner projective integrators to be applied to legacy codes or perform a coarse-grained time integration of microscopic systems to evolve macroscopic behavior, for example.     

High-Order Symplectic Partitioned Lie Group Methods

In this article, a unified approach to obtain symplectic integrators on \(T^{*}G\) from Lie group integrators on a Lie group \(G\) is presented. The approach is worked out in detail for symplectic integrators based on Runge–Kutta–Munthe-Kaas methods and Crouch–Grossman methods. These methods can be interpreted as symplectic partitioned Runge–Kutta methods extended to the Lie group setting in two different ways. In both cases, we show that it is possible to obtain symplectic integrators of arbitrarily high order by this approach.     

Higher order variational integrators in optimal control theory

Higher order variational integrators are analyzed and applied to optimal control problems posed with mechanical systems. First, we derive two different kinds of high order variational integrators based on different dimensions of the underlying approximation space. While the first well-known integrator is equivalent to a symplectic partitioned Runge-Kutta method, the second integrator, denoted as symplectic Galerkin integrator, yields a method which in general, cannot be written as a standard symplectic Runge-Kutta scheme [1]. Furthermore, we use these integrators for the discretization of optimal control problems. By analyzing the adjoint systems of the optimal control problem and its discretized counterpart, we prove that for these particular integrators optimization and discretization commute [2]. This property guarantees that the accuracy is preserved for the adjoint system which is also referred to as the Covector Mapping Principle. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Exponential Rosenbrock integrators for option pricing

In this paper, we are concerned with the time integration of differential equations modeling option pricing. In particular, we consider the Black-Scholes equation for American options. As an alternative to existing methods, we present exponential Rosenbrock integrators. These integrators require the evaluation of the exponential and related functions of the Jacobian matrix. The resulting methods have good stability properties. They are fully explicit and do not require the numerical solution of linear systems, in contrast to standard integrators. We have implemented some numerical experiments in Matlab showing the reliability of the new method.     

Consistent time integration for transient three-dimensional contact using the NTS-method

The present work aims to investigate stable time integrators for large deformation contact problems within the framework of the well known node-to-segment (NTS)-method. For this kind of problem, standard time integrators fail to conserve the total energy of the system. To remedy this drawback, we combine a mixed method with the concept of a discrete gradient applied to the aforementioned NTS-method. In the context of nonlinear elastodynamics stable integrators for ordinary differential equations have been extensively developed and investigated during the last two decades. For contact problems, energy consistent integrators have been developed for the NTS-method (see e.g. Ref. [1]). (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Stochastic integrators indexed by a multi-dimensional parameter

Summary Multi-parameter stochastic integrators are described and classified according to directions of integrability. Sets of directions are distinguished precisely by the combinatorial dimension of corresponding fractional Cartesian products. The main theorem establishes existence of stochastic processes which are integrators in prescribed directions but not others.     

Hamilton–Pontryagin Integrators on Lie Groups Part I: Introduction and Structure-Preserving Properties

In this paper, structure-preserving time-integrators for rigid body-type mechanical systems are derived from a discrete Hamilton–Pontryagin variational principle. From this principle, one can derive a novel class of variational partitioned Runge–Kutta methods on Lie groups. Included among these integrators are generalizations of symplectic Euler and Störmer–Verlet integrators from flat spaces to Lie groups. Because of their variational design, these integrators preserve a discrete momentum map (in the presence of symmetry) and a symplectic form. In a companion paper, we perform a numerical analysis of these methods and report on numerical experiments on the rigid body and chaotic dynamics of an underwater vehicle. The numerics reveal that these variational integrators possess structure-preserving properties that methods designed to preserve momentum (using the coadjoint action of the Lie group) and energy (for example, by projection) lack.     

Addition in a class of nonlinear Stieltjes integrators

A product integral formula is established for the generation of an evolution system by the sum of two Stieltjes integrators. The class of Stieltjes integrators considered is the class recently introduced and studied by J.V. Herod.     

ON ENERGY CONSERVATION BY TRIGONOMETRIC INTEGRATORS IN THE LINEAR CASE WITH APPLICATION TO WAVE EQUATIONS

Trigonometric integrators for oscillatory linear Hamiltonian differential equations are considered. Under a condition of Hairer & Lubich on the filter functions in the method, a modified energy is derived that is exactly preserved by trigonometric integrators. This implies and extends a known result on all-time near-conservation of energy. The extension can be applied to linear wave equations.     

Closed Newton–Cotes trigonometrically-fitted formulae of high order for long-time integration of orbital problems

The connection between closed Newton–Cotes, trigonometrically-fitted differential methods and symplectic integrators is studied in this paper. Several one-step symplectic integrators have been obtained based on symplectic geometry, as is shown in the literature. However, the study of multi-step symplectic integrators is very limited. The well-known open Newton–Cotes differential methods are presented as multilayer symplectic integrators by Zhu et al. [W. Zhu, X. Zhao, Y. Tang, Journal of Chem. Phys. 104 (1996), 2275]. The construction of multi-step symplectic integrators based on the open Newton–Cotes integration methods is investigated by Chiou and Wu [J.C. Chiou, S.D. Wu, Journal of Chemical Physics 107 (1997), 6894]. The closed Newton–Cotes formulae are studied in this paper and presented as symplectic multilayer structures. We also develop trigonometrically-fitted symplectic methods which are based on the closed Newton–Cotes formulae. We apply the symplectic schemes in order to solve Hamilton’s equations of motion which are linear in position and momentum. We observe that the Hamiltonian energy of the system remains almost constant as the integration proceeds. Finally we apply the new developed methods to an orbital problem in order to show the efficiency of this new methodology.     

Weak local linear discretizations for stochastic differential equations: Convergence and numerical schemes

Weak local linear (WLL) discretizations are playing an increasing role in the construction of effective numerical integrators and inference methods for stochastic differential equations (SDEs) with additive noise. However, due to limitations in the existing numerical implementations of WLL discretizations, the resulting integrators and inference methods have either been restricted to particular classes of autonomous SDEs or showed low computational efficiency. Another limitation is the absence of a systematic theoretical study of the rate of convergence of the WLL discretizations and numerical integratos. This task is the main purpose of the present paper. A second goal is introducing a new WLL scheme that overcomes the numerical limitations mentioned above. Additionally, a comparative analysis between the new WLL scheme and some conventional weak integrators is also presented.     

General techniques for constructing variational integrators

The numerical analysis of variational integrators relies on variational error analysis, which relates the order of accuracy of a variational integrator with the order of approximation of the exact discrete Lagrangian by a computable discrete Lagrangian. The exact discrete Lagrangian can either be characterized variationally, or in terms of Jacobi’s solution of the Hamilton-Jacobi equation. These two characterizations lead to the Galerkin and shooting constructions for discrete Lagrangians, which depend on a choice of a numerical quadrature formula, together with either a finite-dimensional function space or a one-step method. We prove that the properties of the quadrature formula, finite-dimensional function space, and underlying one-step method determine the order of accuracy and momentum-conservation properties of the associated variational integrators. We also illustrate these systematic methods for constructing variational integrators with numerical examples.