Application of Extended Messinger Models to Complex Geometries

Since, ice accretion can significantly degrade the performance and the stability of an airborne vehicle, it is imperative to be able to model it accurately. While ice accretion studies have been performed on airplane wings and helicopter blades in abundance, there are few that attempt to model the process on more complex geometries such as fuselages. This paper proposes a methodology that extends an existing in-house Extended Messinger solver to complex geometries by introducing the capability to work with unstructured grids and carry out spatial surface streamwise marching. For the work presented here commercial solvers such as STAR-CCM+ and ANSYS Fluent are used for the flow field and droplet dispersed phase computations. The ice accretion is carried out using an in-house icing solver called GT-ICE. The predictions by GT-ICE are compared to available experimental data, or to predictions by other solvers such as LEWICE and STAR-CCM+. Three different cases with varying levels of complexity are presented. The first case considered is a commercial transport airfoil, followed by a three-dimensional MS(1)-317 swept wing. Finally, ice accretion calculations performed on a Robin fuselage have been discussed. Good agreement with experimental data, where applicable, is observed. Differences between the ice accretion predictions by different solvers have been discussed

Measurement of the Critical Distance Parameter Against Icing Conditions on a NACA 0012 Swept Wing Tip

This work presents the results of three experiments, one conducted in the Icing Research Tunnel (IRT) at NASA Glenn Research Center and two in the Goodrich Icing Wind Tunnel (IWT). The experiments were designed to measure the critical distance parameter on a NACA 0012 Swept Wing Tip at sweep angles of 45deg, 30deg, and 15deg. A time sequence imaging technique (TSIT) was used to obtain real time close-up imaging data during the first 2 min of the ice accretion formation. The time sequence photographic data was used to measure the critical distance at each icing condition and to study how it develops in real time. The effect on the critical distance of liquid water content, drop size, total temperature, and velocity was studied. The results were interpreted using a simple energy balance on a roughness elemen

Overview of Icing Research at NASA Glenn

The aviation industry continues to deal with icing-related incidents and accidents on a regular basis. Air traffic continues to increase, placing more aircraft in adverse icing conditions more frequently and for longer periods. Icing conditions once considered rare or of little consequence, such as super-cooled large droplet icing or high altitude ice crystals, have emerged as major concerns for modern aviation. Because of this, there is a need to better understand the atmospheric environment, the fundamental mechanisms and characteristics of ice growth, and the aerodynamic effects due to icing, as well as how best to protect these aircraft. The icing branch at NASA Glenn continues to develop icing simulation methods and engineering tools to address current aviation safety issues in airframe, engine and rotorcraft icing

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

Collaborative Experiments and Computations in Aircraft Icing

The formation of ice over lifting surfaces can affect aerodynamic performance. The ability to predict ice accumulation and the resulting degradation in vehicle performance is essential to determine the limitations of aircraft in icing encounters. The consequences of underestimating performance degradation can be serious and so it is important to produce accurate predictions, particularly for severe icing conditions. The simulation of ice accretion is a challenging multidisciplinary problem that requires close collaboration between the computational and ground test communities. This paper describes three recent case studies and the lessons learned through collaborative experiments and computations in aircraft icing- one for large commercial transports, one for rotorcraft, and one dealing with icing on regional jets

Evaluation of Scaling Methods for Rotorcraft Icing

This paper reports result of an experimental study in the NASA Glenn Icing Research Tunnel (IRT) to evaluate how well the current recommended scaling methods developed for fixed-wing unprotected surface icing applications might apply to representative rotor blades at finite angle of attack. Unlike the fixed-wing case, there is no single scaling method that has been systematically developed and evaluated for rotorcraft icing applications. In the present study, scaling was based on the modified Ruff method with scale velocity determined by maintaining constant Weber number. Models were unswept NACA 0012 wing sections. The reference model had a chord of 91.4 cm and scale model had a chord of 35.6 cm. Reference tests were conducted with velocities of 76 and 100 kt (39 and 52 m/s), droplet MVDs of 150 and 195 fun, and with stagnation-point freezing fractions of 0.3 and 0.5 at angle of attack of 0deg and 5deg. It was shown that good ice shape scaling was achieved for NACA 0012 airfoils with angle of attack lip to 5deg

