Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userStochastic Parametrization of the Richardson Triple
69
0
0.0
(
0
)
Authors
Darryl D. Holm
Ask ChatGPT about the research
No Arabic abstract
A Richardson triple is an ideal fluid flow map $g_{t/ep,t,ep t} = h_{t/ep}k_t l_{ep t}$ composed of three smooth maps with separated time scales: slow, intermediate and fast; corresponding to the big, little, and lesser whorls in Richardsons well-known metaphor for turbulence. Under homogenisation, as $lim epto0$, the composition $h_{t/ep}k_t $ of the fast flow and the intermediate flow is known to be describable as a single stochastic flow $dd g$. The interaction of the homogenised stochastic flow $dd g$ with the slow flow of the big whorl is obtained by going into its non-inertial moving reference frame, via the composition of maps $(dd g)l_{ep t}$. This procedure parameterises the interactions of the three flow components of the Richardson triple as a single stochastic fluid flow in a moving reference frame. The Kelvin circulation theorem for the stochastic dynamics of the Richardson triple reveals the interactions among its three components. Namely, (i) the velocity in the circulation integrand acquires is kinematically swept by the large scales; and (ii) the velocity of the material circulation loop acquires additional stochastic Lie transport by the small scales. The stochastic dynamics of the composite homogenised flow is derived from a stochastic Hamiltons principle, and then recast into Lie-Poisson bracket form with a stochastic Hamiltonian. Several examples are given, including fluid flow with stochastically advected quantities, and rigid body motion under gravity, i.e., the stochastic heavy top in a rotating frame.
rate research
Read More
We are modelling multi-scale, 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 Hamiltons principle. This is done by coupling an Euler--Poincare {it reduced Lagrangian} for the current flow and a {it 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 Hamiltons 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. The Kelvin circulation theorem stochastic CL model can accept stochasticity in its both its integrand and in the Lagrangian transport velocity of its circulation loop.
We review the role of dual pairs in mechanics and use them to derive particle-like solutions to regularized incompressible fluid systems. In our case we have a dual pair resulting from the action of diffeomorphisms on point particles (essentially by moving the points). We then augment our dual pair by considering the action of diffeomorphisms on Taylor series, also known as jets. The augmented weak dual pairs induce a hierarchy of particle-like solutions and conservation laws with particles carrying a copy of a jet group. We call these augmented particles jetlets. The jet groups serve as finite-dimensional models of the diffeomorphism group itself, and so the jetlet particles serve as a finite-dimensional model of the self-similarity exhibited by ideal incompressible fluids. The conservation law associated to jetlet solutions is shown to be a shadow of Kelvins circulation theorem. Finally, we study the dynamics of infinite time particle mergers. We prove that two merging particles at the zeroth level in the hierarchy yield dynamics which asymptotically approach that of a single particle in the first level in the hierarchy. This merging behavior is then verified numerically as well as the exchange of angular momentum which must occur during a near collision of two particles. The resulting particle-like solutions suggest a new class of meshless methods which work in dimensions $n geq 2$ and which exhibit a shadow of Kelvins circulation theorem. More broadly, this provides one of the first finite-dimensional models of self-similarity in ideal fluids.
The classic evolution equations for potential flow on the free surface of a fluid flow are not closed because the pressure and the vertical velocity dynamics are not specified on the free surface. Moreover, their wave dynamics does not cause circulation of the fluid velocity on the free surface. The equations for free-surface motion we derive here are closed and they are not restricted to potential flow. Hence, true wave-current interaction dynamics can occur. In particular, the Kelvin-Noether theorem demonstrates that wave activity can induce fluid circulation and vorticity dynamics on the free surface. The wave-current interaction equations introduced here open new vistas for both the deterministic and stochastic analysis of nonlinear waves on free surfaces.
Wave--current interaction (WCI) dynamics energizes and mixes the ocean thermocline by producing a combination of Langmuir circulation, internal waves and turbulent shear flows, which interact over a wide range of time scales. Two complementary approaches exist for approximating different aspects of WCI dynamics. These are the Generalized Lagrangian Mean (GLM) approach and the Gent--McWilliams (GM) approach. Their complementarity is evident in their Kelvin circulation theorems. GLM introduces a wave pseudomomentum per unit mass into its Kelvin circulation integrand, while GM introduces a an additional `bolus velocity to transport its Kelvin circulation loop. The GLM approach models Eulerian momentum, while the GM approach models Lagrangian transport. In principle, both GLM and GM are based on the Euler--Boussinesq (EB) equations for an incompressible, stratified, rotating flow. The differences in their Kelvin theorems arise from differences in how they model the flow map in the Lagrangian for the Hamilton variational principle underlying the EB equations. A recently developed approach for uncertainty quantification in fluid dynamics constrains fluid variational principles to require that Lagrangian trajectories undergo Stochastic Advection by Lie Transport (SALT). Here we introduce stochastic closure strategies for quantifying uncertainty in WCI by adapting the SALT approach to both the GLM and GM approximations of the EB variational principle. In the GLM framework, we introduce a stochastic group velocity for transport of wave properties, relative to the frame of motion of the Lagrangian mean flow velocity and a stochastic pressure contribution from the fluctuating kinetic energy. In the GM framework we introduce a stochastic bolus velocity in addition to the mean drift velocity by imposing the SALT constraint in the GM variational principle.
The stochastic variational approach for geophysical fluid dynamics was introduced by Holm (Proc Roy Soc A, 2015) as a framework for deriving stochastic parameterisations for unresolved scales. This paper applies the variational stochastic parameterisation in a two-layer quasi-geostrophic model for a beta-plane channel flow configuration. We present a new method for estimating the stochastic forcing (used in the parameterisation) to approximate unresolved components using data from the high resolution deterministic simulation, and describe a procedure for computing physically-consistent initial conditions for the stochastic model. We also quantify uncertainty of coarse grid simulations relative to the fine grid ones in homogeneous (teamed with small-scale vortices) and heterogeneous (featuring horizontally elongated large-scale jets) flows, and analyse how the spread of stochastic solutions depends on different parameters of the model. The parameterisation is tested by comparing it with the true eddy-resolving solution that has reached some statistical equilibrium and the deterministic solution modelled on a low-resolution grid. The results show that the proposed parameterisation significantly depends on the resolution of the stochastic model and gives good ensemble performance for both homogeneous and heterogeneous flows, and the parameterisation lays solid foundations for data assimilation.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
