On Nonlinear Baroclinic Waves and Adjustment of Pancake Dynamics

Three-dimensional nonhydrostatic Euler–Boussinesq equations are studied for Bu=O(1) flows as well as in the asymptotic regime of strong stratification and weak rotation. Reduced prognostic equations for ageostrophic components (divergent velocity potential and geostrophic departure/thermal wind imbalance) are analyzed. We describe classes of nonlinear anisotropic ageostrophic baroclinic waves which are generated by the strong nonlinear interactions between the quasi-geostrophic modes and inertio-gravity waves. In the asymptotic regime of strong stratification and weak rotation we show how switching on weak rotation triggers frontogenesis. The mechanism of the front formation is contraction in the horizontal dimension balanced by vertical shearing through coupling of large horizontal and small vertical scales by weak rotation. Vertical slanting of these fronts is proportional to μ^−1/2 where μ is the ratio of the Coriolis and Brunt–V?is?l? parameters. These fronts select slow baroclinic waves through nonlinear adjustment of the horizontal scale to the vertical scale by weak rotation, and are the envelope of inertio-gravity waves. Mathematically, this is generated by asymptotic hyperbolic systems describing the strong nonlinear interactions between waves and potential vorticity dynamics. This frontogenesis yields vertical “gluing” of pancake dynamics, in contrast to the independent dynamics of horizontal layers in strongly stratified turbulence without rotation. Received 8 April 1997 and accepted 29 March 1998     

Coherent vorticity dynamics and dissipation in a utility-scale wind turbine wake with uniform inflow

The vorticity dynamics and its relationship to dissipation in the wake of a utility-scale wind turbine are investigated through large-eddy simulation. The vorticity dynamics is assessed through the enstrophy, which is related to the turbulent dissipation. The averaged enstrophy and turbulent dissipation are shown to be quantitatively similar in the wake. Using temporal phase averaging, the vorticity fluctuations are decomposed into coherent and random fluctuations with respect to the frequency of the tip vortices. The enstrophy in the tip vortices is dominated by coherent fluctuations, while the coherent fluctuations of root vortices are immediately saturated by the random vorticity fluctuations of the unstable hub vortex. The coherent strain rate has significant differences compared to the coherent enstrophy within one diameter downwind of blade tip, but the random enstrophy and strain rate are relatively similar. Differences in coherent enstrophy and strain rate decrease further from the rotor.     

The effect of fine structure on the stability of planar vortices

This study considers the linear, inviscid response to an external strain field of classes of planar vortices. The case of a Gaussian vortex has been considered elsewhere, and an enstrophy rebound phenomenon was noted: after the vortex is disturbed enstrophy feeds from the non-axisymmetric to mean flow. At the same time an irreversible spiral wind-up of vorticity fluctuations takes place. A top-hat or Rankine vortex, on the other hand, can support a non-decaying normal mode.In vortex dynamics processes such as stripping and collisions generate vortices with sharp edges and often with bands or rings of fine scale vorticity at their periphery, rather than smooth profiles. This paper considers the stability and response of a family of vortices that vary from a broad profile to a top-hat vortex. As the edge of the vortex becomes sharper, a quasi-mode emerges and vorticity winds up in a critical layer, at the radius where the angular velocity of the fluid matches that of a normal mode on a top-hat vortex. The decay rate of these quasi-modes is proportional to the vorticity gradient at the critical layer, in agreement with theory. As the vortex edge becomes sharper it is found that the rebound of enstrophy becomes stronger but slower.The stability and linear behaviour of coherent vortices is then studied for distributions which exhibit additional fine structure within the critical layer. In particular we consider vorticity profiles with ‘bumps’, ‘troughs’ or ‘steps’ as this fine structure. The modified evolution equation that governs the critical layer is studied using numerical simulations and asymptotic analysis. It is shown that depending on the form of the short-scale vorticity distribution, this can stabilise or destabilise quasi-modes, and it may also lead to oscillatory behaviour.     

Decay of confined,two-dimensional,spatially periodic arrays of vortices: A numerical investigation

