VON NEUMANN STABILITY ANALYSIS OF SYMPLECTIC INTEGRATORS APPLIED TO HAMILTONIAN PDEs

Symplectic integration of separable Hamiltonian ordinary and partial differential equations is discussed. A von Neumann analysis is performed to achieve general linear stability criteria for symplectic methods applied to a restricted class of Hamiltonian PDEs. In this treatment, the symplectic step is performed prior to the spatial step, as opposed to the standard approach of spatially discretising the PDE to form a system of Hamiltonian ODEs to which a symplectic integrator can be applied. In this way stability criteria are achieved by considering the spectra of linearised Hamiltonian PDEs rather than spatial step size.     

A Bifurcation Theory for Three-Dimensional Oblique Travelling Gravity-Capillary Water Waves

This article presents a rigorous existence theory for small-amplitude threedimensional travelling water waves. The hydrodynamic problem is formulated as an infinite-dimensional Hamiltonian system in which an arbitrary horizontal spatial direction is the timelike variable. Wave motions that are periodic in a second, different horizontal direction are detected using a centre-manifold reduction technique by which the problem is reduced to a locally equivalent Hamiltonian system with a finite number of degrees of freedom.Dedicated to Klaus Kirchgässner on the occasion of his seventieth birthday     

Geometric integration of the paraxial equation

Evolution of solitary waves in photovoltaic-photorefractive crystal satisfy the paraxial equation. The paraxial equation is transformed into the symplectic structure of the infinite dimensional Hamiltonian system. The symplectic structure of the paraxial equation is discretizated by the symplectic method. The corresponding symplectic scheme preserves conservation of discrete energy which reflects conservation of energy of the paraxial equation. The symplectic scheme is applied to simulate the solitary wave behaviors of the paraxial equation. Evolution of the solitary waves with the different applied electric field and the different photovoltaic fields are investigated.     

Three‐dimensional internal gravity‐capillary waves in finite depth

We consider three‐dimensional inviscid‐irrotational flow in a two‐layer fluid under the effects of gravity and surface tension, where the upper fluid is bounded above by a rigid lid and the lower fluid is bounded below by a flat bottom. We use a spatial dynamics approach and formulate the steady Euler equations as an infinite‐dimensional Hamiltonian system, where an unbounded spatial direction x is considered as a time‐like coordinate. In addition, we consider wave motions that are periodic in another direction z. By analyzing the dispersion relation, we detect several bifurcation scenarios, two of which we study further: a type of 00(is)(iκ₀) resonance and a Hamiltonian Hopf bifurcation. The bifurcations are investigated by performing a center‐manifold reduction, which yields a finite‐dimensional Hamiltonian system. For this finite‐dimensional system, we establish the existence of periodic and homoclinic orbits, which correspond to, respectively, doubly periodic travelling waves and oblique travelling waves with a dark or bright solitary wave profile in the x direction. The former are obtained using a variational Lyapunov‐Schmidt reduction and the latter by first applying a normal form transformation and then studying the resulting canonical system of equations.     

Symplectic invariants for curves and integrable systems in similarity symplectic geometry

In this paper, similarity symplectic geometry for curves is proposed and studied. Explicit expressions of the symplectic invariants, Frenet frame and Frenet formulae for curves in similarity symplectic geometry are obtained by using the equivariant moving frame method. The relationships between the Euclidean symplectic invariants, Frenet frame, Frenet formulae and the similarity symplectic invariants, Frenet frame, Frenet formulae for curves are established. Invariant curve flows in four-dimensional similarity symplectic geometry are also studied. It is shown that certain intrinsic invariant curve flows in four-dimensional similarity symplectic geometry are related to the integrable Burgers and matrix Burgers equations.     

Invariant tori of locally Hamiltonian systems close to conditionally integrable systems

We study the problem of perturbations of quasiperiodic motions in the class of locally Hamiltonian systems. By using methods of the KAM-theory, we prove a theorem on the existence of invariant tori of locally Hamiltonian systems close to conditionally integrable systems. On the basis of this theorem, we investigate the bifurcation of a Cantor set of invariant tori in the case where a Liouville-integrable system is perturbed by a locally Hamiltonian vector field and, simultaneously, the symplectic structure of the phase space is deformed. __________ Translated from Ukrains’kyi Matematychnyi Zhurnal, Vol. 59, No. 1, pp. 71–98, January, 2007.     

Symplectic method for the construction of ergodic measures on invariant submanifolds of nonautonomous hamiltonian systems: Lagrangian manifolds, their structure, and mather homologies

