# Flow filling a curved pipe

Hervé Michallet, C. Mathis, P. Maïssa, F. Dias

► To cite this version:

Hervé Michallet, C. Mathis, P. Maïssa, F. Dias. Flow filling a curved pipe. Journal of Fluids Engineering, 2001, 123 (3), pp.686-691. 10.1115/1.1374442. hal-00260280

HAL Id: hal-00260280

https://hal.science/hal-00260280v1

Submitted on 18 Feb 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# Flow Filling a Curved Pipe

A small scale experiment was designed to study the propagation of the front of a viscous fluid filling a curved pipe. Several Newtonian fluids with different viscosities and a non-Newtonian fluid have been used. The experiments show that there exists a minimum speed for completely filling the pipe, which depends on the parameters of the experiment (di-ameter d and radius of curvature R of the pipe, kinematic viscosity v of the fluid). Appropriate dimensionless numbers are introduced to characterize the flow and optimal filling conditions.

H. Michallet

C. Mathis<fnref n="1" />

P. Maïssa

F. Dias<fnref n="2" />

Institut Non-Linéaire de Nice,
UMR 6618 - CNRS and UNSA,
1361, route des Lucioles, 06560 Valbonne,
France

## Introduction

Prestressed concrete is commonly used in construction. To prevent the prestressed strands from any corrosion, a cement solution is generally injected in pipes of complex geometry. The pipe has to be completely filled in order to avoid weak zones in the structure. The aim of the present study is to observe experimentally the propagation of the front of a viscous fluid injected in a curved pipe, in order to understand the reasons for the appearance of unfilled zones. The solution used in construction is composed of cement, water, liquefier, and retarder. It is a non-Newtonian fluid with air bubbles and particles in suspension. As a first step, we mainly focus in this paper on the filling by Newtonian fluids.

Permanent flows in pipes have been widely studied [1], particularly in the case of curved pipes [2–4] and in engineering applications [5,6]. However, pipe flows involving a propagating front or interface remain poorly understood. The problem of an air-water interface propagating in straight inclined pipes has been addressed only recently [7,8]. As stated above, in industrial applications such as construction, the fluid is often mixed with air bubbles and the flow ought to be treated as a slug flow (see for example [9,10] for recent work on slug flows in pipes). But the emphasis of the present work is in the propagation of the interface which clearly separates the gas from the liquid.

The front propagation of a perfect fluid in a horizontal pipe was studied by Benjamin [11]. See also the review by Simpson [12] on gravity currents for related studies and Asavenant and Vandenbergroeck [13]. Benjamin analyzed the front in terms of a ‘cavity flow’ displacing a fluid beneath it. The two-dimensional geometry is shown in Fig. 1(a) (Benjamin also considered the case of a circular cross-section). Note that Benjamin studied only the case where the free surface detaches with an angle of 60 degrees. By applying conservation of mass, momentum, and energy, he found that there is a unique solution, characterized by

$$ \frac{U}{\sqrt{gd}} = \frac{1}{2}, \quad \frac{u}{\sqrt{gh}} = \sqrt{2}, \quad \frac{h}{d} = \frac{1}{2}, $$

where the meaning of the various symbols is shown in Fig. 1. In other words, the cavity fills half of the box. An experiment was suggested by Benjamin to realize closely this flow. Liquid initially fills a long rectangular box closed at both ends and fixed horizontally. One end is then opened, and under the action of gravity the liquid flows out freely from this end. It can be expected that, after

the transient effects of starting have disappeared, the air-filled cavity replacing the volume of the ejected fluid will progress steadily along the box. Observed in a frame of reference travelling with the front of the cavity, the motion of the liquid will appear to be steady, as shown in Fig. 1(a). If the effects of viscosity and surface tension are neglected, the velocity of the cavity relative to a stationary observer will be U, and, since h=d/2, the liquid will discharge from the open end with the same velocity. In addition to the solution studied by Benjamin, there is also a one-parameter family of solutions in which the free surface detaches from the box tangentially [14]. Such a solution is plotted in Fig. 1(b). At the exit of the box, the Froude number $u(gh)^{-1/2}$ for these solutions is between 1 and $\sqrt{2}$. The limit $\sqrt{2}$ corresponds to Benjamin’s solution, while the limit 1 corresponds to the vanishing of the cavity. Using conservation of mass, momentum and energy, one finds the relation

$$ \frac{u}{\sqrt{gh}} = \sqrt{\frac{d}{h}}. \qquad (1) $$

Unfortunately, simple arguments based on conservation of mass, momentum and energy do not give as much information in the case of a curved channel.

