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.     

Integral Model for Turbulent Buoyant Jets in Unbounded Stratified Flows Part 2: Plane Jet Dynamics Resulting from Multiport Diffuser Jets

An integral model for the plane buoyant jet dynamics resulting from the interaction of multiple buoyant jet effluxes spaced along a diffuser line is considered as an extension of the round jet formulation that was proposed in Part I. The receiving fluid is given by an unbounded ambient environment with uniform density or stable density stratification and under stagnant or steady sheared current conditions. Applications for this situation are primarily for submerged multiport diffusers for discharges of liquid effluents into ambient water bodies, but also for multiple cooling tower plumes and building air-conditioning. The CorJet model formulation describes the conservation of mass, momentum, buoyancy and scalar quantities in the turbulent jet flow in the plane jet geometry. It employs an entrainment closure approach that distinguishes between the separate contributions of transverse shear and of internal instability mechanisms, and contains a quadratic law turbulent pressure force mechanism. But the model formulation also includes several significant three-dimensional effects that distinguish actual diffuser installations in the water environment. These relate to local merging processes from the individual multiple jets, to overall finite length effects affecting the plume geometry, and to bottom proximity effects given by a “leakage factor” that measures the combined affect of port height and spacing in allowing the ambient flow to pass through the diffuser line in order to provide sufficient entrainment flow for the mixing downstream from the diffuser. The model is validated in several stages: First, comparison with experimental data for the asymptotic, self-similar stages of plane buoyant jet flows, i.e. the plane pure jet, the pure plume, the pure wake, the advected line puff, and the advected line thermal, support the choice of the turbulent closure coefficients contained in the entrainment formulation. Second, comparison with data for many types of non-equilibrium flows with a plane geometry support the proposed functional form of the entrainment relationship, and also the role of the pressure force in the jet deflection dynamics. Third, the observed behavior of the merging process from different types of multiport diffuser discharges in both stagnant and flowing ambient conditions and with stratification appears well predicted with the CorJet formulation. Fourth, a number of spatial limits of applicability, relating to terminal layer formation in stratification or transition to passive diffusion in a turbulent ambient shear flow, have been proposed. In sum, the CorJet integral model appears to provide a mechanistically sound, accurate and reliable representation of complex buoyant jet mixing processes, provided the condition of an unbounded receiving fluid is satisfied.     

Large-eddy simulation of starting buoyant jets

A series of Large Eddy Simulations (LES) are performed to investigate the penetration of starting buoyant jets. The LES code is first validated by comparing simulation results with existing experimental data for both steady and starting pure jets and lazy plumes. The centerline decay and the growth rate of the velocity and concentration fields for steady jets and plumes, as well as the simulated transient penetration rate of a starting pure jet and a starting lazy plume, are found to compare well with the experiments. After validation, the LES code is used to study the penetration of starting buoyant jets with three different Reynolds numbers from 2000 to 3000, and with a wide range of buoyancy fluxes from pure jets to lazy plumes. The penetration rate is found to increase with an increasing buoyancy flux. It is also observed that, in the initial Period of Flow Development, the two penetrative mechanisms driven by the initial buoyancy and momentum fluxes are uncoupled; therefore the total penetration rate can be resolved as the linear addition of these two effects. A fitting equation is proposed to predict the penetration rate by combining the two independent mechanisms.     

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.     

Statistical properties of concentration fluctuations in two merging plumes

