by farhan ahmad [email protected]/uploads/1/4/2/7/... ·...

47
Incompressible Flow in Pipes and Channels By Farhan Ahmad [email protected] Department of Chemical Engineering, University of Engineering & Technology Lahore

Upload: others

Post on 05-Aug-2020

3 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

Incompressible Flow in

Pipes and Channels

By

Farhan Ahmad

[email protected]

Department of Chemical Engineering,

University of Engineering & Technology Lahore

Page 2: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

2

Industrial processes - flow of fluids through pipes, conduits,

and processing equipment.

Circular cross-section

Non-circular cross-section

Flow of fluids in

Totally or partially filled pipes,

Layers down vertically inclined surfaces,

Through beds of solids, and

Agitated vessels.

Significance

Page 3: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

3

Page 4: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,
Page 5: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

5

Page 6: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

6

Page 7: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

7

Consider the steady flow of a viscous fluid at constant

density in fully developed flow through a horizontal tube.

Visualize a disk-shaped element of fluid, concentric with the

axis of the tube.

Flow of Incompressible Fluids in Pipe –

Shear-Stress Distribution

Page 8: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

8

Flow of Incompressible Fluids in Pipe –

Shear-Stress Distribution

Page 9: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

9

At wall

After subtraction

Relation between τ and r

At r =0 , τ = 0

Page 10: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

10

Pressure Drop

Apply the balance

Relation between Skin Friction and Wall Shear

Page 11: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

11

ratio of the wall shear stress to the product of the density and

the velocity head.

Friction Factor

Page 12: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

12

Relations between Skin Friction Parameter

Page 13: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

13

Laminar

Turbulent

Fluid may be

• Newtonian

• Non-Newtonian

Flow in Pipe

Page 14: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

14

Velocity Distribution

Average velocity

Momentum and Kinetic energy correction factors

Laminar Flow of Newtonian Fluids

Page 15: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

15

Circular cross-section

Local velocity u depends on radius r

Consider a thin ring of radius r and width dr

According to Newton's law of viscosity

Laminar Flow of Newtonian Fluids –

Velocity Distribution

Page 16: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

Using both equations

Integrate with boundary condition u = 0, at r = rw

16

Page 17: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

17

Maximum velocity is at the center of pipe i.e., at r = 0

Relation of local to maximum velocity

Maximum velocity

Page 18: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

18

Graphical representation

Page 19: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

19

Average Velocity

𝑺 = 𝝅𝒓𝒘𝟐

Page 20: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

20

Kinetic energy correction factor

For Laminar Flow α = 2

Page 21: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

21

Momentum correction factor

For Laminar Flow β = 4/3

Page 22: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

22

Hagen-Poiseuille Equation

Page 23: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

23

Velocity variations with radius for power law fluids

The pressure difference for power law fluids

Laminar Flow for Non-Newtonian Liquids

Page 24: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

24

Laminar Flow for Non-Newtonian Liquids

Page 25: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

25

Bingham-plastic fluids:

The general shape of the curve of u versus r in case of Bingham-plastic fluids is;

In the central portion - no velocity variation with the radius

the velocity gradient is confined to an annular space between the central portion and tube wall.

The center portion is moving in plug flow.

Laminar Flow for Non-Newtonian Liquids

Page 26: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

26

The velocity distribution is;

The shear diagram is;

Page 27: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

27

Laminar Flow for Non-Newtonian Liquids

Bingham-plastic fluids:

For the velocity variation in the annular space between the tube wall and the plug, the following equation applies;

and

Page 28: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

28

Turbulent Flow in Pipes and Closed Channels

Viscous Sublayer

Buffer layer

Transition core

Page 29: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

29

Velocity Distribution for Turbulent Flow

Newtonian fluid

Turbulent flow at Reynolds No 10000

Smooth pipe

Velocity gradient is zero at centerline

Turbulent core – eddies – large but of low intensity

