Search

This work involves theoretical and numerical analysis of the thermal quasi-geostrophic (TQG) model of submesoscale geophysical fluid dynamics (GFD). Physically, the TQG model involves thermal geostrophic balance, in which the Rossby number, the Froude number and the stratification parameter are all of the same asymptotic order. The main analytical contribution of this paper is to construct local-in-time unique strong solutions for the TQG model. For this, we show that solutions of its regularised version α -TQG converge to solutions of TQG as its smoothing parameter α → 0 and we obtain blow-up criteria for the α -TQG model. The main contribution of the computational analysis is to verify the rate of convergence of α -TQG solutions to TQG solutions as α → 0 , for example, simulations in appropriate GFD regimes., Competing Interests: Conflict of interestThe authors have no competing interests to declare that are relevant to the content of this work. All authors certify that they have no affiliations with or involvement in any organisation or entity with any financial interest or non-financial interest in the subject matter or materials discussed in this manuscript., (© The Author(s) 2023.)

We are modelling multiscale, multi-physics uncertainty in wave-current interaction (WCI). To model uncertainty in WCI, we introduce stochasticity into the wave dynamics of two classic models of WCI, namely the generalised Lagrangian mean (GLM) model and the Craik-Leibovich (CL) model. The key idea for the GLM approach is the separation of the Lagrangian (fluid) and Eulerian (wave) degrees of freedom in Hamilton's principle. This is done by coupling an Euler-Poincaré reduced Lagrangian for the current flow and a phase-space Lagrangian for the wave field. WCI in the GLM model involves the nonlinear Doppler shift in frequency of the Hamiltonian wave subsystem, which arises because the waves propagate in the frame of motion of the Lagrangian-mean velocity of the current. In contrast, WCI in the CL model arises because the fluid velocity is defined relative to the frame of motion of the Stokes mean drift velocity, which is usually taken to be prescribed, time independent and driven externally. We compare the GLM and CL theories by placing them both into the general framework of a stochastic Hamilton's principle for a 3D Euler-Boussinesq (EB) fluid in a rotating frame. In other examples, we also apply the GLM and CL methods to add wave physics and stochasticity to the familiar 1D and 2D shallow water flow models. The differences in the types of stochasticity which arise for GLM and CL models can be seen by comparing the Kelvin circulation theorems for the two models. The GLM model acquires stochasticity in its Lagrangian transport velocity for the currents and also in its group velocity for the waves. However, the CL model is based on defining the Eulerian velocity in the integrand of the Kelvin circulation relative to the Stokes drift velocity induced by waves driven externally. Thus, the Kelvin theorem for the stochastic CL model can accept stochasticity in its both its integrand and in the Lagrangian transport velocity of its circulation loop. In an "Appendix", we also discuss dynamical systems analogues of WCI., (© The Author(s) 2020.)

Suppose the observations of Lagrangian trajectories for fluid flow in some physical situation can be modelled sufficiently accurately by a spatially correlated Itô 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ô noise in the data frame to the equivalent Stratonovich noise. However, the transformation from the Itô 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ô 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., Competing Interests: We declare we have no competing interest., (© 2020 The Authors.)

This paper exploits the theory of geometric gradient flows to introduce an alternative regularization of the thin-film equation valid in the case of large-scale droplet spreading-the geometric diffuse-interface method. The method possesses some advantages when compared with the existing models of droplet spreading, namely the slip model, the precursor-film method and the diffuse-interface model. These advantages are discussed and a case is made for using the geometric diffuse-interface method for the purpose of numerical simulations. The mathematical solutions of the geometric diffuse interface method are explored via such numerical simulations for the simple and well-studied case of large-scale droplet spreading for a perfectly wetting fluid-we demonstrate that the new method reproduces Tanner's Law of droplet spreading via a simple and robust computational method, at a low computational cost. We discuss potential avenues for extending the method beyond the simple case of perfectly wetting fluids., Competing Interests: We declare we have no competing interests., (© 2020 The Author(s).)

