Search

We present a numerical investigation of stochastic transport in ideal fluids. According to Holm (Proc Roy Soc, 2015) and Cotter et al. (2017), the principles of transformation theory and multi-time homogenisation, respectively, imply a physically meaningful, data-driven approach for decomposing the fluid transport velocity into its drift and stochastic parts, for a certain class of fluid flows. In the current paper, we develop new methodology to implement this velocity decomposition and then numerically integrate the resulting stochastic partial differential equation using a finite element discretisation for incompressible 2D Euler fluid flows. The new methodology tested here is found to be suitable for coarse graining in this case. Specifically, we perform uncertainty quantification tests of the velocity decomposition of Cotter et al. (2017), by comparing ensembles of coarse-grid realisations of solutions of the resulting stochastic partial differential equation with the "true solutions" of the deterministic fluid partial differential equation, computed on a refined grid. The time discretization used for approximating the solution of the stochastic partial differential equation is shown to be consistent. We include comprehensive numerical tests that confirm the non-Gaussianity of the stream function, velocity and vorticity fields in the case of incompressible 2D Euler fluid flows., 41 pages, 26 figures Minor changes -- updated figures to improve readability. Corrected typos. Shifted Remark 7 to just after Assumption A1. Added Remark 8

Starting with the exact factorization of the molecular wavefunction, this paper presents the results from the numerical implementation in nonadiabatic molecular dynamics of the recently proposed bohmion method. Within the context of quantum hydrodynamics, we introduce a regularized nuclear Bohm potential admitting solutions comprising a train of δ-functions which provide a finite-dimensional sampling of the hydrodynamic flow paths. The bohmion method inherits all the basic conservation laws from its underlying variational structure and captures electronic decoherence. After reviewing the general theory, the method is applied to the well-known Tully models, which are used here as benchmark problems. In the present case of study, we show that the new method accurately reproduces both electronic decoherence and nuclear population dynamics.

Suppose the observations of Lagrangian trajectories for fluid flow in some physical situation can be modelled sufficiently accurately by a spatially correlated It\^o stochastic process (with zero mean) obtained from data which is taken in fixed Eulerian space. Suppose we also want to apply Hamilton's principle to derive the stochastic fluid equations for this situation. Now, the variational calculus for applying Hamilton's principle requires the Stratonovich process, so we must transform from It\^o noise in the \emph{data frame} to the equivalent Stratonovich noise. However, the transformation from the It\^o process in the data frame to the corresponding Stratonovich process shifts the drift velocity of the transformed Lagrangian fluid trajectory out of the data frame into a non-inertial frame obtained from the It\^o correction. The issue is, "Will non-inertial forces arising from this transformation of reference frames make a difference in the interpretation of the solution behaviour of the resulting stochastic equations?" This issue will be resolved by elementary considerations., Comment: 12 pages, 1 figure, 2nd version, comments welcome!

A famous feature of the Camassa–Holm equation is its admission of peaked soliton solutions known as peakons. We investigate this equation under the influence of stochastic transport. Noting that peakons are weak solutions of the equation, we present a finite-element discretization for it, which we use to explore the formation of peakons. Our simulations using this discretization reveal that peakons can still form in the presence of stochastic perturbations. Peakons can emerge both through wave breaking, as the slope turns vertical, and without wave breaking as the inflection points of the velocity profile rise to reach the summit.

We study the three-dimensional, incompressible, non-hydrostatic Boussinesq fluid equations, which are applicable to the dynamics of the oceans and atmosphere. These equations describe the interplay between velocity and buoyancy in a rotating frame. A hierarchy of dynamical variables is introduced whose members $\Omega_{m}(t)$ ($1 \leq m < \infty$) are made up from the respective sum of the $L^{2m}$-norms of vorticity and the density gradient. Each $\Omega_{m}(t)$ has a lower bound in terms of the inverse Rossby number, $Ro^{-1}$, that turns out to be crucial to the argument. For convenience, the $\Omega_{m}$ are also scaled into a new set of variables $D_{m}(t)$. By assuming the existence and uniqueness of solutions, conditional upper bounds are found on the $D_{m}(t)$ in terms of $Ro^{-1}$ and the Reynolds number $Re$. These upper bounds vary across bands in the $\{D_{1},\,D_{m}\}$ phase plane. The boundaries of these bands depend subtly upon $Ro^{-1}$, $Re$, and the inverse Froude number $Fr^{-1}$. For example, solutions in the lower band conditionally live in an absorbing ball in which the maximum value of $\Omega_{1}$ deviates from $Re^{3/4}$ as a function of $Ro^{-1},\,Re$ and $Fr^{-1}$., Comment: 24 pages, 3 figures and 1 table

