Richard E. Kreeger

Implementation and Validation of 3-D Ice Accretion Measurement Methodology

A research program has been implemented to develop and validate the use of a commercial 3-D laser scanning system to record ice accretion geometry in the NASA Icing Research Tunnel. A main component of the program was the geometric assessment of the 3- D laser scanning system on a 2-D (straight wing) and a 3-D (swept wing) airfoil geometries. This exercise consisted of comparison of scanned ice accretion to castings of the same ice accretion. The scan data were also used to create rapid prototype artificial ice shapes that were scanned and compared to the original ice accretion. The results from geometric comparisons on the straight wing showed that the ice shape models generated through the scan/rapid prototype process compared reasonably well with the cast shapes. Similar results were obtained with the geometric comparisons on the swept wing. It was difficult to precisely compare the scans of the cast shapes to the original ice accretion scans because the cast shapes appear to have shrunk during the mold/casting process by as much as 0.10-inch. However the comparison of the local ice-shape features were possible and produced better results. The rapid prototype manufacturing process was shown to reproduce the original ice accretion scan normally within 0.01-inch.

Assessment of Ice Shape Roughness Using a Self-Orgainizing Map Approach

Self-organizing maps are neural-network techniques for representing noisy, multidimensional data aligned along a lower-dimensional and nonlinear manifold. For a large set of noisy data, each element of a finite set of codebook vectors is iteratively moved in the direction of the data closest to the winner codebook vector. Through successive iterations, the codebook vectors begin to align with the trends of the higher-dimensional data. Prior investigations of ice shapes have focused on using self-organizing maps to characterize mean ice forms. The Icing Research Branch has recently acquired a high resolution three dimensional scanner system capable of resolving ice shape surface roughness. A method is presented for the evaluation of surface roughness variations using high-resolution surface scans based on a self-organizing map representation of the mean ice shape. The new method is demonstrated for 1) an 18-in. NACA 23012 airfoil 2 AOA just after the initial ice coverage of the leading 5 of the suction surface of the airfoil, 2) a 21-in. NACA 0012 at 0AOA following coverage of the leading 10 of the airfoil surface, and 3) a cold-soaked 21-in.NACA 0012 airfoil without ice. The SOM method resulted in descriptions of the statistical coverage limits and a quantitative representation of early stages of ice roughness formation on the airfoils. Limitations of the SOM method are explored, and the uncertainty limits of the method are investigated using the non-iced NACA 0012 airfoil measurements.

Ice Roughness in Short Duration SLD Icing Events

Ice accretion codes depend on models of roughness parameters to account for the enhanced heat transfer during the ice accretion process. While mitigating supercooled large droplet (SLD or Appendix O) icing is a significant concern for manufacturers seeking future vehicle certification due to the pending regulation, historical ice roughness studies have been performed using Appendix C icing clouds which exhibit mean volumetric diameters (MVD) much smaller than SLD clouds. Further, the historical studies of roughness focused on extracting parametric representations of ice roughness using multiple images of roughness elements. In this study, the ice roughness developed on a 21-in. NACA 0012 at 0deg angle of attack exposed to short duration SLD icing events was measured in the Icing Research Tunnel at the NASA Glenn Research Center. The MVDs used in the study ranged from 100 micrometer to 200 micrometers, in a 67 m/s flow, with liquid water contents of either 0.6 gm/cubic meters or 0.75 gm/cubic meters. The ice surfaces were measured using a Romer Absolute Arm laser scanning system. The roughness associated with each surface point cloud was measured using the two-dimensional self-organizing map approach developed by McClain and Kreeger (2013) resulting in statistical descriptions of the ice roughness.

A Finite-Volume Approach to Modeling Ice Accretion

In this work we present a novel, generalized, multiscale physics, unstructured finite-volume, CFD approach for simulating ice accretion on aircraft. A multi-physics solver that evaluates the (1) air flow, (2) droplet trajectories, (3) surface-liquid flow, (4) solidification, and (5) computes the deformed ice shape, is presented. Initial results show promise in the developed methods and solvers, that are expected to later be extended for future rotorcraft ice-accretion analysis. Initial validation cases are presented for the various components of the solver, and compare reasonable well with LEWICE and experiments for simple geometries. This initial capability displays a capability that could be extended, in future efforts, with more detailed models and provide ice shapes of similar quality as the current methodologies, while providing a capability that extends to more complex configurations such as rotorcraft.

Ice Particle Impacts on a Moving Wedge