The statistics of the fluctuating concentration field within a plume is important in the analysis of atmospheric dispersion of toxic, inflammable and odorous gases. Previous work has tended to focus on concentration fluctuations in single plumes released in the surface layer or at ground level and there is a general lack of information about the mixing of two adjacent plumes and how the statistical properties of the concentration fluctuations are modified in these circumstances. In this work, data from wind tunnel experiments are used to analyse the variance, skewness, kurtosis, intermittency, probability density function and power spectrum of the concentration field during the mixing of two identical plumes and results are compared with those obtained for an equivalent single plume. The normalised variance, skewness and kurtosis on the centre-lines of the combined plume increase with distance downwind of the stack and, in the two-source configuration, takes lower values than those found in the single plumes. The results reflect the merging process at short range, which is least protracted for cases in which the sources are in-line or up to 30 \(^{\circ }\) off-line. At angles of 45 \(^{\circ }\) and more, the plumes are effectively side-by-side during the merging process and the interaction between the vortex pairs in each plume is strong. Vertical asymmetry is observed between the upper and the lower parts of the plumes, with the upper part having greater intermittency (i.e. the probability that no plume material is present) and a more pronounced tail to the concentration probability distribution. This asymmetry tends to diminish at greater distances from the source but occurs in both buoyant and neutral plumes and is believed to be associated with the ‘bending-over’ of the emission in the cross-flow and the vortex pair that this generates. The results allowed us to identify three phases in plume development. The first, very near the stack, is dominated by turbulence generated within the plume and characterised by concentration spectra with distinct peaks corresponding to scales comparable with those of the counter-rotating vortex pair. A second phase follows at somewhat greater distances downwind, in which there are significant contributions to the concentration fluctuations from both the turbulence internal to the plume and the external turbulence. The third phase is one in which the concentration fluctuations appear to be controlled by the external turbulence present in the ambient flow.     

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.     

Determination of the entrainment coefficient of a pure plume using the salt-bath technique

Environmental Fluid Mechanics - The entrainment coefficient of the pure plume, the ratio of the radial velocity of the entraining fluid at the edge of the plume to the axial velocity, is...     

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.     

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.     

Experimental Study of the Criteria of Flow Laminarization in Two-Dimensional Dense Gas Plumes

Laminarization of flow in a two-dimensional dense gas plume was experimentally investigated in this study. The plume was created by releasing CO₂ through a ground-level line source into a simulated turbulent boundary layer over an aerodynamically rough surface in a meteorological wind tunnel. The bulk Richardson number (Ri_*), based on negative plume buoyancy, plume thickness, and friction velocity, was varied over a wide range so that the effects of stable stratification on plume laminarization could be observed. A variety of ambient wind speeds as well as three different sizes of roughness arrays were used so that possible effects of roughness Reynolds number (Re_*) on plume laminarization could also be identified. Both flow visualization methods and quantitative measurements of velocity and intermittency of turbulence were used to provide quantitative assessments of plume laminarization.Flow visualization provided an overall picture of how the plume was affected by the negative buoyancy. With increasing Ri_*, both the plume depth and the vertical mixing were significantly suppressed, while upstream propagation of the plume from the source was enhanced. The most important feature of the flow revealed by visualization was the laminarization of flow in the lower part of the plume, which appeared to be closely related to both Ri_* and Re_*.Measurements within the simulated dense gas plumes revealed the influence of the stable stratification on mean velocity and turbulence intensity profiles. Both the mean velocity and turbulence intensity were significantly reduced near the surface; and these reductions systematically depended on Ri_*. The roughness Reynolds number also had considerable influence on the mean flow and turbulence structure of the dense gas plumes.An intermittency analysis technique was developed and applied to the digitized instantaneous velocity signals. It not only confirmed the general flow picture within the dense plume indicated by the flow visualization, but also clearly demonstrated the changes of flow regime with variations in Ri_* and Re_*. Most importantly, based on this intermittency analysis, simple criteria for characterizing different flow regimes are formulated; these may be useful in predicting when plume laminarization might occur.     

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.     

Evaluation of CFD Simulation using RANS Turbulence Models for Building Effects on Pollutant Dispersion

CFD evaluations were performed to examine the applicability of the RANS methods in simulating pollutant dispersion near, within and over three typical building configurations: (1) an isolated building, (2) a building array and (3) an urban intersection. The CFD results are compared with values obtained from wind tunnel tests. In some situations major differences between the wind tunnel tests and the CFD results were observed. The main source of difference between the CFD and wind tunnel results was inadequate modelling of local flow patterns using the RANS turbulence models. Also inappropriate evaluation of high intermittent turbulent mixing in the RANS approach may lead to either over-prediction or under-prediction of the concentration level, by up to a factor of 10, depending on the case investigated.     

Dredged Material Disposal at the Edge of the Florida Current