This paper formulates a variational approach for treating observational uncertainty and/or computational model errors as stochastic transport in dynamical systems governed by action principles under nonholonomic constraints. For this purpose, we derive, analyze and numerically study the example of an unbalanced spherical ball rolling under gravity along a stochastic path. Our approach uses the Hamilton-Pontryagin variational principle, constrained by a stochastic rolling condition, which we show is equivalent to the corresponding stochastic Lagrange-d'Alembert principle. In the example of the rolling ball, the stochasticity represents uncertainty in the observation and/or error in the computational simulation of the angular velocity of rolling. The influence of the stochasticity on the deterministically conserved quantities is investigated both analytically and numerically. Our approach applies to a wide variety of stochastic, nonholonomically constrained systems, because it preserves the mathematical properties inherited from the variational principle. Keywords: Nonholonomic constraints, Stochastic dynamics, Transport noise., Comment: Published version

We study the $G$-strand equations that are extensions of the classical chiral model of particle physics in the particular setting of broken symmetries described by symmetric spaces. These equations are simple field theory models whose configuration space is a Lie group, or in this case a symmetric space. In this class of systems, we derive several models that are completely integrable on finite dimensional Lie group $G$ and we treat in more details examples with symmetric space $SU(2)/S^1$ and $SO(4)/SO(3)$. The later model simplifies to an apparently new integrable $9$ dimensional system. We also study the $G$-strands on the infinite dimensional group of diffeomorphisms, which gives, together with the Sobolev norm, systems of $1+2$ Camassa-Holm equations. The solutions of these equations on the complementary space related to the Witt algebra decomposition are the odd function solutions., To appear in Proceedings of the Royal Society A

Initially, this paper reviews the mathematical issues surrounding hydrostatic primitive equations (HPEs) and non-hydrostatic primitive equations (NPEs) that have been used extensively in numerical weather prediction and climate modelling. A new impetus has been provided by a recent proof of the existence and uniqueness of solutions of viscous HPEs on a cylinder with Neumann-like boundary conditions on the top and bottom. In contrast, the regularity of solutions of NPEs remains an open question. With this HPE regularity result in mind, the second issue examined in this paper is whether extreme events are allowed to arise spontaneously in their solutions. Such events could include, for example, the sudden appearance and disappearance of locally intense fronts that do not involve deep convection. Analytical methods are used to show that for viscous HPEs, the creation of small-scale structures is allowed locally in space and time at sizes that scale inversely with the Reynolds number.

We introduce a set of coupled equations for multilayer water waves that removes the ill-posedness of the multilayer Green-Naghdi (MGN) equations in the presence of shear. The new well-posed equations are Hamiltonian and in the absence of imposed background shear they retain the same travelling wave solutions as MGN. We call the new model the Square Root Depth equations, from the modified form of their kinetic energy of vertical motion. Our numerical results show how the Square Root Depth equations model the effects of multilayer wave propagation and interaction, with and without shear., Comment: 10 pages, 5 figures

The EPDiff equation (or dispersionless Camassa-Holm equation in 1D) is a well known example of geodesic motion on the Diff group of smooth invertible maps (diffeomorphisms). Its recent two-component extension governs geodesic motion on the semidirect product ${\rm Diff}\circledS{\cal F}$, where $\mathcal{F}$ denotes the space of scalar functions. This paper generalizes the second construction to consider geodesic motion on ${\rm Diff} \circledS\mathfrak{g}$, where $\mathfrak{g}$ denotes the space of scalar functions that take values on a certain Lie algebra (for example, $\mathfrak{g}=\mathcal{F}\otimes\mathfrak{so}(3)$). Measure-valued delta-like solutions are shown to be momentum maps possessing a dual pair structure, thereby extending previous results for the EPDiff equation. The collective Hamiltonians are shown to fit into the Kaluza-Klein theory of particles in a Yang-Mills field and these formulations are shown to apply also at the continuum PDE level. In the continuum description, the Kaluza-Klein approach produces the Kelvin circulation theorem., 22 pages, 2 figures. Submitted to Proc. R. Soc. A