This work presents the results of an experimental study of ice particle impacts on a moving wedge. The experiment was conducted in the Adverse Environment Rotor Test Stand (AERTS) facility located at Penn State University. The wedge was placed at the tip of a rotating blade. Ice particles shot from a pressure gun intercepted the moving wedge and impacted it at a location along its circular path. The upward velocity of the ice particles varied from 7 to 12 meters per second. Wedge velocities were varied from 0 to 120 meters per second. Wedge angles tested were 0 deg, 30 deg, 45 deg, and 60 deg. High speed imaging combined with backlighting captured the impact allowing observation of the effect of velocity and wedge angle on the impact and the post-impact fragment behavior. It was found that the pressure gun and the rotating wedge could be synchronized to consistently obtain ice particle impacts on the target wedge. It was observed that the number of fragments increase with the normal component of the impact velocity. Particle fragments ejected immediately after impact showed velocities higher than the impact velocity. The results followed the major qualitative features observed by other researchers for hailstone impacts, even though the reduced scale size of the particles used in the present experiment as compared to hailstones was 4:1.

Experimental Measurement of Frozen and Partially Melted Water Droplet Impact Dynamics

High-speed video of single frozen water droplets impacting a surface was acquired. The droplets diameter ranged from 0.4 mm to 0.9 mm and impacted at velocities ranging from 140 m/sec to 309 m/sec. The techniques used to freeze the droplets and launch the particles against the surfaces is described in this paper. High-speed video was used to quantify the ice accretion area to the surface for varying impact angles (30 deg, 45 deg, 60 deg), impacting velocities, and break-up angles. An oxygen /acetylene cross-flow flame used to ensure partial melting of the traveling frozen droplets is also discussed. A linear relationship between impact angle and ice accretion is identified for fully frozen particles. The slope of the relationship is affected by impact speed. Perpendicular impacts, i.e. 30 deg, exhibited small differences in ice accretion for varying velocities, while an increase of 60% in velocity from 161 m/sec to 259 m/sec, provided an increase on ice accretion area of 96% at an impact angle of 60 deg. The increase accretion area highlights the importance of impact angle and velocity on the ice accretion process of ice crystals. It was experimentally observed that partial melting was not required for ice accretion at the tested velocities when high impact angles were used (45 and 60 deg). Partially melted droplets doubled the ice accretion areas on the impacting surface when 0.0023 Joules were applied to the particle. The partially melted state of the droplets and a method to quantify the percentage increase in ice accretion area is also described in the paper.

Ice Accretion Modeling using an Eulerian Approach for Droplet Impingement

A three-dimensional Eulerian analysis has been developed for modeling droplet impingement on lifting bodes. The Eulerian model solves the conservation equations of mass and momentum to obtain the droplet flow field properties on the same mesh used in CFD simulations. For complex configurations such as a f ull rotorcraft, the Eulerian approach is more efficient because the Lagrangian approach would require a significant amount of seeding for accurate estimates of collection effici ency. Simulations are done for various benchmark cases such as NACA0012 airfoil, MS317 airfoil and oscillating SC2110 airfoil to illustrate its use. The present results are compare d with results from the Lagrangian approach used in an industry standard analysis call ed LEWICE.

Heat transfer from protuberances

An experiment by Henry et al. explored the heat transfer of flows over protuberances in laminar and turbulent flow, simulating conditions during the beginning stages of glaze icing. This paper represents an effort to explain the heat transfer enhancement of roughness elements and protuberances that was observed by Henry et aL In the experiments of Henry et al., a single roughness element was placed on a heated flat plate. The temperature along the flat plate and along the roughness element was measured using an infrared camera to determine the enhancement of heat transfer of the protuberance as opposed to the smooth surface. A one-dimensional extended-surface (fin) analysis was performed to examine the results of Henry et al. Although significant assumptions were made using the extended-surface analysis, the important trends of the Henry et al. data were captured. The extended-surface analysis captured the trends of the apparent enhancement as reported by Henry et al. vs the Reynolds number based on the location from the leading edge of the surface and vs the ratio of the protuberance height to the boundary-layer thickness. Although the absolute magnitudes of the apparent enhancement are overestimated by the extended-surface analysis, the matched trends indicate the importance of the thermal conductivity of the protuberance, the importance of the interaction of the protuberance with the thermal boundary layer, and the importance of radiation into the protuberance.

Heat Transfer over Roughness Elements Larger than the Boundary Layer

*† ‡ The development of an ice feather growth model requires some insight into the effects of roughness on heat transfer. The objective of this study was to find the heat enhancement caused by the acceleration of the flow over an isolated roughness element in laminar flow. The temperature and heat transfer were calculated in 2D along a flat plate and along the roughness element for values of k/δ 0.5 to 3.0 and Re from 200 to 2000. The problem was then extended to 3D turbulent flow.

Ice Particle Impacts on a Flat Plate

An experimental study was conducted at the Ballistic Laboratory of NASA Glenn Research Center to study the impact of ice particles on a stationary flat surface target set at 45 degrees with respect to the direction of motion of the impinging particle (Figure 1). The experiment is part of NASA efforts to study the physics involved in engine power-loss events due to ice-crystal ingestion and ice accretion formation inside engines. These events can occur when aircraft encounter high-altitude convective weather.