Note on the flow of a stratified fluid over a stationary obstacle in a channel

Abstract

A new method is introduced to produce a uniform stratified flow over a stationary obstacle in an open channel. The flow is achieved by discharging the flow from the channel through a sink. The details of the sink are unimportant. The flow speed is limited only by the sink capacity. Selective withdrawal at lower densimetric Froude numbers is effectively eliminated through the use of a contraction. The standing, free-surface, long wave arising from the initiation of the flow is also eliminated by the contraction. Experiments are conducted for flow over a sphere for a range of Reynolds numbers from O(10²) to O(10³) and a range of Richardson numbers from O(10^?1) to O(10). Dye and neutrally buoyant droplets are used for quantitative analysis of the wake structure. The wake is also probed by a hot-film anemometer. The frequency of vortex shedding is obtained. Comparison with data from towed experiments is also presented.     

Internal waves in a randomly stratified fluid

Abstract

We discuss the propagation of internal waves in a rotating stratified unbounded fluid with randomly varying stability frequency, N. The first order smoothing approximation is used to derive the dispersion relation for the mean wave field when N is of the form N ² = N _o ²(1 + ?μ), where μ is a centered stationary random function of either depth (z) or time (t), N _o = constant and O < ?² ≦ 1. Expressions are then derived for the change in phase speed and growth rate due to the random fluctuations μ; in particular, attention is focused on the behaviour of these expressions for short and long correlation lengths (case μ = μ(z)) and times (case μ = μ(t)). For the case μ = μ(z), which represents a model for the temperature and salinity fine-structure in the ocean, the appropriate statistics of the fluctuations observed at station P (50°N, 145°W) have been incorporated into the theory to estimate the actual importance of the effects due to these random fluctuations. It is found that the phase speed of the mean wave decreases significantly if (i) the wavelength is short compared to g/N_o ² or (ii) the wave number vector is essentially horizontal and the wave frequency is very close to N _o. Also, the random fluctuations cause a significant growth (decay) in the amplitude of a wave propagating upwards (downwards) through a depth of a few kilometers. However, in the direction of energy propagation, the kinetic energy is conserved. Finally, it is shown that the average effect of the depth dependent fluctuations at station P is to slightly decrease the stability frequency and the magnitude of the group velocity.     

On the mathematical stability of stratified flow models with local turbulence closure schemes

Occasionally, numerical simulations using local turbulence closure schemes to estimate vertical turbulent fluxes exhibit small-scale oscillations in space, causing the eddy coefficients to vary over several orders of magnitude on short distances. Theoretical developments suggest that these spurious oscillations are essentially due to the way the eddy coefficients depend on the vertical gradient of the model’s variables. An instability criterion is derived based on the assumptions that the artefacts under study are due to the development of small-amplitude, small time- and space-scale perturbations of a smooth solution. The relevance of this criterion is demonstrated by applying it to a series a closure schemes, ranging from the Pacanowski–Philander formulas to the Mellor–Yamada level 2.5 model.     

Instabilities in a staircase stratified shear flow

We study stratified shear flow instability where the density profile takes the form of a staircase of interfaces separating uniform layers. Internal gravity waves riding on density interfaces can resonantly interact due to a background shear flow, resulting in the Taylor-Caulfield instability. The many steps of the density profile permit a multitude of interactions between different interfaces, and a rich variety of Taylor-Caulfield instabilities. We analyse the linear instability of a staircase with piecewise-constant density profile embedded in a background linear shear flow, locating all the unstable modes and identifying the strongest. The interaction between nearest-neighbour interfaces leads to the most unstable modes. The nonlinear dynamics of the instabilities are explored in the long-wavelength, weakly stratified limit (the defect approximation). Unstable modes on adjacent interfaces saturate by rolling up the intervening layer into a distinctive billow. These nonlinear structures coexist when stacked vertically and are bordered by the sharp density gradients that are the remnants of the steps of the original staircase. Horizontal averages remain layer-like.     

Turbulent mixing in a rotating,stratified fluid

Abstract

An experimental study was carried out to investigate the effect of rotation on turbulent mixing in a stratified fluid when the turbulence in the mixed layer is generated by an oscillating grid. Two types of experiments were carried out: one of them is concerned with the deepening of the upper mixed layer in a stable, two-fluid system, and the other deals with the interaction between a stabilizing buoyancy flux and turbulence.

In the first type of experiments, it was found that rotation suppresses entrainment at larger Rossby numbers. As the Rossby number becomes smaller (Ro 0.1), the entrainment rate increases with rotation—the onset of this phenomenon, however, was found to coincide with the appearance of coherent vortices within the mixed layer. The radiation of energy from the mixed layer to the lower non-turbulent layer was found to occur and the magnitude of the energy flux was found to be increased with the rotational frequency. It was also observed that vortices are generated, rather abruptly, in the lower layer as the mixed layer deepens.

