Efficient Approximation of Unsteady Friction Weighting Functions

A simple method is presented for evaluating wall shear stresses from known flow histories in unsteady pipe flows. The method builds on previous work by Trikha, but has two important differences. One of these enables the method to be used with much larger integration time steps than are acceptable with Trikha’s method. The other, a general procedure for determining approximations to weighting functions, enables it to be used at indefinitely small times (high frequencies). The method is applicable to both laminar and turbulent flows.     

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.     

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.     

Timescale Behavior of the Wall Shear Stress in Unsteady Laminar Pipe Flows

Based on two-dimensional (2D) flow model simulations, the effects of the radial structure of the flow (e.g., the nonuniformity of the velocity profile) on the pipe wall shear stress, τw, are determined in terms of bulk parameters such as to allow improved 1D modeling of unsteady contribution of τw. An unsteady generalization, for both laminar and turbulent flows, of the quasi-stationary relationship between τw and the friction slope, J, decomposes the additional unsteady contribution into an instantaneous energy dissipation term and an inertial term (that is, based on the local average acceleration-deceleration effects). The relative importance of these two effects is investigated in a transient laminar flow and an analysis of the range of applicability of this kind of approach of representing unsteady friction is presented. Finally, the relation between the additional inertial term and Boussinesq momentum coefficient, is clarified. Although laminar pipe flows are a special case in engineering practice, solutions in this flow regime can provide some insight into the behavior of the transient wall shear stress, and serve as a preliminary step to the solutions of unsteady turbulent pipe flows.     

Time-Line Interpolation Errors in Pipe Networks

An exact method of assessing numerical errors in analyses of unsteady flows in pipe networks is introduced. The assessment is valid for fixed-grid method of characteristics analyses using time-line interpolations. A pipe polynomial transfer matrix is developed and is analogous to transfer function matrices used in free oscillation theory. The influence of reachback is assessed by comparing exact numerical predictions using a polynomial transfer matrix with exact analytical predictions obtained using free oscillation theory. The investigation is part of a long-term project aimed at automating the selection of numerical grid sizes in unsteady flow analyses. The eventual goal is to enable users of unsteady flow software to prescribe required degrees of accuracy instead of specifying the numerical grid itself. This paper is only a first step toward the long-term aim, but it is a big step toward an intermediate objective of providing exact benchmarking data for the assessment of approximate methods of automatic grid selection.     

Numerical Error in Weighting Function-Based Unsteady Friction Models for Pipe Transients

The accurate simulation of pressure transients in pipelines and pipe networks is becoming increasingly important in water engineering. Applications such as inverse transient analysis for condition assessment, leak detection, and pipe roughness calibration require accurate modeling of transients for longer simulation periods that, in many situations, requires improved modeling of unsteady frictional behavior. In addition, the numerical algorithm used for unsteady friction should be highly efficient, as inverse analysis requires the transient model to be run many times. A popular model of unsteady friction that is applicable to a short-duration transient event type is the weighting function-based type, as first derived by Zielke in 1968. Approximation of the weighting function with a sum of exponential terms allows for a considerable increase in computation speed using recursive algorithms. A neglected topic in the application of such models is evaluation of numerical error. This paper presents a discussion and quantification of the numerical errors that occur when using weighting function-based models for the simulation of unsteady friction in pipe transients. Comparisons of numerical error arising from approximations are made in the Fourier domain where exact solutions can be determined. Additionally, the relative importance of error in unsteady friction modeling and unsteady friction itself in the context of general simulation is discussed.     

Semianalytical Model for Solid Concentration in Turbulent Pipe Flows with Dispersed Particles

This paper presents a semianalytical model for the radial distribution of the solid concentration in a fully developed vertical turbulent pipe two-phase flow. A simplified momentum equation in the radial direction for solid phase in a two-phase flow with dilute suspended particles was first obtained. A linear empirical closure relation for the mean gas and solid velocities along the pipe direction was constructed using published experimental data. By incorporating the closure relation, an analytical solution to the simplified solid momentum equation with the appropriate boundary conditions at the pipe center and wall was obtained. The results from this semianalytical model are able to describe the core-annulus phenomenon commonly occurring in two-phase turbulent pipe flows. Very good agreements were found between the model predictions and published experimental data.     

Efficient Treatment of the Vardy–Brown Unsteady Shear in Pipe Transients