A field data collection project was undertaken to investigate the short-term fate of dredged material discharged in the designated Miami Ocean Dredged Material Disposal Site (ODMDS) before dredging of the Miami River and the Miami Harbor Turning Basin begins. the designated ODMDS is located in relatively deep water for discharge sites with a typical bottom depth of 150 metres and is also located in the western boundary region of the Gulf Stream current off Miami. Acoustical backscattering, current, particulate, temperature and salinity data were gathered over a three day period from April 24, 1990 through April 26, 1990. the major generic features of shallow-water discharge plumes were observed to be present: (a) the presence of a rapid convective descending plume portion; (b) impact of that plume portion with the ocean bottom and concomitant generation of a bottom surge; (c) rapid horizontal width growth of the descending plume through entrainment; and (d) retention of a residual plume portion within the water column. A well-mixed upper water column layer extending to a depth of 40 to 60 metres below the surface of the ocean permitted measurements of the plume entrainment coefficient free from bottom boundary, water column density gradient, and vertical current shear effects which are usually present in relatively shallow, e.g. less than 40 metres bottom depth, coastal ocean discharge studies. Entrainment coefficient estimates obtained in this study were between 0.5 to 0.7. the residual water plume material was tracked over one-half hour during each of eight discharge events and was transported in a north-northeast direction.     

Dredged Material Disposal at the Edge of the Florida Current

A field data collection project was undertaken to investigate the short-term fate of dredged material discharged in the designated Miami Ocean Dredged Material Disposal Site (ODMDS) before dredging of the Miami River and the Miami Harbor Turning Basin begins. the designated ODMDS is located in relatively deep water for discharge sites with a typical bottom depth of 150 metres and is also located in the western boundary region of the Gulf Stream current off Miami. Acoustical backscattering, current, particulate, temperature and salinity data were gathered over a three day period from April 24, 1990 through April 26, 1990. the major generic features of shallow-water discharge plumes were observed to be present: (a) the presence of a rapid convective descending plume portion; (b) impact of that plume portion with the ocean bottom and concomitant generation of a bottom surge; (c) rapid horizontal width growth of the descending plume through entrainment; and (d) retention of a residual plume portion within the water column. A well-mixed upper water column layer extending to a depth of 40 to 60 metres below the surface of the ocean permitted measurements of the plume entrainment coefficient free from bottom boundary, water column density gradient, and vertical current shear effects which are usually present in relatively shallow, e.g. less than 40 metres bottom depth, coastal ocean discharge studies. Entrainment coefficient estimates obtained in this study were between 0.5 to 0.7. the residual water plume material was tracked over one-half hour during each of eight discharge events and was transported in a north-northeast direction.     

Plume rise and spread in buoyant releases from elevated sources in the lower atmosphere

This study focuses on the influence of emission conditions—velocity and temperature—on the dynamics of a buoyant gas release in the atmosphere. The investigations are performed by means of wind tunnel experiments and numerical simulations. The aim is to evaluate the reliability of a Lagrangian code to simulate the dispersion of a plume produced by pollutant emissions influenced by thermal and inertial phenomena. This numerical code implements the coupling between a Lagrangian stochastic model and an integral plume rise model being able to estimate the centroid trajectory. We verified the accuracy of the plume rise model and we investigated the ability of two Lagrangian models to evaluate the plume spread by means of comparisons between experiments and numerical solutions. A quantitative study of the performances of the models through some suitable statistical indices is presented and critically discussed. This analysis shows that an additional spread has to be introduced in the Lagrangian trajectory equation in order to account the dynamical and thermal effects induced by the source conditions.     

Measurement of Odor-Plume Structure in a Wind Tunnel Using a Photoionization Detector and a Tracer Gas

