Next Article in Journal
Three-Dimensional Analytical Solutions for Acoustic Transverse Modes in a Cylindrical Duct with Axial Temperature Gradient and Non-Zero Mach Number
Next Article in Special Issue
Numerical and Experimental Research on Flight Control of a V-Tail Configuration for the Wind Tunnel Model of Aircraft
Previous Article in Journal
Numerical Investigation on Hybrid Rocket Motors with Star-Segmented Rotation Grain
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Investigation on the Effect of Blockage on the Icing of Airfoils

1
Key Laboratory of Unsteady Aerodynamics and Flow Control, Ministry of Industry and Information Technology, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
2
Shanghai Aircraft Design and Research Institute, Shanghai 201210, China
*
Author to whom correspondence should be addressed.
Aerospace 2022, 9(10), 587; https://doi.org/10.3390/aerospace9100587
Submission received: 6 September 2022 / Revised: 5 October 2022 / Accepted: 7 October 2022 / Published: 10 October 2022
(This article belongs to the Special Issue Thermal Fluid Dynamics and Control in Aerospace)

Abstract

:
The blockage is one of the important factors affecting the icing of airfoils in wind tunnel tests. In this paper, numerical simulations are conducted to study the effect of blockage on the icing of different airfoils. By reducing the height of testing wind tunnels, the blockage is increased, and the changes in the height and angle of the ice horn are numerically investigated. The simulation results indicate that as the blockage increases, the flow velocity above the stagnation point of the airfoil increases, leading to larger pressure coefficients distribution and stronger heat transfer capacity. As a result, the position of icing moves forward, and the angle of the upper ice horn becomes smaller. In addition, the increased flow velocity facilitates the collection of water droplets in the area, which improves the icing and increases the height of the upper ice horn. It is also found that the blockage increases the angle of attack of the airfoil, moving the stagnation point backward and decreasing the angle of the upper ice horn. When the blockage is above 15%, the joint influence of the opening angle and height of the upper ice horn significantly reduces the projection height of the upper ice horn in the direction of the incoming flow, leading to unacceptable criticality of the ice shape.

1. Introduction