It is nevertheless important to keep in mind that the Froude number prevails to determine the size of the cavity and thus the ability to fill the pipe. This point will be discussed in the section where we consider various dimensionless numbers. Before that, we first describe the manner in which the experiments have been conducted. Then we present and discuss the experimental results.

## Experimental Setup

The experimental setup is shown in Fig. 2. Gravity acts downwards. A transparent PVC pipe (‘Tubclair’) of inner diameter d is curved with a radius R. Different diameters and radii were considered: $0.6 \le d \le 2$ cm; $5 \le R \le 30$ cm. The special case of a horizontal straight pipe was also investigated. The fluid was injected with a centrifugal pump at a constant mean flow velocity (noted U).

Several fluids were considered in order to investigate the effects of density ρ, surface tension σ, and mainly kinematic viscosity ν (see Table 1). A series of runs has been performed with a non-Newtonian fluid (gel). Its viscosity has been measured with a plane-plane viscosimeter (see Fig. 3).

We assume that capillary effects are negligible, so that $d \gg L_c$ where $L_c = \sigma^{1/2}(\rho g)^{-1/2}$ is the capillary length ($\sigma = 73 \text{ gs}^{-2}$ for water and $\sigma = 21 \text{ gs}^{-2}$ for silicone). In the case of water $L_c = 0.26$ cm and this assumption might be a shortcoming for small diameters. We nevertheless performed a few runs with water contaminated with wetting agents in order to diminish $L_c$: the behavior and the results were the same within measurements errors. Moreover, the results presented in this paper are limited to $d \ge 1$ cm.

<footnote n="1">¹Permanent address: Laboratoire des Ecoulements Géophysiques et Industriels (UJF-INPG-CNRS), BP53, 38041 Grenoble Cedex 9, Fax: +33 4 76 82 52 71; E-mail: herve.michallet@hmg.inpg.fr.</footnote>

<footnote n="2">²Corresponding author, present address: Ecole Normale Supérieure de Cachan - CNRS, CMLA UMR8536, 61 avenue du Président Wilson, 94235 Cachan Cedex; E-mail: dias@cmla.ens-cachan.fr.</footnote>
<figure><img src="image_2.png" /><figcaption>Fig. 1 Steady two-dimensional flow past a cavity. (a) Solution with the free surface leaving the wall with a 60 degree angle. (b) Solution with the free surface leaving the wall tangentially. The point C denotes the detachment point.</figcaption></figure>

<figure><img src="image_3.png" /><figcaption>Fig. 3 Dynamic viscosity η of the gel versus the shear stress ˙γ. Its density is ρ=1.03 g/cm³. We note that this fluid is very viscous in the conditions of our experiments (we can estimate roughly that ˙γ~2U/d is less than 10 s⁻¹). Measurements have been performed with a plane-plane viscosimeter at Ecole des Mines de Paris (CEMEF, Sophia Antipolis).</figcaption></figure>

## Observations

For a very low flow rate (a typical mean flow velocity is $U \sim 0.1$ cm/s) and a sufficiently large pipe section (so that $d \gg L_c$) the interface is more or less horizontal in the upflow part of the pipe. The flow is dominated by gravity. In the downflow part of the pipe the liquid is creeping down the pipe.

If $U$ is slightly increased, the interface becomes more or less perpendicular to the pipe when rising. At the top, the interface deforms as shown in Fig. 2. The intersection between the interface and the upper part of the pipe (point C in Fig. 2) is a stagnation point.

If $U$ is increased, the point C moves (interface in dotted line in Fig. 2). As soon as $U$ is decreased back to zero, the pipe is emptying and the interface moves back to the top.

<figure><img src="image_4.png" /><figcaption>Fig. 2 Diagram of the flow and notation. The interface is shown at two different times in solid and dotted lines.</figcaption></figure>

<table><caption>Table 1 Fluids used in the experiments.</caption><thead><tr><th>Fluid</th><th>Density<br>ρ (g/cm³)</th><th>Viscosity<br>ν (cm²/s)</th><th>R=5 cm</th><th>Symbol<br>R=30 cm</th><th>R=∞</th></tr></thead><tbody><tr><td>Water</td><td>1</td><td>0.01</td><td>●</td><td>○</td><td>○</td></tr><tr><td>Water (40 %)<br>+ glycerol (60 %)</td><td>1.16</td><td>0.086</td><td>▼</td><td>▽</td><td>▽</td></tr><tr><td>Silicone V50</td><td>1</td><td>0.5</td><td>►</td><td>▽</td><td>▽</td></tr><tr><td>Silicone V100</td><td>1</td><td>1</td><td></td><td>▽</td><td>▽</td></tr><tr><td>Silicone V300</td><td>1</td><td>3</td><td></td><td>▽</td><td>▽</td></tr><tr><td>Gel</td><td>1.03</td><td>(See Fig. 3)</td><td></td><td>□</td><td>□</td></tr></tbody></table>