In the second set of experiments a quasi-steady mixed layer was found to develop of which the thickness varies with rotation in a fashion that is consistent with the result of the first experiment. Also the rotation was found to delay the formation of a pycnocline.     

Space time evolution of sand bed topography and associated flow turbulence: experiments with statistical analysis

A series of flume experiments were conducted with varying the flow discharges at the Fluvial Mechanics Laboratory of Indian Statistical Institute (Kolkata) to understand the co-evolution patterns of generating bed forms and the corresponding flow turbulence. Instantaneous bed elevations and velocity components were recorded continuously for sufficient time using high resolution instruments, such as, Ultrasonic Ranging System and acoustic Doppler velocimeter, at some spatial location over the deformed bed for each flow discharge. Increase in mean bed elevations and bed-slopes was found to be increased in discharge. Heavy-tailed nature of the probability density functions of magnitude of bed elevation increments, magnitude of single continuous bed elevation increments and instantaneous Reynolds shear stresses along three planes were analyzed using Pareto and truncated Pareto distributions. The spectral analysis of bed elevations revealed that the slope of log–log linearity increased with increase in flow discharge. Wavelet cross-correlations depicted strong dependence of bed form evolution on the corresponding instantaneous Reynolds shear stress along xz-plane. A Gram–Charlier type of distribution was used to estimate the probability density function of fluctuating velocity components, instantaneous Reynolds shear stresses along three planes and the joint probability density functions of the fluctuating velocity components, which showed good fit with the experimental data.     

The statistics of a passive scalar in field-guided magnetohydrodynamic turbulence

A variety of studies of magnetised plasma turbulence invoke theories for the advection of a passive scalar by turbulent fluctuations. Examples include modelling the electron density fluctuations in the interstellar medium, understanding the chemical composition of galaxy clusters and the intergalactic medium, and testing the prevailing phenomenological theories of magnetohydrodynamic turbulence. While passive scalar turbulence has been extensively studied in the hydrodynamic case, its counterpart in MHD turbulence is significantly less well understood. Herein we conduct a series of high-resolution direct numerical simulations of incompressible, field-guided, MHD turbulence in order to establish the fundamental properties of passive scalar evolution. We study the scalar anisotropy, establish the scaling relation analogous to Yaglom’s law, and measure the intermittency of the passive scalar statistics. We also assess to what extent the pseudo Alfvén fluctuations in strong MHD turbulence can be modelled as a passive scalar. The results suggest that the dynamics of a passive scalar in MHD turbulence is considerably more complicated than in the hydrodynamic case.     

Numerical modelling of stratified flow: a spectral approach

A numerical formulation is developed to solve the three-dimensional hydrodynamic equations which describe flow in a stratified sea.Arbitrary continuous physically realistic variations of density and eddy viscosity can be included in the model, which is sufficiently flexible to be applicable to sea areas of any horizontal extent and depth. A continuous current profile from sea surface to sea bed, is computed with the model. A method for expanding computed current profiles in terms of vertical modes is proposed and the contribution of these modes to the current profiles is considered.The time variation of the wind-induced circulation of a stratified lake in response to a suddenly imposed and maintained wind stress is examined. Calculations show that the wind-driven surface current is modulated by the internal seiche motion of the lake.     

Flow of a double-diffusive stratified fluid in a differentially-rotating cylinder

Abstract

A study has been made of a basic state of axisymmetric flow, at large rotational Reynolds numbers, in a double-diffusive stratified fluid contained in a vertically-mounted, differentially-rotating cylindrical cavity. The aim is to describe the qualitative characteristics of the flow of a fluid, the density of which is stratified by two diffusive effects, i.e., temperature and salinity gradients. Attention is confined to situations in which the temperature and salinity gradients make opposing contributions to the overall density profile, the undisturbed stratification being gravitationally stable. Finite difference numerical solutions of the governing Navier-Stokes equations have been obtained using the Boussinesq approximation. The results are presented in a way that illustrates the explicit effects of double-diffusivity when the cavity aspect ratio, height/radius, is O(1). The principal non-dimensional parameters characterizing the flow field are identified. In the interior core, the primary dynamic balance is between the horizontal density gradient and the vertical shear of the prevailing azimuthal velocity. The effective stratification is seen to decrease as the double-diffusivity increases, even if the overall stratification parameter, St, is held constant. The solute field contains a very thin boundary layer structure at large Lewis numbers. The effective stratification increases with the Prandtl number. Results have been derived for extreme values of the cavity aspect ratio. For small cavity aspect ratios, the dominant dynamic ingredients are viscous diffusion and rotation. For large aspect ratios, the bulk of the flow field is determined by the rotating sidewall. In this case, the direct influence of the double-diffusivity is minor.     

Wind mixing of a stratified shear flow

The response of a shear flow to an imposed wind stress is studied both theoretically and by means of a numerical turbulence model. It is shown that for small initial gradient Richardson numbers (Ri₀ ≲ 4/3) a tail wind causes the slab velocity of the upper mixed layer to decrease. The theory is based on the assumption that during the wind-induced entrainment process the overall Richardson number will adjust to a quasi-constant value (Ri_u ≈ 2/3). The turbulence model is the so-called k-ɛ model. It is calibrated to five conditions by tuning only one constant. The details of the deepening process and the density and velocity distributions of the upper mixed layer during this anomalous behavior are thus made clear. The results imply that the common practice of estimating the total current velocity by vector addition of the original velocity and the wind-induced velocity (calculated from models based on an ocean at rest) may lead to an overestimation of the current speed.     

Upon boundary conditions imposed by a stratified fluid

Abstract

It is shown that in the limit of slow steady linearized flow and infinite stratification a stratified region lying above or below a layer of destabilized fluid produces zero velocity and zero stress boundary conditions. These novel boundary conditions constrain flow more fully than the more common rigid or free conditions, and arise because the penetrative flow in the stabilized region must pump only a finite amount of heat.     

Thermal and magnetically driven instabilities in a non-constantly stratified fluid layer

Abstract

The linear hydromagnetic stability of a non-constantly stratified horizontal fluid layer permeated by an azimuthal non-homogeneous magnetic field is investigated for various widths of the stably stratified part of the layer in the geophysical limit q→0 (q is the ratio of thermal and magnetic diffusivities). The choice of the strength of the magnetic field B_o is as in Soward (1979) (see also Soward and Skinner, 1988) and the equations for the disturbances are treated as in Fearn and Proctor (1983). It was found that convection is developed in the whole layer regardless of the width of its stably stratified part. The thermal instability penetrates essentially from the unstably stratified part of the layer into the stably stratified part for A ～ 1 (A characterises the ratio of the Lorentz and Coriolis forces). When the magnetic field is strong (A>1) the thermal convection is suppressed in the stably stratified part of the layer. However, in this case, it is replaced by the magnetically driven instability; which is fully developed in the whole layer. The thermal instabilities always propagate westward and exist for all the modes m. The magnetically driven instabilities propagate either westward or eastward according to the width of the stably and unstably stratified parts and exist only for the mode m=1.     

Forces on spheres moving horizontally in a rotating stratified fluid

Abstract

Measurements have been made of the net horizontal force F acting on a sphere moving with horizontal velocity U (Reynolds numbers in the range 10²-10⁴) through a stratified fluid rotating about a vertical axis with uniform angular velocity Ω. In both homogeneous and stratified rotating fluids with small Rossby number R(R = U/Ωa ? 1 where a is the radius of the sphere) the force F is of magnitude 2ΩρUV (where ρ is the density of the fluid and V is the volume of the sphere). In a homogeneous fluid the relative directions of F and U were found to depend on the quantity F = 8Ωa ²/UD (where D is the depth of the fluid in which the object is placed (Mason, 1975)). In a rotating stratified fluid the relative directions of F and U are found to depend on the inverse Froude number k(k = Na/U where N ² = (g/δ)?ρ/?z) provided D > 4aΩ/N. In a homogeneous fluid with F ? 1 the force F is mainly in the U direction (a drag force due to inertial wave radiation) and is ～ ?0.4 |MX 2ΩρUV For F ? 1 a “Taylor column” occurs and the force, in correspondence with theoretical expectations, is ～ - 2Ω |MX UρV In a rotating stratified fluid with N ～2Ω and k ? 1 the force F is mainly in the U direction but is roughly one half of that occurring in the homogeneous situation with F ? 1 (tentatively explained as due to the evanescence of inertia-gravity disturbances). In a rotating stratified fluid with k ? 1 the flow should have no vertical motion (as with F ? 1) and again in correspondence with theoretical expectations the drag is ～ ?2 Ω |MX UρV. In a non-rotating stratified fluid the drag coefficient C _D(C _D = F U/½?ρU ²) was measured in the range k = 0.1 to 10 and had a maximum value ～ 1.2 for k ～ 3.     

Closed form solution for a point force in a rotating stratified fluid

Summary The motion generated by an oscillatory point force in an inviscid incompressible rotating stratified fluid is studied, taking into account the effect of stratification on the inertia terms in the equations of motion. The solutions are obtained in closed form using Fourier transforms. The motion is axisymmetric (or not) according as the force acts along (or perpendicular to) the axis of rotation of the fluid. In the hyperbolic case, the motion is confined to the inner region of a cone with vertex at the origin. The solution when there is no rotation-no stratification-neither rotation nor stratification, is deduced and each case is found to differ from the case of rotating stratified fluid in many respects. The solution for a steady force is also deduced in the limit.     

