Integral models for bubble,droplet, and multiphase plume dynamics in stratification and crossflow

We present the development and validation of a numerical modeling suite for bubble and droplet dynamics of multiphase plumes in the environment. This modeling suite includes real-fluid equations of state, Lagrangian particle tracking, and two different integral plume models: an Eulerian model for a double-plume integral model in quiescent stratification and a Lagrangian integral model for multiphase plumes in stratified crossflows. Here, we report a particle tracking algorithm for dispersed-phase particles within the Lagrangian integral plume model and a comprehensive validation of the Lagrangian plume model for single- and multiphase buoyant jets. The model utilizes literature values for all entrainment and spreading coefficients and has one remaining calibration parameter \( \kappa \), which reduces the buoyant force of dispersed phase particles as they approach the edge of a Lagrangian plume element, eventually separating from the plume as it bends over in a crossflow. We report the calibrated form \( \kappa = [(b - r) / b]^4 \), where b is the plume half-width, and r is the distance of a particle from the plume centerline. We apply the validated modeling suite to simulate two test cases of a subsea oil well blowout in a stratification-dominated crossflow. These tests confirm that errors from overlapping plume elements in the Lagrangian integral model during intrusion formation for a weak crossflow are negligible for predicting intrusion depth and the fate of oil droplets in the plume. The Lagrangian integral model has the added advantages of being able to account for entrainment from an arbitrary crossflow, predict the intrusion of small gas bubbles and oil droplets when appropriate, and track the pathways of individual bubbles and droplets after they separate from the main plume or intrusion layer.     

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.     

Verification of the plume rise/dispersion model USPR: plume rise for single stack emissions

Despite proliferation of the use of air pollution models for regulatory application, major discrepancies still occur between models and also between models and observations, especially when oversimplistic models are used. The problem of predicting plume rise (and subsequently ground level concentrations) from a single source is evaluated here in terms of an integral plume rise and dispersion model (USPR) which encompasses both bouyant rise and turbulent spreading; thus avoiding the problems of the concatenation of separate plume rise and dispersion models. The wide range of validity of the USPR model is demonstrated is terms of plume rise by comparison with the highly buoyant GCOS and Kincaid plumes as well as with dense effluents. It is also shown to be in agreement with Briggs' two-thirds law when the restrictions applicable to the latter model are imposed.     

Cooling tower plume abatement and plume modeling: a review

Visible plumes above wet cooling towers are of great concern due to the associated aesthetic and environmental impacts. The parallel path wet/dry cooling tower is one of the most commonly used approaches for plume abatement, however, the associated capital cost is usually high due to the addition of the dry coils. Recently, passive technologies, which make use of free solar energy or the latent heat of the hot, moist air rising through the cooling tower fill, have been proposed to minimize or abate the visible plume and/or conserve water. In this review, we contrast established versus novel technologies and give a perspective on the relative merits and demerits of each. Of course, no assessment of the severity of a visible plume can be made without first understanding its atmospheric trajectory. To this end, numerous attempts, being either theoretical or numerical or experimental, have been proposed to predict plume behavior in atmospheres that are either uniform versus density-stratified or still versus windy (whether highly-turbulent or not). Problems of particular interests are plume rise/deflection, condensation and drift deposition, the latter consideration being a concern of public health due to the possible transport and spread of Legionella bacteria.

     

Lagrangian models of dispersion in marine environment

Turbulent dispersion can be studied successfully by using Lagrangian particle models. In general, the prediction of correct concentration fields is a complex issue when the turbulent field is inhomogeneous and non-stationary. Two classes of Lagrangian dispersion models have been considered in this work, which are based on the Wiener process and the so called “well-mixed” criterion. In order to test the performances of these models and shed light on the underlying physical processes and modeling assumptions, four different numerical models have been compared and tested by means of their long time behavior by considering several study cases concerning idealized marine environment. Furthermore, the coupling of the community model Princeton Ocean Model (POM) with the Lagrangian model LASEMOD (LAgrangian SEa MODel) is used to investigate the temporal and spatial evolution of a passive pollutant released in the vicinity of the coast in the Tyrrhenian Sea basin. The simulation shows with reasonable accuracy the time evolution of both the hydrodynamic and the concentration fields and provides a useful insight into the evaluation of the environmental impact of pollutant releases along the coast.     

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.     

A Numerical Study of the Dynamics of the Riverine Thermal Bar in a Deep Lake

