A continuous/discrete fractional Noether’s theorem

We prove a fractional Noether’s theorem for fractional Lagrangian systems invariant under a symmetry group both in the continuous and discrete cases. This provides an explicit conservation law (first integral) given by a closed formula which can be algorithmically implemented. In the discrete case, the conservation law is moreover computable in a finite number of steps.     

Hamiltonian systems with a discrete symmetry

If $\text{?: M → M}$ is an antisymplectic involution of a symplectic manifold M then the fixed set of ? is a Lagrangian submanifold L ? M. Moreover there exist cotangent bundle coordinates in a neighborhood of L in M such that ? in these coordinates maps a covector into its negative. Thus classical examples which have a discrete symmetry such as the restricted three-body problems are locally like a reversible system.     

Node-pair finite volume/finite element schemes for the Euler equation in cylindrical and spherical coordinates

A numerical scheme is presented for the solution of the compressible Euler equations in both cylindrical and spherical coordinates. The unstructured grid solver is based on a mixed finite volume/finite element approach. Equivalence conditions linking the node-centered finite volume and the linear Lagrangian finite element scheme over unstructured grids are reported and used to devise a common framework for solving the discrete Euler equations in both the cylindrical and the spherical reference systems. Numerical simulations are presented for the explosion and implosion problems with spherical symmetry, which are solved in both the axial–radial cylindrical coordinates and the radial–azimuthal spherical coordinates. Numerical results are found to be in good agreement with one-dimensional simulations over a fine mesh.     

Variational and Geometric Structures of Discrete Dirac Mechanics

In this paper, we develop the theoretical foundations of discrete Dirac mechanics, that is, discrete mechanics of degenerate Lagrangian/Hamiltonian systems with constraints. We first construct discrete analogues of Tulczyjew’s triple and induced Dirac structures by considering the geometry of symplectic maps and their associated generating functions. We demonstrate that this framework provides a means of deriving discrete Lagrange–Dirac and nonholonomic Hamiltonian systems. In particular, this yields nonholonomic Lagrangian and Hamiltonian integrators. We also introduce discrete Lagrange–d’Alembert–Pontryagin and Hamilton–d’Alembert variational principles, which provide an alternative derivation of the same set of integration algorithms. The paper provides a unified treatment of discrete Lagrangian and Hamiltonian mechanics in the more general setting of discrete Dirac mechanics, as well as a generalization of symplectic and Poisson integrators to the broader category of Dirac integrators.     

Lagrangian reduction and the double spherical pendulum

This paper studies the stability and bifurcations of the relative equilibrium of the double spherical pendulum, which has the circle as its symmetry group. The example as well as others with nonabelian symmetry groups, such as the rigid body, illustrate some useful general theory about Lagrangian reduction. In particular, we establish a satisfactory global theory of Lagrangian reduction that is consistent with the classical local Routh theory for systems with an abelian symmetry group.Dedicated to Professor Klaus Kirchgässner on the occasion of his 60th birthdayResearch partially supported by a Humboldt award at the Universität Hamburg and by DOE Contract DE-FGO3-88ER25064.     

Geometric Integrators for Higher-Order Variational Systems and Their Application to Optimal Control

Numerical methods that preserve geometric invariants of the system, such as energy, momentum or the symplectic form, are called geometric integrators. In this paper we present a method to construct symplectic-momentum integrators for higher-order Lagrangian systems. Given a regular higher-order Lagrangian \(L:T^{(k)}Q\rightarrow {\mathbb {R}}\) with \(k\ge 1\), the resulting discrete equations define a generally implicit numerical integrator algorithm on \(T^{(k-1)}Q\times T^{(k-1)}Q\) that approximates the flow of the higher-order Euler–Lagrange equations for L. The algorithm equations are called higher-order discrete Euler–Lagrange equations and constitute a variational integrator for higher-order mechanical systems. The general idea for those variational integrators is to directly discretize Hamilton’s principle rather than the equations of motion in a way that preserves the invariants of the original system, notably the symplectic form and, via a discrete version of Noether’s theorem, the momentum map. We construct an exact discrete Lagrangian \(L_d^e\) using the locally unique solution of the higher-order Euler–Lagrange equations for L with boundary conditions. By taking the discrete Lagrangian as an approximation of \(L_d^e\), we obtain variational integrators for higher-order mechanical systems. We apply our techniques to optimal control problems since, given a cost function, the optimal control problem is understood as a second-order variational problem.     

