Andy P Broeren | Engineering Papers

NASA NTRS2025

Details of Adhering and Testing Artificial Ice Shapes to the High Lift Common Research Model in the National Transonic Facility

There is currently interest in conducting flight Reynolds number aerodynamic performance testing for iced airplane configurations. There has only been limited testing in the past, none of which has resulted in publicly available data. Consequently, there is virtually no documented knowledge for conducting icing related testing at flight Reynolds number where testing in cryogenic conditions is required. The National Transonic Facility (NTF) is a cryogenic wind tunnel at NASA Langley Research Center that is capable of achieving flight Reynolds numbers on subscale models. A large number of aerodynamic tests have been performed in this facility for uniced subscale model configurations. Through the collaborations associated with the High Lift Common Research Model, a test campaign was conducted to fill the existing research gaps with publicly available data. This resulted in a requirement that NASA develop a suitable level of competency for performing icing related tests in cryogenic environments. This document provides information about the challenges faced during the test campaign, how the issues were mitigated, and preserves the institutional knowledge gained by providing recommendations that will be used to reduce risk for future testing.

Iced Aerodynamics↗

NASA NTRSDate not supplied

Ice Shape Classification for Swept Wings

An important consideration for swept-wing ice shapes is the level of geometric fidelity required to accurately capture the aerodynamics. One way to help make the task of addressing this more manageable is to organize the ice shapes into a small number of classifications. These classifications are done by grouping similar ice shape geometries that also have similar iced wing aerodynamics. Ice shape classifications for airfoils are more mature and are reviewed as an example. Swept wing ice shape classifications are based on the well-established airfoil classifications and here five classifications are proposed: leading-edge roughness, streamwise ice, 3D leading-edge horn, highly 3D leading-edge horn and spanwise ridge. Recent data on swept wing icing aerodynamics is used to develop and support this classification. The development, or suppression, of the leading-edge separation vortex on swept wings without a leading-edge device is seen to be key in the classification of the highly 3D leading-edge horn.

Aircraft icing↗

NASA NTRSDate not supplied

Effect of Geometric Fidelity on the Aerodynamics of a Swept Wing with Glaze Ice Accretion

Aerodynamic assessment of icing effects on swept wings is an important component of a larger effort to improve three-dimensional icing simulation capabilities. An understanding of ice-shape geometric fidelity on iced-wing aerodynamics and the associated flowfield features are needed to guide the development and validation of ice-accretion simulation tools. To this end, wind-tunnel testing was carried out for 8.9% and 13.3% scale semispan wing models based upon the Common Research Model airplane configuration. Various levels of geometric fidelity of an artificial ice shape representing a glaze-ice accretion on a swept wing were investigated. The highest fidelity artificial ice shape reproduced all of the three-dimensional features associated with the glaze ice accretion. The lowest fidelity artificial ice shapes were simple, spanwise-varying horn ice geometries intended to represent the maximum ice thickness on the wing upper surface. The results presented in this paper show that the addition of grit roughness to some lower-fidelity artificial ice shapes resulted in favorable lift and pitching moment comparisons to the wing with the highest fidelity artificial ice shape. In the range of 4.3 to 7.4 deg. angle of attack, surface oil flow visualization and pressure data show that the wing with the two lower fidelity simulations clearly demonstrated a leading edge vortex dominated flowfield, referred to as type I. For the wing with the high fidelity ice shape, the flowfield at lower angles of attack was characterized by streamwise-running, counter-rotating vortical flow referred to as type II. Between 6.4 and 7.4 deg. angle of attack, the effect of the type II flow structures was significantly altered and gave way to the type I leading edge vortex. This means that for angles of attack 7.4 deg. and higher, the wing with all three configurations exhibited the same type of flowfield. This helps to explain why there is reasonably good agreement in the lift and pitching moment coefficients among these configurations.

Aircraft icing↗

NASA NTRSDate not supplied

Influence of Airfoil Angle of Attack on Ice Accretion Roughness

The influence of airfoil angle attack on roughness evolution and spatial variation was investigated in the Icing Research Tunnel (IRT) at the NASA Glenn Research Center. Two airfoil models were used for the study: a 53.34-cm (21-in.) NACA 0012 model and a 152.4-cm (60-in.) HAARP-II model. For the NACA 0012, the angle of attack was varied from 0° to 3°, and the ice accretion roughness was characterized while keeping the other icing parameters such as freestream velocity, freestream stagnation temperature, accumulation parameter, and median volumetric diameter (MVD) constant. For the HAARP-II, the angle of attack was varied from -1.9°, which is the non-lifting angle of attack for the airfoil, to 3°. A series of accumulation time progression cases for the HAARP-II at 2.3° angle of attack was also performed in Appendix C and Appendix O conditions and compared to the measurements presented by McClain et al. (2018) for the non-lifting angle of attack condition. The results demonstrate the strong influence of the changing local static pressure along the surface influencing the roughness characteristics. A correlation approach for the maximum equivalent sand-grain roughness along the surface of an airfoil is presented and considers the scaled accumulation time, the stagnation point freezing fraction, and the surface pressure relative to the static and stagnation pressures of the flow. Finally, an attempt to relate the maximum roughness quantities to local pressure coefficients and local icing quantities is presented and discussed.