We derive a new variational principle, leading to a new momentum map and a new multisymplectic formulation for a family of Euler-Poincaré equations defined on the Virasoro-Bott group, by using the inverse map (also called 'back-to-labels' map). This family contains as special cases the well-known Korteweg-de Vries, Camassa-Holm and Hunter-Saxton soliton equations. In the conclusion section, we sketch opportunities for future work that would apply the new Clebsch momentum map with 2-cocycles derived here to investigate a new type of interplay among nonlinearity, dispersion and noise., Competing Interests: We declare we have no competing interests.

The un-reduction procedure introduced previously in the context of classical mechanics is extended to covariant field theory. The new covariant un-reduction procedure is applied to the problem of shape matching of images which depend on more than one independent variable (for instance, time and an additional labelling parameter). Other possibilities are also explored: nonlinear [Formula: see text]-models and the hyperbolic flows of curves.

We derive and study stochastic dissipative dynamics on coadjoint orbits by incorporating noise and dissipation into mechanical systems arising from the theory of reduction by symmetry, including a semidirect product extension. Random attractors are found for this general class of systems when the Lie algebra is semi-simple, provided the top Lyapunov exponent is positive. We study in details two canonical examples, the free rigid body and the heavy top, whose stochastic integrable reductions are found and numerical simulations of their random attractors are shown.

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 non-holonomic constraints. For this purpose, we derive, analyse 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, non-holonomically constrained systems, because it preserves the mathematical properties inherited from the variational principle., Competing Interests: The authors have no competing interests.

Inspired by spatiotemporal observations from satellites of the trajectories of objects drifting near the surface of the ocean in the National Oceanic and Atmospheric Administration's "Global Drifter Program", this paper develops data-driven stochastic models of geophysical fluid dynamics (GFD) with non-stationary spatial correlations representing the dynamical behaviour of oceanic currents. Three models are considered. Model 1 from Holm (Proc R Soc A 471:20140963, 2015) is reviewed, in which the spatial correlations are time independent. Two new models, called Model 2 and Model 3, introduce two different symmetry breaking mechanisms by which the spatial correlations may be advected by the flow. These models are derived using reduction by symmetry of stochastic variational principles, leading to stochastic Hamiltonian systems, whose momentum maps, conservation laws and Lie-Poisson bracket structures are used in developing the new stochastic Hamiltonian models of GFD.

Variational integrators are derived for structure-preserving simulation of stochastic Hamiltonian systems with a certain type of multiplicative noise arising in geometric mechanics. The derivation is based on a stochastic discrete Hamiltonian which approximates a type-II stochastic generating function for the stochastic flow of the Hamiltonian system. The generating function is obtained by introducing an appropriate stochastic action functional and its corresponding variational principle. Our approach permits to recast in a unified framework a number of integrators previously studied in the literature, and presents a general methodology to derive new structure-preserving numerical schemes. The resulting integrators are symplectic; they preserve integrals of motion related to Lie group symmetries; and they include stochastic symplectic Runge-Kutta methods as a special case. Several new low-stage stochastic symplectic methods of mean-square order 1.0 derived using this approach are presented and tested numerically to demonstrate their superior long-time numerical stability and energy behavior compared to nonsymplectic methods.

In Holm (Holm 2015 Proc. R. Soc. A 471 , 20140963. (doi:10.1098/rspa.2014.0963)), stochastic fluid equations were derived by employing a variational principle with an assumed stochastic Lagrangian particle dynamics. Here we show that the same stochastic Lagrangian dynamics naturally arises in a multi-scale decomposition of the deterministic Lagrangian flow map into a slow large-scale mean and a rapidly fluctuating small-scale map. We employ homogenization theory to derive effective slow stochastic particle dynamics for the resolved mean part, thereby obtaining stochastic fluid partial equations in the Eulerian formulation. To justify the application of rigorous homogenization theory, we assume mildly chaotic fast small-scale dynamics, as well as a centring condition. The latter requires that the mean of the fluctuating deviations is small, when pulled back to the mean flow., Competing Interests: We declare we have no competing interests.

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 detail examples with symmetric space SU (2)/ S 1 and SO (4)/ SO -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.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.