Ice accretion occurs immediately on wing, tail, propeller, rotor and even antenna when aircraft fly through the cloud full of supercooled water droplets or encountering precipitations such as freezing rain or drizzle. Ice accumulating on the aerodynamic sensitive surfaces leads safety issue because of its detrimental effect on the performance such as maximum lift penalty, stall angle reduction, and parasite drag increase. Aviation accident data showed that 583 icing accidents occurred and caused more than 800 fatalities from 1982 to 2000 in the United States [1]. For this reason, aircraft icing was even recognized as the “most wanted aviation transportation safety improvement” by the National Transportation Safety Board (NTSB) [2]. The Federal Aviation Regulations (FAR) Part 25 Appendix C define icing envelopes (in terms of air temperature, liquid water content and droplets’ size) for aircraft certification, corresponding to 99.9% of icing conditions found in stratiform clouds [3]. Most light aircraft are not required to pass this certification and do not usually have the required ice-protection systems for flying into icing conditions. Larger aircraft are equipped with a variety of anti-icing or de-icing systems (heaters, pneumatic boots and liquid flows) to help them prevent ice formation [4].
In order to guarantee safe and on schedule operations of transportation aircraft, accurate prediction of ice accretion, related performance degradation and anti/de-icing systems development play an important role for the design of aircraft. Over the years, test approaches including flight and wind-tunnel tests, and numerical simulations have been adopted in the prediction of ice accretion. Following the development and progress of computational science as well as the great advantages of numerical simulation in the economic cost, many ice-prediction codes have been successfully developed and serve as valuable design and certification tools for the aircraft manufacturers such as LEWICE developed by NASA Lewis Research Center [5], ONERA 2-D/3-D developed by ONERA [6], IMPIN3D developed by Italy [7], and the second-generation 3-D icing simulation system FENSAP-ICE developed by Numerical Technologies International [8], etc. Recently, a large number of articles have been published that are devoted to the application of machine learning methods, deep neural networks in the field of fluid dynamics and even the procedure and method for the ice accretion prediction for different airfoils using artificial neural networks (ANNs) are discussed [9]. It is worth expecting that the combination of two approaches: machine learning and numerical simulation, will speed up the prediction of the icing shape and significantly reduce computational costs.
Considering that the natural icing conditions are very difficult and costly to obtain, the icing test in wind tunnels is not only an important means to study aircraft icing (an essential step to obtaining civil aircraft airworthiness certification) [10,11] but also an important method to develop and verify numerical tools used for icing [11,12,13,14]. The civil aircraft authorities have set clear approval requirements for the wind tunnels and numerical tools used for icing in the design of aircraft [15,16]. Previous wind tunnel tests show that wall interference is an important factor affecting the accuracy of testing results [17,18,19]. In a closed wind tunnel, the existence of walls limits the bending and diffusion of streamlines. Therefore, the area through which the airflow between adjacent streamlines is smaller than that in the free space, and the average flow velocity between the airfoil and tunnel wall is larger. This phenomenon is called the blocking effect. In reality, the limitation of the bending of streamlines is equivalent to inducing an upward angle to the airfoil, thereby increasing the angle of attack as compared with that defined using the axis of the wind tunnel [18]. The increased angle of attack will move the icing from the upper surface of the airfoil to the lower surface, leading to overestimated aerodynamic performances of airfoils. Consequently, the criticality of the ice shape is not guaranteed, and the test results will not be recognized.
The blockage is defined as the ratio of the maximum incoming airflow area to the cross-sectional area of the test tunnel. The US Federal Aviation Administration [20] reported that the blockage should not exceed 10% in the icing test of airfoils in wind tunnels. However, in practice, the chord length of the wings of large passenger aircraft is generally long (may exceed 5 m). Even the world’s largest wind tunnel (the FL-16 ice wind tunnel [21] of the China Aerodynamic Research and Development Center, with a size of 4.8 m × 3.2 m) is not suitable for such experiments. On the other hand, the aircraft airworthiness examiners generally do not allow the use of scaled test models in icing wind tunnel tests. Even though the hybrid wing design technology is adopted [22,23], the blockage will still exceed 10%, which makes the reliability of icing test results questionable under certain circumstances (e.g., with a large angle of attack).
Guo Qiling et al. [24,25] of the Key Laboratory of Icing and Deicing of CARDC studied the influence of walls on the impinging characteristics of water droplets and icing of airfoils. Their research shows that the blockage will affect the impinging characteristics of water droplets, and the increased blockage will improve icing on the surface of airfoils. Zocca et al. [26] found that the blockage has a significant impact on the icing of airfoils and brings a large deviation to the ice shape. Qin et al. [27] numerically examined the influence of cavity walls on the impinging characteristics of water droplets. They concluded that the blockage would shift the stagnation point of airflows on the airfoil, which significantly influences the collection of water droplets on the airfoil surface. The NRC wind tunnel in Canada also noticed that the increased blockage is beneficial for the collection of water droplets. They revealed the relationships between the water droplet collection rate and the Reynolds number of water droplets at different blockages [28]. It can be seen from the studies mentioned above that the impinging characteristics of water droplets and resultant ice shape are both affected by the blockage. Nonetheless, there is still no recommended value of the maximum blockage acceptable in wind tunnel tests. In this paper, the NACA0012 and GLC305 airfoils in the literature are used as examples and numerically investigated. By adjusting the height of the testing wind tunnel, the icing characteristics of airfoils at different blockages are simulated. Based on the ice shapes obtained from the tests, the influence of the blockage on the ice shape is studied, and the maximum blockage allowed in the icing tests of airfoils is explored. The outline of this article is organized as follows. Section 2 introduces the numerical methodology employed in this research. Section 3 validates the numerical model. Section 4 presents and discusses the simulation results. Finally, concluding remarks are drawn in Section 5.

2. Numerical Methodology

2.1. Model Description and Mesh Configuration