Icing↗

NASA NTRSDate not supplied

Effect of Water Droplets Crossing the Boundary Layer in a Stagnation Point Configuration

An experimental study was conducted in the Vertical Icing Studies Tunnel at the Icing Physics Flow Laboratory of NASA Glenn Research Center to study the effect of water droplets crossing the boundary layer in a stagnation point configuration. The objective of the experiment was to determine if water droplets crossing a boundary layer create turbulent spots that accelerate the boundary layer transition from laminar to turbulent. Water droplets that crossed the boundary layer were generated with a nozzle installed in the plenum. The turbulence level in the boundary layer was measured with a hot wire system when the nozzle was off, air on, and air and water on. The results indicate that the presence of the nozzle alone did not affect the boundary layer. The activation of the nozzle to eject air only or to generate water droplets affected the turbulence level in the boundary layer. The continued study of this effect is needed because of its implications in the development of heat transfer models for icing codes.

Icing↗

NASA NTRSDate not supplied

Comparison of Ice Shapes on Full-Chord and Truncated Swept Wing Models

A research program was conducted to evaluate the effectiveness of icing tunnel hybrid model design. A hybrid design is where the full-scale leading edge of a wing section is maintained only to a certain percentage of the local chord, while the aft section of the model is redesigned into a shortened or truncated planform. Ice shape geometries on a full-chord length and hybrid versions of the swept CRM65 wing model were obtained in NASA Icing Research Tunnel and compared. For most test conditions, the overall size and shape of the ice shapes compared well. On the main part of the ice shape, the hybrid model had larger chordwise extent, resulting in wider ice shape. For the scallop ice shapes, the ice shape heights at attachment line location were identical. For conditions that produced rime shape, hybrid model ice shapes were larger overall, even at the attachment line location. The only location where there was more ice on the full chord model was on the lower surface, downstream of where the hybrid model no longer had full-scale leading edge. For the test conditions where ice mass measurements were acquired for both models, the full-chord models had 11-22% less ice mass than on the hybrid wing. These differences were larger and opposite of those predicted from LEWICE 3D ice prediction calculations. Further investigation is needed in order determine the cause of the differences between the full-chord and hybrid models.

Icing↗

NASA NTRSDate not supplied

Ice Accretion Roughness Variations on a Hybrid CRM65-Midspan Wing Model

Ice accretion roughness measurements were performed in the Icing Research Tunnel (IRT) at NASA Glenn Research Center for the Hybrid CRM65-Midspan model in a range of icing conditions. The Hybrid CRM65-Midspan model was chosen for this investigation because 1) the model exhibits high sweep relative to models previously explored in the roughness investigations, 2) the model has leading edge characteristics similar to wing shapes currently used in mid-size commercial airliners, and 3) the sweep and thickness ratios relate better to hybrid lifting body designs for N+2 and N+3 vehicles than other models available. The investigation consisted of multiple sets of tests which focused on 1) 0-angle of attack cases replicating the conditions employed by Anderson et al. (1998) using both Appendix C and SLD cloud conditions, 2) cases based on the “Max Scallop” case by Broeren et al. (2016) and a “High Temperature” case with cloud properties similar to the “Max Scallop” case. Additional tests were performed 1) based on the “Max Scallop” case with variations in freestream static temperature and 2) using test section speeds near 10,000-hold flight conditions. The point clouds were characterized using the approach of McClain and Kreeger (2013) for the ice roughness variations and using the approach of McClain (2016) for the mean ice thickness variations. The resulting roughness and mean thickness variations generally follow the temporal scaling previously identified using on airfoil models without sweep, but the collapse of the time progression profiles is not as tight as found for past measurements on models without sweep. LEWICE and modified panel-method predictions were used explore spatial roughness variations and to compare to the roughness correlations developed by McClain et al. (2021) for the “Max Scallop” cases.

Icing↗

NASA NTRSDate not supplied

A Model for Ice Accretion Roughness Evolution and Spatial Variations

Over the past decade, multiple investigations of ice accretion roughness and spatial variations have been performed in the Icing Research Tunnel (IRT) at the NASA Glenn Research Center. The early investigations used models of NACA 0012 airfoils with different chord sizes and focused on temporal scaling and primary cloud scaling parameters such as stagnation point collection efficiency. Subsequent investigations included the effects of model sweep, airfoil shape, airfoil lifting condition, and freestream static temperature. To develop a predictive model for roughness evolution in generalized icing situations, the maximum roughness values for the airfoils and conditions used in the angle of attack and freestream temperature investigations were scaled to eliminate the temporal variations. LEWICE simulations were employed to identify the local collection efficiency at the locations of maximum roughness, and a two-dimensional panel-method code was used to identify the local static pressure at the location of the maximum roughness. Because of the stochastic nature of ice accretion roughness, a physics-directed approach was employed to develop a multi-dimensional correlation based on the local pressure coefficient, the local total temperature, the cloud properties, and the cloud exposure time. The resulting correlation predictions are compared to ice shapes measured in the IRT for both a 21-in. NACA 0012 airfoil model and a 60-in. HAARP-II model. The resulting predictions indicate that the local freezing fraction must be included in the roughness predictive model. Implications regarding icing heat transfer predictions using the resulting roughness model are also discussed.

