Search
 
 

FLUID MECHANICS - THEORY

   

This section begins with a brief review on laminar and turbulent flows followed by an introduction of fully developed flow. The Hagan-Poiseuille and Darcy-Weisbach equations are then presented for circular pipe flow. In this section, the pipe is assumed to be fully occupied by the fluid and the pressure in the fluid is greater than zero gauge. Flows in partially filled pipes are covered in the discussion of open channel flow. At the end of this section, flow in non-circular conduits is discussed.

     
    Entrance and Fully Developed Regions


Entrance and Fully Developed Regions

 

Recall that Reynolds number (Re = ρVD/μ) can be used to classify if a flow is laminar or turbulent. The type of flow not only depends on its velocity (V), but also its density (ρ), viscosity (μ) and length scale (diameter D). For pipe flow, the flow is laminar when Re < 2,100, and it is turbulent when Re > 4,000. The flow is referred to as transition flow when Reynolds number is in between 2,100 and 4,000.

When a fluid enters a pipe, viscous effects due to the pipe wall will develop. The region where viscous effects are important is referred to as the boundary layer. The velocity profile, as shown in the figure, will also vary due to the growth of this boundary layer. When the velocity profile reaches a constant (i.e., velocity profile no longer changes along the pipe), the flow is said to be fully developed. The length required for the flow to reach fully developed conditions is called the entrance length (Le), and it can be determined from the following empirical relations:

     For laminar flow:      Le/D = 0.06Re

     For turbulent flow:     Le/D = 4.4(Re)1/6

     
    Hagan-Poiseuille Flow


Parabolic Velocity Profile

 

Fully developed laminar pipe flow is generally known as Hagan-Poiseuille flow. The velocity profile of a fully developed laminar horizontal circular pipe flow can be derived from the Navier-Stokes equations and is given by

     

The derivation process is similar to that of obtaining the velocity distributions for flow between parallel plates, hence the derivation details are omitted here for brevity. The assumptions used in the derivation are as follows:

     (a) The fluid is incompressible (Dρ/Dt = 0).
     (b) The flow is steady ().
     (c) The flow is in the z-direction only (vr = vθ = 0).
     (d) The flow is axisymmetric ().

The obtained velocity profile is parabolic, and it closely resembles the velocity profile obtained for flow between parallel plates.

The volumetric flow rate (Q) can then be determined by integrating the velocity profile over the cross-sectional area to yield:

     

where Δp is the pressure drop across a pipe length of L and radius R. The above equation is known as Hagan-Poiseuille's law.

     
    Darcy-Weisbach Equation

 

 

From another discussion, it is known that Bernoulli's equation can be applied along a streamline to analyze steady, incompressible and inviscid flow. Most flows in real practice, however, are viscous, and hence losses due to viscous dissipation (friction) are unavoidable. To analyze flow in a viscous pipe, Bernoulli's equation can be modified and applied to two sections of the pipe flow:

     

where p is the static pressure, V is the average velocity, z is the elevation, and hL accounts for the head loss due to friction between section (1) and (2).

     


Moody Chart

 

Table: Typical Values of
Surface Roughness
Pipe (New) Surface
Roughness, ε
ft (10-3) mm
Riveted Steel 3 - 30 0.9-9.0
Concrete 1 - 10 0.3-3.0
Wood Stave 0.6 - 3 0.18-0.9
Cast Iron 0.85 0.26
Galvanized Iron 0.5 0.15
Commercial Steel
or Wrought Iron
0.15 0.046
Drawn Tubing 0.005 0.0015
PVC/Plastic
(smooth)
0.005-
0.023
0.0015-
0.007
 

Based on dimensional analysis, the head loss over the length of the pipe is given by

     

where f is an empirical friction factor. The above expression is also called the Darcy-Weisbach equation. The friction factor is a function of the Reynolds number (Re) and relative roughness of the pipe (ε/D) as shown in the Moody chart to the left. Typical values of the surface roughness of pipe (ε) are summarized in the table.

The Moody Chart actually has three main regions or sections. The first region is for laminar flow for when Re is approximately less than 2100. Then the friction factor is a a direct function of the Re, or

     For laminar flow:     f = 64/Re

The region between laminar and completely turbulent. This is when the chart is required. First determine the relative roughness number, ε/D, and then determine Re. Using these two numbers and the Moody Chart, estimate the friction factor, f.

In general, the friction factor is a function of Re and ε/D for turbulent flow. However, when Re is very high (completely turbulent flow), the friction factor becomes independent of Re, and it depends on the relative roughness only. This is the third region and as an alternative to determining the friction factor for turbulent flow from the Moody chart, it can be calculated using the correlation provided by Colebrook**,

     For turbulent flow:

          

Note that iteration is needed to obtain f since the above equation is not expressed explicitly in f.

[**Reference: Colebrook, C. F., "Turbulent Flow in Pipes with Particular Reference to the Transition between the Smooth and Rough Pipe Laws," Journal of the Institute of Civil Engineering London, Vol. 11, pp. 133-156, 1939.]

     


Hydraulic Diameter (Dh = 4A/P)

  Flow in Non-Circular Conduits

 

To analyze flow in conduits with non-circular shapes, the hydraulic diameter (Dh) can be used in place of the diameter. It is defined as

     Dh = 4A/P

where A is the cross-sectional area, and P is the wetted perimeter.

     
   
 
Practice Homework and Test problems now available in the 'Eng Fluids' mobile app
Includes over 250 free problems with complete detailed solutions.
Available at the Google Play Store and Apple App Store.