We develop a new approach to the study of properties of ergodic measures for nonautonomous periodic Hamiltonian flows on symplectic manifolds, which are used in many problems of mechanics and mathematical physics. Using Mather’s results on homologies of invariant probability measures that minimize some Lagrangian functionals and the symplectic theory developed by Floer and others for the investigation of symplectic actions and transversal intersections of Lagrangian manifolds, we propose an analog of a Mather-type β-function for the study of ergodic measures associated with nonautonomous Hamiltonian systems on weakly exact symplectic manifolds. Within the framework of the Gromov-Salamon-Zehnder elliptic methods in symplectic geometry, we establish some results on stable and unstable manifolds for hyperbolic invariant sets, which are used in the theory of adiabatic invariants of slowly perturbed integrable Hamiltonian systems. __________ Translated from Ukrains’kyi Matematychnyi Zhurnal, Vol. 58, No. 5, pp. 675–691, May, 2006.     

Actions symplectiques de groupes compacts (Symplectic Actions of Compact Groups)

For any symplectic action of a compact connected group on a compact connected symplectic manifold, we show that the intersection of the Weyl chamber with the image of the moment map is a closed convex polyhedron. This extends Atiyah–Guillemin–Sternberg–Kirwan's convexity theorems to non-Hamiltonian actions. As a consequence, we describe those symplectic actions of a torus which are coisotropic (or multiplicity free), i.e. which have at least one coisotropic orbit: they are the product of an Hamiltonian coisotropic action by an anhamiltonian one. The Hamiltonian coisotropic actions have already been described by Delzant thanks to the convex polyhedron. The anhamiltonian coisotropic actions are actions of a central torus on a symplectic nilmanifold. This text is written as an introduction to the theory of symplectic actions of compact groups since complete proofs of the preliminary classical results are given. An erratum to this article is available at .     

Reissner板弯曲的辛求解体系

基于Reissner板弯曲问题的Hellinger-Reissner变分原理,通过引入对偶变量,导出Reissner板弯曲的Hamilton对偶方程组．从而将该问题导入到哈密顿体系,实现从欧几里德空间向辛几何空间,拉格朗日体系向哈密顿体系的过渡．于是在由原变量及其对偶变量组成的辛几何空间内,许多有效的数学物理方法如分离变量法和本征函数向量展开法等均可直接应用于Reissner板弯曲问题的求解．这里详细求解出Hamilton算子矩阵零本征值的所有本征解及其约当型本征解,给出其具体的物理意义．形成了零本征值本征向量之间的共轭辛正交关系．可以看到,这些零本征值的本征解是Saint-Venant问题所有的基本解,这些解可以张成一个完备的零本征值辛子空间．而非零本征值的本征解是圣维南原理所覆盖的部分．新方法突破了传统半逆解法的限制,有广阔的应用前景．     

基于径向基逼近的KdV方程的无网格辛算法

基于径向基逼近理论,本文为KdV方程构造了一个无网格辛算法.首先借助径向基空间离散Hamilton函数以及Poisson括号,把KdV方程转化成一个有限维的Hamilton系统.然后用辛积分子离散有限维系统,得到辛算法.文章进一步讨论了所构造辛算法的收敛性和误差界.数值例子验证了理论分析.     

Backward error analysis for multi-symplectic integration methods

Summary. A useful method for understanding discretization error in the numerical solution of ODEs is to compare the system of ODEs with the modified equations obtained through backward error analysis, and using symplectic integration for Hamiltonian ODEs provides more incite into the modified equations. In this paper, the ideas of symplectic integration are extended to Hamiltonian PDEs, and this paves the way for the development of a local modified equation analysis solely as a useful diagnostic tool for the study of these types of discretizations. In particular, local conservation laws of energy and momentum are not preserved exactly when symplectic integrators are used to discretize, but the modified equations are used to derive modified conservation laws that are preserved to higher order along the numerical solution. These results are also applied to the nonlinear wave equation. Mathematics Subject Classification (1991):65M06, 65P10, 37K05     

Multiple lie-poisson structures,reductions, and geometric phases for the Maxwell-Bloch travelling wave equations

Summary The real-valued Maxwell-Bloch equations on ℝ³ are investigated as a Hamiltonian dynamical system obtained by applying an S¹ reduction to an invariant subsystem of a dynamical system on ℂ³. These equations on ℝ³ are bi-Hamiltonian and possess several inequivalent Lie-Poisson structures parametrized by classes of orbits in the group SL(2, ℝ). Each Lie-Poisson structure possesses an associated Casimir function. When reduced to level sets of these functions, the motion takes various symplectic forms, from that of the pendulum to that of the Duffing oscillator. The values of the geometric (Hannay-Berry) phases obtained in reconstructing the solutions are found to depend upon the choice of Casimir function, that is, upon the parametrization of the reduced symplectic space.     

Comparison theorems for symplectic systems of difference equations

In the present paper, we prove comparison theorems for symplectic systems of difference equations, which generalize difference analogs of canonical systems of differential equations. We obtain general relations between the number of focal points of conjoined bases of two symplectic systems with matrices W _i and $ \hat W_i $ \hat W_i as well as their corollaries, which generalize well-known comparison theorems for Hamiltonian difference systems. We consider applications of comparison theorems to spectral theory and in the theory of transformations. We obtain a formula for the number of eigenvalues λ of a symplectic boundary value problem on the interval (λ ₁, λ ₂]. For an arbitrary symplectic transformation, we prove a relationship between the numbers of focal points of the conjoined bases of the original and transformed systems. In the case of a constant transformation, we prove a theorem that generalizes the well-known reciprocity principle for discrete Hamiltonian systems.     

