DEPARTMENT OF TECHNOLOGY Moment Balance of an axisymmetric jet and the effect of air entrainment from ambient Xue Fei May 2010 Master s Thesis in Building Environment Modeling - CFD, Measuring Techniques and Visualization Supervisor: Mats Sandberg Examiner: Taghi Karimipanah
Preface The author wishes to thank supervisor Mats Sandberg for his patient help and crucial discussions, Claes Blomqvist for his kindly help with the instruments. The author also would like to thank Taghi Karimipanah for helping me access to this project and essential discussions. I
II
Abstract The momentum flux of turbulent impinging jet is not conserved causing a transfer of momentum between the jet and its ambient, which may result in either increasing or decreasing the momentum flux of the jet. In order to study the nature of the momentum variation, the reaction force of an axisymmetric jet was measured by issuing the jet from a round nozzle against a plate which is placed on a digital balance. Tests were carried out in an almost infinite ambient with different heights for Reynolds between 5.5 10 and 6.4 10. Different sizes of ceilings were installed around the nozzle to examine the influence of the air entrainment from ceiling level into the jet. Visualization was performed to understand better the effect of the entrainment. The results show, after a certain distance, the reaction force becomes larger than the momentum flux at the inlet. By using wider ceilings around the jet supply outlet, the reaction force, which is normalized by the inlet momentum, was slightly decreased. III
IV
Content 1 Introduction...1 2 Sensitivity of the balance.........3 3 Method of Test......5 4 Experimental Set-up....9 5 Results.....13 5.1 Jet with a small ceiling on top of the control volume.....13 5.2 Jet with a round sheet on top of the control volume......16 5.3 Jet with an extended sheet on top of the control volume...22 5.4 General comparing....26 6 Conclusion....... 29 7 References......31 Appendix...33 Ⅰ Details for L=0.12 m at air flow 16.7 l/s.......35 Ⅱ Details for L=0.12 m at air flow 18.7 l/s...37 Ⅲ Details for L=2 m at air flow 16.7 l/s...39 Ⅳ Details for L=2 m at air flow 18.7 l/s...41 Ⅴ Details for L=3 m at air flow 16.7 l/s...43 Ⅵ Details for L=3 m at air flow 18.7 l/s...45 V
VI
Nomenclature A: Area of impinging plate (A= D 2 /4) [m 2 ] A f : Free area of supply device [m 2 ] D: Impinging plate diameter [m] d: Diameter of nozzle [m] da: Area differential [m 2 ] F: Flow force (F=M+pA) [N] g: Gravity acceleration [m/s 2 ] h: Impinging height [m] L: Diameter of the Ceiling around nozzle [m] M: Kinematic momentum flux (M U u da) [N] M(0): Exit kinematic momentum flux at nozzle [N] q: Volumetric flow rate [m 3 /s] R: Reaction force [N] U: Instantaneous axial velocity [m/s] : Fluctuating part of axial velocity [m/s] V: Velocity of the air coming in the control volume [m/s] V r : Radial velocity component of flowing in [m/s] x: Coordinate in axial direction [m] θ: Entrainment angle [radian] ρ: Air density [Kg/m 3 ] VII
VIII
1. Introduction A turbulent impinging jet has different functions and is useful in many engineering applications, such as cooling of cast iron, glass tempering and ventilation. It has also been widely studied with different configurations. S.Maurel and C.Solliec [1] studied the plane turbulent jets impinging normally to a flat and smooth plate. Carlo Carcasci [2] investigated air impinging jets by visualisation methods. Karimipanah and Sandberg [3] studied the confinement effect of Jets in ventilated Rooms. Heskestad [4] performed hot wire measurements in a plane turbulent jet. However, the results of kinematic momentum flux of a jet in a still ambient are controversial. Many researchers showed decreasing in the momentum, for instance Brandbury (1965) [5], Goldschmidt (1975) [6], Gutmark (1976) [7], and Kotsovinos (1975) [8]. Others reported increasing momentum, as e.g. Hussain and Clark (1977) [9], Panchapakesan and Lumley (1993) [10]. Result of constancy momentum also reported by Hussein (1994) [11]. Kotsovinos and Angelidis (1991) [12] attributed the increasing or decreasing of the jet momentum flux to the direction of induced flow that enter the jet. Hussein (1994) [11] reported that the previous measurements of increasing or decreasing of momentum are not valid dependent to the model scale. Karimipanah and Sandberg [13] discussed the decay of momentum and velocity in an axisymmetric impinging jet. Karimipanah and Sandberg [14] also studied the momentum balance of an axisymmetric impinging jet both in infinite environment and within enclosure. In this work, the momentum balance of a turbulent axisymmetric jet is studied with special interest on the influence of entrainment of ambient air from ceiling level into the jet. The jets with 2 different air flow rates were investigated and, up to a distance equal to 46 exit diameters, its forces 1
were measured by a digital balance and recorded by a personal computer. 2
2. Sensitivity of the balance To investigate the sensitivity of the digital balance, Mettler ID1 weighting terminal, different activities were performed with different distance from the balance to study how to secure the weighting results of the balance from people s activities around it. Table.1 Sensitivity of the balance with people s activities from different distance Distance (m) 0 0.5 1 2 2.5 3 Weight Range (g) (-13)~14 (-3)~1 (-2)~1 (-1)~2 0 0 Table.1 shows that to secure the weighting results, during the experiment, a protected area should be set with, at lease, 2.5 m as radius. In order to understand how long time common activities can affect the weighting results. The influence of running and walking close to the digital balance has been examined (See Figure.1). 5 times tests were performed for each behavior. Figure.1 Running and walking close to the balance. 3
Figure.2 and figure.3 show that after each adjustment, measurement should be carried out about 20 seconds later. 15 10 Weight (g) 5 0-5 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 test 1 test 2 test 3 test 4 test 5-10 Measuring time (s) Figure.2. Sensitivity test for walking close to the balance 1-1 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 Weight (g) -3-5 -7-9 -11-13 test 1 test 2 test 3 test 4 test 5-15 Measuring time (s) Figure.3 Sensitivity test for running close to the balance 4
3. Method of Test An axisymmetric jet was issued from a round nozzle and perpendicularly impinges on the centre of a black circular plate. The flow force of the axisymmetric jet on the impingement plate was weighted by a balance which is mounted under the impingement plate, See figure.4. The kinematic momentum flux of a free turbulent jet is defined as M = U 2 + u 2 da (1) Where U is the instantaneous axial velocity and uu is the fluctuating part of axial velocity. If one assumes that the air pressures under the plate and at the top of the control volume are equal to the pressure in the ambient, then the reaction force is h 2 AA uuuu RR (h) = MM 0 ππππ ρρvvrr 0 tttttttt(xx)dddd + AAnnnnnnnnnnnn ρρuu 2 dddd (2) Where VV rr is the radial velocity component of the entrainment velocity. One assumes θ is the angle that the entrained air enters the control volume. (See figure.4) 5
Figure.4 Schematic view of the test bench without ceiling on top of the control volume With special interest of the air entrainment from the ceiling level, a round sheet was placed around the nozzle, See figure.5, and the net force becomes RR (h) = MM 0 ππππ 0 h ρρvvrr 2 tttttttt(xx)dddd (3) 6
Figure.5 Schematic view of the test bench with ceiling on top of the control volume This ceiling was used to avoid the entrainment of ambient air from the top. Although the realistic solution which is based on the totally closed ceiling in case of ventilation applications, but the ceiling was not closed properly and this small opening caused some undesired entrainment. (See figure.6) 7
Figure.6 Test bench with ceiling on top of the control volume and the undesired entrainment. 8
4. Experimental Set-up The support for the nozzle was firmly clamped to the ground to avoid vibrations. A fan with variable power was adopted to supply different air flow rates. The device for measuring the flow rate is an orifice plate with 50 mm hole-diameter and a micro manometer, Swema 80, which was used for measuring the pressure difference between the two sides of the orifice plate. Then, with a thermometer, Digitron 2006T, and a barometer, Mechanism LTD No.M2236A, which were used for measuring local temperature and pressure respectively, air flow rate could be calculated. To calculate the air density of the jet, temperature at the nozzle was measured before sampling for each height. A digital balance, Mettler ID1 weighting terminal, with 1 gram resolution and a circular black plate with 1000 mm diameter were placed at downstream of the nozzle for measuring the weight of the jet and the results were recorded by a personal computer. For each measurement, 200 samples were taken. A well designed round nozzle with 25 mm diameter was adopted for issuing a three dimensional axisymmetric jet which was perpendicular to the centre of the impingement plate [15]. The impinging height h was arranged from about 22 to 1144 mm downstream of the nozzle. The distances between measuring points was equal to about 2 AA ff. Three different sizes of ceilings were adopted for the control volume. The diameters are 0.12 m, 2 m, 3 m respectively (See figure.7). Hysteresis appears at a few regions after installing the ceilings with L=2 m and L=3 m. To avoid it, a small board was used to horizontally blow the hysteresis air away before measures. A Philips Focus Generator PCV740 was used to make light sheet for smoke visualization. Canon D100 was used for taking long exposure photos, and the exposure time was set by 10 to 15 seconds. 9
(a) (b) 10
(c) Figure.7 Photos of Experimental Set-Up with different sizes of ceilings. (a) L=0.12 m; (b) L=2 m; (c) L=3 m. 11
12
5. Results 5.1 Jet with a small ceiling on top of the control volume Measurement was performed in an almost infinite environment. Ambient air could entrain the control volume from all the directions. (Details can be found in Appendix Ⅰ & Ⅱ) 5.1.1 Weighting Figure.8 presents the recorded reaction force normalized by the inlet momentum flux MM (0) varying with x/d. It is noted that the normalized reaction force increases rapidly before x/d=4, then it slowly reaches its peak value after x/d=10 approximately, which is in agreement with the findings of Karimipanah and Sandberg [14]. After x/d=40, the force starts to decrease. The force increases slightly with Reynolds number and it becomes larger than the inlet momentum flux after about x/d=2. 1.4 1.2 R/M(0) 1 0.8 0.6 0.4 0.2 q=16.7 l/s Re=57000 L=0.12 m q=18.7 l/s Re=64300 L=0.12 m 0 0 10 20 30 40 50 x/d Figure.8 Measured reaction force versus x/d 13
The results were also analyses in terms of the impinging distance normalized by the diameter of the ceiling of the control volume; the holder of the nozzle (see figure.9) is the smallest ceiling with 0.12 m diameter. Figure.10 shows that an inflection point of the profile of the normalized reaction force happens at about x/l=1. Figure.9 Supply terminal [16] 1.4 R/M(0) 1.2 1 0.8 0.6 0.4 0.2 q=16.7 l/s Re=57000 L=0.12 m q=18.7 l/s Re=64300 L=0.12 m 0 0 2 4 6 8 10 12 x/l Figure.10 Measured reaction force versus x/l 14
5.1.2 Visualization To visualize how the entrainment works, smoke visualization and long exposure camera offers an opportunity to see the view of the entrainment. Figure.11 shows the overall two dimension view. Note that the entrainment happens all the way along the jet. The angle of the entrainment increases as approaching to the plate. Two vortexes are found under the holder of the nozzle, See Figure.12, which may stop the directly air entrainment. The vortexes always appear until the distance between the plate and the ceiling is smaller than the length of the vortex. And the length of the vortex is approximately equal to the distance where the inflection point appears. Figure.11 Photo of air entrainment for Re=57000, q=16.7 l/s, x/d=30, L=0.12 m 15
Figure.12 Photo of air entrainment for Re=57000, q=16.7 l/s, x/d=9, L=0.12 m 5.2 Jet with a round sheet on top of the control volume A round sheet with 2 m diameter was installed to prevent ambient air entrainment from the ambient above the nozzle. (Details can be found in Appendix Ⅲ & Ⅳ) 5.2.1 Weighting Measurements were also performed for two different air flow rates. Hysteresis appears at the region close to the inflection point for both air flow rates; the trajectory of the normalized reaction force did not follow the same path when the experiment performed from different directions (see figure 13). Below the inflection point, ambient air entrainment is stopped by the spreading wall jet (see figure.14). Figure.15 shows the normalized reaction force is quite similar between the two different Reynolds numbers. This shows that the interaction between the jet and its ambient remains almost the same. But the inflection point of the profile moves forward till 16
about x/d=15 to x/d=20. There is a sharply jump at the inflection point, it shows that air entrainment has strong influence on the reaction force cause air starts to entrain the control volume after the inflection point. By the inflection point, it is also able to define two regions which are deentrainment region and entrainment region respectively. Then, the profile becomes relatively stable, it slowly decreases after x/d=40. Under x/d=7, the value of the normalized reaction force is minus. It represents that the spreading wall jet has strong influence in the region between the ceiling and the impinging plate and it may strongly interact with the impinging jet. Figure.16 shows the inflection point appears at about x/l=0.2. Figure.13 Hysteresis at the inflection point, Re=57400, q=16.8 l/s, L=2m Figure.14 Ambient air entrainment stopped by the spreading wall jet 17
Figure.15 Measured reaction force with a round sheet on top of the control volume versus x/d Figure.16 Measured reaction force with a round sheet on top of the control volume versus x/l 18
5.2.2 Visualization According to the visualization, even the round sheet was installed, ambient air is still able to entrain to the control volume from the top of the round sheet. (See Figure.17 and Figure.18) Figure.17 Photo of air entrainment for Re=57400, q=16.7 l/s, x/d=46, L=2 m Figure.18 Photo of air entrainment for Re=57400, q=16.7 l/s, x/d=30, L=2 m 19
Below about x/d=15 (See figure.19), Direct entrainment was stopped by the spreading wall jet, which is due to impingement of the jet on the impinging plate and subsequently the radial deflection of the jet trajectory. Approximately, it is also the distance where the inflection point appears. Figure.19 Photo of jet spreading outward at the edge of the ceiling for Re=57400, q=16.7 l/s, x/d=14, L=2 m Spiral rotation appears at a region close to the nozzle (See figure 20); it is due to the uneven air entrainment from the different sizes of openings at different sides of the ceiling. 20
Figure.20 photo of spiral rotation at the region close to the nozzle for Re=57400, q=16.7 l/s, x/d=30, L=2 m 5.3 Jet with an extended sheet on top of the control volume The round sheet was extended by hard plastic. Approximately the diameter of the sheet increased to 3 meters. ( Details can be found in Appendix Ⅴ & Ⅵ ) 5.3.1 Weighting Figure.21 shows that, with the extended sheet, the results profiles for the two air flow rates are quite similar. It is due to the same magnitude influence from the ambient air. The value of the normalized force becomes minus under x/d=7. Sharply jump appears at the inflection point. Hysteresis appears at the region close to the inflection point. Figure.22 shows the inflection point appears at about x/l=0.14. 21
Figure.21 Measured reaction force with extended sheet on top of the control volume versus x/d Figure.22 Measured reaction force with extended sheet on top of the control volume versus x/l 5.3.2 Visualization According to the visualization(see figure.23), even with the extended 22
sheet of 3 meters diameter. There is still air from the top of the sheet that entrain to the control volume. Figure.23 Air entrainment at the edge of the extended sheet for Re=55300, q=16.7 l/s, x/d=30, L=3 m Figure.24 shows a rotating flow that result from the conflict between the entraining air and spreading wall jet. But when the ceiling was lowered below the distance of inflection point, a weak rotating flow with opposite direction was observed. It is due to the air entrainment from the space between ground and the impinging plate to the spreading wall jet. (See figure.25) 23
Figure.24 A rotating flow at the edge of the extended sheet for Re=55300, q=16.7 l/s, x/d=30, L=3 m Figure.25 A rotating flow at the edge of the extended sheet for Re=55300, q=16.7 l/s, x/d=15, L=3 m 24
5.4 General Comparing Figure.26 and figure.27 show the comparison for air flows rate at 16.7 l/s and 18.7 l/s with different set-up. The inflection point obviously moves forward by installing the round sheet. But the inflection point appears at almost the same place before and after extending the round sheet. It shows, after a certain length, the length of the sheet will not affect the height where ambient air stops entrain the control volume. The results profiles before and after extending the round sheet are similar. It is due to the same magnitude influence from the ambient air, and one can predict the same magnitude influence in case of ventilation application with infinite or totally closed ceiling. The sharply jump at the inflection point only appears after installing the round sheet, and the value of the normalized force is less than 1 before the jump. It shows that the outward spreading jet can stop the entrainment more effectively than the rotating flow that could result in vortexes when L=0.12 m. Note that, with the smallest ceiling, slightly higher normalized reaction force can be achieved. It also shows that the direct air entrainment from the ceiling level of the control volume does not contribute much for the increase of the jet s momentum. 25
1.5 1 0.5 16.7 l/s Re=57000 L=0.12 m R/M(0) 0 0 10 20 30 40 50 16.7 l/s Re=57400 L=2 m -0.5-1 16.7 l/s Re=55300 L=3 m -1.5 x/d Figure.26 Measured reaction force with different Set-Ups versus x/d for q=16.7 l/s 1.5 1 18.7 l/s Re=64300 L=0.12 m R/M(0) 0.5 0-0.5 0 10 20 30 40 50 18.7 l/s Re=64800 L=2 m 18.7 l/s Re=62100 L=3 m -1 x/d Figure.27 Measured reaction force with different Set-Ups versus x/d for q=18.7 l/s Figure.28 and figure.29 show the comparison by weighting results. It is more obviously to note the decreasing of the weighting results by installing bigger ceiling and decreasing Reynolds number. 26
100 Weight (g) 80 60 40 20 0-20 -40-60 0 10 20 30 40 50 1.67 l/s Re=57000 L=0.12 m 1.67 l/s Re=57400 L=2 m 1.67 l/s Re=62100 L=3 m -80-100 x/d Figure.28 Measured weights with different Set-Ups versus x/d for q=16.7 l/s 150 100 1.87 l/s Re=64300 L=0.12 m Weight (g) 50 0-50 0 10 20 30 40 50 18.7 l/s Re=64800 L=2 m 18.7 l/s Re=62100 L=3 m -100 x/d Figure.29 Measured weights with different Set-Ups versus x/d for q=18.7 l/s When RR/MM (0) > 1, the air entrainment from ambient contributes to the increase of the momentum flux. Figure.30 shows the percentages of the contribution from air entrainment. With the smallest ceiling, the highest amount of entrainment is about 0.2. 27
(R-M(0))/R 0.25 0.2 0.15 0.1 0.05 0 0 10 20 30 40 50 x/d q=16.7 l/s Re=57070 L=0.12 m q=18.7 l/s Re=64300 L=0.12 m q=16.7 l/s Re=57400 L=2 m q=18.7 l/s Re=64800 L=2 m q=16.7 l/s Re=55300 L=3 m q=18.7 l/s Re=62100 L=3 m Figure.30 Percentages of the contribution from air entrainment 28
6. Conclusion The change in momentum flux for a jet have been weighted up to the distance x/d= 46. At a given Reynolds number, the momentum flux increases with the impinging distance up to about x/d=30, the air entrainment from ambient could increase the momentum flux of the jet up to 20% approximately. But the direct air entrainment from ceiling level does not contribute much for the momentum increase. From visualization, note that the air entrainment from the top of the control volume still happened at the ceiling edge even the extended round sheet was installed. The profile of the normalized reaction force has an inflection point, which is also approximately represented the distance where ambient air stopping directly entraining the control volume. 29
30
7. Reference [1] S. Maurel, C. Solliec. A turbulent plane jet impinging nearby and far from a flat plate. Experiments in Fluids, 2001 [2] Carlo Carcasci. An experimental investigation on air impinging jets using visualisation methods. Int. J. Thermal Sciences,1999 [3] Karimipanah T, Sandberg M. Confinement Effect of Jets in Ventilated Rooms. ASHRAE. [4] Heskestad, G., Hot wire measurements in a plane turbulent jet. ASME Journal of App. Mech., Vol, 32, P. 721, 1965. [5] Brandbury, L.J.S., The Structure of a self-preserving turbulent plane jet, J. Fluid Mech., Vol. 23, p. 31, 1965. [6] Goldschmidt, V.W., Young, M.F., Energy Spectrum and Turbulent Scales in a Plane Air Jet, Proceedings of 4 th Biennial symposium on Turbulent in Liquids, Rolla-Missori, p,39, 1975. [7] Gutmark, E., and Wygnanski, I., 1976, The planar turbulent Jet, J. Fulid Mech., Vol.73,Part 4, p.465, 1976. [8] Kotsovinos, N.E., A study of the Entrainment and Turbulent in a plane Buoyant Jet, W.M. Keck Lab., Hydraul. Water Resourses, CIT Rep. KHR-32, 1975. [9] Hussain, A.K.M.F., and Clark, A.R., Upstream influence on the Near Field of a Plane Turbulent Jet, Phys. Fluids, Vol. 20., No. 9, pp 1416-1426, 1977 [10] Panchapakesan, N.R., and Lumley, J.L., Turbulence Measurements in Axisymmetric Jets of Air and Helium. Part 1. Air Jet, J.Fluid Mech., Vol. 246, pp 197-223, 1993. [11] Hussein, H.j., Capp, S.P., and George, W.K., Velocity Measurements in a high Reynolds Number, Momentum-Conserving, Axisymmetric, Turbulent Jet, J. Fluid Mech., Vol. 258, pp.31-75, 1994. 31
[12] Kotsovinos, N.E., and Angelidis, P.B., The Momentum and Turbulent in a plane buoyant Jets, J.Fluid Mech., vol.299, pp. 453-470, 1991. [13] Karimipanah T, Sandberg M. Decay of momentum and velocity in an axisymmetric impinging jet. Proceedings of Roomvent 94 Krakow Polen 1994 Volume page 399-410 [14] Karimipanah T, Sandberg M. Momentum balance of an axisymmetric jet and its impingement on a flat surface. Fluid Engineering [15] Karimipanah, T. Turbulent jets in confined saces. Thesis. Center for Built Environment, KTH, Gavle, Sweden, ISBN 91-7170-667-4, 1996 [16] Mateo de guadalfajara Pinilla. Jet kinetic momentum flux in an open ambient recorded with a digital balance, entrain of air from the surrounding. University of Gavle 32
Appendix 33
34
Ⅰ.Measurement and calculation details for Set-up with L=0.12m at air flow 16.7 l/s. 35
36
Ⅱ.Measurement and calculation details for Set-up with L=0.12 m at air flow rate 18.7 l/s. 37
38
Ⅲ.Measurement and calculation details for Set-up with L=2 m at air flow 16.7 l/s. 39
40
Ⅳ.Measurement and calculation details for Set-up with L=2 m at air flow rate 18.7 l/s. 41
42
Ⅴ.Measurement and calculation details for Set-up with L=3 m at air flow rate 16.7 l/s. 43
44
Ⅵ.Measurement and calculation details for Set-up with L=3 m at air flow rate 18.7 l/s. 45