The dynamics of Vlasov kinetic moments is shown to be Lie-Poisson on the dual Lie algebra of symmetric contravariant tensor fields. The corresponding Lie bracket is identified with the symmetric Schouten bracket and the moment Lie algebra is related with a bundle of bosonic Fock spaces, where creation and annihilation operators are used to construct the cold plasma closure. Kinetic moments are also shown to define a momentum map, which is infinitesimally equivariant. This momentum map is the dual of a Lie algebra homomorphism, defined through the Schouten bracket. Finally the moment Lie-Poisson bracket is extended to anisotropic interactions., Comment: 19 pages, no figures. Submitted to Phys. Lett. A

Commutator errors arise in large-eddy simulation of incompressible turbulent flow from the application of non-uniform filters to the continuity -- and Navier-Stokes equations. For non-uniform, high order filters with bounded moments the magnitude of the commutator errors is shown to be of the same order as that of the turbulent stress fluxes. Consequently, one cannot reduce the size of the commutator errors independently of the turbulent stress terms by any judicious construction of such filter operators. Independent control over the commutator errors compared to the turbulent stress fluxes can, instead, be obtained by appropriately restricting the spatial variations of the filter-width and filter-skewness. For situations in which the dynamical consequences of the commutator errors are significant, e.g., near solid boundaries, explicit similarity modelling for the commutator errors is proposed, including the application of Leray regularization. The performance of this commutator error parametrization is illustrated for the one-dimensional Burgers equation. The Leray approach is found to capture the filtered flow with higher accuracy than conventional similarity modelling, which is particularly relevant for large filter-width variations.

In this paper we introduce and study a new model for three–dimensional turbulence, the Leray– α model. This model is inspired by the Lagrangian averaged Navier–Stokes– α model of turbulence (also known Navier–Stokes– α model or the viscous Camassa–Holm equations). As in the case of the Lagrangian averaged Navier–Stokes– α model, the Leray– α model compares successfully with empirical data from turbulent channel and pipe flows, for a wide range of Reynolds numbers. We establish here an upper bound for the dimension of the global attractor (the number of degrees of freedom) of the Leray– α model of the order of ( L / l d ) 12/7 , where L is the size of the domain and l d is the dissipation length–scale. This upper bound is much smaller than what one would expect for three–dimensional models, i.e. ( L / l d ) 3 . This remarkable result suggests that the Leray– α model has a great potential to become a good sub–grid–scale large–eddy simulation model of turbulence. We support this observation by studying, analytically and computationally, the energy spectrum and show that in addition to the usual k −5/3 Kolmogorov power law the inertial range has a steeper power–law spectrum for wavenumbers larger than 1/ α . Finally, we propose a Prandtl–like boundary–layer model, induced by the Leray– α model, and show a very good agreement of this model with empirical data for turbulent boundary layers.

This paper derives stochastic partial differential equations (SPDEs) for fluid dynamics from a stochastic variational principle (SVP). The Legendre transform of the Lagrangian formulation of these SPDEs yields their Lie-Poisson Hamiltonian form. The paper proceeds by: taking variations in the SVP to derive stochastic Stratonovich fluid equations; writing their It\^o representation; and then investigating the properties of these stochastic fluid models in comparison with each other, and with the corresponding deterministic fluid models. The circulation properties of the stochastic Stratonovich fluid equations are found to closely mimic those of the deterministic ideal fluid models. As with deterministic ideal flows, motion along the stochastic Stratonovich paths also preserves the helicity of the vortex field lines in incompressible stochastic flows. However, these Stratonovich properties are not apparent in the equivalent It\^o representation, because they are disguised by the quadratic covariation drift term arising in the Stratonovich to It\^o transformation. This term is a geometric generalisation of the quadratic covariation drift term already found for scalar densities in Stratonovich's famous 1966 paper. The paper also derives motion equations for two examples of stochastic geophysical fluid dynamics (SGFD); namely, the Euler-Boussinesq and quasigeostropic approximations., Comment: 19 pages, no figures, 2nd version. To appear in Proc Roy Soc A. Comments to author are still welcome!

Truncated Taylor expansions of smooth flow maps are used in Hamilton's principle to derive a multiscale Lagrangian particle representation of ideal fluid dynamics. Numerical simulations for scattering of solutions at one level of truncation are found to produce solutions at higher levels. These scattering events to higher levels in the Taylor expansion are interpreted as modeling a cascade to smaller scales., 1st revision -- Comments welcome