A numerical model based on a Finite Volume formulation of the Navier–Stokes equations is used to simulate a range of scenarios leading to a thermal bar formed by a river inflow to an idealised deep lake. The results presented here show that small riverine salinity increases have a profound effect on the dynamics of the thermal bar, suppressing horizontal propagation of the plume and raising the possibility of a thermal bar which is capable of sinking to great depths. This finding is particularly relevant to Lake Baikal in Siberia, where the vigorous deep-water renewal is still not fully understood. An analysis of the buoyancy forces governing the depth of penetration of the thermal bar plume shows that realistic salinity gradients are an important factor in determining the circulation of Baikal waters. Observations of the saline curtailment of the thermal bar's horizontal propagation also reveal a potential for reduced productivity in the ecosystem of any temperate river delta during the Spring renewal period.     

Wave mixing rise inferred from Lyapunov exponents

We study the horizontal surface mixing and the transport induced by waves in a coastal environment. A comparative study is addressed by computing the Lagrangian coherent structures, via Finite Size Lyapunov Exponents, that arise in two different numerical settings: with and without wave coupled to currents. In general, we observe that mixing is increased in the area due to waves. Besides, the methodology presented here is tested by deploying a set of eight Lagrangian drifters at different locations. This dynamical approach is shown as a valuable tool to extract information about transport, mixing and residence embedded in the Eulerian time dependent velocity fields obtained from numerical models.     

Joint use of data and modeling in coastal wave transformation

In the framework of a research project entitled ??BRISA??BReaking waves and Induced SAnd transport??, a methodology was devised to characterize the waves joining together in-situ measurements and numerical wave propagation models. With this goal in mind, a number of in-situ measurements were made, for selected positions in front of Praia de Faro (South Portugal), during four days (25th to 28th March, 2009) by using different types of equipments (e.g., resistive wave gauges, pressure sensors, currentmeters and a new prototype pore pressure sensor using optical fibre). Wave records were obtained simultaneously offshore (at a water depth of 11.7?m below mean sea level, MSL) and at the surf and swash zones. The data processing and analysis were made by applying classical time domain techniques. Numerical simulations of the wave propagation between offshore and inshore for the measurement period were performed with two numerical models, a 1D model based on linear theory and a nonlinear Boussinesq-type model, COULWAVE, both forced by the measured offshore wave conditions of 27th March 2009. Comparisons between numerical results and field data for the pressure sensors placed in the surf and swash zones were made and discussed. This approach enables to evaluate the performance of those models to simulate those specific conditions, but also to validate the models by gaining confidence on their use in other conditions.     

The influence of hydrodynamic processes on zooplankton transport and distributions in the North Western Mediterranean: Estimates from a Lagrangian model

A Lagrangian module has been developed and coupled with the 3D circulation model Symphonie to study the influence of hydrodynamic processes on zooplankton transport and distributions in the North Western Mediterranean (NWM). Individuals are released every 3 days from March to August 2001 in two initial areas: around the DYFAMED sampling station in the central Ligurian Sea and in the Rhône river plume. Then the individuals are tracked for 40 days either as passive particles or with a simple diel vertical migration (DVM) pattern. The simulations suggest strong seasonal patterns in the distributions of the individuals released around the DYFAMED sampling station. Individuals spread all over the NWM basin after 40 days but different patterns occur depending on the season, the initial depths of release and the capacity of DVM. Offshore-shelf transport only occurs in April and May with particles ending up in the Gulf of Lions (GoL) in low concentrations. In other months, the Northern Current (NC) can be considered as a barrier for particles entering the GoL from the offshore sea. A quarter to a half of passive individuals released in the Rhône river plume remain in the GoL. The rest is transported by the NC towards the Catalan Sea. Applying a simple DVM scheme does not increase the retention of particles on the shelf.     

Detrainment of plumes from vertically distributed sources

We present experimental results demonstrating that, for the turbulent plume from a buoyancy source that is vertically distributed over the full area of a wall, detrainment qualitatively changes the shape of the ambient buoyancy profile that develops in a sealed space. Theoretical models with one-way-entrainment predict stratifications that are qualitatively different from the stratifications measured in experiments. A peeling plume model, where density and vertical velocity vary linearly across the width of the plume, so that plume fluid “peels” off into the ambient at intermediate heights, more accurately captures the shape of the ambient buoyancy profiles measured in experiments than a conventional one-way-entrainment model does.     

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.     

Dynamics of the buoyant plume off the Pearl River Estuary in summer