Figure 1a presents the sketch of the airfoil in a wind tunnel under investigation in this study. In the numerical simulations, for simplicity, the airfoil is installed without fixed boundary conditions, and its left and right supporting structures are omitted. The blockage is changed by adjusting the height between the upper and lower tunnel walls, ignoring the interference of the side tunnel walls. For instance, If the height of the teste section is decrease and the installation angle of the wing model remains the same, the blockage is increased. As shown in Figure 1a, two-dimensional unstructured grids are adopted to mesh the airfoil and the surrounding flow field. To save computational cost, denser grids are constructed in the airfoil region as compared with the surrounding flow field that employs relatively coarse grids. As such, triangular prism boundary layer grids are used on the surface of the airfoil, as displayed in Figure 1b. Likewise, triangular pyramid grids are employed near the leading edge of the airfoil, as shown in Figure 1c. The prism layers are 40, the first layer is about 1 × 10−6 m high, and the expansion ratio is 1.2, the number of grid cells along the airfoil is 120.

2.2. Governing Equations

Following the construction of mesh for the airfoil in the wind tunnel, numerical calculations are conducted to simulate the ice accretion. The working fluid is air and follows the Navier-Stokes equations and the equation of state, which are,
ρ ¯ t + x i ( ρ ¯ u ˜ i ) = 0
t ( ρ ¯ u ˜ i ) + x j ( ρ ¯ u ˜ i u ˜ j + p ¯ δ i j τ i j t o t ˜ ) = 0
t ( ρ ¯ e ˜ t o t ) + x j ( ρ ¯ u ˜ j e ˜ t o t + u ˜ j p ¯ + q j t o t ˜ u ˜ i τ i j t o t ˜ ) = 0
p = ρ R g T
where the overbars “~” and “−” denote the density-weighted time averaging (Favre averaging) and classical time averaging (Reynolds averaging), respectively. In this case, p, ρ and T are pressure, density and temperature, ui or uj stands for the velocity component in corresponding direction, etot = e + uiui/2 is the total energy with e being the specific internal energy, τ i j t o t is the total stress tensor, q j t o t is the heat flux, δij is the Kronecker delta, and Rg is the specific gas constant. The Spalart–Allmaras turbulence model is used in the flow field calculation. In the calculation of droplet trajectory, a single average diameter called Median Volumetric Diameter (MVD) is used, and the coupled heat transfer calculation related to the anti-ice system is not involved.
The governing equations stated above are resolved numerically using the commercial CFD (computational fluid dynamics) package ANSYS Fensap-ice R19.0. Fensap-ice R19.0 includes many modules, such as flow field calculation, droplet collision calculation, ice accretion, and heat transfer calculation, which can be used for icing analysis.

2.3. Boundary Conditions

The boundary conditions in the numerical simulations are defined referring to the NACA0012 and GLC305 airfoils. Table 1 lists the parameters of NACA0012 and GLC305 airfoils used in the icing test in NASA Lewis icing Research Tunnel (IRT) [13].

2.4. Sensitivity Studies

Sensitivity studies on the influence of grid number and icing step number on the simulation results are performed. In the multistep icing calculation, the multi-step method based on artificial mesh reconstruction is adopted to ensure that the reconstructed mesh can accurately describe the ice shape and maintain sufficient mesh quality. The time interval of each step of mesh reconstruction is equal: the icing time of each step is defined as the total icing time divided by mesh reconstruction steps.
Figure 2 shows the ice shapes calculated at different grid numbers and mesh reconstruction steps. The “160k 6steps” in the figure represents the case in which the total grid number is 160,000 and the icing step number is 6. First, by comparing the ice shapes at grid numbers of 160,000, 300,000 and 476,000, It has been found that the grid number of 476,000 is sufficient enough to accurately describe the ice shape while minimizing the computational cost simultaneously. Therefore, the optimal grid number is chosen at 476,000. Subsequently, by comparing the ice shapes at icing step numbers of 3, 6 and 7 when the grid number is 476,000, It has been noticed that the mesh reconstruction time has a great impact on the ice horn height. When the icing step number reaches 6, the ice shape becomes less sensitive. Hence, 6 steps of icing are adopted in the subsequent simulations in this paper.

3. Verification of Numerical Results

Comparisons are made between the experimental data of the NACA0012 airfoil in the literature [13] and the ice shape calculated by ANSYS Fensap-ice R19.0 using the same parameters. The results in Figure 3 indicate that the criticality of ice shape from numerical simulations in terms of ice horn height and angle is consistent with that of the experiment, which demonstrates the validity of the numerical methodologies proposed in this research.