The present paper derives systems of partial differential equations that admit a quadratic zero curvature representation for an arbitrary real semisimple Lie algebra. It also determines the general form of Hamilton's principles and Hamiltonians for these systems and analyzes the linear stability of their equilibrium solutions in the examples of $\mathfrak{so}(3)$ and $\mathfrak{sl}(2,\mathbb{R})$., 17 pages, no figures. First version, comments welcome!

International audience; Modifications of the equations of ideal fluid dynamics with advected quantities are introduced that allow selective decay of either the energy h or the Casimir quantities C in the Lie-Poisson (LP) formulation. The dissipated quantity (energy or Casimir, respectively) is shown to decrease in time until the modified system reaches an equilibrium state consistent with ideal energy-Casimir equilibria, namely d(h + C) = 0. The result holds for LP equations in general, independently of the Lie algebra and the choice of Casimir. This selective decay process is illustrated with a number of examples in 2D and 3D magnetohydrodynamics. © 2014 IOP Publishing Ltd & London Mathematical Society.

We propose the viscous Camassa-Holm equations as a closure approximation for the Reynolds-averaged equations of the incompressible Navier-Stokes fluid. This approximation is tested on turbulent channel and pipe flows with steady mean. Analytical solutions for the mean velocity and the Reynolds shear stress are consistent with experiments in most of the flow region.

We propose a new class of models for the mean motion of ideal incompressible fluids in three dimensions, including stratification and rotation. In these models, the amplitude of the rapid fluctuations introduces a length scale, α, below which wave activity is filtered by both linear and nonlinear dispersion. This filtering enhances the stability and regularity of the new fluid models without compromising either their large scale behavior, or their conservation laws. These models also describe geodesic motion on the volume-preserving diffeomorphism group for a metric containing the H1 norm of the fluid velocity.

We present two large families of Silnikov-type homoclinic orbits in a two-mode model that describes second-harmonic generation in a passive optical cavity. These families of homoclinic orbits give rise to chaotic dynamics in the model.

We examine the stability of flows which are the sum of a linear flow with elliptic streamlines and a transverse standing wave. Such flows are exact solutions of the 3D Euler equations for an inviscid incompressible fluid. We examine the stability of such flows under local 3D high frequency perturbations, which can be viewed as secondary perturbations of the elliptic flow itself. We find that these flows are unstable to such perturbations and compute growth rate in several cases. {copyright} {ital 1997} {ital The American Physical Society}

We discuss three examples in which one may extend integrable Euler--Poincar\'e ODEs to integrable Euler--Poincar\'e PDEs in the matrix G-Strand context. After describing matrix G-Strand examples for $SO(3)$ and $SO(4)$ we turn our attention to $SE(3)$ where the matrix G-Strand equations recover the exact rod theory in the convective representation. We then find a zero curvature representation (ZCR) of these equations and establish the conditions under which they are completely integrable. Thus, the G-Strand equations turn out to be a rich source of integrable systems. The treatment is meant to be expository and most concepts are explained in examples in the language of vectors in $\mathbb{R}^3$., Comment: COMMENTS WELCOME!

We introduce a measurement (the discrepancy ) of the minimum energy needed to transform from a standard parametrized planar curve c 0 to an observed curve c 1 . To this end, we say that a curve of transformations in the special Euclidean group SE(2) is admissible if it maps the source curve to the target curve under the point-wise action of SE(2) on the plane. After endowing the group SE(2) with a left-invariant metric, we define a relative geodesic in SE(2) to be a critical point of the energy functional associated to the metric, over all admissible curves. The discrepancy is then defined as the value of the energy of the minimizing relative geodesic. In the first part of the paper, we derive a scalar ODE which is a necessary condition for a curve in SE(2) to be a relative geodesic, and we discuss some of the properties of the discrepancy. In the second part of the paper, we consider discrete curves, and by means of a variational principle, we derive a system of discrete equations for the relative geodesics. We finish with several examples.

