Simulation of and Experimental Research on Rivulet Model on Airfoil Surface
Abstract
:1. Introduction
2. Mathematical Models
2.1. Water Film Flow Model
2.2. Rivulet Flow Model
2.2.1. Two-Dimensional Water Film Breakup Model
2.2.2. Three-Dimensional Mass Flow Allocation Model
2.3. Roughness Model
3. Experimental System
4. Results and Analysis
4.1. Validation of the Simulation Model
- When the flow rate of the water supply from the nozzle is constant, an increase in the air speed delays the breakup of the water film. This is because, when the air speed increases, the shear force on the surface of the wing also increases, resulting in a decrease in the critical thickness of the water film breakup.
- The width and spacing of the rivulets decrease as the air speed increases. However, the velocity of the rivulet flow gradually increases. According to Equations (9) and (23)–(26), the rivulet flow velocity is proportional to τ1/3, and the spacing of the rivulets is proportional to τ−2/3. Thus, the shear stress increases with an increase in air speed, resulting in an increasing rivulet velocity and a decreasing rivulet spacing.
4.2. Effect of Surface Temperature on Water Film Breakup
4.3. Effect of Roughness on Water Film Breakup
- When the surface of the wing has a uniform roughness, an increase in roughness causes the water film to break earlier. This is because, according to Equation (48), an increase in roughness results in a smaller surface shear than that of a smooth surface. When Equation (23) is combined, it is found that the dimensionless critical thickness becomes smaller; therefore, the water film breaks up earlier.
- The thickness, width, and spacing of the rivulets increase with increasing surface roughness. This is because, as the surface roughness increases, the surface shear stress decreases. Combined with the experimental phenomenon, an increase in surface roughness increases the resistance to water flow, which is manifested by the water flowing more slowly and the width and spacing of rivulets becoming larger. From a theoretical point of view, the same conclusion can be drawn from Equations (10) and (25)–(27).
- A comparison of the last two rows in Table 3 shows that the water film breakup location does not change when the roughness increases to a specific value. This is because the water film breakup location moves forward and coincides with the water impingement limit. Since the experiments showed that the water film does not break up within the impingement range of the water droplets, the breakup location does not move forward, even if the roughness increases further. Therefore, it can be concluded that when the water flow rate is relatively high, water film breakup should occur outside of the impingement limit. However, there is a critical value of roughness such that, when the roughness is greater than or equal to it, water film breakup occurs at the impingement limit.
- As the air speed increases, the droplet impingement zone increases, i.e., the impingement limit shifts back. Therefore, it can also be concluded that, as the air speed increases, the critical value of the roughness also increases.
4.4. Effect of Water Film Breakup on Anti-Ice System
5. Conclusions
- At a given water supply flow rate, an increase in air speed delays the breakup of the water film. Meanwhile, the width and spacing of the rivulets decrease, and the rivulet flow velocity increases.
- An increase in surface temperature advances the breakup of the water film. The thickness, spacing, width, and flow velocity of the rivulets decrease with an increasing surface temperature.
- If uniform roughness exists on the surface, the increased roughness causes the water film to break earlier. The thickness, spacing, and width of the rivulets increase, and the rivulet flow velocity decreases with an increasing surface roughness.
- When the water flow rate is relatively large, there is a critical value of roughness, causing water film breakup to occur at the water impingement limit when the roughness is greater than or equal to this value. At this point, the breakup location no longer changes. As the air speed increases, the impingement limit moves downstream; therefore, the critical value of the roughness increases.
Author Contributions
Funding
Data Availability Statement
Conflicts of Interest
References
- Gohardani, O. Progress in Aircraft Icing and Aircraft Erosion Research; Aerospace and System Engineering and Research; Nova Science Publishers: Hauppauge, NY, USA, 2017; ISBN 978-1-5361-2030-1. [Google Scholar]
- Thomas, S.K.; Cassoni, R.P.; MacArthur, C.D. Aircraft Anti-Icing and de-Icing Techniques and Modeling. J. Aircr. 1996, 33, 841–854. [Google Scholar] [CrossRef]
- Carroll, T.; McAvoy, W. Formation of Ice on Airplanes. Airw. Age 1928, 58–59. [Google Scholar]
- Alekseenko, S.; Nakoryakov, V.E.; Pokusaev, B.G. Wave Flow of Liquid Films; Begell House: Danbury, CT, USA, 1994. [Google Scholar]
- Johnson, M.F.G.; Schluter, R.A.; Miksis, M.J.; Bankoff, S.G. Experimental Study of Rivulet Formation on an Inclined Plate by Fluorescent Imaging. J. Fluid Mech. 1999, 394, 339–354. [Google Scholar] [CrossRef]
- Hoffmann, A.; Ausner, I.; Repke, J.-U.; Wozny, G. Fluid Dynamics in Multiphase Distillation Processes in Packed Towers. Comput. Chem. Eng. 2005, 29, 1433–1437. [Google Scholar] [CrossRef]
- Hagemeier, T.; Hartmann, M.; Kühle, M.; Thévenin, D.; Zähringer, K. Experimental Characterization of Thin Films, Droplets and Rivulets Using LED Fluorescence. Exp. Fluids 2012, 52, 361–374. [Google Scholar] [CrossRef]
- Zhang, K.; Hu, H. An Experimental Study of the Wind-Driven Water Droplet/Rivulet Flows over an Airfoil Pertinent to Wind Turbine Icing Phenomena. In Proceedings of the Fluids Engineering Division Summer Meeting, Chicago, IL, USA, 3–7 August 2014; Volume 46247, p. V01DT39A001. [Google Scholar]
- Bonart, H.; Marek, A.; Repke, J.-U. Experimental Characterization of Stable Liquid Rivulets on Inclined Surfaces: Influence of Surface Tension, Viscosity and Inclination Angle on the Interfacial Area. Chem. Eng. Res. Des. 2017, 125, 226–232. [Google Scholar] [CrossRef]
- Liu, Y.; Zhang, K.; Tian, W.; Hu, H. An Experimental Study to Characterize the Effects of Initial Ice Roughness on the Wind-Driven Water Runback over an Airfoil Surface. Int. J. Multiph. Flow 2020, 126, 103254. [Google Scholar] [CrossRef]
- Messinger, B.L. Equilibrium Temperature of an Unheated Icing Surface as a Function of Air Speed. J. Aeronaut. Sci. 1953, 20, 29–42. [Google Scholar] [CrossRef]
- Hobler, T. Heat and Mass Transfer Bibliography: Polish Works. Int. J. Heat Mass Transf. 1965, 8, 841–844. [Google Scholar] [CrossRef]
- Mikielewicz, J.; Moszynski, J.R. Breakdown of a Shear Driven Liquid Film; Polish Academy of Sciences: Gdansk, Poland, 1975. [Google Scholar]
- Al-Khalil, K.M. Numerical Simulation of an Aircraft Anti-Icing System Incorporating a Rivulet Model for the Runback Water; The University of Toledo: Toledo, Spain, 1991. [Google Scholar]
- Rothmayer, A.; Tsao, J. Water Film Runback on an Airfoil Surface. In Proceedings of the 38th Aerospace Sciences Meeting and Exhibit, Reno, NV, USA, 10–13 January 2000; p. 237. [Google Scholar]
- Saber, H.H.; El-Genk, M.S. On the Breakup of a Thin Liquid Film Subject to Interfacial Shear. J. Fluid Mech. 2004, 500, 113–133. [Google Scholar] [CrossRef]
- Silva, G.; Silvares, O.; Zerbini, E. Water Film Breakdown and Rivulets Formation Effects on Thermal Anti-Ice Operation Simulation. In Proceedings of the 9th AIAA/ASME Joint Thermophysics and Heat Transfer Conference, San Francisco, CA, USA, 5–8 June 2006; p. 3785. [Google Scholar]
- Silva, G.; Silvares, O.; Zerbini, E.; Hefazi, H.; Chen, H.-H.; Kaups, K. Differential Boundary-Layer Analysis and Runback Water Flow Model Applied to Flow around Airfoils with Thermal Anti-Ice. In Proceedings of the 1st AIAA Atmospheric and Space Environments Conference, San Antonio, TX, USA, 22–25 June 2009; p. 3967. [Google Scholar]
- Wang, G.; Rothmayer, A. Air Driven Water Flow Past Small Scale Roughness. In Proceedings of the 43rd AIAA Aerospace Sciences Meeting and Exhibit, Reno, NV, USA, 10–13 January 2005; p. 653. [Google Scholar]
- Wang, G. Thin Water Films Driven by Air through Surface Roughness; Iowa State University: Ames, IA, USA, 2008. [Google Scholar]
- Dong, W.; Zheng, M.; Zhu, J.; Lei, G. Calculation and Analysis of Runback Water Flow on Anti-Icing Airfoil Surface. J. Aircr. 2016, 53, 1597–1605. [Google Scholar] [CrossRef]
- Myers, T.G.; Thompson, C.P. Modeling the Flow of Water on Aircraft in Icing Conditions. AIAA J. 1998, 36, 1010–1013. [Google Scholar] [CrossRef]
- Zhao, H.; Zhu, C.; Zhao, N.; Zhu, C.; Wang, Z. Experimental Research on the Influence of Roughness on Water Film Flow. Aerospace 2021, 8, 225. [Google Scholar] [CrossRef]
- Bu, X.; Lin, G.; Shen, X.; Hu, Z.; Wen, D. Numerical Simulation of Aircraft Thermal Anti-Icing System Based on a Tight-Coupling Method. Int. J. Heat Mass Transf. 2020, 148, 119061. [Google Scholar] [CrossRef]
- LI, H.; HUANG, P. Numerical Simulation of Mixed Phase Icing on Two-Dimensional Airfoil. Acta Aeronaut. Astronaut. Sin. 2021, 42, 124085. [Google Scholar]
- Shen, X.; Wang, H.; Lin, G.; Bu, X.; Wen, D. Unsteady Simulation of Aircraft Electro-Thermal Deicing Process with Temperature-Based Method. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2020, 234, 388–400. [Google Scholar] [CrossRef]
Contact Angle (deg) | A (mm) | λ (mm) | |
---|---|---|---|
Simulation | 2 | 2.75 | 2.51 |
4 | 1.57 | 1.89 | |
16 | 0.41 | 0.55 | |
32 | 0.20 | 0.27 | |
64 | 0.09 | 0.13 | |
Experiment | - | 1.74 ± 0.30 | 3.22 ± 0.30 |
Surface Temperature (K) | A (mm) | λ (mm) | Thickness at the Breakup Location (10−5 m) | Breakup Location x/c (%) | Velocity (mm/s) |
---|---|---|---|---|---|
275.15 | 2.46 | 2.88 | 3.14 | 16.93 | 113.13 |
278.15 | 2.45 | 2.85 | 3.11 | 16.03 | 112.48 |
283.15 | 2.41 | 2.81 | 3.07 | 15.15 | 110.82 |
293.15 | 2.30 | 2.71 | 2.95 | 12.71 | 105.82 |
Rz (μm) | A (mm) | λ (mm) | Thickness at the Breakup Location (10−5 m) | Breakup Location x/c (%) | Velocity (mm/s) |
---|---|---|---|---|---|
150 | 4.67 | 5.44 | 5.93 | 42.27 | 107.19 |
200 | 4.68 | 5.46 | 5.96 | 20.91 | 86.69 |
250 | 5.82 | 6.56 | 7.16 | 12.71 | 68.09 |
300 | 11.26 | 11.11 | 12.19 | 12.71 | 45.87 |
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Lou, Y.; Bu, X.; Shen, X.; Lin, G.; Zhang, R.; Zeng, F.; Jin, H.; Ma, K.; Wen, D. Simulation of and Experimental Research on Rivulet Model on Airfoil Surface. Aerospace 2022, 9, 570. https://doi.org/10.3390/aerospace9100570
Lou Y, Bu X, Shen X, Lin G, Zhang R, Zeng F, Jin H, Ma K, Wen D. Simulation of and Experimental Research on Rivulet Model on Airfoil Surface. Aerospace. 2022; 9(10):570. https://doi.org/10.3390/aerospace9100570
Chicago/Turabian StyleLou, Yanxia, Xueqin Bu, Xiaobin Shen, Guiping Lin, Ruchen Zhang, Feixiong Zeng, Haichuan Jin, Kuiyuan Ma, and Dongsheng Wen. 2022. "Simulation of and Experimental Research on Rivulet Model on Airfoil Surface" Aerospace 9, no. 10: 570. https://doi.org/10.3390/aerospace9100570
APA StyleLou, Y., Bu, X., Shen, X., Lin, G., Zhang, R., Zeng, F., Jin, H., Ma, K., & Wen, D. (2022). Simulation of and Experimental Research on Rivulet Model on Airfoil Surface. Aerospace, 9(10), 570. https://doi.org/10.3390/aerospace9100570