Icing simulation: A survey of computer models and experimental facilities

A survey of the current methods for simulation of the response of an aircraft or aircraft subsystem to an icing encounter is presented. The topics discussed include a computer code modeling of aircraft icing and performance degradation, an evaluation of experimental facility simulation capabilities, and ice protection system evaluation tests in simulated icing conditions. Current research focussed on upgrading simulation fidelity of both experimental and computational methods is discussed. The need for increased understanding of the physical processes governing ice accretion, ice shedding, and iced airfoil aerodynamics is examined

Numerical simulation of ice growth on a MS-317 swept wing geometry

An effort to develop a 3-D ice accretion modeling method was initiated. This first step towards creation of a complete aircraft icing simulation code builds on previously developed methods for calculating 3-D flow fields and particle trajectories combined with a 2-D ice accretion calculation along coordinate locations corresponding to streamlines. The types of calculations necessary to predict 3-D ice accretion is demonstrated. Results of calculations using 3-D method for a MS-317 swept wing geometry are projected onto a 2-D plane parallel to the free stream direction and compared to experimental results for the same geometry. It is anticipated that many modifications will be made to this approach, however this effort will lay the groundwork for future modeling efforts. Results indicate that rime ice shapes indicate a difficulty in accurately calculating the ice shape in the runback region

The low frequency oscillation in the flow over a NACA0012 airfoil with an iced leading edge

The unusually low frequency oscillation in the wake of an airfoil is explored experimentally as well as computationally for a NACA0012 airfoil with a glaze ice accretion at the leading edge. Experimentally, flow oscillations were observed at low frequencies that correspond to a Strouhal number of about 0.02. This occurred in the angle of attack range of 8 to 9 deg, near the onset of static stall for this airfoil. With a Navier-Stokes computation, limit-cycle oscillations in the flow and in the aerodynamic forces were also observed at low Strouhal numbers. However, the occurrence of the oscillation is found to depend on the turbulence model in use as well as the Reynolds number

Numerical Modeling of Anti-icing Systems and Comparison to Test Results on a NACA 0012 Airfoil

A series of experimental tests were conducted in the NASA Lewis IRT on an electro-thermally heated NACA 0012 airfoil. Quantitative comparisons between the experimental results and those predicted by a computer simulation code were made to assess the validity of a recently developed anti-icing model. An infrared camera was utilized to scan the instantaneous temperature contours of the skin surface. Despite some experimental difficulties, good agreement between the numerical predictions and the experiment results were generally obtained for the surface temperature and the possibility for each runback to freeze. Some recommendations were given for an efficient operation of a thermal anti-icing system

Modeling of surface roughness effects on glaze ice accretion

A series of experimental investigations focused on studying the cause and effect of roughness on accreting glaze ice surfaces were conducted. Detailed microvideo observations were made of glaze ice accretions on 1 to 4 inch diameter cylinders in three icing wind tunnels (the Data Products of New England six inch test facility, the NASA Lewis Icing Research Tunnel, and the B. F. Goodrich Ice Protection Research Facility). Infrared thermal video recordings were made of accreting ice surfaces in the Goodrich facility. Distinct zones of surface water behavior were observed; a smooth wet zone in the stagnation region with a uniform water film; a rough zone where surface tension effects caused coalescence of surface water into stationary beads; a horn zone where roughness elements grow into horn shapes; a runback zone where surface water ran back as rivulets; and a dry zone where rime feathers formed. The location of the transition from the smooth to the rough zone was found to migrate with time towards the stagnation point. The behavior of the transition appeared to be controlled by boundary layer transition and bead formation mechanisms at the interface between the smooth and rough zones. Regions of wet ice growth and enhanced heat transfer were clearly visible in the infrared video recordings of glaze ice surfaces. A simple multi-zone modification to the current glaze ice accretion model was proposed to include spatial variability in surface roughness

Ice Accretion and Performance Degradation Calculations with LEWICE/NS

The LEWICE ice accretion computer code has been extended to include the solution of the two-dimensional Navier-Stokes equations. The code is modular and contains separate stand-alone program elements that create a grid, calculate the flow field parameters, calculate the droplet trajectory paths, determine the amount of ice growth, calculate aeroperformance changes, and plot results. The new elements of the code are described. Calculated results are compared to experiment for several cases, including both ice shape and drag rise

Simulation of iced wing aerodynamics

