Integral Model for Turbulent Buoyant Jets in Unbounded Stratified Flows. Part I: Single Round Jet

The mechanics of buoyant jet flows issuing with a general three-dimensional geometry into an unbounded ambient environment with uniform density or stable density stratification and under stagnant or steady sheared current conditions is investigated. An integral model is formulated for the conservation of mass, momentum, buoyancy and scalar quantities in the turbulent jet flow. The model employs an entrainment closure approach that distinguishes between the separate contributions of transverse shear (leading to jet, plume, or wake internal flow dynamics) and of azimuthal shear mechanisms (leading to advected momentum puff or thermal flow dynamics), respectively. Furthermore, it contains a quadratic law turbulent drag force mechanism as suggested by a number of recent detailed experimental investigations on the dynamics of transverse jets into crossflow. The model is validated in several stages: First, comparison with basic experimental data for the five asymptotic, self-similar stages of buoyant jet flows, i.e., the pure jet, the pure plume, the pure wake, the advected line puff, and the advected line thermal, support the choice and magnitude of the turbulent closure coefficients contained in the entrainment formulation. Second, comparison with many types of non-equilibrium flows support the proposed transition function within the entrainment relationship, and also the role of the drag force in the jet deflection dynamics. Third, a number of spatial limits of applicability have been proposed beyond which the integral model necessarily becomes invalid due to its parabolic formulation. These conditions, often related to the breakdown of the boundary layer nature of the flow, describe features such as terminal layer formation in stratification, upstream penetration in jets opposing a current, or transition to passive diffusion in a turbulent ambient shear flow. Based on all these comparisons, that include parameters such as trajectories, centerline velocities, concentrations and dilutions, the model appears to provide an accurate and reliable representation of buoyant jet physics under highly general flow conditions.     

Vertical Buoyant Jets in a Linearly Stratified Ambient Cross-Stream

In this study a numerical simulation is performed to investigate the effect of ambient density stratification on the characteristic of a vertical buoyant jet in a stably linearly stratified ambient cross-stream. Based on the ensemble integral method, the theoretical formulation for such a flow field consists of a set of elliptic Reynolds-averaged equations incorporating with the k– transport equations for the turbulence closure. An oscillating motion can be observed in the computed jet trajectory, and the corresponding alternative variation of dominant quantities for the induced momentum and buoyancy of the jet are examined by direct integration on a cross-section along the jet axis. The influences on the jet development both by the ambient cross-stream and the stratification are investigated. The oscillation characteristic shows that a linear relation holds between the wavenumber of jet trajectory, crossflow velocity and the Brunt–Väisälä frequency of ambient stratification. Computational results indicate that the formation of the secondary and a third pairs of vortices, which are not induced in the unstratified environment, causes the jet flow oscillation from its maximum height-of-rise in the flowing direction. The ambient stratification prohibits the growth of the plume radius and reduces the mixing rate as well as the plume rise. The developed flow indicates the transformation of entrainment mechanism in stratified crossflow.     

Integrated modeling of single port brine discharges into unstratified stagnant ambient

An integrated model is presented for the calculation of the characteristics in the intermediate field region of brine discharges from reverse osmosis desalination plants into unstratified stagnant coastal waters. The model consists of the near field model Modified CorJet Model and the far field model, which are interconnected via a coupling algorithm. This algorithm has been developed to simulate the flow and concentration characteristics of negatively buoyant jets (NBJ) after their impingement on the bottom. The coupling method was developed to be active according to literature, however further work and investigation is needed to be applicable for NBJ discharged into other ambient environments and especially in cases where the background values of ambient flow and concentrations affect the NF values and vice versa. The integrated model was validated with data from the literature as well as with data from experiments conducted in this study showing a good agreement. The coupling algorithm was also compared to other coupling techniques used in the literature for NBJ discharges showing better estimations of the experimental data.     

Influence of Downstream Control and Limited Depth on Flow Hydrodynamics of Impinging Buoyant Jets

