Experimental and Numerical Investigation on the

0 downloads 0 Views 4MB Size Report
Results indicated that when flow drag and buoyancy forces are summed the observed. 13 ... Martin (1976), which was complemented by Zhou et al (2002) among others, indicated that compressed air. 43 .... using lumped inertia analysis and the ideal gas law. ... which was computed using the isothermal Euler equations.
1

Experimental and Numerical Investigation on the Motion of Discrete Air

2

Pockets in Pressurized Water Flows

3

Carmen D. Chosie1 Thomas M. Hatcher2 Jose G. Vasconcelos3 , A.M. ASCE

4

ABSTRACT

5

Entrapped air pockets are linked to various operational issues in closed pipe flows. Tracking of entrapped

6

air pockets is a relevant task in numerical simulation of hydraulic systems subjected to air pocket entrapment.

7

However, few studies focused on the kinematics of discrete air pockets and the motion of these at various

8

combinations of pipe flow velocities, air pocket volumes and pipe slopes. Such information is relevant for

9

the improvement of current numerical models for the simulation of transition between pressurized and free

10

surface flows, and related air pocket entrapment. This work explores the link between pipe flow velocity, pipe

11

slope, and entrapped air pocket celerities. Systematic experimental studies were performed and provided

12

insights on the nature of entrapped air pocket motion in cases when buoyancy forces are either summed or

13

opposed to flow drag. Results indicated that when flow drag and buoyancy forces are summed the observed

14

air pocket celerity can be approximated by the celerity in quiescent conditions plus the pipe velocity. This

15

approximation fails in cases when drag force opposed buoyancy.This work also proposes a novel approach

16

for the simulation of the air pocket motion as non-Boussinesq gravity currents. Modeling results compared

17

favorably with experimental measurements in horizontal slope and pipe flow velocity cases. Such results

18

are an indication of the feasibility of integrating this approach into more complex hydraulic models able to

19

simulate flow regime transitions.

20

Keywords: Closed pipe flows; Air pockets; Laboratory experiments; Numerical modeling

21

INTRODUCTION

22

A key problem for closed conduit flows in which two-phase flows may occur is the limited information

23

on the kinematics of entrapped air pockets. The ability of describing the motion of entrapped pockets is

24

important as it can indicate the likelihood of uncontrolled air releases through ventilation shafts in stormwater

25

systems. Such information could also help assess the effectiveness of air pocket removal in transmission 1 Graduate student, Dept. of Civil Engrg., Auburn University, 238 Harbert Engineering Center, Auburn, AL, 36849. E-mail: [email protected] 2 Ph.D. candidate, Dept. of Civil Engrg., Auburn University, 238 Harbert Engineering Center, Auburn, AL, 36849. E-mail: [email protected] 3 A.M. ASCE, Assistant Professor, Dept. of Civil Engrg., Auburn University, 238 Harbert Engineering Center, Auburn, AL, 36849. E-mail: [email protected]

1

26

mains. Finally, the ability to track entrapped air pockets is relevant for the simulation of flow regime

27

transition between free surface and pressurized flows, with the goal of detecting a priori deleterious air-

28

water interactions. Past investigations dealing with topics related to the motion of air pockets in closed

29

conduit flows yielded important insights in this process. Yet, as is discussed, some knowledge gaps remained

30

both in terms of experimental and numerical approaches on how to predict this motion.

31

Most water conveyance systems based on closed conduits are not designed to operate in two-phase flows,

32

yet such undesirable flows have been observed in various transmission pipelines as well in in stormwater sys-

33

tems. Various mechanisms for large air pocket entrapment have been identified to date, and a comprehensive

34

summary is presented in Lauchlan et al (2005) in the context of transmission mains. Other mechanisms have

35

been identified in the context of the rapid filling of stormwater systems. The pioneer study by Hamam

36

and McCorquodale (1982) indicated that interface instability caused by the relative motion between air and

37

water can lead to surface waves in water that can grow and eventually promote air pocket entrapment. Other

38

mechanisms for air pocket entrapment in the context of stormwater systems have been studied experimentally

39

by Vasconcelos and Wright (2006).

40

The presence of air pockets in closed conduit water flows is generally undesirable because it can cause

41

various operational problems. One of the clearest adverse impacts is the increased potential for structural

42

damage, caused by the great compressibility of air that leads to increased surging. The pioneer study by

43

Martin (1976), which was complemented by Zhou et al (2002) among others, indicated that compressed air

44

pressure can greatly exceed the pressure that caused the flow prior to air compression. Zhou et al (2002)

45

points out that such compression may be behind an episode in which an entire manhole structure was blown

46

off the pipe. Issues related to potential structural damage in rapid filling of closed conduits have also been

47

explored by Song et al (1983) and Guo and Song (1991) in the context of pressurization of large stormwater

48

tunnels. A more recent work by Vasconcelos and Leite (2012) confirmed these earlier findings, but also

49

showed that when relief is provided next to the location where air is compressed, as would be anticipated in

50

some real applications, the pressure rise can be significantly mitigated.

51

Another adverse impact related to air pocket entrapment is loss of conveyance and storage capacity

52

in closed conduits. As shown in Falvey (1980), entrapped air pockets act as flow constrictions generating

53

additional turbulence and hence energy losses. Also, large entrapped air pockets will occupy a volume in

54

closed conduits that would otherwise be filled with water. This could be a relevant issue in the context of

55

stormwater storage tunnels since early pressurization of these systems caused by large entrapped air pockets

56

have the potential of increasing frequency of combined sewer overflows (CSO) episodes.

57

Release of entrapped air pockets can be a problematic issue as well. In the context of water transmission

58

mains one concern is the possibility of air slamming, defined as the waterhammer-type pressures observed

2

59

at the moment when air is completely evacuated through air valves. The study by Lingireddy et al (2004)

60

presents a formulation to calculate such peaks in terms of the diameter of the pipeline, the diameter of the

61

air valve and the air pocket pressure. In context of large collection sewers and tunnels, release of large air

62

pockets through water-filled ventilation shafts have been linked to geysering phenomena. This phenomenon

63

was observed in physical modeling studies by Vasconcelos and Wright (2011), as well as field observations of

64

large trunk sewers Wright et al (2011) and in penstocks Nielsen and Davis (2009). An interesting aspect in

65

Wright et al (2011) studies is that while water may be ejected violently above ground, the hydraulic grade

