Navier-Stokes-alpha model: LES equations with nonlinear dispersion

We present a framework for discussing LES equations with nonlinear dispersion. In this framework, we discuss the properties of the nonlinearly dispersive Navier-Stokes-alpha model of incompressible fluid turbulence --- also called the viscous Camassa-Holm equations and the LANS equations in the literature --- in comparison with the corresponding properties of large eddy simulation (LES) equations obtained via the approximate-inverse approach. In this comparison, we identify the spatially filtered NS-alpha equations with a class of generalized LES similarity models. Applying a certain approximate inverse to this class of LES models restores the Kelvin circulation theorem for the defiltered velocity and shows that the NS-alpha model describes the dynamics of the defiltered velocity for this class of generalized LES similarity models. We also show that the subgrid scale forces in the NS-alpha model transform covariantly under Galilean transformations and under a change to a uniformly rotating reference frame. Finally, we discuss in the spectral formulation how the NS-alpha model retains the local interactions among the large scales, retains the nonlocal sweeping effects of large scales on small scales, yet attenuates the local interactions of the small scales amongst themselves.Comment: 15 pages, no figures, Special LES volume of ERCOFTAC bulletin, to appear in 200

Continuous and discrete Clebsch variational principles

The Clebsch method provides a unifying approach for deriving variational principles for continuous and discrete dynamical systems where elements of a vector space are used to control dynamics on the cotangent bundle of a Lie group \emph{via} a velocity map. This paper proves a reduction theorem which states that the canonical variables on the Lie group can be eliminated, if and only if the velocity map is a Lie algebra action, thereby producing the Euler-Poincar\'e (EP) equation for the vector space variables. In this case, the map from the canonical variables on the Lie group to the vector space is the standard momentum map defined using the diamond operator. We apply the Clebsch method in examples of the rotating rigid body and the incompressible Euler equations. Along the way, we explain how singular solutions of the EP equation for the diffeomorphism group (EPDiff) arise as momentum maps in the Clebsch approach. In the case of finite dimensional Lie groups, the Clebsch variational principle is discretised to produce a variational integrator for the dynamical system. We obtain a discrete map from which the variables on the cotangent bundle of a Lie group may be eliminated to produce a discrete EP equation for elements of the vector space. We give an integrator for the rotating rigid body as an example. We also briefly discuss how to discretise infinite-dimensional Clebsch systems, so as to produce conservative numerical methods for fluid dynamics

A study of human performance in a rotating environment

Consideration is given to the lack of sufficient data relative to the response of man to the attendant oculovestibular stimulations induced by multi-directional movement of an individual within the rotating environment to provide the required design criteria. This was done to determine the overall impact of artificial gravity simulations on potential design configurations and crew operational procedures. Gross locomotion and fine motor performance were evaluated. Results indicate that crew orientation, rotational rates, vehicle design configurations, and operational procedures may be used to reduce the severity of the adverse effects of the Coriolis and cross-coupled angular accelerations acting on masses moving within a rotating environment. Results further indicate that crew selection, motivation, and short-term exposures to the rotating environment may be important considerations for future crew indoctrination and training programs

Two-component {CH} system: Inverse Scattering, Peakons and Geometry

An inverse scattering transform method corresponding to a Riemann-Hilbert problem is formulated for CH2, the two-component generalization of the Camassa-Holm (CH) equation. As an illustration of the method, the multi - soliton solutions corresponding to the reflectionless potentials are constructed in terms of the scattering data for CH2.Comment: 22 pages, 3 figures, draft, please send comment

Variational Principles for Lagrangian Averaged Fluid Dynamics

The Lagrangian average (LA) of the ideal fluid equations preserves their transport structure. This transport structure is responsible for the Kelvin circulation theorem of the LA flow and, hence, for its convection of potential vorticity and its conservation of helicity. Lagrangian averaging also preserves the Euler-Poincar\'e (EP) variational framework that implies the LA fluid equations. This is expressed in the Lagrangian-averaged Euler-Poincar\'e (LAEP) theorem proven here and illustrated for the Lagrangian average Euler (LAE) equations.Comment: 23 pages, 3 figure

Phase diagram for a single flexible Stockmayer polymer at zero field

The equilibrium conformations of a flexible permanent magnetic chain that consists of a sequence of linked magnetic colloidal nanoparticles with short-ranged Lennard-Jones attractive interactions (Stockmayer polymer) are thoroughly analysed via Langevin dynamics simulations. A tentative phase diagram is presented for a chain of length $N=100$. The phase diagram exhibits several unusual conformational phases when compared with the non-magnetic chains. These phases are characterised by a large degree of conformational anisotropy, and consist of closed chains, helicoidal-like states, partially collapsed states, and very compact disordered states. The phase diagram contains several interesting features like the existence of at least two 'triple points'

The Hamiltonian structure and Euler-Poincar\'{e} formulation of the Vlasov-Maxwell and gyrokinetic systems

We present a new variational principle for the gyrokinetic system, similar to the Maxwell-Vlasov action presented in Ref. 1. The variational principle is in the Eulerian frame and based on constrained variations of the phase space fluid velocity and particle distribution function. Using a Legendre transform, we explicitly derive the field theoretic Hamiltonian structure of the system. This is carried out with a modified Dirac theory of constraints, which is used to construct meaningful brackets from those obtained directly from Euler-Poincar\'{e} theory. Possible applications of these formulations include continuum geometric integration techniques, large-eddy simulation models and Casimir type stability methods. [1] H. Cendra et. al., Journal of Mathematical Physics 39, 3138 (1998)Comment: 36 pages, 1 figur

An Optimal Control Formulation for Inviscid Incompressible Ideal Fluid Flow

In this paper we consider the Hamiltonian formulation of the equations of incompressible ideal fluid flow from the point of view of optimal control theory. The equations are compared to the finite symmetric rigid body equations analyzed earlier by the authors. We discuss various aspects of the Hamiltonian structure of the Euler equations and show in particular that the optimal control approach leads to a standard formulation of the Euler equations -- the so-called impulse equations in their Lagrangian form. We discuss various other aspects of the Euler equations from a pedagogical point of view. We show that the Hamiltonian in the maximum principle is given by the pairing of the Eulerian impulse density with the velocity. We provide a comparative discussion of the flow equations in their Eulerian and Lagrangian form and describe how these forms occur naturally in the context of optimal control. We demonstrate that the extremal equations corresponding to the optimal control problem for the flow have a natural canonical symplectic structure.Comment: 6 pages, no figures. To appear in Proceedings of the 39th IEEEE Conference on Decision and Contro

Lattice Models of Quantum Gravity

Standard Regge Calculus provides an interesting method to explore quantum gravity in a non-perturbative fashion but turns out to be a CPU-time demanding enterprise. One therefore seeks for suitable approximations which retain most of its universal features. The $Z_2$-Regge model could be such a desired simplification. Here the quadratic edge lengths $q$ of the simplicial complexes are restricted to only two possible values $q=1+\epsilon\sigma$, with $\sigma=\pm 1$, in close analogy to the ancestor of all lattice theories, the Ising model. To test whether this simpler model still contains the essential qualities of the standard Regge Calculus, we study both models in two dimensions and determine several observables on the same lattice size. In order to compare expectation values, e.g. of the average curvature or the Liouville field susceptibility, we employ in both models the same functional integration measure. The phase structure is under current investigation using mean field theory and numerical simulation.Comment: 4 pages, 1 figure