4. Results and Discussion

4.1. Criticality of Ice Shape

Gray [29] demonstrated that the height, position, and angle of the upper ice horn are the critical parameters that affect the formation of ice shape which in return influences the aerodynamic performance of the airfoil. The angle of the upper ice horn has a direct impact on the size and position of bubbles separating on the surface of the airfoil while the parameters of the lower ice horn generally affect the range of negative angle of attack, having little influence on the aerodynamic characteristics of the airfoil. Figure 4a displays the height of the upper ice horn h, the chord length of the upper ice horn c, and the opening angle of the upper ice horn θ. The characteristic position x of the ice shape is determined by the intersection of the trailing edge of the ice horn and the surface of the airfoil. The characteristic height h is the distance between the intersection described above and the highest point of the ice horn. In practical work, the projection height H of the upper ice horn (see Figure 4b) in the incoming flow direction is generally used to judge the criticality of the ice shape. Specifically, the larger H is, the more critical the ice shape is [30]. The parameter H reflects the joint influence of h, x, and θ.

4.2. Effect of Blockage on the Criticality of Ice Shape

To simulate the installation of supports (left and right) of the airfoil in practical work, the upper and lower tunnel walls are moved toward the airfoil, thus increasing the blockage. Three blockages (10%, 15.2%, and 20%) are investigated in the simulations and the pressure coefficient distribution of each blockage is calculated. The results show that when the angle of attack is unchanged, as the blockage increases, the peak of the pressure coefficient distribution increases (see Figure 5a), and the stagnation point moves backward (see Figure 5b). This phenomenon is also reported in reference [17] which concludes that the average flow velocity between the airfoil and the tunnel wall increases after the blockage increases, leading to more restrictions on the streamline bending by the tunnel wall. Thus, the angle of attack increases because of the increase of the angle (upward) of the incoming flow.
In the simulations, the ice shapes at different blockages are calculated and presented in Figure 6. By comparing the calculated ice shape without a fixed boundary with that measured in the experiment, it has been found that as the blockage increases, the height of the upper ice horn increases, and the opening angle of the upper ice horn decreases. The increase in the height of the upper ice corner is due to the increase of local velocity around the head of the airfoil, thereby facilitating the collection of water droplets.
The local droplet collection efficiency β is defined as the ratio, for a given mass of water, of the area of impingement to the area through which the water passes at some distance upstream of the airfoil. Taking a unit width as one dimension of both area terms, the local droplet collection efficiency can then be defined as,
β = d y d s = Δ y 0 Δ s
where Δ y 0 is the spacing between water droplets at the release plane and Δ s is the distance along the airfoil surface between the impact locations of the same two droplets.
As can be seen from Figure 7, the local droplet collection coefficient β increases with the increase of blockage, and the peak of β gradually moves towards the leading edge. Meanwhile, due to the increase of flow velocity in this area, the heat transfer rate is enhanced which results in faster icing of water film, smaller flow distance, and reduced angle of the upper ice horn. Figure 6b shows that when the blockage reaches 20% or even 30%, the change in the opening angle of the upper ice horn becomes more obvious.
Numerical simulations are also performed on the airfoil of GLC305 reported in the literature [13]. As can be seen in Figure 8, as the blockage increases, the opening angle of the upper ice horn decreases. In particular, when the blockage increases to 17.6%, there is an obvious decline in the opening angle of the upper ice horn. Figure 9 shows the relationship between the blockage and the dimensionless projection height H/c of the upper ice horn in the direction of incoming flow, where c represents the chord length of the airfoil. It is found that H/c is small when the blockage falls between 10% and 15%. When the blockage is greater than 15%, H/c decreases remarkably, which significantly affects the criticality of the ice shape.

5. Conclusions