Icing↗

NASA NTRSDate not supplied

Analysis of Ice Mass Growth Over Time on the CRM65 Midspan Hybrid Model

The Aeronautics Research Mission Directorate at NASA is developing and applying tools to enable future technologies towards sustainable flight. Aircraft icing has been identified as a potential barrier to entry into service for innovative designs necessitating improvements to computational ice accretion tools. NASA is developing the Glenn Icing Computational Environment (GlennICE) to address deficiencies in the computational modeling capabilities of previously developed ice accretion solvers. To benchmark and improve the ability to model highly three-dimensional ice accretion, high quality validation data against experimental data is required. The CRM65 Midspan Hybrid geometry was previously tested at the NASA Icing Research Tunnel to generate experimental data for swept wing geometries typical for commercial transport aircraft. As a part of a 2018 icing test campaign, experimental data characterizing the relationship between ice accretion time and ice mass growth was obtained and can be leveraged for use in validation of computational tools. The desire for computational ice accretion solvers to predict ice shapes profiles accreted experimentally has often overshadowed the comparison to the mass and bulk volume of ice accreted. To address this deficiency, an analysis is presented in which GlennICE is applied to simulations of the CRM65 Midspan Hybrid model tested in the NASA Icing Research Tunnel. Results from the computational fluid dynamics simulations compared favorably to the experimental pressure coefficient data, thus validating the modeling setup. The experimental data showed excellent repeatability for the 15.0 minute accretion time. The comparisons between the experimental and computational ice mass over time showed good agreement up to 10.0 minutes after which the ice mass was underpredicted. The experimental ice mass was largely linear with some nonlinear data. The bulk volume of ice accreted experimentally compared well to GlennICE for the scanned ice shapes and mean combined cross section ice shapes, but was underpredicted for the maximum combined cross section ice shapes at longer accretion times. The experimental minimum combined cross section, mean combined cross section, and maximum combined cross section profiles when compared to GlennICE show good agreement for the mean combined cross section up to 15.0 minutes. The analyses show that with a single-shot method, GlennICE currently underpredicts the ice mass for longer accretion times, is not able to match the bulk volume of the maximum combined cross section due to dominating scallop features, and future work is required to generate a more generalized ice bulk density model.

Icing↗

NASA NTRSDate not supplied

Analysis of Ice Mass Growth Over Time on the CRM65 Midspan Hybrid Model

The Aeronautics Research Mission Directorate at NASA is developing and applying tools to enable future technologies towards sustainable flight. Aircraft icing has been identified as a potential barrier to entry into service for innovative designs necessitating improvements to computational ice accretion tools. NASA is developing the Glenn Icing Computational Environment (GlennICE) to address deficiencies in the computational modeling capabilities of previously developed ice accretion solvers. To benchmark and improve the ability to model highly three-dimensional ice accretion, high quality validation data against experimental data is required. The CRM65 Midspan Hybrid geometry was previously tested at the NASA Icing Research Tunnel to generate experimental data for swept wing geometries typical for commercial transport aircraft. As a part of a 2018 icing test campaign, experimental data characterizing the relationship between ice accretion time and ice mass growth was obtained and can be leveraged for use in validation of computational tools. The desire for computational ice accretion solvers to predict ice shapes profiles accreted experimentally has often overshadowed the comparison to the mass and bulk volume of ice accreted. To address this deficiency, an analysis is presented in which GlennICE is applied to simulations of the CRM65 Midspan Hybrid model tested in the NASA Icing Research Tunnel. Results from the computational fluid dynamics simulations compared favorably to the experimental pressure coefficient data, thus validating the modeling setup. The experimental data showed excellent repeatability for the 15.0 minute accretion time. The comparisons between the experimental and computational ice mass over time showed good agreement up to 10.0 minutes after which the ice mass was underpredicted. The experimental ice mass was largely linear with some nonlinear data. The bulk volume of ice accreted experimentally compared well to GlennICE for the scanned ice shapes and mean combined cross section ice shapes, but was underpredicted for the maximum combined cross section ice shapes at longer accretion times. The experimental minimum combined cross section, mean combined cross section, and maximum combined cross section profiles when compared to GlennICE show good agreement for the mean combined cross section up to 15.0 minutes. The analyses show that with a single-shot method, GlennICE currently underpredicts the ice mass for longer accretion times, is not able to match the bulk volume of the maximum combined cross section due to dominating scallop features, and future work is required to generate a more generalized ice bulk density model.

Icing↗