An accurate, simple, and efficient approximation to the Vardy–Brown unsteady friction equation is derived and shown to be easily implemented within a one-dimensional characteristics solution for unsteady pipe flow. For comparison, the exact Vardy–Brown unsteady friction equation is used to model shear stresses in transient turbulent pipe flows and the resulting waterhammer equations are solved by the method of characteristics. The approximate Vardy–Brown model is more computationally efficient (i.e., requires one-sixth the execution time and much less memory storage) than the exact Vardy–Brown model. Both models are compared with measured data from different research groups and with numerical data produced by a two-dimensional turbulence waterhammer model. The results show that the exact Vardy–Brown model and the approximate Vardy–Brown model are in good agreement with both laboratory and numerical experiments over a wide range of Reynolds number and wave frequencies. The proposed approximate model only requires the storage of flow variables from a single time step while the exact Vardy–Brown model requires the storage of flow variables at all previous time steps and the two-dimensional model requires the storage of flow variables at all radial nodes.     

Bubble Movement in Downward-Inclined Pipes

An experimental investigation with large-scale model tests of bubbles moving upward and downward in downward-inclined pipes is presented. The shape, velocity, and drag coefficient of single nonspherical air bubbles in continuous air-water flows are discussed. The bubble height depends mainly on the approach flow water velocity and the pipe slope. For stagnant bubbles, the bubble height is determined depending on these two parameters. Equilibrium of the drag and buoyancy forces is applied on single air bubbles in downward-inclined pipes. In pipes with pipe slope ranging from 0.052–0.087, the bubble drag coefficient is independent of the bubble Reynolds number. However, the bubble drag coefficient depends on the pipe slope and the approach flow water velocity. Using the approach of the equilibrium of the main forces the volume of stagnant bubbles can be predicted.     

Wall Shear Stress in the Early Stage of Unsteady Turbulent Pipe Flow

Conventionally, wall shear stress in an unsteady turbulent pipe flow is decomposed into a quasi-steady component and an “unsteady wall shear stress” component. Whereas the former is evaluated by using “standard” steady flow correlations, extensive research has been carried out to develop methods to predict the latter leading to various unsteady friction models. A different approach of decomposition is used in the present paper whereby the wall shear in an unsteady flow is split into the initial steady value and perturbations from it. It is shown that in the early stages of an unsteady turbulent pipe flow, these perturbations are well described by a laminar-flow formulation. This allows simple expressions to be derived for unsteady friction predictions, which are in good agreement with experimental and computational results.     

Numerical Calculation of Turbulence Structure in Depth-Varying Unsteady Open-Channel Flows

Unsteady depth-varying open-channel flows are really observed in flood rivers. Owing to highly accurate laser Doppler anemometers (LDA), some valuable experimental databases of depth-varying unsteady open-channel flows are now available. However, these LDA measurements are more difficult to conduct in open-channel flows at higher unsteadiness, in comparison with unsteady wall-bounded flows such as oscillatory boundary layers and duct flows. Therefore, in this study, a low-Reynolds-number k–ε model involved with a function of unsteadiness effect was developed and some numerical calculations were conducted using the volume of fluid method as a free-surface condition. The present calculated values were in good agreement with the existing LDA data in the whole flow depth from the wall to the time-dependent free surface. These values were also compared with those of unsteady wall-bounded flows. The present calculations were able to predict the distributions of turbulence generation and its dissipation, and consequently the unsteadiness effect on turbulence structure was discussed on the basis of the outer-variable unsteadiness parameter α, which is correlated with the inner-variable unsteadiness parameter ω+ in unsteady wall-bounded flows.     

Explicit Design of Commercial Pipes with No Secondary Losses

This paper presents a direct solution technique derived by using a power law for the friction factor for three cases of pipe design. The parameters used in the form of the power law presented herein are allowed to vary with pipe size and Reynolds number as well as pipe type, thus providing wide coverage with good accuracy for the computation of friction factors of commercial pipes.     

Development of Experimental Test Setup to Evaluate Abrasion Characteristics of Drainage Pipes

The longevity of drainage pipes is dependent on the material properties of the pipe, the surrounding load, the abrasive characteristics of the effluent, and the frequency, velocity, and quantity of the effluent that flows through the pipe. Each of these parameters has impacts on the age-related abrasive decrease in the wall thickness of the drainage pipe. Traditionally, the rotating pipe test is used to assess the vulnerability of drainage pipes to abrasion. To simulate more realistically the natural environment of drainage pipes, we designed an abrasion test that keeps the drainage pipe at variable slope with respect to the horizontal and continuously pumps abrasive slurry through the pipe for 30?days. The time dependent change in wall thickness of a cellulose fiber reinforced cement pipe was simultaneously compared to the change in wall thickness of a steel-reinforced concrete pipe. The wall thickness of the fiber reinforced cement pipe decreased approximately 5 times faster than the thickness of the steel-reinforced cement pipe.     

Unsteady Friction in Rough Pipes