Vortex blob methods are typically characterized by a regularization length scale, below which the dynamics are trivial for isolated blobs. In this article, we observe that the dynamics need not be trivial if one is willing to consider distributional derivatives of Dirac delta functionals as valid vorticity distributions. More specifically, a new singular vortex theory is presented for regularized Euler fluid equations of ideal incompressible flow in the plane. We determine the conditions under which such regularized Euler fluid equations may admit vorticity singularities which are stronger than delta functions, e.g., derivatives of delta functions. We also describe the symplectic geometry associated with these augmented vortex structures, and we characterize the dynamics as Hamiltonian. Applications to the design of numerical methods similar to vortex blob methods are also discussed. Such findings illuminate the rich dynamics which occur below the regularization length scale and enlighten our perspective on the potential for regularized fluid models to capture multiscale phenomena.

This paper derives stochastic partial differential equations (SPDEs) for fluid dynamics from a stochastic variational principle (SVP). The paper proceeds by taking variations in the SVP to derive stochastic Stratonovich fluid equations; writing their Itô 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ô representation, because they are disguised by the quadratic covariation drift term arising in the Stratonovich to Itô transformation. This term is a geometric generalization 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; namely, the Euler-Boussinesq and quasi-geostropic approximations.

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.

We formulate the quasi-Lagrangian fluid transport dynamics of mass density ρ and the projection q=ω·∇ρ of the vorticity ω onto the density gradient, as determined by the three-dimensional compressible Navier-Stokes equations for an ideal gas, although the results apply for an arbitrary equation of state. It turns out that the quasi-Lagrangian transport of q cannot cross a level set of ρ. That is, in this formulation, level sets of ρ (isopycnals) are impermeable to the transport of the projection q.

A quantum spline is a smooth curve parametrized 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 minimized. 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.

In the framework of large deformation diffeomorphic metric mapping (LDDMM), we present a practical methodology to integrate prior knowledge about the registered shapes in the regularizing metric. Our goal is to perform rich anatomical shape comparisons from volumetric images with the mathematical properties offered by the LDDMM framework. We first present the notion of characteristic scale at which image features are deformed. We then propose a methodology to compare anatomical shape variations in a multi-scale fashion, i.e., at several characteristic scales simultaneously. In this context, we propose a strategy to quantitatively measure the feature differences observed at each characteristic scale separately. After describing our methodology, we illustrate the performance of the method on phantom data. We then compare the ability of our method to segregate a group of subjects having Alzheimer's disease and a group of controls with a classical coarse to fine approach, on standard 3D MR longitudinal brain images. We finally apply the approach to quantify the anatomical development of the human brain from 3D MR longitudinal images of pre-term babies. Results show that our method registers accurately volumetric images containing feature differences at several scales simultaneously with smooth deformations., (© 2011 British Crown Copyright)

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.

In this paper, we present a fine and coarse approach for the multiscale registration of 3D medical images using Large Deformation Diffeomorphic Metric Mapping (LDDMM). This approach has particularly interesting properties since it estimates large, smooth and invertible optimal deformations having a rich descriptive power for the quantification of temporal changes in the images. First, we show the importance of the smoothing kernel and its influence on the final solution. We then propose a new strategy for the spatial regularization of the deformations, which uses simultaneously fine and coarse smoothing kernels. We have evaluated the approach on both 2D synthetic images as well as on 3D MR longitudinal images out of the Alzheimer's Disease Neuroimaging Initiative (ADNI) study. Results highlight the regularizing properties of our approach for the registration of complex shapes. More importantly, the results also demonstrate its ability to measure shape variations at several scales simultaneously while keeping the desirable properties of LDDMM. This opens new perspectives for clinical applications.