The disarrangement of a perturbed lattice of vortices was studied numerically. The basic state is an exponentially decaying, exact solution of the Navier-Stokes equations. Square arrays of vortices with even numbers of vortex cells along each side were perturbed and their evolution was investigated. Whether the energy in the perturbation grows somewhat before it decays or decays monotonically depends on the initial strength of the vortices of the basic state, the extent of lateral confinement and the structure of the perturbation. The critical condition for temporally local instability, i.e. the critical amplitude of the basic state that must be exceeded to allow energy transfer from the basic state to the perturbation, is discussed. In the strongly confined case of a square lattice of four vortices the appearance of enchancement of global rotation is the result of energy transfer from the basic state to a temporally local unstable mode. Energy is transferred from the basic state to larger-scaled structures (inverse cascade) only if the scales of the larger structures are inherently contained in the initial structure of the perturbation. The initial structure of the double array of vortices is not maintained except for a very special form of perturbation. The facts that large scales decay more slowly than small scales and that, when non-linearities are sufficiently strong, energy is transferred from one scale to another explain the differences in the disarrangement process for different initial strengths of the vortices of the basic state. The stronger vortices, i.e. the vortices perturbed in a manner that increases their strength, tend to dominate the weaker vortices. The pairing and subsequent merging (or capture) of vortices of like sense into larger-scale vortices are described in terms of peaks in the evolution of the square root of the palinstrophy divided by the enstrophy.     

Direct numerical simulation of homogeneous stratified rotating turbulence

The effects of the Prandtl number on stratified rotating turbulence have been studied in homogeneous turbulence by using direct numerical simulations and a rapid distortion theory. Fluctuations under strong stable-density stratification can be theoretically divided into the WAVE and the potential vorticity (PV) modes. In low-Prandtl-number fluids, the WAVE mode deteriorates, while the PV mode remains. Imposing rotation on a low-Prandtl-number fluid makes turbulence two-dimensional as well as geostrophic; it is found from the instantaneous turbulent structure that the vortices merge to form a few vertically-elongated vortex columns. During the period toward two-dimensionalization, the vertical vortices become asymmetric in the sense of rotation. Communicated by S. Obi PACS 47.55.Hd     

Flow structures,nonlinear inertial waves and energy transfer in rotating spheres

We investigate flow structures, nonlinear inertial waves and energy transfer in a rotating fluid sphere, using a Galerkin spectral method based on helical-wave decomposition (HWD). Numerical simulations of flows in a sphere are performed with different system rotation rates, where a large-scale forcing is employed. For the case without system rotation, the intense vortex structures are tube-like. When a weak rotation is introduced, small-scale structures are reduced and vortex tubes tend to align with the rotation axis. As the rotation rate increases, a large-scale anticyclonic vortex structure is formed near the rotation axis. The structure is shown to be led by certain geostrophic modes. When the rotation rate further increases, a cyclone and an anticyclone emerge from the top and bottom of the boundary, respectively, where two quasi-geostrophic equatorially symmetric inertial waves dominate the flow. Based on HWD, effects of spherical confinement on rotating turbulence are systematically studied. It is found that the forward cascade becomes weaker as the rotation increases. When the rotation rate becomes larger than some critical value, dual energy cascades emerge, with an inverse cascade at large scales and a forward cascade at small scales. Finally, the flow behavior near the boundary is studied, where the average boundary layer thickness gets smaller when system rotation increases. The flow behavior in the boundary layer is closely related to the interior flow structures, which create significant mass flux between the boundary layer and the interior fluid through Ekman pumping.     

Decaying homogeneous turbulence under strong density stratification at low-Prandtl number

The characteristics of decaying homogeneous turbulence under strong density stratification have been studied using direct numerical simulations. While our previous study dealt with rotating stratified turbulence, here we investigate the detailed flow structure of stratified turbulence without rotation especially at low-Prandtl number. By assuming a low-Prandtl-number fluid, e.g. liquid sodium: Pr ≈ 0.01, gallium: Pr ≈ 0.025, internal gravity waves are markedly attenuated due to the large thermal conductivity, and turbulence soon reaches a two-component state, where vertical energy, coupled with potential energy, significantly decays, and becomes negligible as observed experimentally (Praud et al. in J Fluid Mech 522:1–33, 2005). In the horizontal plane, there appear large-scale vortices with vertical vorticity, and those with the same sign of vorticity increase their horizontal length scale by merging with each other. In the vertical plane, highly sheared regions represented by horizontal vorticity also tend to horizontally increase their length scale and become layered structures by the combined effects of vortex coalescence and energy cascade into higher vertical wavenumbers.      

