Large-eddy simulations of turbidity plumes in crossflow
Introduction
Our specific field of interest is the highly-concentrated sediment-water plume released through the overflow shaft of a trailing suction hopper dredging vessel, e.g. Smith and Friedrichs [1]. The latter type of dredger is widely applied worldwide and allows dredging while sailing. By means of a suction pipe lowered to the sea- or riverbed, a sediment-water mixture is pumped into the hopper. Whereas the coarser sediment particles settle in the hopper, the finest sediment fractions flow overboard with the excess water through a vertical overflow dropshaft, the outlet of which is flush with the keel of the vessel’s hull. As a consequence, a negatively buoyant plume of water and fine sediment particles is released vertically below the vessel. Because of the sailing speed of the vessel and/or the ambient currents, the generated plumes are subjected to a crossflow. To enable proper assessment of the environmental impact of the plumes it is important to be able to predict its fate and dispersion in the water column [2]. The subsequent interaction of the plume with the bed, is outside the scope of this paper.
The objective of this paper is to develop a numerical model capable of accurately predicting the mean trajectory and the turbulent dispersion of a negatively buoyant sediment plume in a crossflow. Both the flow field of a plume in crossflow as well as the flow field of a plume in crossflow influenced by a backward-facing step (BFS) are analysed thoroughly. To this end, a Large Eddy Simulation (LES) approach is selected. The model will be set up and validated based upon general plume data reported in literature, as well as scaled-down laboratory experiments with sediment plumes by the authors (Decrop et al.) [3]. In the latter experiments, the hull of the dredging vessel was highly schematised as a quasi prismatic box with a given draft. The numerical model will be validated both without and with the schematised hull. It was intended to set up a model with acceptable process times to enable future upscaling of the model to prototype scale, including realistic dredging vessel geometry and propeller action.
A buoyant jet or forced plume is a plume in which both an initial momentum flux and a buoyancy flux govern the flow up to a certain distance from the release point [4], [5]. Pure plumes and pure jets, to the contrary, are only forced by a mass density difference and an initial momentum flux, respectively.
For round buoyant jets with top-hat velocity profile, initial volume, momentum and buoyancy fluxes can be written, respectively, as: where is the exit pipe diameter, is the (uniform) exit velocity and is the reduced gravity of the mixture being discharged in an ambient fluid.
In the case presented in this paper, the buoyancy force is generated by the presence of particles in a mixture. Assuming the mass density of the fluid phase in the mixture is equal to the mass density of the ambient fluid, can be written as where is the initial particle volume concentration (), is the mass density of the particles and is the mass density of the fluid phase in the mixture and of the receiving ambient fluid (both tap water in the experiments carried out by the authors).
The quantities , and are considered the primary variables in the study of turbulent jets, plumes and forced plumes [5]. In the remainder of this article, the forced plumes studied will be referred to as ‘plumes’.
Plumes released from a circular pipe or orifice exhibit self-similarity in the region past the so-called zone of flow establishment, typically at about seven pipe diameters [4], [5]. The profiles of flow velocity and tracer concentration collapse when normalised by the appropriate parameters and they follow a Gaussian distribution. This allows researchers to compare experimental and numerical results of plumes of different , and . In this work, the self-similarity property will be used to validate the LES model by simulating a vertical plume and comparing the results with measurements by Decrop et al. [3], which are shown to be consistent with single-phase plume experiments by Papanicolaou and List [6], Shabbir et al. [7] and Dai et al. [8].
The flow in plumes released vertically in a crossflow is characterised by large-scale turbulent structures. These are related to the blending of initial vertical momentum with crossflow horizontal momentum (counter-rotating vortex pair, CRVP) as well as Kelvin–Helmholtz instabilities resulting in leading-edge vortices and associated convection cells ejecting in the direction of the buoyancy force vector (e.g. Tian and Roberts [9], Cambonie et al. [10] and Diez et al. [11]). Additionally, a pattern of wake vortices similar to the vortices in the wake of a cylinder has been observed by Fric and Roshko [12] and Muldoon and Acharya [13] in the wake of a jet in cross flow. These vortices are, however, clearly distinct from a von Karman vortex street since they originate from the crossflow boundary layer. The two-liquid case of these flows has been studied intensively, whereas turbulent plumes in which negative buoyancy is created by fine particles are not very well studied. Much of the large scale flow structure is expected to be similar to the intermiscible liquids case, however, turbulent diffusion of mass is expected to behave differently in dispersed two-phase flows. Different authors have shown that the turbulent Schmidt number, has a different average value for particulate suspensions compared to the standard value of 0.7 for intermiscible liquids diffusion [14], [15], [16], [17]. In principle, a two-phase LES model using the mixture model for the dispersed phase of the spreading of a sediment laden jet should resolve enough of the turbulent motions to provide a spreading rate associated with the diffusivity in a sediment plume. The sub-grid scale turbulent Schmidt number is derived from test-filtered turbulent fluxes [18], as described in Section 2.1.2.
Dimensional analysis leads to two main dimensionless numbers characterising plumes released in a cross-flow, the densimetric Froude number and the velocity ratio .
Trajectories followed by plumes in crossflow are determined by several possible flow regimes: jet regime, bent jet regime, plume regime and the bent plume regime. In some literature, e.g. Wright [19], these are referred to as momentum-dominated near field (MDNF), momentum-dominated far field (MDFF), buoyancy-dominated near field (BDNF) and buoyancy-dominated far field (BDFF) respectively, e.g. Wright [19].
The plumes studied here range within (, ) limits occurring in sediment plumes released from dredging vessels. Even though the initial relative density difference is usually in the order of 1–10%, the buoyancy is relatively weak compared to the crossflow in these cases, with the velocity ratio usually in the range . Therefore, the momentum length scale (Eq. (6)) is larger than the buoyancy length scale (Eq. (7)) in most cases, leading to a plume trajectory sequence MDNF–MDFF–BDFF. However, in strong cross-flow cases both and are around or less than unity, due to which the plume transforms very rapidly to the BDFF regime. The plume trajectory in the two main regimes occurring in the studied flows (MDFF, BDFF) scale with and , respectively.
As mentioned above, many types of turbulent structures occur in plumes in crossflow. Some of the large-scale structures have a steady nature and can also be solved with a Reynolds-Averaged Navier–Stokes (RANS) model. Therefore, one could wonder why a more time-consuming LES model is used in that case. Many of the larger turbulent structures are only resolved in a time-domain model, for example buoyancy bursts referred to as cauliflowers and tracer entrained in a wake vortices (observed by Smith and Mungal [20] in a jet). Moreover, it is known that the turbulence in jets and plumes in cross-flow is anisotropic [21]. This makes the use of RANS model with an isotropic turbulence model less suited. The influence of the geometry of the surrounding walls on these structures can therefore ideally be studied using time-domain solutions. Also, future study of pulsed plumes using the presented model would require time-domain resolution in any case [22], [23]. Frequency analysis of the turbulent fluctuations of velocity and tracer concentration is possible using LES results. Experiments in the past show both and exponential spectral energy cascade laws [8].
LES simulations of jets and plumes in a still environment have been reported repeatedly in the past [24], [25], [26]. LES simulations of a dispersed two-phase plume in still environment were described by Dimitrova et al. [27]. A limited number of studies, though, reports LES simulations of jets or single phase plumes in a crossflow, e.g. by Yuan et al. [28], Recker et al. [29] and Coussement et al. [22]. Muppidi and Mahesh [30], [31] performed Direct Numerical Simulations (DNS) simulations of transverse jets and passive scalar transport therein. To the best of our knowledge, LES of two-phase, small Stokes number particle plumes in crossflow did not receive much attention yet. In order to set up an LES model of good quality, two requirements are defined. The first one is to ensure that the percentage of resolved turbulent kinetic energy (TKE) is sufficient. Pope [32] defined a criterion for the minimum amount of turbulence to be resolved in an LES model. The idea of complete LES was defined by Pope [32] as an LES simulation in which solution adapted gridding is applied in order to ensure that at each time and at each location the criterion is met. The second requirement is that the time-averaged result of the LES simulation should be grid-independent. The flow fields obtained by the numerical simulations presented in this work are evaluated using both quality requirements.
In the next section of this paper, the governing equations and sub-gridscale models selected for the development of the plume model are presented, as well as the discretisation of the equations on unstructured grids. In Section 3, a test case on a vertical plume is described. In order to manage the large number of relevant processes, complexity is added in a step-wise manner: firstly crossflow is added to the vertical plume, and the results are described and analysed in Section 4. Secondly, the effect of a schematised ship hull on the plume dispersion is studied by adding a backward-facing step to the flow (Section 5). Finally, the results are discussed and conclusions are drawn in Section 6.
Section snippets
Model setup
A negatively-buoyant, vertically downward injected fine sediment plume in a horizontal crossflow is considered in this paper. The buoyancy vector and initial velocity have the same direction, pointing both in nadir direction. Apart from the particles’ settling velocity, the flow is dynamically equivalent to an upward injection of positively buoyant fluids in a horizontal crossflow such as in a smokestack plume. The negatively-buoyant sediment plume occurring when for example dredging vessels
Discussion and conclusions
A two-phase LES model with a mixture model approach for the dispersed phase is presented. The objective was to develop a numerical model capable of simulating the plume trajectory and turbulent dispersion properly in time and in three spatial dimensions. The LES model was designed to execute simulations with acceptable processing times so that in the future upscaling from laboratory to prototype scale becomes possible, as well as implementation of a realistic dredging vessel geometry. The LES
Acknowledgements
The authors wish to thank the Agentschap voor Innovatie door Wetenschap en Technologie (IWT) and International Marine and Dredging Consultants (IMDC) for funding this work. The authors are also grateful to the anonymous reviewers for their numerous instructive and constructive suggestions to improve the quality of this article.
References (60)
- et al.
Size and settling velocities of cohesive flocs and suspended sediment aggregates in a trailing suction hopper dredge plume
Cont. Shelf Res.
(2011) - et al.
PLIF and PIV measurements of the self-preserving structure of steady round buoyant turbulent plumes in crossflow
Int. J. Heat Fluid Flow
(2005) - et al.
Direct numerical simulation of pulsed jets-in-crossflow
Comput. & Fluids
(2010) - et al.
Effects of crossflow on puff and oscillation modes of a pulsed elevated transverse jet
Eur. J. Mech. B Fluids
(2012) - et al.
Large-eddy simulation of a turbulent forced plume
Eur. J. Mech. B Fluids
(2001) Large eddy simulation using an unstructured mesh based finite-volume solver
Vortex identification: New requirements and limitations
Int. J. Heat Fluid Flow
(2007)Turbulent flow over a swept backward-facing step
Eur. J. Mech. B Fluid
(2004)Environmental Aspects of Dredging
(2008)- B. Decrop, T. De Mulder, P. Troch, E. Toorman, M. Sas, 2012. Experimental investigation of negatively buoyant sediment...
Turbulent gravitational convection from maintained and instantaneous sources
Proc. R. Soc. Lond. Ser. A Math. Phys. Sci.
Mixing in Inland and Coastal Waters
Investigations of round vertical turbulent buoyant jets.
J. Fluid Mech.
Structure of round, fully developed, buoyant turbulent plumes
ASME Trans. J. Heat Transfer
A 3D LIF system for turbulent buoyant jet flows
Exp. Fluids
Experimental study of counter-rotating vortex pair trajectories induced by a round jet in cross-flow at low velocity ratios
Exp. Fluids
Vortical structure in the wake of a transverse jet
J. Fluid Mech.
Vertical mixing in the fully developed turbulent layer of sediment-laden open-channel flow
J. Hydraul. Eng.
Principles of Sediment Transport in Rivers, Estuaries and Coastal Seas
Cohesive sediment transport. I: Process description
J. Hydraul. Eng.
Modeling suspended sediment transport in nonequilibrium situations
J. Hydraul. Eng.
A proposed modification of the Germano subgrid-scale closure method
Phys. Fluids A
Buoyant jets in density-stratified crossflow
J. Hydraul. Eng.
Mixing, structure and scaling of the jet in crossflow
J. Fluid Mech.
The interaction of jets with crossflow
Annu. Rev. Fluid Mech.
Large eddy simulation of a pulsed jet in cross-flow
J. Fluid Mech.
Towards LES models of jets and plumes
Center for Turbulence Research, Annu. Res. Briefs
Large Eddy Simulation of buoyant plumes
On the reliability of large-eddy simulation for dispersed two-phase flows
Cited by (15)
Fluid-particle interaction regimes during the evolution of turbidity currents from a coupled LES/DEM model
2022, Advances in Water ResourcesCitation Excerpt :Depth-resolving models can be further subdivided into quasi-single-phase and two-phase models. In quasi-single-phase models, the particle phase is represented by a continuous particle concentration field and assumed to exhibit the same velocity as the fluid phase, or a velocity mapped to the fluid phase velocity using an empirical relationship (e.g., Decrop et al., 2015), while the mass and momentum balances of the particle-fluid mixture control the TC evolution. In two-phase models, particle and fluid phase are solved separately from one another, with separate mass and momentum conversation equations for each phase.
New plume rise modeling in a turbulent atmosphere via hybrid RANS-LES numerical simulation
2018, Journal of Wind Engineering and Industrial AerodynamicsCitation Excerpt :They validated the results with field data. Decrop et al. (2015) studied momentum-dominated and buoyancy-dominated regimes’ feature of a two-phase plume. They also modified the coefficients of the integral plume rise models using numerical simulation results.
Direct effect of atmospheric turbulence on plume rise in a neutral atmosphere
2017, Atmospheric Pollution ResearchCitation Excerpt :Pournazeri et al. (2012) modified the conservation equations to account for updrafts and downdrafts in the integral plume rise models and verified the results with field data. Decrop et al. (2015) studied momentum dominated and buoyancy dominated regimes’ feature of a two phase plume experimentally and numerically. They also modified the coefficients of integral plume rise models using numerical simulation results.