An exact algorithm for the capacitated facility location problems with single sourcing

Facility location problems are often encountered in many areas such as distribution, transportation and telecommunication. We describe a new solution approach for the capacitated facility location problem in which each customer is served by a single facility. An important class of heuristic solution methods for these problems are Lagrangian heuristics which have been shown to produce high quality solutions and at the same time be quite robust. A primal heuristic, based on a repeated matching algorithm which essentially solves a series of matching problems until certain convergence criteria are satisfied, is incorporated into the Lagrangian heuristic. Finally, a branch-and-bound method, based on the Lagrangian heuristic is developed, and compared computationally to the commercial code CPLEX. The computational results indicate that the proposed method is very efficient.     

A cell-centered lagrangian scheme in two-dimensional cylindrical geometry

A new Lagrangian cell-centered scheme for two-dimensional compressible flows in planar geometry is proposed by Maire et al.The main new feature of the algorithm is that the vertex velocities and the numerical fluxes through the cell interfaces are all evaluated in a coherent manner contrary to standard approaches.In this paper the method introduced by Maire et al.is extended for the equations of Lagrangian gas dynamics in cylindrical symmetry.Two different schemes are proposed,whose difference is that one uses volume weighting and the other area weighting in the discretization of the momentum equation.In the both schemes the conservation of total energy is ensured,and the nodal solver is adopted which has the same formulation as that in Cartesian coordinates.The volume weighting scheme preserves the momentum conservation and the area-weighting scheme preserves spherical symmetry.The numerical examples demonstrate our theoretical considerations and the robustness of the new method.     

A sampling theorem on homogeneous manifolds

We consider a generalization of entire functions of spherical exponential type and Lagrangian splines on manifolds. An analog of the Paley-Wiener theorem is given. We also show that every spectral entire function on a manifold is uniquely determined by its values on some discrete sets of points.

The main result of the paper is a formula for reconstruction of spectral entire functions from their values on discrete sets using Lagrangian splines.

     

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.     

Locally symmetric minimal affine Lagrangian surfaces in C2

One of the basic facts known in the theory of minimal Lagrangian surfaces is that a minimal Lagrangian surface of constant curvature in C ² must be totally geodesic. In affine geometry the constancy of curvature corresponds to the local symmetry of a connection. In Opozda (Geom. Dedic. 121:155–166, 2006), we proposed an affine version of the theory of minimal Lagrangian submanifolds. In this paper we give a local classification of locally symmetric minimal affine Lagrangian surfaces in C ². Only very few of surfaces obtained in the classification theorems are Lagrangian in the sense of metric (pseudo-Riemannian) geometry.     

A connection-theoretic approach to reduction of second-order dynamical systems with symmetry

We deal with reduction of Lagrangian systems that are invariant under the action of the symmetry group. Unlike the bulk of the literature we do not rely on methods coming from the calculus of variations. Our method is based on the geometrical analysis of regular Lagrangian systems, where solutions of the Euler-Lagrange equations are interpreted as integral curves of the associated second-order differential equation field. In particular, we explain so-called Lagrange-Poincaré reduction [1] and Routh reduction [3] from the viewpoint of that vector field. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Lagrangian bounds in multiextremal polynomial and discrete optimization problems