Field measurements of salinity, wind and river discharge and numerical simulations of hydrodynamics from 1978 to 1984 are used to investigate the dynamics of the buoyant plume off the Pearl River Estuary (PRE), China during summer. The studies have shown that there are four major horizontal buoyant plume types in summer: Offshore Bulge Spreading (Type I), West Alongshore Spreading (Type II), East Offshore Spreading (Type III), and Symmetrical Alongshore Spreading (Type IV). River mouth conditions, winds and ambient coastal currents have inter-influences to the transport processes of the buoyant plume. It is found that all of the four types are surface-advected plumes by analysing the vertical characteristic of the plumes, and the monthly variations of the river discharge affect the plume size dominantly. The correlation coefficient between the PRE plume size and the river discharge reaches 0.85 during the high river discharge season. A wind strength index has been introduced to examine the wind effect. It is confirmed that winds play a significant role in forming the plume morphology. The alongshore wind stress and the coastal currents determine the alongshore plume spreading. The impact of the ambient currents such as Dongsha Current and South China Sea (SCS) Warm Current on the plume off the shelf has also assessed. The present study has demonstrated that both the river discharge and wind conditions affect the plume evolution.     

Numerical analysis on the Brisbane River plume in Moreton Bay due to Queensland floods 2010–2011

During the Queensland floods in the summer of 2010–2011, a flood-driven Brisbane River plume extended into Moreton Bay, Queensland, Australia, and then seaward, travelling in a northward direction. It covered approximately 500 km $^{2}$ . This paper presents a three- dimensional hydrodynamic numerical model investigation into the behaviour of the Brisbane River plume. The model was verified by using satellite observations and field measurement data. The present study concludes that the high river discharge was the primary factor determining the plume size and its seaward extensions. A notable finding was that the plume was a bottom-trapped type rather than a buoyant type. Further, the southerly winds were found to have moderately confined the alongshore extension of the plume, and had caused the plume to mix thoroughly with the ocean water.     

Filament-Based Atmospheric Dispersion Model to Achieve Short Time-Scale Structure of Odor Plumes

This article presents the theoretical motivation, implementation approach, and example validation results for a computationally efficient plume simulation model, designed to replicate both the short-term time signature and long-term exposure statistics of a chemical plume evolving in a turbulent flow. Within the resulting plume, the odor concentration is intermittent with rapidly changing spatial gradient. The model includes a wind field defined over the region of interest that is continuous, but which varies with location and time in both magnitude and direction. The plume shape takes a time varying sinuous form that is determined by the integrated effect of the wind field. Simulated and field data are compared. The motivation for the development of such a simulation model was the desire to evaluate various strategies for tracing odor plumes to their source, under identical conditions. The performance of such strategies depends in part on the instantaneous response of target receptors; therefore, the sequence of events is of considerable consequence and individual exemplar plume realizations are required. Due to the high number of required simulations, computational efficiency was critically important.     

Laboratory measurements and multi-block numerical simulations of the mean flow and turbulence in the non-aerated skimming flow region of steep stepped spillways

We present and discuss the results of a comprehensive study addressing the non-aerated region of the skimming flow in steep stepped spillways. Although flows in stepped spillways are usually characterized by high air concentrations concomitant with high rates of energy dissipation, the non-aerated region becomes important in small dams and/or spillways with high specific discharges. A relatively large physical model of such spillway was used to acquire data on flow velocities and water levels and, then, well-resolved numerical simulations were performed with a commercial code to reproduce those experimental conditions. The numerical runs benefited from the ability of using multi-block grids in a Cartesian coordinate system, from capturing the free surface with the TruVOF method embedded in the code, and from the use of two turbulence models: the k-e{k{-}\varepsilon} and the RNGk-e{k{-}\varepsilon} models. Numerical results are in good agreement with the experimental data corresponding to three volumetric flow rates in terms of the time-averaged velocities measured at diverse steps in the spillway, and they are in very satisfactory agreement for water levels along the spillway. In addition, the numerical results provide information on the turbulence statistics of the flow. This work also discusses important aspects of the flow, such as the values of the exponents of the power-law velocity profiles, and the characteristics of the development of the boundary layer in the spillway.     

Warm discharges in cold fresh water: 2. Numerical simulation of laminar line plumes