We study a trajectory-planning problem whose solution path evolves by means of a Lie group action and passes near a designated set of target positions at particular times. This is a higher-order variational problem in optimal control, motivated by potential applications in computational anatomy and quantum control. Reduction by symmetry in such problems naturally summons methods from Lie group theory and Riemannian geometry. A geometrically illuminating form of the Euler-Lagrange equations is obtained from a higher-order Hamilton-Pontryagin variational formulation. In this context, the previously known node equations are recovered with a new interpretation as Legendre-Ostrogradsky momenta possessing certain conservation properties. Three example applications are discussed as well as a numerical integration scheme that follows naturally from the Hamilton-Pontryagin principle and preserves the geometric properties of the continuous-time solution., Comment: 32 pages, 3 figures, comments welcome!

A variational framework is defined for vertical slice models with three-dimensional velocity depending only on x and z . The models that result from this framework are Hamiltonian, and have a Kelvin–Noether circulation theorem that results in a conserved potential vorticity in the slice geometry. These results are demonstrated for the incompressible Euler–Boussinesq equations with a constant temperature gradient in the y -direction (the Eady–Boussinesq model), which is an idealized problem used to study the formation and subsequent evolution of weather fronts. We then introduce a new compressible extension of this model. Unlike the incompressible model, the compressible model does not produce solutions that are also solutions of the three-dimensional equations, but it does reduce to the Eady–Boussinesq model in the low Mach number limit. Hence, the new model could be used in asymptotic limit error testing for compressible weather models running in a vertical slice configuration.

International audience; The problem of parameterizing the interactions of larger scales and smaller scales in fluid flows is addressed by considering a property of two-dimensional (2D) incompressible turbulence. The property we consider is selective decay, in which a Casimir of the ideal formulation (enstrophy in 2D flows, helicity in three-dimensional flows) decays in time, while the energy stays essentially constant. This paper introduces a mechanism that produces selective decay by enforcing Casimir dissipation in fluid dynamics. This mechanism turns out to be related in certain cases to the numerical method of anticipated vorticity discussed in Sadourny and Basdevant (1981 C. R. Acad. Sci. Paris 292 1061-4, 1985 J. Atm. Sci. 42 1353-63). Several examples are given and a general theory of selective decay is developed that uses the Lie-Poisson structure of the ideal theory. A scale-selection operator allows the resulting modifications of the fluid motion equations to be interpreted in several examples as parametrizing the nonlinear, dynamical interactions between disparate scales. The type of modified fluid equation systems derived here may be useful in modelling turbulent geophysical flows where it is computationally prohibitive to rely on the slower, indirect effects of a realistic viscosity, such as in large-scale, coherent, oceanic flows interacting with much smaller eddies. © 2013 IOP Publishing Ltd & London Mathematical Society.

The Kirchhoff-Kida family of elliptical vortex columns in flows with uniform strain and rotation displays a rich variety of dynamical behaviours, even in a purely two-dimensional setting. In this paper, we address the stability of these columns with respect to three-dimensional perturbations via the geometrical optics method. In the case when the external strain is equal to zero, the analysis reduces to the stability of a steady elliptical vortex in a rotating frame. When the external strain is non-zero, the stability analysis reduces to the theory of a Schrodinger equation with quasi-periodic potential. We present stability results for a variety of different Kirchhoff-Kida flows. The vortex columns are typically unstable except when the interior vorticity is approximately the negative of the background vorticity, so that the flow in the inertial frame is nearly a potential flow.

The dynamics of an ensemble of two-level atoms in a single-mode resonant laser cavity with external pumping and a weak coherent probe are investigated. Homoclinic orbits connected to the completely inverted atomic state are found analytically using linear perturbation theory (the Melnikov method) in the near-integrable limit, and are continued in the parameter space to O(1) parameter values using the bifurcation code AUTO. The breakup of these homoclinic orbits is believed to be a source of chaos in the system.

We develop a variational method of deriving stochastic partial differential equations whose solutions follow the flow of a stochastic vector field. As an example in one spatial dimension we numerically simulate singular solutions (peakons) of the stochastically perturbed Camassa-Holm (CH) equation derived using this method. These numerical simulations show that peakon soliton solutions of the stochastically perturbed CH equation persist and provide an interesting laboratory for investigating the sensitivity and accuracy of adding stochasticity to finite dimensional solutions of stochastic partial differential equations (SPDE). In particular, some choices of stochastic perturbations of the peakon dynamics by Wiener noise (canonical Hamiltonian stochastic deformations, or CH-SD) allow peakons to interpenetrate and exchange order on the real line in overtaking collisions, although this behaviour does not occur for other choices of stochastic perturbations which preserve the Euler-Poincar\'e structure of the CH equation (parametric stochastic deformations, or P-SD), and it also does not occur for peakon solutions of the unperturbed deterministic CH equation. The discussion raises issues about the science of stochastic deformations of finite-dimensional approximations of evolutionary PDE and the sensitivity of the resulting solutions to the choices made in stochastic modelling., Comment: 21 pages, 15 figures -- 2nd version

