The Energy-Momentum Method

This paper develops the energy momentum methodJor studying stability and bifurcation of Lagrangian and Hamiltonian systems with symmetry. The method was specifically designed to deal with the stability of rotating structures. The relation with the energy-Casimir method is given and the energy-momentum method is shown to be more general. Stability of rigid body motion is given 10 illustrate the method. Some discussion of its applicability to general rotating systems and block diagonalization is also given

Stress tensors, Riemannian metrics and the alternative descriptions in elasticity

N/

Symplectic-energy-momentum preserving variational integrators

The purpose of this paper is to develop variational integrators for conservative mechanical systems that are symplectic and energy and momentum conserving. To do this, a space–time view of variational integrators is employed and time step adaptation is used to impose the constraint of conservation of energy. Criteria for the solvability of the time steps and some numerical examples are given

Hamiltonian reductions of the one-dimensional Vlasov equation using phase-space moments

We consider Hamiltonian closures of the Vlasov equation using the phase-space moments of the distribution function. We provide some conditions on the closures imposed by the Jacobi identity. We completely solve some families of examples. As a result, we show that imposing that the resulting reduced system preserves the Hamiltonian character of the parent model shapes its phase space by creating a set of Casimir invariants as a direct consequence of the Jacobi identity

Variational Multisymplectic Formulations of Nonsmooth Continuum Mechanics

This paper develops the foundations of the multisymplectic formulation of nonsmooth continuum mechanics. It may be regarded as a PDE generalization of previous techniques that developed a variational approach to collision problems. These methods have already proved of value in computational mechanics, particularly in the development of asynchronous integrators and efficient collision methods. The present formulation also includes solid-fluid interactions and material interfaces and, in addition, lays the groundwork for a treatment of shocks

Optimal trajectory generation in ocean flows

In this paper it is shown that Lagrangian Coherent Structures (LCS) are useful in determining near optimal trajectories for autonomous underwater gliders in a dynamic ocean environment. This opens the opportunity for optimal path planning of autonomous underwater vehicles by studying the global flow geometry via dynamical systems methods. Optimal glider paths were computed for a 2-dimensional kinematic model of an end-point glider problem. Numerical solutions to the optimal control problem were obtained using Nonlinear Trajectory Generation (NTG) software. The resulting solution is compared to corresponding results on LCS obtained using the Direct Lyapunov Exponent method. The velocity data used for these computations was obtained from measurements taken in August, 2000, by HF-Radar stations located around Monterey Bay, CA

Lagrangian coherent structures in n-dimensional systems

Numerical simulations and experimental observations reveal that unsteady fluid systems can be divided into regions of qualitatively different dynamics. The key to understanding transport and stirring is to identify the dynamic boundaries between these almost-invariant regions. Recently, ridges in finite-time Lyapunov exponent fields have been used to define such hyperbolic, almost material, Lagrangian coherent structures in two-dimensional systems. The objective of this paper is to develop and apply a similar theory in higher dimensional spaces. While the separatrix nature of these structures is their most important property, a necessary condition is their almost material nature. This property is addressed in this paper. These results are applied to a model of Rayleigh-Bénard convection based on a three-dimensional extension of the model of Solomon and Gollub

Is SGR 1900+14 a Magnetar?

We present RXTE observations of the soft gamma--ray repeater SGR 1900+14 taken September 4-18, 1996, nearly 2 years before the 1998 active period of the source. The pulsar period (P) of 5.1558199 +/- 0.0000029 s and period derivative (Pdot) of (6.0 +/- 1.0) X 10^-11 s/s measured during the 2-week observation are consistent with the mean Pdot of (6.126 +/- 0.006) X 10^-11 s/s over the time up to the commencement of the active period. This Pdot is less than half that of (12.77 +/- 0.01) X 10^-11 s/s observed during and after the active period. If magnetic dipole radiation were the primary cause of the pulsar spindown, the implied pulsar magnetic field would exceed the critical field of 4.4 X 10^13 G by more than an order of magnitude, and such field estimates for this and other SGRs have been offered as evidence that the SGRs are magnetars, in which the neutron star magnetic energy exceeds the rotational energy. The observed doubling of Pdot, however, would suggest that the pulsar magnetic field energy increased by more than 100% as the source entered an active phase, which seems very hard to reconcile with models in which the SGR bursts are powered by the release of magnetic energy. Because of this, we suggest that the spindown of SGR pulsars is not driven by magnetic dipole radiation, but by some other process, most likely a relativistic wind. The Pdot, therefore, does not provide a measure of the pulsar magnetic field strength, nor evidence for a magnetar.Comment: 14 pages, aasms4 latex, figures 1 & 2 changed, accepted by ApJ letter

Stability Analysis of a Rigid Body with Attached Geometrically Nonlinear Rod by the Energy-Momentum Method

This paper applies the energy-momentum method to the problem of nonlinear stability of relative equilibria of a rigid body with attached flexible appendage in a uniformly rotating state. The appendage is modeled as a geometrically exact rod which allows for finite bending, shearing and twist in three dimensions. Application of the energy-momentum method to this example depends crucially on a special choice of variables in terms of which the second variation block diagonalizes into blocks associated with rigid body modes and internal vibration modes respectively. The analysis yields a nonlinear stability result which states that relative equilibria are nonlinearly stable provided that; (i) the angular velocity is bounded above by the square root of the minimum eigenvalue of an associated linear operator and, (ii) the whole assemblage is rotating about the minimum axis of inertia

Frictional Collisions Off Sharp Objects

This work develops robust contact algorithms capable of dealing with multibody nonsmooth contact geometries for which neither normals nor gap functions can be defined. Such situations arise in the early stage of fragmentation when a number of angular fragments undergo complex collision sequences before eventually scattering. Such situations precludes the application of most contact algorithms proposed to date