66

line was far below the ground surface elevation.

67

Another context where very important contributions in this field originates are investigations on the

68

motion of air cavities as gravity currents. In such studies water typically discharges in one end of an initially

69

filled pipe. An exchange flow condition is established at the discharge point as an air cavity is admitted

71

at the top of the pipe and moves toward the opposite end of the pipe. The work by Benjamin (1968) was √ one of the pioneers in establishing that the celerity of the air cavity leading edge scales with gH, with

72

g the gravitational acceleration and H the pipe characteristic dimension (e.g. diameter, D). It was also

73

determined that the gravity current celerity generally decreased with its thickness for a given pipe cross-

74

sectional dimension. The study from Wilkinson (1982) expanded the work by Benjamin by incorporating

75

surface tension effects, which have an effect on the location of the stagnation point located at the leading

76

edge of the air-water interface at the pipe crown. Baines (1991) work followed these studies and studied an

77

interesting air gulping flow feature observed when the pipe was set in a downward slope, in which the outlet

70

79

was at a lower elevation. This work showed that the observed celerity increased with slope (tested up to 8 √ degrees) and that the discrete air pocket celerity formed after gulping also scaled with gH. A limitation

80

to these studies was that there was no pipe flow prior to the arrival of the air cavity. A related work by

81

Vasconcelos and Wright (2008) presented the advance of an air cavity against an accelerating pressurized

82

flow column and indicated that the increasing water flow velocity could eventually reverse the advance of the

83

cavity. At this point in time, however, the authors determined that the air-water shape is no longer strongly

84

curved but was strongly affected by the brink depth discharge conditions at the far-downstream end.

78

85

The “hydraulic clearing” of transmission mains presents the problem of air pocket motion in a different

86

prism with the very practical objective of removing air from mains. Because of the potential issues in closed

87

conduits, various studies have been presented to date attempting to determine a minimum value for the

88

pipe flow velocity that would ensure that drag is sufficiently large enough to overcome anticipated air pocket

89

buoyancy forces and clear pipelines of pockets. Works of Falvey (1980), Little et al (2008), and Pozos et

90

91

al (2010) provide examples of formulas to compute the velocity magnitude, which depends on pipe angle √ and gD. While very useful, these studies have not focused on predicting the displacement and celerity of

3

92

discrete air pockets for various combinations of pipe flow velocities and pipe slopes.

93

Studies performed by Pothof and Clemens (2010) focused on the motion of elongated air pockets in

94

favorable slope pipelines. These studies were useful in providing insights as to the controlling aspects of the

95

pocket motion, with potential application in mitigating block back occurrences. A separate work presented

96

by Pothof and Clemens (2011) focused on the the co-current flow of air and water in downward sloping pipes.

97

These studies determined the water discharge rate required to prevent air accumulation as a function of air

98

flow rate, pipeline diameter, and pipeline slope. These experimental results pointed various dimensionless

99

parameters that are useful in characterizing air-water flow regimes. Such studies, however, have not focused

100

in predicting the displacement and celerity of discrete air pockets for various combinations of pipe flows

101

velocities and slopes, including cases when buoyancy and drag forces are not in opposition.

102

Glauser and Wickenhauser (2009) provided a relevant experimental study on the dynamics of an air

103

cavity advancing in a pressurized pipe. The focus of that work was on the shape and movement of single air

104

bubbles in continuous air-water flow with the intention to determine the bubble volume that represents the

105

stagnation velocity based on a balance of buoyancy and drag forces. It was noted that in favorably sloped

106

pipes, volumes larger than this critical volume move against this flow due to buoyancy and volumes smaller

107

than this critical volume would be dragged by the flow against buoyancy. Similar to studies in hydraulic

108

clearing of water mains, it was determined that pipeline slope, background water velocity, and bubble volume

109

controlled the observed pocket celerity. This investigation involved the use of slopes ranging from 1.7% to

110

8.7%, so effects such as air pocket spreading observed in horizontal slopes could not be evaluated.

111

Numerical studies have been presented that attempt to describe the motion of entrapped air pockets in

112

closed conduit flow. The majority of these studies are in the realm of multi-phase flow applications, and

113

exemplified by the works of Barnea and Taitel (1993); DeHenau and Raithby (1995) and Issa and Kempf

114

(2003). In these applications there is a continuous injection of gas and liquid in closed conduits, and the

115

mechanisms for pocket formation are diverse (e.g. terrain induced slugging) from the ones that are more

116

frequently observed in water, wastewater and stormwater systems. Two studies were to numerically describe

117

the motion of entrapped air pockets in water systems. The study by Li and McCorquodale (1999) presents

118

the mathematical framework for the simulation of flow regime transition (also referred to as mixed flows)

119

using lumped inertia analysis and the ideal gas law. This works also considered the mechanism for air

120

pocket formation based on shear flow instabilities presented in Hamam and McCorquodale (1982). In the

121

formulation the location of the air-water interface is calculated explicitly, providing a means to compute the

122

advance of an entrapped pocket. Yet, the velocity of the air bubble, rather than being calculated by the

123

model, appears as a calibration parameter in the simulations.

124

More recently Trindade and Vasconcelos (2013) presented a study that included a numerical framework

4

125

to simulate pipeline priming operations considering the possibility of air pressurization using the Two-

126

component pressure approach following the work by Vasconcelos and Wright (2009). With regards to air

127

pressurization, this work included the possibility of air pressure variation along the length of the air pocket,

128

which was computed using the isothermal Euler equations. The model updates the air-water interface using

129

a source term for the Euler equations presented in Toro (2009). The trajectory of air-water interface was

130

accurately predicted including complex flow interactions such as interface breakdown of the pressurization

131

interface. A model limitation was that the only mechanism that accounted for the motion of the air-water

132

interface was the displacement of the air phase caused by the water inflow during the priming event. Also,

133

it was assumed that one of the boundaries of the air pocket was fixed at the air release valve. An important

134

difficulty in the simulation of entrapped pocket motion is that such one-dimensional models are generally

135

constructed with the hypothesis of hydrostatic pressure distribution at the flow cross-sections, which will

136

not hold at strongly curved air-water interfaces.

137

Understanding of air-water interactions is key to a more accurate description of unsteady closed conduit

138

flows. Previous studies have provided significant insights for air-water interactions in closed conduit flows,