The Camassa-Holm (CH) equation is a well-known integrable equation describing the velocity dynamics of shallow water waves. This equation exhibits spontaneous emergence of singular solutions (peakons) from smooth initial conditions. The CH equation has been recently extended to a two-component integrable system (CH2), which includes both velocity and density variables in the dynamics. Although possessing peakon solutions in the velocity, the CH2 equation does not admit singular solutions in the density profile. We modify the CH2 system to allow a dependence on the average density as well as the pointwise density. The modified CH2 system (MCH2) does admit peakon solutions in the velocity and average density. We analytically identify the steepening mechanism that allows the singular solutions to emerge from smooth spatially confined initial data. Numerical results for the MCH2 system are given and compared with the pure CH2 case. These numerics show that the modification in the MCH2 system to introduce the average density has little short-time effect on the emergent dynamical properties. However, an analytical and numerical study of pairwise peakon interactions for the MCH2 system shows a different asymptotic feature. Namely, besides the expected soliton scattering behavior seen in overtaking and head-on peakon collisions, the MCH2 system also allows the phase shift of the peakon collision to diverge in certain parameter regimes.

We investigate the emergence of singular solutions in a nonlocal model for a magnetic system. We study a modified Gilbert-type equation for the magnetization vector and find that the evolution depends strongly on the length scales of the nonlocal effects. We pass to a coupled density-magnetization model and perform a linear stability analysis, noting the effect of the length scales of nonlocality on the system's stability properties. We carry out numerical simulations of the coupled system and find that singular solutions emerge from smooth initial data. The singular solutions represent a collection of interacting particles (clumpons). By restricting ourselves to the two-clumpon case, we are reduced to a two-dimensional dynamical system that is readily analyzed, and thus we classify the different clumpon interactions possible.

We compute solutions of the Lagrangian-averaged Navier-Stokes alpha - (LANS alpha ) model for significantly higher Reynolds numbers (up to Re approximately 8300 ) than have previously been accomplished. This allows sufficient separation of scales to observe a Navier-Stokes inertial range followed by a second inertial range specific to the LANS alpha model. Both fully helical and nonhelical flows are examined, up to Reynolds numbers of approximately 1300. Analysis of the third-order structure function scaling supports the predicted l3 scaling; it corresponds to a k-1 scaling of the energy spectrum for scales smaller than alpha. The energy spectrum itself shows a different scaling, which goes as k1. This latter spectrum is consistent with the absence of stretching in the subfilter scales due to the Taylor frozen-in hypothesis employed as a closure in the derivation of the LANS alpha model. These two scalings are conjectured to coexist in different spatial portions of the flow. The l3 [E(k) approximately k-1] scaling is subdominant to k1 in the energy spectrum, but the l3 scaling is responsible for the direct energy cascade, as no cascade can result from motions with no internal degrees of freedom. We demonstrate verification of the prediction for the size of the LANS alpha attractor resulting from this scaling. From this, we give a methodology either for arriving at grid-independent solutions for the LANS alpha model, or for obtaining a formulation of the large eddy simulation optimal in the context of the alpha models. The fully converged grid-independent LANS alpha model may not be the best approximation to a direct numerical simulation of the Navier-Stokes equations, since the minimum error is a balance between truncation errors and the approximation error due to using the LANS alpha instead of the primitive equations. Furthermore, the small-scale behavior of the LANS alpha model contributes to a reduction of flux at constant energy, leading to a shallower energy spectrum for large alpha. These small-scale features, however, do not preclude the LANS alpha model from reproducing correctly the intermittency properties of the high-Reynolds-number flow.

New model equations are derived for dynamics of aggregation of finite-size particles. The differences from standard Debye-Hückel and Keller-Segel models are that the mobility of particles depends on the configuration of their neighbors and linear diffusion acts on locally averaged particle density. The evolution of collapsed states in these models reduces exactly to finite-dimensional dynamics of interacting particle clumps. Simulations show these collapsed (clumped) states emerge from smooth initial conditions, even in one spatial dimension. Extensions to two and three dimensions are also discussed.

We consider the wide class of systems modeled by an integrable approximation to the 3 degrees of freedom elastic pendulum with 1:1:2 resonance, or the swing-spring. This approximation has monodromy which prohibits the existence of global action-angle variables and complicates the dynamics. We study the quantum swing-spring formed by bending and symmetric stretching vibrations of the CO2 molecule. We uncover quantum monodromy of CO2 as a nontrivial codimension 2 defect of the three dimensional energy-momentum lattice of its quantum states.

