O. A. Druzhinin

Regular and chaotic dynamics of a fountain in a stratified fluid

In the present paper, we study by direct numerical simulation (DNS) and theoretical analysis, the dynamics of a fountain penetrating a pycnocline (a sharp density interface) in a density-stratified fluid. A circular, laminar jet flow of neutral buoyancy is considered, which propagates vertically upwards towards the pycnocline level, penetrates a distance into the layer of lighter fluid, and further stagnates and flows down under gravity around the up-flowing core thus creating a fountain. The DNS results show that if the Froude number (Fr) is small enough, the fountain top remains axisymmetric and steady. However, if Fr is increased, the fountain top becomes unsteady and oscillates in a circular flapping (CF) mode, whereby it retains its shape and moves periodically around the jet central axis. If Fr is increased further, the fountain top rises and collapses chaotically in a bobbing oscillation mode (or B-mode). The development of these two modes is accompanied by the generation of different patterns of internal waves (IW) in the pycnocline. The CF-mode generates spiral internal waves, whereas the B-mode generates IW packets with a complex spatial distribution. The dependence of the amplitude of the fountain-top oscillations on Fr is well described by a Landau-type two-mode-competition model.

Internal wave generation by a fountain in a stratified fluid

The purpose of the study is the direct numerical and theoretical modeling of fountain dynamics in a fluid with density stratification in the form of a pycnocline. The fountain is formed as a vertical jet penetrates through the pycnocline. In numerical simulation the jet flow is initiated by means of preassigning a boundary condition in the form of an upward-directed laminar flow of a neutral-buoyancy fluid with an axisymmetric Gaussian velocity profile. The calculations show that at a Froude number Fr greater than a certain critical value the flow becomes unstable and the fountain executes self-oscillations accompanied by internal wave generation in the pycnocline. Depending on Fr, two self-oscillation modes can be distinguished. At fairly low Fr the fountain executes circular motion in the horizontal plane, in the vicinity of the center of jet, its shape remaining almost invariant. In this case, internal waves in the form of unwinding spirals are radiated. At fairly high Fr another mode predominates, when the fountain top chaotically “strays” in the vicinity of the center of the jet and, periodically breaking down, generates wave packets propagating toward the periphery of the computation domain. In both cases, the main peak in the frequency spectrum of the internal waves coincides with the fountain top oscillation frequency which monotonically decreases with increase in Fr. In numerical simulation the Fr-dependence of the fountain top oscillation amplitude is in good agreement with that predicted by the theoretical model of the concurrence of the interacting modes in the soft self-excitation regime.

Internal Wave Radiation by a Turbulent Fountain in a Stratified Fluid

Large eddy simulation is applied to model a fountain in a density-stratified fluid. The fountain is formed, as a vertical turbulent jet penetrates through a pycnocline. The jet flow is initiated by the formulation of a boundary condition in the form of an upward neutral-buoyancy fluid flow with the Gaussian axisymmetric mean-velocity profile and a given fluctuation level. It is shown that at a Froude number Fr higher than a certain critical value the fountain executes self-oscillations accompanied by internal wave generation within the pycnocline. The predominant self-oscillation mode is axisymmetric, when the fountain top periodically breaks down generating internal wave packets traveling toward the periphery of the computation domain. The characteristic frequency of the internal waves coincides with that of the fountain top oscillations and monotonically decreases with increase in Fr. The Fr-dependence of the fountain top oscillation amplitude obtained in the numerical solution is in good agreement with the predictions of the theoretical Landau model for the instability mode in the soft self-excitation regime.

On the onset of the instability of a three-dimensional jet in a stratified fluid

The onset of a three-dimensional jet flow in a stratified fluid is studied with the aid of a direct numerical simulation. An initially cylindrical jet with a Gaussian velocity profile is considered in a fluid with stable linear density stratification. The results indicate that, if an initial small perturbation of the velocity field has a wide spectrum, an exponential growth of the isolated quasi-two-dimensional mode occurs and its spectral maximum is shifted toward smaller wave numbers in comparison with the maximum of the helical mode of the instability of a nonstratified jet. The growth rate is proportional to Ri0.5, where Ri is the global Richardson number. The onset of the instability leads to the formation of the flow’s vortex structure, which consists of a collection of different-polarity quasi-two-dimensional vortices located in a horizontal plane near the longitudinal axis of the jet. At sufficiently long times (Nt > 100, where N is the buoyancy frequency and t is time), the growth of instability reaches the saturation stage and further fluctuations in velocity and density decay under the effect of viscous diffusion. At this stage, the flow becomes self-similar and the time dependences of the transverse and vertical widths of the jet are consistent with the asymptotic behaviors of integral parameters of the flow that are observed experimentally in the far stratified wake. The results suggest that the onset of the instability of a quasitwo-dimensional mode can play the determining role in the dynamics of flow in the far stratified wake.

Laboratory, numerical, and theoretical modeling of the flow in a far wake in a stratified fluid