This study investigates the effect of blockage on the icing of airfoils through computational fluid dynamics. Simulations are performed on the airfoils of NACA0012 and GLC305 tested in the wind tunnels, and comparisons are made between the calculated results and experimental data in the literature. The criticality of ice shape from numerical simulations in terms of ice horn height and angle is consistent with that of the experiment, which demonstrates the validity of the numerical methodologies proposed in this study. The key conclusions of this research are summarized as follows.
(1)
The increase of ice horn height will increase the criticality of ice shape, while the decrease of ice horn angle will reduce the criticality of ice shape
(2)
As the blockage increases, the peak of the pressure coefficient distribution increases, and the stagnation point moves backward. This is because the average flow velocity between the airfoil and the tunnel wall increases after the blockage increases, leading to more restrictions of the streamline bending by the tunnel wall.
(3)
With the increase of the blockage of the ice tunnel, the joint influence of the opening angle and height of the upper ice horn significantly reduces the projection height of the upper ice horn in the direction of the incoming flow. As a result, the criticality of the ice shape is reduced. To ensure the criticality of the ice shape, the blockage should be below 15%.

Author Contributions

Conceptualization and methodology, Z.L. and D.L.; methodology and software Z.H.; validation, X.X. and Y.H.; formal analysis and investigation, D.L.; resources and data curation, Z.H.; writing—original draft preparation, D.L.; writing—review and editing, Z.L.; visualization, Z.H.; supervision, Z.L.; project administration, D.L.; funding acquisition, Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Fundation of China (grant nos. 12102187 and 11872212), the National Key R&D Program of China (Project No. 2020YFA0712000), and a project funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions.

Conflicts of Interest

