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.