A quantum spline is a smooth curve parameterised by time in the space of unitary transformations, whose associated orbit on the space of pure states traverses a designated set of quantum states at designated times, such that the trace norm of the time rate of change of the associated Hamiltonian is minimised. The solution to the quantum spline problem is obtained, and is applied in an example that illustrates quantum control of coherent states. An efficient numerical scheme for computing quantum splines is discussed and implemented in the examples., Comment: Published version: 5 pages, 3 figures

We study the dynamics of the classical homogeneous SU(2)⊗U(1) gauge theory in the spontaneously broken phase. A transition from regular to chaotic dynamics is shown to occur for energy densities e>0.3 m W 2 / G F ≈2×10 10 GeV/fm 3 . The regular regime is characterized by a sharp frequency spectrum for the Yang-Mills fields. The frequency spectrum in the chaotic regime is continuous and the sharp frequency lines of the regular regime are dispersed. In the early universe or at high energy densities in a Yang-Mills plasma this transition may have observable effects that can be studied by analyzing the frequency spectrum and the time correlation functions. In the chaotic regime these time correlation functions decay considerably faster than in the regular regime.

We discuss the problem of designing a quantum computer based on one-dimensional “alternating” Ising systems (linear chains with periodically recurring spin groups) in an external magnetic field which exceeds the interaction between spins. Loading and processing of information in the alternating Ising system may be accomplished by using a scheme suggested recently by Lloyd for heteropolymer systems. The detailed operation of a simple quantum logical device is described in the framework of a binary Ising system. Estimates of physical parameters are presented that show that the experimental realization of such quantum computer elements would be feasible as a research task. However, many difficulties remain to be addressed, before the approach discussed here would be applicable in real devices.

Quantum effects are considered in the dynamics of a system of N paramagnetic atoms in a resonant cavity interacting with a constant magnetic field and with a resonant external magnetic field. In the semi-classical limit (classical radiation field in the cavity) this resonantly driven system shows developed (global) chaos. Expectation-value dynamics shows however that quantum corrections cause a departure from the semi-classical chaotic dynamics on a time-scale τ ℏ ∼ln N and that quantum correlation functions grow exponentially in time. The possibility of experimentally observing this effect is discussed.

We derive a new completely integrable dispersive shallow water equation that is biHamiltonian and thus possesses an infinite number of conservation laws in involution. The equation is obtained by using an asymptotic expansion directly in the Hamiltonian for Euler's equations in the shallow water regime. The soliton solution for this equation has a limiting form that has a discontinuity in the first derivative at its peak., Comment: LaTeX file. Figure available from authors upon request

The relations between smooth and peaked soliton solutions are reviewed for the Camassa-Holm (CH) shallow water wave equation in one spatial dimension. The canonical Hamiltonian formulation of the CH equation in action-angle variables is expressed for solitons by using the scattering data for its associated isospectral eigenvalue problem, rephrased as a Riemann-Hilbert problem. The momentum map from the action-angle scattering variables $T^*({\mathbb{T}^N})$ to the flow momentum ($\mathfrak{X}^*$) provides the Eulerian representation of the $N$-soliton solution of CH in terms of the scattering data and squared eigenfunctions of its isospectral eigenvalue problem. The dispersionless limit of the CH equation and its resulting peakon solutions are examined by using an asymptotic expansion in the dispersion parameter. The peakon solutions of the dispersionless CH equation in one dimension are shown to generalize in higher dimensions to peakon wave-front solutions of the EPDiff equation whose associated momentum is supported on smoothly embedded subspaces. The Eulerian representations of the singular solutions of both CH and EPDiff are given by the (cotangent-lift) momentum maps arising from the left action of the diffeomorphisms on smoothly embedded subspaces., First version -- comments welcome! Submitted to JPhysA

Competition between two degenerate spatial modes in a ring-cavity laser is shown to lead to chaotic time-evolution of their amplitudes. Analysis of this temporal chaos implies that the laser output will show randomly intermittent bursts of apparent spatial complexity due to interference between the two competing modes.