Application of a High-Fidelity Icing Analysis Method to a Model-Scale Rotor in Forward Flight

An icing analysis process involving the loose coupling of OVERFLOW-RCAS for rotor performance prediction and with LEWICE3D for thermal analysis and ice accretion is applied to a model-scale rotor for validation. The process offers high-fidelity rotor analysis for the noniced and iced rotor performance evaluation that accounts for the interaction of nonlinear aerodynamics with blade elastic deformations. Ice accumulation prediction also involves loosely coupled data exchanges between OVERFLOW and LEWICE3D to produce accurate ice shapes. Validation of the process uses data collected in the 1993 icing test involving Sikorsky's Powered Force Model. Non-iced and iced rotor performance predictions are compared to experimental measurements as are predicted ice shapes

Effects of Hydrogen Bonding and Molecular Chain Flexibility of Substituted n-Alkyldimethylsilanes On Impact Ice Adhesion Shear Strength

The effects of hydrogen bonding and molecular flexibility upon ice adhesion shear strength were investigated using aluminum substrates coated with substituted n-alkyldimethylalkoxysilanes. The location of the chemical group substitution was on the opposing end of the linear n-alkyl chain with respect to silicon. Three hydrogen-bonding characteristics were evaluated: 1) non-hydrogen bonding, 2) donor/acceptor, and 3) acceptor. Varying the length of the n-alkyl chain provided an assessment of molecular chain flexibility. Coated and uncoated aluminum surfaces were characterized by receding water contact angle and surface roughness. Ice adhesion shear strength was determined in the Adverse Environment Rotor Test Stand facility from -16 to -8C that simulated aircraft in-flight icing conditions within the FAR Part 25/29 Appendix C icing envelope. Surface roughness of the coatings was similar allowing for comparison of the test results. An adhesion reduction factor, based on the ice adhesion shear strength data with respect to uncoated aluminum obtained at the same temperature, was calculated to compare the data. The results revealed complex interactions with impacting supercooled water droplets that were interdependent upon ice accretion temperature, surface energy characteristics of water and ice, hydrogen bonding characteristic of the substituent, and length of the n-alkyl chain. To aid in explaining the results, 1) changes in the surface energy component (i.e., non-polar and polar) values that water undergoes during its phase change from liquid to solid that arise from the freezing of impacting supercooled water droplets on the surface depended upon the temperature during accretion were taken into account and 2) the physical properties (i.e., water solubility and melting point) of small compounds analogous to the substituted n-alkyldimethylalkoxysilanes used in this study were compared

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

A Non-Cut Cell Immersed Boundary Method for Use in Icing Simulations

This paper describes a computational fluid dynamic method used for modelling changes in aircraft geometry due to icing. While an aircraft undergoes icing, the accumulated ice results in a geometric alteration of the aerodynamic surfaces. In computational simulations for icing, it is necessary that the corresponding geometric change is taken into consideration. The method used, herein, for the representation of the geometric change due to icing is a non-cut cell Immersed Boundary Method (IBM). Computational cells that are in a body fitted grid of a clean aerodynamic geometry that are inside a predicted ice formation are identified. An IBM is then used to change these cells from being active computational cells to having properties of viscous solid bodies. This method has been implemented in the NASA developed node centered, finite volume computational fluid dynamics code, FUN3D. The presented capability is tested for two-dimensional airfoils including a clean airfoil, an iced airfoil, and an airfoil in harmonic pitching motion about its quarter chord. For these simulations velocity contours, pressure distributions, coefficients of lift, coefficients of drag, and coefficients of pitching moment about the airfoil's quarter chord are computed and used for comparison against experimental results, a higher order panel method code with viscous effects, XFOIL, and the results from FUN3D's original solution process. The results of the IBM simulations show that the accuracy of the IBM compares satisfactorily with the experimental results, XFOIL results, and the results from FUN3D's original solution process