The behaviour of a discharge of warm water upwards into a homogeneous body of cold fresh water was investigated by means of a numerical model. The discharge has a parabolic velocity profile, with Reynolds number \(Re=50\), Prandtl number \(Pr=7\) and Froude number varied over the range \(0.2 \le {\rm Fr} \le 2.5\). Water density is taken to be a quadratic function of temperature, so that an initially positively buoyant discharge will experience buoyancy reversal as it mixes with an ambient below the temperature of maximum density. The resulting plume has some similarities to a fountain resulting from injection of negatively buoyant fluid upward into a less dense ambient. The plume is initially symmetric, but then its head detaches as it approaches its maximum height. The detached head is denser than the fluid in the plume below it, and the interaction between the sinking head and the rising plume causes a sideways deflection; as this cycle is repeated, the plume displays side-to-side flapping motion and vertical bobbing. As Froude number is increased (i.e. buoyancy reduced) the growth of the plume becomes slower, but the plume eventually reaches a greater height. We obtain empirical power-law scalings for maximum height and time taken to reach that height as functions of Froude number; these scalings are simlar to those for fountains with a linear dependence of density on temperature in the very weak regime.     

A stochastic Lagrangian micromixing model for the dispersion of reactive scalars in turbulent flows: role of concentration fluctuations and improvements to the conserved scalar theory under non-homogeneous conditions

Several reaction schemes, based on the conserved scalar theory, are implemented within a stochastic Lagrangian micromixing model to simulate the dispersion of reactive scalars in turbulent flows. In particular, the formulation of the reaction-dominated limit (RDL) reaction scheme is here extended to improve the model performance under non-homogeneous conditions (NHRDL scheme). The validation of the stochastic model is obtained by comparison with the available measurements of reactive pollutant concentrations in a grid-generated turbulent flow. This test case describes the dispersion of two atmospheric reactant species (NO and O₃) and their reaction product (NO₂) in an unbounded turbulent flow. Model inter-comparisons are also assessed, by considering the results of state-of-the-art models for pollutant dispersion. The present validation shows that RDL reaction scheme provides a systematic overestimation (relative error of ca. 85% around the centreline) in computing the local reactant consumption/production rate, whereas the NHRDL scheme drastically reduces this gap (relative error lower than 5% around the centreline). In terms of NO₂ production (or reactant consumption), neglecting concentration fluctuations determines overestimations of the product mean of around 100% and a NO₂ local production of one order of magnitude higher than the reference simulation. In terms of standard deviations, the concentration fluctuations of both the passive and reactive scalars are generally of the same order of magnitude or up to 1 or 2 orders of magnitudes higher than the corresponding ensemble mean values, except for the background reactant close to the plume edges. The study highlights the importance of modelling pollutant reactions depending on the instantaneous instead of the mean concentrations of the reactants, thus quantifying the role of the turbulent fluctuations of concentration, in terms of scalar statistics (mean, standard deviation, intensity of fluctuations, skewness and kurtosis of concentration, segregation coefficient, simulated reaction rate). This stochastic particle method represents an efficient numerical technique to solve the convection–diffusion equation for reactive scalars and involves several application fields: micro-scale air quality (urban and street-canyon scales), accidental releases, impact of odours, water quality and fluid flow industrial processes (e.g. combustion).     

Influence of planform geometry and momentum ratio on thermal mixing at a stream confluence with a concordant bed

The effects of planform geometry and momentum flux ratio on thermal mixing at a stream confluence with concordant bed morphology are investigated based on numerical simulations that can capture the dynamics of large-scale turbulence. In two simulations, the bathymetry and asymmetrical planform geometry are obtained from field experiments and the momentum flux ratio is set at values of one and four. These two conditions provide the basis for studying differences in thermal mixing processes at this confluence when the wake mode and the Kelvin–Helmholtz mode dominate the development of coherent structures within the mixing interface (MI). The effects of channel curvature and angle between the two incoming streams on thermal mixing processes are investigated based on simulations conducted with modified planform geometries. Two additional simulations are conducted for the case where the upstream channels are parallel but not aligned with the downstream channel and for the zero-curvature case where the upstream channels are parallel and aligned with the downstream channel. The simulations highlight the influence of large-scale coherent structures within the MI and of streamwise-oriented vortical (SOV) cells on thermal mixing processes within the confluence hydrodynamics zone. Simulation results demonstrate the critical role played by the SOV cells in promoting large-scale thermal mixing for cases when such cells form in the immediate vicinity of the MI and in modifying the shape of the thermal MI within cross sections of the downstream channel—predictions consistent with empirical measurements of thermal mixing at the confluence. The set of numerical simulations reveal that the degree of thermal mixing occurring within the confluence hydrodynamic zone varies dramatically with planform geometry and incoming flow conditions. In some cases thermal mixing at the downstream end of the confluence hydrodynamic zone is limited to the MI and its immediate vicinity, whereas in others substantial thermal mixing has occurred over most of the cross-sectional area of the flow. Overall, the simulations highlight the flow conditions and the controls of these conditions that influence mixing within the immediate vicinity of a confluence.