We consider that the pipe is completely filled when the point C has reached the point B (Fig. 2). We define $U_c$ as the critical speed beyond which the interface moves until the pipe is completely filled. The speed $v$ is the mean velocity of the point C moving from A to B (Fig. 2): if $U$ is larger than $U_c$, then $v$ is greater than zero.

For the special case of the non-Newtonian fluid gel, the viscosity is smaller where the shear stress is larger, that is near the boundaries. Away from the boundaries, the behavior of the fluid is similar to a solid core pushed by the flow. In other words, the characteristic time of deformation of the interface is very long.

## Measurements

The mean flow velocity $U$ was obtained by measuring the flow rate in the steady-state regime, i.e., when the pipe was completely filled. The flow rate and $v$ were estimated by measuring the filling time of a graded beaker and the travelling time of point C from A to B (Fig. 2) with a stopwatch. We may consider that these simple techniques do not lead to any bias error. By repeating the same run (at least five times) for different flow rates, we have estimated that the precision limit and thus the uncertainty was always better than 10 percent on $U$ and $v$.

We measured $v$ as a function of $U$ for different experiments using various fluids, pipes, and curvatures (see Table 1). In Fig. 4, $v$ is plotted versus $U$ for several runs in four different configurations. The line $v=U$ delimits the region where physical solutions exist: $v$ cannot be larger than $U$. In all our experiments $v$ apparently varies linearly with $U$. We may thus determine the critical velocity $U_c$ (for $v=0$), the slope $\alpha=dv/dU$ and the interfacial velocity $v(U=0)$, by rms fitting a straight line into the data for each fluid, pipe diameter and radius of curvature.

Let us first consider the slope $\alpha$ of the four straight lines in Fig. 4. The analysis presented in the Introduction considered a perfect fluid. It was assumed that the flow is identical in a reference frame moving with the front. In this framework, the filling of a pipe is similar to the emptying of the pipe. This would imply that $\alpha$ is equal to one: adding $\delta U$ to the mean entrance flow velocity would add the same value to the front propagation velocity. It is therefore not surprising that the slopes of the straight lines in Fig. 4 related to water are close to one. When the fluid is far from perfect
<figure><img src="image_5.png" /><figcaption>Fig. 4 Speed of the front propagation v versus the mean entrance velocity U (d=1.5 cm; see list of symbols in Table 1). The straight lines are fitted curves by quadratic means.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>Fig. 5 Critical speed versus the pipe diameter (see list of symbols in Table 1)</figcaption></figure>

but still Newtonian (case of the silicone oil V100 in Fig. 4), α is clearly smaller than one: when the viscosity is large, the deformation of the front strongly depends on U. For a low viscosity, the deformation of the interface takes place on a very short time scale and the propagation of the front apparently does not depend on the reference frame. In the other limit, that is for a very large viscosity, the deformation of the front takes place on a very long time scale and the slope is again close to one (this is the case of the non-Newtonian fluid).

The velocity $v(0)$ corresponds to the draining of the pipe. The measurements have been performed by bringing $U$ to zero after filling the pipe: the flow is then from right to left in Fig. 2. It is important to emphasize that $v(0)$ is very small in the case of the non-Newtonian fluid. We have estimated that the draining velocity is about 1 cm per week. The deformation of the front is visible only when the fluid is sheared and $U_c$ is therefore also close to zero for the gel. The agreement in $v(0)$ between the measured values and the values deduced from the fitted straight lines in Fig. 4 is rather good except in the case of water for $R=30$ cm. The discrepancy is probably due to measurement errors but also to wetting effects. The contact between the fluid and the PVC is indeed not the same when filling or emptying the pipe. It was not our purpose to investigate wetting effects and we did not test other materials for the pipe. We therefore have to consider that wetting effects may induce a bias error over $U_c$ of about 10 percent in the worst case, i.e., an uncertainty of 20 percent in the water case.

The critical speed $U_c$ (mean entrance velocity for which $v=0$) is plotted against the pipe diameter $d$ in Fig. 5. $U_c$ clearly tends to increase with $d$. Considering gravitational effects, the larger is $d$, the more efficient are gravity effects to deform the interface. It is thus not surprising that inertial effects have to be increased to compensate for it. Gravity effects could also reasonably depend on the radius of curvature. The interface deforms more rapidly when the pipe is curved: for a same diameter, $U_c$ is therefore smaller when $R=\infty$. Otherwise, if $R$ is large but not infinite, the travelling time of the front in the pipe is large and the mean entrance velocity has to be large for the front to propagate. In other words, for a very small radius, the front has not enough time to deform before the exit of the pipe. This would explain the fact that $U_c$ is larger for $R=30$ cm than for $R=5$ cm.