139

yet significant knowledge gaps remain with respect to air pocket kinematics. These gaps exist both in terms

140

of tested experimental conditions and in terms of numerical modeling approaches to describe motion of

141

entrapped air pockets. The proposed research attempts to address some of these gaps, as is outlined below

142

in the Objectives.

143

OBJECTIVES

144

This research program has two main objectives. The first is to further explore the link between pipe flow

145

velocity, pipeline slope and the celerity of entrapped air pockets of various volumes, including effects of air

146

pocket spreading. This work presents results from an experimental program on the kinematics of entrapped

147

air pockets in steady pressurized water flows and shallow pipe slopes.

148

The second objective is to propose and assess the accuracy of an innovative framework to describe

149

the motion of air pockets in horizontal slopes based on gravity current theory. In particular, this proposed

150

framework intends to provide a simpler method to model entrapped air pocket motion that can be embedded

151

in unsteady flow models in the future. Numerical model and experimental results are compared and discussed

152

along with conclusions and recommendations for future work.

153

METHODOLOGY

154

Experimental methods

155

Figure 1 presents a schematic of the apparatus used in the experimental procedure. These experiments

156

were performed with a variable slope, clear PVC pipeline supplied upstream by a fixed head reservoir

5

157

with a throttled discharge downstream to ensure pressurized flow throughout the pipeline. The pipeline

158

has a length L = 11.1 m and a diameter D = 101.6 mm. The upstream supply reservoir is attached to

159

the pipeline by a 50 mm diameter ball valve. Two knife gate valves with the same diameter as that of

160

the pipeline were positioned 3 m downstream of the upstream supply reservoir with a 2.1 m separation for

161

adverse and horizontal sloped conditions. For positive sloped experiments, the intermediate knife gate valves

162

were relocated to a position 3 m upstream of the downstream end of the pipe. Pre-determined volumes of

163

air at atmospheric pressure were injected into the region between the two knife gate valves forming an air

164

pocket. These partially closed knife gate valves prevented the motion of the air pockets while still allowing

165

flow underneath.

166

Air pockets of different volumes (V ol) were injected via a graduated acrylic tank located at the bottom

167

of the clear PVC pipeline, which is separated from the water flow by a control valve. This supply tank was

168

equipped with valves that ensured the air injected into the water filled pipeline was at atmospheric pressure.

169

170

171

The measurement devices used in the experimental procedure include: • Up to four high definition camcorders (1080 pixels, 30 FPS), providing a total view of 7.0 m upstream and downstream of the intermediate knife gate valves;

172

• Two piezoresistive pressure transducers, MEGGIT Endevco 8510B-5, located at the upstream end

173

(L1 = Lcoord /L = 0.0) and at approximately 50% of the pipeline length (L1 = 0.5), recording with a

174

frequency of 100 Hz for each channel;

175

• National Instruments data acquisition board NI-USB 6210 with SignalExpress logging software;

176

• Digital manometers located at both ends of the apparatus and at 50% of the pipeline length (L1 = 0.5);

177

178

and • Paddle-wheel flow meter to gage pipe flow rate.

179

Conditions that merited testing included those that represented a wide range of air pocket volumes,

180

pipeline slopes, and background flow rates that stormwater storage tunnels could encounter during rapid

181

filling after an intense rain event.

182

The experimental procedure was performed as follows:

183

1. Set the selected slope for the pipeline apparatus;

184

2. Create steady flow conditions in the PVC pipeline by turning on submersible pumps at desired flow

185

186

187

rate and opening the intermediate knife gate valves to desired opening percentage; 3. Inject air volume at atmospheric pressure by opening the valve connecting the graduated acrylic cylinder and the water filled PVC pipeline;

6

188

189

4. Allow system to achieve a new steady state condition due to the increase in energy losses created by the addition of air into the system;

190

5. Initiate pressure measurements and collect readings on calibration manometers;

191

6. Open the intermediate knife gate valves simultaneously and rapidly (t < 0.6 s) to allow sudden release

192

193

194

195

196

of air pocket; 7. Record motion and spreading of the air pocket(s) with camcorders until moving outside field of view; and 8. Close valves and stop submersible pumps after considerable time and make final readings on the calibration manometers.

197

The results presented in this paper include 99 differing conditions. Seven slopes have been tested:

198

horizontal slope, 0.5%, 1.0%, and 2.0% adverse slope(downstream end at higher elevation than inlet) and

199

0.5%, 1.0%, and 2.0% positive slope (downstream end at lower elevation than inlet). For each of these

200

slopes, three different flow rates ranging from 1.4 L/s to 4.0 L/s and up to four air pocket volumes were

201

tested. These pocket volumes ranged between 1.3 L and 4.3 L for horizontal slopes and between 0.5 L and

202

4.0 L for adverse and positive slopes.

203

It was noted visually that a pronounced drag created by larger flow rates caused bits of some air pockets

204

to be dragged underneath the knife gate valves intended to keep the pocket in place as a whole until the start

205

of the experiment. In order to keep this phenomenon from happening, certain runs with larger air pocket

206

volumes and larger flow rates were not able to be performed. Along with the flow rates and air pockets

207

volumes mentioned above, experiments were run to study the thinning of the air pocket when a single knife

208

gate valve was opened with no background flow.

209

Pressure monitoring throughout the system allowed the following objectives to be completed: gain insight

210

into the pressure variation in the system during motion of released air pockets, assess amount of energy loss

211

caused by the partially closed knife gate valves holding the pocket in its initial position, and measure

212

thickness of air pockets in selected cases since there is evidence in literature of the relationship between

213

pocket thickness and celerity during motion. Each case was repeated at least once in order to ensure that

214

data was consistent. Table 1 presents a summary of the parameters tested in this investigation.

215

Numerical modeling

216

According to modeling objectives, simulation of gravity current flows can be performed with three primary

217

approaches: 1) accurate CFD simulations solving Navier-Stokes equations; 2) approximations based on

218

shallow water equation models; and 3) simple integral models. As this work is not focused on air-water

219

flow details and interactions but instead aims to track the endpoint coordinates of air pockets, a simple

7

220

integral model was formulated (this approach is outlined in Ungarish (2009)). For this numerical model, a

221

relationship between the gravity current thickness and its propagation speed is utilized to update the location

222

of the gravity current fronts/leading edges. Integral models assume uniform flow depth at each time step

