Approximation of Turbulent Wall Shear Stresses in Highly Transient Pipe Flows

Theoretical predictions of wall shear stresses in unsteady turbulent flows in pipes are developed for all flow conditions from fully smooth to fully rough and for Reynolds numbers from 103 to 108. A weighting function approach is used, based on a two-region viscosity distribution in the pipe cross section that is consistent with the Colebrook–White expression for steady-state wall friction. The basic model is developed in an analytical form and the resulting weighting function is then approximated as a sum of exponentials using a modified form of an approximation due to Trikha. A straightforward method is presented for the determination of appropriate values of coefficients for any particular Reynolds number and pipe roughness ratio. The end result is a method that can be used relatively easily by analysts seeking to model unsteady flows in pipes and ducts.     

Curvilinear Flow Profiles Based on Reynolds Averaging

A new theoretical approach is presented for the derivation of free surface profiles of two-dimensional steady and unsteady flows by solving the Reynolds-averaged Navier-Stokes equations applied to the turbulent flow regime. This approach enables us to compute the steady and unsteady curvilinear flows having small curvatures, such as free overfall and constant velocity surge. In addition, the applications of the theory to the second-order waves are illustrated through the problems of small height bore and second-order tide.     

Well-Balanced Scheme between Flux and Source Terms for Computation of Shallow-Water Equations over Irregular Bathymetry

Although many numerical techniques such as approximate Riemann solvers can be used to analyze subcritical and supercritical flows modeled by hyperbolic-type shallow-water equations, there are some difficulties in practical applications due to the numerical unbalance between source and flux terms. In this study, a revised surface gradient method is proposed that balances source and flux terms. The new numerical model employs the MUSCL–Hancock scheme and the HLLC approximate Riemann solver. Several verifications are conducted, including analyses of transcritical steady-state flows, unsteady dam break flows on a wet and dry bed, and flows over an irregular bathymetry. The model consistently returns accurate and reasonable results comparable to those obtained through analytical methods and laboratory experiments. The revised surface gradient method may be a simple but robust numerical scheme appropriate for solving hyperbolic-type shallow-water equations over an irregular bathymetry.     

Modeling Depth-Averaged Velocity and Boundary Shear in Trapezoidal Channels with Secondary Flows

The Shiono and Knight method (SKM) offers a new approach to calculating the lateral distributions of depth-averaged velocity and boundary shear stress for flows in straight prismatic channels. It accounts for bed shear, lateral shear, and secondary flow effects via 3 coefficients—f,λ, and Γ—thus incorporating some key 3D flow feature into a lateral distribution model for streamwise motion. The SKM incorporates the effects of secondary flows by specifying an appropriate value for the Γ parameter depending on the sense of direction of the secondary flows, commensurate with the derivative of the term Hρ(UV)d. The values of the transverse velocities, V, have been shown to be consistent with observation. A wide range of boundary shear stress data for trapezoidal channels from different sources has been used to validate the model. The accuracy of the predictions is good, despite the simplicity of the model, although some calibration problems remain. The SKM thus offers an alternative methodology to the more traditional computational fluid dynamics (CFD) approach, giving velocities and boundary shear stress for practical problems, but at much less computational effort than CFD.     

Evaluation of Unsteady Wall Shear Stress by Zielke’s Method

The method used in the classical paper by Zielke to estimate the unsteady component of shear stress in unsteady pipe flows is revisited. It is found that the method is undesirably sensitive to the size of the integration time step. The sensitivity is shown to be caused dominantly by the first term in the integration when inadequate allowance is made for the infinite value of the weighting function. A simple method of avoiding the error without requiring the use of small grid sizes is presented.     

Transcritical Flow due to Channel Contraction

The singular point is a physical phenomenon consistent with critical flow conditions, and, consequently, the real control section of a water surface profile. A general method to study the location, type, and water surface slope of a singular point is described. This method, in addition to improving the classic gradually-varied flow theory, can be used for the design of channel contractions involving transcritical flows, as in the case of chute spillways and inlets to river diversion tunnels. The method is explained in the case of a chute spillway and verified against experimental data recorded in a Venturi channel.     

Simulation of Transcritical Flow in Pipe/Channel Networks

Using finite difference methods in conjunction with the reduced momentum equation and applying boundary condition structure inherent to subcritical flow to all regimes, is an approach that enables efficient numerical simulation of supercritical and transcritical flows in pipe/channel systems. However, as well as certain errors within a single channel due to incomplete equations, this technique also may introduce unwanted effects propagating across a network in both upstream and downstream directions. These may include: unrealistic backwater effects due to improper boundary conditions, nonamplifying oscillations due to jerky jump movement, and other computational instabilities. Practical implications of these are analyzed in detail and are illustrated using a set of examples. Sensitivity analyzes and comparisons with analytical solutions and laboratory experiments are made. The measures to reduce the inaccuracies inevitable in simulation of transcritical flows are discussed.     