An energy dissipation model is presented for the computation of unsteady friction losses adapted to smooth-to-rough transition and fully rough pipes. The eddy viscosity model used to compute the Reynolds stresses in turbulent flow is modified to include the effect of roughness, which is considered in the computation of the velocity profiles. The model is tested by comparing the computed transient pressures with measured data from a laboratory test facility and a prototype test. Details of the experimental setup, pressure-head measurements, and valve characteristics during transient flow conditions are presented. The quasi-steady approximation gives an inaccurate prediction of the pressure head history; however, significantly better results are obtained if the unsteady friction effects are included.     

Transient Friction in Pressurized Pipes. II: Two-Coefficient Instantaneous Acceleration–Based Model

The goal in the field of modeling of hydraulic transients is a comprehensive model for pipe networks that is computationally fast and accurate. The fastest models are the one-dimensional (1D) models that use instantaneous acceleration–based (IAB) properties, but unfortunately these models are not as accurate as the more demanding 1D convolution-based (CB) models or quasi two-dimensional models. Focusing on a single pipe, this paper investigates the fundamental behavior of the much more accurate 1D CB model to find two coefficients for use with the two-coefficient formulation of the much-used modified IAB (MIAB) model for complete closing of a downstream valve. Two coefficients are found based on the weighting function used in the CB model, and these coefficients vary along the pipe length. Simulations are compared with two experimental results from tests performed at University of Adelaide in Australia in 1995. The experimental results are for different initial Reynolds numbers of approximately 2,000 and 5,800. The results show very good agreement between simulations and experiments. The improvement of the MIAB model is not general, and for the time being, only complete closure of a downstream valve in a single pipeline at low Reynolds numbers has been investigated.     

Nonconservative Formulation of Unsteady Pipe Flow Model

A new approach to numerical modeling of water hammer is proposed. An unsteady pipe flow model incorporating Brunone’s unsteady friction model is used, but in contrast to the standard treatment of the unsteady friction term as a source term, the writers propose a nonconservative formulation of source term. Second-order flux limited and high order weighted essentially nonoscillating numerical schemes were applied to the proposed formulation, and results are in better agreement with measurements when compared with results obtained with standard form.     

Experimental Study and 3D Numerical Simulations for a Free-Overflow Spillway

The main objectives of the present work were to investigate the flow field over a spillway and to simulate the flow by means of a three-dimensional (3D) numerical model. Depending on the wall curvature, the boundary layer parameters decreased or increased with increasing distance along the spillway. The growth of the boundary layer along the spillway is better described as a function of Reynolds number than the normalized streamwise length. A simplified form of the 3D momentum equation can be used to obtain a rough estimate of the skin friction. The velocity profile in the boundary layer along the spillway is described by a velocity–defect relationship. Numerical models provide a cost-effective means of simulating spillway flows. In this study, the water surface profiles and the discharge coefficients for a laboratory spillway were predicted within an accuracy range of 1.5–2.9%. The simulations were sensitive to the choice of the wall function, grid spacing, and Reynolds number. A nonequilibrium wall function with a grid spacing equal to a distance of 30 wall units gave good results.     

Testing Block Probes for Wall Shear Stress Measurement in Water Flows

The authors developed and tested two types of sublayer pressure probe or blocks, to indirectly measure the average modulus and direction of the wall parallel component of the shear stress tensor, usually referred to as wall shear stress, or boundary shear, in water flows. The local shear stress has been investigated by means of a Preston tube for different Reynolds numbers, in a smooth rectangular section duct. The observed data have been used to calibrate the blocks. The blocks’ design is a direct outcome of the Dexter yaw meter, developed to function in an air flow. The observed pressure field associated with a given shear stress is compared to its theoretical counterpart, based on the assumption that the velocity profile satisfies a two-dimensional law of the wall, for smooth boundary.     

Role of Artificial Dissipation Scaling and Multigrid Acceleration in Numerical Solutions of the Depth-Averaged Free-Surface Flow Equations

A numerical model is developed for solving the depth-averaged, open-channel flow equations in generalized curvilinear coordinates. The equations are discretized in space in strong conservation form using a space-centered, second-order accurate finite-volume method. A nonlinear blend of first- and third-order accurate artificial dissipation terms is introduced into the discrete equations to accurately model all flow regimes. Scalar- and matrix-valued scaling of the artificial dissipation terms are considered and their effect on the accuracy of the solutions is evaluated. The discrete equations are integrated in time using a four-stage explicit Runge–Kutta method. For the steady-state computations, local time stepping, implicit residual smoothing, and multigrid acceleration are used to enhance the efficiency of the scheme. The numerical model is validated by applying it to calculate steady and unsteady open-channel flows. Extensive grid sensitivity studies are carried out and the potential of multigrid acceleration for steady depth-averaged computations is demonstrated.