Computational anatomy (CA) has introduced the idea of anatomical structures being transformed by geodesic deformations on groups of diffeomorphisms. Among these geometric structures, landmarks and image outlines in CA are shown to be singular solutions of a partial differential equation that is called the geodesic EPDiff equation. A recently discovered momentum map for singular solutions of EPDiff yields their canonical Hamiltonian formulation, which in turn provides a complete parameterization of the landmarks by their canonical positions and momenta. The momentum map provides an isomorphism between landmarks (and outlines) for images and singular soliton solutions of the EPDiff equation. This isomorphism suggests a new dynamical paradigm for CA, as well as new data representation.

The statistics of the energy and helicity fluxes in isotropic turbulence are studied using high resolution direct numerical simulation. The scaling exponents of the energy flux agree with those of the transverse velocity structure functions through refined similarity hypothesis, consistent with Kraichnan's prediction. The helicity flux is even more intermittent than the energy flux. Consistent with this observation, the spatial helicity-flux structures are finer than those of energy flux and more tubelike in geometry.

Elliptic instability in fluids is discussed in the context of the Lagrangian-averaged Navier-Stokes-alpha (LANS-alpha) turbulence model. This model preserves the Craik-Criminale (CC) family of solutions consisting of a columnar eddy and a Kelvin wave. The LANS-alpha model is shown to preserve elliptic instability. However, the model shifts the critical stability angle. This shift increases (decreases) the maximum growth rate for long (short) waves. It also introduces a band of stable CC solutions for short waves.

A vorticity surge event that could be a paradigm for a wide class of bursting events in turbulence is studied. The coherent mechanism is characterized by locally transverse vortex configurations that are intrinsically helical in both physical and Fourier space when there is a peak of the maximum vorticity parallel omega parallel(infinity)(t). At no time are nonhelical, antiparallel vorticity elements observed. This event precedes the appearance of the traditional signatures of an energy cascade such as strong growth of the dissipation, spectra approaching -5/3, and strongly Beltramized vortex tubes. Comparing how different large-eddy simulations reproduce these properties demonstrates the importance of properly modeling nonlinear transport of both energy and helicity.

We begin by placing the generalized Lagrangian mean (GLM) equations for a compressible adiabatic fluid into the Euler-Poincare (EP) variational framework of fluid dynamics, for an averaged Lagrangian. This is the Lagrangian averaged Euler-Poincare (LAEP) theorem. Next, we derive a set of approximate small amplitude GLM equations (glm equations) at second order in the fluctuating displacement of a Lagrangian trajectory from its mean position. These equations express the linear and nonlinear back-reaction effects on the Eulerian mean fluid quantities by the fluctuating displacements of the Lagrangian trajectories in terms of their Eulerian second moments. The derivation of the glm equations uses the linearized relations between Eulerian and Lagrangian fluctuations, in the tradition of Lagrangian stability analysis for fluids. The glm derivation also uses the method of averaged Lagrangians, in the tradition of wave, mean flow interaction. Next, the new glm EP motion equations for incompressible ideal fluids are compared with the Euler-alpha turbulence closure equations. An alpha model is a GLM (or glm) fluid theory with a Taylor hypothesis closure. Such closures are based on the linearized fluctuation relations that determine the dynamics of the Lagrangian statistical quantities in the Euler-alpha equations. Thus, by using the LAEP theorem, we bridge between the GLM equations and the Euler-alpha closure equations, through the small-amplitude glm approximation in the EP variational framework. We conclude by highlighting a new application of the GLM, glm, and alpha-model results for Lagrangian averaged ideal magnetohydrodynamics. (c) 2002 American Institute of Physics.

We use asymptotic analysis and a near-identity normal form transformation from water wave theory to derive a 1+1 unidirectional nonlinear wave equation that combines the linear dispersion of the Korteweg-deVries (KdV) equation with the nonlinear/nonlocal dispersion of the Camassa-Holm (CH) equation. This equation is one order more accurate in asymptotic approximation beyond KdV, yet it still preserves complete integrability via the inverse scattering transform method. Its traveling wave solutions contain both the KdV solitons and the CH peakons as limiting cases.