Evolution of solitary vortices in a rotating fluid

The evolution of cyclones and anticyclones with a characteristic scale substantially exceeding the radius of deformation is investigated numerically. An equation obtained by an asymptotic method for small values of the Kibel'-Rossby number is employed. The model ensures conservation of the potential vorticity in the fluid particles for motions of finite amplitude (when the thickness of the layer H deviates considerably from the undisturbed value H₀). It is shown that anticyclones of a certain type adapt themselves to a steady shape and travel in a direction opposite to the direction of rotation of the sphere (westwards). It is found that the existence of a steadily migrating anticyclone requires the presence of a region of closed isolines of the potential vorticity (trapping region), in which the fluid is transported together with the anticyclone. Cyclones drift westwards more slowly than anticyclones and migrate towards the poles, losing energy by emitting Rossby waves. It is found that in cyclones the trapping region retains an almost circular shape, the decay of the eddy, associated with wave emission, showing as its initial intensity increases.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 52–59, July–August, 1986.     

Direct simulation of breakdown to turbulence following oblique instability waves in a supersonic boundary layer

The late stages of transition to turbulence in a Mach two boundary layer are investigated by direct numerical simulation of the compressible Navier-Stokes equations. The primary instability at this Mach number consists of oblique waves, which are known to form a pattern of quasi-streamwise vortices. It is found that breakdown does not follow immediately from these vortices, which decay in intensity. The generation of new vortices is observed by following the evolution of the pressure and vorticity in the simulation, and analysed by consideration of vorticity stretching. It is found that the slight inclined and skewed nature of the quasi-streamwise vortices leads to a production of oppositely signed streamwise vorticity, which serves as a strong localised forcing of the shear layer alongside the original vortices, formed by convection and stretching of spanwise vorticity. The shear layer rolls up into many new vortices, and is followed by a sharp increase in the energy of higher frequencies and in the skin friction.     

A New Class of Equations for Rotationally Constrained Flows

The incompressible Navier–Stokes equation is considered in the limit of rapid rotation (small Ekman number). The analysis is limited to horizontal scales small enough so that both horizontal and vertical velocities are comparable, but the horizontal velocity components are still in geostrophic balance. Asymptotic analysis leads to a pair of nonlinear equations for the vertical velocity and vertical vorticity coupled by vertical stretching. Statistically stationary states are maintained against viscous dissipation by boundary forcing or energy injection at larger scales. For thermal forcing direct numerical simulation of the reduced equations reveals the presence of intense vortical structures spanning the layer depth, in excellent agreement with simulations of the Boussinesq equations for rotating convection by Julien et al. (1996). Received 30 May 1997 and accepted 4 January 1998     

Lagrangian Covariance Analysis of β-Plane Turbulence

The effects of Rossby wave–turbulence interactions on particle dispersion are investigated in a Lagrangian analysis of the potential vorticity equation. The analysis produces several exact statistical results for fluid particle dispersion in barotropic turbulence on a β-plane. In the inviscid problem the first integral time scale of the meridional velocity is found to be zero, as might occur in pure wave processes, and the meridional particle dispersion is bounded. The second integral time scale, which determines the magnitude of the bound, is shown to depend explicitly on β, the enstrophy and the energy of the meridional velocity. Expressions relating the autocorrelation of the vorticity to the autocorrelation of the meridional velocity are derived and the Lagrangian integral time scale of the relative vorticity is diagnostically related to the meridional velocity correlation. The applicability of these predictions is verified in a series of numerical simulations. For a range of β values, the meridional extent of quasisteady alternating zonally averaged jets occurring in the numerical solutions scales with a length scale given by the the standard deviation of the meridional particle dispersion. Received 8 March 1999 and accepted 12 December 1999     

Vortex structures in geophysical convection

An intriguing variety of vortex structures arise during buoyant convection, especially in the presence of background stratification and rotation. These vortices play an important role in environmental fluid motions, bearing upon small-scale turbulence to planetary-scale circulation. A brief review of vortex motions associated with buoyant convection is presented in this paper, emphasizing the sources of vorticity, evolution of vortex structures and their role in oceanic and atmospheric dynamics. The genesis of a variety of vortices, for example, mushroom vortices, geostrophic and ageostrophic vortices, dipolar structures and hetons in buoyant convection flows is described, and parameterizations to represent their properties are discussed. New laboratory and numerical simulation results on vortex-related phenomena in stratified and rotating fluids and their implications in geophysical convective flows are also presented.     