Godunov-Type Solutions for Transient Flows in Sewers

This work is part of a long term project which aims at developing a hydraulic model for real-time simulation of unsteady flows in sewers ranging from gravity flows, to partly gravity–partly surcharged flows to fully surcharged flows. The success of this project hinges on the ability of the hydraulic model to handle a wide range of complex boundaries and to provide accurate solutions with the least central processing unit time. This first paper focuses on the development and assessment of two second-order explicit finite-volume Godunov-type schemes (GTS) for unsteady gravity flows in sewers, but with no surcharging. Traditionally, hydraulic transients have been modeled using the method of characteristics (MOC), which is noted for its ability to handle complex boundary conditions (BCs). The two GTS described herein incorporate BCs in a similar manner to the MOC. The accuracy and efficiency of these GTS schemes are investigated using problems whose solution contains features that are relevant to transient flows in sewers such as shock, expansion, and roll waves. The results show that these GTS schemes are significantly faster to execute than the fixed-grid MOC scheme with space-line interpolation, and in some cases, the accuracy produced by the two GTS schemes cannot be matched by the accuracy of the MOC scheme, even when a Courant number close to one and a large number of grids is used. Furthermore, unlike the MOC solutions, which exhibit increasing numerical dissipation with decreasing Courant numbers, the resolution of the shock fronts was maintained by the GTS schemes even for very low Courant numbers (0.001).     

Numerical Sensitivity Study of Unsteady Friction in Simple Systems with External Flows

This paper investigates the importance of unsteady friction effects when performing water hammer analyses for pipe systems with external fluxes due to demands, leaks, and other system elements. The transient energy equation for a system containing an orifice-type external flow is derived from the two-dimensional, axial momentum equation. A quasi-two-dimensional flow model is used to evaluate the relative energy contribution of total friction, unsteady friction, and the external flow, in a 1,500?m pipeline, with orifice flows ranging from steady-state flows of 2–70% of the mean pipe flow, and a Reynolds number of 600,000. It is found that for initial lateral flows larger than around 30% of the mean flow, unsteady friction effects can probably be neglected, whereas for external flows smaller than this, unsteady friction should generally be considered. Overall, the relative role of unsteady friction is found to diminish as the external flux increases, implying that unsteady friction is not critical for systems with large external flows. These results imply that unsteady friction may have a significant impact on the validity of transient leak detection techniques that have been derived assuming quasi-steady friction. To demonstrate this point, an existing transient leak detection method, originally derived under quasi-steady conditions, is tested with unsteady friction included.     

Hydraulic Elements Chart for Pipe Flow Using New Definition of Hydraulic Radius

The Manning formula is used routinely to calculate the mean velocity of uniform flow. Although this empirical formula is effective when applied to uniform flow in simple rectangular or trapezoidal cross sections, the roughness coefficient of the formula is variable when examining flow in a pipe that is partially full. Thus, the coefficient must be altered depending on the relative depth of fluid in the pipe. As this seems to be due to the definition of the hydraulic radius, a new definition of hydraulic radius is proposed here that was used to calculate a hydraulic elements chart for flow in pipes with a constant roughness coefficient. The results of the calculations showed very good agreement with Camp’s chart. Furthermore, with adjustment of the “free-surface weight factor,” this method was also capable of expressing other hydraulic elements charts reported previously. This new definition of hydraulic radius can also be applied to flow in simple cross sections and may be developed further for use with compound channel flows in future studies.     

Junction and Drop-Shaft Boundary Conditions for Modeling Free-Surface, Pressurized, and Mixed Free-Surface Pressurized Transient Flows

A junction and drop-shaft boundary conditions (BCs) for one-dimensional modeling of transient flows in single-phase conditions (pure liquid) are formulated, implemented and their accuracy are evaluated using two computational fluid dynamics (CFD) models. The BCs are formulated in the case when mixed flows are simulated using two sets of governing equations, the Saint-Venant equations for the free-surface regions and the compressible water hammer equations for the pressurized regions. The proposed BCs handle all possible flow regimes and their combinations. The flow in each pipe can range from free surface to pressurized flow and the water depth at the junction or drop shaft can take on all possible levels. The BCs are applied to the following three cases: (1) a three-way merging flow; (2) a three-way dividing flow; and (3) a drop shaft connected to a single-horizontal pipe subjected to a rapid variation of the water surface level in the drop shaft. The flow regime for the first two cases range from free surface to pressurized flows, while for the third case, the flow regime is pure pressurized flow. For the third case, laboratory results as well as CFD results were used for evaluating its accuracy. The results suggest that the junction and drop-shaft BCs can be used for modeling transient free-surface, pressurized, and mixed flow conditions with good accuracy.     