An optimization method used in image-processing (metamorphosis) is found to imply Euler's equations for incompressible flow of an inviscid fluid, without requiring that the Lagrangian particle labels exactly follow the flow lines of the Eulerian velocity vector field. Thus, an optimal control problem and an optimization problem for incompressible ideal fluid flow both yield the \emph {same} Euler fluid equations, although their Lagrangian parcel dynamics are \emph{different}. This is a result of the \emph{gauge freedom} in the definition of the fluid pressure for an incompressible flow, in combination with the symmetry of fluid dynamics under relabeling of their Lagrangian coordinates. Similar ideas are also illustrated for SO(N) rigid body motion., Comment: 12 pages

We present a tri-Hamiltonian formulation of the self-induced transparency (SIT) equations.

Hamiltonian Silnikov orbits are shown to be produced when a system of two-level atoms interacting via a self-consistent electric field in a single-model lossless cavity is perturbed by a small-amplitude probe laser (a weak, externally imposed, monochromatic electromagnetic field). This result is obtained by utilizing the two constants of motion which exist in the perturbed problem and by calculating a Melnikov function whose nondegenerate zeros imply transverse intersections of the stable and unstable manifolds of a hyperbolic fixed point set. The dynamics of the system on the Silnikov orbits contains a Smale horsehoe construction, leading to intermittent switching of the sign of the electric field which has extreme sensitivity to initial conditions.

The Green Nagdhi equations are frequently used as a model of the wave-like behaviour of the free surface of a fluid, or the interface between two homogeneous fluids of differing densities. Here we show that their multilayer extension arises naturally from a framework based on the Euler Poincare theory under an ansatz of columnar motion. The framework also extends to the travelling wave solutions of the equations. We present numerical solutions of the travelling wave problem in a number of flow regimes. We find that the free surface and multilayer waves can exhibit intriguing differences compared to the results of single layer or rigid lid models.

A framework for statistical-mechanical analysis of quantum Hamiltonians is introduced. The approach is based upon a gradient flow equation in the space of Hamiltonians such that the eigenvectors of the initial Hamiltonian evolve toward those of the reference Hamiltonian. The nonlinear double-bracket equation governing the flow is such that the eigenvalues of the initial Hamiltonian remain unperturbed. The space of Hamiltonians is foliated by compact invariant subspaces, which permits the construction of statistical distributions over the Hamiltonians. In two dimensions, an explicit dynamical model is introduced, wherein the density function on the space of Hamiltonians approaches an equilibrium state characterised by the canonical ensemble. This is used to compute quenched and annealed averages of quantum observables., Comment: 8 pages, 2 figures, references added

This paper applies Hamiltonian methods to the Stokes representation of the one-beam and two-beam problems of polarized optical pulses propagating as travelling waves in nonlinear media. We treat these two dynamical systems as follows. First, we use the reduction method of Marsden and Weinstein to map each of the systems to the two-dimensional sphere, S 2 . The resulting reduced systems are then analyzed from the viewpoints of their stability properties and of bifurcations with symmetry; in particular, several degenerate bifurcations are found and described. We also establish the presence of chaotic dynamics in these systems by demonstrating the existence of Smale horseshoe maps in the three- and four-dimensional cases, as well as Arnold diffusion in the higher-dimensional cases. The method we use to establish such complex dynamics is the Mel'nikov technique, as extended by Holmes and Marsden, and Wiggins for the higher-dimensional cases. These results apply to perturbations of homoclinic and heteroclinic orbits of the reduced integrable problems for static, as well as travelling-wave, solutions describing either a single opt ical beam, or two such beams counterpropagating. Thus, we show that these optics problems exhibit complex dynamics and predict the experimental consequences of this dynamics.

An asymptotic Green's function in homogeneous anisotropic viscoelastic media is derived. The Green's function in viscoelastic media is formally similar to that in elastic media, but its computation is more involved. The stationary slowness vector is, in general, complex valued and inhomogeneous. Its computation involves finding two independent real-valued unit vectors which specify the directions of its real and imaginary parts and can be done either by iterations or by solving a system of coupled polynomial equations. When the stationary slowness direction is found, all quantities standing in the Green's function such as the slowness vector, polarization vector, phase and energy velocities and principal curvatures of the slowness surface can readily be calculated. The formulae for the exact and asymptotic Green's functions are numerically checked against closed-form solutions for isotropic and simple anisotropic, elastic and viscoelastic models. The calculations confirm that the formulae and developed numerical codes are correct. The computation of the P-wave Green's function in two realistic materials with a rather strong anisotropy and absorption indicates that the asymptotic Green's function is accurate at distances greater than several wavelengths from the source. The error in the modulus reaches at most 4% at distances greater than 15 wavelengths from the source.