Much study has been performed on the mixing properties of submerged, turbulent buoyant jets. It is safe to say that the problem of estimating dilution rates in vertical buoyant jets spreading in an `infinitely deep' ambient water has been more than adequately resolved by previous researchers. However, the majority of environmental applications involve discharges into ambient waters of finite depths in which a bounding surface serves to re-direct the impinging buoyant jet horizontally into a radial spreading layer. Previous research indicates that this impinging jet undergoes additional mixing before buoyancy stabilizes vertical mixing and confines the spreading layer to the vicinity of the bounding surface. Unfortunately, the conceptualization and subsequent mathematical modeling of this additional mixing phenomenon is surrounded by considerable amount of disagreement between researchers. The purpose of this study is to provide, by means of velocity and concentration profile measurements, independent experimental evidence for the existence of a critical flow state immediately downstream of the active mixing zone in the horizontally flowing, radial flow that forms after impingement. It is further shown that this critical flow state must be expressed in terms of a composite Froude Number that takes into account the possibility of a non-zero exchange layer flow. Finally, the influence of the presence of a sill-like topographic downstream control on the criticality of the radial flow immediately downstream of the active mixing zone is also investigated.     

Stability and Mixing of a Vertical Axisymmetric Buoyant Jet in Shallow Water

The stability, mixing and effect of downstream control on axisymmetric turbulent buoyant jets discharging vertically into shallow stagnant water is studied using 3D Reynolds-averaged Navier–Stokes equations (RANS) combined with a buoyancy-extended k –ε model. The steady axisymmetric turbulent flow, temperature (or tracer concentration) and turbulence fields are computed using the finite volume method on a high resolution grid. The numerical predictions demonstrate two generic flow patterns for different turbulent heated jet discharges and environmental parameters (i) a stable buoyant discharge with the mixed fluid leaving the vertical jet region in a surface warm water layer; and (ii) an unstable buoyant discharge with flow recirculation and re-entrainment of heated water. A stratified counterflow region always appears in the far-field for both stable and unstable buoyant discharges. Provided that the domain radius L exceeds about 6H, the near field interaction and hence discharge stability is governed chiefly by the jet momentum length scale to depth ratio l_M/H, regardless of downstream control. The near field jet stability criterion is determined to be l_M/H = 3.5. A radial internal hydraulic jump always exists beyond the surface impingement region, with a 3- to 6-fold increase in dilution across the jump compared with vertical buoyant jet mixing. The predicted stability category, velocity and temperature/concentration fields are well-supported by experiments of all previous investigators.     

A simplified computational analysis of turbulent plumes and jets

The application of computational fluid dynamics (CFD), particularly Large Eddy Simulation, for the modelling of buoyant turbulent plumes, has been demonstrated to be very accurate, but computationally expensive. Here a more basic, and therefore more generally practicable, approach is presented. Commercial CFD software is used to model such plumes using Reynolds-Averaged Navier-Stokes (RANS) turbulence models. A careful comparison is made between the numerical predictions and well-established results regarding the bulk properties of plumes. During this process, we are able to observe the well-known approximate Gaussian nature of the plume and achieve quantitative agreement with empirical plume spread coefficients. The use of numerical modelling allows for the investigation of the flow field and turbulence in those regions of the plume of most interest—the plume edge and near source regions. A comprehensive sensitivity study is conducted to identify the limits of applicability of this modelling approach. It is shown that the standard modelling approach of Morton, Taylor and Turner, which introduced the well-known entrainment assumption, pertains in a region well above the source region. At the plume edge, the levels of turbulence are contrasted with the value of the entrainment parameter. Finally, the effects of forcing the plumes with additional momentum at the source are considered, including the case of a pure jet. We show how these forced plumes eventually lose their momentum excess and tend to the behaviour of a pure, buoyant plume.     

Initial mixing of inclined dense jet in perpendicular crossflow

A comprehensive experimental investigation for an inclined ( $60^{\circ }$ to vertical) dense jet in perpendicular crossflow—with a three-dimensional trajectory—is reported. The detailed tracer concentration field in the vertical cross-section of the bent-over jet is measured by the laser-induced fluorescence technique for a wide range of jet densimetric Froude number $Fr$ and ambient to jet velocity ratios $U_r$ . The jet trajectory and dilution determined from a large number of cross-sectional scalar fields are interpreted by the Lagrangian model over the entire range of jet-dominated to crossflow-dominated regimes. The mixing during the ascent phase of the dense jet resembles that of an advected jet or line puff and changes to a negatively buoyant thermal on descent. It is found that the mixing behavior is governed by a crossflow Froude number $\mathbf{F} = U_r Fr$ . For $\mathbf{F} < 0.8$ , the mixing is jet-dominated and governed by shear entrainment; significant detrainment occurs and the maximum height of rise $Z_{max}$ is under-predicted as in the case of a dense jet in stagnant fluid. While the jet trajectory in the horizontal momentum plane is well-predicted, the measurements indicate a greater rise and slower descent. For $\mathbf{F} \ge 0.8$ the dense jet becomes significantly bent-over during its ascent phase; the jet mixing is dominated by vortex entrainment. For $\mathbf{F} \ge 2$ , the detrainment ceases to have any effect on the jet behavior. The jet trajectory in both the horizontal momentum and buoyancy planes are well predicted by the model. Despite the under-prediction of terminal rise, the jet dilution at a large number of cross-sections covering the ascent and descent of the dense jet are well-predicted. Both the terminal rise and the initial dilution for the inclined jet in perpendicular crossflow are smaller than those of a corresponding vertical jet. Both the maximum terminal rise $Z_{max}$ and horizontal lateral penetration $Y_{max}$ follow a $\mathbf{F}^{-1/2}$ dependence in the crossflow-dominated regime. The initial dilution at terminal rise follows a $S \sim \mathbf{F}^{1/3}$ dependence.     

Tracer distributions of line advected thermals

A Light Attenuation system is employed to study the behaviour of strongly advected buoyant jets, focusing in particular on line advected thermals. A double-Gaussian approximation is employed to characterize the cross-sectional form of the tracer distributions. The new data highlights distinct differences in the form of line advected thermals, when compared to advected line momentum puffs. Thermals have a wider peak separation and a relatively large cross-sectional distortion that decays as the flow becomes fully established. It is therefore possible to improve predictions of mean minimum dilution, spread and flow path through separate parameterization of these strongly-advected flows. However, there are notable inconsistencies in recently obtained trajectory data, which require further investigation before such improvements can be made.     

2D brine sewage after impinging on a shallow sea free surface

We consider the problem of the vertically upwards disposal of heavy brine sewage from a two-dimensional diffuser in a lighter, homogeneous, motionless and shallow ambient sea. The rejected high salinity water of seawater desalination plants for urban and agricultural uses is such a case of a two dimensional fountain. The disposal of brine sewage produces a negative buoyant jet due to its initial momentum, which impinges on the free surface, spreads laterally on it and then sinks downwards, because of the negative buoyancy. Laboratory experiments and dimensional considerations are used in this paper in order to investigate the spreading behavior (width) of the vertical fountain which impinges on the free surface of the shallow ambient fluid. The experimental results have been used to derive an equation relating the width at the free surface with the initial parameters of the flow. In addition, the experimentally measured dilution of the heavier brine sewage on the recipient’s surface is compared with the dilution which was calculated by a numerical simulation of a well-known commercial software package, CORJET (a CORMIX sub model).     

Small-scale entrainment in inclined gravity currents

We investigate the effect of buoyancy on the small-scale aspects of turbulent entrainment by performing direct numerical simulation of a gravity current and a wall jet. In both flows, we detect the turbulent/nonturbulent interface separating turbulent from irrotational ambient flow regions using a range of enstrophy iso-levels spanning many orders of magnitude. Conform to expectation, the relative enstrophy isosurface velocity \(v_n\) in the viscous superlayer scales with the Kolmogorov velocity for both flow cases. We connect the integral entrainment coefficient E to the small-scale entrainment and observe excellent agreement between the two estimates throughout the viscous superlayer. The contribution of baroclinic torque to \(v_n\) is negligible, and we show that the primary reason for reduced entrainment in the gravity current as compared to the wall-jet are 1) the reduction of \(v_n\) relative to the integral velocity scale \(u_T\); and 2) the reduction in the surface area of the isosurfaces.     

Surface signatures of submerged heated jet

The flow induced at the surface of a water body by a submerged heated horizontal turbulent jet was investigated experimentally with the aim of developing parameterizations for surface mean temperature/velocity fields. The jet nozzle diameter was fixed, the depth of the jet beneath the free surface was varied, and two jet Reynolds numbers (5020, 11300) were considered. The surface temperature was measured using a highly sensitive infrared camera, and the near-surface horizontal velocity field was measured using particle image velocimetry. The experimental results were explained using a model based on similarity solutions with variable turbulent viscosity. While classical Schlichting’s solution with constant turbulent viscosity predicts complete similarity for transverse velocity/temperature distributions only in a plane that coincides with the flow axis, the present solution predicts similarity in an arbitrary plane parallel to the flow axis, which was confirmed using data collected at the surface. Comparisons of present data with available previous results also showed general agreement.     

Modeling the buoyant flow of heated water discharged from surface and submerged side outfalls in shallow and deep water with a cross flow

In this study, a three-dimensional model was used to numerically study the buoyant flow, along with its mixing characteristics, of heated water discharged from the surface and submerged side outfalls in shallow and deep water with a cross flow. Hydraulic experimental data were used to evaluate the applicability of the model. The simulation results agree well with the experimental results, particularly, the jet trajectories, the dimensions of the recirculating zone, and the distribution of the dimensionless excess temperature. The level of accuracy of the simulation results of the present study is nearly identical to that of the results conducted by McGuirk and Rodi (1978). If the heated water is discharged into shallow water where the momentum flux ratio and the discharge densimetric Froude number are high, the submerged discharge method is better than the surface discharge method in terms of the scale of the recirculating zone and the minimum dilution. In deep water, where the momentum flux ratio and discharge densimetric Froude number are low, however, the submerged discharge method had few advantages. In shallow water, the discharge jet is deflected by the ambient cross flow, while forcing the ambient flow to bend towards the far bank for the full depth. For a submerged discharge in shallow water, the recirculating zone is the largest in the lowest layer but becomes smaller in the upper layer. As the water depth increases, the ambient flow goes over the jet and diminishes the blocking effect, thereby decreasing the bending of the jet.     

Advanced integral model for groups of interacting round turbulent buoyant jets

An integral model that combines all advantages of Superposition Method (SM), Entrainment Restriction Approach (ERA) and Second Order Approach (SOA) is proposed to predict the mean axial velocity and concentration fields of a group of N interacting vertical round turbulent buoyant jets. SM is successful in predicting the fields of mean axial velocity and mean concentration for a group of N interacting jets or plumes and ERA is advantageous in predicting the above fields for either two or large number (N → ∞) of interacting buoyant jets in the whole range of buoyancy. SOA takes into consideration in a dynamic way the turbulent contribution to the momentum and buoyancy fluxes and provides better accuracy than the usual procedures. A novelty of the proposed model is the production and utilisation of advanced profile distributions, convenient for the mean axial velocities and concentrations in a cross-section of the entire group of buoyant jets. These profiles are developed on the basis of flux conservation of momentum, buoyancy and kinetic energy for the mean motion. They enhance dynamic adaptation of the individual buoyant jet axes to the group centreline. Due to these profile distributions, the present model owns generality of application and better accuracy of predictions compared to usual integral models using simple Gaussian or top-hat profiles; thus it conferred the name Advanced Integral Model (AIM). AIM is herein applied to predict the mean flow properties of two different arrangement types of any number of buoyant jets: (a) linear diffusers and (b) rosette-type risers. Present results are compared to available experimental data and traditional solutions based on Gaussian profiles. Findings may be useful for design purposes and environmental impact assessment.     

Laboratory studies defining flow regimes for negatively buoyant surface discharges into crossflow

Surface discharges of negatively buoyant jets into moving ambient water create a range of complex flow patterns. These complexities arise through the interplay between the discharge’s initial fluxes and the motion of the ambient current. In this study a series of laboratory experiments were conducted for negatively buoyant surface discharges into crossflow to investigate flow patterns under different discharge and ambient conditions. The results compared with simulations of the CORMIX model, an expert system for ocean outfall design. In CORMIX, the simulation module DHYDRO for dense discharges has been used. Finally the flow different patterns were arranged in a dimensionless diagram to propose a modified flow classification system with new criteria.     

Characterising strongly advected discharges in the initial dilution zone

Mean concentration fields of strongly advected non-buoyant discharges are characterised with a double-Gaussian assumption. Comparisons with experimental data show that the approximation provides a reasonable representation of the cross-sectional profiles. The self-similarity of these profiles enables their form to be represented by two additional parameters, one describing the relative separation of the peaks and the other the ratio of the cross-sectional spreads. Values for these additional parameters are determined from experimental data. This systematic approach to characterising the strongly advected flows provides a consistent framework for determining spreading rates and concentration ratios, such as the peak to centreline maximum and the peak to top hat. The double-Gaussian framework also provides a basis for comparisons with the CorJet and VisJet numerical models. In addition the double-Gaussian assumption is employed to interpret data obtained using the Light Attenuation technique. This is a relatively simple measuring system, which provides depth integrated concentration information. The data obtained using this technique is shown to be generally consistent with that from previous studies.     

Evaluating the dimensionality and significance of “active periods” in turbulent environmental flows defined using Lipshitz/Hölder regularity

Active periods within perturbed boundary-layer flows are considered in terms of the local roughness of measured velocity time series and defined in terms of Hölder/Lipshitz exponents. Such events are associated with the passage of energetic, coherent flow structures and are responsible for exerting high turbulent stresses because of the rapid changes in velocity that occur at such times. A method is proposed for assessing the effective dimensionality of such active periods, as well as their significance to the flow field, for a particular choice of flow metric. The method is applied to the turbulent flow through a confluence flow geometry, with velocity samples acquired close to the bed of the channel in a zone of complex mixing. The dimensionality of the active periods is consistent with the observed patterns of sediment entrainment from the bed, with the significance of the active periods decaying away from the erosional zone.     

PIV experimental study on the flow characteristics upstream of a floating intake in nonlinear stratified ambient conditions

Selective withdrawal is commonly implemented in nonlinearly stratified ambient, which typically has stratified ambient conditions, for purposes of controlling quality. A floating intake is applied as an effective facility of selective withdrawal. However, the outflow dynamics of a floating intake in a nonlinearly stratified ambient have been disregarded, which has a significant effect on the outflow water quality of a reservoir. Experiments were conducted to investigate the effect of thermal stratification on the flow characteristics using particle image velocimetry at three temperature distributions (no stratification, weak stratification and strong stratification). The flow fields upstream of the floating intake showed that the withdrawal layer was formed inhibited by the thermal stratification. And strong stratification produced the thinner withdrawal layer thickness, leading to a larger nonuniform coefficient of the velocity profile. To quantitatively describe the velocity profiles, formulas of dimensionless velocity profiles were proposed. The flow developments were analysed, and the virtual control points located 0.56d above the floating intake (where d is the straight pipe diameter of the floating intake) were obtained. The positions of virtual control points mainly depended on the withdrawal discharge. The decay rate of the velocity along the horizontal line passing through the virtual control point was inversely proportional to the stratification intensity.

     

Modelling hydraulic jump using the bubbly two-phase flow method

Hydraulic jumps have complex flow structures, characterised by strong turbulence and large air contents. It is difficult to numerically predict the flows. It is necessary to bolster the existing computer models to emphasise the gas phase in hydraulic jumps, and avoid the pitfall of treating the phenomenon as a single-phase water flow. This paper aims to improve predictions of hydraulic jumps as bubbly two-phase flow. We allow for airflow above the free surface and air mass entrained across it. We use the Reynolds-averaged Navier–Stokes equations to describe fluid motion, the volume of fluid method to track the interface, and the k–ε model for turbulence closure. A shear layer is shown to form between the bottom jet flow and the upper recirculation flow. The key to success in predicting the jet flow lies in formulating appropriate bottom boundary conditions. The majority of entrained air bubbles are advected downstream through the shear layer. Predictions of the recirculation region’s length and air volume fraction within the layer are validated by available measurements. The predictions show a linear growth of the shear layer. There is strong turbulence at the impingement, and the bulk of the turbulence kinetic energy is advected to the recirculation region via the shear layer. The predicted bottom-shear-stress distribution, with a peak value upstream of the toe of the jump and a decaying trend downstream, is realistic. This paper reveals a significant transient bottom shear stress associated with temporal fluctuations of mainly flow velocity in the jump. The prediction method discussed is useful for modelling hydraulic jumps and advancing the understanding of the complex flow phenomenon.     

Three-Dimensional Modelling of Concentration Fluctuations in Complicated Geometry

The strong fluctuating component in the measured concentration time series of a dispersing gaseous pollutant in the atmospheric boundary layer, and the hazard level associated to short-term concentration levels, demonstrate the necessity of calculating the magnitude of turbulent fluctuations of concentration using computational simulation models. Moreover the computation of concentration fluctuations in cases of dispersion in realistic situations, such as built-up areas or street canyons, is of special practical interest for hazard assessment purposes. In this paper, the formulation and evaluation of a model for concentration fluctuations, based on a transport equation, are presented. The model is applicable in cases of complex geometry. It is included in the framework of a computational code, developed for simulating the dispersion of buoyant pollutants over complex geometries. The experimental data used for the model evaluation concerned the dispersion of a passive gas in a street canyon between 4 identical rectangular buildings performed in a wind tunnel. The experimental concentration fluctuations data have been derived from measured high frequency concentrations. The concentration fluctuations model is evaluated by comparing the model's predictions with the observations in the form of scatter plots, quantile-quantile plots, contour plots and statistical indices as the fractional bias, the geometrical mean variance and the factor-of-two percentage. From the above comparisons it is concluded that the overall model performance in the present complex geometry case is satisfactory. The discrepancies between model predictions and observations are attributed to inaccuracies in prescribing the actual wind tunnel boundary conditions to the computational code.     

Combined integral and particle model for describing the dispersion,dilution, terminal layer formation and influence area from a point source discharge into a water body

Contamination of coastal water is a persistent threat to ecosystems around the world. In this study, a novel model for describing the dispersion, dilution, terminal layer formation and influence area from a point source discharge into a water body is presented and compared with field measured data. The model is a Combined Integral and Particle model (CIPMO). In the initial stage, the motion, dispersion and dilution of a buoyant jet are calculated. The output from the buoyant jet model is then coupled with a Lagrangian Advection and Diffusion model describing the far-field. CIPMO ensures that both the near- and far-field processes are adequately resolved. The model either uses empirical data or collects environmental forcing data from open source hydrodynamic models with high spatial and temporal resolution. The method for coupling the near-field buoyant jet and the particle tracking model is described and the output is discussed. The model shows good results when compared with measurements from a field study.