Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: FLUID MECHANICS TUTORIAL No. 3 BOUNDARY LAYER THEORY In order to complete this tutorial you should already have completed tutorial 1 and 2 in this series. This tutorial examines boundary layer theory in some depth. When you have completed this tutorial, you should be able to do the following. Discuss the drag on bluff objects including long cylinders and spheres. Explain skin drag and form drag. Discuss the formation of wakes. Explain the concept of momentum thickness and displacement thickness. Solve problems involving laminar and turbulent boundary layers. Throughout there are worked examples, assignments and typical exam questions. You should complete each assignment in order so that you progress from one level of knowledge to another. Let us start by examining how drag is created on objects. 1 1. DRAG When a fluid flows around the outside of a body, it produces a force that tends to drag the body in the direction of the flow. The drag acting on a moving object such as a ship or an aeroplane must be overcome by the propulsion system. Drag takes two forms, skin friction drag and form drag. 1.1 SKIN FRICTION DRAG Skin friction drag is due to the viscous shearing that takes place between the surface and the layer of fluid immediately above it. This occurs on surfaces of objects that are long in the direction of flow compared to their height. Such bodies are called STREAMLINED. When a fluid flows over a solid surface, the layer next to the surface may become attached to it (it wets the surface). This is called the ‘no slip condition’. The layers of fluid above the surface are moving so there must be shearing taking place between the layers of the fluid. The shear stress acting between the wall and the first moving layer next to it is called the wall shear stress and denoted τw. The result is that the velocity of the fluid u increases with height y. The boundary layer thickness δ is taken as the distance required for the velocity to reach 99% of uo. This layer is called the BOUNDARY LAYER and δ is the boundary layer thickness. Fig. 1.1 Shows how the velocity "u" varies with height "y" for a typical boundary layer. Fig.1.1 In a pipe, this is the only form of drag and it results in a pressure and energy lost along the length. A thin flat plate is an example of a streamlined object. Consider a stream of fluid flowing with a uniform velocity uo. When the stream is interrupted by the plate (fig. 1.2), the boundary layer forms on both sides. The diagram shows what happens on one side only. Fig. 1.2 The boundary layer thickness δ grows with distance from the leading edge. At some distance from the leading edge, it reaches a constant thickness. It is then called a FULLY DEVELOPED BOUNDARY LAYER. 2 ρu o x µ x is the distance from the leading edge. At low Reynolds numbers, the boundary layer may be laminar throughout the entire thickness. At higher Reynolds numbers, it is turbulent. This means that at some distance from the leading edge the flow within the boundary layer becomes turbulent. A turbulent boundary layer is very unsteady and the streamlines do not remain parallel. The boundary layer shape represents an average of the velocity at any height. There is a region between the laminar and turbulent section where transition takes place The Reynolds number for these cases is defined as: (R e ) x = The turbulent boundary layer exists on top of a thin laminar layer called the LAMINAR SUB LAYER. The velocity gradient within this layer is linear as shown. A deeper analysis would reveal that for long surfaces, the boundary layer is turbulent over most of the length. Many equations have been developed to describe the shape of the laminar and turbulent boundary layers and these may be used to estimate the skin friction drag. Note that for this ideal example, it is assumed that the velocity is the undisturbed velocity uo everywhere outside the boundary layer and that there is no acceleration and hence no change in the static pressure acting on the surface. There is hence no drag force due to pressure changes. CALCULATING SKIN DRAG The skin drag is due to the wall shear stress τw and this acts on the surface area (wetted area). The drag force is hence: R = τw x wetted area. The dynamic pressure is the pressure resulting from the conversion of the kinetic energy of the stream into pressure and is defined by the 2 ρu o expression .The drag coefficient is defined as 2 Drag force C Df = dynamic pressure x wetted area 2τ 2R C Df = 2 =w 2 ρu 0 x wetted area ρu 0 Note that this is the same definition for the pipe friction coefficient Cf and it is in fact the same thing. It is used in the Darcy formula to calculate the pressure lost in pipes due to friction. For a smooth surface, it can be shown that CDf = 0.074 (Re)x-1/5 ρu L (Re)x is the Reynolds number based on the length. (R e ) x = o µ 3 WORKED EXAMPLE 1.1 Calculate the drag force on each side of a thin smooth plate 2 m long and 1 m wide with the length parallel to a flow of fluid moving at 30 m/s. The density of the fluid is 800 kg/m3 and the dynamic viscosity is 8 cP. SOLUTION ρu L 800 x 30 x 2 (R e ) x = o = = 6 x 10 6 µ 0.008 C Df = 0.074 x (6 x10 6 ) 2 ρu 0 1 5 = 0.00326 800 x 30 2 = 360 kPa 2 2 x dynamic pressure = 0.00326 x 360 x 10 3 = 1173.6 Pa Dynamic pressure = τ w = C Df − = R = τ w x Wetted Area = 1173.6 x 2 x 1 = 2347.2 N On a small area the drag is dR = τw dA. If the body is not a thin plate and has an area inclined at an angle θ to the flow direction, the drag force in the direction of flow is τw dA cosθ. Fig.1.3 The drag force acting on the entire surface area is found by integrating over the entire area. R = τ w cos θ dA ∫ Solving this equation requires more advanced studies concerning the boundary layer and students should refer to the classic textbooks on this subject. SELF ASSESSMENT EXERCISE No. 1 1. A smooth thin plate 5 m long and 1 m wide is placed in an air stream moving at 3 m/s with its length parallel with the flow. Calculate the drag force on each side of the plate. The density of the air is 1.2 kg/m3 and the kinematic viscosity is 1.6 x 10-5 m2/s. (0.128 N) 2. A pipe bore diameter D and length L has fully developed laminar flow throughout the entire length with a centre line velocity uo. Given that the drag coefficient is given as ρu D CDf = 16/Re where Re = o , show that the drag force on the inside of the pipe is µ given as R=8πµuoL and hence the pressure loss in the pipe due to skin friction is pL = 32µuoL/D2 4 1.2 FORM DRAG and WAKES Form or pressure drag applies to bodies that are tall in comparison to the length in the direction of flow. Such bodies are called BLUFF BODIES. Consider the case below that could for example, be the pier of a bridge in a river. The water speeds up around the leading edges and the boundary layer quickly breaks away from the surface. Water is sucked in from behind the pier in the opposite direction. The total effect is to produce eddy currents or whirl pools that are shed in the wake. There is a build up of positive pressure on the front and a negative pressure at the back. The pressure force resulting is the form drag. When the breakaway or separation point is at the front corner, the drag is almost entirely due to this effect but if the separation point moves along the side towards the back, then a boundary layer forms and skin friction drag is also produced. In reality, the drag is always a combination of skin friction and form drag. The degree of each depends upon the shape of the body. Fig.1.4 The next diagram typifies what happens when fluid flows around a bluff object. The fluid speeds up around the front edge. Remember that the closer the streamlines, the faster the velocity. The line representing the maximum velocity is shown but also remember that this is the maximum at any point and this maximum value also increases as the stream lines get closer together. Fig.1.5 5 Two important effects affect the drag. Outside the boundary layer, the velocity increases up to point 2 so the pressure acting on the surface goes down. The boundary layer thickness δ gets smaller until at point S it is reduced to zero and the flow separates from the surface. At point 3, the pressure is negative. This change in pressure is responsible for the form drag. Inside the boundary layer, the velocity is reduced from umax to zero and skin friction drag results. Fig.1.6 In problems involving liquids with a free surface, a negative pressure shows up as a drop in level and the pressure build up on the front shows as a rise in level. If the object is totally immersed, the pressure on the front rises and a vacuum is formed at the back. This results in a pressure force opposing movement (form drag). The swirling flow forms vortices and the wake is an area of great turbulence behind the object that takes some distance to settle down and revert to the normal flow condition. Here is an outline of the mathematical approach needed to solve the form drag. Form drag is due to pressure changes only. The drag coefficient due to pressure only is denoted CDp and defined as Drag force C Dp = dynamic pressure x projected area 2R C Dp = 2 ρu 0 x projected area The projected area is the area of the outline of the shape projected at right angles to the flow. The pressure acting at any point on the surface is p. The force exerted by the pressure on a small surface area is p dA. If the surface is inclined at an angle θ to the general direction of flow, the force is p cosθ dA. The total force is found by integrating all over the surface. R = p cos θ dA ∫ The pressure distribution over the surface is often expressed in the form of a pressure coefficient defined as follows. 2( p - p o ) Cp = 2 ρu o po is the static pressure of the undisturbed fluid, uo is the velocity of the undisturbed fluid and 2 ρu o is the dynamic pressure of the stream. 2 6 Consider any streamline that is affected by the surface. Applying Bernoulli between an undisturbed point and another point on the surface, we have the following. 2 ρu 2 ρu o p0 + =p + 2 2 ρ2 p - p o = (u o - u 2 ) 2 ρ 2 2 (u o - u 2 ) 2 2 2 2(p - p o ) 2 = (u o - u ) = 1 − u Cp = = 2 2 2 2 ρu o ρu o uo uo In order to calculate the drag force, further knowledge about the velocity distribution over the object would be needed and students are again recommended to study the classic textbooks on this subject. The equation shows that if u<uo then the pressure is positive and if u>uo the pressure is negative. 1.3 TOTAL DRAG It has been explained that a body usually experiences both skin friction drag and form drag. The total drag is the sum of both. This applies to aeroplanes and ships as well as bluff objects such as cylinders and spheres. The drag force on a body is very hard to predict by purely theoretical methods. Much of the data about drag forces is based on experimental data and the concept of a drag coefficient is widely used. The DRAG COEFFICIENT is denoted CD and is defined by the following expression. CD = Resistance force Dynamic pressure x projected Area CD = 2R ρu x projected Area 2 o 7 WORKED EXAMPLE 1.2 A cylinder 80 mm diameter and 200 mm long is placed in a stream of fluid flowing at 0.5 m/s. The axis of the cylinder is normal to the direction of flow. The density of the fluid is 800 kg/m3. The drag force is measured and found to be 30 N. Calculate the drag coefficient. At a point on the surface the pressure is measured as 96 Pa above the ambient level. Calculate the velocity at this point. SOLUTION Projected area = 0.08 x 0.2 =0.016 m2 R = 30 N uo = 0.5 m/s ρ = 800 kg/m3 dynamic pressure = ρuo2/2 = 800 x 0.52/2 = 100 Pa CD = Resistance force 30 = = 18.75 Dynamic pressure x projected Area 100 x 0.016 ρ2 p - p o = (u o - u 2 ) 2 800 96 = (0.5 2 − u 2 ) 2 96 x 2 = (0.5 2 − u 2 ) 800 0.24 = 0.25 - u 2 u 2 = 0.01 u = 0.1 m/s 8 1.4 APPLICATION TO A CYLINDER 2R The projected Area is LD x projected Area where L is the length and D the diameter. The drag around long cylinders is more predictable than for short cylinders and the following applies to long cylinders. Much research has been ρu d carried out into the relationship between drag and Reynolds number. Re = o and d is the µ diameter of the cylinder. At very small velocities, (Re <0.5) the fluid sticks to the cylinder all the way round and never separates from the cylinder. This produces a streamline pattern similar to that of an ideal fluid. The drag coefficient is at its highest and is mainly due to skin friction. The pressure distribution shows that the dynamic pressure is achieved at the front stagnation point and vacuum equal to three dynamic pressures exists at the top and bottom where the velocity is at its greatest. The drag coefficient is defined as : CD = 2 ρu o Fig.1.7 As the velocity increases the boundary layer breaks away and eddies are formed behind. The drag becomes increasingly due to the pressure build up at the front and pressure drop at the back. Fig.1.8 9 Further increases in the velocity cause the eddies to elongate and the drag coefficient becomes nearly constant. The pressure distribution shows that ambient pressure exists at the rear of the cylinder. Fig.1.9 At a Reynolds number of around 90 the vortices break away alternatively from the top and bottom of the cylinder producing a vortex street in the wake. The pressure distribution shows a vacuum at the rear. Fig.1.10 Up to a Reynolds number of about 2 x 105, the drag coefficient is constant with a value of approximately 1. The drag is now almost entirely due to pressure. Up to this velocity, the boundary layer has remained laminar but at higher velocities, flow within the boundary layer becomes turbulent. The point of separation moves back producing a narrow wake and a pronounced drop in the drag coefficient. When the wake contains vortices shed alternately from the top and bottom, they produce alternating forces on the structure. If the structure resonates with the frequency of the vortex shedding, it may oscillate and produce catastrophic damage. This is a problem with tall chimneys and suspension bridges. The vortex shedding may produce audible sound. Fig. 1.12 shows an approximate relationship between CD and Re for a cylinder and a sphere. 10 SELF ASSESSMENT EXERCISE No. 2 1. Calculate the drag force for a cylindrical chimney 0.9 m diameter and 50 m tall in a wind blowing at 30 m/s given that the drag coefficient is 0.8. The density of the air is 1.2 kg/m3. (19.44 N) 2 Using the graph (fig.1.12) to find the drag coefficient, determine the drag force per metre length acting on an overhead power line 30 mm diameter when the wind blows at 8 m/s. The density of air may be taken as 1.25 kg/m3 and the kinematic viscosity as 1.5 x 10-5 m2/s. (1.8 N). 1.5 APPLICATION TO SPHERES The relationship between drag and Reynolds number is roughly the same as for a cylinder but ρu d it is more predictable. The Reynolds number is Re = o where d is the diameter of the µ sphere. The projected area of a sphere of diameter d is ¼ πd2. In this case, the expression for 8R the drag coefficient is as follows. C D = 2 ρu x πd 2 . At very small Reynolds numbers (less than 0.2) the flow stays attached to the sphere all the way around and this is called Stokes flow. The drag is at its highest in this region. As the velocity increases, the boundary layer separates at the rear stagnation point and moves forward. A toroidal vortex is formed. For 0.2<Re<500 the flow is called Allen flow. Fig.1.11 The breakaway or separation point reaches a stable position approximately 80o from the front stagnation point and this happens when Re is about 1000. For 500<Re the flow is called Newton flow. The drag coefficient remains constant at about 0.4. Depending on various factors, when Re reaches 105 or larger, the boundary layer becomes totally turbulent and the separation point moves back again forming a smaller wake and a sudden drop in the drag coefficient to about 0.3. There are two empirical formulae in common use. 24 6 For 0.2 < Re < 105 + + 0 .4 CD = Re 1+ Re 11 For Re <1000 CD = [ 24 1 + 0.15 Re 0.687 Re Fig. 1.12 shows this approximate relationship between CD and Re. Fig.1.12 WORKED EXAMPLE 1.3 A sphere diameter 40 mm moves through a fluid of density 750 kg/m3 and dynamic viscosity 50 cP with a velocity of 0.6 m/s. Note 1 cP = 0.001 Ns/m2. Calculate the drag on the sphere. Use the graph to obtain the drag coefficient. SOLUTION ρud 750 x 0.6 x 0.04 = = 360 µ 0.05 from the graph C D = 0.8 Re = CD = R= 2R 2 ρu x projected Area Projected area = π d2 0.04 2 =π = 1.2566 x 10 -3 m 2 4 4 C D ρu 2 A 0.8 x 750 x 0.6 2 x 1.2566 x 10 -3 = = 0.136 N 2 2 12 1.6 TERMINAL VELOCITY When a body falls under the action of gravity, a point is reached, where the drag force is equal and opposite to the force of gravity. When this condition is reached, the body stops accelerating and the terminal velocity reached. Small particles settling in a liquid are usually modelled as small spheres and the preceding work may be used to calculate the terminal velocity of small bodies settling in a liquid. A good application of this is the falling sphere viscometer described in earlier work. For a body immersed in a liquid, the buoyant weight is W and this is equal to the viscous resistance R when the terminal velocity is reached. πd 3 g (ρ s − ρ f ) R = W = volume x gravity x density difference = 6 ρs = density of the sphere material ρf = density of fluid d = sphere diameter STOKES’ FLOW For Re<0.2 the flow is called Stokes flow and Stokes showed that R = 3πdµut For a falling sphere viscometer, Stokes flow applies. Equating the drag force and the buoyant weight we get πd 3 g(ρ s − ρ f ) 3πdµu t = 6 2 d g(ρ s − ρ f ) µ= 18u t d 2 g(ρ s − ρ f ) 18µ This formula assumes a fluid of infinite width but in a falling sphere viscometer, the liquid is squeezed between the sphere and the tube walls and additional viscous resistance is produced. The Faxen correction factor F is used to correct the result. The terminal velocity for Stokes flow is u t = WORKED EXAMPLE 1.4 The terminal velocity of a steel sphere falling in a liquid is 0.03 m/s. The sphere is 1 mm diameter and the density of the steel is 7830 kg/m3. The density of the liquid is 800 kg/m3. Calculate the dynamic and kinematic viscosity of the liquid. SOLUTION Assuming Stokes’ flow the viscosity is found from the following equation. d 2 g(ρ s − ρ f ) 0.0012 x 9.81 x (7830 - 800) µ= = = 0.1277 Ns/m 2 = 127.7 cP 18u t 18 x 0.03 µ 0.1277 = = 0.0001596 m 2 / s = 159.6 cSt ρs 800 ρ ud 800 x 0.03 x 0.001 Check the Reynolds number. R e = f = = 0.188 µ 0.0547 As this is smaller than 0.2 the assumption of Stokes’ flow is correct. ν= 13 ALLEN FLOW For 0.2 < Re < 500 the flow is called Allen flow and the following expression gives the empirical relationship between drag and Reynolds number. CD=18.5Re-0.6 Equating for CD gives the following result. C D = Substitute R = CD = πd 3 g (ρ s − ρ f ) 6 8dg(ρ s − ρ f ) 6ρ f u 2 t 8dg(ρ s − ρ f ) 6ρ f u 2 t ρ u d − = 18.5R e 0.6 = 18.5 f t µ ρ u d = 18.5 f t µ 8R − = 18.5R e 0.6 2 2 ρ f u t πd −0.6 − 0.6 From this equation the velocity ut may be found. NEWTON FLOW For 500 < Re < 105 CD takes on a constant value of 0.44. Equating for CD gives the following. C D = Substitute R = 8dg(ρ s − ρ f ) 6ρ f u 2 t πd 3 g (ρ s − ρ f ) 6 8R = 0.44 ρ f u 2 πd 2 t = 0.44 29.73dg(ρ s − ρ f ) ρf When solving the terminal velocity, you should always check the value of the Reynolds number to see if the criterion used is valid. ut = 14 WORKED EXAMPLE 1.5 Small glass spheres are suspended in an up wards flow of water moving with a mean velocity of 1 m/s. Calculate the diameter of the spheres. The density of glass is 2630 kg/m3. The density of water is 1000 kg/m3 and the dynamic viscosity is 1 cP. SOLUTION First, try the Newton flow equation. This is the easiest. ut = d= 29.73d g(ρ s − ρ f ) ρf u 2ρ f 12 x 1000 t = = 0.0021 m or 2.1 mm 29.73 g(ρ s − ρ f ) 29.73 x 9.81 x (2630 - 1000) Check the Reynolds number. ρ f u t d 1000 x 1 x 0.0021 = = 2103 µ 0.001 The assumption of Newton flow was correct so the answer is valid. Re = WORKED EXAMPLE 1.6 Repeat the last question but this time with a velocity of 0.05 m/s. Determine the type of flow that exists. SOLUTION If no assumptions are 24 6 formula C D = + + 0 .4 Re 1+ Re made, we should use the general ρ f u t d 1000 x 0.05 x d = = 50 000d µ 0.001 24 6 CD = + + 0.4 Re 1+ Re Re = 24 6 + + 0.4 50000d 1 + 50000d 6 C D = 0.00048d −1 + + 0.4 1 + 223.6d 0.5 CD = CD = 8dg(ρ s − ρ f ) 6ρ f u 2 = 8d x 9.81 x (2630 - 1000) = 8528.16d 6 x 1000 x 0.05 2 8528.16d = 0.00048d −1 + 6 + 0.4 1 + 223.6d 0.5 15 This should be solved by any method known to you such as plotting two functions and finding the point of interception. f1(d) = 8528.16d f2(d) = 0.00048d −1 + 6 + 0.4 1 + 223.6d 0.5 The graph below gives an answer of d = 0.35 mm. Fig. 1.13 ρ f u t d 1000 x 0.05 x 0.00035 = = 17.5 µ 0.001 This puts the flow in the Allen's flow section. Checking the Reynolds’ number R e = 16 ANOTHER METHOD OF SOLUTION It has been shown previously that the drag coefficient for a sphere is given by the formula CD = 8R . R is the drag force. One method of solving problems is to arrange the πd 2 ρu 2 formula into the form CD Re2 as follows. CD = 8R πd ρ f u 2 C D Re2 = 2 x ρ f µ2 ρfµ = 2 8Rρ f πµ 2 x µ2 ρud 2 f 2 2 = 8Rρ f πµ 2 x 1 Re2 8Rρ f πµ 2 If the sphere is falling and has reached its terminal velocity, R = buoyant weight. R= πd 3 g (ρ s − ρ f ) 6 CD R = 2 e C D Re2 = πd 3 g (ρ s − ρ f )8 ρ f 6πµ 2 4d 3 gρ f (ρ s − ρ f 3µ 2 ) ...............( A) The drag coefficient for a sphere is related to the Reynolds number as described previously. There are two equations commonly used for this relationship as follows. CD = and CD = 24 6 + + 0.4..................( B) Re 1 + Re [ 24 0.687 1 + 0.15Re ..................(C ) Re Either B or C may be used in the solution of problems. The general method is to solve ReCD2 from equation A. Next compose a table of values of Re, CD, and ReCD2. Plot ReCD2 vertically and Re horizontally. Find the value of Re that gives the required value of ReCD2. From this the velocity may be deduced. 17 WORKED EXAMPLE 1.7 A sphere 1.5 mm diameter falls in water. The density of the sphere is 2500 kg/m3. The density and dynamic viscosity of water is 997 kg/m3 and 0.89 x 10-3 Ns/m2 respectively. The drag coefficient is given by the formula C D = [ 24 0.687 1 + 0.15Re . Determine the terminal Re velocity. SOLUTION 4d 3 gρ f (ρ s − ρ f ) 4(0.0015) 3 x9.81x997(2500 − 997 ) = = 83513 CD R = 3µ 2 3(0.89 x10 −3 ) 2 24 1 + 0.15 Re 0.687 . Next compile a table using the formula C D = Re 2 e [ Re CD CDRe2 0.1 24.7 2.47 1 27.6 27.6 10 4.15 415 100 1.09 109 1000 0.44 438288 We are looking for a value of CDRe2 = 83513 and it is apparent that this occurs when Re is between 100 and 1000. By plotting or by narrowing down the figure by trial and error we find that the correct value of Re is 356. Re = 356 = ρfud/µ 356 = 997 x u x 0.0015/0.89 x 10-3 u = 0.212 m/s and this is the terminal velocity. 18 SELF ASSESSMENT EXERCISE No. 3 1. a. Explain the term Stokes flow and terminal velocity. b. Show that the terminal velocity of a spherical particle with Stokes flow is given by the formula u = d2g(ρs - ρf)/18µ Go on to show that CD=24/Re 2. Calculate the largest diameter sphere that can be lifted upwards by a vertical flow of water moving at 1 m/s. The sphere is made of glass with a density of 2630 kg/m3. The water has a density of 998 kg/m3 and a dynamic viscosity of 1 cP. (20.7 mm) 3. Using the same data for the sphere and water as in Q2, calculate the diameter of the largest sphere that can be lifted upwards by a vertical flow of water moving at 0.5 m/s. (5.95 mm). 4. Using the graph (fig. 1.12) to obtain the drag coefficient of a sphere, determine the drag on a totally immersed sphere 0.2 m diameter moving at 0.3 m/s in sea water. The density of the water is 1025 kg/m3 and the dynamic viscosity is 1.05 x 10-3 Ns/m2. (0.639 N). 5. A glass sphere of diameter 1.5 mm and density 2 500 kg/m3 is allowed to fall through water under the action of gravity. The density of the water is 1000 kg/m3 and the dynamic viscosity is 1 cP. Calculate the terminal velocity assuming the drag coefficient is CD = 24 Re -1 (1+ 0.15Re 0.687) (Ans. 0.215 m/s 6. Similar to part of Q1 1990 A glass sphere of density 2 690 kg/m3 falls freely through water. Find the terminal velocity for a 4 mm diameter sphere and a 0.4 mm diameter sphere. The drag coefficient is CD = 8F/{πd2ρu2} This coefficient is related to the Reynolds number as shown. Re CD 15 3.14 20 2.61 25 2.33 30 2.04 35 1.87 The density and viscosity of the water is 997 kg/m3 and 0.89 x 10-3 N s/m2. Answer 0.45 m/s and 0.06625 m/s. 7. Similar to part of Q4 1988. A glass sphere of diameter 1.5 mm and density 2 500 kg/m3 is allowed to fall through water under the action of gravity. Find the terminal velocity assuming the drag coefficient is CD = 24 Re-1(1+ 0.15Re0.687) (Ans. 0.215 m/s) 19 8. Similar to Q1 1986 The force F on a sphere of diameter d moving at velocity um in a fluid is given by F = CD {πd2ρum2}/8 For Reynolds numbers less than 1 000, CD is given by CD= 24 Re-1(1+ 0.15Re0.687) Estimate the terminal velocity of a glass sphere 1 mm diameter and density 2 650 kg/m3 in water of density 997 kg/m3 and viscosity 0.89 x 10-3 N s/m2. Answer 0.15 m/s 20 2. BOUNDARY LAYERS In this tutorial we will look at the shape of various types of boundary layers. We will look at the mathematical equations for the shape of the boundary layer and use them to solve problems. You may recall that the definition of a BL is the thickness of that layer next to a surface in which the velocity grows from zero to a maximum value (or so close to a maximum as to be of no practical difference). This thickness is usually given the symbol δ (small delta). The boundary layer, once established may have a constant thickness but, for example, when a flow meets the leading edge of a surface, the boundary layer will grow as shown (fig.2.1). Fig.2.1 When the flow enters a pipe the BL builds up from all around the entrance and a cross section shows the layer meets at the centre (fig.2.2) Fig.2.2 The symbol u1 is used to designate the maximum velocity in the fully developed layer. The fully developed layer may be laminar or turbulent depending on the Reynolds’ Number. 21 The velocity profile for a typical case is shown on fig.2.3. Fig.2.3 The shear stress between any two horizontal layers is τ. For a Newtonian Fluid the relationship between shear stress, dynamic viscosity (µ) and rate of shear strain (du/dy) is τ = µ du/dy At the wall the shear stress is called the WALL SHEAR STRESS, τo and occurs at y = 0. Note that the gradient du/dy is the rate of shear strain and it is steeper for turbulent flow than for laminar flow giving a greater shear resistance. The solution of problems is simplified by the concepts of DISPLACEMENT THICKNESS AND MOMENTUM THICKNESS which we will now examine. 2.1. DISPLACEMENT THICKNESS δ* The flow rate within a boundary layer is less than that for a uniform flow of velocity u1. The reduction in flow is equal to the area under the curve in fig.2.3. If we had a uniform flow equal to that in the boundary layer, the surface would have to be displaced a distance δ* in order to produce the reduction. This distance is called the displacement thickness and it is given by : δ flow redution = ∫ [u 1 − u ]dy = u1δ ∗ 0 If this is equivalent to a flow of velocity u1 in a layer δ* thick then : δ u 1 δ ∗ = ∫ 1 − dy u 0 2.2. MOMENTUM THICKNESS θ The momentum in a flow with a BL present is less than that in a uniform flow of the same thickness. The momentum in a uniform layer at velocity u1 and height h is ρhu12. When a BL exists this is reduced by ρu12θ. Where θ is the thickness of the uniform layer that contains the equivalent to the reduction. Using the same reasoning as before we get : δ u u θ = ∫ 1 − dy u u1 0 1 22 2.3. BOUNDARY LAYER LAWS The velocity at any distance y above a surface is a function of the wall shear stress, the dynamic viscosity and the density. u =φ(y, τo, ρ ,µ) If you are familiar with the method of dimensional analysis you may wish to show for yourself that : u(ρ/τo)½ = φ{y(τo)½} Generally the law governing the growth of a BL is of the form u = φ(y) and the limits must be that u = 0 at the wall and u = u1 in the fully developed flow. There are many ways in which this is expressed according to the Reynolds’ Number for the flow. The important boundary conditions that are used in the formulation of boundary layer laws are: 1. The velocity is zero at the wall ( u = 0 @ y = 0). 2. The velocity is a maximum at the top of the layer ( u = u1 @ y = δ). 3. The gradient of the b.l. is zero at the top of the layer (du/dy = 0 @ y = δ). 4. The gradient is constant at the wall (du/dy = C @ y = 0). 5. Following from (4) d2u/dy2 = 0 @ y = 0). Let us start by considering LAMINAR BOUNDARY LAYERS. 23 LAMINAR BOUNDARY LAYERS 2.3.1 One of the laws which seem to work for laminar flow is u = u1 sin(π y/2δ) WORKED EXAMPLE No.2.1 Find the displacement thickness δ* for a Laminar BL modelled by the equation u = u1 sin(π y/2δ) δ δ u πy δ ∗ = ∫ 1 − dy = ∫ 1 − sin u 2δ 0 1 0 δ 2δ 2δ πy δ = y + cos = {δ + 0} − 0 + = 0.364δ π π 2δ 0 ∗ Another way of expressing the shape of the laminar BL is with a power law. The next example is typical of that used in the examination. WORKED EXAMPLE No.2.2 The velocity distribution inside a laminar BL over a flat plate is described by the cubic law : u/u1 = a0 + a1y + a2y2 + a3y3 Show that the momentum thickness is 39δ/280 SOLUTION At y = 0, u = 0 so it follows that a0 = 0 d2u/dy2 = 0 @ y = 0 so a2= 0. Show for yourself that this is so. The law is reduced to at y = δ, u = u1 so u/u1 = a1y + a3y3 1 = a1δ + 3a3δ2 hence a1 =(1-a3δ3)/δ Now differentiate and du/dy = u1(a1+3a3y2) at y = δ, du/dy is zero so 0 =a1+ 3a3δ2 so a1= -3a3δ2 Hence by equating a1= 3/2δ and a3 = -1/2δ3 Now we can write the velocity distribution as and u/u1=3y/2δ -(y/δ)3/2 du/dy = u1{3/2δ + 3y2/2δ3} 24 If we give the term y/δ the symbol η we may rewrite the equation as: u/u1=3η/2 - η3/2 The momentum thickness θ is given by : δ u u θ = ∫ 1 − dy but dy = δ dη u u1 0 1 1 3η η 3 3η η 3 θ = ∫ − 1 − − dη 0 2 2 2 2 Integrating gives : 3η 2 η 4 9η 3 η 7 3η 5 θ =δ − − − + 8 12 28 10 4 between the limits η=0 and η=1 this evaluates to θ = 39δ/280 25 WORKED EXAMPLE No.2.3 Show that δ/x=4.64Re-0.5 for the same case as before. SOLUTION We must first go back to the basic relationship. From the previous page du/dy = u1{3/2δ + 3y2/2δ3} At the wall where y=0 the shear stress is τo=µdu/dy = µu1{3/2δ + 3y2/2δ3} =(µu1/δ) δ [(3/2δ) +3y2/2δ3] Putting y/δ= η we get : τo = (µu1/δ) δ [(3/2δ) +3δ2/2δ] τo = (µu1/δ) [(3/2) +3δ2/2] at the wall η=0 τo = (µu1/δ) (3/2)………........(2.1) The friction coefficient Cf is always defined as Cf = τo/(ρ u12/2).............(2.2) It has been shown elsewhere that Cf = 2dθ/dx. The student should search out this information from test books. Putting θ = 39δ/280 (from the last example) then Cf = 2dθ/dx = (2x39/280) dδ/dx ..............(2.3) equating (2.2) and (2.3) gives τo=( ρu12)(39/280)dδ/dx .............(2.4) equating (2.1) and (2.4) gives (ρu12)(39/280)dδ/dx = (µu/δ)(3/2) hence (3 x 280)/(2 x 39)(µdx)/ρu) = δdδ Integrating 10.77(µx/ ρu1) = δ2/2 + C Since δ = 0 at x = 0 (the leading edge of the plate) then C=0 hence δ= {21.54µx/ρu1}½ dividing both sides by x gives NB δ/x = 4.64(µ/ρu1x)-½ = 4.64Re-½ Rex= ρu1x/µ and is based on length from the leading edge. 26 SELF ASSESSMENT EXERCISE No. 4 1. The BL over a plate is described by u/u1=sin(πy/2δ). Show that the momentum thickness is 0.137δ. 2. The velocity profile in a laminar boundary layer on a flat plate is to be modelled by the cubic expression u/u1=a0+a1y + a2y2+a3y3 where u is the velocity a distance y from the wall and u1 is the main stream velocity. Explain why a0 and a2 are zero and evaluate the constants a1 and a3 in terms of the boundary layer thickness δ. Define the momentum thickness θ and show that it equals 39δ/280 Hence evaluate the constant A in the expression δ/x = A (Rex)-0.5 where x is the distance from the leading edge of the plate. It may be assumed without proof that the friction factor Cf = 2 dθ/dx 3. (a) The velocity profile in a laminar boundary layer is sometimes expressed in the form u/u1=a0+a1(y/δ)+a2(y/δ)2+a3 (y/δ)3+a4(y/δ)4 where u1 is the velocity outside the boundary layer and δ is the boundary layer thickness. Evaluate the coefficients a0 to a4 for the case when the pressure gradient along the surface is zero. (b) Assuming a velocity profile u/u1=2(y/δ) - (y/δ)2 obtain an expression for the mass and momentum fluxes within the boundary layer and hence determine the displacement and momentum thickness. 4. When a fluid flows over a flat surface and the flow is laminar, the boundary layer profile may be represented by the equation 2 u/u1= 2(η) - (η) where η = y/δ y is the height within the layer and δ is the thickness of the layer. u is the velocity within the layer and u1 is the velocity of the main stream. Show that this distribution satisfies the boundary conditions for the layer. Show that the thickness of the layer varies with distance (x) from the leading edge by the equation δ=5.48x(Rex)-0.5 2 It may be assumed that τo = ρu1 dθ/dx 27 5. Define the terms displacement thickness δ * and momentum thickness θ. Find the ratio of these quantities to the boundary layer thickness δ if the velocity profile within the boundary layer is given by u/u1=sin(πy/2δ) Show, by means of a momentum balance, that the variation of the boundary layer -0.5 thickness δ with distance (x) from the leading edge is given by δ = 4.8(Rex) It may be assumed that τo = ρu12 dθ/dx Estimate the boundary layer thickness at the trailing edge of a plane surface of length 0.1 m when air at STP is flowing parallel to it with a free stream velocity u1 of 0.8 m/s. It may be assumed without proof that the friction factor Cf is given by Cf = 2 dθ/dx N.B. standard data µ = 1.71 x 10-5 N s/m2. ρ = 1.29 kg/m3. 6. In a laminar flow of a fluid over a flat plate with zero pressure gradient an approximation to the velocity profile is u/u1=(3/2)(η) - (1/2)(η)3 η = y/δand u is the velocity at a distance y from the plate and u1 is the mainstream velocity. δ is the boundary layer thickness. Discuss whether this profile satisfies appropriate boundary conditions. Show that the local skin-friction coefficient Cf is related to the Reynolds’ number (Rex) based on distance x from the leading edge by the expression Cf =A (Rex)-0.5 and evaluate the constant A. It may be assumed without proof that Cf = 2 dθ/dx and that θ is the integral of (u/u1)(1 - u/u1)dy between the limits 0 and δ 28 2.3.2 TURBULENT BOUNDARY LAYERS When a fluid flows at high velocities, the boundary layer becomes turbulent and the gradient at the wall becomes smaller so the wall shear stress is larger and the drag created on the surface increases. Fig. 2.4 Prandtl found that a law which fits the turbulent case well for Reynolds’ numbers below 107 is: u = u1(y/δ)1/7 This is called the 1/7th law. The gradient of the B.L. is du/dy = u1δ1/7y-6/7/7 This indicates that at the wall where y=0, the gradient is infinite (horizontal). This is obviously incorrect and is explained by the existence of a laminar sub-layer next to the wall. In this layer the velocity grows very quickly from zero and merges with the turbulent layer. The gradient is the same for both at the interface of laminar and turbulent flow. The drag on the surface is due to the wall shear stress in the laminar sub-layer. 29 WORKED EXAMPLE No.2.4 Show that the mean velocity in a pipe with fully developed turbulent flow is 49/60 of the maximum velocity. Assume the 1/7th law. For a pipe, the B.L. extends to the centre so δ = radius = R. Consider an elementary ring of flow. Fig.2.5 The velocity through the ring is u. The volume flow rate through the ring is 2πrudr The volume flow rate in the pipe is Q = 2π∫rudr 1/7 Since δ = R then u = u1(y/R) also r = R-y -1/7 1/7 Q = 2π∫ (R-y)udr = 2π∫u1R (R-y)y dy -1/7 1/7 8/7 Q =2πu1R [Ry -y ] -1/7 8/7 15/7 Q =2πu1R [(7/8)Ry - (7/15)y 2 Q =(49/60)πR u1. 2 The mean velocity is defined by um=Q/πR hence um=(49/60)u1 30 2.4 FRICTION COEFFICIENT AND BOUNDARY LAYERS Earlier it was explained that the friction coefficient Cf is the ratio of the wall shear stress to the dynamic pressure so : Cf = 2τo/(ρum2) ..........................(2.4.1) For smooth walled pipes, Blazius determined that Cf = 0.079Re-0.25 ........................(2.4.2) 2τo/(ρum2) = 0.079Re-0.25 Equating (2.4.1) and (2.4.2) gives : Note that um is the mean velocity and u1 is the maximum velocity. Research shows that um=0.8u1 Also Note that Re=ρu1D/µ and D = 2δ. Hence τo=0.02125ρu12(µ/ρδu1) 0.25.............(2.4.3) 2.5 FORCE BALANCE IN THE BOUNDARY LAYER The student should refer to textbooks for finer details of the following work. Consider again the growth of the B.L. as the fluid comes onto a flat surface. A stream line for the flow is not parallel to the B.L. Now consider a control volume A, B, C, D. Fig.2.6 Balancing pressure force and shear force at the surface with momentum changes gives : δ δ δ δu τ o = ρ ∫ [u − u1 ]udy + ρ ∫ [u − u1 ]udy...........(2.4.4) δx 0 δx 0 Using equations (2.4.2), (2.4.3) and (2.4.4) gives (4/5)δ5/4= 0.231{µ/ρu}1/5x = Re-1/5 Fs= τo x surface area The surface skin friction coefficient is Cf=2Fs/(ρu12) = 0.072Re-1/5 Experiments have shown that a more accurate figure is : Cf= 0.074Re-1/5 The shear force on the surface is 31 SELF ASSESSMENT EXERCISE No. 5 1. Under what circumstances is the velocity profile in a pipe adequately represented by 1/7 the 1/7 th power law u/u1=(y/R) where u is the velocity at distance y from the wall, R is the pipe radius and u1 is the centre-line velocity ? The table shows the measured velocity profile in a pipe radius 30 mm. Show that these data satisfy the 1/7 th power law and hence evaluate (i) the centre-line velocity (ii) the mean velocity um (iii) the distance from the wall at which the velocity equals um. 1.0 2.0 1.54 1.70 1.94 2. 5.0 2.14 10.0 2.26 15.0 2.36 20.0 u (m/s) y (mm) 1/7 (a) Discuss the limitations of the 1/7th power law u/u1=(y/R) for the velocity profile in a circular pipe of radius R, indicating the range of Reynolds numbers for which this law is applicable. (b) Show that the mean velocity is given by 49u1/60. (c) Water flows at a volumetric flow rate of 1.1 x 10-3 m3/s in a tube of diameter 25 mm. Calculate the centre-line velocity and the distance from the wall at which the velocity is equal to the mean velocity. -0.25 (d) Assuming that Cf=0.079(Re) evaluate the wall shear stress and hence estimate the laminar sub-layer thickness. µ = 0.89 x 10-3 N s/m2. ρ=998 kg/m3. 32 ...
View Full Document

This note was uploaded on 01/05/2011 for the course CSU 3 taught by Professor Handsome during the Spring '10 term at CSU Pueblo.

Ask a homework question - tutors are online