Many dynamical systems, such as the Lotka-Volterra predator-prey model and the Euler equations for the free rotation of a rigid body, are PT symmetric. The standard and well-known real solutions to such dynamical systems constitute an infinitessimal subclass of the full set of complex solutions. This paper examines a subset of the complex solutions that contains the real solutions, namely, those having PT symmetry. The condition of PT symmetry selects out complex solutions that are periodic., 13 pages, 12 figures, to appear in Fast Track Communications, Journal of Physics A: Mathematical and Theoretical

The Vlasov equation for the collisionless evolution of the single-particle probability distribution function (PDF) is a well-known Lie-Poisson Hamiltonian system. Remarkably, the operation of taking the moments of the Vlasov PDF preserves the Lie-Poisson structure. The individual particle motions correspond to singular solutions of the Vlasov equation. The paper focuses on singular solutions of the problem of geodesic motion of the Vlasov moments. These singular solutions recover geodesic motion of the individual particles., 16 pages, no figures. Submitted to Phys. Lett. A

The motion of a classical pendulum in a gravitational field of strength g is explored. The complex trajectories as well as the real ones are determined. If g is taken to be imaginary, the Hamiltonian that describes the pendulum becomes PT-symmetric. The classical motion for this PT-symmetric Hamiltonian is examined in detail. The complex motion of this pendulum in the presence of an external periodic forcing term is also studied., 11 pages, 12 figures

We derive equations of motion for the dynamics of anisotropic particles directly from the dissipative Vlasov kinetic equations, with the dissipation given by the double bracket approach (Double Bracket Vlasov, or DBV). The moments of the DBV equation lead to a nonlocal form of Darcy's law for the mass density. Next, kinetic equations for particles with anisotropic interaction are considered and also cast into the DBV form. The moment dynamics for these double bracket kinetic equations is expressed as Lie-Darcy continuum equations for densities of mass and orientation. We also show how to obtain a Smoluchowski model from a cold plasma-like moment closure of DBV. Thus, the double bracket kinetic framework serves as a unifying method for deriving different types of dynamics, from density--orientation to Smoluchowski equations. Extensions for more general physical systems are also discussed., Comment: 19 pages; no figures. Submitted to Proc. Roy. Soc. A

Vorticity dynamics of the three-dimensional incompressible Euler equations is cast into a quaternionic representation governed by the Lagrangian evolution of the tetrad consisting of the growth rate and rotation rate of the vorticity. In turn, the Lagrangian evolution of this tetrad is governed by another that depends on the pressure Hessian. Together these form the basis for a direction of vorticity theorem. Moreover, in this representation, fluid particles carry ortho-normal frames whose Lagrangian evolution in time are shown to be directly related to the Frenet-Serret equations for a vortex line. The frame dynamics suggest an elegant Lagrangian relation regarding the pressure Hessian tetrad. The equations for ideal MHD are similarly considered., Comment: 14 pages, 3 figures, revised abstract, To appear in Nonlinearity

After a review of the isentropic compressible magnetohydrodynamics (ICMHD) equations, a quaternionic framework for studying the alignment dynamics of a general fluid flow is explained and applied to the ICMHD equations., Comment: 12 pages, 2 figures, submitted to a Focus Issue of New Journal of Physics on "Magnetohydrodynamics and the Dynamo Problem" J-F Pinton, A Pouquet, E Dormy and S Cowley, editors

Experimentalists now measure intense rotations of Lagrangian particles in turbulent flows by tracking their trajectories and Lagrangian-average velocity gradients at high Reynolds numbers. This paper formulates the dynamics of an orthonormal frame attached to each Lagrangian fluid particle undergoing three-axis rotations, by using quaternions in combination with Ertel's theorem for frozen-in vorticity. The method is applicable to a wide range of Lagrangian flows including the three-dimensional Euler equations and its variants such as ideal MHD. The applicability of the quaterionic frame description to Lagrangian averaged velocity gradient dynamics is also demonstrated., Comment: 9 pages, one figure, revised