The authors declare that they have no known competing financial interest or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Petty, K.R.; Floyd, C.D.J. A Statistical Review of Aviation Airframe Icing Accidents in the US. In Proceedings of the 11th Conference on Aviation, Range, and Aerospace, National Trans-portation Safety Board, Hyannis, MA, USA, 7 October 2004; Available online: https://ams.confex.com/ams/11aram22sls/techprogram/paper_81425.htm (accessed on 8 October 2022).
  2. Ratvasky, T.P.; Barnhart, B.P.; Lee, S. Current Methods Modeling and Simulating Icing Effects on Aircraft Performance, Stability, Control. J. Aircr. 2010, 47, 201–211. [Google Scholar] [CrossRef]
  3. Heinrich, A.; Ross, R.; Zumwalt, G. Aircraft Icing Handbook; Civil Aviation Authority: Lower Hutt, New Zealand, 2000; pp. 13–14. [Google Scholar]
  4. Shinkafi, A.; Lawson, C.; Seresinhe, R.; Quaglia, D.; Madani, I. An Intelligent Ice Protection System for Next Generation Aircraft Trajectory Optimisation. In Proceedings of the 29th Congress of the International Council of Aeronautical Sciences, ICAS 2014, St. Petersburg, Russia, 7–12 September 2014; Available online: https://www.researchgate.net/publication/290354048_An_intelligent_ice_protection_system_for_next_generation_aircraft_trajectory_optimisation (accessed on 10 March 2019).
  5. Wright, W.B. User’s Manual for LEWICE 3.2; NASA/CR-2008-214255; NASA: Washington, DC, USA, 2008.
  6. Hedde, T.; Guffond, D. Improvement of the ONERA 3D icing code, comparison with 3D experiment shapes. In Proceedings of the 31st Aerospace Sciences Meeting, Reno, NV, USA, 11–14 January 1993; AIAA 93-0169. [Google Scholar]
  7. Mingione, G.; Brandi, V.; Saporiti, A. A 3D ice accretion simulation code. In Proceedings of the 37th Aerospace Sciences Meeting and Exhibit, Reno, NV, USA, 11–14 January 1999; AIAA 99-0247. [Google Scholar]
  8. Habashi, W.G.; Morency, F.; Beaugendre, H. FENASP-ICE: A Second Generation 3D CFD-Based In-Flight Icing Simulation System; SAE Technical Paper 2003-01-2157, 2003. Available online: https://saemobilus.sae.org/content/2003-01-2157/ (accessed on 5 September 2022). [CrossRef]
  9. Strijhak, S.; Ryazanov, D.; Koshelev, K.; Ivanov, A. Neural Network Prediction for Ice Shapes on Airfoils Using iceFoam Simulations. Aerospace 2022, 9, 96. [Google Scholar] [CrossRef]
  10. Zhou, F.; Feng, L.; Xu, C.; Zhao, K.; Han, Z. Determination and verification of critical ice shape for the certification of civil aircraft. J. Exp. Fluid Mech. 2016, 30, 8–13. [Google Scholar]
  11. Zhan, P. Review on the system of icing facilities in NASA. Aeronaut. Sci. Technol. 2021, 32, 1–6. [Google Scholar]
  12. Hedde, T.; Guffond, D. Development of a three-dimensional icing code—Comparison with experimental shapes. In Proceedings of the 30th Aerospace Sciences Meeting and Exhibit, Reno, NV, USA, 6–9 January 1992. [Google Scholar]
  13. Weight, W.B.; Rutkowski, A. Validation Results for LEWICE2.0; NASA/CR-1999-208690; NASA: Washington, DC, USA, 1999.
  14. Weight, W.B. Validation Results for LEWICE 3.0. In Proceedings of the 43rd AIAA Aerospace Sciences Meeting and Exhibit, Reno, NV, USA, 10–13 January 2005. [Google Scholar]
  15. SAE ARP5903-2003; Droplet Impingement and Ice Accretion Computer Codes. SAE International: Warrendale, PA, USA, 2003. [CrossRef]
  16. SAE ARP5905-2003; Calibration and Acceptance of Icing Wind Tunnels. SAE International: Warrendale, PA, USA, 2005. [CrossRef]
  17. Li, Z. Wind Tunnel Test Manual; Aviation Industry Press: Beijing, China, 2015; pp. 853–854. [Google Scholar]
  18. Editorial Committee of military training materials of PLA General equipment department. Low Speed Wind Tunnel Test; National Defense Industry Press: Beijing, China, 2002; pp. 332–333. [Google Scholar]
  19. Barlow, J.B., Jr.; Rae, W.H.; Pope, A. Low-Speed Wind Tunnel Testing, 3rd ed.; Wiley: New York, NY, USA, 1999; pp. 234–366. [Google Scholar]
  20. Aircraft Ice Protection, AC20-73A; Federal Aviation Administration: Washington, DC, USA, 2006.
  21. Guo, X.; Zhang, P.; Zhao, X.; Yang, S.; Lin, W. Validation of thermal flow field in large icing wind tunnel. Exp. Fluid Dyn. 2020, 34, 82–91. [Google Scholar]
  22. Zhao, K.L.; Lu, Z.L.; Ding, L.; Tan, G.; Feng, L.; Guo, T. A design method of hybrid airfoil applied in icing wind tunnel test. Acta Aerodyn. Sin. 2013, 31, 718–722. [Google Scholar]
  23. Han, Z.R.; Zhao, K.L.; Yan, W. Application of gurney flap in hybrid wing design for icing wind tunnel. Acta Aeronaut. Et Astronaut. Sin. 2019, 40, 522387, In Chinese. [Google Scholar] [CrossRef]
  24. Guo, Q.; Liu, S.; Guo, X.; Liu, Y.; Yang, S. Influence of sidewall interference on the impact characteristics of water droplets on wing surface in icing wind tunnel. In Proceedings of the 11th National Conference on Hydrodynamics, Shenzhen, China, 3 December 2020. [Google Scholar]
  25. Guo, Q.; Liu, S.; Yi, X.; GUO, X.; Liu, Y. Influence of wind tunnel sidewall interference on NACA0012 airfoil icing growth. Key Laboratory of icing and anti icing, China Aerodynamics Research and development center. In Proceedings of the Fourth National Symposium on Icing and Anti Icing, Chongqing, China, 12–14 November 2020. [Google Scholar]
  26. Zocca, M.; Gori, G.; Guardone, A. Blockage and three-dimensional effects in wind-tunnel testing of ice accretion over wings. J. Aircr. 2016, 54, 1–9. [Google Scholar] [CrossRef]
  27. Qin, C.; Loth, E. Numerical Study of Droplet Trajectory and Collection Efficiency in IRT with Large Blockage Effects. In Proceedings of the 9th AIAA Atmospheric and Space Environments Conference, Denver, CO, USA, 5–9 June 2017. [Google Scholar]
  28. Catherine, C. Effects of Aerodynamic Blockage on Stagnation Collection Efficiency in a Wind Tunnel Icing Environment. In Proceedings of the 2018 Atmospheric and Space Environments Conference, Atlanta, GA, USA, 25–29 June 2018. [Google Scholar]
  29. Gray, V.H. Correlations among Ice Measurements, Impingement Rates, Icing Conditions and Drag of a 65A004 Airfoil; NACA TN 4151; NASA: Washington, DC, USA, 1958.
  30. Report of the 12A Working Group on Determination of Critical Ice Shapes for the Certification of Aircraft; DOT/FAA/AR-00/37; Federal Aviation Administration: Washington, DC, USA, 2000.
