Search

Stratified wakes past an isolated conical seamount are simulated at a Froude number of $Fr = 0.15$ and Rossby numbers of $Ro = 0.15$ , 0.75 and $\infty$ . The wakes exhibit a Kármán vortex street, unlike their unstratified, non-rotating counterpart. Vortex structures are studied in terms of large-scale global modes, as well as spatially localised vortex evolution, with a focus on rotation effects. The global modes are extracted by spectral proper orthogonal decomposition (SPOD). For all three studied $Ro$ ranging from mesoscale, submesoscale and non-rotating cases, the frequency of the SPOD modes at different heights remains coupled as a global constant. However, the shape of the SPOD modes changes from slanted ‘tongues’ at zero rotation ( $Ro=\infty$ ) to tall hill-height columns at strong rotation ( $Ro=0.15$ ). A novel method for vortex centre tracking shows that, in all three cases, the vortices at different heights advect uniformly at approximately $0.9U_{\infty }$ beyond the near wake, consistent with the lack of variability of the global modes. Under system rotation, cyclonic vortices and anticyclonic vortices (AVs) present considerable asymmetry, especially at $Ro = 0.75$ . The vorticity distribution as well as the stability of AVs are tracked downstream using statistics conditioned to the identified vortex centres. At $Ro=0.75$ , intense AVs with relative vorticity up to $\omega _z/f_{c}=-2.4$ (where $\omega_z$ is the vertical vorticity and $f_c$ is the Coriolis frequency) are seen with small regions of instability and they maintain large $\omega _z/f_{c}$ magnitude in the far wake. Recent stability analysis that accounts for stratification and viscosity is found to improve on earlier criteria and show that these intense AVs are stable.

We conducted experimental investigations on the effect of stable thermal conditions on rough-wall boundary layers, with a specific focus on their response to abrupt increases in surface roughness. For stably stratified boundary layers, a new analytical relation between the skin-friction coefficient, $C_f$ , and the displacement thickness was proposed. Following the sharp roughness change, the overshoot in $C_f$ is slightly enhanced in stably stratified layers when compared with that of neutral boundary layers. Regarding the velocity defect law, we found that the displacement thickness multiplied by $\sqrt{2/C_f}$ , performs better than the boundary layer thickness alone when describing the similarity within internal boundary layers for both neutral and stable cases. A non-adjusted region located just beneath the upper edge of the internal boundary layer was observed, with large magnitudes of skewness and kurtosis of streamwise and wall-normal velocity fluctuations for both neutral and stable cases. At a fixed wall-normal location, the greater the thermal stratification, the greater the magnitudes of skewness and kurtosis. Quadrant analysis revealed that the non-adjusted region is characterised by an enhancement/reduction of ejection/sweep events, particularly for stably stratified boundary layers. Spatially, these ejections correspond well with peaks of kurtosis, exhibit stronger intensity and occur more frequently following the abrupt change in surface conditions. [ABSTRACT FROM AUTHOR]

Air curtain devices (ACD) are commonly installed in domestic and commercial buildings to suppress the buoyancy-driven exchange flow through a doorway opening. Generally, the operating density of an ACD is equal to that of the indoor space making it neutrally buoyant. In the present study, we evaluate the performance of heavy air curtains where the operating density of the ACD is higher than that of the ambient fluid. The primary objective is to quantify the air curtain effectiveness, E, that determines the thermal comfort of building occupants based on the mean temperature inside the interrogated region. Experiments and numerical simulations are conducted and validated for various values of deflection modulus, D m , that compare the relative magnitude of the jet momentum and transverse stack effect due to buoyancy. The other important non-dimensional parameter is the density ratio, S, which compares the extent of added buoyancy in ACD to that of across the doorway. In addition, the velocity dynamics of the air curtains are compared with an isothermal jet to understand the underlying effects that the buoyancy causes on the jet development. The general structure of air curtains that characterize the jet inclination and penetration is visualized through injecting a dye, and it agrees very well with the buoyancy distribution obtained using simulations at different D m . Upon introduction of an assisting buoyancy, it has been found that the infiltration reduces by 25% compared to a neutrally buoyant air curtain for practical values of D m . Article Highlights: Heavier air curtains are more effective as compared to the neutrally buoyant air curtains. The deceleration of mean centreline velocity is more prominent at smaller deflection modulus. Stability may be achieved at slightly smaller deflection modulus for heavier air curtains. [ABSTRACT FROM AUTHOR]

Resistive tearing instabilities are common in fluids that are highly electrically conductive and carry strong currents. We determine the effect of stable stratification on the tearing instability under the Boussinesq approximation. Our results generalise previous work that considered only specific parameter regimes, and we show that the length scale of the fastest-growing mode depends non-monotonically on the stratification strength. We confirm our analytical results by solving the linearised equations numerically, and we discuss whether the instability could operate in the solar tachocline. [ABSTRACT FROM AUTHOR]

The linear stability of plane Couette–Poiseuille flow (CPF) is studied with the physical effects of stratification, rotation and viscosity all included for the first time together. With no stratification, two instability mechanisms are present due to the shear and rotation which, for the most part, do not interact as they favour different forms of two-dimensionality. However, there are some small parts of parameter space where new three-dimensional instability appears indicating that Rayleigh's criterion is also violated in parameter space beyond where shear instability is expected. No fully localised centrifugal instabilities can be found for CPF, but they are shown to exist if the base flow shear has a maximum in the domain (the base flow needs to be at least cubic in the cross-stream variable rather than just quadratic as in CPF). With stable stratification present, new instabilities are found due to the combined effects of stratification and rotation, but only some appear to be of the resonance-type associated with the strato-rotational instability. The other unstable branches are more accurately interpreted as a stratification-modified centrifugal instability. Three-dimensional versions of this violate Rayleigh's criterion even when this is extended to include stratification. When stratification is stronger than rotation, the resonance-type instabilities are only dominant for cyclonic flows. [ABSTRACT FROM AUTHOR]

A comprehensive and fine-grained transient analysis of the kinetic budget equation in the spectral and spatial domains for channel flow subjected to stable stratification is performed, through a series of large eddy simulations (LES) at a bulk Reynolds number of 2200 and for Froude numbers of 0.5, 2.0, 4.5 and 6.8, thereby covering a range of stratification strengths from high to low. The dynamic Smagorinsky eddy viscosity-based sub-grid scale (SGS) model is used for closure. Validation of the methodology is carried out by comparing the turbulent flow LES results to previously published channel flow results. The analysis is focussed on the effect of stratification on the energy cascade process of an initially turbulent flow in a channel configuration. The terms in the kinetic budget equation, in the homogeneous and inhomogeneous directions, are compared for different Froude numbers to assess the effects of weak and strong stratification on turbulent statistics. Transient results are plotted to highlight the role played by different terms in the budget equation and to understand the formation of the gravity waves in the core or outer regions of the wall-bounded flow in the presence of stratification. The energy cascade reveals some interesting observations about the functions of viscous diffusion, pressure diffusion, and inhomogeneous energy transfer in the initial phase of the development of flow, which in turn helps to explain the stationary waves observed in the presence of strong stratification. [ABSTRACT FROM AUTHOR]

Emptying-filling boxes have been studied in a wide range of configurations for decades, but the flow created in the box by two plumes rising from sources of arbitrary strength and elevation was previously unsolved. Guided by experiments and simplified analytical modelling, we reveal a rich array of two- and three-layer stratifications across seven possible flow regimes. The governing equations for these regimes show how the prevailing regime and stratification properties vary with three key parameters: the relative strength of the plumes, the height difference between their sources and a parameter characterising the resistance of the box to emptying. We observe and explain new behaviours not described in previous studies that are crucial to understanding emptying-filling boxes with multiple plumes. In particular, we demonstrate that the oft-assumed premise that n plumes leads to a stratification with n + 1 layers is not necessarily true, even in the absence of mixing. Two emptying-filling box models are developed: an analytical model addressing all combinations of the governing parameters and an extended model for three-layer stratifications that incorporates two mixing processes observed in the experiments. The predictions of these two models are in generally excellent agreement with measurements from the experimental campaign covering 69 combinations of the governing parameters. This study improves our understanding of emptying-filling boxes and could facilitate improvements to natural ventilation building design, as demonstrated by an example scenario in which occupants feel cooler upon the addition of a second source of heat. [ABSTRACT FROM AUTHOR]

The linear stability of the Stuart vortices, which is a model of arrays of vortices often observed in the atmosphere and the oceans, in rotating stratified fluids is investigated by local and modal stability analysis. As in the case of the two-dimensional (2-D) Taylor-Green vortices, five types of instability appear in general: the pure-hyperbolic instability, the strato-hyperbolic instability, the rotational-hyperbolic instability, the centrifugal instability and the elliptic instability. The condition for each instability and the estimate of the growth rate derived by Hattori & Hirota (J. Fluid Mech., vol. 967, 2023, A32) are shown to also be useful for the Stuart vortices, which supports their applicability to general flows. The properties of each instability depend on stratification and rotation in a way similar to the case of the 2-D Taylor-Green vortices. For the Stuart vortices, however, the centrifugal instability and the elliptic instability become more dominant than the three hyperbolic instabilities in comparison to the 2-D Taylor-Green vortices; this is explained by the larger ratios of the maximum vorticity and the strain rate at the elliptic stagnation points to the strain rate at the hyperbolic stagnation points. Direct correspondence between the modal and local stability results is further established by comparing unstable modes to solutions to the local stability equations; this is useful for identifying the types of modes since the mechanism of instability is readily known in the local stability analysis. This helps us to discover the modes of the ring-type elliptic instability, which have been predicted only theoretically. [ABSTRACT FROM AUTHOR]

Wave turbulence provides a conceptual framework for weakly nonlinear interactions in dispersive media. Dating from five decades ago, applications of wave turbulence theory to oceanic internal waves assigned a leading-order role to interactions characterized by a near equivalence between the group velocity of high-frequency internal waves with the phase velocity of near-inertial waves. This scale-separated interaction leads to a Fokker-Planck (generalized diffusion) equation. More recently, starting four decades ago, this scale-separated paradigm has been investigated using ray tracing methods. These ray methods characterize spectral transport of energy by counting the amplitude and net velocity of wave packets in phase space past a high-wavenumber gate prior to 'breaking'. This explicitly advective characterization is based on an intuitive assignment and lacks theoretical underpinning. When one takes an estimate of the net spectral drift from the wave turbulence derivation and makes the corresponding assessment, one obtains a prediction of spectral transport that is an order of magnitude larger than either observations or reported ray tracing estimates. Motivated by this contradiction, we report two parallel derivations for transport equations describing the refraction of high-frequency internal waves in a sea of random inertial waves. The first uses standard wave turbulence techniques and the second is an ensemble-averaged packet transport equation characterized by the dispersion of wave packets about a mean drift in the spectral domain. The ensemble-averaged transport equation for ray tracing differs in that it contains the intuitively motivated advective term. We conclude that the aforementioned contradiction between theory, numerics and observations needs to be taken at face value and present a pathway for resolving this contradiction. [ABSTRACT FROM AUTHOR]

The concept of added mass is generalized to stratified fluids, accounting for the presence of internal waves. Once the added mass of a moving body is known, so is the hydrodynamic force exerted on it by the fluid, and the energy imparted by it to the fluid. As a function of frequency, added mass is complex. Its real part is associated with inertia and its imaginary part, only present in the frequency range of propagating waves, with wave damping. Owing to causality, these two parts satisfy Kramers-Kronig relations. The added masses of an elliptic cylinder of horizontal axis, typical of two-dimensional bodies, and a spheroid of vertical axis, typical of three-dimensional bodies, are deduced from their dipole strengths, themselves deduced from their representations as single layers. The wave power is shown to be a maximum, for fixed oscillation amplitude, at approximately 0.8 times the buoyancy frequency. In the temporal domain, added mass appears as a new memory force taking the form of a convolution integral. The kernel of this integral combines algebraically decaying oscillations at the buoyancy frequency on the one hand; and an exponentially damped oscillation for the horizontal motion of the spheroid, implying short-termmemory, an aperiodic algebraic decay for its vertical motion, implying long-term memory, and a constant for the motion of the cylinder, implying everlasting memory, on the other hand. A limitation of the study is its restriction to translational motion. [ABSTRACT FROM AUTHOR]

We conducted experimental investigations on the effect of stable thermal conditions on rough-wall boundary layers, with a specific focus on their response to abrupt increases in surface roughness. For stably stratified boundary layers, a new analytical relation between the skin-friction coefficient, Cf, and the displacement thickness was proposed. Following the sharp roughness change, the overshoot in Cf is slightly enhanced in stably stratified layers when compared with that of neutral boundary layers. Regarding the velocity defect law, we found that the displacement thickness multiplied by x 2/Cf, performs better than the boundary layer thickness alone when describing the similarity within internal boundary layers for both neutral and stable cases. A non-adjusted region located just beneath the upper edge of the internal boundary layer was observed, with large magnitudes of skewness and kurtosis of streamwise and wall-normal velocity fluctuations for both neutral and stable cases. At a fixed wall-normal location, the greater the thermal stratification, the greater the magnitudes of skewness and kurtosis. Quadrant analysis revealed that the non-adjusted region is characterised by an enhancement/reduction of ejection/sweep events, particularly for stably stratified boundary layers. Spatially, these ejections correspond well with peaks of kurtosis, exhibit stronger intensity and occur more frequently following the abrupt change in surface conditions. [ABSTRACT FROM AUTHOR]

Through a combination of laboratory experiments and theoretical models, we investigate the interaction of a mean upwelling through a closed basin with a vertical buoyancy flux. The fluid is mixed by a horizontally oscillating rake, which either traverses the whole basin or which oscillates just near one vertical boundary. We first review the steady state and demonstrate that, in both mixing regimes, the vertical density profile across the basin is controlled by the steady-state balance between the upward advective and diffusive fluxes of salinity as described by the classical model introduced by Munk (Deep-Sea Res. , vol. 13, issue 4, 1966, pp. 707–730). However, with boundary mixing, we show that both the upwelling and the buoyancy transport are localised to the mixing zone near the boundary, and the interior fluid is stagnant. We then develop a model to describe the transient evolution of the system if there is either a discrete increase or gradual decrease to the buoyancy flux. In the boundary mixing case, the change in the buoyancy flux at the lower boundary leads to a change in the buoyancy of the fluid in the boundary mixing region, and this induces a transient, buoyancy-driven flow in the boundary region in addition to the steady upwelling. In turn, an equal and opposite vertical flow develops in the interior, and this leads to a change in the density stratification of the interior fluid as the system adjusts to a new equilibrium. However, in our experiments, there is no vertical mixing in the interior and interior fluid may upwell or downwell dependent on the change to the buoyancy forcing. We discuss the implications of our results for the transport and mixing in the deep ocean, and the associated interpretation of field experiments. [ABSTRACT FROM AUTHOR]

This paper addresses the stability of plane Couette flow in the presence of strong density and viscosity stratifications. It demonstrates the existence of a generalised inflection point that satisfies the generalised Fjørtoft criterion of instability when a minimum of kinematic viscosity is present in the base flow. The characteristic scales associated with this minimum are identified as the primary controlling parameters of the associated instability, regardless of the type of stratification. To support this finding, analytical stability models are derived in the long-wave approximation using piecewise linear base flows. Numerical stability calculations are carried out to validate these models and to provide further information on the production of disturbance vorticity. All instabilities are interpreted as arising from the interaction between two vorticity waves. Depending on the type of stratification, these two waves are produced by different physical mechanisms. When both strong density and viscosity stratifications are present, we show that they result from the concurrent action of shear and inertial baroclinic effects. The stability models developed for simple fluid models ultimately shed light on a recently observed unstable mode in supercritical fluids (Ren et al. , J. Fluid Mech. , vol. 871, 2019, pp. 831–864), providing a quantitative prediction of the stability diagram and identifying the dominant mechanisms at play. Furthermore, our study suggests that the minimum of kinematic viscosity reached at the Widom line in these fluids is the leading cause of their instability. The existence of similar instabilities in different fluids and flows (e.g. miscible fluids) is finally discussed. [ABSTRACT FROM AUTHOR]

The oscillations of buoyant bodies in stratified fluids are deduced from the variations of their added mass. Three configurations are considered: a body displaced from its neutral level then released; a Cartesian diver set into oscillation by a modulation of the hydrostatic pressure, then released; and a body attached to a pendulum to which an impulse is applied. The first configuration is related to the dynamics of Lagrangian floats in the ocean. Two particular bodies are considered: an elliptic cylinder of horizontal axis, typical of two-dimensional bodies; and a spheroid of vertical axis, typical of three-dimensional bodies. The ultimate motion of the body consists of oscillations at the buoyancy frequency with an amplitude decaying algebraically with time. Before that, the resonant response of the system is observed, either aperiodic exponential decay when the system has no intrinsic dynamics, or exponentially damped oscillation otherwise. Comparison with available measurements demonstrates the need to include viscous dissipation in the analysis. At high Stokes number, dissipation comes from the Basset–Boussinesq memory force and is affected negligibly by the stratification; at low Stokes number, dissipation comes from Stokes resistance and exhibits a significant effect of the stratification. [ABSTRACT FROM AUTHOR]

Detailed observations of the velocities of Jovian vortices exist at only one height in the atmosphere, so their vertical structures are poorly understood. This motivates this study that computes stable three-dimensional, long-lived planetary vortices that satisfy the equations of motion. We solve the anelastic equations with a high-resolution pseudo-spectral method using the observed Jovian atmospheric temperatures and zonal flow. We examine several families of vortices and find that constant-vorticity vortices, which have nearly uniform vorticity as a function of height, and horizontal areas that go to zero at their tops and bottoms, converge to stable vortices that look like the Great Red Spot (GRS) and other Jovian anticyclones. In contrast, the constant-area vortices proposed in previous studies, which have nearly uniform areas as a function of height, and vertical vorticities that go to zero at their tops and bottoms, are far from equilibrium, break apart, and converge to constant-vorticity vortices. Our late-time vortices show unexpected properties. Vortices that are initially non-hollow become hollow (i.e. have local minima of vertical vorticity at their centres), which is a feature of the GRS that cannot be explained with two-dimensional simulations. The central axes of the final vortices align with the planetary spin axis even if they align initially with the local direction of gravity. We present scaling laws for how vortex properties change with the Rossby number and other non-dimensional parameters. We prove analytically that the horizontal mid-plane of a stable vortex must lie at a height above the top of the convective zone. [ABSTRACT FROM AUTHOR]

We develop a physics-informed neural network (PINN) to significantly augment state-of-the-art experimental data of stratified flows. A fully connected deep neural network is trained using time-resolved experimental data in a salt-stratified inclined duct experiment, consisting of three-component velocity fields and density fields measured simultaneously in three dimensions at Reynolds number = O(103) and at Prandtl or Schmidt number = 700. The PINN enforces incompressibility, the governing equations for momentum and buoyancy, and the boundary conditions at the duct walls. These physics-constrained, augmented data are output at an increased spatio-temporal resolution and demonstrate five key results: (i) the elimination of measurement noise; (ii) the correction of distortion caused by the scanning measurement technique; (iii) the identification of weak but dynamically important three-dimensional vortices of Holmboe waves; (iv) the revision of turbulent energy budgets and mixing efficiency; and (v) the prediction of the latent pressure field and its role in the observed asymmetric Holmboe wave dynamics. These results mark a significant step forward in furthering the reach of experiments, especially in the context of stratified turbulence, where accurately computing three-dimensional gradients and resolving small scales remain enduring challenges. [ABSTRACT FROM AUTHOR]

The stratified inclined duct (SID) sustains an exchange flow in a long, gently sloping duct as a model for continuously forced density-stratified flows such as those found in estuaries. Experiments have shown that the emergence of interfacial waves and their transition to turbulence as the tilt angle is increased appears to be linked to a threshold in the exchange flow rate given by inviscid two-layer hydraulics. We uncover these hydraulic mechanisms by (i) using recent direct numerical simulations (DNS) providing full flow data in the key flow regimes (Zhu et al., J. Fluid Mech., vol. 969, 2023, A20), (ii) averaging these DNS into two layers, and (iii) using an inviscid two-layer shallow-water and instability theory to diagnose interfacial wave behaviour and provide physical insight. The laminar flow is subcritical and stable throughout the duct and hydraulically controlled at the ends of the duct. As the tilt is increased, the flow becomes supercritical everywhere and unstable to long waves. An internal jump featuring stationary waves first appears near the centre of the duct, then leads to larger-amplitude travelling waves, and to stronger jumps, wave breaking and intermittent turbulence at the largest tilt angle. Long waves described by the (nonlinear) shallow-water equation are interpreted locally as linear waves on a two-layer parallel base flow described by the Taylor--Goldstein equation. This link helps us to interpret long-wave instability and contrast it with short-wave (e.g. Kelvin--Helmholtz) instability. Our results suggest a transition to turbulence in SID through long-wave instability relying on vertical confinement by the top and bottom walls. [ABSTRACT FROM AUTHOR]

We introduce a geometric analysis of turbulent mixing in density-stratified flows based on the alignment of the density gradient in two orthogonal bases that are locally constructed from the velocity gradient tensor. The first basis connects diapycnal mixing to rotation and shearing motions, building on the recent 'rortex--shear decomposition' in stratified shear layers (Jiang et al., J. Fluid Mech., vol. 947, 2022, A30), while the second basis connects mixing to the principal axes of the viscous dissipation tensor. Applying this framework to datasets taken in the stratified inclined duct laboratory experiment reveals that density gradients in locations of high shear tend to align preferentially (i) along the direction of minimum dissipation and (ii) normal to the plane spanned by the rortex and shear vectors. The analysis of the local alignment across increasingly turbulent flows offers new insights into the intricate relationship between the density gradient and dissipation, and thus diapycnal mixing. [ABSTRACT FROM AUTHOR]

Diapycnal mixing plays a key role in the thermohaline circulation of the deep ocean. Field observations have suggested that this mixing is intensified over the rough topography along the boundaries of the ocean. In this study, we experimentally explore the transport of salinity and tracer in a horizontally stirred, stratified fluid with a steady vertical buoyancy flux. The mechanical mixing occurs either uniformly across the tank or near a sidewall. To compare uniform and boundary mixing, first we explore the steady state dynamics and find that, in both cases, the surfaces of constant density are horizontal. With uniform mixing, vertical transport of buoyancy occurs uniformly throughout the tank while, with boundary mixing, transport is confined to the turbulent region and the interior space remains quiescent and does not mix or have any net movement. We then explore the transient evolution of the stratification when the source of buoyancy is removed from the base of the system. In the boundary mixing case, the resulting divergence of the turbulent diffusive flux in the boundary region leads to a reduction in the buoyancy of the fluid in the boundary region, and a net upflow develops in the boundary region. In turn, this drives a downwelling in the interior. Vertical gradients in the rate of downwelling lead to stretching of isopycnals and, together, these processes enable the interior stratification to evolve. The experiments highlight that, with boundary mixing, the main transport of buoyancy occurs near the boundaries even though the interior is stratified; independent measurements of tracer mixing in the interior show that this fluid may be quiescent and stratified even though there is a large flux being transported in the boundary region. This has important implications for the interpretation of mixing data in the ocean. In a companion paper, (Li and Woods Part 2), we explore the interaction of a net upwelling with such boundary mixing. [ABSTRACT FROM AUTHOR]

Stratified flows are commonly observed two-phase flow regime in nuclear and chemical reactors. In the present work, stratified and wavy-stratified flows occurring in circular tubes of size 1.8–4.5 mm are studied experimentally and numerically. Free surface numerical modeling based on volume of fluids approach is used to study the changes occurring in stratified flow pattern, and the model is validated with adiabatic experimental results. The transition boundaries of smooth to wavy patterns are presented in the form of a map. The liquid film height distribution in spatial and time domain is also presented. The effect of constant heat flux on the flow patterns and liquid film thickness is detailed. It is found that, above a certain heat flux (33,000 W/m2), wavy appearance becomes predominant. A correlation is developed to determine the maximum and average liquid film thickness for various heat loads. [ABSTRACT FROM AUTHOR]

Direct numerical simulations of stratified open channel flows subject to a varying surface heat flux are performed. The influence of the diurnal heating time on the spatial and temporal variation of mixing in the flow and the characteristics of the mean flow state are examined. The control parameters are the bulk stability parameter λ B , defined through the ratio of the channel height δ and a bulk Obukhov length scale L B , and the diurnal time scale t ^ , defined as the ratio of the heating time to an eddy turnover time. The Prandtl number Pr and Reynolds number R e τ have values of 1 and 400. Simulations are performed over t ^ = 1 to 24 and λ B = 0.6 to 26. Two key flow features are used to classify the flow regimes observed, namely the laminar layer depth (LLD) and stratified layer depth (SLD) where the LLD is defined as the depth from the free surface when the buoyancy Reynolds number R e B ≈ 7 and the SLD is the depth from the free surface when the turbulent Froude number F r ≈ 1 . This study attempts to characterise how these length scales vary across the diel cycle. The LLD is a viscous length scale and a regime map of a viscous parameter, the bulk Obhukov Reynolds number R e L , and t ^ is presented to classify the LLD behaviour. A regime map of λ B and t ^ is presented to classify the behaviour of the SLD. Three classifications for each layer depth behaviour within a diel cycle form the basis of the regime maps for this paper: a neutral flow where the LLD or SLD does not exist (denoted by NL and NS), a stratified flow where the LLD or SLD are diurnally varying (denoted as DL and DS) and a persistent layer of the LLD or SLD (denoted as PL and PS). The transition between the NL to DL is t ^ ∝ R e L 4.5 , DL to PL is t ^ ∝ R e L - 0.5 , NS to DS is t ^ ∝ λ B 0 and DS to PS is t ^ ∝ λ B 1 . The regime maps may be used as a predictive tool to determine when suppressed mixing regimes occur in rivers. At each flow depth, the flow sweeps though a range of mixing states across the diel cycle. The local mixing efficiency are briefly assessed and found to scale well with the instantaneous Fr number according to the regimes proposed by Garanaik and Venayagamoorthy (J. Fluid Mech., vol. 867, 2019, pp. 323-333). Article Highlights: This paper reports on direct numerical simulations of stratified open channel flows subject to a varying surface heat flux. The results have found that: Increasing the diurnal time scale allows the flow to sweep through a wider range of flow states from turbulent to strongly stratified, Simulation data of the mixing efficiency and turbulent Froude number from this temporally varying and spatially inhomogenous flow that undergoes strong temporal forcing collapses well onto the parameterisation scheme of Garanaik & Venayagamoorthy (J. Fluid Mech., vol. 867, 2019, pp. 323–333) found for homogenous stratified flows, and Three distinct classifications (a persistent layer, a diurnal layer and one where the layer does not exist) of the laminar layer depth (LLD) and stratified layer depth (SLD) behaviour persists throughout a diel cycle and form a regime map given a λ B , R e τ and t ^ value. The relationship between each transitions are: from no LLD (NL) to a diurnal LLD (DL) t ^ ∝ R e L 4.5 , DL to a persistent LLD (PL) t ^ ∝ R e L - 0.5 , no SLD (NS) to a diurnal SLD (DS) t ^ ∝ λ B 0 and DS to a persistent SLD (PS) is t ^ ∝ λ B 1 . The regime maps may used as a predictive model to calculate when suppressed mixing transpires in rivers. [ABSTRACT FROM AUTHOR]

The effect of the initial condition upon the transport dynamics of miscible flowing fluids in a porous medium is investigated under viscosity and density contrasts. Such flows have attracted significant attention due to their importance in many fields of science and engineering, such as $\mathrm {CO}_2$ sequestration and aquifer remediation. Using high-resolution two-dimensional numerical simulations, we illustrate the impact of viscosity and density contrasts on the temporal evolution of the spreading and mixing quantities. We show that such impact depends on the initial shape of the source distribution where the solute is injected and on the intensity of the horizontal background flux. We find that rates of mixing are dependent on whether the solute is more or less viscous than the ambient fluid, a result usually not taken into consideration in studies on gravity fingering. At higher background flux, the effects due to horizontal viscous fingering dominate over gravitational fingering. Our computational analysis also suggests a non-trivial relationship between mixing and the length of the plume's interface under fingering instabilities. Finally, we show how a stratified permeability field can interact with these sources of instabilities and affect the transport behaviour of the plume. [ABSTRACT FROM AUTHOR]

Surface differential heating on a stably stratified fluid body drives an overturning circulation confined to the upper fluid region – here coined stratified horizontal convection (SHC). In this manuscript, we investigate the dynamics of SHC via laboratory experiments, exploring local and global flow properties. By considering the available potential energy of the system, we derive a unique length scale of SHC and introduce the Péclet number $Pe$ that captures both the stabilising effect of stratification and the destabilising effect of the baroclinic adjustment. We found that $Pe$ characterises local and global flow properties, including the fluid transport of the overturning circulation, the available mechanical energy and the flow dimensionality. Our study provides insights into the fluid dynamics of stratified environments that experience horizontal convection, such as lakes, oceans and atmospheres. [ABSTRACT FROM AUTHOR]

The stratified inclined duct (SID) experiment consists of a zero-net-volume exchange flow in a long tilted rectangular duct, which allows the study of realistic stratified shear flows with sustained internal forcing. We present the first three-dimensional direct numerical simulations (DNS) of SID to explore the transitions between increasingly turbulent flow regimes first described by Meyer & Linden (J. Fluid Mech. , vol. 753, 2014, pp. 242–253). We develop a numerical set-up that faithfully reproduces the experiments and sustains the flow for arbitrarily long times at minimal computational cost. We recover the four qualitative flow regimes found experimentally in the same regions of parameter space: laminar flow, waves, intermittent turbulence and fully developed turbulence. We find good qualitative and quantitative agreement between DNS and experiments and highlight the added value of DNS to complement experimental diagnostics and increase our understanding of the transition to turbulence, both temporally (laminar/turbulent cycles) and parametrically (as the tilt angle of the duct and the Reynolds number are increased). These results demonstrate that numerical studies of SID – and deeper integration between simulations and experiments – have the potential to lead to a better understanding of stratified turbulence. [ABSTRACT FROM AUTHOR]

Understanding the dynamics of staircases in salt fingering convection presents a long-standing theoretical challenge to fluid dynamicists. Although there has been significant progress, particularly through numerical simulations, there are a number of conflicting theoretical explanations as to the driving mechanism underlying staircase formation. The Phillips effect proposes that layering in stirred stratified flow is due to an antidiffusive process, and it has been suggested that this mechanism may also be responsible for salt fingering staircases. However, the details of this process, as well as mathematical models to predict the evolution and merger dynamics of staircases, have yet to be established. We generalise the theory of the Phillips effect to a three-component system (e.g. temperature, salinity, energy) and demonstrate a regularised nonlinear model of layering based on mixing length parametrisations. The model predicts both the inception of layering and its long-term evolution through mergers, while generalising, and remaining consistent with, previous results for double-diffusive layering based on flux ratios. Our model of salt fingering is formulated using spatial averaging processes, and closed by a mixing length parametrised in terms of the kinetic energy and the ratio of the temperature and salt gradients. The model predicts a layering instability for a bounded range of parameter values in the salt fingering regime. Nonlinear solutions show that an initially unstable linear buoyancy gradient develops into layers, which proceed to merge through a process of stronger interfaces growing at the expense of weaker ones. Our results indicate that these mergers are responsible for the characteristic increase of buoyancy flux through thermohaline staircases. [ABSTRACT FROM AUTHOR]

The linear stability of the two-dimensional Taylor–Green vortices, which is a spatially periodic array of vortices, in rotating stratified fluids is investigated by local and modal stability analysis. Five types of instability appear in general: the pure hyperbolic instability, the strato-hyperbolic instability, the rotational-hyperbolic instability, the centrifugal instability and the elliptic instability. The condition for each instability and the estimate of the growth rate, which are useful in interpreting numerical results, are obtained in the framework of local stability analysis. Realizability of an instability is introduced to predict whether an unstable mode corresponding to an unstable region found in the local stability analysis exists at finite Reynolds numbers. In the absence of stratification, the pure hyperbolic instability is dominant for weak rotation; it is stabilized for strong rotation. For strong anti-cyclonic rotation, the elliptic instability or the centrifugal instability becomes dominant depending on the parameter values; further stronger rotation stabilizes both instabilities. For strong cyclonic rotation, the rotational-hyperbolic instability or the elliptic instability becomes dominant, although the growth rate is smaller than the anti-cyclonic cases. Strong stratification changes the stability properties. The strato-hyperbolic instability occurs for weak rotation. The rotational-hyperbolic instability and the elliptic instability are weakened under cyclonic rotation, while the latter survives and extends the unstable range under anti-cyclonic rotation. The pure hyperbolic instability and the centrifugal instability are less affected by stratification. The mode structures of each instability are in good agreement with the corresponding solution to local stability equations, confirming the physical mechanism of the instability. [ABSTRACT FROM AUTHOR]

When an immiscible oil drop is immersed in a stably stratified ethanol-water mixture, the Marangoni flow on the surface of the drop can experience an oscillatory instability, so that the drop undergoes a transition from levitating to bouncing. The onset of the instability and its mechanisms have been studied previously (Li et al., Phys. Rev. Lett., vol. 126, 2021, 124502; Li et al., J. Fluid Mech., vol. 932, 2022, A11), yet the bouncing motion of the drop itself, which is a completely different problem, has not yet been investigated. Here we study how the bouncing characteristics (jumping height, rising and sinking time) depend on the control parameters (drop radius, stratification strength, drop viscosity). We first record experimentally the bouncing trajectories of drops of different viscosities in different stratifications. Then a simplified dynamical analysis is performed to get the scaling relations of the jumping height and the rising and sinking times. The rising and sinking time scales are found to depend on the drag coefficient CDS of the drop in the stratified liquid, which is determined empirically for the current parameter space (Zhang et al., J. Fluid Mech., vol. 875, 2019, 622-656). For low-viscosity (5 cSt) oil drops, the results on the drag coefficient match those from the literature (Yick et al., J. Fluid Mech., vol. 632, 2009, pp. 49-68; Candelier et al., J. Fluid Mech., vol. 749, 2014, pp. 184-200). For high-viscosity (100 cSt) oil drops, the parameter space had not been explored and the drag coefficients are not readily available. Numerical simulations are therefore performed to provide external verification for the drag coefficients, which well match with the experimental results. [ABSTRACT FROM AUTHOR]

In this study, we model the evolution of near-inertial waves (NIWs) under the β -effect by imposing initial zonal velocity in the mixed layer. We consider various background stratification values and initial velocity magnitudes, to investigate their effect on the wave characteristics and the decay rate of mixed layer kinetic energy. Increasing the interior stratification strength led to increased energy content in the higher vertical modes and also resulted in a faster decay of mixed layer kinetic energy. Double-inertial waves are observed in almost all the simulations, and the energy content in these waves was found to increase with an increase in the interior stratification strength as well as with the initial velocity amplitude. The mixed layer energy decay rate obtained from our simulations agrees well with the theoretical estimate by (Moehlis and Llewellyn Smith , 2001); however, we found a significant disagreement at higher velocity amplitudes. A strong shear zone is observed at the base of the mixed layer and the spatial extent of this zone increases with an increase in the strength of the stratification, initial velocity amplitude, and with the inclusion of the pycnocline. This study also highlights the importance of accurate parameterization of NIWs in the global models, as we found that the decay rate of mixed-layer kinetic energy is very sensitive to the vertical eddy viscosity parametrization in the model. [ABSTRACT FROM AUTHOR]

Abstract:

The reflection and transmission of internal waves at a horizontal interface between layers of constant buoyancy frequency are treated. This interface is a simple model of the atmospheric tropopause. The waves are weakly nonlinear and obey the nonlinear Schrödinger equation away from the interface. The waves are horizontally periodic and vertically confined to a long packet. The interfacial conditions are formally taken to third order, and include higher-order linear as well as nonlinear effects. The higher-order linear affects show an instability at the interface for steep waves. Numerical results for shorter wave packets show that the higher-order linear terms result in non-physical wave generation, and are ultimately neglected for such cases. Nonlinear effects result in a decrease in the reflected wave amplitudes and an increase in the transmitted wave amplitudes when compared to linear theory. Overall, the nonlinear effects make the interface more transparent to upwardly propagating internal waves. [ABSTRACT FROM AUTHOR]

In this paper we study the self-induced low-Reynolds-number flow generated by a cylinder immersed in a stratified fluid. In the low Péclet limit, where the Péclet number is the ratio of the radius of the cylinder and the Phillips length scale, the flow is captured by a set of linear equations obtained by linearising the governing equations with respect to the prescribed far field conditions. We specifically focus on the low Péclet regime and develop a Green's function approach to solve the linearised equations governing the flow over the cylinder. We cross check our analytical solution against numerical solution of the nonlinear equations to obtain the range of the Péclet numbers for which the linear solution is valid. We then take advantage of the analytical solution to find explicit far-field decay rates of the flow. Our detailed analysis points out that the streamfunction and the velocity field decays algebraically in the far field. Intriguingly, this algebraic decay of the flow is much slower when compared with the exponential decay of the flow generated by a slow moving cylinder in the homogeneous Stokes regime, in the absence of stratification. Consequently, the flow generated by a cylinder in the stratified Stokes regime will have a larger domain of influence when compared with the flow generated by a cylinder in the homogeneous Stokes regime. [ABSTRACT FROM AUTHOR]

The strongly nonlinear Miyata–Choi–Camassa model under the rigid lid approximation (MCC-RL model) can describe accurately the dynamics of large-amplitude internal waves in a two-layer fluid system for shallow configurations. In this paper, we apply the MCC-RL model to study the internal waves generated by a moving body on the bottom. For the case of the moving body speed $U=1.1c_{0}$ , where ${c_0}$ is the linear long-wave speed, the accuracy of the MCC-RL results is assessed by comparing with Euler's solutions, and very good agreement is observed. It is found that when the moving body speed increases from $U=0.8c_{0}$ to $U=1.241c_{0}$ , the amplitudes of the generated internal solitary waves in front of the moving body become larger. However, a critical moving body speed is found between $U=1.241c_{0}$ and $U=1.242c_{0}$. After exceeding this critical speed, only one internal wave right above the body is generated. When the moving body speed increases from $U=1.242c_{0}$ to $U=1.5c_{0}$ , the amplitudes of the internal waves become smaller. [ABSTRACT FROM AUTHOR]

We consider a special type of hydraulic jumps (internal bores) which, in the vertically bounded system of two immiscible fluids with slightly different densities, conserve not only the mass and impulse but also the circulation and energy. This is possible only at specific combinations of the upstream and downstream states. Two such combinations are identified with arbitrary upstream and downstream interface heights. The first has a cross-symmetry between the interface height and shear on both sides of the jump. This symmetry, which is due to the invariance of the two-layer shallow-water system with swapping the interface height and shear, ensures the automatic conservation of the impulse and energy as well as the continuity of characteristic velocities across the jump. The speed at which such jumps propagate is defined uniquely by the conservation of the mass and circulation. The other possibility is a marginally stable shear flow that can have fully conservative jumps with discontinuous characteristic velocities. Both types of conservative jumps are shown to represent a long-wave approximation to the so-called solibores that appear as smooth permanent-shape solutions in a weakly non-hydrostatic model. A new analytical solution for solibores is obtained and found to agree very well with the previous DNS results for partial-depth lock release flow. The finding that certain large-amplitude hydraulic jumps can be fully conservative, while most are not such even in the inviscid approximation, points towards the wave dispersion as a primary mechanism behind the lossy nature of internal bores. [ABSTRACT FROM AUTHOR]

Understanding how turbulence enhances irreversible scalar mixing in density-stratified fluids is a central problem in geophysical fluid dynamics. While isotropic overturning regions are commonly the focus of mixing analyses, we here investigate whether significant mixing may arise in anisotropic statically stable regions of the flow. Focusing on a single forced direct numerical simulation of stratified turbulence, we analyse spatial correlations between the vertical density gradient $\partial \rho /\partial z$ and the dissipation rates of kinetic energy $\epsilon$ and scalar variance $\chi$ , the latter quantifying scalar mixing. The domain is characterized by relatively well-mixed density layers separated by sharp stable interfaces that are correlated with high vertical shear. While static instability is most prevalent within the mixed layers, much of the scalar mixing is localized to the intervening interfaces, a phenomenon not apparent if considering local static instability or $\epsilon$ alone. While the majority of the domain is characterized by the canonical flux coefficient $\varGamma \equiv \chi /\epsilon =0.2$ , often assumed in ocean mixing parametrizations, extreme values of $\chi$ within the statically stable interfaces, associated with elevated $\varGamma$ , strongly skew the bulk statistics. Our findings suggest that current parametrizations of turbulent mixing may be biased by undersampling, such that the most common, but not necessarily the most significant, mixing events are overweighted. Having focused here on a single simulation of stratified turbulence, it is hoped that our results motivate a broader investigation into the role played by stable density interfaces in mixing, across a wider range of parameters and forcing schemes representative of ocean turbulence. [ABSTRACT FROM AUTHOR]

Copyright of Journal of Fluid Mechanics is the property of Cambridge University Press and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

Linear and nonlinear stability analyses are performed to determine critical Rayleigh numbers (${Ra}_{cr}$) for a Rayleigh–Bénard convection configuration with an imposed bottom boundary heat flux that varies harmonically in time with zero mean. The ${Ra}_{cr}$ value depends on the non-dimensional frequency $\omega$ of the boundary heat-flux modulation. Floquet theory is used to find ${Ra}_{cr}$ for linear stability, and the energy method is used to find ${Ra}_{cr}$ for two different types of nonlinear stability: strong and asymptotic. The most unstable linear mode alternates between synchronous and subharmonic frequencies at low $\omega$ , with only the latter at large $\omega$. For a given frequency, the linear stability ${Ra}_{cr}$ is generally higher than the nonlinear stability ${Ra}_{cr}$ , as expected. For large $\omega$ , ${Ra}_{cr} \omega ^{-2}$ approaches an $O(10)$ constant for linear stability but zero for nonlinear stability. Hence the domain for subcritical instability becomes increasingly large with increasing $\omega$. The same conclusion is reached for decreasing Prandtl number. Changing temperature and/or velocity boundary conditions at the modulated or non-modulated plate leads to the same conclusions. These stability results are confirmed by selected direct numerical simulations of the initial value problem. [ABSTRACT FROM AUTHOR]

We study the dynamics of capillary waves at the interface of a two-layer stratified turbulent channel flow. We use a combined pseudo-spectral/phase field method to solve for the turbulent flow in the two liquid layers and to track the dynamics of the liquid–liquid interface. The two liquid layers have same thickness and same density, but different viscosity. We vary the viscosity of the upper layer (two different values) to mimic a stratified oil–water flow. This allows us to study the interplay between inertial, viscous and surface tension forces in the absence of gravity. In the present set-up, waves are naturally forced by turbulence over a broad range of scales, from the larger scales, whose size is of order of the system scale, down to the smaller dissipative scales. After an initial transient, we observe the emergence of a stationary capillary wave regime, which we study by means of temporal and spatial spectra. The computed frequency and wavenumber power spectra of wave elevation are in line with previous experimental findings and can be explained in the frame of the weak wave turbulence theory. Finally, we show that the dispersion relation, which gives the frequency ($\omega$) as a function of the wavenumber ($k$), is in good agreement with the well-established theoretical prediction, $\omega (k) \sim k^{3/2}$. [ABSTRACT FROM AUTHOR]