Transition zone – eddies – small but intense

Kinetic energy

At centerline - isentropic turbulence – anisotropic in turbulence core

Page 30: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

30

Velocity Distribution for Turbulent Flow

It is customary to express the velocity distribution in turbulent flow not as velocity vs. distance but in terms of dimensionless parameters defined by the following eqns;

Page 31: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

31

Velocity Distribution for Turbulent Flow

For the velocity distribution in the laminar sublayer;

An empirical equation for the so-called buffer layer is;

An equation proposed by Prandtl for the turbulent core is;

Page 32: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

32

Velocity Distribution for Turbulent Flow

Page 33: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

33

Flow Quantities for Turbulent Flow

Average Velocity:

Page 34: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

34

Flow Quantities for Turbulent Flow

Reynolds Number – Friction Factor Law for Smooth Pipe:

Von Karman equation

Page 35: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

35

Flow Quantities for Turbulent Flow

Kinetic Energy and Momentum Correction Factors:

For turbulent flow f is of the order of 0.004, and for this value

and both are assumed to be unity in case of turbulent flow.

Page 36: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

36

Flow Quantities for Turbulent Flow

Relation between Maximum velocity and Average Velocity:

Page 37: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

37

Flow Quantities for Turbulent Flow

Effect of Roughness:

In turbulent flow, a rough pipe leads to a larger friction factor for

a given Reynolds number than a smooth pipe does.

If a rough pipe is smoothed, the friction factor is reduced.

When further smoothing brings about no further reduction in

friction factor for a given Reynolds number, the tube is said to be

hydraulically smooth.

Page 38: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

38

Flow Quantities for Turbulent Flow

Effect of Roughness:

Roughness parameter k

f is a function of both NRe and the relative roughness k/D, where D is the diameter of the pipe.

Page 39: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

39

Flow Quantities for Turbulent Flow

Effect of Roughness:

All clean, new commercial pipes seem to have the same type of

roughness.

Each material of construction has its own characteristic

roughness parameter.

Old, foul and corroded pipe can be very rough, and the character

of the roughness differs from that of clean pipe.

Roughness has no appreciable effect on the friction factor for

laminar flow unless k is so large that the measurement of the

diameter becomes uncertain.

Page 40: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

40

Friction Factor Chart

Page 41: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

41

Friction Factor Chart

For Laminar flow – straight line with slope -1

For turbulent flow the lowest line represents the friction factor

for smooth tubes. A much more convenient empirical equation

for this line is the relation;

Over a range of Reynolds number from about 50,000 to 1 × 106

Over a range of Reynolds number from about 3000 to 3 × 106

Page 42: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

42

Friction Factor Chart

Comparing the above two equations

For Power Law Fluids

Page 43: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

43

Friction Factor Chart

Page 44: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

44

Drag Reduction in Turbulent Flow

Page 45: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

45

Effect of Heat transfer / Non-isothermal flow

When the fluid is either heated or cooled by a conduit wall hotter or colder than the fluid, the velocity gradient is changed.

The effect on the velocity gradients is especially pronounced with liquids where viscosity is a strong function of temperature.

1. The Reynolds number is calculated on the assumption that the fluid temperature equals the mean bulk temperature, which is defined as the arithmetic average of the inlet and outlet temperatures.

2. The friction factor corresponding to the mean bulk temperature is divided by a factor

Page 46: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

46

Effect of Heat transfer / Non-isothermal flow

<

Page 47: By Farhan Ahmad farhanahmad@uet.educhemicalengineeringdocuments.weebly.com/uploads/1/4/2/7/... · 2018-09-07 · Farhan Ahmad farhanahmad@uet.edu.pk Department of Chemical Engineering,

47

Flow through Channels of Non-Circular cross-sections

Equivalent Diameter:

It is four times the hydraulic radius.

Hydraulic Radius:

It is the ratio of the cross-sectional area of the channel to the wetted perimeter of the channel.