Many polynomial and discrete optimization problems can be reduced to multiextremal quadratic type models of nonlinear programming. For solving these problems one may use Lagrangian bounds in combination with branch and bound techniques. The Lagrangian bounds may be improved for some important examples by adding in a model the so-called superfluous quadratic constraints which modify Lagrangian bounds. Problems of finding Lagrangian bounds as a rule can be reduced to minimization of nonsmooth convex functions and may be successively solved by modern methods of nondifferentiable optimization. This approach is illustrated by examples of solving polynomial-type problems and some discrete optimization problems on graphs.     

Locally symmetric minimal affine Lagrangian surfaces in C<Superscript>2</Superscript>

One of the basic facts known in the theory of minimal Lagrangian surfaces is that a minimal Lagrangian surface of constant curvature in C ² must be totally geodesic. In affine geometry the constancy of curvature corresponds to the local symmetry of a connection. In Opozda (Geom. Dedic. 121:155–166, 2006), we proposed an affine version of the theory of minimal Lagrangian submanifolds. In this paper we give a local classification of locally symmetric minimal affine Lagrangian surfaces in C ². Only very few of surfaces obtained in the classification theorems are Lagrangian in the sense of metric (pseudo-Riemannian) geometry. The research supported by the KBN grant 1 PO3A 034 26.     

Conservation laws via the partial Lagrangian and group invariant solutions for radial and two-dimensional free jets

The conservation laws for Prandtl’s boundary layer equations for an incompressible fluid governing the flow in radial and two-dimensional jets are investigated. For both radial and two-dimensional jets the partial Lagrangian method is used to derive conservation laws for the system of two differential equations for the velocity components. The Lie point symmetries are calculated for both cases and a symmetry is associated with the conserved vector that is used to establish the conserved quantity for the jet. This associated symmetry is then used to derive the group invariant solution for the system governing the flow in the free jet.     

The inverse problem of the calculus of variations for scalar fourth-order ordinary differential equations

A simple invariant characterization of the scalar fourth-order ordinary differential equations which admit a variational multiplier is given. The necessary and sufficient conditions for the existence of a multiplier are expressed in terms of the vanishing of two relative invariants which can be associated with any fourth-order equation through the application of Cartan's equivalence method. The solution to the inverse problem for fourth-order scalar equations provides the solution to an equivalence problem for second-order Lagrangians, as well as the precise relationship between the symmetry algebra of a variational equation and the divergence symmetry algebra of the associated Lagrangian.

     

A nondisplaceable Lagrangian torus in T*S2

We show that the Lagrangian torus in the cotangent bundles of the 2‐sphere obtained by applying the geodesic flow to the unit circle in a fiber is not displaceable by computing its Lagrangian Floer homology. The computation is based on a symmetry argument. © 2007 Wiley Periodicals, Inc.     

Oscillatory solutions of fourth order conservative systems via the Conley index

In this paper we investigate periodic solutions of second order Lagrangian systems which oscillate around equilibrium points of center type. The main ingredients are the discretization of second order Lagrangian systems that satisfy the twist property and the theory of discrete braid invariants developed by Ghrist et al. (2003) [5]. The problem with applying this topological theory directly is that the braid types in our analysis are so-called improper. This implies that the braid invariants do not entirely depend on the topology: the relevant braid classes are non-isolating neighborhoods of the flow, so that their Conley index is not universal. In first part of this paper we develop the theory of the braid invariant for improper braid classes and in the second part this theory is applied to second order Lagrangian system and in particular to the Swift–Hohenberg equation.     

Chaotic Attractors with Discrete Planar Symmetries

Chaotic behavior is known to be compatible with symmetry and illustrations are constructed using functions equivariant with respect to the desired symmetries. Earlier investigations determined families of equivariant functions for a few of the discrete symmetry groups in the plane; those results are extended to all the discrete symmetry groups of the plane. This includes consideration of the all the frieze and two-dimensional crystallographic groups. © 1998 Elsevier Science Ltd. All rights reserved.