We conduct three-dimensional direct numerical simulations to investigate the mixing, entrainment and energy budgets of gravity currents emerging from two-layer stratified locks. Depending on the density and layer thickness ratios, we find that either the upper layer or lower layer fluid can propagate faster, and that the density structure of the overall gravity current can range from strongly stratified to near-complete mixing. We furthermore observe that intermediate values of the density ratio can maximise mixing between the gravity current layers. Based on the vorticity budget, we propose a theoretical model for predicting the overall gravity current height, along with the front velocity of the two layers, for situations in which the lower layer moves faster than the upper layer. The model identifies the role of the height and thickness ratios in determining the velocity structure of the current, and it clarifies the dynamics of the ambient counter-current. A detailed analysis of the energy budget quantifies the conversion of potential into kinetic energy as a function of the governing parameters, along with the energy transfer between the different layers of the gravity current and the ambient fluid. Depending on the values of the density and layer thickness ratios, we find that the lower lock layer can gain or lose energy, whereas the upper layer always loses energy. [ABSTRACT FROM AUTHOR]

The formation of a layered structure in the form of vertically separated density steps (staircases) in stably stratified fluids has been reported in many laboratory and oceanic studies as well as in the terrestrial atmospheric boundary layer (ABL) to a lesser extent, with attribution to different dynamical mechanisms. This paper presents observations of layered structures in fog-laden marine ABL, where both fog and density steps appear almost simultaneously following a turbulent mixing event under nocturnal conditions. The observations were made during the C-FOG (2018) field campaign aboard a research vessel using rawinsonde launches, aided by a suite of supporting onboard instruments. This is a case of great practical interest because of the impediment by fog-laden staircases to optical and near-infrared wave propagation in the ABL due to enhanced beam jitter by density steps and beam attenuation by fog. A new mechanism is proposed to explain the genesis of density layering, wherein steps appear when fluid parcels with significant buoyancy differences ( Δ b ) osculate in regions of weak turbulence (local length and velocity scales, L H and u H , respectively) devoid of adequate inertial forces ( ∼ u H 2 / L H ) to cause fluid parcels to stir past each other. This is expressed in terms of a local bulk Richardson number criterion R i = Δ b L H / u H 2 > R i c , where Ri c is a critical value. A simple laboratory experiment with an idealized (three layer) density stratification and a known turbulence source (oscillating grid) was performed to demonstrate the proposed mechanism, and through a combination of measurements and modeling it was found Ri c ≈ 1.5. The proposed criterion was consistent with C-FOG field observations as well as representative previous layering observations in the atmosphere and ocean due to localized turbulence mixing events caused by Kelvin–Helmholtz billowing. [ABSTRACT FROM AUTHOR]

Analysis of satellite altimetry and Argo float data leads Ni et al. (J. Geophys. Res. , 125, 2020, e2020JC016479) to argue that mesoscale dipoles are widespread features of the global ocean having a relatively uniform three-structure that can lead to strong vertical exchanges. Almost all the features of the composite dipole they construct can be derived from a model for multipoles in the surface quasi-geostrophic equations for which we present a straightforward novel solution in terms of an explicit linear algebraic eigenvalue problem, allowing simple evaluation of the higher radial modes that appear to be present in the observations and suggesting that mass conservation may explain the observed frontogenetic velocities. [ABSTRACT FROM AUTHOR]

The present paper investigates theoretically and experimentally the boundary layer generated by a stably stratified fluid flowing horizontally along a surface tilted in the transverse direction and deformed by sinusoidal undulations with crests perpendicular to the flow direction. In the absence of undulations, a weak transverse velocity proportional to the normal velocity is created such that the flow remains purely horizontal. In the presence of undulations of amplitude $h$ , a stronger transverse flow is generated that exhibits a singular behaviour at the critical altitude where the frequency of the perturbation matches the buoyancy frequency of the fluid. This baroclinic critical layer was previously analysed by Passaggia et al. (J. Fluid Mech. , vol. 751, 2014, pp. 663–684) for a boundary layer flow with a small sliding velocity on the surface. Here, the no-slip boundary condition of the experimental flow is applied. For this purpose, we solve the viscous sub-layer to obtain a complete theoretical model for the solution in the critical layer without any adjusting parameter. The theoretical predictions for the transverse velocity are compared with experimental measurements, and a good quantitative agreement is demonstrated. Compared with the sliding case, the no-slip boundary condition on the surface reduces the amplitude of the critical layer solution by a factor $Re^{-1/3}$ , where the Reynolds number Re is defined using the velocity at infinity and the thickness of the boundary layer. As a consequence, the transverse velocity has a maximum in the critical layer of order $h$ , but it still induces a shear rate of order $h\,Re^{1/3}$. [ABSTRACT FROM AUTHOR]

We present an investigation of the resonance conditions governing triad interactions of cylindrical internal waves, i.e. Kelvin modes, described by Bessel functions. Our analytical study, supported by experimental measurements, is performed both in confined and unconfined axisymmetric domains. We are interested in two conceptual questions: can we find resonance conditions for a triad of Kelvin modes? What is the impact of the boundary conditions on such resonances? In both the confined and unconfined cases, we show that sub-harmonics can be spontaneously generated from a primary wave field if they satisfy at least a resonance condition on their frequencies of the form $\omega _0 = \pm \omega _1 \pm \omega _2$. We demonstrate that the resulting triad is also spatially resonant, but that the resonance in the radial direction may not be exact in confined geometries due to the prevalence of boundary conditions – a key difference compared with Cartesian plane waves. [ABSTRACT FROM AUTHOR]

In a three-layer system with equal upper and lower layer thicknesses that are sufficiently thin and with the same density difference across each interface, breathers have been shown to exist using fully nonlinear governing equations. These breathers are well modelled by theoretical solutions of the mKdV equation, provided the interfaces between the layers do not cross a critical depth. The soliton-like characteristics of fully nonlinear breathers, in particular how two breathers interact, have yet to be explored. Using numerical simulations, this study addresses this shortcoming by studying fully nonlinear overtaking collisions of two breathers in a three-layer symmetric stratification. We apply the fully nonlinear and strongly dispersive FDI-3s internal wave equations, based on a variational principle, in a three-layer system. When the amplitude is small, the analytic breathers fit the wave shapes of the overtaking collision breathers. We find that the larger the upper and lower layer thicknesses are, provided they are below the critical thickness, the more the breathers behave like solitons. We show that an overtaking collision of two breathers is close to elastic. [ABSTRACT FROM AUTHOR]