5 ENERGY EQUATION OF FLUID MOTION 5.1 Eulerian Approach & Control Volume In order to develop the equations that describe a flow, it is assumed that fluids are subject to certain fundamental laws of physics. The pertinent laws are: 1. Conservation of mass;. Conservation of energy; 3. Conservation of momentum. These principles were initially developed for the case of a solid body and the application of these laws to a solid body is relatively straightforward since the body will be of measurable size and mass. However, it is not for a flowing fluid. To describe the displacement, velocity and acceleration of a fluid, control volume approach (Eulerian approach) is commonly used. A control volume is a purely imaginary region within a body of flowing fluid. The region is usually at a fixed location and of fixed size. Inside the region, all of the dynamic forces cancel each other. Attention may therefore be focused on the forces acting externally on the control volume. The control volume may be of any shape. Therefore, a shape may be selected which is most convenient for any particular application. 5. Continuity Equation The principle of conservation of mass can be applied to a flowing fluid. Considering any fixed region in the flow constituting a control volume, Mass of fluid entering per unit time = Mass of fluid leaving per unit time + Increase of mass of fluid in the control volume per unit time P.5-1
mass of fluid entering Control Volume mass of fluid leaving Continuity of Flow For steady flow, the mass of fluid in the control volume remains constant and the relation reduces to Mass of fluid entering per unit time = Mass of fluid leaving per unit time Apply this principle to steady flow in a streamtube as shown below, Area = da1 Velocity = V1 Pressure = p1 1 Area = d A Velocity = V Pressure = p Continuous flow through a streamtube If there is no flow being evacuated from the stream tube except at the outlet section, then Q = ρv 1 A 1 = ρv A Q = V 1 A 1 = V A where V 1, V are the mean velocities at sections 1 and respectively. The continuity equation for liquids many be generally expressed in the form Q = V*A = constant P.5-
where V and A are velocity of flow and the cross-sectional area normal to the flow respectively at any section of the liquid stream. The continuity equation can also be applied to determine the relation between the flows into and out of a junction. For steady condition, Total inflow to junction = total outflow from junction A V Q A1 V1 Q1 A3 V3 Q3 ρ 1 Q 1 = ρ Q + ρ 3 Q 3 For an incompressible flow, ρ 1 = ρ = ρ 3, then Q 1 = Q + Q 3 or A 1 V 1 = A V + A 3 V 3 In general, if we consider flow towards the junction as positive and flow away for the junction as negative, then for steady flow at any junction the algebraic sum of all the mass flows must be zero: Σ ρq = 0 P.5-3
Worked example: Water flows through a pipe AB of diameter d 1 = 50 mm, which is in series with a pipe BC of diameter d = 75 mm in which the mean velocity V = m/s. At C the pipe forks and one branch CD is of diameter d 3 such that the mean velocity V 3 is 1.5 m/s. The other branch CE is of diameter d 4 = 30 mm and conditions are such that the discharges Q from BC divides so that Q 4 = ½ Q 3. Calculate the values of Q 1, V 1, Q, Q 3, d 3, Q 4 and V 4. Answer A B C d1=50 mm d=75 mm V= m/s V3=1.5 m/s d4=30 mm D Q3=Q4 Since pipes AB and BC are in series and water is incompressible, the volume flow rate will be the same in each pipe, Q 1 = Q. But E Q = A * V = (π/4)*d *V Q 1 = Q = (π/4)*(0.075) * = 8.836*10-3 m 3 /s Q1 V 1 = Q 1 /A 1 = ( π / 4) *d 8.836 *10 ( π / 4) * 0.05 = 4.5 m/s 3 = Considering pipes BC, CD and DE, the discharge from BC must be equal to the sum of the discharges through CE and CE. Therefore 1 And since Q = Q 3 + Q 4 Q 4 = ½ Q 3 P.5-4
Q = 1.5 Q 3 Hence Q 3 = Q / 1.5 = 8.836*10-3 /1.5 = 5.891*10-3 m 3 /s & Q 4 = ½ Q 3 =.945*10-3 m 3 /s Also, since Q 3 = (π/4)*d 3 *v 3 4Q3 Hence d 3 = πv 3 4*5.891*10 = π*1.5 = 0.071 m Q4 V 4 = ( π / 4) *d 4 3.945*10 ( π / 4) * 0.03 = 4.17 m/s = 3 P.5-5
5.3 Bernoulli s Equation Figure below shows a small streamtube element of cylindrical section with uniform cross-sectional area. The fluid accelerates in the direction of the flow. Taking this to be a control volume at an instant in the time, ds p' (p+dp).a z p.a α mg p' Forces acting on the element in the s direction include: (i) Pressure forces expressed as: upstream end p. A downstream end -(p + p s s) A (pressure may vary with respect to space, s and time, t, therefore the partial derivative is used). (ii) Gravity force W.sinα where W = γ. s. A (iii) sinα = z s, then as s approaches zero, then sinα = z s By applying Newton s second law in the s direction, the following equation is reached. ΣF s = M*a s where M = mass of the cylindrical element (ρ. s. A) a s = acceleration in the s direction P.5-6
Therefore, we arrive at the following equation by substituting the above expression into the above equation. γ ρ s ( p + z) =.a s - Euler s equation By rewriting the acceleration term in a general form as follows: a s = = dv s dt Vs t V + s s ds dt V s For steady flow, = 0 t Then a s = V s s V s By substituting the above expression into Euler s equation we obtain the following (p + γz) = ρ.v s s V s s Assuming the density of the fluid remains unchanged, z ρ + γ + (Vs s s s p ) = 0 By integrating the above equation with respect to s we have, p + γz + ρ V s = constant or p γ + z + V g = constant The subscript s is usually omitted and v is used to represent the velocity in the direction of fluid flow. P.5-7
The term p/γ - pressure head z - potential head (elevation head) V /g - velocity head The constant at the right hand side of the equation is the total head (total energy per unit weight) of the flow field. The equation is applicable when the flow is steady, non-viscous, and incompressible with the constant density. Under these conditions the total head at any point along a streamline of the flow field is the same. Bernoulli s equation is an expression of the principle of conservation of energy. In most applications, Bernoulli s equation is applied at two points in the flow field with z 1 and z are referred to the same datum. p1 γ + z 1 + V1 g = p γ + z + V g P.5-8
Worked examples: 1. Determine the velocity of the water flowing at point 1 and of the pipe shown below, with the following data. Neglect all losses in heads. p 1 = 10 kpa, p = 00 kpa, d 1 = 0.4 m and d = 0.6 m 1 Answer Apply Continuity equation to point 1 and V 1 A 1 = V A π(0.4) Since A 1 = m, A = 4 (0.6) V 1 = V * (0.4) =.5 V π(0.6) 4 m Apply Bernoulli s equation along the central streamline to points 1 and p1 V1 p V + + z1 = + + z γ g γ g p 1 = 10 kpa, p = 00 kpa V 1 =.5 V & z 1 = z (same level) Hence 10 (.5V ) + 9.81 *9.81 V = 6.3 m/s 00 = 9.81 + V *9.81 and V 1 =.5 v = 14. m/s P.5-9
. For a 50 mm diameter siphon drawing oil (S.G. = 0.8) from the oil reservoir as shown below, the head loss from point 1 to point is 1.5 m and from point to point 3 is.40 m. Find the discharge of oil from the siphon and the oil pressure at point. 1 oil m 5m dia = 50 mm 3 Answer Apply Bernoulli s equation to points 1 & 3. p γ 1 V1 p3 V3 + + z1 = + + z3 + loss1-3 g γ g p 1 = p 3 = 0 (p atm = 0) V 1 = 0 (Assume the oil vessel is very large, the rate of the oil level drop is negligible) z 3 = 0, z 1 = 5 m loss 1-3 = loss 1- + loss -3 = 1.5 +.4 m = 3.9 m 0 + 0 + 5 = 0 + V 3 /(*9.81) + 0 + 3.9 V 3 = 4.646 m/s Q = A 3 * V 3 = [π(0.05) /4]*4.646 = 0.0091 m 3 /s P.5-10
Again, apply Bernoulli s equation to points 1 & p1 V1 p V + + z1 = + + z + loss 1- γ g γ g z = 7 m V = V 3 (same pipe diameter) = 4.646 m loss 1- = 1.5 m 0 + 0 + 5 = p /γ + (4.646) /(*9.81) + 7 + 1.5 p /γ = - 4.6 m p = - 4.6 * 9.81 * 0.8 kpa = - 36.9 kpa P.5-11
5.4 Flow Measurement The three parameters most often involved in flow measurement are pressure, velocity and discharge of flow. The relationships between these parameters have already been defined through the Continuity and Bernoulli s equation. In this section, the principles and applications of common flow measurement devices are discussed. 5.4.1 Stagnation Tube and Pitot Tube When an L shaped tube is placed in a flow field, it can be used to measure the velocity of flow. h va,pa A z vb,pb B flow (a) (b) Since the flow is brought to a standstill at the front end of the tube, this simple device is called a stagnation tube (fig. a). The velocity of flow measured by the stagnation tube is expressed as V = gh - Torricelli s formula where V = velocity of flow h = the level of fluid in the tube above the free surface. This expression can be derived by applying Bernoulli s equation to the stagnation point and any other point upstream of it. For simplicity, choose a point (point A) at the same level as the stagnation point (point B, i.e. the velocity of fluid flow at this point is zero). P.5-1
From Bernoulli s equation, p γ A + A V g + z A = p γ B + B V g + z B Since z A = z B, and V B = 0 (stagnation point), therefore VA g = p γ B - p γ A which can be applied to flow in closed or open conduits. Since also p A = γz and p B = γ(z + h), then V A = gh, where tube V A = the velocity of flow which is measured by the stagnation Since the application of the stagnation tube is limited by the pressure of the flow, pitot tubes are often used in pressure pipes. P.5-13
Worked examples: 1. Water flows through the pipe contraction shown in figure below. For the given 0. m difference in manometer level, determine the flow rate when the diameter of the small pipe, D is 0.05m. 0. m 0.1 m (1) h1 water Q h () D Answer Again, apply Bernoulli s equation to points 1 & p1 V1 p V + + z1 = + + z γ g γ g V 1 = 0 (stagnation point) z 1 = z = 0 (same level) p1 p V + 0 + 0 = + + 0 γ γ g V = (p g p ) γ 1 but p 1 = γh 1 and p = γh p 1 p = γ(h 1 h ) = 0. γ Thus V = 0.γ g * = g *0. m/s γ = 1.98 m/s Q = A * V = [π(0.05) /4]*1.98 = 3.888*10-3 m 3 /s = 3.888 L/s P.5-14
. Water is siphoned from a tank as shown below. Determine the flow rate and the pressure at point A, a stagnation point. 0.04 m diameter 1 Answer A 3m water Apply Bernoulli s equation to points 1 & p1 V1 p V + + z1 = + + z γ g γ g p 1 = p = 0 (p atm ) V 1 = 0 z 1 = 3m, z = 0 0 + 0 + 3 = 0 + V /(*9.81) + 0 V = *9.81* 3 m/s = 7.67 m/s Hence Q = A * V = [π(0.04) /4]*7.67 m 3 /s = 9.64 * 10-3 m 3 /s Again, apply Bernoulli s equation to points 1 & A p1 V1 pa VA + + z1 = + + z A γ g γ g V A = 0 & z A = 0 0 + 0 + 3 = p A /γ + 0 + 0 or p A = γ* 3 = 9.81 * 3 kpa = 9.43 kpa P.5-15
5.4. Orifice An orifice is a geometric opening in the side of a wall and tank, through which fluid can flow. A circular sharp edged orifice at the side of a water tank is shown. A va H B vb Orifice area, A The volume rate of flow discharged through an orifice will depend upon the head of the fluid above the level of the orifice and it can therefore be used as a means of flow measurement. Applying Bernoulli s equation to A and B, assuming that there is no loss of energy, pa VA pb VB + + z A = + + z B γ g γ g Putting z A - z B = H, V A = 0, V B = V and p A = p B, velocity of jet, V = gh Theoretically, if A is the cross-sectional area of the orifice, Discharge, Q = Area * Velocity = A gh In practice, the actual discharge is considerably less than the theoretical discharge and is modified by introducing a coefficient of discharge, C d, so that Actual discharge, Q actual = C d * Q theoretical = C d A gh P.5-16
There are two reasons for this phenomenon. First, the velocity of the jet is less than that the theoretical velocity because there is a loss of energy between A and B: Actual velocity at B = C v * V = C v gh where C v = coefficient of velocity 0.97 Second, the paths of the particles of the fluid converge on the orifice and the area of the issuing jet at B is less than the area of the orifice A at C. In the plane of the orifice, the particles have a component of velocity towards the centre and the pressure at C is greater than atmospheric pressure. The streamlines will converge at a short distance downstream of the orifice and the minimum flow area is called the vena contracta. This is smaller than the cross-sectional area of the orifice. vena contracta C v B vb Contraction of issuing jet The ratio of the cross-sectional area of the vena contracta to that of the orifice is called the coefficient of contraction, C c. Actual area of the jet at B = C c A Hence Actual Discharge = Actual area at B * Actual velocity at B = C c A * C v gh = C c * C v A gh C d = C c * C v P.5-17
Evaluation of the actual velocity from the displacement of the jet vena contracta centre line of jet y x For horizontal motion, S = x, u = V, and a = 0 Hence x = V * t where V = actual velocity at vena contracta For vertical motion, u = 0, y = a = g 1 g * t By eliminating t, V = gx y P.5-18
Worked examples: 1. Oil of specific gravity 0.8 discharges from an open tank through an orifice of diameter 14 mm. The coefficient of velocity is 0.88 and the coefficient of contraction is 0.6. The centre of the orifice is at a depth of 0.9m from the surface of the oil. Determine the diameter of the vena contracta and the discharge of oil through the orifice. Answer 0.9m orifice hole dia = 14mm A = π * (14/1000) /4 = 1.54 * 10-4 m Since C c = A c /A 0.6 = A c / 1.54 * 10-4 hence A c = 0.955 * 10-4 m or d c = 0.0110 m = 11.0 mm Theoretical velocity, V = gh = * 981. * 09. = 4.00 m/s Theoretical discharge, Q = V*A = 4.0 * 1.54 * 10-4 = 6.47 * 10-4 m 3 /s Coefficient of discharge, C d = C v * C c = 0.88 * 0.6 = 0.546 Actual discharge = C d * Theoretical discharge = 0.546 * 6.47 * 10-4 = 3.53 * 10-4 m 3 /s = 0.353 L/s P.5-19
. If viscous effects are neglected and the tank is large, determine the flow rate from the tank as shown below. (Neglect all losses) orifice dia = 50 mm 3 1 oil m S.G.=0.81 0.7m water Answer Apply Bernoulli s equation between to & 3 p V p3 V3 + + z = + + z3 γ g γ g p = p 1 + γ oil h = γ oil (p 1 = p atm = 0) p 3 = 0 (p atm ) V = 0 z = 0.7 m & z 3 = 0 Thus γ o h/γ w + 0 + 0.7 = 0 + V 3 /(g) + 0 γ oil V 3 = g *( + 0.7) m/s γ w Since S.G. oil = γ oil /γ w = 0.81 V 3 = g * ( * 0.81+ 0.7) m/s = 6.747 m/s Q = A 3 * V 3 = [π(0.05) /4]*6.747 m 3 /s = 0.013 m 3 /s P.5-0
5.4.3 Venturi Meter The venturi meter is the most commonly used device for flow measurement and is designed with a streamlined throat which can, 1. reduce the overall head loss, and. reduce the normal mechanical wear of the device. The discharge equation for a venturi meter is the same as for an orifice meter, except that the flow coefficient value of a venturi meter is generally higher than that of an orifice meter and a flow nozzle. This follows since the coefficient of contraction of a venturi meter is unity. converging cone entry throat diverging section 1 h Applying Bernoulli s equation with p 1 = ρgh 1, p = ρgh and z 1 = z, and also by continuity equation, V 1 A 1 = V A V 1 = gh ( h) A 1 ( ) A 1 1 Hence discharge Q = A 1 V 1 A1 = A1 ( ) A 1 gh ( h) 1 Actual discharge = C d Q CA = A1 ( ) A d 1 1 gh ( h) 1 P.5-1
Worked examples: 1. The flow of kerosene is measure using a venturi meter. The diameter of the pipe and the throat are 50 mm and 5 mm respectively. A differential manometer shows a deflection, h of 55 mm of mercury. The coefficient of discharge is 0.96. Determine the flow rate of kerosene. Take the S.G. of kerosene as 0.8. Answer 1 y ρ A h' B ρ m Area of cross-section of the pipe, A 1 = π * 0.05 /4 = 0.00196 m Area of cross-section of the throat, A = π * 0.05 /4 = 0.00049 m By continuity equation, A 1 0.00196 V = ( )*V1 = ( )*V 1 = 4 V 1 A 0.00049 By Bernoulli s equation, z 1 = z p 1 = p + h *(γ m - γ k ) (p 1 - p )/ γ k S. G. m = h *( S. G. k - 1) 13.6 = 0.055*( -1) m 0.8 = 0.857 m p1 V1 p V + + z1 = + + z γ g γ g V1 (4V1) 0.857 + = g g V 1 = 1.059 m/s or V = 4*1.059 = 4.35 m/s Q actual = 0.96 * 4.35 * 0.00049 = 0.00199 m 3 /s = 1.99 L/s P.5-
. The water supply to a gas water heater contracts from 10 mm in diameter at 1 to 7 mm in diameter at (figure below). If the pipe is horizontal, calculate the difference in pressure between 1 and when the velocity of the water at 1 is 4.5 m/s. The pressure difference operates the gas control through connections which are taken to horizontal cylinder in which a piston of 0 mm diameter moves. Ignoring friction and the area of the piston connecting rod, what is the force on the piston? 1 Answer For Continuity of flow, V 1 A 1 = V A Since π(0.1) A 1 = 4 m, A = V 1 = 4.5 m/s (0.1) V = 4.5 * (0.07) m/s = 9.18 m/s Apply Bernoulli s equation to points 1 & p1 V1 p V + + z1 = + + z γ g γ g z 1 = z (same level) V 1 = 4.5 m/s V = 9.18 m/s π(0.07) 4 m p 4.5 p 9.18 hence 1 + + 0 = + 0 γ *9.81 γ *9.81 + P.5-3
p 1 p γ p 1 p 9.18 4.5 = *9.81 = 3.67 m of water = 3.67 * 9.81 kpa = 3.045 kpa π(0.0) Area of piston, A p = m 4 = 0.000314 m Force on piston = (p 1 p )*A p = 3.045 * 0.000314 kn = 0.01 kn P.5-4
3. Air flows through the device shown below. If the flow rate is large enough, the pressure within the construction will be low enough to draw the water up into the tube. Determine the flow rate, Q and the pressure needed at section 1 to draw the water into section. Neglect compressibility and viscous effects. dia = 5 mm 1 air 3 dia = 50 mm 0.3m dia = 50 mm free jet water Answer For Continuity of flow, V A = V 3 A 3 (d3) V = V 3 * (d ) (50) = V 3 * (5) = 4V 3 Apply Bernoulli s equation between to & 3 p V p3 V3 + + z = + + z3 γ g γ g p 3 = 0 z = z 3 (same level) p = -γ w h p γ w = - h γ γ 9.81*1000 = - *0. 3 1 = - 45.5 m m P.5-5
Thus - 45.5 + (4V 3 ) /(g) + 0 = 0 + V 3 /(g) + 0 (4 1) V 3 = 45.5 * * 9.81 = 17.91 m/s V 3 Q = A 3 *V 3 = π(0.05) /4*17.91 m 3 /s = 0.0351 m 3 /s Also apply Bernoulli s equation between to 1 & 3 p1 V1 p3 V3 + + z1 = + + z3 γ g γ g V 1 = V 3 (same diameter) z 1 = z 3 (same level) therefore p 1 = p 3 = 0 (P atm ) P.5-6
Class Exercise 5.1: Water flows down the ramp in the channel as shown below. The channel width decreases from 4.5 m at section 1 to.5 m at section. For the conditions shown, determine the flow rate. 1 Q m 0.5m width = 4.5 m 1m width =.5 m (14.1 m 3 /s) P.5-7
Class Exercise 5.: The water in a tank is 1.8 m deep and over the surface is air at pressure 70 kpa (gauge). Find the flow rate from an orifice of 50 mm in the bottom of the tank if the C d = 0.6 pressure 70 kpa 1 water 1.8 m orifice dia = 50 mm (0.0156 m 3 /s) P.5-8
Class Exercise 5.3: A horizontal venturi tube, 80 mm diameter at the entrance and 140 mm diameter at the throat, has a discharge coefficient of 0.97. A U-tube manometer filled with mercury is connected between the entrance and the throat to measure the water flow between them. Calculate the flowrate when the difference in the mercury level is 50 mm. (0.054 m 3 /s) P.5-9
Tutorial: Energy Equation 1. Water flows from the faucet on the first floor of the building shown below with a maximum velocity of 6.5 m/s. For steady inviscid flow, determine the maximum water velocity from the basement faucet and from the faucet on the second floor. Assume only two taps will be opened at each time and each floor is 4 m tall. 1.5m v=6.5m/s 4m.5m 1.5m 1.0m 1.5m inlet. For the water shooting out of the pipe and nozzle under the conditions shown in figure below, find the height, h about the nozzle to which the water jet will shoot. Assume negligible head loss. water h 1.1m nozzle dia=0.1m A point A dia=0.m PA=55kPa 3 Oil (S.G = 0.86) flows through the system as shown in figure below. A manometer is used to measure the flowrate of the oil. By neglecting any losses, determine the flowrate of the oil. A=0.01m 1.5m oil 0.5m Q A=0.030m SG=.5 P.5-30