The N-vortex problem on a sphere: geophysical mechanisms that break integrability

We describe the dynamical system governing the evolution of a system of point vortices on a rotating spherical shell, highlighting features which break what would otherwise be an integrable problem. The importance of the misalignment of the center-of-vorticity vector associated with a cluster of point vortices with the axis of rotation is emphasized as a crucial factor in the interpretation of dynamical features for many flow configurations. We then describe two important physical mechanisms which break what would otherwise be an integrable problem—the interactions between the local center-of-vorticity vectors of more than one region of concentrated vorticity, and the coupling between the center-of-vorticity vector and the background vorticity field which supports Rossby waves. Focusing on the Polar vortex splitting event of September 2002, we describe simple (i.e., low dimensional) mechanisms that can trigger instabilities whose subsequent development cause the onset of chaotic advection and global particle transport. At the linear level, eigenvalues that oscillate between elliptic and hyperbolic configurations initiate the pinch-off process of a passive patch representing the Polar vortex. At the nonlinear level, the evolution and topological bifurcations of the streamline patterns are responsible for its further splitting, stretching, and subsequent transport over the sphere. We finish by briefly describing how to incorporate conservation of potential vorticity and the development of a model governing the probability density function associated with the point vortex system.     

On the Asymptotic Regimes and the Strongly Stratified Limit of Rotating Boussinesq Equations

Asymptotic regimes of geophysical dynamics are described for different Burger number limits. Rotating Boussinesq equations are analyzed in the asymptotic limit of strong stratification in the Burger number of order one situation as well as in the asymptotic regime of strong stratification and weak rotation. It is shown that in both regimes the horizontally averaged buoyancy variable is an adiabatic invariant (approximate conservation law) for the full Boussinesq system. Spectral phase shift corrections to the buoyancy time scale associated with vertical shearing of this invariant are deduced. Statistical dephasing effects induced by turbulent processes on inertial-gravity waves are evidenced. The “split” of the energy transfer of the vortical and the wave components is established in the Craya–Herring cyclic basis. As the Burger number increases from zero to infinity, we demonstrate gradual unfreezing of energy cascades for ageostrophic dynamics. This property is related to the nonlinear geostrophic adjustment mechanism which is the capacity of ageostrophic dynamics to transfer energy to small scales. The energy spectrum and the anisotropic spectral eddy viscosity are deduced with an explicit dependence on the anisotropic rotation/stratification time scale which depends on the vertical aspect ratio parameter. Intermediate asymptotic regime corresponding to strong stratification and weak rotation is analyzed where the effects of weak rotation are accounted for by an asymptotic expansion with full control (saturation) of vertical shearing. The regularizing effect of weak rotation differs from regularizations based on vertical viscosity. Two scalar prognostic equations for ageostrophic components (divergent velocity potential and geostrophic departure) are obtained. Received 23 January 1997 and accepted 11 July 1997     

On the correlation between enstrophy and energy dissipation rate in a turbulent wake

All three components of the vorticity fluctuation have been measured simultaneously in a turbulent wake using a new eight-sensor vorticity probe. The vorticity fluctuation spectra agree reasonably well with those from a direct numerical simulation of a turbulent channel flow at high wavenumbers. The similarity between the instantaneous energy dissipation rate ε and the instantaneous enstrophy ω² is examined using spectra and probability density functions. The correlation between ω² and ε is evaluated in some detail. The homogeneous value of ε is strongly correlated with ω². The full value of ε and, more especially its isotropic value, are less well correlated with the enstrophy. Conditional averaging indicates that high enstrophy regions are associated with high energy dissipation rate regions.     

Nonlinear interaction of a vortex pair with clean and surfactant-covered free surfaces