Ultrasonic Lift Measurement Technique for Flow-Induced Structural Vibrations

Deficiencies exist in our current ability to measure lift forces in wind-tunnel experiments on vibrating structures in a fluid flow. The ultrasonic lift measurement (ULM) technique has been previously developed to measure time-averaged fluid circulation and lift on stationary structures. The ULM technique is based on measuring transit times of acoustic pulses along paths enclosing the structure. A quasi-steady method based on the Kutta–Joukowski theorem has been used in the past to convert fluid circulation to lift values in ULM studies. In this paper, the largely unstudied extension of the ULM technique to measure unsteady lift forces in flows involving structural vibration is considered. Analytic methods are developed that can be used to properly convert the instantaneous circulation measurements (attainable from ULM experiments) to lift values. These unsteady methods are validated by using numerical simulations of flows over flat plates undergoing oscillating motion. It is shown that the addition of unsteady terms provides a method with improved accuracy over the previous quasi-steady assumption. The methods are also applied to ULM data from an oscillating airfoil experiment in a preliminary study.     

Noise of Acoustic Doppler Velocimeter Data in Bubbly Flows

Acoustic Doppler velocimeter (ADV) measurements are useful and powerful for measurements of mean and turbulent components of fluids in both hydraulic experimental facilities and fields. However, it is difficult to use the ADV in bubbly flows because air bubbles generate spike noise in the ADV velocity data. This study describes the validity of the ADV measurements in bubbly flows. The true three-dimensional phase space method originally developed by Goring and Nikora in 2002 is significantly useful to eliminate spike noise of ADV recorded data in bubbly flow in comparison with the classical low correlation method.     

Role of Energy Loss on Discharge Characteristics of Sluice Gates

A theoretical method was used to derive an equation for the discharge coefficient of sluice gates in rectangular channels under orifice-flow (both free and submerged) conditions. The proposed equation allows for the effects of energy dissipation between the upstream section of the gate and the vena contracta. The hydraulic energy loss in the upstream pool is attributable to the induced turbulence by the recirculating region and to the growth of bottom boundary layer. For the submerged-flow condition, turbulent shear-layer entrainment is also responsible for the energy loss. This energy loss is introduced into the equations through a coefficient k that has been conventionally assumed to be negligible. Experimental data from the literature were used to validate the equation, which showed good agreement with the measured values. It is also shown that the magnitude of the energy-loss factor is a function of the geometry of the gate and can modify the discharge coefficient. An equation for the distinguishing condition between free and submerged flows is also presented. The new equations can be used to predict the performance of sluice gates with different edge shapes under free- and submerged-flow situations.     

Formulas for Friction Factor in Transitional Regimes

This note concerns variations of the friction factor in the two transitional regimes, one between laminar and turbulent flows and the other between fully smooth and fully rough turbulent flows. An interpolation approach is developed to derive a single explicit formula for computing the friction factor in all flow regimes. The results obtained for pipe flows give a better representation of Nikuradse’s experimental data, in comparison with other implicit formulas available in the literature. Certain modifications are also made for applying the obtained friction formula to open-channel flows.     

Hydraulic Models of the Flow Distribution in a Four Branch Open Channel Junction with Supercritical Flow

Intense rainfall on urban areas can generate severe flooding in the city, and if the conditions are right, the flow in the streets can be supercritical. The redistribution of the flow in street intersections determines the flow rates and water levels in the street network. We have investigated the flow that occurs when two supercritical flows collide in a 90° junction formed by streets of identical cross section. Several flow configurations within the intersection are possible, depending on the position of the hydraulic jumps that form in and upstream of the intersection. Previous work has identified three flow types, with Type II flows being further classified into three subregimes. Hydraulic models have been developed, based on the principles of the conservation of flow and momentum flux in the intersection, which predict the angles at which the jumps will form. These models can be used to determine the flow type that will occur. Moreover, additional models have been developed for computing the outflow discharge distribution. For Type I flows, it has not been possible to develop such a hydraulic model for the discharge distribution, but some data are provided for one configuration to indicate the influence of different parameters. For Type II and Type III flows, such models are developed, and their predictions agree with data obtained from the channel intersection facility at the Laboratory of Fluid Mechanics and Acoustics in Lyon.     