223

while providing the average gravity current depth and celerity. Because of this uniform depth assumption,

224

integral models are also referred to as “box models” Ungarish (2009).

225

For the experiments presented in this work, the integral model is developed for horizontal (no slope) flow

226

conditions. Figure 2 presents a schematic diagram of this modeling approach where x0 and y0 are respectively

227

the initial length and thickness (measured from the pipe crown) of the air pocket. After the knife gate valves

228

are fully opened, the air pocket spreads in both directions at unique velocities according to the leading edge

229

celerity and relative pipe flow velocity. Front conditions relate the air pocket leading edge celerity and its

230

thickness and offer a closure to the integral model along with the air phase continuity equation. The front

231

conditions used in integral model framework are adapted to account for pipe flow velocities following the

232

rationale by Hallworth et al (1998). The resulting front condition is:

uf =

p dxf = F gD + Vf low dt

(1)

233

where xf is the distance traveled by the respective air pocket front, F is the local Froude number at the

234

upstream and/or downstream air pocket front, D is the pipe diameter and Vf low is the background flow

235

velocity.

236

Each front condition uses a different formulation to determine the value of F at the leading edge of the

237

air pocket. The first front condition applies the theory by Benjamin (1968) in which there is a momentum

238

balance solved at the air pocket leading edge. The control volume (CV) encloses the leading edge so that the

239

flow is uniform and hydrostatic at the CV boundaries. This formulation considers surface tension following

240

Wilkinson (1982). After integrating Equation 1 and solving for xf , the result is the first integral model

241

considered in this work (referred to as M 1): √ uf = (F gD + Vf low )(tnew − told )   1/2 2πDA2 σ 1 2 3  πD A2 − ρgrc − 4A2 A2 cos α + 6 D sin α     F =   1 2 2 πD πD − A2 2 

2

(2)

242

where tnew is the updated time step; told is the previous time step; α = is the pipe half-angle, with 2α as

243

the angle subtended from the center of the pipe to the free surface; σ is the air-water surface tension; rc is

244

the leading edge radius of curvature; and A2 is the cross-sectional flow area beneath the air pocket equal to

8

245

π−α+

1 2

sin α

 D2 4

.

246

The second integral model approach (M 2) utilizes the front condition presented in Vasconcelos and

247

Wright (2008), which in turn was based on the work by Townson (1991). This front condition accounts

248

for the curvature at the air pocket front and has the advantage of accounting for the depth change at

249

the hydraulic jump that trails the leading edge of the gravity current. A limitation is that while this

250

front condition formulation does not assume uniform air pocket thickness, the integral model makes this

251

assumption. This has impacts in the accuracy of model predictions, as is shown in the Results section. The

252

expression for the 2nd proposed integral model (M 2) is: √ uf = (F gD + Vf low )(tnew − told )  1/2 1  3 1− 3 sin α + 3(π − α) cos α − sin α  π − α − sin α cos α  3π     F =   A2 A2 π 2− π π 

(3)

253

Since there are initially two air intrusion fronts for the cases studied in this work, Equation 2 or 3 are

254

solved twice at each time step. With the front velocities known, the locations of each front are updated

255

(xf U and xf D ) according to Figure 2. The methodology used to solve each leading edge location is identical

256

except that the sign of U changes depending on the orientation of the background flow with respect to the

257

front motion. One of the difficulties in this procedure is that xf (xf U or xf D ) is a function of α and vice

258

versa. For the implementation utilized in this work, xf was first computed from parameters determined at

259

the previous time step. Then, α was updated using the following volume/continuity expression: Lp D2 V ol = 4

  1 α − sin 2α 2

(4)

260

in which V ol and Lp = x0 + xf U + xf D are the air pocket volume and length, respectively. This nonlinear

261

expression was solved at each time step using the iterative Newton-Raphson method with a tolerance for α

262

of 10−5 . A schematic diagram for this integral model formulation is displayed in Figure 2.

263

Albeit relatively simple, the integral models formulated in this work account for buoyancy, drag, and

264

background flows, but there are important limitations. The model formulation utilizes steady-state theory

265

for the air cavity motion in which zero displacement is assumed at the stagnation point. The steady-

266

state assumption is valid during the initial slumping stage, but the air pocket motion will diverge from

267

steady-state theory in later stages of the propagation Simpson (1997). In addition, as previously stated, the

268

spatial variation in depth is neglected. Finally, the adopted radius of curvature at the front leading edge

269

(rc = 0.0034 m) followed the work of Wilkinson (1982) in rectangular cross-sections, while results in this

9

270

work were obtained in circular pipes.

271

RESULTS AND DISCUSSION

272

Measurements of air pocket spreading with no pipe flows

273

Figure 3 presents the path of the air pocket leading edge location over time and corresponding celerity

274

for four different air pocket volumes for horizontal slopes without pipe flows. Such paths are defined in

275

this work as air pocket front trajectories, and one may notice the curved nature of the trajectories. This

276

curvature indicates a reduction in the leading edge celerity, caused by the air pocket spreading and reducing

277

of the air pocket thickness. The slope of the air pocket’s leading edge trajectory decreased for larger air

278

pocket volumes, indicating an increase in celerity, an expected result. It can be noticed that the trajectory

279

of the back-propagating front is symmetrical to the forward propagating front in the absence of flow in the

280

pipe, also as anticipated. In all results, air pocket volumes are presented in a normalized format, with V ol∗

281

defined as V ol/D3 .

282

Results for adverse slopes (downstream end at a higher elevation than the upstream end) and no pipe

283

flows are presented Figure 4. Trajectories and celerity values of the leading edge of the moving air pockets

284

are shown for various volumes as well as various adverse pipeline slopes. The results do not present the

285

curved trajectory from the air pocket leading edge that was present in horizontal cases even for the smallest

286

slope tested of 0.5%. This indicates that the effect of air pocket spreading at the pocket leading edge was

287

not significant in the experiments. Another observation of this figure indicates the relative celerity difference

288

between pocket of various volumes decreases for larger pipeline slopes. That may indicate that larger slopes

289

lead to stronger concentration of air at the leading edge of the air pocket as it propagates, contributing to

290

the reduction of these celerity values. Results for the smallest air pocket volume of V ol∗ = 0.48 and 2% slope

292