Two-dimensional unsteady viscous-flow problem associated with the normal incidence of a counter-rotating vortex pair on a free surface is analyzed. Effects of surface tension and insoluble surfactants on the generation of free-surface vorticity and surface waves are investigated. A recently developed finite-difference method based on boundary-fitted coordinates is used to solve the fully-nonlinear problem. Results show that in the absence of surfactants and at low Froude number (based on circulation strength and initial separation distance of the vortex pair), waves of short lengths are generated. However, secondary vorticity generated in this case is not strong enough to affect the outward translation of the primary vortices. At intermediate Froude number, a transient wave developing outboard of the primary vortex becomes steep, and eventually breaks because of local instability. Consequently, free-surface vorticity inhibits the outward translation of the primary vortices. Surface tension in a clean free surface dampens the steep short waves, hence also the generation of free-surface vorticity. However, variation in surface tension induced by surfactants intensifies the generation of surface vorticity, thereby causing the primary vortices to rebound. The increase in the rotational part of wave motion results in the dampening of overall free-surface deformations. However, it is found that the shear stress associated with a large gradient of surfactant concentration could cause local steepening of the short wave generated outboard of the primary vortex.     

Velocity field after wave breaking

The velocity field in breaking water waves is considered in this paper. A numerical simulation describes in detail the transition from a primary overturning and consequent rebounding jets into a bore front, where the vorticity in the coherent large‐scale eddy structures devolves into turbulence. Spatial changes in the frequency spectra of the kinetic energy and the enstrophy are associated with the production, transport and dissipation of the Reynolds stress and the various wave and turbulent mixing length scales. Mean velocity fields and the wave and kinetic energy in a surf zone are evaluated. Fourier and wavelet spectral analysis is applied to study both the surface elevation and energy changes, and the distinction that must be made between spilling and plunging breakers is clarified in this paper. Copyright © 2002 John Wiley & Sons, Ltd.     

Pressure drop for flow in channels subjected to strong system rotation

A disc stack centrifuge is an industrial example of a fluid machine in which all the internal flow takes place in a rapidly rotating frame. The present report gives a survey of the experimental and theoretical work performed at Alfa-Laval in order to estimate the pressure drops in the different internal passages in the centrifuge, including both laminar and turbulent flow.For the laminar flow between the discs, a theory has been developed using the concept of a rotating Hele Shaw cell and conformal mapping. The theory is valid in the limit of very small Rossby numbers. For moderately large Rossby numbers, this model overestimates the pressure drop. The linear theory was extended by introducing advecting vortices in a computer model. The vortices cause vertical fluid transport between the Ekman and geostrophic layers by Ekman pumping, an effect which decreases the pressure drop in the disc stack. The linear model and the enhanced model have both been confirmed by experiments.The flow is turbulent in most parts of the centrifuge, except in the disc stack. The theoretical or numerical modelling for rotating turbulent flows is very difficult and no reliable models exist so far. We therefore have to rely on measurements, which show that the pressure is significantly influenced by rotation for Rossby numbers below unity.     

Vortex stretching and filaments

This paper presents a new technique to produce controlled stretched vortices. Intense elliptical vortices are created by stretching of an initial vorticity sheet. The initial vorticity comes from a laminar boundary layer flow and the stretching is parallel to the vorticity vectors. This low velocity flow enables direct observation of the formation and destabilization of vortices. Visualizations are combined with quasi-instantaneous measurements of a full velocity profile. The velocity profile is obtained with an ultrasonic pulsed Doppler velocimeter. The evolution of the central diameter of the vortices is related to the stretching. It is observed that destabilization occurs by pairing of two vortices, by hairpin deformation, and by breakdown of vortices into a “coil shape”.     

Understanding coherent vortices through computational fluid dynamics

A new definition of coherent vortices in turbulence is proposed, where the vorticity equation reduces to a cyclostrophic balance. Afterward, we describe five fundamental vortex interactions, the sheet, the spiral, the pairing, the even longitudinal, and the odd longitudinal modes. Numerous examples of these interactions are provided from direct numerical or large-eddy simulations. The resulting vortices are responsible for the internal intermittent character of turbulence, with highly nongaussian tails for the probability density functions of vorticity, passive scalar, and low pressure. In a mixing layer, the combination of the odd longitudinal and the pairing modes (helical pairing) is inhibited by compressibility, above a convective Mach number of 0.7. When turbulence is submitted to a solid-body rotation, anticyclonic vortices of local Rossby number of the order of 1 transform into intense perpendicular Görtler-type alternate longitudinal vortices.The support of CCVR, CEA, CNRS, Dassault/CNES, DRET, LHF, and Région Rhône-Alpes is acknowledged. This paper is the text of an invited lecture given at the IUTAM Symposium on Eddy Structure Identification in Free Turbulent Shear Flows, Poitiers, 12–14 October 1992.