The patterns of stimulus available to moths flying along pheromone plumes in a 3-m-long wind tunnel were characterized using a high frequency photoionization detector in conjunction with an inert tracer gas. Four contrasting flow regimes and source conditions were produced: odor released in pulses from a vertical and horizontal array of four sources, odor released continuously from a point source, and odor released continuously from a point source into an oscillatory wake. Although the four flow regimes produced plumes of intermittent and fluctuating concentration, there were considerable differences in the structure of the signal presented to the sensor. Pulses of tracer gas released at 10 Hz retained most of their longitudinal and lateral separation. The plume growing in the disturbed flow (`oscillatory'), was broader in its lateral extent than the plume growing in an undisturbed flow (`continuous'), and the concentrations in the former were lower at each downstream position. The signal recorded in the disturbed flow had higher intermittency, but the ratio between the peak concentration and the signal mean was lower than in the continuous plume. Time scales were typically longer in the tunnel than in a field setting, but length scales and the main features of intermittency and fluctuation were similar. Moths flying along plumes of pheromone in this and similar wind tunnels typically slow their velocity and narrow the lateral excursions of their flight track as they approach a pheromone source. Which features of the plumes measured in this study account for these behavioral reactions remains to be determined.     

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.     

New analytical formulations for calculation of dispersion parameters of Gaussian model using parallel CFD

New analytical formulations are presented for calculation of most effective parameters in the Gaussian plume dispersion model; the standard deviations of concentration for horizontal and vertical dispersion in neutral atmosphere conditions. Employing parallel Computational Fluid Dynamics (CFD) as a powerful tool, some well-known analytical generations of Pasquill–Gifford–Turner experimental data are modified. To achieve this aim, CFD simulations are carried out for single stack dispersion on flat terrain surface and ground level concentrations are determined in different distances. An inverse procedure in Gaussian plume dispersion model is then applied and standard deviations of horizontal and vertical dispersions are obtained. The values are compared with those of the well-known methods of Doury, Briggs and Hanna in two cases: the experimental data for release of krypton-85 from 100 m high and pollution dispersion from three 28 m high stacks of Besat power plant near Tehran. The comparison indicates that new formulations for plume dispersion are more accurate than other well-known formulations.     

Mass Transport in Vegetated Shear Flows

Submerged aquatic vegetation has the potential to greatly improve water quality through the removal of nutrients, particulates and trace metals. The efficiency of this removal depends heavily upon the rate of vertical mixing, which dictates the timescale over which these constituents remain in the canopy. Continuous dye injection experiments were conducted in a flume with model vegetation to characterize vertical mass transport in vegetated shear flows. Through the absorbance–concentration relationship of the Beer–Lambert Law, digital imaging was used to provide high-resolution concentration profiles of the dye plumes. Vertical mass transport is dominated by the coherent vortices of the vegetated shear layers. This is highlighted by the strong periodicity of the transport and its simple characterization based on properties of the shear layer. For example, the vertical turbulent diffusivity is directly proportional to the shear and thickness of the layer. The turbulent diffusivity depends upon the size of the plume, such that the rate of plume growth is lower near the source. In the far-field, mass is mixed more than twice as rapidly as momentum. Finally, plume size is dictated predominantly by X, a dimensionless distance that scales upon the number of vortex rotations experienced by the plume.     

Mixing at a sediment concentration interface in turbulent open channel flow

In this work we address the role of turbulence on mixing of clear layer of fluid with sediment-laden layer of fluid at a sediment concentration interface. This process can be conceived as the entrainment of sediment-free fluid into the sediment-laden layer, or alternatively, as the transport of sediment into the top sediment-free flow. This process is governed by four parameters—Reynolds number of the flow \(Re_\tau\), non-dimensional settling velocity of the sediment (proxy for sediment size) \(\tilde{V}\), Richardson number \(Ri_\tau\) and Schmidt number Sc. For this work we have performed direct numerical simulations for fixed Reynolds and Schmidt numbers while varying the values of Richardson number and particle settling velocity. In the simple model considered here, the flow’s momentum and turbulence pre-exists over the entire layer of fluid, while the sediment is initially confined to a layer close to the bed. Mixing of sediment-free fluid with the sediment-laden layer is associated primarily with upward transport of sediment and buoyancy. There is no simultaneous upward transport of fluid momentum and turbulence into the sediment-free fluid layer, which is already in motion and turbulent. The analysis performed shows that the ability of the flow to transport a given sediment size decreases with the distance from the bottom, and thus only fine enough sediment particles are transported across the sediment concentration interface. For these cases, the concentration profiles evolve to a final steady state in good agreement with the well-known Rouse profile. The approach towards the Rouse profile happens through a transient self-similar state. This behavior of the flow is not seen for larger particles. Detailed analysis of the three dimensional structure of the sediment concentration interface shows the mechanisms by which sediment particles are lifted up by tongues of sediment-laden fluid with positive correlation between vertical velocity and sediment concentration. Finally, the mixing ability of the flow is addressed by monitoring the time evolution of the center of mass of the sediment-laden layer and the vertical location of the sediment-free/sediment-laden interface.