In addition, above $v=1$ cm$^2$/s, viscous effects tend to maintain the shape of the interface and thus decrease $U_c$. This is noticeable

for the silicone V100 and clearly visible for the silicone V300. Again the case of the gel requires specific comments. $U_c$ is very close to zero because the front deforms only when the fluid is sheared as mentioned above. The graphical estimation of $U_c$ gives a nonzero value only for $d=1.5$ cm. This configuration will allow us to estimate dimensionless numbers for the gel in the next section.

The data in Figs 4 and 5 are available in tabular form in Table 2.

## Interpretation in Terms of Dimensionless Numbers

We may consider various dimensionless numbers that could characterize the flow in our special configuration. Note that in our experiments both the Bond number ($Bo=d^2\rho g/\sigma$) and the Weber number ($We=\rho dU^2/\sigma$) are very large: capillary effects are a lot weaker than inertial and gravitational effects.

Considering inertial and viscous effects in the framework of a curved pipe, we introduce a modified Reynolds number $Re$, which partly takes into account centrifugal effects, and is in fact the classical Reynolds number plus the Dean number [15]:

$$Re = \frac{Ud}{\nu} \left( 1 + \sqrt{\frac{d}{R}} \right). \qquad (2)$$

To balance inertial and gravity effects, we build a modified Froude number $Fr$, such as:

$$Fr = \frac{U}{\sqrt{gd}} \left( 1 + \sqrt{\frac{d}{R}} \right). \qquad (3)$$

The modification of the classical dimensionless numbers may be understood by considering that there is a need to add centrifugal effects to maintain the fluid against the concave side of the pipe. Note that (2) and (3) tend toward the classical Reynolds and Froude numbers in the limiting case of a straight pipe, i.e., when $R$ tends toward infinity.

Let us finally define $K$, which measures the balance between gravity and viscosity:

$$K = \frac{Fr}{Re} = \frac{\nu}{\sqrt{gd^3}}. \qquad (4)$$
<table><caption>**Table 2 Measured values**</caption><thead><tr><th>Fluid</th><th>Radius<br>of curvature<br>R (cm)</th><th>Pipe<br>diameter<br>d (cm)</th><th>Speed<br>of front<br>v (cm/s)</th><th>Mean entrance<br>velocity<br>U (cm/s)</th></tr></thead><tbody><tr><td rowspan="10">► Silicone V50</td><td>5</td><td>1.0</td><td>-8.6</td><td>0.0</td></tr><tr><td>5</td><td>1.0</td><td>0.6</td><td>10.2</td></tr><tr><td>5</td><td>1.0</td><td>3.2</td><td>14.1</td></tr><tr><td>5</td><td>1.0</td><td>11.8</td><td>28.7</td></tr><tr><td>5</td><td>1.0</td><td>18.9</td><td>36.6</td></tr><tr><td>5</td><td>1.2</td><td>-10.7</td><td>0.0</td></tr><tr><td>5</td><td>1.2</td><td>4.4</td><td>19.8</td></tr><tr><td>5</td><td>1.2</td><td>6.7</td><td>22.8</td></tr><tr><td>5</td><td>1.2</td><td>12.2</td><td>31.4</td></tr><tr><td rowspan="10">Δ</td><td>30</td><td>1.0</td><td>-8.8</td><td>0.0</td></tr><tr><td>30</td><td>1.0</td><td>1.0</td><td>13.0</td></tr><tr><td>30</td><td>1.0</td><td>4.9</td><td>18.9</td></tr><tr><td>30</td><td>1.0</td><td>8.9</td><td>23.6</td></tr><tr><td>30</td><td>1.0</td><td>19.5</td><td>35.6</td></tr><tr><td>30</td><td>1.0</td><td>22.1</td><td>38.4</td></tr><tr><td>30</td><td>1.2</td><td>-10.8</td><td>0.0</td></tr><tr><td>30</td><td>1.2</td><td>1.4</td><td>19.1</td></tr><tr><td>30</td><td>1.2</td><td>4.5</td><td>24.0</td></tr><tr><td>30</td><td>1.2</td><td>8.9</td><td>31.0</td></tr><tr><td rowspan="10">Δ silicone V100</td><td>30</td><td>1.2</td><td>-13.6</td><td>0.0</td></tr><tr><td>30</td><td>1.5</td><td>2.2</td><td>24.2</td></tr><tr><td>30</td><td>1.5</td><td>7.3</td><td>34.2</td></tr><tr><td>30</td><td>1.5</td><td>16.4</td><td>49.2</td></tr><tr><td>30</td><td>1.5</td><td>22.5</td><td>57.0</td></tr><tr><td>30</td><td>2.0</td><td>-18.0</td><td>0.0</td></tr><tr><td>30</td><td>2.0</td><td>2.0</td><td>34.6</td></tr><tr><td>30</td><td>2.0</td><td>3.4</td><td>38.4</td></tr><tr><td>30</td><td>2.0</td><td>4.8</td><td>41.6</td></tr><tr><td>30</td><td>2.0</td><td>6.4</td><td>46.3</td></tr><tr><td rowspan="10">◄ Silicone V300</td><td>30</td><td>1.2</td><td>-8.5</td><td>0.0</td></tr><tr><td>30</td><td>1.2</td><td>2.3</td><td>14.7</td></tr><tr><td>30</td><td>1.2</td><td>2.3</td><td>14.7</td></tr><tr><td>30</td><td>1.2</td><td>3.9</td><td>17.1</td></tr><tr><td>30</td><td>1.2</td><td>5.0</td><td>18.0</td></tr><tr><td>30</td><td>1.5</td><td>-11.0</td><td>0.0</td></tr><tr><td>30</td><td>1.5</td><td>1.3</td><td>19.6</td></tr><tr><td>30</td><td>1.5</td><td>6.8</td><td>29.4</td></tr><tr><td>30</td><td>1.5</td><td>10.3</td><td>35.3</td></tr><tr><td>30</td><td>1.5</td><td>11.8</td><td>38.5</td></tr><tr><td rowspan="10">● water</td><td>30</td><td>1.5</td><td>-6.2</td><td>0.0</td></tr><tr><td>30</td><td>1.5</td><td>0.9</td><td>10.5</td></tr><tr><td>5</td><td>1.0</td><td>-9.5</td><td>0.0</td></tr><tr><td>5</td><td>1.0</td><td>3.0</td><td>12.4</td></tr><tr><td>5</td><td>1.0</td><td>7.2</td><td>15.5</td></tr><tr><td>5</td><td>1.0</td><td>10.8</td><td>18.1</td></tr><tr><td>5</td><td>1.2</td><td>-10.3</td><td>0.0</td></tr><tr><td>5</td><td>1.2</td><td>0.7</td><td>13.7</td></tr><tr><td>5</td><td>1.2</td><td>1.4</td><td>14.0</td></tr><tr><td>5</td><td>1.2</td><td>7.0</td><td>19.1</td></tr><tr><td rowspan="10">○</td><td>5</td><td>1.5</td><td>-11.9</td><td>0.0</td></tr><tr><td>30</td><td>1.0</td><td>-9.8</td><td>0.0</td></tr><tr><td>30</td><td>1.0</td><td>3.2</td><td>14.8</td></tr><tr><td>30</td><td>1.0</td><td>9.6</td><td>19.4</td></tr><tr><td>30</td><td>1.0</td><td>13.5</td><td>23.1</td></tr><tr><td>30</td><td>1.0</td><td>41.3</td><td>47.1</td></tr><tr><td>30</td><td>1.2</td><td>-12.0</td><td>0.0</td></tr><tr><td>0</td><td>1.2</td><td>0.5</td><td>16.9</td></tr><tr><td>30</td><td>1.2</td><td>4.6</td><td>19.5</td></tr><tr><td>30</td><td>1.2</td><td>5.6</td><td>22.5</td></tr><tr><td>30</td><td>1.2</td><td>9.2</td><td>23.2</td></tr><tr><td rowspan="10">◎</td><td>30</td><td>1.2</td><td>14.1</td><td>29.0</td></tr><tr><td>30</td><td>1.2</td><td>28.0</td><td>42.3</td></tr><tr><td>30</td><td>1.2</td><td>28.8</td><td>44.9</td></tr><tr><td>30</td><td>1.2</td><td>52.7</td><td>66.6</td></tr><tr><td>30</td><td>1.5</td><td>-15.0</td><td>0.0</td></tr><tr><td>30</td><td>1.5</td><td>9.9</td><td>28.9</td></tr><tr><td>30</td><td>1.5</td><td>15.0</td><td>33.1</td></tr><tr><td>30</td><td>1.5</td><td>14.9</td><td>33.1</td></tr><tr><td>30</td><td>1.5</td><td>22.5</td><td>41.1</td></tr><tr><td>30</td><td>1.5</td><td>32.6</td><td>49.2</td></tr><tr><td>30</td><td>1.5</td><td>55.2</td><td>72.2</td></tr><tr><td>∞</td><td>1.5</td><td>-12.0</td><td>0.0</td></tr><tr><td>∞</td><td>1.5</td><td>0.0</td><td>9.4</td></tr><tr><td>∞</td><td>1.5</td><td>0.0</td><td>13.3</td></tr><tr><td>∞</td><td>1.5</td><td>7.1</td><td>17.8</td></tr><tr><td>∞</td><td>1.5</td><td>9.2</td><td>20.6</td></tr><tr><td>∞</td><td>1.5</td><td>13.0</td><td>25.2</td></tr></tbody></table>

