Professional Documents
Culture Documents
Water Hammer With Column Separation A Historical Review
Water Hammer With Column Separation A Historical Review
Received 9 June 2005; accepted 4 August 2005 Available online 13 October 2005
Abstract Column separation refers to the breaking of liquid columns in fully lled pipelines. This may occur in a water-hammer event when the pressure in a pipeline drops to the vapor pressure at specic locations such as closed ends, high points or knees (changes in pipe slope). The liquid columns are separated by a vapor cavity that grows and diminishes according to the dynamics of the system. The collision of two liquid columns, or of one liquid column with a closed end, may cause a large and nearly instantaneous rise in pressure. This pressure rise travels through the entire pipeline and forms a severe load for hydraulic machinery, individual pipes and supporting structures. The situation is even worse: in one water-hammer event many repetitions of cavity formation and collapse may occur. This paper reviews water hammer with column separation from the discovery of the phenomenon in the late 19th century, the recognition of its danger in the 1930s, the development of numerical methods in the 1960s and 1970s, to the standard models used in commercial software packages in the late 20th century. A comprehensive survey of laboratory tests and eld measurements is given. The review focuses on transient vaporous cavitation. Gaseous cavitation and steam condensation are beyond the scope of the paper. r 2005 Elsevier Ltd. All rights reserved.
Keywords: Water hammer; Fluid transients; Column separation; Cavitation; History; Review
1. Introduction 1.1. Water hammer and column separation Modern hydraulic systems operate over a broad range of operating regimes. Any change of ow velocity in the system induces a change in pressure. The sudden shutdown of a pump or closure of a valve causes uid transients which may involve large pressure variations, local cavity formation, distributed cavitation (bubble ow), hydraulic and structural vibrations and excessive mass oscillations. In particular, the occurrence of column separation may have a signicant impact on subsequent transients in the system.
Corresponding author. Tel.: +31 40 247 2755; fax: +31 40 244 2489.
E-mail addresses: [email protected] (A. Bergant), [email protected] (A.R. Simpson), [email protected] (A.S. Tijsseling). 0889-9746/$ - see front matter r 2005 Elsevier Ltd. All rights reserved. doi:10.1016/j.juidstructs.2005.08.008
ARTICLE IN PRESS
136 A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171
Nomenclature Abbreviations in References ASCE ASME AWWA BHRA BHR FED IAHR IChemE IMechE JSCE JSME MISI PVP SMiRT American Society of Civil Engineers American Society of Mechanical Engineers American Water Works Association British Hydromechanics Research Association Group British Hydromechanics Research Group Fluids Engineering Division International Association of Hydraulic Research Institution of Chemical Engineers Institution of Mechanical Engineers Japan Society of Civil Engineers Japan Society of Mechanical Engineers Moscow Institute of Civil Engineering Pressure Vessels and Piping Structural Mechanics in Reactor Technology
Abbreviations DGCM DVCM FEM FSI GIVCM HGL MOC Scalars A a B C D E e f g H K L l p p* Q R S T Tc T cs t V 8 x cross-sectional pipe area (m2) pressure wave speed (m/s) constant in water-hammer compatibility equations (s/m2) constant in water-hammer compatibility equations (m) inner diameter of pipe (m) Youngs modulus of pipe wall material (Pa) pipe wall thickness (m) Darcy-Weisbach friction factor gravitational acceleration (m/s2) (piezometric) head (m) bulk modulus of liquid (Pa) pipe length (m) cavity length (m) gage uid pressure (Pa) absolute uid pressure (Pa) discharge (m3/s) pipeline resistance coefcient (s2/m5) cavitation severity index absolute temperature (K) valve closure time (s) duration of the rst column separation (s) time (s) ow velocity (m/s) volume (m3) distance along the pipeline (m) discrete gas cavity model discrete vapor cavity model nite element method uidstructure interaction generalized interface vaporous cavitation model hydraulic grade line method of characteristics
ARTICLE IN PRESS
A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171 137
a g D y r c
void fraction specic weight (kg/m2/s2) numerical step size, change in magnitude pipe slope (rad) mass density (kg/m3) numerical weighting factor
Subscripts A B b cs f in j m max mc P Pu pc q RV s sv u v 0 point in xt plane point in xt plane barometric column separation nal inception, initial computational section index mixture maximum negative (minus) characteristic point P in xt plane, downstream side of point P upstream side of point P positive (plus) characteristic quasi-steady reservoir-valve (see Eq. (5)) shock shock-vapor upstream, unsteady vapor initial (steady) state
Column separation is like the breaking of a solid rod or ropeGalilei already described this analogy (Rouse and Ince, 1957, 1963, p. 57)and the phenomenon can be nicely demonstrated in a simple toy apparatus (Wylie, 1999) that emits light ashes when lled with glycerin (Schmid, 1959; Chakravarty et al., 2004). Large pressures with steep wave fronts may occur when column separations collapse and the practical implications are therefore signicant. As an outcome, uid transients may lead to severe accidents (Jaeger, 1948; Bonin, 1960; Parmakian, 1985; Kottmann, 1989; De Almeida, 1991; Galante and Pointer, 2002). The general policy in hydraulic design involving uid transients is to avoid column separation. 1.2. Previous reviews Jaeger et al. (1965), Martin (1973) and Thorley (1976) provided valuable summaries with extensive bibliographies of the historical development of many aspects of water hammer including column separation. Streeter and Wylie, in their three textbooks (Streeter and Wylie, 1967; Wylie and Streeter, 1978a, 1993), give account of previous work on column separation. De Almeida (1987) reviewed the period 19781987. Beuthe (1997) gave an extensive general review with emphasis on steam condensation, omitting perhaps the rst contribution to this subject by Stromeyer (1901). During the period 19711991 an international Working Group of the IAHR carried out a major research effort with respect to column separation in industrial systems. One of the main aims of the Group was the development of computer codes and the validation of these against experimental results. The outcomes of this work have been summarized in an extensive synthesis report (Fanelli, 2000). Current collaborative research in Europe is carried on through Surge-Net. The present paper is a deliverable of the Surge-Net Work Packages 2 (Multiphase Flow) and 4 (Collection of Data). 1.3. Scope and outline of the paper Cavitation is a broad eld of research. This review is conned to the macroscopic aspects of transient vaporous cavitation with focus on the important case of column separation. The aim is to give a historical account, to summarize
ARTICLE IN PRESS
138 A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171
the state of the art, and to have a list of references as complete as practically possible. Microscopic aspects of cavitation, cavitation erosion, cavitation noise and vibration, and advanced cavitation theorysubjects not considered hereinare covered by Knapp et al. (1970), Hammitt (1980), Trevena (1987), Young (1989), Brennen (1995) and Li (2000). Valuable information may also be obtained from the review by Arndt (1981). The presence of air in water may have a signicant impact on uid transients. Air presence is considered briey herein, but the focus is on systems with very little air, as these situations will, in most circumstances, result in the highest pressures. Following this introduction, water hammer and the different types of cavitation (including column separation) are explained in Section 2. Section 3 considers mathematical models and their numerical implementation. Experimental studies of column separation and eld measurements are surveyed in Section 4 and Section 5 gives the conclusion to this paper.
2. Water hammer, column separation and cavitation 2.1. Water hammer in a historical context The conception of the theory of water hammer can, amongst others, be traced to Menabrea (1858), translated by Anderson (1976), Menabrea (1862), Michaud (1878), Von Kries (1883), Frizell (1898), Joukowsky (1900), Gibson (1908) and Allievi (1902, 1913). Joukowsky conducted experiments in Moscow in 1897/1898 when he derived his famous law for instantaneous water hammer. This law states that the (piezometric) head rise DH resulting from a fast (Tco2L/a) closure of a downstream valve is given by DH aV 0 , g (1)
in which a is the pressure of wave speed, V0 the initial ow velocity, g the gravitational acceleration, L the pipe length and Tc the valve closure time. The wave speed is estimated from Kortewegs (1878) formula s K=r a , (1a) 1 K=ED=e where K is the bulk modulus, r the mass density, E the Youngs modulus of the pipe wall material, D the inner diameter of the pipe, and e the wall thickness. The period in a pipe, 2L/a, is the time for a pressure wave to travel from the valve to a reection point (e.g. a reservoir) and back. Formula (1) is generally referred to as the Joukowsky equation, but Tijsseling and Anderson (2004) pointed out that Von Kries (1883) was actually the rst one to derive and validate it. 2.2. Column separation Column separation can have a devastating effect on a pipeline system. A well-known example is the accident at Oigawa hydropower station in 1950 in Japan (Bonin, 1960). Three workers died. The plant was designed in the early 20th century. A fast valve-closure due to the draining of an oil control system during maintenance caused an extreme high-pressure wave that split the penstock open. The resultant release of water generated a low-pressure wave resulting in substantial cavitation that caused crushing (pipe collapse) of a signicant portion of the pipeline due to the external atmospheric pressure load. Jaeger (1948) reviewed a number of the most serious accidents due to water hammer in pressure conduits. Many of the failures described were related to vibration, resonance and auto-oscillation. Two of the failures involved column separation. In one case a governor caused a valve to open suddenly and thus produced a lowpressure wave that resulted in a column separation at a change in the pipe prole. When the liquid columns rejoined, strong pressure rises caused cracking of a concrete section of the penstock. Kottmann (1989) described two accidents related to column separation in which two workers died. List et al. (1999) reported severe damage to the lining of a pump discharge pipeline, nally resulting in leaks. The cause was vapor cavity collapse. 2.3. First observations of sub-atmospheric pressures during water-hammer events Carpenter (1894) and Carpenter and Barraclough (1894) were the rst to record sub-atmospheric pressures in a water-hammer event, as noted by Joukowsky (1900, pp. 35). Joukowsky (1900, pp. 3132) was the rst to observe (see Fig. 1) and understand column separation. He explained the events in his experimental main-pipeline-gate system
ARTICLE IN PRESS
A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171 139
Fig. 1. Pressure record exhibiting column separation. Horizontal axis: time (each dot indicates half a second). Vertical axis: pressure (Joukowsky pressure of 15.3 bar). The upper horizontal line is the static pressure and the lower horizontal line is the atmospheric pressure. [Adapted from Joukowsky (1900, Fig. 17); also shown by Simin (1904, Fig. 13), and by Moshnin and Timofeeva (1965, Fig. 1)].
literally as follows: Starting at the moment of closure of the gate the water in the pipe is continuously being stopped, whereby it is being compressed, the pipe expands and the pressure increases with Dp. When this state travels with the celerity a up to the main, the latter transmits back along the pipe the pressure of the main (a little raised due to water hammer in the main itself) and a velocity of the water, which is directed in the direction of the main. This phase rst passes the cabins II and III (measuring points), as a result of which the pressure in the indicators (gages) in these cabins falls to the pressure of the main. When however the mentioned phase reaches the gate, this instantaneously causes, since the velocity of the water is directed away from the gate, a decrease of the pressure at the gate. If thereby the velocity V is so large, that according to the theory the reduced pressure would be negative, a break of the water mass will occur. The water column will be separated from the gate, ahead of which a small rareed void develops. Similar separations can also form in other parts of a liquid column, the parts towards which the reduced pressure propagated. and The condition, that the water column is separated from the gate, prolongs the duration of the reduced pressure and makes the second impact stronger than the rst, because it takes place with the velocity, at which the liquid column speeds into the rareed void. See also Simin (1904, pp. 381382). Gibson (1908) performed water-hammer experiments with closure and opening of a downstream valve in a laboratory pipeline apparatus. He indicated that a low-pressure wave initiated gas release. Strowger and Kerr (1926) warned that load rejection could cause a full column-separation in the draft tube of a reaction water turbine. Thorley (1976) attributed the rst work on vapor cavities to Hogg and Traill (1926) and Langevin and Boullee (1928). Mostowsky (1929) presented the rst theoretical analysis of column separation in an explanation of his laboratory measurements. Billings et al. (1933) presented a paper at the rst symposium on water hammer in 1933 in Chicago (Proceedings, 1933) that dealt with parting of the water column. The authors noted that dangerous instantaneous pressure rises often originated in the upper portion of a penstock, when the liquid column parted and rejoined abruptly. 2.4. Vaporous cavitation Two types of vaporous cavitation in pipelines are distinguished. The magnitude of the void fraction of the vapor within the liquid is the basis for identifying the two types. It is dened as the ratio of the volume of the vapor, 8v , to the total volume of the liquid/vapor mixture, 8m : av 8v . 8m (2)
The symbol a was introduced by Wallis (1969) in his classic textbook on two-phase ow. The two types are: (i) local column-separation (large a, a % 1) and (ii) distributed vaporous cavitation (small a, a % 0). In contrast to column separations, distributed vaporous cavitation occurs over an extended length of the pipe. All cavities maintain the vapor pressure of the liquid. 2.4.1. Cavitation inception and tensile stress The vapor pressure of the liquid is adopted as the cavitation inception pressure in most mathematical models for transient cavitation. However, there are a number of reported experiments with cavitation inception pressures (negative absolute pressure peaks) much lower than the liquid vapor pressure (Lee et al., 1985; Takenaka, 1987; Fan and Tijsseling, 1992; Simpson and Bergant, 1996). Washio et al. (1994) even observed traveling tensile waves (negative
ARTICLE IN PRESS
140 A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171
absolute pressure waves) in a dead-end branch of a main oil pipe. The magnitude of the negative absolute pressure is the dynamic tensile strength of the liquid. It is determined by the ow conditions and the cavitational properties of liquid and pipe walls. Intense pre-pressurization of the liquid tends to increase its tensile strength. Plesset (1969), Overton et al. (1984) and Trevena (1984, 1987) give in-depth treatments of tensile stresses in liquids. Recent results have been presented by Williams et al. (1999), Williams and Williams (2000) and Brown and Williams (2000). Tensile stress is a metastable condition for the liquid, which in a transient event is governed by nonequilibrium thermodynamics. Shinada and Kojima (1987) measured small negative absolute pressures of longer duration in transient cavitation tests. They attributed this to the effect of surface tension, which can be important when the distributed cavitation zone consists of many miniscule bubbles. 2.4.2. Local large vapor cavities After the studies of column separation by Joukowsky (1900) and Mostowsky (1929), it was LeConte (1937) who presented some of the rst experimental results for a local liquid column separation. He also proposed an analytical model based on rigid-column analysis. An arbitrary coefcient had to be introduced to provide a match between the analytical and experimental results. Knapp (1939) attributed many failures in the upper portion of a penstock to the pressure rise that resulted from column separation at a knee in the prole. Previously cited causes for these failures had been described as of obscure origin. Lupton (1953) presented a summary of the graphical method, with one section devoted to separation of water columns. He introduced the possibility of the formation of an intermediate vapor cavity not located adjacent to a hydraulic device (valve, turbine, etc.) or at a high point and he presented an example that exhibited the sequence of events leading to the formation of such an intermediate cavity. Moores (1953) initiated consideration of the topic in a discussion of Luptons paper, when he posed a question as to what happened when pressure waves meet. Luptons reply stated that: if the meeting of the waves were negative, and their sum exceeded the absolute head H0 existing initially, a gap would be formed. Jordan (1961) developed an analytical method for the computation of the exact locations along a pipeline of intermediate cavities. ONeill (1959) noted that most previous graphical studies overlooked the formation of intermediate cavities, which impose an internal boundary condition in the pipe. Experimental results with visualization studies of the growth and decay of intermediate cavities were presented. Sharp (1960, 1965a, b), continuing ONeills work, considered the growth and collapse of small vapor cavities produced by a rarefaction wave. An ideal spherical cavity was analyzed theoretically (Sharp, 1965a). Experimental results, which included high-speed photography of single intermediate cavities, were presented. Sharp proposed that another type of cavity also existed during the rst and succeeding rupture phases for an entirely different reason. A succession of a small number of cavities at regular intervals was asserted to form moving away from the valve (Sharp, 1965a, 1977). Reversal of the transient ow in the pipeline caused these cavities to collapse and a series of regular pulses resulted when the liquid columns rejoined. This phenomenon appears to be similar to the occurrence of distributed vaporous cavitation as discussed below. 2.4.3. Distributed small vapor cavities The difference between a local column separation and distributed cavitation was rst described by Knapp (1937b) in a paper at the second water hammer symposium in New York (Proceedings, 1937), which expanded the work in an earlier reference (Knapp, 1937a). He presented an example in which a pressure drop was produced by a rupture of the pipe just below a shutoff valve. The negative wave traveled up the pipeline unchanged until the wave front intersected the zero pressure line. The liquid between the point of intersection and the downstream reservoir cavitated partially, but column separation did not occur. Reection of the pressure wave from the reservoir brought the liquid back to its original state (without cavitation). He stated that further investigations were necessary to clear up completely such water-hammer conditions with cavitation. Knapp (1939) further developed the concept of distributed cavitation in a discussion of LeContes (1937) paper and concluded that this could not be solved by the graphical method. De Haller and Bedue (1951) presented an analytical treatment of column separation. They suggested that cavities could occur along long sections of pipeline rather than forming a local gap occupying the entire pipe cross-section. Lupton (1953) provided a description of the sequence of events associated with distributed vaporous cavitation in the transmission of a negative wave along an upward sloping frictionless pipe. Bunt (1953) presented ndings of a laboratory investigation on column separation at a valve and at high points, where the possibility of the formation of a distributed vaporous cavitation region was recognized. Jordan (1965) investigated column separation and distributed cavitation in pumping systems with horizontal, upward and downward sloping pipe sections and he studied analytically the effect of hydraulic grade line (HGL) and pipe slope
ARTICLE IN PRESS
A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171 141
(a)
(d)
(b)
(e)
(c)
(f )
Fig. 2. Dotted lines: vapor head; thin dashed lines: steady state head; dashed lines: unsteady head. (Left) Fast-closing valve generating cavitation in (a) downward sloping pipe, (b) horizontal pipe, and (c) upward sloping pipe. [Adapted from Zielke and Perko (1985, Fig. 3)]. (Right) Cavity formation at a knee with pipe slope (d) increasing, (e) decreasing but still upward, and (f) decreasing to downward [adapted from Zielke and Perko (1985, Fig. 4)].
on the formation of cavitation zones. A cavitation zone may result from the passage of a negative wave traveling in the direction of decreasing steady state pressure. In contrast, if the steady state pressure increases in the direction of wave propagation a vaporous cavitation zone cannot occur (Zielke and Perko, 1985; Simpson, 1986; Simpson and Wylie, 1989). This is portrayed in Fig. 2. 2.5. Gaseous cavitation Many papers have addressed the effects of free gas, dissolved gas and gas release on transients in liquid-lled pipelines (Swafeld, 19691970, 1972a, b; Enever, 1972; Kranenburg, 1974b; Martin and Padmanabhan, 1975; Sundquist, 1977; Wiggert and Sundquist, 1979; Wylie, 1980; Baasiri and Tullis, 1983; Chaudhry et al., 1990; Hadj-Taieb and Lili, 1998, 1999, 2000, 2001; Kessal and Amaouche, 2001). One of the main features of liquids is their capability of absorbing a certain amount of gas with which they come into contact through a free surface. In contrast to rupture and vaporization, which takes only a few microseconds, gas release is a diffusive process, which takes several seconds (in water at room temperature). Gas absorption is much slower than gas release (order of minutes) (Zielke et al., 1989). Dissolved gas is an important consideration in sewage water lines, aviation fuel and oil lines. Gas comes out of solution when the pressure drops. It may be assumed that released gas stays in cavities and does not immediately redissolve following a rise in pressure. Kobori et al. (1955), Pearsall (1965) and Raiteri and Siccardi (1975) showed that the presence of free gas reduces the wave speed and, according to Eq. (1), the transient pressure variations.
ARTICLE IN PRESS
142 A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171
2.6. Short-duration pressure peaks following cavity collapse Short-duration pressure peaks exceeding the Joukowsky value given by Eq. (1) may occur after cavity collapse. The phenomenon is explained for a frictionless reservoir-pipe-valve system, Fig. 3(a). Instantaneous valve closure causes a head rise DH a=gV 0 . The head rise travels towards the reservoir, where it reects negatively, and returns at the closed valve at time 2L/a. Fig. 3(b) shows the pattern of propagation and reection of pressure waves in the distancetime plane. Fig. 3(c) shows the heads at the valve calculated with a water-hammer model with (thick solid line) and without (thin solid line) a cavity forming at the valve. The liquid ows in the reverse direction (towards the reservoir) at time 2L/a. Complete stoppage of the ow at the valve now requires a head drop of (a/g)V0 relative to H0. This head drop (thin solid line in Fig. 3(c)) would result in a head less than the vapor head (horizontal dashed line in Fig. 3(c)). When the head is not allowed to drop below the vapor head (thick solid line in Fig. 3(c)), the reverse velocity at the valve is not reduced to zero at time 2L/a, but to (Mostowsky, 1929, Fig. 7(b) herein with DH in H 0 H b ) V 0 DV vc V 0 where DH in H 0 H b p =g; v DH in oDH (3a) gDH in , a (3)
is the head drop (relative to H0) which starts vaporization. Here, H0 is the steady state head at the valve, Hb the barometric head, p the absolute vapor pressure at temperature T, and g rg the specic weight of the liquid. Now, at v the closed valve, the ow is towards the reservoir, the liquid detaches from the solid and a cavity begins to grow. The cavity acts as a xed-head boundary condition for the pressure waves. At each successive time interval (4L/a, 6L/a, etc.) a pressure wave reects off the vapor cavity and the reverse liquid velocity decreases by an amount of 2DV vc . At a certain time this velocity changes direction, the cavity then begins to shrink and it nally collapses at point A in Fig. 3(b) and (c). The head directly caused by the cavity collapse is less than the head H0+DH generated by valve closure, but the maximum head (larger than Joukowsky pressure head rise!) occurs, in this example, at a time of about 6L/a (point B in Fig. 3(b) and (c)) in the form of a short-duration peak. If the cavity collapse would have taken place exactly at the arrival of a pressure wave front, at times that are multiples of 2L/a, then a short-duration peak would not have occurred.
(b)
(c)
(a)
Fig. 3. Formation of a short-duration pressure peak. (a) Reservoir-pipe-valve system. (b) Wave paths in distance-time plane. (c) Piezometric-head history at valve.
ARTICLE IN PRESS
A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171 143
In a similar consideration, Wylie and Streeter (1993, pp. 192196) showed that the time of existence of the rst vapor cavity approximately is T cs 2a L DH 2L V0 , gDH in a DH in a (4)
thus conrming the formula found by Mostowsky (1929) and shown in Fig. 7(a) herein. Angus (1935, 1937) presented one example in which the collision of a liquid column with a closed valve resulted in a short-duration pressure peak of large magnitude. The duration of the pressure peak was about one-tenth of a 2L/a period. Pressure traces measured by Binnie and Thackrah (1951) indicated the existence of a short-duration pressure peak, which exceeded the main shock pressure. The authors attributed the higher pressures to the existence of additional pressure waves caused by reections from bends, sockets and other discontinuities in the pipeline. Their analysis was based on rigid-column theory and on water-hammer theory. Some of ONeills (1959) analytical examples exhibited short-duration pressure peaks, but his experimental pressure records did not. Gayed and Kamel (1959) clearly explained the occurrence of short-duration peaks and they claimed their appearance in one of the measurements (Fig. 11(a) in their paper). Heaths (1962) graphical-method water-hammer analysis predicted short-duration pressure peaks, but these were not apparent in his experimental results. Sharp (1960) published in Nature an example with shortduration pressure peaks and intermediate cavities. Walsh (1964) and Li and Walsh (1964) estimated the maximum possible pressure rise following the collapse of the rst cavity at an upstream valve as a DH max V f 2H RV , g (5)
where Vf is the velocity of the liquid column at the valve just before cavity collapse and HRV the difference of reservoir head and vapor head at the valve. Walsh (1964) presented results for a number of experimental runs, a couple of which possibly show short-duration pressure peaks. The time of existence of the rst cavity was quite long with respect to 2L/ a. The velocity Vf was estimated from experimental results using Eq. (5). Equations similar to Eq. (5) were presented by Moshnin (1961), Moshnin and Timofeeva (1965), De Almeida (1983), Simpson and Wylie (1985) and Wylie and Streeter (1993, p. 194 with DH max DH 2DH in ) for the case of instantaneous closure of upstream and downstream valves, where DHmax can be more than two times the Joukowsky value (1). The situation can even be worse when intermediate cavities form (Sharp, 1960). Kottmann (1989) used rigid-column theory to show that the collapse of a midpoint cavity caused pressure rises of three times Joukowsky. Yamaguchi and Ichikawa (1976) and Yamaguchi et al. (1977) presented oscilloscope traces of experimental results exhibiting what appear to be the rst results in the literature that clearly show short-duration pressure peaks. Tarasevich (1980) considered column separation at an upstream valve. He predicted maximum pressures as a function of initial ow velocity. His analysisbased on rigid column theorywas conrmed by the experiments of Smirnov and Zubov (1972). Martins (1983) measurements concerned limited cavitation where the time of existence of the cavity at the valve is relatively short. The experimental results showed that the maximum pressure may exceed the Joukowsky pressure rise in the form of a short-duration pressure peak as shown in Fig. 4(a). Unfortunately, the reservoir pressure was rising during the experiment because the tank was too small (Martin, 1989). The measured reservoir pressure shown in Fig. 4(b) was thought to be typical for all experimental runs. Kojima et al. (1984) presented a mathematical model for predicting both the pressure rise associated with cavity collapse and the duration of the column separation. Their experimental results exhibited short-duration pressure peaks. Simpson (1986) showed a range of short-duration pressure peaks measured in a reservoirupward sloping pipeline valve system. Due to the upward slope of the pipe, the vapor cavity was conned to be adjacent to the valve with no distributed cavitation along the pipe (at least until the collapse of the rst vapor cavity).
2.7. Severity of cavitation Martin (1983) inferred the severity of cavitation from the duration, Tcs, of the rst, mostly largest, column separation at the valve. He introduced the cavitation severity index S T cs a=2L, which according to Eq. (4) is equivalent with S DH=DH in . Paredes et al. (1987), Carmona et al. (1987, 1988) and Anderson et al. (1991) proposed a somewhat different index. Bergant and Simpson (1999a) made a classication based on maximum heads (shown in Fig. 5) whether or not exceeding the Joukowsky value (Eq. (1)). The steepness of wave fronts is another important parameter in the assessment of dynamic loads on pipe systems and their supporting structures.
ARTICLE IN PRESS
144 A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171
(a)
(b)
Fig. 4. Experimental result showing short-duration pressure peaks [adapted from Martin (1983, 1989)]. Absolute HGL variation with time, (a) at valve, and (b) in reservoir.
(a)
(b)
Fig. 5. Computed maximum head at valve as function of initial velocity for varying reservoir-head [adapted from Bergant and Simpson (1999a)]: (a) downward sloping pipe, (b) upward sloping pipe.
3. Mathematical models and numerical methods A comprehensive investigation of column separation and cavitation in pipe systems was not possible until the 1960s [see Hager (2001)] due to the unavailability of computers. Most studies used graphical and arithmetic procedures originally set forth by Gibson (19191920), Schnyder (1932) [see Hager (2001)], Angus (1935), Bergeron (1935, 1950), Stepanoff (1949) and Parmakian (1955). The rst computer-oriented procedures for the treatment and analysis of water hammer include work by Thibessard (1961), Lai (1961), Streeter and Lai (1962), Streeter (1963, 1964, 1965), Van De Riet (1964), Vreugdenhil (1964) and Contractor (1965). The rst computer models for column separation were ` developed by Thibessard (1961) at Liege in Belgium, Streeter and Wylie (1967), Baltzer (1967a, b) and Weyler (1969) at the University of Michigan, and Vreugdenhil (1964) and Siemons (1967) at Delft in The Netherlands.
ARTICLE IN PRESS
A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171 145
3.1. Water-hammer equations The water-hammer equations are applied when the pressure is above vapor pressure. They comprise the continuity equation and the equation of motion: @H @H a2 @V V V sin y 0, @t @x g @x @H @V @V fV jV j V 0, @x @t @x 2D (6)
(7)
where H is the piezometric head (HGL), t the time, V the ow velocity, x the distance along the pipeline, y the pipe slope, a the pressure wave speed, g the gravitational acceleration, f the Darcy-Weisbach friction factor [its history is nicely described by Brown (2002)] and D the inner diameter of the pipe. For most engineering applications, the convective terms V(qH/qx), V(qV/qx) and V sin y are very small compared to the other terms and therefore neglected. Streeter and Wylies (1967) classic textbook led the world to the method of characteristics (MOC) as the numerical method to solve Eqs. (6) and (7) on a digital computer. The MOC uses the compatibility relations that are valid along the positive (C+: Dx=Dt a) and negative (C: Dx=Dt a) characteristic lines. These are for the xt grid shown in Fig. 6, H P C pc Bpc QPu Rpc QPu QA , (8) H P C mc Bmc QP Rmc QP QB , (9)
where Cpc, Bpc and Rpc are constants for the positive or plus characteristic (pc), and Cmc, Bmc and Rmc are constants for the negative or minus characteristic (mc). In anticipation of discrete cavity models, the discharge may be discontinuous with up- and downstream values QPu and QP, respectively. In nite-difference form Eqs. (8) and (9) for the computational section P with index j become H t H tDt j j1 a f Dx tDt fQu t Qj1 g Qu t QtDt 0, j j j1 2 gA 2gDA a f Dx tDt fQt Qu j1 g Qt Qu tDt 0, j1 j 2 j gA 2gDA (10)
H t H tDt j j1
(11)
where Dx is the reach length and Dt is the time step as indicated in Fig. 6.
Fig. 6. The method of characteristics (Dx aDt) staggered grid for a reservoir-pipe-valve system.
ARTICLE IN PRESS
146 A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171
(a)
(b)
(c)
Fig. 7. Theoretical and experimental analysis of column separation in a reservoir-pipe-valve system. (a) Theoretical pressure as function of time [adapted from Mostowsky (1929, Fig. 7)]. (b) Theoretical ow velocity as function of time [adapted from Mostowsky (1929, Fig. 8)]. (c) Experimental pressure as function of dimensionless time [adapted from Mostowsky (1929, Fig. 16)].
3.2. Discrete single-cavity models Single-cavity numerical models deal with local column separations as described in Section 2.4.2. A single cavity is allowed to form either at a closed end, at a high point or at a change in pipe slope. All graphical solutions of water hammer employ this modeling approach. Rigid column theory has also been used to compute the behavior of systems with single cavities. Mostowsky (1929) analyzed column separation at a downstream valve. Fig. 7(a) shows the pressuretime diagram for a frictionless pipe where the duration Tcs of the column separation is an integer multiple of the wave travel time 2L/a. This gure is a conrmation of Eq. (4). Fig. 7(b) is the corresponding velocitytime diagram that determines the time of separation. Mostowsky (1929) performed measurements in a 29.5 m long, 2 inch diameter pipe. The measured pressure history in Fig. 7(c) shows that, unlike Fig. 7(a), the second pressure rise is lower than the rst one, and the measured number of 11-times 2L/a is not the 12-times predicted by Eq. (4). Mostowsky (1929), attributing these discrepancies to friction, performed a rigid-column analysis with a quadratic friction law. He found closed-form solutions, for example a formula relating the effective friction coefcient to p1 and p2 (these are indicated in Fig. 7(c)), and an improved formula for Tcs giving a value of 10.6-times 2L/a, which is in close agreement with the experiment. Angus (1935, 1937) included a single vapor cavity at a boundary in the graphical method (Fig. 8). For a pump failure on a discharge line, the possibility was investigated of a cavity forming on the pipe side of a check valve near the pump.
ARTICLE IN PRESS
A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171 147
After the cavity formed and the liquid column returned to the closed valve, the pressure rise was found to be about four times the normal value, because of a large short-duration pressure peak of the type discussed in Section 2.6. Bergerons (1939, 1950) treatment of column separation is similar to that of Angus (1935). He indicates that the underpressure in a cavity is not the barometric vacuum, but rather the vapor pressure at the temperature of the liquid. He neglects the inuence of gravity on the shape of the liquidvapor interface. In his application of the graphical method, friction was lumped at the upstream end by using an orice. A number of cycles of successive formation and collapse of a column separation were considered. Friction loss resulted in damping and the mean value of head rise decreased with every cycle, as did the size of the cavity. Gayed and Kamels (1959) analysis included variable-length rigid-column theory and water-hammer theory with noninstantaneous excitation. Kephart and Davis (1961) used rigid-column theory to determine the magnitude of the pressure rise when liquid columns rejoined in pump discharge lines equipped with check valves at the pump outlet. This method served as a preliminary design technique. Escande (1962) included cavities and lumped nonlinear friction in his graphical calculations. Ruus et al. (1984) performed an extensive numerical study of column separation located at a high point. 3.3. Discrete multiple-cavity models 3.3.1. The discrete vapor cavity model (DVCM) The discrete (concentrated, lumped) vapor cavity model (DVCM) is the most commonly used model for column separation and distributed cavitation at the present time. It is easily implemented within standard water-hammer software and it reproduces many of the physical features of column-separation events. Wylie and Streeter (1978a, 1993) have described the DVCM in detail and they provided a FORTRAN computer code for its implementation. We will use their notation. Cavities are allowed to form at any of the computational sections (grid points) if the pressure is computed to be below the vapor pressure. Vapor cavities are thus conned to computational sections as sketched in Fig. 9 and pure liquid with a constant pressure wave speed is assumed in between. The absolute pressure in a cavity is set equal to the vapor pressure: p p . v (12)
ARTICLE IN PRESS
148 A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171
Fig. 9. Denition sketch for discrete vapor cavity model (adapted from Tijsseling, 1993, Fig. 4.9).
Fig. 10. Numerical oscillations in DVCM results [adapted from De Vries (1972, Fig. 42)]. Pressure-head histories at different positions after pump shutdown.
The upstream and downstream discharges QPu and QP at a cavity are computed from the compatibility relations (8) and (9). The cavity volume 8 follows then from Z t d8 QP QPu dt. (13) QP QPu or 8 dt tin The numerical integration of Eq. (13) in the MOC with staggered grid (Fig. 6) is usually given by
t2Dt t2Dt 8t 8t2Dt fcQt Qt 1 cQP QPu g2Dt P Pu
(14)
in which 8t and 8t2Dt are the volumes at the current time and at 2Dt earlier, and c is a numerical weighting factor. If the cavity volume becomes zero or negative, the cavity disappears and the computation returns to the standard waterhammer procedure. Provoost (1976) used a cavity closure condition that exactly satises the mass balance. Although the DVCM was correctly formulated by Thibessard (1961), Vreugdenhil (1964), Streeter (1969) and Tanahashi and Kasahara (1969), these investigators did not apply the model to regions of distributed cavitation. This was done by De Vries (1972, 1973, 1974) when he simulated the experiments performed by Vreugdenhil et al. (1972) and Kloosterman et al. (1973). De Vries (1972) was the rst to report on numerical oscillations induced by the condensation of a region of distributed cavitation (Fig. 10). To suppress these oscillations he added small amounts of free gas to the discrete cavities. In fact he then used the discrete free gas cavity model developed by Brown (1968), but with nonzero vapor pressure. This model is described in Section 3.3.2. Evans and Sage (1983) had condence in the DVCM and used it for the water-hammer analysis of a 110 km long gravity line of 0.91.6 m diameter. Bergant and Simpson (1999b) incorporated cavitation inception with negative absolute pressure (tensile stress) in the DVCM. Numerical and experimental results, both with negative pressure spikes, were compared. The negative pressure spike at cavitation inception did not signicantly affect the results.
3.3.2. The discrete gas cavity model (DGCM) The discrete gas cavity model (DGCM) is similar to the DVCM, with a quantity of free air assumed to be concentrated at each computational section. The pressure in a cavity satises the ideal gas law: p p 8 p p 80 constant; v 0 v (15)
where the free gas is assumed to behave isothermally, which is valid for tiny bubbles. Large bubbles and column separations tend to behave adiabatically. The only difference with the DVCM is in the 8p* curve as shown in Fig. 11,
ARTICLE IN PRESS
A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171 149
Fig. 11. p*8 curves for DVCM (L-shaped solid line, 80 0 mm3) and DGCM (thin solid lines, 80 0.01 mm3, 80 0.1 mm3, 80 1 mm3, 80 5 mm3, 80 10 mm3) with p0* 1 bar and pv* 0.2 bar. The value of pv* (water at 60 1C) is chosen large to show its inuence.
Fig. 12. Wave-speed ratios [adapted from Liou (2000, Fig. 5)]. Dots: numerical (DGCM) wave speed/gas-free wave-speed. Lines: theoretical mixture wave-speed/gas-free wave-speed.
which also tells us that the DVCM is a limit case of the DGCM. The DGCM exhibits dispersion of rarefaction waves and steepening of compressive waves and thus allows shock fronts to form. Brown (1968) presented the rst attempt at describing column separation with the effects of entrained air. Entrained air was assumed to be evenly distributed in concentrated pockets at equal distances along the pipeline. The presence of air decreased the overall pressure wave speed. A pressurevolume ideal-gas relation was assumed to describe the behavior of the concentrated air pockets, whose locations were assumed to be permanent with no change due to the prevailing direction of ow. The presence of entrained air was neglected above a certain head, where the solution reverted back to normal water-hammer computations. After Brown (1968), this model was further developed by De Vries (1973), Provoost (1976), Provoost and Wylie (1981) and Wylie (1984). In a rigorous mathematical treatise of the DGCM, applying Von Neumann analysis to a linearized set of equations, Liou (1999, 2000) nicely showed that the numerical wave speed converges to the theoretical (physical) one (see Fig. 12). He thus explained why the DCGM exhibits nonlinear variable wave-speed features. The DGCM model was successfully used by Enever (1972), Tullis et al. (1976), Aga et al. (1980), Ewing (1980), Suda (1990) and Barbero and Ciaponi (1991). Ewing (1980) discussed various damping mechanisms in liquidgas mixtures. Capozza (1986) applied the method to condenser cooling circuits. Wylie (1992) and Wylie and Streeter (1993,
ARTICLE IN PRESS
150 A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171
pp. 188192, 202205) compared DGCM results with analytical results, with bubbly ow experimental data from Akagawa and Fujii (1987) and with experimental column-separation data from Simpson (1986). 3.3.3. Features and limitations of discrete cavity models When the absolute pressure reaches the vapor pressure, cavities or bubbles will develop in the liquid. In discrete cavity models (like DVCM and DGCM) the cavities are concentrated at the grid points. Between the grid points pure liquid is assumed for which the basic water-hammer equations remain valid. This means that the pressure wave speed a is maintained (and convective terms neglected) between grid points in distributed cavitation regions. However, in bubble ow the actual pressure wave speed is very low and pressure-dependent. These matters are implicit in the model (Liou, 1999, 2000). Complete and clear descriptions of discrete cavity models are given by Wylie (1984) and Zielke and Perko (1985). In discrete cavity models the cavities do not move. This is consistent with the acoustic approximation: since the overall time scale is acoustic (water hammer), the displacements of vapor bubbles are small. However, Vreugdenhil (1964) took into account, within the DVCM, the motion of liquidvapor boundaries. The maximum length, l 8=A, of a cavity must be small compared to the spatial mesh size. Simpson and Bergant (1994a) recommended l o0:1. (16) Dx For distributed cavitation regions, condition (16) is mostly fullled. If not, the discrete cavity model is not valid and the application of models for two-phase plug ow, slug ow or open-channel ow should be considered. For column separations, condition (16) may sometimes be violated, which is acceptable in the opinion of the authors, since column separation is a local phenomenon, and only a few grid points are concerned. However, care should be taken when l/Dx becomes larger than 1. In that case liquidvapor boundaries moving from grid point to grid point should be explicitly modeled. From a physical microscopic point of view the macroscopic DVCM is not correct; a two-phase ow (Wallis, 1969) approach would be better. However, during vaporous cavitation, Eq. (12) is physically a strong condition, if there is no free gas or gas release involved. Furthermore, the continuity equation, that is the mass balance, is satised throughout. Note that the mass of vapor has been neglected, just as the inuence of radial pipe displacements on the cavity volume. The DVCM is a relatively simple model, which is able to cover the essential phenomena in transient cavitation. It ts in with the standard MOC approach, so that it can be used in general water-hammer computer-codes. Its main deciency is in the numerical oscillations and unrealistic spikes appearing in the calculated pressure histories, when regions of distributed cavitation occur (De Vries, 1972; Bergant and Simpson, 1999a). One way of partly suppressing the oscillations and spikes is by assuming small amounts of initial free gas in the grid points (De Vries, 1973; Provoost, 1976; Wylie, 1984; Zielke and Perko, 1985; Simpson, 1986; Barbero and Ciaponi, 1991; Bergant, 1992; Simpson and Bergant, 1994a). Condition (12) is then replaced by Eq. (15), so that DVCM becomes DGCM. The recommended free gas void fractions to be used are of the order of 107 at standard atmospheric conditions. The numerical integration of Eq. (15) may also affect the amount of oscillations and spikes (Provoost and Wylie, 1981; Simpson and Wylie, 1985; Simpson and Bergant, 1994a; Liou 1999, 2000). The application of a numerical lter may be considered (Vliegenthart, 1970; Kranenburg, 1974a). Discrete cavity models are able to predict accurately the maximum pressures in a system, but they are usually poor in the prediction of the timing of repeated cavity formation and collapse. The inclusion of unsteady friction (Section 3.8) might help in this respect. 3.4. Shallow-water ow models Shallow-water (open-channel) modeling provided the rst real attempt at a more realistic description of cavitating ow. Bubbles were assumed to form, rise quickly and agglomerate into a single long thin cavity, when the pressure reached the vapor pressure. Li (1962, 1963) and Li and Walsh (1964) presented a study of the mechanics of pipe ow following a column separation at an upstream closing valve. Spreading of the vaporliquid interface was described by shallow-water theory with two quasi-linear partial differential equations. The study revealed that the spreading of the surface may be neglected in computing the pressure resulting from the cavity collapse. Rigid-column theory was used for the liquid, which required the assumption that the time of existence of the cavity is long compared to 2L/a. Baltzer (1967a, b) used a shallow-water ow model to calculate the shape, movement and collapse of a vapor cavity formed at the upstream side of a valve. A detailed description was presented of the sequence of events. This is illustrated in Fig. 13.
ARTICLE IN PRESS
A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171 151
Fig. 13. Simulated transient pressures at gate valve and concurrent, free-surface proles of vapor cavity (adapted from Baltzer, 1967b, Fig. 4).
Siemons (1966, 1967) assumed that the thickness of the cavity layer was small in comparison to the diameter of the pipe. He stressed that the rise in pressure at cavity collapse would not be great. Kalkwijk and Kranenburg (1971) noted that Siemons results did not maintain a mass balance at the boundary of the cavity and therefore they questioned the validity of his conclusion concerning the generation of high pressures. The transition from the water-hammer region to the cavitating region is one of the major problems in the shallow-water approach. For this reason, Kranenburg (1974a) concluded that the description of cavitating ow and column separation by shallow-water theory did not seem attractive. Another problem was the appearance of gravity waves. Furthermore, the model is physically incorrect for vertical pipes. Marsden and Fox (1976) used a similar approach to cavitation as that taken by Baltzer (1967a, b). High pressures were not predicted by their model and the authors concluded that uneconomic overdesign of pipelines resulted if the cavity collapse was computed by normal techniques. Fox and McGarry (1983) presented a variant of the discrete cavity model with a cavity assumed to occupy the upper portion of the pipe and to be spread over a Dx length. Thermodynamic effects were included, but the authors concluded that their inuence was insignicant if the vapor pressure of the liquid was small. 3.5. Two-phase ow models Kalkwijk and Kranenburg (1971, 1973) presented two theoretical bubble ow approaches to describe distributed vaporous cavitation in a horizontal pipeline. The rst approach was entirely based on the dynamic behavior of gas bubbles. The method failed at the point where the radii of the bubbles exceeded a critical value. At this size the bubbles became unstable and the MOC characteristic lines became imaginary. The incorporation of added liquid mass might solve this problem (Geurst, 1985; Thorley and Wiggert, 1985; Gale and Tiselj, 2003). The second approach distinguished between regions with and without cavitation. Different systems of equations held for the water-hammer and the vaporous-cavitation regions (Kranenburg, 1972). The pressure wave speed in the cavitation region, relative to the uid particles, is reduced to zero, because the pressure in such a region is constant at vapor pressure. Analytical expressions were developed for the velocity and void fraction in horizontal vaporous zones. When a cavitation region stopped growing, a shock formed at the transition from the water-hammer to the vaporous region, which penetrated into the cavitation region.
ARTICLE IN PRESS
152 A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171
Kranenburg (1972, 1974b) developed a simplied one-dimensional mathematical model, referred to as the simplied bubble ow model. The slope of the pipe and the inuence of gas release were both considered (see also Section 3.7). Kranenburg (1974a) found that there was considerable difculty in using the MOC due to the pressure dependence of the wave speed as a result of free gas. He asserted that discontinuities or shocks between the water-hammer and vaporous region should be tted in the continuous solution only for simple cases. To simplify his modeling approach, the bubble ow regime was assumed for the whole pipe, even for the water-hammer regions. A modied surface-tension term was used to achieve this simplication. This model did not show explicit transitions between the water-hammer and vaporous cavitation regions. Kranenburg (1972) used a Lax-Wendroff two-step numerical scheme despite the occurrence of shock waves. An articial viscosity was used to suppress the numerical instability that spread the developing shock front over a number of grid points. Column separation was explicitly taken into account at mesh points where it was expected to occur. Wylie and Streeter (1978b, 1993) developed a similar case-specic analytic model involving vaporous cavitation. 3.5.1. Two-phase ow equations for distributed vaporous cavitation The two-phase equations describing a vaporous cavitation region are the following equations of continuity and motion (Streeter, 1983; Bergant and Simpson, 1992; Wylie and Streeter, 1993): @av @av @V m Vm 0, @t @x @x @V m @V m fV m jV m j Vm g sin y 0, @t @x 2D (17) (18)
where av is the void fraction of vapor and Vm the mixture velocity. The pressure is assumed constant at vapor pressure, so that only gravitational and friction forces act. The two equations are valid for small void fractions (av oo1) and up to a temperature of about 330 K (Hatwin et al., 1970), so that thermodynamic effects are not important. The conventional Darcy-Weisbach friction factor f for fully liquid ow is assumed in Eq. (18), because the effect of bubbles on friction loss can be ignored for small void fractions (Grifth, 1987). The solution of Eq. (18) for Vm depends on the pipe slope being upward, downward or horizontal and it depends on the velocity of the liquidvapor mixture at the time of cavitation inception. Assuming that Vm is a uniform velocity independent of xin each individual cavitation zone, the different results of integration for Vm can be found in Streeter (1983), Bergant (1992), and Bergant and Simpson (1999a). 3.5.2. Shock equations for the condensation of distributed vaporous cavitation A vaporous cavitation region expands in size due to a propagating low-pressure wave. Eventually, the cavitation region stops expanding and condensation starts. The movement of the interface (shock-wave front) separating the onephase (liquid) and the two-phase (liquidvapor mixture) uids is described by shock equations. Isothermal conditions across the interface of innitesimal width are assumed (Campbell and Pitcher, 1958). The shock equations are (Bergant, 1992; Bergant and Simpson, 1992): hg i as 2 H s H sv av V V m 0, (19) a gH s H sv V V m V V m as 0, (20)
where as is the shock-wave speed, Hs the head on the liquid side of the wave front and Hsv the head on the cavitation side of the wave front. Relations (19) and (20) couple the water-hammer Eqs. (6) and (7) and the two-phase ow Eqs. (17) and (18). 3.6. Interface models This type of model combines discrete cavities with two-phase ows. Flow regions with different characteristics (that is: water hammer, distributed vaporous cavitation, end cavities and intermediate cavities) are modeled separately, while the region interfaces are tracked. Kranenburg (1974b) modeled column separation at a valve in combination with the description of vaporous cavitation regions using his bubble ow model. The model was applied to the inclined-pipe experiments of Baltzer (1967a, b). Streeter (1983) was the rst to develop a combined analysis for modeling column separations at high points and a number of distributed vaporous cavitation regions, while retaining the shock-tting approach to explicitly compute the
ARTICLE IN PRESS
A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171 153
locations of transitions between water-hammer and vaporous-cavitation regions. Gas release was not considered, thereby removing the problem associated with the variable wave speed due to the presence of free gas. The model was applicable to pipes at any angle with the horizontal. Many separate cavitation zones could be modeled, as well as their collapse and reforming. The equations developed for the vaporous regions were for various combinations of slope and initial velocity at cavity inception. A computational section became vaporous once a pressure less than the vapor pressure had been computed from the water-hammer MOC. A time-line interpolation was used to estimate the rst occurrence of vapor at a computational section. The development by Streeter (1983) also considered the formation of discrete vapor cavities. It was asserted that a discrete cavity may only form if the angle with the horizontal between two adjacent pipe sections decreased in the downstream direction, such as at a high point. If vapor pressure occurred at such a section, a column separation was assumed to form. For all other pipe slope conditions a distributed cavitation region was assumed. This approach did not account for the possibility of intermediate cavities due to the interaction of two low-pressure waves. Bergant (1992) and Bergant and Simpson (1992) extended a standard DVCM with vaporous cavitation zones, shock waves and various types of discrete cavities. The model, referred to as a generalized interface vaporous cavitation model (GIVCM), handles a number of pipeline congurations (sloping and horizontal) and various interactions between water-hammer regions, vaporous-cavitation regions, and intermediate and boundary column-separations. For example, Fig. 14 shows a typical sequence of transient events in a horizontal pipe including growth and collapse of a discrete cavity, propagation of a vaporous cavitation zone and two shock-wave fronts. In essence, the GIVCM algorithm maintains the same basic structure as the DVCM and it is therefore simpler than previous interface models (Simpson, 1986). A loop for the shock treatment at appropriate computational sections was added to the basic DVCM loop and a module for combined discrete-cavity and distributed-cavitation computation was incorporated. Flags to control the correct physical behavior of various phase interactions and to identify possible new interactions strengthened the algorithm. 3.7. Modeling of gas release Dijkman (1968) and Dijkman and Vreugdenhil (1969) extended Siemons (1967) shallow-water ow model by considering gas release into a single cavity at a high point. The pressure rise, after compression of the cavity incorporating released gas, was concluded to be less serious than for the vapor-only case. Gas ow equations were solved in combination with the shallow-water equations. The MOC was applied to solve the resulting fourth-order hyperbolic system. The authors concluded that it was uncertain how the collapse pressures may be computed and suggested an approach of attempting to prevent the occurrence of cavitation. Kranenburg (1974a) presented an extensive work on the effect of free gas and gas release on cavitation in pipelines. Gas release in the vaporous region was concluded to cause damping of the pressure peaks caused by the collapses of
Fig. 14. Cavitation and shock formation in a horizontal pipe [adapted from Bergant and Simpson (1992)].
ARTICLE IN PRESS
154 A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171
column separations. To support this contention, a hydrodynamic energy balance was computed for a column separation in a reservoir/horizontal-pipeline/valve system. Contributions to the energy balance included elastic energy of the liquid and pipe wall, elastic energy of the free gas, work done at the reservoir, dissipation caused by shock waves, and dissipation caused by wall friction. The elastic energy of the free gas was small compared with the dissipation terms. This explained the relatively small inuence of gas release at the column separation void. There was only slight damping due to wall friction. The conclusion was drawn that the marked energy loss may be attributed to the dissipation at the shock-wave fronts due to heat-transfer and viscosity effects. In summary, Kranenburg (1974b) concluded that the inclusion of gas release had no effect when only cavitating ow occurred, whereas the inuence was considerable when column separation occurred in combination with cavitating ow. Gas release in the cavitating ow region adjacent to a column separation diminished the duration of subsequent separations and thus the maximum pressures upon collapse. Barbero and Ciaponi (1991) examined the inuence of initial free gas and gas release in their calculations. 3.8. Modeling of friction In horizontal regions of distributed cavitation the pressure is constant and the ow accelerations are governed by friction forces only, see Eq. (18). Therefore, an accurate modeling of friction is a prerequisite. To do so, the friction factor f used in Eqs. (7) and (18) can be expressed as the sum of a quasi-steady part fq and an unsteady part fu (Zielke, 1968; Vardy, 1980; Vardy and Brown, 2000, 2003, 2004; Bergant et al., 2001). In this way a number of unsteady friction models have been incorporated in standard discrete cavity models (Safwat et al., 1986; Dudlik, 1999; Axworthy and Chabot, 2004). Shuy and Apelt (1983) performed numerical analyses with ve different friction models (steady, quasi-steady and the three unsteady friction models developed by Carstens and Roller (1959), Trikha (1975) and Hino et al. (1977)) that had been incorporated into a standard DVCM. The authors studied slow transients in two long pipelines (2.3 and 9 km). For the case of pure water-hammer they found only small differences in the results of the ve models, but for the case with column-separation large discrepancies occurred. Brunone and Greco (1990), Golia and Greco (1990) and Brunone et al. (1991) used the DVCM in combination with Golias (1990) unsteady friction model. Numerical results were compared with results of measurements of rapid waterhammer and column separation. Signicant discrepancies between experiment and theory were found for all runs when using a quasi-steady friction term. Results obtained with Golias friction model showed an improved agreement between the computed and measured results. However, the agreement was better for the water-hammer case than for the column-separation case. Kojima et al. (1984) developed a gas non-bubbly ow model using unsteady pipe friction that accurately predicted their experimental results. The inuence of unsteady laminar friction was studied. Bergant and Simpson (1994a) investigated the performance of the quasi-steady and three distinct types of unsteady friction models, namely those of Zielke (1968), Hino et al. (1977) and Brunone et al. (1991). The Zielke and the Brunone et al. unsteady friction models gave the best t with experimental data for the case of pure water-hammer. Bughazem and Anderson (1996, 2000) extended an earlier DVCM study by Anderson and Arfaie (1991) with Brunones unsteady friction model and found excellent agreement between theory and experiment. Shu (2003b) extended a two-phase ow model with Zielkes (1968) friction model. 3.9. Alternative methods for modeling column separation 3.9.1. Jordan algebraic method Jordan (1965, 1975) developed an analytical method for the treatment of distributed vaporous cavitation zones. He asserted that pressure waves cannot propagate through an established mixture of liquid and vapor and that therefore the SchnyderBergeron graphical method could not be used in this region. Jordan developed equations for distributed vaporous cavitation and for the penetration of liquid columns into cavitating regions. He thus calculated the duration of cavitation. Tarasevich (1975, 1997) developed a similar analytical method. 3.9.2. Finite element method Howlett (1971) modeled the liquid contained within a pipe system by means of solid beams without bending stiffness in a nite element (FEM) solution procedure. Watt et al. (1980), Bach and Spangenberg (1990), Jovic (1995) and Shu (2003a) applied the FEM to the classical water-hammer equations. Giesecke (1981) used the FEM and mentioned the discrete cavity model, but did not show any results. Axisa and Gibert (1982) and Schwirian (1982, 1984) employed the DVCM within the context of the FEM; gaps were allowed to form between the axial beam elements simulating the liquid as shown in Fig. 15. They compared numerical results obtained with and without cavitation.
ARTICLE IN PRESS
A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171 155
Fig. 15. DVCM in FEM context [adapted from Schwirian (1984, Fig. 4)]. Dynamic element with gap, which allows the formation of a cavity.
3.9.3. Other methods Mansour (1996) used the DVCM in combination with a special nite-difference scheme for the water-hammer equations. Fanelli (2000) presented unsteady pipe ow equations for a simple condenser cooling water system, including variable wave speed and friction losses due to the presence of gas bubbles, and equations for column separation. In addition, a number of boundary conditions were fully described, representing pumps, valves, condensers, siphons, surge tanks and nonreturn valves. 3.10. A comparison of models De Vries et al. (1971), Kalkwijk et al. (1972), Vreugdenhil et al. (1972) and Provoost (1976) compared Siemons (1966, 1967) shallow-water ow model with Kalkwijk and Kranenburgs (1971, 1973) bubble ow model. Experimental results were presented for a horizontal 1450 m test circuit. The construction of some large water supply pipelines in The Netherlands prompted these studies. Allowing the occurrence of pipeline cavitation was considered as an alternative to expensive water-hammer control devices such as surge tanks, air vessels and ywheels. The authors concluded that the results obtained from both computer programs exhibited adequate agreement with the experimental results for the horizontal test circuit. Provoost (1976) found that the shallow-water model did not reproduce the eld measurements for a pipeline system with two high points and the simplied bubble model of Kranenburg (1972, 1974b) was not suited to describe column separations at high points. As a result, the DGCM was used by De Vries (1973) and Provoost (1976). Simpson (1986), Bergant (1992) and Bergant and Simpson (1999a) compared numerical results from discrete vapor (DVCM), discrete gas (DGCM) and (generalized) interface vaporous cavitation models (GIVCM) with results of measurements. The principal source of discrepancies between the computed and measured results was found to originate from the method of physical description of distributed vaporous cavitation. Simpson and Bergant (1994a) and Bergant and Simpson (1994b) compared a number of cavitation models. They found that within the MOC the staggered grid is preferred above the rectangular grid, which actually comprises two independent staggered grids that may cause chessboard instability. Miwa et al. (1990) validated numerical results from the DVCM, the DGCM and a variable wave-speed model with consideration of gas release against experimental data. Dudlik (1999) and Dudlik et al. (2003) compared the results of commercial codes with new experimental data. Dudlik et al. (2000) compared the DVCM with a three-phase model that allowed for the calculation of sudden changes of gas content in the liquid. Shu (2003b) compared the DVCM with a two-phase model and with experimental data from Sanada et al. (1990).
4. Laboratory experiments and eld tests In all the experiments described in this section, unless otherwise stated, water hammer and column separation are initiated by the rapid closure of a valve and the test liquid is water. 4.1. Visualization of cavity formation Escande and Nougaro (1953) reported an early ow visualization study. They conducted column-separation experiments in a laboratory apparatus comprised of a 25 m long horizontal pipeline of 200 mm internal diameter. A
ARTICLE IN PRESS
156 A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171
transparent section was positioned next to the valve for ow visualization using a high-speed camera. A large vapor cavity was observed following the closure of the valve. Bunt (1953), Smirnov (1954), Kamel (1954) and Blind (1956) made similar observations; they presented photographs exhibiting the shape of the cavity at the valve during its growth and collapse. The vaporliquid interface sloped gently over quite a long distance. Duc (1959) investigated column separation, for three different piping congurations at a high point, due to the shutdown of pumps in a 1 km long discharge line at a eld installation. A clear piece of plexiglas pipe at the high point allowed a sequence of photographs to be taken, which exhibited the changes in the cavity during a column separation. These tests contradicted the generally accepted supposition of a complete separation of the liquid columns at an elevated point. Liquid remained at the separation point until return-ow lled the void. ONeill (1959) also concluded from a photographic study that intermediate cavities did not form over the entire pipe cross section, as ideally assumed in the theory; most cavities appeared in the upper portion of the cross section only. Li and Walsh (1964), Baltzer (1967b) and Safwat (1972a) presented photographs of a discrete cavity at the downstream side of a closing valve. The appearance of tiny bubbles was observed across the whole cross-section, which extended along a large portion of the pipe. Tanahashi and Kasahara (1970) studied photographically the formation and collapse of a discrete cavity at the high point of a pumping system. Nonoshita et al. (1991, 1992, 1999) presented photographs of column separation in a laboratory draft tube of a water turbine following wicket gates closure. They studied the effect of the draft tube inlet swirl on column-separation events. The swirl ow generated gas release and subsequent attenuation of maximum pressures following cavity collapse. Bergant (1992) and Bergant and Simpson (1996) presented photographs of column separation in a sloping pipeline. A discrete vapor cavity and a distributed vaporous cavitation zone were observed following rapid valve closure. Dudlik et al. (1997, 1999) and Dudlik (1999) presented photographs of column separation in a 230 m long test rig with complex pipe geometry. Swafeld (19691970) and Kojima et al. (1984) employed photography to visualize column separation in liquids other than water (kerosene, mineral oil). Differences in the column-separation mechanism for different liquids were not observed. Yang et al. (1996) and Adam et al. (1998) used a 100 Hz nonintrusive electrical-capacitance tomography technique to observe cross-sectional images from changes in uid capacitance in a 32 m long pipeline of 42 mm diameter. Cavities forming at the top of the pipe were of approximately circular cross-section rather than having a horizontal lower surface. Jaworek and Krupa (2004) developed a similar technique, but working at 80 MHz. Dudlik et al. (1997, 1999) and Dudlik (1999) used the 1 kHz wire-mesh visualization technique developed by Prasser et al. (1998, 2001) which is based on the uid conductivity. The ow-disturbing intrusive effect of such a device was studied by Wangjiraniran et al. (2003). Tabei et al. (2003) added small amounts of the noble gases xenon and argon to initially degassed water to provoke light emission at bubble collapse. In this way, they were able to accurately determine the speed of a shock wave entering a region of bubbly ow. Their theoretical study showed that local temperatures up to 7000 K may occur when small bubbles collapse. Su et al. (2003) observed the strongest light signals when water with xenon was cooled to near its freezing point. Chakravarty et al. (2004) conducted tests with many different liquids. They found the best cavitation sonoluminescence in liquids with a low vapor pressure and a moderately high viscosity. Mechanical energy could be converted into visible light with an efciency as high as 1%.
4.2. Laboratory experiments with column separation At the end of 19th century Joukowsky (1900) conducted water-hammer and column-separation experiments in Aleksejews water supply system in Moscow. He had three pipelines of lengths and diameters {(L, D) ((320, 51); (320, 102); (325, 152)) (m, mm)}, each with a downstream valve, connected to the main supply pipeline of diameter 610 mm. Pressure records were taken at the valve and along the pipeline. Joukowsky was the rst to qualitatively explain the phenomenon of column separation (see Section 2.3). Langevin and Boullee (1928)as presented in Bergeron (1950)were probably the rst persons using piezoelectric pressure transducers for the measurement of column separation. Their record shown in Fig. 16 is typical for all later measurements. Mostowskys (1929) measurement in a 30 m long, 51 mm diameter pipe is shown in Fig. 7(c). Binnie and Thackrah (1951) tested a pipeline with an automatic air-inlet valve. When the air-valve was removed, a series of violent impacts took place due to cavity formation and collapse. Maximum recorded pressures were actually somewhat greater than the theoretical estimates. Lupton (1953), in discussing Binnie and Thackrahs (1951) paper, did not agree with the use of air valves because of the unpredictably high pressures that may result from the presence of air in a pipeline. Gayed and Kamel (1959) showed pressure histories obtained in a 2.6 m long pipe for different initial ow velocities.
ARTICLE IN PRESS
A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171 157
Fig. 16. Typical pressure history for repeated column separation and collapse [adapted from Bergeron (1950, Fig. 56)].
Carstens and Hagler (1964, 1966) considered water hammer in a pipeline transporting phosphate-ore slurry. Results were presented from both a laboratory model and a theoretical analysis that described the column-separation phenomenon. Jordan (1965) investigated column separation and distributed vaporous cavitation in a laboratory apparatus comprising a 202 m long upward-sloping pipeline of 52 mm diameter. Baltzer (1967a, b) found that the measured pressure rises in a coiled-pipe apparatus were appreciably smaller than those predicted by his shallow-water ow model. Tanahashi and Kasahara (1970) presented a comparison of experimental and analytical results of water hammer with column separations at the valve and at the pipe midpoint. Raiteri and Siccardis (1970) experimental records were later on used by De Bernardinis et al. (1975) to validate their mathematical models. Safwat (1972ac), Safwat and De Kluyver (1972) and Safwat and Van Den Polder (1973) performed measurements in 4046 m long horizontal plexiglas pipelines of 90 mm diameter with and without a condenser. Vreugdenhil et al. (1972) conducted tests in a 1450 m long horizontal pipeline of 100 mm diameter to which Provoost (1976) added two high points. Piga and Sambiago (1974) tested in a laboratory cooling water system with lengths 2530 m and diameter 100 mm, and Thorley and Chohan (1976) in a horizontal pipe of length 16 m and diameter 38 mm. Krivchenko et al. (1975) presented results of column-separation measurements in a draft tube of the Kaplan turbine test rig at MISI, Moscow. A large cavity formed after rapid closure of the wicket gates. Column separation rst occurred in the space between the guide vanes and the runner, the cavity then grew into the draft tube inlet cone. A large pressure peak occurred after cavity collapse. The authors provided measurements of the pressures under the turbine head cover and at the draft tube inlet, as well as measurements of the axial hydraulic force acting on the turbine runner. Nonoshita et al. (1991, 1992, 1999) performed similar experimental tests. Katz and Chai (1978) carried out experiments in 0.3 m long tubes of 5 or 6 mm diameter with 2 ms valve closures. They showed one diagram of column-separation duration as a function of initial velocity (ranging from 1 to 20 m/s). This is the only paper (known to the authors) describing a benecial use of column separation, namely as feedback mechanism in liquid oscillators that can be used for industrial cleaning. Their theory, based on rigid-column theory, included the effect of gas release. Kot and Youngdahl (1978a, b) gave a clear explanation of the DVCM and they used experiments in a 9 m long closed tube for its validation. Aga et al. (1980) measured cavitating oil ow in a 250 m long, 90 mm diameter, test rig. Van De Sande and Belde (1981) tested in a U-tube of 9 m length and 45 mm diameter. Graze and Horlacher (1983) built two horizontal test rigs (120 and 86 m, and diameters of 82 and 200 mm, respectively). The authors stressed the importance of using adequate pressure transducers and recording equipment to avoid unrealistic pressure spikes, which is discussed in Section 4.6. Gottlieb et al. (1981) recorded extremely high pressure peaks immediately upon the collapse of vapor cavities. These could be unrealistic because of the type of pressure transducer used. Fox and McGarry (1983) developed a test rig for the study of pressure transients in pipelines carrying volatile liquids (18 m length, 55 mm diameter). Martin (1983) tested in a coiled tube (102 m length, 13 mm diameter). He presented experimental results (Fig. 4) with limited, moderate and severe cavitation, the severity being dened in Section 2.7. Golia and Greco (1990) found excellent agreement between DVCM computations and Martins (1983) experimental data. Borga and De Almeida (1985) tested in a horizontal pipeline of 105 m length and studied the inuence of an in-line nonreturn valve on reducing the pressure pulses. Simpson (1986) studied the short-duration pressure peaks explained in Section 2.6. Measurements were done for eight levels of cavitation severity. The test rig consisted of a 36 m long upward-sloping pipeline of 19 mm diameter connecting two reservoirs. The difference in elevation between the reservoirs was 1 m. Paredes et al. (1987) and Carmona et al. (1987) constructed a 1460 m long and 104 mm diameter horizontal pipeline. The IAHR working group (Fanelli, 2000) has a data base with 13 experimental records of the laboratory tests by Paredes et al. (1987) with different initial ow velocities (0.350.82 m/s) and static heads (1680 m). Barbero and Ciaponi (1991) reported on 23 experiments performed in a 500 m long, 110 mm diameter, test circuit. Anderson et al.
ARTICLE IN PRESS
158 A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171
(1991) and Anderson and Arfaie (1991) discussed several aspects of the DVCM and they showed results of laboratory measurements for three levels of cavitation severity. The laboratory apparatus of Shinada and Kojima (1987, 1989, 1995) was a small-scale physical model of a hydraulic press. The oil-lled test pipe was 5 m long with a diameter of 19 mm. The test results were compared with the results of a single vapor-cavity model that included laminar unsteady friction and a dynamic equation of motion for the valve. Miwa et al. (1990) performed column separation tests in a 500 m long steel pipeline of 0.10 m diameter. Pressures were recorded at seven positions along the pipeline. Tijsseling and Fan (1991a, b, 1992) carried out measurements in a 4.5 m long, 52 mm diameter, closed pipe. Transients were generated through the impact of a solid steel rod at one end of the pipe, which led to very steep wave fronts (pressure rises in microseconds). Two complications encountered in the conventional reservoir-pipe-valve system were absent: the initial steady state pressure gradient and the nonlinear valve closure. The water was stored under pressure in a closed container so that the amount of free gas was negligible. Because of the time scale of the experiment (in milliseconds) there was no gas release and negligible wall friction. The experiment isolated vaporous cavitation (in combination with uidstructure interaction, FSI), the severity of which was regulated by the static pressure of the water and the impact velocity of the rod. The experimental results were used to validate the DVCM. Sayir and Hausler (1991) also performed cavitation experiments in a closed tube. The tube was made of transparent PVC, 20 m long and 110 mm in diameter. Bergant (1992), Bergant and Simpson (1995) and Simpson and Bergant (1996) performed a series of tests in an adjustable experimental apparatus comprising a 37 m long pipeline of 22 mm diameter and sloping at about 5%. The following quantities were varied: initial ow velocity, heads in the tanks, and position (downstream end, midpoint, upstream end) and closure time of the valve. Repeat tests and uncertainty analyses were performed. The results and documentation of all 116 measurements can be obtained through the School of Civil and Environmental Engineering at the University of Adelaide. It is difcult to directly measure ow rates in transient ows. Washio et al. (1996) developed a procedure to accurately deduce ow rates in transient laminar ow from two adjacent pressure transducers. The procedure is based on frequency-dependent friction theory (Zielke, 1968) and relies on a specially constructed differential amplier. The method was shown to work very well for oil column-separation in the test rig of Washio et al. (1994). The typical changes of ow rates during column separation and rejoining were successfully measured for the rst time (Fig. 17). Mitosek (1997, 1998, 2000) made cavitation measurements in plastic pipes. He recorded extremely high pressure peaks that might be attributed to the strain-gage type pressure transducers used (see Section 4.6). Greenshields and Leevers (2000) studied the brittle fracture behavior of plastic pipes. Their laboratory apparatus consisted of a closed vertically falling pipe where column separation occurred at the tapered bottom end. Surge and cavity collapse were able to fracture pipes with an initial defect. Nakagawa and Takenaka (1995), Lai et al. (2000) and Gale and Tiselj (2003) carried out tests with initial vapor voids of different air content. Nakagawa and Takenaka (1995) created an initial void by cooling down hot liquid in a closed
(a)
(b)
Fig. 17. Measured (rst and third column) and calculated (second and fourth column) results in column-separation experiment. The second row depicts ow rates [adapted from Washio et al. (1996, Fig. 8)].
ARTICLE IN PRESS
A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171 159
pipe. Lai et al. (2000) performed tests in a pipeline with one vertical branch. Deaerated water was the test liquid. The test results were used to validate their DGCM code, which had a polytropic gas law, but no gas release mechanism. 4.3. Laboratory experiments with gas release Keller and Zielke (1976) measured free gas variations subsequent to a rapid drop in pressure in a 32 m long pipe with a diameter of 125 mm that was connected to a cavitation tunnel. Wiggert and Sundquist (1979) conducted experiments using 129 and 295 m long coiled-tube apparatuses with a diameter of 25 mm. They investigated gas release during transients at different initial gas concentrations. The effects of gas release, cavitation nuclei and turbulence were studied. Martin (1981) used a plexiglas pipe apparatus (length 32 m, diameter 26 mm) where the water was saturated with injected air. Kazama (1983) performed experiments with watermethanol mixtures. A test-section of 2 m length contained an initial cavity of vapor and foam at its vertical closed end. Sudden valve opening created a pressure wave collapsing the cavity. The pressure at one location and the volume of remaining air were measured. Shinada (1994) studied experimentally and theoretically column separation with gas release. In a 2.5 m long, 19 mm diameter pipe he tested with saturated and deaerated oil. He measured air content, surface tension and diffusion rate. On the basis of experimental results, the proposed bubble-diffusion model allowed for gas release only at the column separation. 4.4. Laboratory experiments with uidstructure interaction The repeated collapse of column separations and the almost instantaneous pressure rises associated with them (see Fig. 16) form a severe load for pipelines and their supporting structures. Structural vibration is likely to occur. Fluidinduced structural motion, structure-induced uid motion and the underlying coupling mechanisms are commonly referred to as FSI. Most of the experimental researchers mentioned herein prevented unwanted FSI effects by rigidly anchoring their pipes. Fan and Tijsseling (1992), however, focused on the simultaneous occurrence of cavitation and FSI. They performed experiments in a closed pipe, the vibrating ends of which interacted with transient columnseparations. They observed distributed vaporous cavitation caused by a stress wave in the pipe wall. Tijsseling et al. (1996) investigated experimentally and theoretically column separation in an unrestrained one-elbow pipe system. Tijsseling and Vardy (2006) presented cavitation tests in an unrestrained one-branch system. Review papers by Tijsseling (1996) and Wiggert and Tijsseling (2001) provide information on combined cavitation/FSI models and on FSI in general. 4.5. Field tests Apelt (1956) measured pressures in the eld for water hammer in pump discharge lines that undergo liquid column separation. His investigations veried that water-hammer theory may be applied with condence, but that it could not account for the phenomenon of column separation. Richards (1956) also presented some eld test data for pumping plants. He contended that it was virtually impossible to analytically analyze such systems. OBrien (1956) took issue with this point and cited some extensive calculations that involved six or seven separate column-separations. Whiteman and Pearsall (1962) conducted eld tests at a pumping station on reux-valve characteristics and pressure rises after pump shutdown. Heavy ywheels were used to lengthen the rundown time, thereby preventing column separation. Duc (1959, 1965) measured column separation in a eld installation. In addition to recording pressure peaks due to the rejoining of liquid columns, the phenomena were observed visually. The results exhibited some narrow high-pressure peaks of short duration. These may have resulted from the type of pressure gage being used (see discussion in Section 4.6). Brown (1968) presented eld measurements of transients in two pump discharge lines with distinct knees in their proles. The pressures measured in the eld were greater than those predicted by the graphical method in the design stage. This was attributed to the presence of free and dissolved air in the liquid. He noted that the presence of air may result in large pressure surges and higher reverse speeds of pumps due to the prolonging of column separation. De Vries (1975a, b) carried out measurements in a thermal power plant cooling water system (pipeline at the upstream end of the condenser: length of 76 m and diameter of 1.8 m; pipeline at the downstream end of the condenser: length of 182 m and diameter of 1.8 m). Provoost (1976) conducted an investigation of a 28 km long pipeline with a diameter of 1.8 m. The pipeline was carefully deaerated before each test. Six cases of pump shutdown were presented. Sharp (1977) presented eld measurements in a 12 km long, 250 mm diameter, pipeline. Siccardi (1979) presented a
ARTICLE IN PRESS
160 A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171
comprehensive summary of measurements and computations of water hammer and column separation in a cooling water system at La Casella thermal power plant (pipeline at the upstream end of the condenser: length of 440 m and diameter of 2 m; pipeline at the downstream end of the condenser: length of 234 m and diameter of 1.6 m). The data base of the IAHR working group (Fanelli, 2000) contains data measured in the La Casella power plant; Daco and Meregalli (1981) provided 25 experimental records of pump start-up and pump power-failure at different operating conditions. De Almeida and Hipolito (1981), Jolas (1981), Enever (1983) and Yow et al. (1985) conducted extensive tests in the cooling water systems of thermal and nuclear power stations. Wang and Locher (1991) found surprisingly good agreement between simulations and eld data obtained in a 47 km long cross-country pipeline with a diameter of 0.84 m. Axworthy and Chabot (2004) took measurements in a 0.25 m diameter sewage main consisting of 720 m of ductile-iron pipes and 915 m of PVC pipes. 4.6. Selection of pressure transducers Care must be taken in the selection (and way of mounting) of pressure transducers to be used in cavitation measurements. Pressure transducers may easily be damaged during cavitation tests, because the nearby explosions and implosions of small bubbles are too severe a load (Chen and Israelachvili, 1991; Broos, 1993, p. 13). The natural frequency of the pressure transducer must be high compared to the frequencies in the measured signal. It is noted that exploding and imploding cavitation bubbles may lead to pressure signals with a frequency spectrum up to 1 MHz (Oldenziel and Teijema, 1976, p. 14). Graze and Horlacher (1983) and Simpson and Bergant (1991, 1994b) reported unrealistic pressure spikes and oscillations for inductive and strain-gage type pressure transducers and attributed these to the fact that the natural frequency of the transducer was too low. Borga and De Almeida (1985) concluded that the pressure transducer should be a sophisticated low-inertia instrument. Sayir and Hausler (1991) had problems with a too low natural frequency of piezoelectric transducers. Le et al. (1989, p. 3) developed their own special transducers with a natural frequency of 1.7 MHz to overcome this problem. Arndt et al. (1995) described how to make your own pressure transducer. Mitosek (1997, 1998, 2000) found unrealistically high-pressure peaks with strain-gage type pressure transducers. Greenshields and Leevers (2000) discussed the frequencies in measured signals in relation to the natural frequencies of pipe-wall and piezoelectric pressure transducer, noting that the latter are reduced by added liquid mass.
5. Conclusion During the 20th century there was considerable research conducted into water hammer with column separation. This report attempts to span all of the signicant research that has been carried out during this period. The occurrence of low pressures and associated column separation during water-hammer events has been a concern for much of the 20th century in the design of pipe systems for distribution and cooling. The closure of a valve or shutdown of a pump may cause pressures so low that the liquid will cavitate. The collapse of vapor cavities and rejoinder of water columns can generatenearly instantaneouslyextremely large pressure that may cause signicant damage or ultimately failure of the pipe system. As early as 1900, Joukowsky had identied the physical occurrence of column separation. The 1930s produced the rst mathematical models of vapor cavity formation and collapse based on the graphical method. The identication of the various physical attributes of column separation occurred in the mid-20th century (distributed or vaporous cavitation in the 1930s; intermediate vapor cavities in the 1950s). These both led to a better physical understanding of the process of column separation and ultimately laid out the groundwork for the development of computer-based numerical models. The 1960s saw the development of the rst computer models of column separation within the framework of the method-of-characteristics solution of the water-hammer equations. A variety of alternative numerical models were developed from the late 1960s to the early 1990s. The most signicant models that have been developed include: the DVCM, the DGCM and the generalized interface vapor cavity model (GIVCM). The rst two are the easiest to implement. The DVCM is the most popular model used in currently available commercial computer codes for water hammer analysis. Advanced two-phase ow models are still under development. Numerous experimental studies have been carried out during the 20th century. These have beneted from the technological advances in measuring equipment (the introduction of the piezoelectric pressure transducer in the late 1920s was a milestone) and from the development of novel ow visualization technologies in the 1990s including highspeed video imaging, electrical-capacitance tomography and wire-mesh sensors. From all the validation tests presented in the literature it may be concluded that, despite its simplicity, the DVCM reproduces the essential features of transient cavitation. The versatility of the model has been demonstrated by the
ARTICLE IN PRESS
A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171 161
variety of pipe systems used in the tests. The major deciency of the model is the appearance of nonphysical oscillations in the results. The DGCM is better in this respect and therefore recommended in developing and revising industrial engineering water-hammer computer codes. The interface models give reliable results, but they still are quite complicated for general use.
Acknowledgments The Surge-Net project (www.surge-net.info) is supported by funding under the European Commissions Fifth Framework Growth Programme via Thematic Network Surge-Net contract reference: G1RT-CT-2002-05069. The authors of this paper are solely responsible for the content and it does not represent the opinion of the Commission. The Commission is not responsible for any use that might be made of data herein.
References
Adam, M.S., Yang, W.Q., Watson, R., 1998. A capacitance tomograc system for the measurement of void fraction in transient cavitation. IAHR Journal of Hydraulic Research 36, 707719. Aga, J., Karterud, T.J., Nielsen, T.K., 1980. Testing of transient ow and column separation in crude oil pipelines. In: Proceedings of the Third International Conference on Pressure Surges. BHRA, Canterbury, UK, pp. 113126. Akagawa, K., Fujii, T., 1987. Development of research on water hammer phenomena in two phase ow. In: Proceedings of the ASMEJSME Joint Technology Conference. Honolulu, Hawaii, USA, pp. 333349. Allievi, L., 1902. Teoria generale del moto perturbato dellacqua nei tubi in pressione. (General theory of the variable motion of water ` in pressure conduits.). Annali della Societa degli Ingegneri ed Architetti Italiani, Milan, Italy (in Italian). (French translation by Allievi himself, in Revue de Mecanique, Paris, 1904; German translation by R. Dubs and V. Bataillard, Springer, Berlin, 1909.) Allievi, L., 1913. Teoria del colpo dariete. (Theory of water-hammer.). Atti del Collegio degli Ingegneri ed Architetti Italiani, Milan, Italy (in Italian). (French translation by D. Gaden, 1921, Dunod, Paris; English translation by E.E. Halmos, 1925, Riccardo Garroni, Rome.) Anderson, A., 1976. Menabreas note on waterhammer: 1858. ASCE Journal of the Hydraulics Division 102, 2939. Anderson, A., Arfaie, M., 1991. Variable waterhammer wavespeed in column separation. In: Proceedings of the International Meeting on Hydraulic Transients with Water Column Separation. 9th Round Table of the IAHR Group, Valencia, Spain, pp. 183199. Anderson, A., Sandoval-Pena, R., Arfaie, M., 1991. Column separation behaviour modes in a simple test rig. In: Proceedings of the International Meeting on Hydraulic Transients with Water Column Separation. 9th Round Table of the IAHR Group, Valencia, Spain, pp. 3350. Angus, R.W., 1935. Simple graphical solution for pressure rise in pipes and discharge lines. Journal of the Engineering Institute of Canada 18 (2), 7281, 264273. Angus, R.W., 1937. Water hammer in pipes, including those supplied by centrifugal pumps: graphical treatment. Proceedings of the Institution of Mechanical Engineers 136, 245331 (Also: Bulletin No. 152, University of Toronto, Canada, 1938). Apelt, C.J., 1956. Investigation of water-hammer at the University of Queensland. Journal of Institution of Engineers, Australia 28 (3), 7581. Arndt, R.E.A., 1981. Cavitation in uid machinery and hydraulic structures. Annual Review of Fluid Mechanics 13, 273328. Arndt, R.E.A., Ellis, C.R., Paul, S., 1995. Preliminary investigation of the use of air injection to mitigate cavitation erosion. ASME Journal of Fluids Engineering 117, 498504. Axisa, F., Gibert, R.J., 1982. Non-linear analysis of uid-structure coupled transients in piping systems using nite elements. ASMEPVP 63, Flow-Induced Vibration of Circular Cylindrical Structures, pp. 151165. Axworthy, D.H., Chabot, N., 2004. Pressure transients in a Canadian sewage force main. Canadian Journal of Civil Engineering 31, 10391050. Baasiri, M., Tullis, J.P., 1983. Air release during column separation. ASME Journal of Fluids Engineering 105, 113118. Bach, P., Spangenberg, S., 1990. A numerical method for simulation of liquid and gas ow in pipeline networks. Rohre Rohrleitungsbau Rohrleitungstransport 3R (international) 29, 185190. Baltzer, R.A., 1967a. A study of column separation accompanying transient ow of liquids in pipes. Ph.D. Thesis, The University of Michigan, Ann Arbor, Michigan, USA. Baltzer, R.A., 1967b. Column separation accompanying liquid transients in pipes. ASME Journal of Basic Engineering 89, 837846. Barbero, G., Ciaponi, C., 1991. Experimental validation of a discrete free gas model for numerical simulation of hydraulic transients with cavitation. In: Proceedings of the International Meeting on Hydraulic Transients with Water Column Separation. 9th Round Table of the IAHR Group, Valencia, Spain, pp. 5169. Bergant, A., 1992. Kavitacijski tok med prehodnimi rezimi v cevnih sistemih. (Transient cavitating ow in pipelines.). Ph.D. Thesis, University of Ljubljana, Ljubljana, Slovenia (in Slovene).
ARTICLE IN PRESS
162 A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171 Bergant, A., Simpson, A.R., 1992. Interface model for transient cavitating ow in pipelines. In: Bettess, R., Watts, J. (Eds.), Unsteady Flow and Fluid Transients. A.A. Balkema, Rotterdam, pp. 333342. Bergant, A., Simpson, A.R., 1994a. Estimating unsteady friction in transient cavitating pipe ow. In: Miller, D.S. (Ed.), Water Pipeline Systems. Mechanical Engineering Publications, London, pp. 316. Bergant, A., Simpson, A.R., 1994b. Range of validity of the discrete vapor cavity model for pipeline water-hammer. Report No. R99, Deptartment of Civil and Environmental Engineering, University of Adelaide, Adelaide, Australia. Bergant, A., Simpson, A.R., 1995. Water hammer and column separation measurements in an experimental apparatus. Report No. R128, Dept. of Civil and Environmental Engineering, University of Adelaide, Adelaide, Australia. Bergant, A., Simpson, A.R., 1996. Vizualizacija kavitacijskega toka med prehodnimi rezimi v cevnih sistemih. (Visualization of transient cavitating ow in piping systems.). Strojnis ki VestnikJournal of Mechanical Engineering 42 (12), 316 (in Slovene and in English). Bergant, A., Simpson, A.R., 1999a. Pipeline column separation ow regimes. ASCE Journal of Hydraulic Engineering 125, 835848. Bergant, A., Simpson, A.R., 1999b. Cavitation inception in pipeline column separation. In: Proceedings of the 28th IAHR Congress, Graz, Austria, CD-ROM, 7pp. Bergant, A., Simpson, A.R., V tkovsky, J., 2001. Developments in unsteady pipe ow friction modelling. IAHR Journal of Hydraulic Research 39, 249257. Bergeron, L., 1935. Etude des variations de regime dans les conduites deau. Solution graphique generale. (Study on the steady-state variations in water-lled conduits. General graphical solution.). Revue Generale de lHydraulique 1 (1), 1225 (in French). Bergeron, L., 1939. Discussion of Experiments and calculations on the resurge phase of water hammer by J.N. LeConte and discussion of Air chambers for discharge lines by L. Allievi. Transactions of the ASME 61, 441445. Bergeron, L., 1950. Du Coup de Belier en HydrauliqueAu Coup de Foudre en Electricite. (Waterhammer in hydraulics and wave surges in electricity.). Dunod, Paris (in French). (English translation by ASME Committee, Wiley, New York, 1961.) Beuthe, T.G., 1997. Review of two-phase water hammer. In: Proceedings of the 18th Canadian Nuclear Society Conference, Toronto, Canada, 20pp. Billings, A.W.K., Dodkin, O.K., Knapp, F., Santos, A., 1933. High head penstock design. In: Proceedings of ASME Water Hammer Symposium, Chicago, USA, pp. 2961. Binnie, A.M., Thackrah, D.G., 1951. Water hammer in a pumping main and its prevention. Proceedings of the Institution of Mechanical Engineers 165, 4352, 6573. Blind, H., 1956. Nichtstationare Stromungen in Unterwasserstollen. (Unsteady ow in tailrace tunnels.). In: Tolke, F. (Ed.), Veroffentlichungen zur Forschung der DruckstoXprobleme in Wasserkraftanlagen und Rohrleitungen, Second ed. Springer, Berlin, pp. 67108 (in German). Bonin, C.C., 1960. Water-hammer damage to Oigawa Power Station. ASME Journal of Engineering for Power 82, 111119. Borga, A.A., De Almeida, A.B., 1985. Transient cavitation modeling. In: Proceedings of the Seventh International Round Table on Hydraulic Transients with Water Column Separation. IAHR, Alvkarleby, Sweden. Brennen, C.E., 1995. Cavitation and Bubble Dynamics. Oxford University Press Inc., New York. Broos, Ph., 1993. Knocking and gravel in vibrating pipelines predictable. Delft Outlook 93(2), 1013. ISSN: 0926-7212. Brown, G.O., 2002. The history of the Darcy-Weisbach equation for pipe ow resistance. In: Proceedings of the 150th Anniversary Conference of ASCE. Washington DC, USA, pp. 3443. Brown, R.J., 1968. Water-column separation at two pumping plants. ASME Journal of Basic Engineering 90, 521531. Brown, S.W.J., Williams, P.R., 2000. The tensile behaviour of elastic liquids under dynamic stressing. Journal of Non-Newtonian Fluid Mechanics 90, 111. Brunone, B., Greco, M., 1990. Un modello per la ricostruzione di fenomeni di colpo dariete anche in presenza di cavitazione. (A model for the reconstruction of the phenomena of waterhammer also in the presence of cavitation.). In: Proceedings of the 22nd Italian Congress of Hydraulics and Hydraulic Construction, Cosenza, Italy, vol. 4, pp. 147160 (in Italian). Brunone, B., Golia, U.M., Greco, M., 1991. Modelling of fast transients by numerical methods. In: Proceedings of the International Meeting on Hydraulic Transients with Column Separation. 9th Round Table, IAHR, Valencia, Spain, pp. 215222. Bughazem, M.B., Anderson, A., 1996. Problems with simple models for damping in unsteady ow. In: Boldy, A. (Ed.), Pressure Surges and Fluid Transients in Pipelines and Open Channels. Mechanical Engineering Publications, Burry St. Edmunds, pp. 537548. Bughazem, M.B., Anderson, A., 2000. Investigation of an unsteady friction model for waterhammer and column separation. In: Anderson, A. (Ed.), Pressure Surges. Safe Design and Operation of Industrial Pipe Systems. Professional Engineering Publishing Ltd., Burry St. Edmunds, pp. 483498. Bunt, E.A., 1953. Preliminary study of valve cavitation in pipelines. Journal of the South African Institution of Mechanical Engineers 2, 235260. Campbell, I.J., Pitcher, A.S., 1958. Shock waves in a liquid containing gas bubbles. Proceedings of the Royal Society, Series A 243, 535545. Capozza, A., 1986. TRANSID: a computer program for the study of hydraulic transients (involving cavitation and waterhammer). Energia Nucleare 3 (2), 105122. Carmona, R., Sanchez, A., Sanchez, J.L., 1987. Experimental relation between the highest transient pressure and the severity of water column separation. In: Proceedings of the Eighth International Round Table on Hydraulic Transients in Power Stations. IAHR, Madeira, Portugal Paper D2.
ARTICLE IN PRESS
A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171 163 Carmona, R., Sanchez, J.L., Carmona, L., 1988. A simplied procedure to evaluate liquid column separation phenomena. Water Power and Dam Construction 40 (12), 4246. Carpenter, R.C., 1894. Some experiments on the effect of water hammer. The Engineering Record 30, 173175. Carpenter, R.C., Barraclough, S.H., 1894. Some experiments on the effect of water hammer. Transactions of the ASME 15, 510535. Carstens, M.R., Hagler, T.W., 1964. Water hammer resulting from cavitating pumps. ASCE Journal of the Hydraulics Division 90 (HY6), 161184. Carstens, M.R., Hagler, T.W., 1966. Closure to Water hammer resulting from cavitating pumps. ASCE Journal of the Hydraulics Division 92 (HY2), 357359. Carstens, M.R., Roller, J.E., 1959. Boundary-shear stress in unsteady turbulent pipe ow. ASCE Journal of the Hydraulics Division 85 (HY2), 6781. Chakravarty, A., Georghiou, T., Phillipson, T.E., Walton, A.J., 2004. Stable sonoluminescence within a water hammer tube. Physical Review E 69 (066317), 18. Chaudhry, M.H., Bhallamudi, S.M., Martin, C.S., Naghash, M., 1990. Analysis of transient pressures in bubbly, homogeneous, gasliquid mixtures. ASME Journal of Fluids Engineering 112, 225231. Chen, Y.L., Israelachvili, J., 1991. New mechanism of cavitation damage. Science 252, 11571160. Contractor, D.N., 1965. The reection of water hammer pressure waves from minor losses. ASME Journal of Basic Engineering 87, 445452. Daco, E., Meregalli, E., 1981. Data bank of in-situ tests performed on the condenser cooling circuits of La Casella thermal power plant (Piacenza, Italy). In: Proceedings of the Fifth International Round Table on Hydraulic Transients with Water Column Separation. IAHR, Obernach, Germany. De Almeida, A.B., 1983. Cavitation and water-column separation by the method of characteristics. In: Proceedings of the Sixth International Symposium on Water Column Separation. IAHR, Gloucester, UK. De Almeida, A.B., 1987. Some contributions to the transient macrocavitations knowledge and modelling (19781987). In: Proceedings of the Eighth International Round Table on Hydraulic Transients in Power Stations. IAHR, Madeira, Portugal (Paper D1). De Almeida, A.B., 1991. Accidents and incidents: a harmful/powerful way to develop expertise on pressure transients. In: Proceedings of the International Meeting on Hydraulic Transients with Water Column Separation. 9th Round Table of the IAHR Group, Valencia, Spain, pp. 379401. De Almeida, A.B., Hipolito, H.J., 1981. Transient macrocavitation simulation for hydraulic pumping systems. In: Proceedings of the Fifth International Round Table on Hydraulic Transients with Water Column Separation. IAHR, Obernach, Germany. De Bernardinis, B., Federici, G., Siccardi, F., 1975. Transient with liquid column separation: numerical evaluation and comparison with experimental results. LEnergia Elettrica 52, 471477. De Haller, P., Bedue, A., 1951. The break-away of water columns as a result of negative pressure shocks. Sulzer Technical Review 43 (4), 1825. De Vries, A.H., 1972. Cavitatie door waterslag in horizontale leidingen. (Cavitation due to waterhammer in horizontal pipelines.). Delft Hydraulics Laboratory, Report M 1116, Delft, The Netherlands (in Dutch). De Vries, A.H., 1973. Cavitatie door waterslag in horizontale leidingen met enige hoge punten. (Cavitation due to waterhammer in horizontal pipelines with several high points.). Delft Hydraulics Laboratory, Report M 1152, Delft, The Netherlands (in Dutch). De Vries, A.H., 1974. Hydraulic aspects of cooling water systems for thermal power plants. In: Transactions of the Seventh IAHR Symposium, Section for Hydraulic Machinery, Equipment and Cavitation, Vienna, Austria. Also: Delft Hydraulics Laboratory, Publication No. 136. De Vries, A.H., 1975a. Research on cavitation due to waterhammer in the Netherlands. In: Proceedings of the Second International Round Table on Hydraulic Transients with Water Column Separation. IAHR, Vallombrosa, Italy. De Vries, A.H., 1975b. Research on cavitation due to water hammer in the Netherlands. LEnergia Elettrica 52, 478485. De Vries, A.H., Vreugdenhil, C.B., Kranenburg, C., 1971. Cavitatie door waterslag in horizontale leidingen. (Cavitation due to waterhammer in horizontal pipelines.). Delft Hydraulics Laboratory, Report S 143, Delft, The Netherlands (in Dutch). Dijkman, H.K.M., 1968. Cavitatieverschijnselen in een horizontale leiding met water waarin gas is opgelost. (Cavitation phenomena in a horizontal conduit lled with water in which gas is dissolved.). Delft Hydraulics Laboratory, Report S 103-III, Delft, The Netherlands (in Dutch). Dijkman, H.K.M., Vreugdenhil, C.B., 1969. The effect of dissolved gas on cavitation in horizontal pipe-lines. IAHR Journal of Hydraulic Research 7, 301314 (Also: Delft Hydraulics Laboratory, Publication No. 70). Duc, J., 1959. Negative pressure phenomena in a pump pipeline. Sulzer Technical Review 41 (3), 311. Duc, J., 1965. Negative pressure phenomena in pump pipe lines. In: Proceedings of the International Symposium on Water Hammer in Pumped Storage Projects. ASME, Chicago, USA, pp. 154167. Dudlik, A., 1999. Vergleichende Untersuchungen zur Beschreibung von transienten Stromungsvorgangen in Rohrleitungen. (Comparing investigations on the description of transient ow processes in pipelines.). Ph.D. Thesis, Universitat Dortmund, UMSICHTSchriftenreihe Band 20, ISBN: 3-8167-5550-X, Fraunhofer IRB Verlag, Dortmund (in German). Dudlik, A., Schluter, S., Hoyer, N., Prasser, H.M., 1997. Transiente Stromungsvorgange in Rohrleitungen. (Transient ow processes in pipelines.). VDI-Bericht Nr. 1359, pp. 353368, VDI-Verlag, Dusseldorf (in German). Dudlik, A., Schluter, S., Weinspach, P.M., 1999. Water hammer and cavitation in long-distance energy pipeworksmeasurement, simulation and prevention. The inuence of valve characteristics and ow-rate dependant tting positions on uiddynamic
ARTICLE IN PRESS
164 A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171 calculations of pipeworks. In: Proceedings of the Seventh International Symposium on District Heating and Cooling, Lund, Sweden, Paper 29. Dudlik, A., Schluter, S., Hoyer, N., Prasser, H.-M., 2000. Pressure surgesexperimental investigations and calculations with software codes using different physical models and assumptions. In: Anderson, A. (Ed.), Pressure Surges. Safe Design and Operation of Industrial Pipe Systems. Professional Engineering Publishing Ltd., Burry St. Edmunds, pp. 279289. Dudlik, A., Schonfeld, S.B.H., Hagemann, O., Fahlenkamp, H., 2003. Water hammer and cavitational hammer in process plant pipe systems. Kerntechnik 68 (3), 9196. Enever, K.J., 1972. Surge pressures in a gasliquid mixture with low gas content. In: Proceedings of the First International Conference on Pressure Surges. BHRA, Canterbury, UK (Paper C1, 11pp). Enever, K.J., 1983. Comparison of theory and practice at Grain power station. In: Proceedings of the Sixth International Round Table on Hydraulic Transients with Water Column Separation. IAHR, Gloucester, UK. Escande, L., 1962. Arret instantane du debit dune conduite forcee avec cavitations. (Instantaneous closing of the discharge in a penstock with cavitations.). In: Proceedings of the IAHR Symposium on Cavitation and Hydraulic Machinery, Sendai, Japan, pp. 113124 (in French). Escande, L., Nougaro, J., 1953. Etude theorique et experimentale du fonctionnement en charge des canaux de fuite en labsence de cheminee dequilibre. (Theoretical and experimental study of the operation of tailraces in the absence of a surge shaft.). La Houille Blanche 8, 607639 (in French). Evans, E.P., Sage, P.V., 1983. Surge analysis of a large gravity pipeline. In: Proceedings of the Fourth International Conference on Pressure Surges. BHRA, Bath, UK, pp. 447460. Ewing, D.J.F., 1980. Allowing for free air in waterhammer analysis. In: Proceedings of the Third International Conference on Pressure Surges. BHRA, Canterbury, UK, pp. 127146. Fan, D., Tijsseling, A., 1992. Fluidstructure interaction with cavitation in transient pipe ows. ASME Journal of Fluids Engineering 114, 268274. Fanelli, M., 2000. Hydraulic Transients with Water Column Separation, IAHR Working Group 19711991 Synthesis Report, IAHR, Delft and ENEL-CRIS, Milan. Fox, P., McGarry, M., 1983. Pressure transients in pipelines transporting volatile liquids. In: Proceedings of the Fourth International Conference on Pressure Surges. BHRA, Bath, UK, pp. 3957. Frizell, J.P., 1898. Pressures resulting from changes of velocity of water in pipes. Transactions of the ASCE 39, Paper No. 819, pp. 118. Galante, C., Pointer, S., 2002. Catastrophic water hammer in a steam dead leg. IChemE Loss Prevention Bulletin 167, 1620. Gale, J., Tiselj, I., 2003. Modelling of cold water hammer with WAHA code. In: Proceedings of the International ConferenceNuclear Energy for New Europe, Portoroz, Slovenia, pp. 214.1214.8. Gayed, Y.K., Kamel, M.Y.M., 1959. Mechanics of secondary water-hammer waves. Proceedings of the Institution of Mechanical Engineers 173 (26), 675686. Geurst, J.A., 1985. Virtual mass in two-phase bubbly ow. Physica A 129, 233261. Gibson, A.H., 1908. Water Hammer in Hydraulic Pipelines. Archibald Constable and Co. Ltd., London. Gibson, N.R., 19191920. Pressures in penstocks caused by gradual closing of turbine gates. Transactions of ASME 83, 707775. Giesecke, H.D., 1981. Calculation of piping response to uid transients including effects of uid/structure interaction. In: Transactions of SMiRT6, Paris, France, Paper B4/4. Golia, U.M., 1990. Sulla valutazione delle forze resistenti nel colpo dariete. (On the valuation of friction forces in waterhammer.). Report No. 639, Dept. of Hydraulic and Environmental Engineering, University of Naples, Naples, Italy (in Italian). Golia, U.M., Greco, M., 1990. Cavitation during water-hammer: quick closure of a downstream valve. In: Proceedings of the Third International Conference on Hydraulic Engineering Software, Hydrosoft 90, Boston, USA, pp. 121129. Gomez, D.M., Langevin, A., 1937. La piezographie directe et instantanee. (Direct and instantaneous piezography.). Actualites Scientiques et Industrielles 512, Recherches dHemodynamique et Cardiologie. Hermann et Cie, in Paris (in French). Gottlieb, L., Larns, G., Vasehus, J., 1981. Transient cavitation in pipelineslaboratory tests and numerical calculations. In: Proceedings of the Fifth International Symposium on Water Column Separation. IAHR, Obernach, Germany, pp. 487508. Graze, H.R., Horlacher, H.B., 1983. Pressure transients following the collapse of vapor cavities. In: Proceedings of the Sixth International Symposium on Hydraulic Transients in Power Stations. IAHR, Gloucester, UK. Greenshields, C.J., Leevers, P.S., 2000. Failure of plastic water pipes by surge and cavitation. In: Proceedings of the Eighth International Conference on Pressure Surges. BHR Group, The Hague, The Netherlands, pp. 377391. Grifth, P., 1987. Characteristics of bubbly ow. In: Hewitt, G.F., Delhaye, J.J., Zuber, N. (Eds.), Multiphase Science and Technology 3. Hemisphere Publishing Corporation, Washington, pp. 155162. Hadj-Taieb, E., Lili, T., 1998. Transient ow of homogeneous gasliquid mixtures in pipelines. International Journal of Numerical Methods for Heat & Fluid Flow 8, 350368. Hadj-Taieb, E., Lili, T., 1999. The numerical solution of the transient two-phase ow in rigid pipelines. International Journal for Numerical Methods in Fluids 29, 501514. Hadj-Taieb, E., Lili, T., 2000. Validation of hyperbolic model for water-hammer in deformable pipes. ASME Journal of Fluids Engineering 122, 5764. Hadj-Taieb, E., Lili, T., 2001. Ecoulements transitoires dans les conduites deformables avec degazage de lair dissous. (Transient ows in plastic pipes with dissolved air degasication.). La Houille Blanche 56, 99107 (in French).
ARTICLE IN PRESS
A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171 165 Hager, W.H., 2001. Swiss contribution to water hammer theory. IAHR Journal of Hydraulic Research 39, 310. Hammitt, F.G., 1980. Cavitation and Multiphase Flow Phenomena. McGraw-Hill, New York. Hatwin, P., Henwood, G.A., Huber, R., 1970. On the collapse of water vapor cavities in a bubble analogue apparatus. Chemical Engineering Science 25, 11971209. Heath, W.E., 1962. Vapor-cavity formation in a pipe after valve closure. M.Sc. Thesis, Georgia Institute of Technology, Atlanta, USA. Hino, M., Sawamoto, M., Takasu, S., 1977. Study on the transition to turbulence and frictional coefcient in an oscillatory pipe ow. Transactions of JSCE 9, 282284. Hogg, T.H., Traill, J.J., 1926. Discussion of Speed changes of hydraulic turbines for sudden changes of load by E.B. Strowger and S.L. Kerr. Transactions of the ASME 48, 252257. Howlett, J.T., 1971. Applications of NASTRAN to Coupled Structural and Hydrodynamic Responses in Aircraft Hydraulic Systems. NASTRAN: Users Experiences, NASA TM X-2378. National Aeronautics and Space Administration, Washington, DC, pp. 407419. Jaeger, C., 1948. Water hammer effects in power conduits. (4 Parts). Civil Engineering and Public Works Review 23, 7476, 138140, 192194, 244246. Jaeger, C., Kerr, L.S., Wylie, E.B., 1965. Selected Bibliography. In: Proceedings of the International Symposium on Water Hammer in Pumped Storage Projects. ASME Winter Annual Meeting, Chicago, USA, pp. 233241. Jaworek, A., Krupa, A., 2004. Gas/liquid ratio measurements by rf resonance capacitance sensor. Sensors and Actuators A 113, 133139. Jolas, C., 1981. Hydraulic transients in closed cooling water systems. In: Proceedings of the Fifth International Round Table on Hydraulic Transients with Water Column Separation. IAHR, Obernach, Germany. Jordan, V., 1961. Vpliv povratne lopute na hidravlicni udar pri izklopu crpalke. (The inuence of check valves on water hammer at pump failure.). Strojnis ki Vestnik 7 (4, 5), 1921 (in Slovene). Jordan, V., 1965. Odredivanje hidraulickog udarca pri iskljucenju crpke bez ublazivaca udarca pod uslovima raskidanja vodenog stuba. (Prediction of water hammer at pump failure without surge protection under water column separation conditions.). Ph.D. Thesis, University of Belgrade, Belgrade, Yugoslavia (in Serbian). Jordan, V., 1975. Neue Ermittlungen uber den DruckstoX in Pumpenleitungen ohne DruckstoXdampfung. (New investigations on waterhammer in pump pipelines without waterhammer-damping.). GWFWasser/Abwasser 116 (12), 540548 (in German). Joukowsky, N., 1900. Uber den hydraulischen Stoss in Wasserleitungsrohren. (On the hydraulic hammer in water supply pipes.). Memoires de lAcademie Imperiale des Sciences de St.-Petersbourg, Series 8, vol. 9, No. 5 (in German). (English translation, partly, by Simin, 1904.) Jovic, V., 1995. Finite elements and the method of characteristics applied to water hammer modelling. Engineering Modelling 8 (34), 5158. Kalkwijk, J.P.Th., Kranenburg, C., 1971. Cavitation in horizontal pipelines due to water hammer. ASCE Journal of the Hydraulics Division 97(HY10), 15851605. (Discussion by Vreugdenhil, C.B., De Vries, A.H., in 98(HY9), 17231725.). Also part of: Delft University of Technology, Dept. of Civil Engineering, Laboratory of Fluid Mechanics, Report No. B/71/3. Kalkwijk, J.P.Th., Kranenburg, C., 1973. Closure to Cavitation in horizontal pipelines due to water hammer. ASCE Journal of the Hydraulics Division 99 (HY3), 529530. Kalkwijk, J.P.Th., Kranenburg, C., Vreugdenhil, C.B., De Vries, A.H., 1972. Cavitation caused by water hammer in horizontal pipelines. Delft Hydraulics Laboratory, Publication No. 97, Delft, The Netherlands. Kamel, M.Y.M., 1954. Cavitation produced by fast velocity changes in pipes. M.Sc. Thesis, Cairo University, Faculty of Engineering, Cairo, Egypt. Katz, S., Chai, D., 1978. Column separation as a feedback mechanism in liquid oscillators. Fluidics Quarterly 10 (1), 4359. Kazama, T., 1983. Column separation and waterhammer in binary mixture. M.Sc. Thesis, Massachusetts Institute of Technology, Dept. of Mechanical Engineering, Boston, USA. Keller, A., Zielke, W., 1976. Variation of free gas content in water during pressure uctuations. In: Proceedings of the Second Conference on Pressure Surges. BHRA, London, UK. Kephart, J.T., Davis, K., 1961. Pressure surges following water-column separation. ASME Journal of Basic Engineering 83, 456460. Kessal, M., Amaouche, M., 2001. Numerical simulation of transient vaporous and gaseous cavitation in pipelines. International Journal of Numerical Methods for Heat & Fluid Flow 11, 121137. Kloosterman, A.H., Kranenburg, C., De Vries, A.H., Wijdieks, J., 1973. Niet-stationaire hydrodynamische verschijnselen in transportleidingen voor vloeistoffen. (Non-stationary hydrodynamic phenomena in uid-conveying conduits.). Bijdragen KIvI Symposium, Delft, June 1972. Also: Delft Hydraulics Laboratory, Publication No. 109N, May 1973, Delft, The Netherlands (in Dutch). Knapp, F., 1937a. Discussion of Water hammer in pipes, including those supplied by centrifugal pumps: graphical treatment by R.W. Angus. Proceedings of the Institution of Mechanical Engineers 136, 304309. Knapp, F., 1937b. Operation of emergency shutoff valves in pipelines. Transactions of the ASME 59, 679682. Knapp, F., 1939. Discussion of Experiments and calculations on the resurge phase of water hammer by J.N. LeConte. Transactions of the ASME 61, 440441. Knapp, R.T., Daily, J.W., Hammitt, F.G., 1970. Cavitation. McGraw-Hill, New York. Kobori, T., Yokoyama, S., Miyashiro, H., 1955. Propagation velocity of pressure wave in pipe line. Hitachi Hyoron 37 (10).
ARTICLE IN PRESS
166 A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171 Kojima, E., Shinada, M., Shindo, K., 1984. Fluid transient phenomena accompanied with column separation in uid power pipeline. Bulletin of JSME 27 (233), 24212429. Korteweg, D.J., 1878. Ueber die Fortpanzungsgeschwindigkeit des Schalles in elastischen Rohren. (On the velocity of propagation of sound in elastic tubes.). Annalen der Physik und Chemie, New Series 5, 525542 (in German). Kot, C.A., Youngdahl, C.K., 1978a. Transient cavitation effects in uid piping systems. Nuclear Engineering and Design 45, 93100. Kot, C.A., Youngdahl, C.K., 1978b. The analysis of uid transients in piping systems, including the effects of cavitation. In: Proceedings of the Winter Annual Meeting of the ASME. Fluid Transients and Acoustics in the Power Industry, San Francisco, USA, pp. 4552. Kottmann, A., 1989. Vorgange beim AbreiXen einer Wassersaule. (Phenomena during breakaway of a water column.). Rohre Rohrleitungsbau Rohrleitungstransport 3R (international) 28, 106110 (in German). Kranenburg, C., 1972. The effect of free gas on cavitation in pipelines induced by water hammer. In: Proceedings of the First International Conference on Pressure Surges. BHRA, Canterbury, UK, pp. 4152. Kranenburg, C., 1974a. Transient cavitation in pipelines. Ph.D. Thesis, Delft University of Technology, Dept. of Civil Engineering, Laboratory of Fluid Mechanics, Delft, The Netherlands. Also: Communications on Hydraulics, Delft University of Technology, Dept. of Civil Engineering, Report No. 73-2, 1973. Kranenburg, C., 1974b. Gas release during transient cavitation in pipes. ASCE Journal of the Hydraulics Division 100(HY10), 13831398. Also part of: The effect of gas release on column separation. Communications on Hydraulics, Delft University of Technology, Dept. of Civil Engineering, Report No. 74-3. Krivchenko, G.I., Arshenevskii, N.N., Kvjatkovskaja, E.V., Klabukov, V.M., 1975. !nspomexahn#eckne &epexoshFIe &poueccFI b !nspoDhepletn#ecknx Yctahobkax. (Hydromechanical Propagation Processes in Hydropower Installations.). Energiya, Moscow (in Russian). Lai, A., Hau, K.F., Noghrehkar, R., Swartz, R., 2000. Investigation of waterhammer in piping networks with voids containing noncondensable gas. Nuclear Engineering and Design 197, 6174. Lai, C., 1961. A study of waterhammer including effect of hydraulic losses. Ph.D. Thesis, The University of Michigan, Ann Arbor, Michigan, USA. Langevin, A., Boullee, M., 1928. In: Bulletin de lUnion Technique du Batiment, p. 81. (Reference in (Gomez and Langevin, 1937) and (Bergeron, 1939, 1950)) (in French.) Le, Q., Franc, J.P., Michel, J.M., 1989. Analyse experimentale des pressions statique et dynamique dans la zone de fermeture des poches de cavitation partielle. (Experimental analysis of static and dynamic pressures in the region of pockets of collapse of partial cavitation.). Journees de lHydraulique, Societe Hydrotechnique de France, Machines Hydrauliques, Conception et Exploitation, 20(1), l.15.11.15.9, Lyon, France (in French). LeConte, J.N., 1937. Experiments and calculations on the resurge phase of water hammer. Transactions of the ASME 59, Paper HYD59-12, pp. 691694. Lee, I.Y., Kitagawa, A., Takenaka, T., 1985. On the transient behaviour of oil ow under negative pressure. Bulletin of JSME 28 (240), 10971104. Li, S.C., 2000. Cavitation of Hydraulic Machinery. Imperial College Press, London. Li, W.H., 1962. Mechanics of pipe ow following column separation. ASCE Journal of the Engineering Mechanics Division 88 (4), 97118. Li, W.H., 1963. Mechanics of pipe ow following column separation. Transactions of the ASCE 128 (Part I), 12331254. Li, W.H., Walsh, J.P., 1964. Pressure generated by cavitation in a pipe. ASCE Journal of the Engineering Mechanics Division 90 (6), 113133. Liou, J.C.P., 1999. Numerical properties of the discrete gas cavity model for transients. In: Proceedings of the Third ASME-JSME Joint Fluids Engineering Conference, Symposium S-290 Water Hammer, San Francisco, USA, ASME-FED 248, Paper FEDSM99-6901, 9pp. Liou, J.C.P., 2000. Numerical properties of the discrete gas cavity model for transients. ASME Journal of Fluids Engineering 122, 636639. List, E.J., Burnam, J., Solbrig, R., Hogatt, J., 1999. Vapor cavity formation and collapse: eld evidence for major pipeline damage. In: Proceedings of the Third ASME-JSME Joint Fluids Engineering Conference, Symposium S-290 Water Hammer, San Francisco, USA, July 1999, ASME-FED 248, Paper FEDSM99-6886, 7pp. Lupton, H.R., 1953. Graphical analysis of pressure surges in pumping systems. Journal of the Institution of Water Engineers 7, 87125. Mansour, S.G.S., 1996. The application of a space-compact high-order implicit scheme for cavitation computation in waterhammer simulation. In: Proceedings of the Second International Conference on Hydroinformatics, Zurich, Switzerland, pp. 315319. Marsden, N.J., Fox, J.A., 1976. An alternative approach to the problem of column separation in an elevated section of pipeline. In: Proceedings of the Second International Conference on Pressure Surges. BHRA, London, UK, pp. 113. Martin, C.S., 1973. Status of uid transients in Western Europe and the United Kingdom. Report on laboratory visits by Freeman scholar. ASME Journal of Fluids Engineering 95, 301318. Martin, C.S., 1981. Gas release in transient pipe ow. In: Proceedings of the Fifth International Round Table on Hydraulic Transients with Water Column Separation. IAHR, Obernach, Germany. Martin, C.S., 1983. Experimental investigation of column separation with rapid closure of downstream valve. In: Proceedings of the Fourth International Conference on Pressure Surges. BHRA, Bath, UK, pp. 7788.
ARTICLE IN PRESS
A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171 167 Martin, C.S., 1989. Private communication to C.S.W. Lavooij at Delft Hydraulics, Delft, The Netherlands. Martin, C.S., Padmanabhan, M., 1975. The effect of free gases on pressure transients. LEnergia Elettrica 52, 262267. Menabrea, L.-F., 1858. Note sur les effets du choc de leau dans les conduites. (Note on the effects of water shock in conduits.). Comptes Rendus Hebdomadaires des Seances de lAcademie des Sciences (Paris) 47, 221224 (in French). (English translation by Anderson, 1976.). Menabrea, L.-F., 1862. Note sur leffet du choc de leau dans les conduites. (Note on the effect of water shock in conduits.). Annales du Genie Civil (Paris) I, 269275 (in French). Michaud, J., 1878. Coup de belier dans les conduites: etude des moyens employes pour en attenuer les effects. (Water hammer in conduits: study of the means used for diminishing the effects.). Bulletin de la Societe Vaudoise des Ingenieurs et des Architectes, Lausanne, Switzerland 4 (3/4), 5664/6577 (in French). Mitosek, M., 1997. Study of cavitation due to water hammer in plastic pipes. Plastics, Rubber and Composites Processing and Applications 26 (7), 324329. Mitosek, M., 1998. Cavitation in pipes at water hammer conditions. Rohre Rohrleitungsbau Rohrleitungstransport 3R (international) 37, 806810. Mitosek, M., 2000. Study of transient vapor cavitation in series pipe systems. ASCE Journal of Hydraulic Engineering 126, 904911. Miwa, T., Sano, M., Yamamoto, K., 1990. Experimental studies on water hammer phenomenon including column separation. Water Supply 8, 430438. Moores, R.G., 1953. Discussion of Graphical analysis of pressure surge in pumping systems by H.R. Lupton. Journal of the Institution of Water Engineers 7, 139140. Moshnin, L.F., 1961. YkaDahnr fo Dawnte bosobosob ot lnspabjn#eckolo ysapa. (Instructions for the protection of conduits against waterhammer.). USSR Scientic Research Institute (WNII) Wodgeo, State Publisher of Literature on Building, Architecture and Building Materials, Moscow, pp. 192200 (in Russian). Moshnin, L.F., Timofeeva, E.T., 1965. &ob FIvehne sabjehnr fpn lnspabjn#ecknx ysapax, cofpobo|saiwnxcr paDpFIbamn cfjovhoctn fotoka. (Pressure-rise in water-hammer accompanied with column-separation.). Bosochad|ehne n Cahntaphar Texhnka (Water Supply and Sanitary Technology) 7, 35 (in Russian). Mostowsky, A.F., 1929. Nccjesobahne lnspabjn#eckoto ysapa b tpydax fpn majFIx hafopax. (Investigation of water hammer in pipes at low pressures.). Bulletins of the Moscow Transport Engineering Institutein memory of F.E. Dzerjinsky 11, 263304 (in Russian). Nakagawa, M., Takenaka, T., 1995. Damping and peak period shortening of a water hammer in collapse of a cooling cavity. JSME International Journal, Series B 38 (2), 300307. Nonoshita, T., Matsumoto, Y., Ohashi, H., Kubota, T., 1991. Model analysis of water column separation accompanied with swirl ow. In: Proceedings of the International Meeting on Hydraulic Transients with Water Column Separation. 9th Round Table of the IAHR Group, Valencia, Spain, pp. 235249. Nonoshita, T., Matsumoto, Y., Ohashi, H., Kubota, T., 1992. Water column separation in an elbow draft tube. In: Proceedings of the 16th IAHR Symposium on Hydraulic Machinery and Cavitation, Sao Paolo, Brazil, vol. 1, pp. 141150. Nonoshita, T., Matsumoto, Y., Ohashi, H., Kubota, T., 1999. Water column separation in a straight draft tube. In: Proceedings of the Third ASME-JSME Joint Fluids Engineering Conference, Symposium S-290 Water Hammer, San Francisco, USA, ASME-FED 248, Paper FEDSM99-6902, 9pp. OBrien, E., 1956. Discussion of Water column separation in pump discharge lines by R.T. Richards. Transactions, ASME 78, 1305. Oldenziel, D.M., Teijema, J., 1976. Cavitation on valves in correlation to liquid properties. In: Proceedings of the IAHR Symposium on Problems of Hydraulic MachineHydraulic Structure Interaction, Leningrad, USSR, pp. 221240. Also: Delft Hydraulics Laboratory, Publication No. 173. ONeill, I.C., 1959. Water-hammer in simple pipe systems. M.Sc. Thesis, University of Melbourne, Melbourne, Australia. Overton, G.D.N., Williams, P.R., Trevena, D.H., 1984. The inuence of cavitation history and entrained gas on liquid tensile strength. Journal of Physics D 17, 979987. Paredes, R., Huerta, A., Sanchez, J.L., 1987. Experimental relation between the highest transient pressure and the severity of water column separation. In: Proceedings of the Eighth International Round Table on Hydraulic Transients with Water Column Separation. IAHR, Madeira, Portugal. Parmakian, J., 1955. Waterhammer Analysis. Prentice-Hall, New York (Reprint in 1963, Dover Publications, New York). Parmakian, J., 1985. Water column separation in power and pumping plants. Hydro Review 4 (2), 8589. Pearsall, I.S., 1965. The velocity of water hammer waves. Proceedings of the Institution of Mechanical Engineers 180 (3E), 1220. Piga, E., Sambiago, G., 1974. Contributo sperimentale allo studio dei transitori in presenza di cavitazione negli impianti di raffreddamento delle centrali termoelettriche. (Contribution to the experimental study of transients in the presence of cavitation in cooling systems of thermal power stations.). In: Proceedings of the Second International Round Table on Hydraulic Transients with Water Column Separation. IAHR, Vallombrosa, Italy (in Italian). Plesset, M.S., 1969. The tensile strength of liquids. In: Proceedings of ASME Fluids Engineering and Applied Mechanics Conference. Cavitation State of Knowledge, Evanston, Illinois, USA, pp. 1525. Prasser, H.-M., Bottger, A., Zschau, J., 1998. A new electrode-mesh tomograph for gasliquid ows. Flow Measurement and Instrumentation 9, 111119. Prasser, H.-M., Scholz, D., Zippe, C., 2001. Bubble size measurement using wire-mesh sensors. Flow Measurement and Instrumentation 12, 299312.
ARTICLE IN PRESS
168 A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171 Proceedings, 1933. Joint ASCE-ASME Symposium on Water Hammer, Chicago, USA (reprinted in 1949 and 1961). Proceedings, 1937. Joint ASCE-ASME-AWWA Second Symposium on Water Hammer, New York, USA. Transactions of ASME 59, 651713, Transactions, ASME 60, 605610, 675682. Provoost, G.A., 1976. Investigation into cavitation in a prototype pipeline caused by water hammer. In: Proceedings of the Second International Conference on Pressure Surges. BHRA, London, UK, pp. 1329. Also: Delft Hydraulics Laboratory, Publication No. 170. Provoost, G.A., Wylie, E.B., 1981. Discrete gas model to represent distributed free gas in liquids. In: Proceedings of the Fifth International Symposium on Water Column Separation. IAHR, Obernach, Germany, 8pp. Also: Delft Hydraulics Laboratory, Publication No. 263, 1982. Raiteri, E., Siccardi, F., 1970. Cavitation overpressure in pressing pipelines. In: Proceedings of the XII Congresso di Idraulica e Costruzioni Idrauliche, Bari, Italy. Raiteri, E., Siccardi, F., 1975. Transients in conduits conveying a two-phase bubbly ow: experimental measurements of celerity. LEnergia Elettrica 52, 256261. Richards, R.T., 1956. Water-column separation in pump discharge lines. Transactions of ASME 78, 12971306. Rouse, H., Ince, S., 1957, 1963. History of Hydraulics. Iowa Institute of Hydraulic Research, University of Iowa, Ames, Iowa, USA. (Reprint in 1963, Dover Publications, New York.) Ruus, E., Karney, B., El-Fitiany, F.A., 1984. Charts for water hammer in low head pump discharge lines resulting from water column separation and check valve closure. Canadian Journal of Civil Engineering 11, 717742. Safwat, H.H., 1972a. Photographic study of water column separation. ASCE Journal of the Hydraulics Division 98 (HY4), 739746. Safwat, H.H., 1972b. Transients in cooling water systems of thermal power plants. Ph.D. Thesis, Delft University of Technology, Dept. of Mechanical Engineering, Delft, The Netherlands. Also: Delft Hydraulics Laboratory, Publication No. 101. Safwat, H.H., 1972c. Experimental study of pressure surges in condenser cooling water systems. In: Proceedings of the First International Conference on Pressure Surges. BHRA, Canterbury, UK, pp. 1732. Safwat, H.H., De Kluyver, J.P., 1972. Digital computations for water hammer-column separation. In: Proceedings of the First International Conference on Pressure Surges. BHRA, Canterbury, UK, pp. 5368. Safwat, H.H., Van Den Polder, J., 1973. Experimental and analytic data correlation study of water column separation. ASME Journal of Fluids Engineering 95, 9197. Safwat, H.H., Arastu, A.H., Husaini, S.H., 1986. Generalized applications of the method of characteristics for the analysis of hydraulic transients involving empty sections. In: Proceedings of the Fifth Conference on Pressure Surges. BHRA, Hanover, Germany, pp. 157167. Sanada, K., Kitagawa, A., Takenaka, T., 1990. A study on analytical methods by classication of column separations in a water pipeline. Transactions of the JSME, Series B 56 (523), 585593 (in Japanese). Sayir, M., Hausler, K., 1991. Cavitation and pipeline failure: a case study. ASME-PVP 219, Transient thermal-hydraulics and coupled vessel and piping system responses, pp. 4750. Schmid, J., 1959. Kinematographische Untersuchung der Einzelblasen-Kavitation. (Cinematographic investigation of single-bubble cavitation.). Acustica 9, 321326 (in German). Schnyder, O., 1932. Uber DruckstoXe in Rohrleitungen. (On water hammer in pipe lines.). Wasserkraft und Wasserwirtschaft 27 (5), 4954; 27(6), 6470; 27(8), 96 (in German). Schwirian, R.E., 1982. Methods for simulating uidstructure interaction and cavitation with existing nite element formulations. ASME-PVP 64, Fluid transients and uidstructure interaction, pp. 261285. Schwirian, R.E., 1984. On the use of structural nite element computer codes to perform uid-dynamical and uid-structure interactive analyses. In: Proceedings of the International Conference on Numerical Methods for Transient and Coupled Problems, Venice, Italy, pp. 186197. Sharp, B.B., 1960. Cavity formation in simple pipes due to rupture of the water column. Nature 185 (4709), 302303. Sharp, B.B., 1965a. The growth and collapse of cavities produced by a rarefaction wave with particular reference to rupture of water column. Ph.D. Thesis, The University of Melbourne, Melbourne, Australia. Sharp, B.B., 1965b. Rupture of the water column. In: Proceedings of the Second Australasian Conference on Hydraulics and Fluid Mechanics. Auckland, New Zealand, pp. A169A176. Sharp, B.B., 1977. A simple model for water column rupture. In: Proceedings of the 17th IAHR Congress, vol. 5. Baden-Baden, Germany, pp. 155161. Shinada, M., 1994. Inuence of gas diffusion on uid transient phenomena associated with column separation generated during decompression operation. JSME International Journal, Series B 37 (3), 457466. Shinada, M., Kojima, E., 1987. Fluid transient phenomena associated with column separation in the return line of a hydraulic machine press. JSME International Journal 30 (268), 15771586. Shinada, M., Kojima, E., 1989. Fluid transient phenomena associated with column separation in the prell line of a hydraulic machine press. JSME International Journal, Series II 32 (4), 550558. Shinada, M., Kojima, E., 1995. Fluid transient phenomena associated with column separation in uid power pipelines. Transactions of the JSME, Series B 61 (587), 25792585 (in Japanese). Shu, J.-J., 2003a. A nite element model and electronic analogue of pipeline pressure transients with frequency-dependent friction. ASME Journal of Fluids Engineering 125, 194199. Shu, J.-J., 2003b. Modelling vaporous cavitation on uid transients. International Journal of Pressure Vessels and Piping 80, 187195.
ARTICLE IN PRESS
A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171 169 Shuy, E.B., Apelt, C.J., 1983. Friction effects in unsteady pipe ows. In: Proceedings of the Fourth International Conference on Pressure Surges. BHRA, Bath, England, pp. 147164. Siccardi, F., 1979. Centrali termoelettriche: transitori idraulici nei circuiti di raffreddamento. (Thermal power stations: hydraulic transients in cooling circuits.). In: Proceedings of the Eighth International Round Table on Hydraulic Transients with Water Column Separation. IAHR, Cagliari, Italy (in Italian). Siemons, J., 1966. Cavitatieverschijnselen in een buisleiding. (Cavitation phenomena in a conduit.). Delft Hydraulics Laboratory, Report S 103-I, Delft, The Netherlands (in Dutch). Siemons, J., 1967. The phenomenon of cavitation in a horizontal pipe-line due to a sudden pump-failure. IAHR Journal of Hydraulic Research 5, 135152 (Also: Delft Hydraulics Laboratory, Publication No. 53). Simin, O., 1904. Water hammer. In: Proceedings of the 24th Annual Convention of the American Water Works Association, St. Louis, USA, pp. 341424. (Including a partial translation of Joukowsky, 1900). Simpson, A.R., 1986. Large water hammer pressures due to column separation in sloping pipes. Ph.D. Thesis, The University of Michigan, Dept. of Civil Engineering, Ann Arbor, USA. Simpson, A.R., Bergant, A., 1991. Column separation research at the University of Adelaide, South Australia. In: Proceedings of the International Meeting on Hydraulic Transients with Water Column Separation. 9th Round Table of the IAHR Group, Valencia, Spain, pp. 253269. Simpson, A.R., Bergant, A., 1994a. Numerical comparison of pipe column-separation models. ASCE Journal of Hydraulic Engineering 120, 361377. Simpson, A.R., Bergant, A., 1994b. Developments in pipeline column separation experimentation. IAHR Journal of Hydraulic Research 32, 183194. Simpson, A.R., Bergant, A., 1996. Interesting lessons from column separation experiments. In: Boldy, A. (Ed.), Pressure Surges and Fluid Transients in Pipelines and Open Channels. Mechanical Engineering Publications Ltd., Burry St. Edmunds, pp. 8397. Simpson, A.R., Wylie, E.B., 1985. Problems encountered in modeling vapor column separation. In: Proceedings of Symposium on Fluid Transients in FluidStructure Interaction. ASME Winter Annual Meeting, Miami Beach, Florida, USA, pp. 103107. Simpson, A.R., Wylie, E.B., 1989. Towards an improved understanding of waterhammer column separation in pipelines. Civil Engineering Transactions 1989, The Institution of Engineers, Australia, CE 31 (3), 113120. Smirnov, D.N., 1954. O lnspabjn#eckom ysape ha hacochFI x yctahobkax hnDkolo sabjehnr. (On water hammer in pump systems at low pressure.). !nspotexhnka n Mejnopaunr (Hydraulics and Irrigation) 6 (9), 4855 (in Russian). Smirnov, D.N., Zubov, L.B., 1972. Experimental investigation of water hammer accompanied with cavitation. Tr. Lab. Inzh. Gidrav., USSR Scientic Research Institute (WNII) Wodgeo, Moscow, Coll. 13 (in Russian). Stepanoff, A.J., 1949. Elements of graphical solution of water-hammer problems in centrifugal-pump systems. Transactions of ASME 83, 19191920. Streeter, V.L., 1963. Water hammer analysis with non-linear frictional resistance. In: Proceedings of the First Australian Conference on Hydraulics and Fluid Mechanics. Pergamon Press, New York. Streeter, V.L., 1964. Water hammer analysis of pipelines. ASCE Journal of the Hydraulics Division 90 (HY4), 151172. Streeter, V.L., 1965. Computer solution of surge problems. Proceedings of the Institution of Mechanical Engineers 80 (3E), 6282. Streeter, V.L., 1969. Water hammer analysis. ASCE Journal of the Hydraulics Division 95 (HY6), 19591972. Streeter, V.L., 1983. Transient cavitating pipe ow. ASCE Journal of Hydraulic Engineering 109 (HY11), 14081423. Streeter, V.L., Lai, C., 1962. Water hammer analysis including friction. ASCE Journal of the Hydraulics Division 88 (HY3), 79112. Streeter, V.L., Wylie, E.B., 1967. Hydraulic Transients. McGraw-Hill, New York. Stromeyer, C.E., 1901. On explosions of steam pipes due to water-hammers. Memoirs and Proceedings of the Manchester Literary and Philosophical Society 46 (3), 116. Strowger, E.B., Kerr, S.L., 1926. Speed changes of hydraulic turbines for sudden changes of load. Transactions, ASME 48, 209262. Su, C.-K., Camara, C., Kappus, B., Putterman, S.J., 2003. Cavitation luminescence in a water hammer: Upscaling sonoluminescence. Physics of Fluids 15, 14571461. Suda, M., 1990. Berechnung von DruckstoXen in einer ussigkeitsgefullten Rohrleitung bei Kavitationsprozessen. (Computersimula tion of water hammer and cavitation process in a liquid pipeline.). Rohre Rohrleitungsbau Rohrleitungstransport 3R (international) 29, 378381 (in German). Sundquist, M.J., 1977. Pipeline transients in highly saturated uids with gaseous and vaporous cavitation. M.Sc. Thesis, Michigan State University, East Lansing, Michigan, USA. Swafeld, J.A., 19691970. A study of column separation following valve closure in a pipeline carrying aviation kerosene. Proceedings of the Institution of Mechanical Engineers 184 (3G), 5764. Swafeld, J.A., 1972a. A study of the inuence of air release on column separation in an aviation kerosine pipeline. Proceedings of the Institution of Mechanical Engineers 186 (56/72), 693703. Swafeld, J.A., 1972b. Column separation in an aircraft fuel system. In: Proceedings of the First International Conference on Pressure Surges. BHRA, Canterbury, UK, Paper C2, 16pp. Tabei, K., Mashiko, S., Shirai, H., 2003. Study of cavitation light emission generated by a waterhammer. In: Proceedings of the Fourth ASME-JSME Joint Fluids Engineering Conference, Honolulu, Hawaii, USA, Paper FEDSM2003-45274, 6pp. Takenaka, T., 1987. Some problems on uid transient phenomena. JSME International Journal, Series B 32 (266), 12001206. Tanahashi, T., Kasahara, E., 1969. Analysis of the waterhammer with water column separation. Bulletin of the JSME 12 (50), 206214.
ARTICLE IN PRESS
170 A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171 Tanahashi, T., Kasahara, E., 1970. Comparisons between experimental and theoretical results of the waterhammer with column separation. Bulletin of the JSME 13 (61), 914925. Tarasevich, V.V., 1975. Metos kohtpojFhFIx to#ek sjr pac#eta lnspabjn#eckolo ysapa, cofpobo|saiwelocr kabntaune{ fotoka |nskoctn. (Method of control points for the calculation of water hammer with cavitating uid ow.). !nspabjnka Coopy|ehn{ OpocntejFhFIx Cnctem (Hydraulic Design of Irrigation Systems) 17 (5), 4760 (NIMI, Institute of Irrigation Engineering, Novocherkassk, Russia (in Russian)). Tarasevich, V.V., 1980. Maximum pressure during water hammer accompanied by discontinuity of the ow. Hydrotechnical Construction 14, 784790. Tarasevich, V.V., 1997. Pacfpoctpahehne n lpahcfopmaunr boDmywehn{ fpn tnspabjn#eckom ysape, cofpobo|saiwemcr kabntaune{. (Propagation and transformation of perturbations of water hammer with cavitation.). In: Kerdinski, V.K. (Ed.), Akyctnka Heoshoposhnx Cpes (Acoustics of Non-Homogeneous Media). Institute of Hydrodynamics, Novosibirsk, Russia, pp. 226234 (in Russian). Thibessard, G., 1961. La simulation du coup de belier sur calculateur numerique. (The simulation of water hammer on a numerical calculator.). In: Proceedings of the Ninth IAHR Convention, Dubrovnik, Yugoslavia, pp. 881890 (in French) Thorley, A.R.D., 1976. A Survey of Investigations into Pressure Surge Phenomena. Research Memorandum ML83. The City University, Department of Mechanical Engineering, London, UK. Thorley, A.R.D., Chohan, R.K., 1976. Development of an experimental facility for studying transient pressure waves in airwater mixtures. Research Memorandum ML 92. The City University, Department of Mechanical Engineering, London, UK. Thorley, A.R.D., Wiggert, D.C., 1985. The effect of virtual mass on the basic equations for unsteady one-dimensional heterogeneous ows. International Journal of Multiphase Flow 11, 149160. Tijsseling, A.S., 1993. Fluidstructure interaction in case of waterhammer with cavitation. Ph.D. Thesis, Delft University of Technology, Faculty of Civil Engineering, Communications on Hydraulic and Geotechnical Engineering, Report No. 93-6, ISSN: 0169-6548, Delft, The Netherlands. Tijsseling, A.S., 1996. Fluidstructure interaction in liquid-lled pipe systems: a review. Journal of Fluids and Structures 10, 109146. Tijsseling, A.S., Anderson, A., 2004. A precursor in waterhammer analysisrediscovering Johannes von Kries. In: Proceedings of the Ninth International Conference on Pressure Surges. BHR Group, Chester, UK, pp. 739751. Tijsseling, A.S., Fan, D., 1991a. The response of liquid-lled pipes to vapor cavity collapse. In: Transactions of SMiRT11, Tokyo, Japan, Paper J10/2, pp. 183188. Tijsseling, A.S., Fan, D., 1991b. The concentrated cavity model validated by experiments in a closed tube. In: Proceedings of the International Meeting on Hydraulic Transients with Water Column Separation. 9th Round Table of the IAHR Group, Valencia, Spain, pp. 145155. Tijsseling, A.S., Fan, D., 1992. Fluidstructure interaction and column separation in a closed pipe. In: Proceedings of the Second National Mechanics Congress, Kerkrade, The Netherlands. Kluwer Academic Publishers, Dordrecht, pp. 205212. Tijsseling, A.S., Vardy, A.E., 2006. Fluidstructure interaction and transient cavitation tests in a T-piece pipe. Journal of Fluids and Structures 20, 753762. Tijsseling, A.S., Vardy, A.E., Fan, D., 1996. Fluidstructure interaction and cavitation in a single-elbow pipe system. Journal of Fluids and Structures 10, 395420. Trevena, D.H., 1984. Cavitation and the generation of tension in liquids. Journal of Physics D 17 (11), 21392164. Trevena, D.H., 1987. Cavitation and Tension in Liquids. Adam Hilger, Bristol, UK and Philadelphia, USA. Trikha, A.K., 1975. An efcient method for simulating frequency-dependent friction in transient liquid ow. ASME Journal of Fluids Engineering 97, 97105. Tullis, J.P., Streeter, V.L., Wylie, E.B., 1976. Waterhammer analysis with air release. In: Proceedings of the Second International Conference on Pressure Surges. BHRA, London, UK, pp. 3547. Van De Riet, R.P., 1964. A computational method for the water hammer problem. Mathematisch Centrum, Report TW 95, Amsterdam, The Netherlands. Van De Sande, E., Belde, A.P., 1981. The collapse of vaporous cavities. An experimental study. In: Proceedings of the Fifth International Round Table on Hydraulic Transients with Water Column Separation. IAHR, Obernach, Germany. Vardy, A.E., 1980. Unsteady ows: fact and friction. In: Proceedings of the Third International Conference on Pressure Surges. BHRA, Canterbury, England, pp. 1526. Vardy, A.E., Brown, J.M.B., 2000. Friction-dependent wavefront evolution. In: Proceedings of the Eighth International Conference on Pressure Surges. BHR Group, The Hague, The Netherlands, pp. 337349. Vardy, A.E., Brown, J.M.B., 2003. Transient turbulent friction in smooth pipe ows. Journal of Sound and Vibration 259 (5), 10111036. Vardy, A.E., Brown, J.M.B., 2004. Transient turbulent friction in fully rough pipe ows. Journal of Sound and Vibration 270 (12), 233257. Vliegenthart, A.C., 1970. The Shuman ltering operator and the numerical computation of shock waves. Journal of Engineering Mathematics 4, 341348. Von Kries, J., 1883. Ueber die Beziehungen zwischen Druck und Geschwindigkeit, welche bei der Wellenbewegung in elastischen Schlauchen bestehen. (On the relations between pressure and velocity, which exist in the wavelike motion in elastic tubes.). Festschrift der 56. Versammlung deutscher Naturforscher und Arzte gewidmet von der Naturforschenden Gesellschaft zu Freiburg i. B., Freiburg, Germany, pp. 6788 (in German). (see: Tijsseling and Anderson, 2004.)
ARTICLE IN PRESS
A. Bergant et al. / Journal of Fluids and Structures 22 (2006) 135171 171 Vreugdenhil, C.B., 1964. Digital computations of water-hammer. Delft Hydraulics Laboratory, Report S 103-I, Delft, The Netherlands. Vreugdenhil, C.B., De Vries, A.H., Kalkwijk, J.P.Th., Kranenburg, C., 1972. Investigation into cavitation in long horizontal pipelines caused by water hammer. In: Transactions of the sixth IAHR Symposium, Section for Hydraulic Machinery, Equipment and Cavitation, Rome, Italy, Paper J3. Also: Kalkwijk et al., Delft Hydraulics Laboratory, Publication No. 115, 1974. Wallis, G.B., 1969. One-dimensional Two-phase Flow. McGraw-Hill, New York. Walsh, J.P., 1964. Pressure generated by cavitation in a pipe. M.Sc. Thesis, Syracuse University, New York, USA. Wang, J.S., Locher, F.A., 1991. Verication of modelling water column separation. In: Proceedings of the International Meeting on Hydraulic Transients with Water Column Separation. 9th Round Table of the IAHR Group, Valencia, Spain, pp. 343354. Wangjiraniran, W., Motegi, Y., Richter, S., Kikura, H., Aritomi, M., Yamamoto, K., 2003. Intrusive effect of wire mesh tomography on gasliquid ow measurement. Journal of Nuclear Science and Technology 40, 932940. Washio, S., Takahashi, S., Konishi, T., Morikawe, H., 1994. Creation and observation of tensile waves in oil column. JSME International Journal, Series B 37 (2), 342348. Washio, S., Takahashi, S., Yamaguchi, S., 1996. Measurement of transiently changing ow rates in oil hydraulic column separation. JSME International Journal, Series B 39 (1), 5156. Watt, C.S., Boldy, A.P., Hobbs, J.M., 1980. Combination of nite difference and nite element techniques in hydraulic transient problems. In: Proceedings of the Third International Conference on Pressure Surges. BHRA, Canterbury, UK, pp. 4362. Weyler, M.E., 1969. An investigation of the effect of cavitation bubbles on momentum loss in transient pipe ow. Ph.D. Thesis, The University of Michigan, Ann Arbor, USA. Whiteman, K.J., Pearsall, I.S., 1962. Reux valve and surge tests at a pumping station. Fluid Handling 13, 248250, 282286. Wiggert, D.C., Sundquist, M.J., 1979. The effect of gaseous cavitation on uid transients. ASME Journal of Fluids Engineering 101, 7986. Wiggert, D.C., Tijsseling, A.S., 2001. Fluid transients and uidstructure interaction in exible liquid-lled piping. ASME Applied Mechanics Reviews 54, 455481. Williams, P.R., Williams, R.L., 2000. On anomalously low values of the tensile strength of water. Proceedings of the Royal Society (London) A: Mathematical and Physical Sciences 456, 13211332. Williams, P.R., Williams, P.M., Brown, S.W.J., Temperley, H.N.V., 1999. On the tensile strength of water under pulsed dynamic stressing. Proceedings of the Royal Society: Mathematical, Physical and Engineering Sciences A 455, 33113323. Wylie, E.B., 1980. Free air in liquid transient ow. In: Proceedings of the Third International Conference on Pressure Surges. BHRA, Cantenbury, UK, pp. 2742. Wylie, E.B., 1984. Simulation of vaporous and gaseous cavitation. ASME Journal of Fluids Engineering 106, 307311. Wylie, E.B., 1992. Low void fraction two-component two-phase transient ow. In: Bettess, R., Watts, J. (Eds.), Unsteady Flow and Fluid Transients. A.A. Balkema, Rotterdam, pp. 39. Wylie, E.B., 1999. Simple device to demonstrate vaporization. In: Proceedings of the Third ASME-JSME Joint Fluids Engineering Conference, Symposium S-290 Water Hammer, San Francisco, USA, ASME-FED 248, Paper FEDSM99-6876, 7pp. Wylie, E.B., Streeter, V.L., 1978a. Fluid Transients. McGraw-Hill, New York (Republished with minor corrections by FEB Press, Ann Arbor, Michigan, USA, 1983.). Wylie, E.B., Streeter, V.L., 1978b. Column separation in horizontal pipelines. In: Proceedings of the Joint Symposium on Design and Operation of Fluid Machinery, vol. 1. IAHR/ASME/ASCE, Colorado State University, Fort Collins, USA, pp. 313. Wylie, E.B., Streeter, V.L., 1993. Fluid Transients in Systems. Prentice-Hall, Englewood Cliffs. Yamaguchi, K., Ichikawa, T., 1976. Pressure surge due to oil column separation. In: Proceedings of the Second International Conference on Pressure Surges. BHRA, London, UK, pp. 2946. Yamaguchi, K., Ichikawa, T., Suzuki, S., 1977. Transient characteristics of oil pipelines with column separation. Bulletin of the JSME 20 (143), 630637. Yang, W.Q., Adam, M.S., Watson, R., Beck, M.S., 1996. Monitoring water hammer by capacitance tomography. Electronic Letters 32 (19), 17781779. Young, F.R., 1989. Cavitation. McGraw-Hill, London. Yow, W., Van Duyne, D.A., Budlong, L.A., Quinlan, P.J., 1985. Column separation and rejoining eld study. In: Proceedings of Forum on Unsteady Flow. ASME Winter Annual Meeting, Miami, USA. ASME-FED 27, pp. 1820. Zielke, W., 1968. Frequency-dependent friction in transient pipe ow. ASME Journal of Basic Engineering 90, 109115. Zielke, W., Perko, H.-D., 1985. Unterdruckerscheinungen und DruckstoXberechnung (Low pressure phenomena and water hammer analysis.). Rohre Rohrleitungsbau Rohrleitungstransport 3R (international) 24, 348355 (in German). Zielke, W., Perko, H.-D., Keller, A., 1989. Gas release in transient pipe ow. In: Proceedings of the Sixth International Conference on Pressure Surges. BHRA, Cambridge, UK, pp. 313.