Double-Averaged Open-Channel Flows with Small Relative Submergence

We investigate the turbulent structure of shallow open channel flows where the flow depth is too small (compared with the roughness height) to form a logarithmic layer but large enough to develop an outer layer where the flow is not directly influenced by the roughness elements. Since the log layer is not present, the displacement height d, which defines the position of the zero plane, and the shear velocity u* cannot be found by fitting the velocity data to the log law. However, these parameters are still very important because they are used for scaling flow statistics for the outer and roughness layers. In this paper we propose an alternative procedure for evaluating d in laboratory conditions, where d is found from additional experiments with the fully developed log layer. We also point out the appropriate procedure for evaluating the shear velocity u* for flows with low submergence. These procedures are applied to our own laboratory flume experiments with uniform sphere roughness, where velocities were measured using Particle Image Velocimetry. Results were interpreted within the framework of the double-averaged Navier–Stokes equations and include mean velocities, turbulence intensities, Reynolds stresses, and form-induced normal and shear stresses. The data collapse well and show that in flows without a developed log layer the structure of turbulence in the outer layer remains similar to that of flows with a log layer. This means that even though the roughness layer in the experiments reported herein was sufficiently high to prevent the development of the log layer, influence of the bed roughness did not spread further up into the outer layer. Furthermore, the results show that flow statistics do not depend on relative submergence except for the form-induced stresses which increase when relative submergence decreases.     

Extensive Development of Leak Detection Algorithm by Impulse Response Method

The oscillatory flows in pipeline systems due to excitation by valve operation are efficiently analyzed by the impulse response method. The impact of leakage is incorporated into the transfer functions of the complex head and discharge. Frequency-dependent friction is used to consider the impact of unsteady friction for laminar condition. Extensive development of the impulse response method was made by considering the sources of friction associated with the local and convective acceleration of velocity for turbulent flow. The genetic algorithm was integrated into the impulse response method to calibrate the location and the quantity of leakage. The calibration function for leakage detection can be made using the pressure-head response at the valve, or the pressure-head and flow response at the section upstream from the valve. The proposed leak detection algorithm shows the potentials for being applied to a simple pipeline system with a single leak or multiple leaks.     

Two-Dimensional Depth-Averaged Flow Modeling with an Unstructured Hybrid Mesh

An unstructured hybrid mesh numerical method is developed to simulate open channel flows. The method is applicable to arbitrarily shaped mesh cells and offers a framework to unify many mesh topologies into a single formulation. A finite-volume discretization is applied to the two-dimensional depth-averaged equations such that mass conservation is satisfied both locally and globally. An automatic wetting-drying procedure is incorporated in conjunction with a segregated solution procedure that chooses the water surface elevation as the main variable. The method is applicable to both steady and unsteady flows and covers the entire flow range: subcritical, transcritical, and supercritical. The proposed numerical method is well suited to natural river flows with a combination of main channels, side channels, bars, floodplains, and in-stream structures. Technical details of the method are presented, verification studies are performed using a number of simple flows, and a practical natural river is modeled to illustrate issues of calibration and validation.     

Three-Dimensional Numerical Modeling of Initial Mixing of Thermal Discharges at Real-Life Configurations

A three-dimensional Reynolds-averaged Navier–Stokes computational fluid dynamics (CFD) model is developed for simulating initial mixing in the near field of thermal discharges at real-life geometrical configurations. The domain decomposition method with multilevel embedded overset grids is employed to handle the complexity of real-life diffusers as well as to efficiently account for the large disparity in length scales arising from the relative size of the ambient river reach and the typical diffuser diameter. An algebraic mixing length model with a Richardson-number correction for buoyancy effects is used for the turbulence closure. The governing equations are solved with a second-order-accurate, finite-volume, artificial compressibility method. The model is validated by applying it to simulate thermally stratified shear flows and negatively buoyant wall jet flows and the computed results are shown to be in good overall agreement with the experimental measurements. To demonstrate the potential of the numerical model as a powerful engineering simulation tool we apply it to simulate turbulent initial mixing of thermal discharges loaded from both single-port and multiport diffusers in a prismatic channel and a natural river. Comparisons of the CFD model results with those obtained by applying two widely used empirical mixing zone models show that the results are very similar in terms of both the rate of dilution and overall shape of the plumes. The CFD model further resolves the complex three-dimensional features of such flows, including the complex interplay of the ambient flow and thermal discharges as well as the interaction between each of discharges loaded from multiple ports, which are obviously not accessible by the simpler empirical models.