<table><caption>**Table 2 (Continued)**</caption><thead><tr><th>Fluid</th><th>Radius<br>of curvature<br>R (cm)</th><th>Pipe<br>diameter<br>d (cm)</th><th>Speed<br>of front<br>v (cm/s)</th><th>Mean entrance<br>velocity<br>U (cm/s)</th></tr></thead><tbody><tr><td>-</td><td>∞</td><td>1.5</td><td>14.1</td><td>25.2</td></tr><tr><td>-</td><td>∞</td><td>1.5</td><td>18.3</td><td>29.8</td></tr><tr><td>-</td><td>∞</td><td>1.5</td><td>20.0</td><td>29.8</td></tr><tr><td>-</td><td>∞</td><td>1.5</td><td>29.0</td><td>41.8</td></tr><tr><td>-</td><td>∞</td><td>1.5</td><td>30.0</td><td>41.8</td></tr><tr><td>▼ water (40%)<br>+ glyc. (60%)</td><td>5</td><td>1.2</td><td>1.6</td><td>16.7</td></tr><tr><td>-</td><td>5</td><td>1.2</td><td>4.8</td><td>20.2</td></tr><tr><td>-</td><td>5</td><td>1.2</td><td>12.6</td><td>27.5</td></tr><tr><td>▽</td><td>30</td><td>1.2</td><td>6.1</td><td>26.7</td></tr><tr><td>-</td><td>30</td><td>1.2</td><td>15.0</td><td>38.8</td></tr><tr><td>-</td><td>30</td><td>1.2</td><td>14.9</td><td>38.8</td></tr><tr><td>-</td><td>30</td><td>1.2</td><td>29.5</td><td>63.2</td></tr><tr><td>-</td><td>30</td><td>1.2</td><td>51.6</td><td>95.7</td></tr><tr><td>-</td><td>30</td><td>2.0</td><td>5.7</td><td>39.3</td></tr><tr><td>-</td><td>30</td><td>2.0</td><td>3.3</td><td>42.4</td></tr><tr><td>-</td><td>30</td><td>2.0</td><td>10.2</td><td>56.8</td></tr><tr><td>-</td><td>30</td><td>2.0</td><td>15.2</td><td>64.2</td></tr><tr><td>-</td><td>∞</td><td>2.0</td><td>0.0</td><td>6.9</td></tr><tr><td>▽</td><td>∞</td><td>2.0</td><td>0.0</td><td>8.2</td></tr><tr><td>-</td><td>∞</td><td>2.0</td><td>16.4</td><td>24.8</td></tr><tr><td>-</td><td>∞</td><td>2.0</td><td>22.7</td><td>32.1</td></tr><tr><td>-</td><td>∞</td><td>2.0</td><td>34.9</td><td>47.1</td></tr><tr><td>-</td><td>∞</td><td>2.0</td><td>45.3</td><td>69.2</td></tr><tr><td>-</td><td>∞</td><td>2.0</td><td>57.0</td><td>80.4</td></tr><tr><td>□ gel</td><td>30</td><td>1.2</td><td>0.0</td><td>0.0</td></tr><tr><td>-</td><td>30</td><td>1.2</td><td>1.8</td><td>2.0</td></tr><tr><td>-</td><td>30</td><td>1.2</td><td>20.3</td><td>23.1</td></tr><tr><td>-</td><td>30</td><td>1.2</td><td>22.3</td><td>27.8</td></tr><tr><td>-</td><td>30</td><td>1.5</td><td>0.0</td><td>0.0</td></tr><tr><td>-</td><td>30</td><td>1.5</td><td>6.2</td><td>9.3</td></tr><tr><td>-</td><td>30</td><td>1.5</td><td>6.5</td><td>9.6</td></tr><tr><td>-</td><td>30</td><td>1.5</td><td>11.1</td><td>15.1</td></tr><tr><td>-</td><td>30</td><td>1.5</td><td>13.7</td><td>18.5</td></tr><tr><td>-</td><td>30</td><td>1.5</td><td>42.3</td><td>46.0</td></tr><tr><td>-</td><td>30</td><td>2.0</td><td>0.0</td><td>0.0</td></tr><tr><td>-</td><td>30</td><td>2.0</td><td>38.7</td><td>49.6</td></tr><tr><td>-</td><td>30</td><td>2.0</td><td>35.9</td><td>53.6</td></tr><tr><td>-</td><td>30</td><td>2.0</td><td>31.0</td><td>55.6</td></tr><tr><td>-</td><td>30</td><td>2.0</td><td>44.1</td><td>64.0</td></tr><tr><td>-</td><td>30</td><td>2.0</td><td>55.9</td><td>88.9</td></tr></tbody></table>

