Regular Reduction of Controlled Hamiltonian System with Symplectic Structure and Symmetry

In this paper, our goal is to study the regular reduction theory of regular controlled Hamiltonian (RCH) systems with symplectic structure and symmetry, and this reduction is an extension of regular symplectic reduction theory of Hamiltonian systems under regular controlled Hamiltonian equivalence conditions. Thus, in order to describe uniformly RCH systems defined on a cotangent bundle and on the regular reduced spaces, we first define a kind of RCH systems on a symplectic fiber bundle. Then introduce regular point and regular orbit reducible RCH systems with symmetry by using momentum map and the associated reduced symplectic forms. Moreover, we give regular point and regular orbit reduction theorems for RCH systems to explain the relationships between RpCH-equivalence, RoCH-equivalence for reducible RCH systems with symmetry and RCH-equivalence for associated reduced RCH systems. Finally, as an application we regard rigid body and heavy top as well as them with internal rotors as the regular point reducible RCH systems on the rotation group $textmd{SO}(3)$ and on the Euclidean group $textmd{SE}(3)$,as well as on their generalizations, respectively, and discuss their RCH-equivalence. We also describe the RCH system and RCH-equivalence from the viewpoint of port Hamiltonian system with a symplectic structure.

Space-time FLAVORS: finite difference, multisymlectic, and pseudospectral integrators for multiscale PDEs

We present a new class of integrators for stiff PDEs. These integrators are generalizations of FLow AVeraging integratORS (FLAVORS) for stiff ODEs and SDEs introduced in [Tao, Owhadi and Marsden 2010] with the following properties: (i) Multiscale: they are based on flow averaging and have a computational cost determined by mesoscopic steps in space and time instead of microscopic steps in space and time; (ii) Versatile: the method is based on averaging the flows of the given PDEs (which may have hidden slow and fast processes). This bypasses the need for identifying explicitly (or numerically) the slow variables or reduced effective PDEs; (iii) Nonintrusive: A pre-existing numerical scheme resolving the microscopic time scale can be used as a black box and easily turned into one of the integrators in this paper by turning the large coefficients on over a microscopic timescale and off during a mesoscopic timescale; (iv) Convergent over two scales: strongly over slow processes and in the sense of measures over fast ones; (v) Structure-preserving: for stiff Hamiltonian PDEs (possibly on manifolds), they can be made to be multi-symplectic, symmetry-preserving (symmetries are group actions that leave the system invariant) in all variables and variational.

Multi-Dirac Structures and Hamilton-Pontryagin Principles for Lagrange-Dirac Field Theories

The purpose of this paper is to define the concept of multi-Dirac structures and to describe their role in the description of classical field theories. We begin by outlining a variational principle for field theories, referred to as the Hamilton-Pontryagin principle, and we show that the resulting field equations are the Euler-Lagrange equations in implicit form. Secondly, we introduce multi-Dirac structures as a graded analog of standard Dirac structures, and we show that the graph of a multisymplectic form determines a multi-Dirac structure. We then discuss the role of multi-Dirac structures in field theory by showing that the implicit field equations obtained from the Hamilton-Pontryagin principle can be described intrinsically using multi-Dirac structures. Furthermore, we show that any multi-Dirac structure naturally gives rise to a multi-Poisson bracket. We treat the case of field theories with nonholonomic constraints, showing that the integrability of the constraints is equivalent to the integrability of the underlying multi-Dirac structure. We finish with a number of illustrative examples, including time-dependent mechanics, nonlinear scalar fields and the electromagnetic field.

Temperature and Friction Accelerated Sampling of Boltzmann-Gibbs Distribution

This paper is concerned with tuning friction and temperature in Langevin dynamics for fast sampling from the canonical ensemble. We show that near-optimal acceleration is achieved by choosing friction so that the local quadratic approximation of the Hamiltonian is a critical damped oscillator. The system is also over-heated and cooled down to its final temperature. The performances of different cooling schedules are analyzed as functions of total simulation time.

From efficient symplectic exponentiation of matrices to symplectic integration of high-dimensional Hamiltonian systems with slowly varying quadratic stiff potentials

We present a multiscale integrator for Hamiltonian systems with slowly varying quadratic stiff potentials that uses coarse timesteps (analogous to what the impulse method uses for constant quadratic stiff potentials). This method is based on the highly-non-trivial introduction of two efficient symplectic schemes for exponentiations of matrices that only require O(n) matrix multiplications operations at each coarse time step for a preset small number n. The proposed integrator is shown to be (i) uniformly convergent on positions; (ii) symplectic in both slow and fast variables; (iii) well adapted to high dimensional systems. Our framework also provides a general method for iteratively exponentiating a slowly varying sequence of (possibly high dimensional) matrices in an efficient way.

Structure preserving Stochastic Impulse Methods for stiff Langevin systems with a uniform global error of order 1 or 1/2 on position

Impulse methods are generalized to a family of integrators for Langevin systems with quadratic stiff potentials and arbitrary soft potentials. Uniform error bounds (independent from stiff parameters) are obtained on integrated positions allowing for coarse integration steps. The resulting integrators are explicit and structure preserving (quasi-symplectic for Langevin systems).

The Dynamics of a Rigid Body in Potential Flow with Circulation

We consider the motion of a two-dimensional body of arbitrary shape in a planar irrotational, incompressible fluid with a given amount of circulation around the body. We derive the equations of motion for this system by performing symplectic reduction with respect to the group of volume-preserving diffeomorphisms and obtain the relevant Poisson structures after a further Poisson reduction with respect to the group of translations and rotations. In this way, we recover the equations of motion given for this system by Chaplygin and Lamb, and we give a geometric interpretation for the Kutta-Zhukowski force as a curvature-related effect. In addition, we show that the motion of a rigid body with circulation can be understood as a geodesic flow on a central extension of the special Euclidian group SE(2), and we relate the cocycle in the description of this central extension to a certain curvature tensor.