The sectional and total aerodynamic load characteristics of moderate aspect ratio wings with and without simulated glaze leading edge ice were studied both computationally, using a three dimensional, compressible Navier-Stokes solver, and experimentally. The wing has an untwisted, untapered planform shape with NACA 0012 airfoil section. The wing has an unswept and swept configuration with aspect ratios of 4.06 and 5.0. Comparisons of computed surface pressures and sectional loads with experimental data for identical configurations are given. The abrupt decrease in stall angle of attack for the wing, as a result of the leading edge ice formation, was demonstrated numerically and experimentally

An experimental investigation of multi-element airfoil ice accretion and resulting performance degradation

An investigation of the ice accretion pattern and performance characteristics of a multi-element airfoil was undertaken in the NASA Lewis 6- by 9-Foot Icing Research Tunnel. Several configurations of main airfoil, slat, and flaps were employed to examine the effects of ice accretion and provide further experimental information for code validation purposes. The text matrix consisted of glaze, rime, and mixed icing conditions. Airflow and icing cloud conditions were set to correspond to those typical of the operating environment anticipated tor a commercial transport vehicle. Results obtained included ice profile tracings, photographs of the ice accretions, and force balance measurements obtained both during the accretion process and in a post-accretion evaluation over a range of angles of attack. The tracings and photographs indicated significant accretions on the slat leading edge, in gaps between slat or flaps and the main wing, on the flap leading-edge surfaces, and on flap lower surfaces. Force measurments indicate the possibility of severe performance degradation, especially near C sub Lmax, for both light and heavy ice accretion and performance analysis codes presently in use. The LEWICE code was used to evaluate the ice accretion shape developed during one of the rime ice tests. The actual ice shape was then evaluated, using a Navier-Strokes code, for changes in performance characteristics. These predicted results were compared to the measured results and indicate very good agreement

Experimental Investigation of Ice Accretion Effects on a Swept Wing

An experimental investigation was conducted to study the effects of 2-, 5-, 10-, and 22.5-min ice accretions on the aerodynamic performance of a swept finite wing. The ice shapes tested included castings of ice accretions obtained from icing tests at the NASA Glenn Icing Research Tunnel (IRT) and simulated ice shapes obtained with the LEWICE 2.0 ice accretion code. The conditions used for the icing tests were selected to provide five glaze ice shapes with complete and incomplete scallop features and a small rime ice shape. The LEWICE ice shapes were defined for the same conditions as those used in the icing tests. All aerodynamic performance tests were conducted in the 7- x 10-ft Low-Speed Wind Tunnel Facility at Wichita State University. Six component force and moment measurements, aileron hinge moments, and surface pressures were obtained for a Reynolds number of 1.8 million based on mean aerodynamic chord and aileron deflections in the range of -15o to 20o. Tests were performed with the clean wing, six IRT ice shape castings, seven smooth LEWICE ice shapes, and seven rough LEWICE ice shapes. Roughness for the LEWICE ice shapes was simulated with 36-size grit. The experiments conducted showed that the glaze ice castings reduced the maximum lift coefficient of the clean wing by 11.5% to 93.6%, while the 5-min rime ice casting increased maximum lift by 3.4%. Minimum iced wing drag was 133% to 3533% greater with respect to the clean case. The drag of the iced wing near the clean wing stall angle of attack was 17% to 104% higher than that of the clean case. In general, the aileron remained effective in changing the lift of the clean and iced wings for all angles of attack and aileron deflections tested. Aileron hinge moments for the iced wing cases remained within the maximum and minimum limits defined by the clean wing hinge moments. Tests conducted with the LEWICE ice shapes showed that in general the trends in aerodynamic performance degradation of the wing with the simulated ice shapes were similar to those obtained with the IRT ice shape castings. However, in most cases, the ice castings resulted in greater aerodynamic performance losses than those obtained with the LEWICE ice shapes. For the majority of the LEWICE ice shapes, the addition of 36-size grit roughness to the smooth ice shapes increased aerodynamic performance losses

Progress toward the development of an airfoil icing analysis capability

The NASA-Lewis aircraft icing analysis program is composed of three major sub-programs. These sub-programs are ice accretion simulation, performance degradation evaluation, and ice protection system evaluation. These topics cover all areas of concern related to the simulation of aircraft icing and its consequences. The motivation for these activities is twofold, reduction of time and effort required in experimental programs and the ability to provide reliable information for aircraft certification in icing, over the complete range of environmental conditions. In addition to the analytical activities associated with development of these codes, several experimental programs are underway to provide verification information for existing codes. These experimental programs are also used to investigate the physical processes associated with ice accretion and removal for improvement of present analytical models. The NASA-Lewis icing analysis program is thus striving to provide a full range of analytical tools necessary for evaluation of the consequences of icing and of ice protection systems