Note that *K* does not depend on the mean flow velocity *U*. Re and Fr are computed for *U* = *U*<sub>c</sub> and plotted against *K* in Figs. 6 and 7.

We may consider two kinds of estimation of the dimensionless numbers for the gel. First we may use the graphical estimation of *U*<sub>c</sub> from Fig. 4 for *d* = 1.5 cm: *U*<sub>c</sub> ~ 1.2 cm/s, leading to *v* ~ 200 cm²/s from Fig. 3 (assuming that *v* ~ 2*U*<sub>c</sub>/*d*). This allows us to plot the squares in Figs. 6 and 7. Second, we may estimate from our measurement that *U*<sub>c</sub> ~ *v*(0) is nonzero but very small (six orders of magnitude less than the other values of velocity). In that case, the squares would be far off the bottom right-hand corners of Figs. 6 and 7. In any case, the gel is possibly following the trend of Newtonian fluids with very large viscosity.

In Fig. 6, the data can be separated into two regions corresponding to partially (below the data) and completely (above the data) filled pipes. In a situation where the flow is characterized by a point belonging to the border line between the two regions, the filling would be effective over an infinite time, using an infinite amount of fluid. In the frame of a linear approximation, the data follow approximately the (−1) slope, showing that Fr(*U* = *U*<sub>c</sub>) is constant from Eq. (4).

Another way to represent the same results is given in Fig. 7. One sees that Fr is of the same order of magnitude for all experiments while *K* varies over four orders of magnitude, but Fr is not rigorously constant. In the range of diameters we investigated, we can say roughly that Fr characterizes the ability to fill the pipes
<figure><img src="image_7.png" /><figcaption>Fig. 6 Modified Reynolds number (2) computed for $U=U_c$ versus K (4) (see list of symbols in Table 1.)</figcaption></figure>

(Fr=1 corresponds to the theoretical limit (1) where the flow de-taches tangentially at the exit of a horizontal box, as suggested above in the introduction). Setting a flow rate leading to Fr≥1 would certainly fill any straight pipe.

A careful examination of Fig. 7 leads to several remarks. For Newtonian fluids, the dotted symbols, which correspond to straight pipes, lie below the points corresponding to curved pipes. The latter are gathered with no obvious dependence on the radius of curvature R. Therefore the modified Froude number, which takes into account centrifugal effects, can predict the effective filling of the curved pipes, but is of less interest in the limit of straight pipes. Each group of points corresponds to a different fluid, i.e., viscosity. The tendency for Fr(U=Uc) to increase with d is clearly visible in each group of points. For K≥10⁻², Fr(U = Uc) is apparently decreasing with K and the flow would eventually tend to follow a non Newtonian behavior.

<figure><img src="image_8.png" /><figcaption>Fig. 7 Modified Froude number (3) computed for $U=U_c$ versus K (4) (see list of symbols in Table 1.)</figcaption></figure>

In almost all our experiments, Re is larger than 10 (see Fig. 6) indicating that viscous effects are smaller than inertial ones. Meanwhile, the data on Fr indicate that both inertial and gravita-tional effects are comparable.

## Conclusion

In order to understand the failure in filling curved pipes, we performed a series of experiments using various fluids. We measured the velocity of propagation of the front v and we deduced from its linear dependence with the mean flow velocity U a criti-cal value Uc, below which there is no longer a complete filling of the pipe. This critical velocity Uc is increasing with d for each fluid.

Results are presented in a synthetic way as relations between dimensionless numbers. We have shown, in the scope of our work, that a modified Froude number determines the ability of the flow to fill a curved pipe. In all our experiments and for Newtonian fluids of different viscosities, the filling was complete for Fr ≥0.6±0.2.

We also performed experiments with a non-Newtonian fluid whose rheological behavior is similar to the one of cement sus-pensions. The front deforms very slowly when the shear is weak because this increases the viscosity of the fluid. In that case, the complete filling may be obtained for Fr very small. Following [16,17], a more extensive study would be needed in order to un-derstand the detail of the non-Newtonian effects.

For industrial applications we advise considering pipes of small diameters and injecting the solution at a relatively low flow rate in order to increase its effective viscosity. However, a future study of interest will be to repeat the present experiments in the context of slug flows.

## Acknowledgments

The authors thank C. Peiti at Ecole des Mines de Paris (CEMEF, Sophia Antipolis) for his help, and R. Leroy and O. Coussy for providing helpful comments. The authors are grateful to J.Ch. Bery for his help in the conception and realization of the experimental apparatus. This study was part of a contract with the Laboratoire Central des Ponts et Chaussées.

## Nomenclature

$d$ = diameter of pipe

$R$ = radius of curvature

$g$ = acceleration of gravity

$v$ = kinematic viscosity

$\rho$ = density

$U$ = mean entrance velocity

$u$ = mean exit velocity

$v$ = velocity of the upper point of the front

$U_c = U(v=0)$ = critical mean velocity

$h$ = height of the exit flow

$Fr = U(1+(d/R)^{1/2})/(gd)^{1/2}$ = modified Froude number

$Re = Ud(1+(d/R)^{1/2})/v$ = modified Reynolds number

$K$ = Fr/Re

$\sigma$ = surface tension

$L_c = \sigma^{1/2}(\rho g)^{-1/2}$ = capillary length

$Bo = d^2 \rho g / \sigma$ = Bond number

$We = \rho d U^2 / \sigma$ = Weber number

## References

[1] Berger, S. A., Talbot, L., and Yao, L.-S., 1983, "Flow in curved pipes," Annu. Rev. Fluid Mech., 15, pp. 461-512.

[2] Soh, W. Y., and Berger, S. A., 1984, "Laminar entrance flow in a curved pipe," J. Fluid Mech., 148, pp. 109-135.

[3] Olson, D. E., and Snyder, B., 1985, "The upstream scale of flow development in curved circular pipes," J. Fluid Mech., 150, pp. 139-158.

[4] Ishigaki, H., 1994, "Analogy between laminar flows in curved pipes and orthogonally rotating pipes," J. Fluid Mech., 268, pp. 133-145.
[5] Popiel, C. O., and van der Merwe, D. F., 1996, "Friction factor in sine-pipe flow," ASME J. Fluids Eng., **118**, pp. 341–345.

[6] Majumdar, B., Mohan, R., Singh, S. M., and Agrawal, D. P., 1998, "Experim-ental study of flow in a high aspect ratio 90 deg curved diffuser," ASME J. Fluids Eng., **120**, pp. 83–89.

[7] Chanson, H., 1997, "Air-water flows in partially-filled conduits," J. Hydraul. Res. **35**, No. 5, pp. 591–602.

[8] Capart, H., Sillen, X., and Zech, Y., 1997, "The design of sewage sludge pumping systems," J. Hydraul. Res., **35**, No. 5, pp. 659–672.

[9] Woods, B. D., and Hanratty, T. J., 1999, "Influence of Froude number on physical processes determining frequency of slugging in horizontal gas-liquid flows," Int. J. Multiphase Flow, **25**, No. 6-7, pp. 1195–1223.

[10] Taitel, Y., Sarica, C., and Brill, J. P., 2000, "Slug flow modeling for down-ward inclined pipe flow: theoretical considerations," Int. J. Multiphase Flow, **26**, No. 5, pp. 833–844.

[11] Benjamin, T. B., 1968, "Gravity currents and related phenomena," J. Fluid Mech., **31**, pp. 209–248.

[12] Simpson, J. E., 1982, "Gravity currents in the laboratory, atmosphere, and ocean," Annu. Rev. Fluid Mech., **14**, pp. 213–234.

[13] Asavenant, J., and Vanden-Broeck, J.-M., 1996, "Nonlinear free-surface flows emerging from vessels and flows under a sluice gate," J. Austral. Math. Soc., **B 38**, pp. 63–86.

[14] Vanden-Broeck, J.-M., and Keller, J. B., 1987, "Weir flows," J. Fluid Mech., **176**, pp. 283–293.

[15] Dean, W. R., 1927, "Note on the motion of fluid in a curved pipe," Philos. Mag., **20**, pp. 208–223.

[16] Chilton, R., Stainsby, R., and Thompson, S., 1996, "The design of sewage sludge pumping systems," J. Hydraul. Res., **34**, No. 3, pp. 395–408.

[17] Chilton, R., and Stainsby, R., 1998, "Pressure loss equations for laminar and turbulent non-Newtonian pipe flow," J. Hydraul. Eng., **124**, No. 5, pp. 522–529.