Figure 1. Model description and mesh configuration. (a) Sketch of the airfoil in a wind tunnel. (b) Boundary layer grids. (c) Grids near the edge of the airfoil.
Figure 1. Model description and mesh configuration. (a) Sketch of the airfoil in a wind tunnel. (b) Boundary layer grids. (c) Grids near the edge of the airfoil.
Aerospace 09 00587 g001
Figure 2. Calculated ice shapes at different grid numbers and icing steps.
Figure 2. Calculated ice shapes at different grid numbers and icing steps.
Aerospace 09 00587 g002
Figure 3. Comparison of ice shapes from the experiment and numerical simulations.
Figure 3. Comparison of ice shapes from the experiment and numerical simulations.
Aerospace 09 00587 g003
Figure 4. Characteristic parameters of the ice horn. (a) Ice horn dimensions. (b) Projection of upper ice horn in the direction of the incoming flow.
Figure 4. Characteristic parameters of the ice horn. (a) Ice horn dimensions. (b) Projection of upper ice horn in the direction of the incoming flow.
Aerospace 09 00587 g004
Figure 5. (a) Pressure coefficients distribution at different blockages. (b) Positions of stagnation point at different blockages.
Figure 5. (a) Pressure coefficients distribution at different blockages. (b) Positions of stagnation point at different blockages.
Aerospace 09 00587 g005
Figure 6. Ice shapes at different blockages. (a) Ice shape with 10% blockage. (b) Ice shape of NACA0012 airfoil at different blockages.
Figure 6. Ice shapes at different blockages. (a) Ice shape with 10% blockage. (b) Ice shape of NACA0012 airfoil at different blockages.
Aerospace 09 00587 g006
Figure 7. Droplet collection coefficients at different blockages of NACA0012 airfoil.
Figure 7. Droplet collection coefficients at different blockages of NACA0012 airfoil.
Aerospace 09 00587 g007
Figure 8. Ice shape of GLC305 airfoil at different blockages.
Figure 8. Ice shape of GLC305 airfoil at different blockages.
Aerospace 09 00587 g008
Figure 9. Relationship between the blockage and the dimensionless projection height.
Figure 9. Relationship between the blockage and the dimensionless projection height.
Aerospace 09 00587 g009
Table 1. Parameters of airfoils used in the numerical simulations.
Table 1. Parameters of airfoils used in the numerical simulations.
NACA0012 AirfoilGLC305 Airfoil
Chord length21 in36 in
Velocity102.8 m/s90 m/s
Angle of attack4.5°
Time7 min16.7 min
LWC0.55 g/m30.405 g/m3
MVD20 μm20 μm
Temperature265.37 K267.40 K
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lu, D.; Lu, Z.; Han, Z.; Xu, X.; Huang, Y. Numerical Investigation on the Effect of Blockage on the Icing of Airfoils. Aerospace 2022, 9, 587. https://doi.org/10.3390/aerospace9100587

AMA Style

Lu D, Lu Z, Han Z, Xu X, Huang Y. Numerical Investigation on the Effect of Blockage on the Icing of Airfoils. Aerospace. 2022; 9(10):587. https://doi.org/10.3390/aerospace9100587

Chicago/Turabian Style

Lu, Daixiao, Zhiliang Lu, Zhirong Han, Xian Xu, and Ying Huang. 2022. "Numerical Investigation on the Effect of Blockage on the Icing of Airfoils" Aerospace 9, no. 10: 587. https://doi.org/10.3390/aerospace9100587

APA Style

Lu, D., Lu, Z., Han, Z., Xu, X., & Huang, Y. (2022). Numerical Investigation on the Effect of Blockage on the Icing of Airfoils. Aerospace, 9(10), 587. https://doi.org/10.3390/aerospace9100587

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop