On the front conditions for gravity currents in channels of general cross-section

We consider the propagation of a high-Reynolds-number gravity current in a horizontal channel with general cross-section whose width is \(f(z), 0 \le z\le H\), and the gravity acceleration g acts in \(-z\) direction. (The classical rectangular cross-section is covered by the particular case \(f(z) =\) const.) We assume a two-layer system of homogeneous fluids of constant densities \(\rho _{c}\) (current, of height \(h < H \)) and smaller \(\rho _{a}\) (ambient, filling the remaining part of the channel). We focus attention on the calculation and assessment of the nose Froude-number condition \(Fr = U/(g' h)^{1/2}\); here U is the speed of propagation of the current and \(g' = (\rho _{c}/\rho _{a}-1) g\) is the reduced gravity. We first revisit the steady-state current, and derive compact insightful expressions of Fr and energy dissipation as a function of \(\varphi \) (\(=\) area fraction occupied by the current in the cross-section). We show that the head loss \(\delta _0\) on the stagnation line is formally a degree of freedom in the determination of \(Fr(\varphi )\), and we clarify the strong connections with the head loss \(\delta \) in the ambient fluid, and with the overall rate of dissipation \(\dot{{\mathcal{D}}}\). We demonstrate that the closure \(\delta _0 = 0\) [suggested by Benjamin (J Fluid Mech 31, 209–248, 1968) for the rectangular cross-section] produces in general the smallest Fr for a given \(\varphi \); the results are valid for a significant range \([0, \varphi _{\max }]\), in which the current is dissipative, except for the point \(\varphi _{\max }\) where \(\delta = \dot{{\mathcal{D}}} = 0\). We show that imposing the closure \(\delta = \dot{{\mathcal{D}}} = 0\), which corresponds to an energy-conserving or non-dissipative current, produces in general unacceptable restrictions of the range of validity, and large values of Fr; in particular, deep currents (\(\varphi < 0.3\) say) must be excluded because they are inherently dissipative. On the other hand, the compromise closure \(\delta (\varphi ) =\delta _0(\varphi )\) produces the simple \(Fr(\varphi ) = \sqrt{2}(1 - \varphi )\) formula whose values and dissipation properties are very close, and the range of validity is identical, to these obtained with Benjamin’s closure (moreover, we show that this corresponds to circulation-conservation solutions). The results are illustrated for practical cross-section geometries (rectangle, \(\Delta \) and \(\nabla \) triangle, circle, and the general power-law \(f(z) = b z ^\alpha \) (\(b>0, \alpha \ge 0, 0< z \le H\)). Next, we investigate the connection of the steady-state results with the time-dependent current, and show that in a lock-released current the rate of dissipation of the system is equal to, or larger than, that obtained for Fr corresponding to the conditions at the nose of the current. The results and insights of this study cover a wide range of cross-section geometry and apply to both Boussinesq and non-Boussinesq systems; they reveal a remarkable robustness of Fr as a function of \(\varphi \).     

Penetrative convection in slender containers

An experimental study was conducted to investigate the penetration of a convective mixed layer into an overlying stably (solutally) stratified layer contained in a narrow, tall vessel when the fluid is subjected to a destabilizing heat flux from below. The interest was the evolution of the bottom mixed-layer height (\(h\)) with time (\(t\)) in the presence of side-wall effects, but without the formation of conventional double-diffusive layers. The side-wall effects are expected at small mixed-layer aspect ratios, \(\varGamma_{h} = (W/h)\), where \(W\) is the container width. This case has not been studied hitherto, although there are important environmental and industrial applications. The mixed-layer growth laws for low aspect ratio convection were formulated by assuming a balance between the vertical kinetic energy flux at the interface and the rate of change of potential energy of the fluid system due to turbulent entrainment. The effects of sidewalls were considered using similarity arguments, by taking characteristic rms velocities to be a function of \(\varGamma_{h}\), in addition to buoyancy flux (\(q_{0}\)) and \(h\). In all stages of evolution, the similarity variables \(\xi = h/W\) and \(t^{\prime } = Nt/A\), where \(A = N^{3} W^{2} /4q_{0}\) and \(N\) is the buoyancy frequency, scaled the mixed-layer evolution data remarkably well. Significant wall effects were noted when \(\varGamma_{h} < 1\), and for this case the interfacial vertical turbulent velocity and length scales were identified via scaling arguments and experimental data.     

Experiments on gravity currents propagating on unbounded uniform slopes

Gravity currents propagating on \(12^\circ \), \(9^\circ \), \(6^\circ \), \(3^\circ \) unbounded uniform slopes and on an unbounded horizontal boundary are reported. Results show that there are two stages of the deceleration phase. In the early stage of the deceleration phase, the front location history follows \({(x_f+x_0)}^2 = {(K_I B)}^{1/2} (t+t_{I})\), where \((x_f+x_0)\) is the front location measured from the virtual origin, \(K_I\) an experimental constant, B the total buoyancy, t time and \(t_I\) the t-intercept. In the late stage of the deceleration phase for the gravity currents on \(12^\circ \), \(9^\circ \), \(6^\circ \) unbounded uniform slopes, the front location history follows \({(x_f+x_0)}^{8/3} = K_{VS} {{B}^{2/3} V^{2/9}_0 }{\nu }^{-1/3} ({t+t_{VS}})\), where \(K_{VS}\) is an experimental constant, \(V_0\) the initial volume of heavy fluid, \(\nu \) the kinematic viscosity and \(t_{VS}\) the t-intercept. In the late stage of the deceleration phase for the gravity currents on a \(3^\circ \) unbounded uniform slope and on an unbounded horizontal boundary, the front location history follows \({(x_f+x_0)}^{4} = K_{VM} {{B}^{2/3} V^{2/3}_0 }{\nu }^{-1/3} ({t+t_{VM}})\), where \(K_{VM}\) is an experimental constant and \(t_{VM}\) the t-intercept. Two qualitatively different flow morphologies are identified in the late stage of the deceleration phase. For the gravity currents on \(12^\circ \), \(9^\circ \), \(6^\circ \) unbounded uniform slopes, an ‘active’ head separates from the body of the current. For the gravity currents on a \(3^\circ \) unbounded uniform slope and on an unbounded horizontal boundary, the gravity currents maintain an integrated shape throughout the motion. Results indicate two possible routes to the final stage of the gravity currents on unbounded uniform slopes.

     

Identifying demand parameters in the presence of unobservables: A combined revealed and stated preference approach

• ALTHOUGH MOST ACADEMIC programs use hospital-based sites for perioperative clinical experiences, the exponential growth in the number of ambulatory surgery centers (ASCs) and the increased variety of procedures performed in these settings make ASCs excellent venues for student learning experiences.• NURSING STUDENTS in one program in Florida have the opportunity to take a perioperative practicum designed to allow them to learn about all aspects of patient care and administrative functions in ASCs.• THE ASC PRACTICUM EXPERIENCE is popular with students and staff members alike, and ASC administrators are seeking ways to increase the number of students in the program in an effort to attract new graduates to ASCs. AORN J 84 (August 2006) 219–232. © AORN, Inc, 2006.

Thin-layer models for gravity currents in channels of general cross-section area,a review

We present a brief review of the recent investigations on gravity currents in horizontal channels with non-rectangular cross-section area (such as triangle, \(\bigvee \)-valley, circle/semi-circle, trapezoid) which occur in nature (e.g., rivers) and constructed environment (tunnels, reservoirs, canals). To be specific, we discuss the propagation of a gravity current (GC) in a horizontal channel along the horizontal coordinate x, with gravity g acting in the \(-z\) direction, and y the horizontal–lateral coordinate. The bottom and top of the channel are at \(z=0,H\). The “standard” problem is concerned with 2D flow in a channel with rectangular (or laterally unbounded) cross-section area (CSA). Recent investigations have successfully extended the standard knowledge to the channels of CSA given by the quite general \(-f_1(z)\le y \le f_2(z)\) for \(0 \le z \le H\). This includes the practical \(\bigvee \)-valley, triangle, circle/semi-circle and trapezoid; these geometries may be in “up” or “down” setting with respect to gravity, e.g., \(\bigtriangleup \) and \(\bigtriangledown \). The major objective of the extended theory is to predict the height of the interface \(z=h(x,t)\) and the velocity (averaged over the CSA) u(x, t), where t is time; the prediction includes the speed and position of the nose \(u_N(t), x_N(t)\). We show that the motion is governed by a set of simplified equations, called “model,” that provides versatile and insightful solutions and trends. The emphasis in on a high-Reynolds-number current whose motion is dominated by buoyancy–inertia balance; in particular a GC released from a lock, which also contains general effects such as front and internal jumps (shocks), and reflected bore. We discuss two-layer, one-layer, and box models; Boussinesq and non-Boussinesq systems; compositional and particle-driven cases; and the effect of stratification of the ambient fluid. The models are self-contained, and admit realistic initial and boundary conditions. The governing equations are amenable to analytical solutions in some special circumstances. Some salient features of the buoyancy-viscous regime, and the estimate for the length at which transition to this regime takes place, are also presented. Some experimental support to the theory, and open questions for further investigations, are also mentioned. The major conclusions are (1) The CSA geometry has significant influence on the motion of the GC; and (2) The new theory is a useful, very significant, extension of the standard two-dimensional GC problem. The standard current is just a particular case, \(f_{1,2} =\) constants, among many other covered by the new theory.     

Variations of bed elevations due to turbulence around submerged cylinder in sand beds

This paper presents the spatio-temporal variations in bed elevations and the near-bed turbulence statistics over the deformed bed generated around the submerged cylindrical piers embedded vertically on loose sediment bed at a constant flow discharge. Experiments were carried out in a laboratory flume for three blockage ratios in the range of 0.04–0.06 using three different sizes of submerged cylinders individually placed vertically at the centerline of the flume. Clear-water experimental conditions were maintained over the smooth sediment bed surface with a constant flow discharge (\(Q = 0.015\,{\rm m}^3/{\rm sec}\)), thereby giving three different cylinder Reynolds numbers \(Re_{D_c} = \frac{U_mD_c}{\nu }\) (=10200, 12750, 15300) away from the cylinder locations, where \(U_m\) is the maximum mean velocity, \(D_c\) is the cylinder diameter and \(\nu\) is the kinematic viscosity of fluid. Instantaneous sand bed elevations around the cylinders were recorded using a SeaTek 5MHz ultrasonic ranging system of net 24 transducers to estimate bed form migration, and the near-bed velocity data at transducer locations over the stable deformed bed around the pier-like structures were collected using down-looking three-dimensional (3D) Micro-acoustic Doppler velocimeter to estimate the bottom Reynolds shear stresses and the contributions of bursting events to the dominant shear stress component. The flow perturbation generated due to relatively lower flow blockage ratio favored to achieve the stable bed condition more rapidly than the others, and larger upstream scour-depth and deformed areas were noticed for greater flow blockage ratio due to larger cylinder diameter. For larger blockage ratio in the upstream of scour-hole near the bed, occurrences of probabilities of both boundary-ward interactions (Q1 and Q3) were the dominant; whereas in the downstream of the scoured region, occurrences of probabilities of second and third quadrant events (Q2 and Q4) were dominant. On the other hand, for the lower blockage ratio, quadrant (Q2) was dominant over Q4 in the downstream of scour-hole, and in the upstream of scour-hole, quadrant Q4 was the dominant.     

Numerical simulations of intrusive gravity currents interacting with a bottom-mounted obstacle in a continuously stratified ambient

In this study, the flow dynamics of intrusive gravity currents past a bottom-mounted obstacle were investigated using highly resolved numerical simulations. The propagation dynamics of a classic intrusive gravity current was first simulated in order to validate the numerical model with previous laboratory experiments. A bottom-mounted obstacle with a varying non-dimensional height of \(\tilde{D}=D/H\), where D is the obstacle height and H is the total flow depth, was then added to the problem in order to study the downstream flow pattern of the intrusive gravity current. For short obstacles, the intrusion re-established itself downstream without much distortion. However, for tall obstacles, the downstream flow was found to be a joint effect of horizontal advection, overshoot-springback phenomenon, and associated Kelvin-Helmholtz instabilities. Analysis of the numerical results show that the relationship between the downstream propagation speed and the obstacle height can be subdivided into three regimes: (1) a retarding regime (\(\tilde{D}\) \(\approx \) 0–0.3) where a 30 % increase in obstacle height leads to a 20 % reduction in propagation speed, simply due to the obstacle’s retarding effect; (2) an impounding regime (\(\tilde{D}\) \(\approx \) 0.3–0.6) where the additional 30 % increase in obstacle height only leads to a further (negligible) 5 % reduction in propagation speed, due to the accelerating effect of upstream impoundment and downstream enhanced mixing; and (3) a choking regime (\(\tilde{D}\) \(\approx \) 0.6–1.0) where the propagation speed is dramatically reduced due to the dominance of the obstacle’s blocking effect. The obstacle thickness was found to be irrelevant in determining the downstream propagation speed at least for the parameter range explored in this study. The present work highlights the significance of topographic effects in stratified flows with horizontal pressure forcing.     

Flow and drag force around a free surface piercing cylinder for environmental applications

This paper investigates flows around a free surface piercing cylinder with Froude number F > 0.5 and Reynolds number around Re = 50,000. The aim of this work is to gain a better understanding of the flow behaviour in environmental systems such as fishways. The advances are based upon experimental and numerical results. Several flow discharges and slopes are tested to obtain both subcritical and supercritical flows. The drag force exerted on the cylinder is measured with the help of a torque gauge while the velocity field is obtained using particle velocimetry. For the numerical part, two URANS turbulence models are tested, the k-\(\omega\) SST and the RNG k-\(\varepsilon\) models using the OpenFOAM software suite for subcritical cases, and then compared with the corresponding experimental results. With fishways applications in mind, the changes in drag coefficient \(C_d\) versus Froude number and water depth are studied and experimental correlations proposed. We conclude that the most suitable URANS turbulence model for reproducing this kind of flow is the k-\(\omega\) SST model.     

Drag coefficient for rigid vegetation in subcritical open-channel flow

Drag coefficient has been commonly used as a quantifying parameter to represent the vegetative drag, i.e., resistance to the flow by vegetation. In this study, the measured data on the drag coefficient for rigid vegetation in subcritical open-channel flow reported in previous studies are collected and preprocessed for multi-parameter analysis. The effect of Froude number (Fr) on the drag coefficient for rigid vegetation in subcritical flow cannot be ignored, especially when \(Fr < 0.12\). The drag coefficient is observed to exponentially decrease with the stem Reynolds number (R _d ) and logarithmically decreased with the vegetation density (λ) when \(0.012 < \lambda < 0.12\). The relative submergence (h ^* ) has a significant effect on the drag coefficient, and a positive logarithmic relationship is summarized. A simplified three-stage empirical formula is obtained based on the divisions of Fr. Laboratory tests (with \(Fr < 0.02\)) prove that the present empirical model has higher precision compared with existing models.     

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.     

Passive scalar roughness lengths for atmospheric boundary layer flow over complex, fractal topographies

When modeling atmospheric boundary layer flow over rough landscapes, surface fluxes of flow quantities (momentum, temperature, etc.) can be described with equilibrium logarithmic law expressions, all of which require specification of a roughness length that is, physically, the elevation at which the flow quantity equals its surface value. In high Reynolds number flows, such as the atmospheric boundary layer, inertial forces associated with turbulent eddy motions are responsible for surface momentum fluxes (form, or pressure drag). Surface scalar fluxes, on the other hand, occur exclusively via diffusion in the immediate vicinity of the topography—the interfacial region—before being advected by turbulent eddy motions into the bulk of the flow. Owing to this difference in surface transfer mechanism, the passive scalar roughness length, $z_{0S}$ , is known to be less than the momentum roughness length, $z_0$ . In this work, classical relations are used to specify $z_{0S}$ during large-eddy simulation of atmospheric boundary layer flow over aerodynamically rough, synthetic, fractal topographies which exhibit power-law height energy spectrum, $E_h (k) \sim k^{\beta _s}$ , where $\beta _s$ is a (predefined) spectral exponent. These topographies are convenient since they resemble natural landscapes and $\beta _s$ can be varied to change the topography’s aerodynamic roughness (the study considers a suite of topographies with $-2.4 \le \beta _s \le -1.2$ , where $-2.4$ and $-1.2$ are the “most smooth” and “most rough” cases, respectively, corresponding with roughness Reynolds number, $Re_0 \approx 10$ and $300$ ). It is often assumed that $z_{0S}/z_{0} \approx 10^{-1}$ for all $Re_0$ . But results from this work show that the roughness length ratio, $z_{0S}/z_{0}$ , depends strongly on $Re_0$ , ranging between $10^{-3}$ and $10^{-1}$ .     

Fitting the truncated negative binomial distribution to count data

Modeling empirical distributions of repeated counts with parametric probability distributions is a frequent problem when studying species abundance. One must choose a family of distributions which is flexible enough to take into account very diverse patterns and possess parameters with clear biological/ecological interpretations. The negative binomial distribution fulfills these criteria and was selected for modeling counts of marine fish and invertebrates. This distribution depends on a vector \(\left( K,\mathfrak {P}\right) \) of parameters, and ranges from the Poisson distribution (when \(K\rightarrow +\infty \)) to Fisher’s log-series, when \(K\rightarrow 0\). Moreover, these parameters have biological/ecological interpretations which are detailed in the literature and in this study. We compared three estimators of K, \(\mathfrak {P}\) and the parameter \(\alpha \) of Fisher’s log-series, following the work of Rao CR (Statistical ecology. Pennsylvania State University Press, University Park, 1971) on a three-parameter unstandardized variant of the negative binomial distribution. We further investigated the coherence underlying parameter values resulting from the different estimators, using both real count data collected in the Mauritanian Exclusive Economic Zone (MEEZ) during the period 1987–2010 and realistic simulations of these data. In the case of the MEEZ, we first built homogeneous lists of counts (replicates), by gathering observations of each species with respect to “typical environments” obtained by clustering the sampled stations. The best estimation of \(\left( K,\mathfrak {P}\right) \) was generally obtained by penalized minimum Hellinger distance estimation. Interestingly, the parameters of most of the correctly sampled species seem compatible with the classical birth-and-dead model of population growth with immigration by Kendall (Biometrika 35:6–15, 1948).     

Throughflow analysis: A stochastic approach

Ecological network analysis (ENA), predicated on systems theory and Leontiev input–output analysis, is a method widely used in ecology to reveal ecosystem properties. An important ecosystem property computed in ENA is throughflows, the amount of matter/energy leaving each compartment of the ecosystem. Throughflows are analyzed via a matrix representing their relationships to the driving input at the boundary. Network particle tracking (NPT) builds on ENA to offer a Lagrangian particle method that describes the activity of the ecosystem at the microscopic level. This paper introduces a Lagrangian throughflow analysis methodology using NPT and shows that the NPT throughflow matrix, , agrees with the conventional ENA throughflow matrix, , for ecosystems at steady-state with donor-controlled flows. The matrix is computed solely from the pathways (particles’ histories) generated by NPT simulations and its average over multiple runs of the algorithm with longer simulation time agrees with the Eulerian matrix (Law of Large Numbers). While the traditional NEA throughflow analysis is mostly used with steady-state ecosystem models, the Lagrangian throughflow analysis that we propose can be used with non-steady-state models and paves the way for the development of dynamic throughflow analysis.     

Evaluation of existing equations for temporal scour depth around circular bridge piers

Scour is defined as the processes of removal of sediment particles from water stream bed by the erosive action of activated water, and also carries sediment away from the hydraulic structures. Scour is the main cause of pier failure. Numerous equations are available for estimating temporal and equilibrium scour depth. The present study describes the phenomenon of temporal scour depth variation at bridge piers and deals with the methods for its estimation. The accuracy of six temporal scour depth equations are also checked in this study. After graphical and statistical analysis, it was found that the relationship proposed by Oliveto and Hager (J Hydraul Eng (ASCE) 128(9):811–520, 2002) predicts temporal scour depth better than other equations. Three equations of equilibrium time of scour are also used for computing equilibrium time. Equilibrium time equation proposed by Choi and Choi (Water Environ J 30(1–2):14–21, 2016) gives better agreements with observed values.