LEWICE/E: An Euler based ice accretion code

A new version of the LEWICE ice accretion computer code was developed which calculates the ice growth on two dimensional surfaces, incorporating the effects of compressibility through the solution of the Euler 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, and plot results. This code increases the applicability of ice accretion predictions by allowing calculations at higher Mach numbers. The new elements of the code are described. Calculated results are compared to experiment for several cases, including a LEWICE example case and a thin airfoil section at a Mach number of 0.58

Navier-Stokes analysis of airfoils with leading edge ice accretions

A numerical analysis of the flowfield characteristics and the performance degradation of an airfoil with leading edge ice accretions was performed. The important fluid dynamic processes were identified and calculated. Among these were the leading edge separation bubble at low angles of attack, complete separation on the low pressure surface resulting in premature shell, drag rise due to the ice shape, and the effects of angle of attack on the separated flow field. Comparisons to experimental results were conducted to confirm these calculations. A computer code which solves the Navier-Stokes equations in two dimensions, ARC2D, was used to perform the calculations. A Modified Mixing Length turbulence model was developed to produce grids for several ice shape and airfoil combinations. Results indicate that the ability to predict overall performance characteristics, such as lift and drag, at low angles of attack is excellent. Transition location is important for accurately determining separation bubble shape. Details of the flowfield in and downstream of the separated regions requires some modifications. Calculations for the stalled airfoil indicate periodic shedding of vorticity that was generated aft of the ice accretion. Time averaged pressure values produce results which compare favorably with experimental information. A turbulence model which accounts for the history effects in the flow may be justified

Swept wing ice accretion modeling

An effort to develop a three-dimensional modeling method was initiated. This first step towards creation of a complete aircraft icing simulation code builds on previously developed methods for calculating three-dimensional flow fields and particle trajectories combined with a two-dimensional ice accretion calculation along coordinate locations corresponding to streamlines. This work is a demonstration of the types of calculations necessary to predict a three-dimensional ice accretion. Results of calculations using the 3-D method for a MS-317 swept wing geometry are projected onto a 2-D plane normal to the wing leading edge and compared to 2-D 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 the flow field over the surface and the particle trajectories differed for the two calculations. This led to lower collection efficiencies, convective heat transfer coefficients, freezing fractions, and ultimately ice accumulation for the 3-D calculation

Software Released by LEWICE 2.0 Ice Accretion Software Development Project

Computational icing simulation methods are making the transition from the realm of research to commonplace use in design and certification. As such, standards of software management, design, validation, and documentation must be adjusted to accommodate the increased expectations of the user community with respect to accuracy, reliability, capability, and usability. With this in mind, in collaboration with Glenn's Engineering Design and Analysis Division, the Icing Branch of the NASA Glenn Research Center at Lewis Field began a software improvement project focused on the two-dimensional ice accretion simulation tool LEWICE. This project is serving as an introduction to the concepts of software management and is intended to serve as a pilot project for future icing simulation code development. The LEWICE 2.0 Software Development Project consisted of two major elements: software management and software validation. The software management element consisted of identifying features of well-designed and well-managed software that are appropriate for an analytical prediction tool such as LEWICE and applying them to a revised version of the code. This element included tasks such as identification of software requirements, development and implementation of coding standards, and implementation of software revision control practices. With the application of these techniques, the LEWICE ice accretion code became a more stable and reliable software product. In addition, the lessons learned about software development and maintenance can be factored into future software projects at the outset. The software validation activity was an integral part of our effort to make LEWICE a more accurate and reliable analysis tool. Because of the efforts taken to extensively validate this software, LEWICE 2.0 is more robust than previous releases and can reproduce results accurately across several computing platforms. It also differs from previous versions in the extensive quantitative comparisons of its results with a data base of ice shapes that were generated in Glenn's Icing Research Tunnel. The results of the shape comparisons were analyzed to determine the range of meteorological conditions under which LEWICE 2.0 is within the repeatability found for ice shapes obtained from experiments. These comparisons show that the average variation of LEWICE 2.0 from experimental data is 7.2 percent, whereas the overall variability of the experimental data is 2.5 percent. The LEWICE 2.0 Software Development Project resulted in the release, at the recent Icing Branch LEWICE Workshop, of a more useful and robust analysis tool than was previously available to the icing community. The lessons learned from this project will be applied to the Icing Branch s future software development efforts

Unstructured grid research and use at NASA Lewis Research Center

Computational fluid dynamics applications of grid research at LRC include inlets, nozzles, and ducts; turbomachinery; propellers - ducted and unducted; and aircraft icing. Some issues related to internal flow grid generation are resolution requirements on several boundaries, shock resolution vs. grid periodicity, grid spacing at blade/shroud gap, grid generation in turbine blade passages, and grid generation for inlet/nozzle geometries. Aircraft icing grid generation issues include (1) small structures relative to airfoil chord must be resolved; (2) excessive number of grid points in far-field using structured grid; and (3) grid must be recreated as ice shape grows

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

Users manual for the NASA Lewis three-dimensional ice accretion code (LEWICE 3D)

A description of the methodology, the algorithms, and the input and output data along with an example case for the NASA Lewis 3D ice accretion code (LEWICE3D) has been produced. The manual has been designed to help the user understand the capabilities, the methodologies, and the use of the code. The LEWICE3D code is a conglomeration of several codes for the purpose of calculating ice shapes on three-dimensional external surfaces. A three-dimensional external flow panel code is incorporated which has the capability of calculating flow about arbitrary 3D lifting and nonlifting bodies with external flow. A fourth order Runge-Kutta integration scheme is used to calculate arbitrary streamlines. An Adams type predictor-corrector trajectory integration scheme has been included to calculate arbitrary trajectories. Schemes for calculating tangent trajectories, collection efficiencies, and concentration factors for arbitrary regions of interest for single droplets or droplet distributions have been incorporated. A LEWICE 2D based heat transfer algorithm can be used to calculate ice accretions along surface streamlines. A geometry modification scheme is incorporated which calculates the new geometry based on the ice accretions generated at each section of interest. The three-dimensional ice accretion calculation is based on the LEWICE 2D calculation. Both codes calculate the flow, pressure distribution, and collection efficiency distribution along surface streamlines. For both codes the heat transfer calculation is divided into two regions, one above the stagnation point and one below the stagnation point, and solved for each region assuming a flat plate with pressure distribution. Water is assumed to follow the surface streamlines, hence starting at the stagnation zone any water that is not frozen out at a control volume is assumed to run back into the next control volume. After the amount of frozen water at each control volume has been calculated the geometry is modified by adding the ice at each control volume in the surface normal direction

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

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

The Influence of SLD Drop Size Distributions on Ice Accretion in the NASA Icing Research Tunnel

An ice shape database has been created to document ice accretions on a 21-inch chord NACA0012 model and a 72-inch chord NACA 23012 airfoil model resulting from an exposure to a Supercooled Large Drop (SLD) icing cloud with a bimodal drop size distribution. The ice shapes created were documented with photographs, laser scanned surface measurements over a section of the model span, and measurement of the ice mass over the same section of each accretion. The icing conditions used in the test matrix were based upon previously used conditions on the same models but with an alternate approach to evaluation of drop distribution effects. Ice shapes resulting from the bimodal distribution as well as from equivalent monomodal drop size distributions were obtained and compared. Results indicate that the ice shapes resulting from the monomodal and bimodal drop size distributions had similar shapes, but the bimodal distributions had greater mass and volume measurements and icing limits that extended further back on the chord of the model