show a much smaller celerity than corresponding cases with shallower slopes, a result which could not be √ explained at this point. However, for most tested cases, the relative celerity C ∗ = C/ gD is limited by 0.47,

293

which is consistent with results presented by Baines (1991). Results obtained when slopes were reversed to

294

positive and no flow are consistent with the ones with the cases with adverse slope, and are not shown here

295

for brevity.

296

Measurements of air pocket spreading with pipe flows

291

297

Results are significantly affected by the presence of pipe flows. Figure 5 presents the trajectories for

299

the cases with air pocket volumes of V ol∗ =1.2 and V ol∗ =2.2 and horizontal slope for various pipe flows. √ Celerity values C ∗ are complemented by celerity values relative to the pipe flow flow (C - Vf low / gD. The

300

symmetry between the two air pocket fronts previously observed for horizontal case is no longer observed due

301

to drag forces caused by the introduction of pipe flow. The results in Figure 5 indicate that larger flow rates

298

10

302

result in an increase in the celerity of the forward moving air pocket front, and that this gain in velocity is

303

proportional to the flow velocity. Results also indicate that in the initial stages following the pocket release

305

there is an increase in the celerity that is generally finished at the instant t∗ ≈ 2, with the normalized time p defined as t∗ = (t/X0 )/( D/g) . The previously observed air pocket front spreading in the horizontal slope

306

and no flow, which led to decreasing celerity values, is not observed is when there is pipe flows. Results in

307

Figure 5 support the assumption that for horizontal pipes the pocket celerity within flowing pipes can be

308

estimated as the summation of its value in quiescent conditions plus the pipe flow velocity.

304

309

An interesting result is that is even for the smallest pipe flow, the backward propagating air pocket leading

310

edge is eventually sheared and no longer observed for values of L1 < 0.3 and t∗ ≈ 1. This indicate that in

311

the absence of buoyancy forces the propagation of discrete air pockets against an non-zero pipe flow velocity

312

is short lived. Figure 6 presents seven snapshots of the pocket propagation with 1 second time lapse between

313

the images, for a condition corresponding to a pocket with V ol∗ =4.1 and flow rate Q∗ =0.27. Air-water

314

interfaces in the figure are artificially enhanced for clarity. The initial shape of the front resembles a typical,

315

dissipative gravity current with a curved nose with a trailing a hydraulic jump, as described in Benjamin

316

(1968). As the front moved upstream, the trailing hydraulic jump entrained air at the leading edge of the

317

pocket, and these air bubbles were carried with the pipe flow. The pocket volume thus decreased over time

318

and within 7-8 seconds later it effectively vanished. In the process, the observed pocket celerity decreased as

319

the thickness of the backward propagating front decreased, a result consistent with gravity current theory.

320

The other extremity of the front did propagate as a curved gravity current front. One speculates that this

321

type of pocket shearing may also be observed in deep stormwater storage tunnels laid in shallow slopes. This

322

behavior was not observed for downward slopes, when buoyancy forces oppose drag forces, as is explained

323

below.

325

√ Figure 7 presents the trajectories, observed celerity C normalized by gD, and the relative celerity √ accounting for pipe flow velocity (C - Vf low ) also normalized by gD, for the cases with V ol∗ =0.95 and

326

V ol∗ =1.9 and 1% adverse slope. As anticipated, the observed air pocket front celerity increases with the

327

flow rate, but there is no significant changes in celerity between the two tested pocket volumes. In such

328

slopes buoyancy and drag forces are summed and the air pocket celerity can also be approximated by the

329

summation of the pipe flow velocity and the celerity observed in quiescent conditions. An exception to this

330

observation was the case with V ol∗ =0.95 and the highest pipe flow rate Q∗ =0.37. An explanation of this

331

observation is not available at this point.

324

332

The most complex cases studied in this work involved positive slopes, when buoyancy forces opposed drag

333

forces. A sample of these results is presented in Figure 8 that corresponds to the case where pocket volumes

334

V ol∗ =0.95 and V ol∗ =1.9 and 1% positive slope. Negative celerity values are due to the propagation direction

11

335

of the pocket toward the upstream end of the apparatus. Visually, the leading edge of the pockets advanced

336

upstream against the pipe flow due to buoyancy in all tested slopes. Strong turbulence was observed at

337

the trailing hydraulic jump behind the air pocket front. Air was sheared from the air pocket leading edge,

338

but consistently with observations are consistent with those by Pozos (2007) these smaller bubbles gathered

339

downstream, and eventually grew in size and also moved against the flow and rejoin with the main pocket.

340

The observed celerity C ∗ values for these positive slopes, as anticipated, are smaller than the case with

341

no pipe flow and ranged from -0.2 to -0.4 for all tested conditions. However, the assumption that these

342

celerity values corresponded to the summation of celerity in quiescent conditions and the pipe flow velocity

343

conditions is no longer valid, a result that is likely linked to these complex air-water interactions observed

344

in positive slopes.

345

Comparison with numerical model predictions

346

The proposed gravity current integral model was tested for four different air pocket volumes and all

347

background flow rates, all involving horizontal pipe. The limitation of slope was due to the numerical model

348

inability to incorporate slope effects. Two alternatives (M 1 and M 2) discussed in the Methodology were

349

compared to these experimental results to better understand the role of front conditions in the predictions of

350

air cavity motion. Both alternatives were simple to implement, involved relatively inexpensive calculations,

351

and thus simulations were concluded in a fraction of a second. Another limitation of the numerical modeling

352

is that comparisons only involved downstream moving fronts when pipe flows were present. This is mostly

353

due to the fact that the upstream moving fronts lasted briefly as they were quickly sheared by pipe flows.

354

Such conditions yielded a limited amount of experimental data points and prevented a comparison between

355

experimental and model results.

356

The air cavity front trajectories for horizontal flow experiments are compared to both numerical model

358

alternatives (M1 and M2) in Figure 9. The pocket front coordinate (X) was normalized by the initial √ air pocket length (X0 ) and the time by (X0 / gD). The experimental results are clearly affected by the

359

background flow rate in which the air cavity front velocity is governed by Equation 1 as described in Hallworth

360

et al (1998). The effect of air pocket volume on front velocity in the tested conditions was not as important

361

as the pipe flow velocity.

357

362

Figure 9 indicates that in the absence of background flow, M 1 consistently compared better to experi-

363

mental results, and this agreement was considered good. On the other hand, M 2 consistently underestimated

364

the air pocket front velocity, and this error increased for larger air pocket sizes. M 2 predictions were more

365

accurate or comparable to M 1 predictions only for cases involving smaller air pocket volumes (V ol∗ = 1.2

366

and V ol∗ = 2.2) in the presence of higher background flows. The accuracy of M 1 appears unaffected by

12

367

different background flows and pocket sizes. These results indicate that the surface tension parameters

368

presented in Wilkinson (1982) for rectangular pipe cross sections were adequate for circular pipes.

369

While experimental observations indicated that the initial stages (t∗ 2, results from M 1 are fairly close to the observations from tested conditions.

375

Results yielded by the M 2 model compared well only with the the smallest air pocket V ol∗ = 1.2.

377

The accuracy of air pocket velocity predictions are reflected in a comparison between measured and √ predicted cavity Froude numbers (F = uf / gD), which are provided in Table 2. The smallest (V ol∗ = 1.2)

378

and largest (V ol∗ = 4.1) air pocket volumes from experiments were selected for the comparison. The front

379

velocity (uf ) was computed from the average velocity between t∗ = 0 and t∗ ≈ 2.5. M 1 was more accurate of

380

the two integral models for each one of the experiments tested in this comparison. This model underpredicted

381

F values by only 0.1% for V ol∗ = 1.2 experiments while overpredicting results in V ol∗ = 4.1 experiments by

382

4.2%. M 2 underpredicted F values by an average of 11.5% with larger discrepancies occuring for larger air

383

pockets and smaller background flows.

384

CONCLUSIONS AND FUTURE WORK

376

385

This work presented results on an investigation of the motion of entrapped air pockets in water pipes.

386

The focus was on the kinematic aspects of air pocket motion and air-water interactions related to the motion

387

of these discrete, finite volume pockets. A long term goal of this research is to incorporate in flow regime

388

transition models the ability to track entrapped pockets and prevent negative impacts related to these

389

events. One relevant context in which the findings of this study can be of practical used is in the rapid

390

filling of stormwater storage tunnels. Entrapped air is a source of various operational issues which range

391

from surging, loss in storage, and uncontrolled air release through water-filled shafts (e.g. geysering). The

392

ability of tracking the motion of entrapped air pockets can be useful in determining the likelihood of such

393

operational problems, and adjust designs that would avoid such adverse issues.

394

An experimental program was proposed where different conditions of pipe slope and velocities rates were

395

tested. Use of clear PVC pipe enabled tracking the coordinates of entrapped pockets over time for each

396

tested condition. In addition to the experimental studies, a non-Boussinesq, integral model was proposed to

397

simulate some of the conditions tested in the experiments. The purpose of this gravity current model was the

398

possibility of such a simple, computationally inexpensive modeling approach being included as a sub-model

13

399

400

in a more complex flow regime transition model. The findings of this research may be summarized as follows:

401

• The celerity of discrete air pockets in stagnant flow conditions depended on air pocket volume, par-

402

ticularly for shallower slopes (below 0.5%), in agreement with Benjamin (1968) and Wilkinson (1982)

403

observations. For stronger slopes, celerity values were not as dependent on volumes as air accumula-

404

tion at the leading edge of the air pocket reduced differences on the pocket thickness;

405

• For horizontal slope and stagnant conditions, the air pocket celerity values slightly decreased over

406

time, indicating that the air pocket spreading decreased the leading edge thickness, resulting in

407

smaller celerity values;

408

• It was observed that the propagation of air pockets in horizontal slope pipes against pipe flows is

409

short lived due to the shearing of the air pockets in the hydraulic jumps that trailed the leading edge

410

of the air pocket;

411

• The celerity of the air pocket leading edge in horizontal and adverse slopes was well approximated

412

by the summation of values observed in quiescent conditions and the pipe flow velocity. Such an

413

approximation did not hold for air pocket propagation in positive slopes where buoyancy forces oppose

414

drag forces;

415

• Air pocket propagation in positive slopes was slowed down by pipe flows, and shearing of the leading

416

edge of the main air pocket by the trailing hydraulic jump was also observed. However, these sheared

417

bubbles regrouped in larger bubbles/pockets downstream from the main air pocket. At a certain

418

size, they would have enough velocity to overcome the drag forces and eventually rejoin the main air

419

pocket.

420

• Predictions of the air pockets leading edge coordinate yielded by proposed numerical models compared

421

well with experiments involving horizontal slopes and pipe flows. The M 1 model alternative (based on

422

the work of Benjamin (1968) and Wilkinson (1982)) yielded better results for large air pocket volumes

423

and smaller pipe flows. The M 2 model alternative results were more accurate for the smallest air

424

pockets volumes (V ol∗ = 1.2) and higher pipe flows.

425

This work also indicates various points that should be addressed in future investigations. One is the

426

development of similar studies using other pipe diameters to assess scale effects on these findings. Exper-

427

iments with 50 mm and 200 mm pipes are planned in the future to complement findings presented in this

428

work. Another important area for future work is to assess the performance of the proposed gravity current

429

integral models when they are embedded in hydraulic models able to simulate both flow regime transition

430

and entrapment of air pockets. Experimental studies that involve both occurrences are planned to help

14

431

assess the effectiveness of such a combination in controlled laboratory conditions.

432

ACKNOWLEDGMENTS

433

434

The authors would like to acknowledge the support provided by LimnoTech Inc., and Auburn University which has funded part of this research.

15

435

APPENDIX I. NOTATION The following symbols are used in this paper.

436

A =

1 2

2

sin α) D4

A2

=

Cross-sectional area below air pocket=(π − α +

C

=

Celerity

C∗

=

Normalized celerity= √CgD

CSO

=

Combined sewer overflow

D

=

Pipeline diameter

F

=

Local Froude number at upstream and/or downstream air pocket front

g

=

Gravitational acceleration

L =

Pipeline length

L1

=

Relative pipeline length coordinate= X L

Lp

=

Air pocket length

Q =

437

Pipeline

area= π4 D2

Flow rate in the pipeline

Q∗

=

Normalized flow rate in the pipeline= √Q

rc

=

Radius of curvature

t =

gD 5

Time

tnew

=

Updated time parameter

told

=

Previous time parameter

t∗

=

q 0 Normalized time= t/X D g

uf

=

Front velocity

Vf low

=

Average flow velocity

Vf∗low

=

low Normalized flow velocity= √fgD

V ol

=

Air pocket volume

V ol∗

=

Normalized air pocket volume=V ol/D3

X

=

Location along pipeline length

X0

=

Initial air pocket length=2.1 m

xf

=

Distance traveled by respective cavity front

X∗

=

Normalized location along pipeline length= XX0

α

=

Pipe half-angle: 2α is the angle subtended from the center of the pipe to the free surface

σ

=

Surface tension

v

16

438

REFERENCES

439

Baines, W.F. (1991). Air cavities as gravity currents on slope. J. Hydr. Eng.,117(12): 1600–1615, 1991.

440

Barnea, D. and Taitel, Y. (1993). Model for slug length distribution in gas-liquid slug flow. Int. J. Multiphase

441

Flow,19(5): 829–838, 1993.

442

Benjamin, T.B. (1968). Gravity currents and related phenomena. J. Fluid Mech., 31: 209–248, 1968.

443

DeHenau, V. and Raithby, G.D. (1995). Transient two-fluid model for the simulation of slug flow in pipelines.

444

Int. J. Multiphase Flow,21(3): 335–349, 1995.

445

Falvey, H.T. (1980). Air-water flow in hydraulic structures. USBR Engineering Monograph, n. 41

446

Glauser, S. and Wickenhauser, M. (2009). Bubble movement in downward-inclined pipes. J. Hydr. Eng.,

447

448

449

450

451

452

453

454

455

456

457

458

459

460

461

462

463

135(11): 1012–1015, 2009. Guo, Q. and Song, C.S.S. (1991). Dropshaft Hydrodynamics under Transient Conditions. J. Hydr. Eng., 117(8): 1042–1055, 1991. Hallworth, M.A., Hogg, A.J., and Huppert, H.E. (1998). Effects of external flow on compositional and particle gravity currents. J. Fluid Mech., 359: 109–142, 1998. Hamam, M.A. and McCorquodale, A. (1982). Transient conditions in the transition from gravity to surcharged sewer flow. Can. J. Civ. Engrg, (9): 189–196, 1982. Issa, R.I. and Kempf, M.H.W. (2003). Simulation of slug flow in horizontal and nearly horizontal pipes with the two-fluid model. Int. J. Multiphase Flow, 29(1): 69–95, 2003. Lauchlan, C.S., Escarameia, M., May, R.W.P., Burrows, R. and Gahan, C.(2005). Air in pipelines - A literature review. HR Wallingford Technical Report SR., 649. Li, J. and McCorquodale, A. (1999). Modeling Mixed Flow in Storm Sewers. J. Hydr. Eng., 125(11): 1170–1180, 1999. Lingireddy, S., Wood, D.J. and Zloczower, N. (2004). Pressure surges in pipeline systems resulting from air releases. J. AWWA, 96(7): 88–94, 2004. Little, M.J., Powell, J.C. and Clark, P.B. (2008). Air movement in water pipelines - some new developments. Proc., 10th Int. Conf. on Pressure Surges, BHR, Edinburgh, UK, 111–122, 2008.

17

464

465

466

467

468

469

470

471

472

473

474

475

Martin, C.S. (1976). Entrapped air in pipelines. Proc., 2nd Int. Conf. on Pressure Surges, BHR, Bedford, England, 15–28, 1976. Nielsen, K.D., and Davis, A.L. (2009). Air migration analysis of the Terror Lake tunnel. Proc., 33rd IAHR Congress, IAHR, Vancouver, BC, 262–268, 2009. Pothof, I. and Clemens, F.H.L.R. (2010). On elongated air pockets in downward sloping pipes. J. Hydr. Res., 48(4): 499–503, 2010. Pothof, I. and Clemens, F..H.L.R. (2011). Experimental study of air-water flow in downward sloping pipes. Int. J. Multiphase Flow, 37: 278–292, 2011. Pozos, O.E. (2007). Investigation on the effects of entrained air in pipelines. PhD Thesis, University of Stuttgart, 2007. Pozos, O., Gonzalez, C., Giesecke, J., Marx, W. and Rodal, E. (2010). Air entrapped in gravity pipeline systems. J. Hydr. Res., 49(3): 394–397, 2010.

476

Simpson, J.E. (1997). Gravity currents in the environment and the laboratory. Cambridge University Press..

477

Song, C.S.S., Cardle, J.A. and Leung, K.S. (1983). Transient mixed-flow models for storm sewers. J. Hydr.

478

479

480

Eng., 109(11): 1487–1504, 1983. Toro, E. (2009). Riemann solvers and numerical methods for fluid dynamics: a practical introduction. Springer Verlag.

481

Townson, J.M.(1991) Free-surface hydraulics. Cambridge University Press.

482

Trindade, B.C. and Vasconcelos, J.G. (2013) Modeling of water pipelines filling events accounting for air

483

phase effects. J. Hydr. Eng., (): 1943–7900, 2013.

484

Ungarish, M.(2009) An introduction to gravity currents and intrusions. CRC Press, 2009.

485

Vasconcelos, J.G. and Leite, G.M. (2012). Pressure surges following sudden air pocket entrapment in

486

487

488

489

490

stormwater tunnels. J. Hydr. Eng., 138(12): 1080–1089, 2012. Vasconcelos, J.G. and Wright, S.J. (2006). Mechanisms for air pocket entrapment in stormwater storage tunnels. Proc., 2006 ASCE EWRI Congress, Omaha, NE, 2006. Vasconcelos, J.G. and Wright, S.J. (2008). Rapid flow startup in filled horizontal pipelines. J. Hydr. Eng., 134(7): 984–992, 2008.

18

491

492

493

494

Vasconcelos, J.G. and Wright, S.J. (2009). Rapid filling of poorly ventilated stormwater storage tunnels. J. Hydr. Res., 47(5): 547–558, 2009. Vasconcelos, J.G. and Wright, S.J. (2011). Geysering generated by large air pockets released through waterfilled ventilation shafts. J. Hydr. Eng., 135(5): 543–555, 2011.

495

Wilkinson, D.L. (1982). Motion of air cavities in long horizontal ducts. J. Fluid Mech., 118: 109–122, 1982.

496

Wright, S.J., Lewis, J.W. and Vasconcelos, J.G. (2011). Geysering in rapidly filling storm-water tunnels. J.

497

498

499

Hydr. Eng., 137(5): 543–555, 2011. Zhou, F., Hicks, F.E. and Steffler, P.M. (2002). Transient flow in a rapidly filling horizontal pipe containing trapped air. J. Hydr. Eng., 128(6): 625–634, 2002.

19

500 501 502 503

List of Tables 1 Experimental variables considered √ with respective tested ranges . . . . . . . . . . . . . . . . . 2 Cavity Froude number (F = uf / gD) comparison for the largest and smallest air pocket volumes tested with various background flows. . . . . . . . . . . . . . . . . . . . . . . . . . . .

20

21 22

TABLE 1.

Experimental variables considered with respective tested ranges

Experimental variable p Flow rate (normalized by Q/ gD5 ) Pipeline slope

Air pocket volume (normalized by D3 )

Range Three values ranging from Q∗ =0 up to 0.388 Horizontal slope 0.5%, 1.0%, and 2.0% adverse slopes 0.5%, 1.0%, and 2.0% positive slopes Horizontal: up to 4 values in the range V ol∗ =1.2-4.1 Adverse: up to 4 values in the range V ol∗ =0.48-3.8 Positive: up to 4 values in the range V ol∗ =0.48-3.8

21

√ TABLE 2. Cavity Froude number (F = uf / gD) comparison for the largest and smallest air pocket volumes tested with various background flows. p Q∗ = Q/ gD5 0 0.13 0.27 0.39 V ol∗ = 1.2 Experiments M1 prediction M2 prediction

0.284 0.268 0.253

Experiments M1 prediction M2 prediction

0.419 0.427 0.323

0.451 0.467 0.431 V ol∗ 0.594 0.612 0.496

22

0.672 0.664 0.617 = 4.1 0.751 0.809 0.683

0.804 0.826 0.773 N/A 0.972 0.840

504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532

List of Figures 1 Sketch of apparatus used in experiments. For positive slopes, the length of the reaches upstream and downstream from the valve were adjusted to 6.0 m and 3.0 m respectively. . . . . 2 Schematic Diagram of an air pocket propagating in a circular pipe using the integral model assumption. The direction of uf D is always positive while the sign of uf U will change depending on the magnitude of U and y/D. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 Trajectory of the air pocket leading edge and observed celerity for horizontal slope and various air pocket volumes, no pipe flow. Negative coordinates indicate pockets are propagating toward upstream. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4 Trajectory of the air pocket leading edge and observed celerity for various adverse slopes and air pocket volumes, no pipe flows. Negative coordinates indicate pockets are propagating toward upstream. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 Trajectory of the air pocket leading edge, observed celerity and relative celerity for conditions with horizontal slope and pipe flows for V ol∗ =1.2 and V ol∗ =2.2. . . . . . . . . . . . . . . . . 6 Sequence of snapshots illustrating the trajectory of the backward moving pocket trajectory for Q∗ = 0.26, V ol∗ =4.1 and horizontal slope, illustrating the shearing process. Interface between air and water is enhanced. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7 Trajectory of the air pocket leading edge, observed celerity and relative celerity for conditions with 1% adverse slope and pipe flows for V ol∗ =0.95 and V ol∗ =1.9. . . . . . . . . . . . . . . . 8 Trajectory of the air pocket leading edge, observed celerity and relative celerity for conditions with 1% positive slope and pipe flows for V ol∗ =0.95 and V ol∗ =1.9. . . . . . . . . . . . . . . . 9 Air pocket front trajectory comparison between experiments and both integral models for various background flows and pocket volumes: a) V ol∗ = 1.2, b) V ol∗ = 2.2, c) V ol∗ = 3.1, and d) V ol∗ = 4.1. The solid lines represent the integral model predictions, and the data markers represent experimental values. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10 Air pocket front celerity comparison between experiments and both integral models for various background flows and pocket volumes: a) V ol∗ = 1.2, b) V ol∗ = 2.2, c) V ol∗ = 3.1, and d) V ol∗ = 4.1. The solid lines represent the integral model predictions, and the data markers represent experimental values. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