The far-wake flow past a sphere towed in a fluid with high Reynolds and Froude numbers and with a pycnocline-form salt-density stratification is studied in a laboratory experiment based on particle image velocimetry and in numerical and theoretical modeling. In the configuration under consideration, the axis of sphere towing is located under a pycnocline. Flow parameters, the profiles of density and average velocity, and the initial field of velocity fluctuation in numerical modeling are specified from the data of the laboratory experiment. The fields of fluid velocity at different times and the time dependences of integral parameters of wake flow, such as the average velocity at the axis and the transverse width of the flow, are obtained. The results of numerical modeling are in good qualitative and quantitative agreement with the data of the laboratory experiment. The results of the laboratory experiment and numerical modeling are compared to the predictions of a quasi-linear and quasi-two-dimensional theoretical model. The time evolution of both the average velocity at the axis and the transverse width of the wake is obtained with the model and is in good agreement with the experimental data. The results of numerical modeling also show that, under the effect of velocity fluctuation in the wake, internal waves whose spatial period is equal to the characteristic period of the wake’s vortex structure are excited efficiently in the pycnocline.

The study of the effects of sea-spray drops on the marine atmospheric boundary layer by direct numerical simulation

The detailed knowledge of turbulent exchange processes occurring in the atmospheric marine boundary layer are of primary importance for their correct parameterization in large-scale prognostic models. These processes are complicated, especially at sufficiently strong wind forcing conditions, by the presence of sea-spray drops which are torn off the crests of sufficiently steep surface waves by the wind gusts. Natural observations indicate that mass fraction of sea-spray drops increases with wind speed and their impact on the dynamics of the air in the vicinity of the sea surface can become quite significant. Field experiments, however, are limited by insufficient accuracy of the acquired data and are in general costly and difficult. Laboratory modeling presents another route to investigate the spray-mediated exchange processes in much more detail as compared to the natural experiments. However, laboratory measurements, contact as well as Particle Image Velocimetry (PIV) methods, also suffer from inability to resolve the dynamics of the near-surface air-flow, especially in the surface wave troughs. In this report, we present a first attempt to use Direct Numerical Simulation (DNS) as tool for investigation of the drops-mediated momentum, heat and moisture transfer in a turbulent, droplet-laden air flow over a wavy water surface. DNS is capable of resolving the details of the transfer processes and do not involve any closure assumptions typical of Large-Eddy and Reynolds Averaged Navier-Stokes (LES and RANS) simulations. Thus DNS provides a basis for improving parameterizations in LES and RANS closure models and further development of large-scale prognostic models. In particular, we discuss numerical results showing the details of the modification of the air flow velocity, temperature and relative humidity fields by multidisperse, evaporating drops. We use Eulerian-Lagrangian approach where the equations for the air-flow fields are solved in a Eulerian frame whereas the drops dymanics equations are solved in a Largangain frame. The effects of air flow and drops on the water surface wave are neglected. A point-force approximation is employed to model the feed-back contributions by the drops to the air momentum, heat and moisture transfer.

The “Bag Breakup” Spume Droplet Generation Mechanism at High Winds. Part II: Contribution to Momentum and Enthalpy Transfer

AbstractIn Part I of this study, we used high-speed video to identify “bag breakup” fragmentation as the dominant mechanism by which spume droplets are generated at gale-force and hurricane wind sp...

Numerical simulation of small-scale mixing processes in the upper ocean and atmospheric boundary layer

The processes of turbulent mixing and momentum and heat exchange occur in the upper ocean at depths up to several dozens of meters and in the atmospheric boundary layer within interval of millimeters to dozens of meters and can not be resolved by known large- scale climate models. Thus small-scale processes need to be parameterized with respect to large scale fields. This parameterization involves the so-called bulk coefficients which relate turbulent fluxes with large-scale fields gradients. The bulk coefficients are dependent on the properties of the small-scale mixing processes which are affected by the upper-ocean stratification and characteristics of surface and internal waves. These dependencies are not well understood at present and need to be clarified. We employ Direct Numerical Simulation (DNS) as a research tool which resolves all relevant flow scales and does not require closure assumptions typical of Large-Eddy and Reynolds Averaged Navier-Stokes simulations (LES and RANS). Thus DNS provides a solid ground for correct parameterization of small-scale mixing processes and also can be used for improving LES and RANS closure models. In particular, we discuss the problems of the interaction between small-scale turbulence and internal gravity waves propagating in the pycnocline in the upper ocean as well as the impact of surface waves on the properties of atmospheric boundary layer over wavy water surface.

Direct numerical simulation of a turbulent stably stratified air flow above a wavy water surface

The influence of the roughness of the underlaying water surface on turbulence is studied in a stably stratified boundary layer (SSBL). Direct numerical simulation (DNS) is conducted at various Reynolds (Re) and Richardson (Ri) numbers and the wave steepness ka. It is shown that, at constant Re, the stationary turbulent regime is set in at Ri below the threshold value Ric depending on Re. At Ri > Ric, in the absence of turbulent fluctuations near the wave water surface, three-dimensional quasiperiodical structures are identified and their threshold of origin depends on the steepness of the surface wave on the water surface. This regime is called a wave pumping regime. The formation of three-dimensional structures is explained by the development of parametric instability of the disturbances induced by the surface water in the air flow. The DNS results are quite consistent with prediction of the theoretical model of the SSBL flow, in which solutions for the disturbances of the fields of velocity and temperature in the wave pumping regime are found to be a solution of a two-dimensional linearized system with the heterogeneous boundary condition, which is caused by the presence of the surface wave. In addition to the turbulent fluctuations, the three-dimensional structures in the wave pumping regime provide for the transfer of impulse and heat, i.e., the increase in the roughness of the water–air boundary caused by the presence of waves intensifies the exchange in the SSBL.