Internal waves in a rotating stratified fluid in an arbitrary gravitational field

Abstract

Small amplitude oscillations of a uniformly rotating, density stratified, Boussinesq, non-dissipative fluid are examined. A mathematical model is constructed to describe timedependent motions which are small deviations from an initial state that is motionless with respect to the rotating frame of reference. The basic stable density distribution is allowed to be an arbitrary prescribed function of the gravitational potential. The problem is considered for a wide class of gravitational fields. General properties of the eigenvalues and eigenfunctions of square integrable oscillations are demonstrated, and a bound is obtained for the magnitude of the frequencies. The modal solutions are classified as to type. The eigenfunctions for the pressure field are shown to satisfy a second-order partial differential equation of mixed type, and the equation is obtained for the critical surfaces which delineate the elliptic and hyperbolic regions. The nature of the problem is examined in detail for certain specific gravitational fields, e.g., a radially symmetric field. Where appropriate, results are compared with those of other investigations of waves in a rotating fluid of spherical configuration and the novel aspects of the present treatment are emphasized. Explicit modal solutions are obtained in the specific example of a fluid contained in a rigid cylinder, stratified in the presence of vertical gravity, with the buoyancy frequency N being an arbitrary prescribed function of depth.     

Lagrangian approach to geostrophic adjustment of frontal anomalies in a stratified fluid

Geostrophic adjustment of frontal anomalies in a rotating continuously stratified fluid is studied in the standard framework of strictly rectilinear fronts and jets. Lagrangian approach to this problem is developed allowing to analyze, in a conceptually and technically simple way, both major problems of the nonlinear adjustment: the existence of a smooth adjusted state for a given set of initial conditions and the attainability of the adjusted state during the adjustment process. Dynamical splitting into balanced (adjusted state) and unbalanced (inertia-gravity waves) motions becomes transparent in the Lagrangian approach. Conditions of existence of the balanced state in the unbounded domain are established. It is shown that nonexistence of a smooth adjusted state in the vertically bounded domains is generic and a parallel with the classical scenario of deformation frontogenesis is developed. Small perturbations around smooth adjusted states are then studied with special emphasis on the wave-trapping inside the jet/front. Trapped modes with horizontal scales comparable to the width of the jet are explicitly constructed for a barotropic jet and their evolution is studied with the help of the WKB-approximation for weakly baroclinic jets. Modifications of the standard scenario of adjustment due to subinertial (quasi-) trapped modes and implications for data analysis are discussed.     

Effects of dissipation on internal waves in a contained rotating stratified fluid

Abstract

Boundary layer techniques are used to examine the modifications due to dissipation in the normal modes of a uniformly rotating, density stratified, Boussinesq fluid in a rigid container. Arbitrary relative influence of rotation and stratification is considered. The existence of critical regions of the container boundary is discussed. In cylindrical geometry a formula is derived for the decay factor on the homogeneous “spin-up” time scale which reveals how the dominant dissipation varies as a function of several parameters. For the situation where the buoyancy and inertial frequency are exactly equal, all boundaries are everywhere critical. In this case the method of multiple time-scales is employed to investigate the confluence inertial-gravity mode which is shown to persist until the diffusive time-scale is achieved.     

Critical-level crossing of rossby-gravity wave-packets in a linear-shear,stratified flow

Abstract

The propagation of Rossby-gravity waves, which are of astronomical interest near the critical level (or corotation point), is considered in a linear-shear, exponentially stratified flow by means of a ray tracing method. It is shown, using the analytic solutions for the stream function, that for a Richardson number J > 1/4 a wave-packet can cross the critical level in a finite time. In the unstable stratified case (J < 0), it cannot cross the critical level for large |J| (J < 0), but may do so for some intermediate |J| values.

Based on the above results, the possible existence of the regular normal mode with a discrete point eigenvalue in the continuous spectrum is discussed for bounded systems.     

Physical and numerical experiments on layered convection in a density-Stratified fluid

Abstract

The formation and growth of horizontal layered convection cells in a density stratified solution of salt water subject to an impulsively applied lateral temperature gradient is investigated with physical and numerical experiments. Results indicate that lyers are induced by two mechanisms. One is the successive formation of layers due to the presence of the top and bottom boundaries. The other is the spontaneous occurrence of layers when a suitably defined Rayleigh number exceeds a critical value. It is found that well established layers are homogeneous in temperature and salinity and are separated by sharp gradients in density. Lateral heat transfer is of a periodic nature. Numerical experiments were carried out for finite and infinite geometry cases. For the finite geometry case, convection cells are generated successively inward from the horizontal boundaries. For the infinite geometry case, periodic conditions in the vertical direction are assumed. With continuous input of small perturbations, simultaneous occurrence of the convection cells is obtained at supercritical Rayleigh numbers. Criteria for determining the onset of spontaneous cells numerically are explored.