23

24

25

26

27 28

29 30 31

32

33

FIG. 1. Sketch of apparatus used in experiments. For positive slopes, the length of the reaches upstream and downstream from the valve were adjusted to 6.0 m and 3.0 m respectively.

24

FIG. 2. Schematic Diagram of an air pocket propagating in a circular pipe using the integral model assumption. The direction of uf D is always positive while the sign of uf U will change depending on the magnitude of U and y/D.

25

FIG. 3. Trajectory of the air pocket leading edge and observed celerity for horizontal slope and various air pocket volumes, no pipe flow. Negative coordinates indicate pockets are propagating toward upstream.

26

FIG. 4. Trajectory of the air pocket leading edge and observed celerity for various adverse slopes and air pocket volumes, no pipe flows. Negative coordinates indicate pockets are propagating toward upstream.

27

FIG. 5. Trajectory of the air pocket leading edge, observed celerity and relative celerity for conditions with horizontal slope and pipe flows for V ol∗ =1.2 and V ol∗ =2.2.

28

FIG. 6. Sequence of snapshots illustrating the trajectory of the backward moving pocket trajectory for Q∗ = 0.26, V ol∗ =4.1 and horizontal slope, illustrating the shearing process. Interface between air and water is enhanced.

29

FIG. 7. Trajectory of the air pocket leading edge, observed celerity and relative celerity for conditions with 1% adverse slope and pipe flows for V ol∗ =0.95 and V ol∗ =1.9.

30

FIG. 8. Trajectory of the air pocket leading edge, observed celerity and relative celerity for conditions with 1% positive slope and pipe flows for V ol∗ =0.95 and V ol∗ =1.9.

31

FIG. 9. Air pocket front trajectory comparison between experiments and both integral models for various background flows and pocket volumes: a) V ol∗ = 1.2, b) V ol∗ = 2.2, c) V ol∗ = 3.1, and d) V ol∗ = 4.1. The solid lines represent the integral model predictions, and the data markers represent experimental values.

32

FIG. 10. Air pocket front celerity comparison between experiments and both integral models for various background flows and pocket volumes: a) V ol∗ = 1.2, b) V ol∗ = 2.2, c) V ol∗ = 3.1, and d) V ol∗ = 4.1. The solid lines represent the integral model predictions, and the data markers represent experimental values.

33