光伏光折变晶体中孤立波的数值模拟

高斯光束在光伏光折变晶体中孤立波的演化满足傍轴方程.傍轴方程可以看作无限维Hamil-tonian系统并可以利用辛几何算法进行计算.数值结果表明外加电场和光伏场的强弱和入射高斯光束的振辐对形成稳定的孤立波有显著的影响.傍轴方程的辛几何差分格式能很好地模拟傍轴方程中孤立波的演化行为.     

A multisymplectic framework for classical field theory and the calculus of variations II: space + time decomposition

In a previous paper I laid the foundations of a covariant Hamiltonian framework for the calculus of variations in general. The purpose of the present work is to demonstrate, in the context of classical field theory, how this covariant Hamiltonian formalism may be space + time decomposed. It turns out that the resulting “instantaneous” Hamiltonian formalism is an infinite- dimensional version of Ostrogradski 's theory and leads to the standard symplectic formulation of the initial value problem. The salient features of the analysis are: (i) the instantaneous Hamiltonian formalism does not depend upon the choice of Lepagean equivalent; (ii) the space + time decomposition can be performed either before or after the covariant Legendre transformation has been carried out, with equivalent results; (iii) the instantaneous Hamiltonian can be recovered in natural way from the multisymplectic structure inherent in the theory; and (iv) the space + time split symplectic structure lives on the space of Cauchy data for the evolution equations, as opposed to the space of solutions thereof.     

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.     

Decomposing a new nonlinear differential-difference system under a Bargmann implicit symmetry constraint

Firstly, a hierarchy of integrable lattice equations and its bi-Hamilt-onian structures are established by applying the discrete trace identity. Secondly, under an implicit Bargmann symmetry constraint, every lattice equation in the nonlinear differential-difference system is decomposed by an completely integrable symplectic map and a finite-dimensional Hamiltonian system. Finally, the spatial part and the temporal part of the Lax pairs and adjoint Lax pairs are all constrained as finite dimensional Liouville integrable Hamiltonian systems.     

Variational and Hamiltonian models in the dynamics of magnetizable conducting fluid

The variational model and the Hamiltonian canonical equation of motion are updated using the Lagrangian invariant for three-dimensional unsteady adiabatic flows of magnetizable, ideally conducting, compressible inviscid fluid. The results are applied to derive Hamiltonian noncanonical equations of motion in physically defined variables. Translated from Nelineinye Dinamicheskie Sistemy: Kachestvennyi Analiz i Upravlenie — Sbornik Trudov, No. 2, pp. 44–46, 1994.     

Global geometrical optics method for vector-valued Schrödinger problems

We extend the theory of global geometrical optics method, proposed originally for the linear scalar high-frequency wave-like equations in [Commun. Math. Sci., 2013, 11(1): 105-140], to the more general vector-valued Schrödinger problems in the semi-classical regime. The key ingredient in the global geometrical optics method is a moving frame technique in the phase space. The governing equation is transformed into a new equation but of the same type when expressed in any moving frame induced by the underlying Hamiltonian flow. The classical Wentzel-Kramers-Brillouin (WKB) analysis benefits from this treatment as it maintains valid for arbitrary but fixed evolutionary time. It turns out that a WKB-type function defined merely on the underlying Lagrangian submanifold can be obtained with the help of this moving frame technique, and from which a uniform first-order approximation of the wave field can be derived, even around caustics. The general theory is exemplified by two specific instances. One is the two-level Schrödinger system and the other is the periodic Schrödinger equation. Numerical tests validate the theoretical results.     

On the orbital stability of pendulum-like oscillations and rotations of a symmetric rigid body with a fixed point

We deal with the problem of orbital stability of planar periodic motions of a dynamically symmetric heavy rigid body with a fixed point. We suppose that the center of mass of the body lies in the equatorial plane of the ellipsoid of inertia. Unperturbed periodic motions are planar pendulum-like oscillations or rotations of the body around a principal axis keeping a fixed horizontal position. Local coordinates are introduced in a neighborhood of the unperturbed periodic motion and equations of the perturbed motion are obtained in Hamiltonian form. Regions of orbital instability are established by means of linear analysis. Outside the above-mentioned regions, nonlinear analysis is performed taking into account terms up to degree 4 in the expansion of the Hamiltonian in a neighborhood of unperturbed motion. The nonlinear problem of orbital stability is reduced to analysis of stability of a fixed point of the symplectic map generated by the equations of the perturbed motion. The coefficients of the symplectic map are determined numerically. Rigorous results on the orbital stability or instability of unperturbed motion are obtained by analyzing these coefficients. The orbital stability is investigated analytically in two limiting cases: small amplitude oscillations and rotations with large angular velocities when a small parameter can be introduced.