Search NASA⌕ Search

SEARCH · Search NASA

Results for “ice accretion”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

Flight Testing an Iced Business Jet for Flight Simulation Model Validation

A flight test of a business jet aircraft with various ice accretions was performed to obtain data to validate flight simulation models developed through wind tunnel tests. Three types of ice accretions were tested: pre-activation roughness, runback shapes that form downstream of the thermal wing ice protection system, and a wing ice protection system failure shape. The high fidelity flight simulation models of this business jet aircraft were validated using a software tool called "Overdrive." Through comparisons of flight-extracted aerodynamic forces and moments to simulation-predicted forces and moments, the simulation models were successfully validated. Only minor adjustments in the simulation database were required to obtain adequate match, signifying the process used to develop the simulation models was successful. The simulation models were implemented in the NASA Ice Contamination Effects Flight Training Device (ICEFTD) to enable company pilots to evaluate flight characteristics of the simulation models. By and large, the pilots confirmed good similarities in the flight characteristics when compared to the real airplane. However, pilots noted pitch up tendencies at stall with the flaps extended that were not representative of the airplane and identified some differences in pilot forces. The elevator hinge moment model and implementation of the control forces on the ICEFTD were identified as a driver in the pitch ups and control force issues, and will be an area for future work.

Ratvasky, Thomas P.↗

Passive infrared ice detection for helicopter applications

A technique is proposed to remotely detect rotor icing on helicopters by using passive IR thermometry to detect the warming caused by latent heat release as supercooled water freezes. During icing, the ice accretion region will be warmer than the uniced trailing edge, resulting in a characteristic chordwise temperature profile. Preliminary tests were conducted on a static model in the NASA Icing Research Tunnel for a variety of wet (glaze) and dry (rime) ice conditions. The chordwise temperature profiles were confirmed by observation with an IR thermal video system and thermocouple observations. The IR observations were consistent with predictions of the LEWICE ice accretion code, which was used to extrapolate the observations to rotor icing conditions. Based on the static observations, the passive IR ice detection technique appears promising; however, further testing or rotating blades is required.

Dershowitz, Adam L.↗

Computational Aerodynamic Analysis of Three-Dimensional Ice Shapes on a NACA 23012 Airfoil

The present study identifies a process for performing computational fluid dynamic calculations of the flow over full three-dimensional (3D) representations of complex ice shapes deposited on aircraft surfaces. Rime and glaze icing geometries formed on a NACA23012 airfoil were obtained during testing in the NASA Glenn Research Center's Icing Research Tunnel (IRT). The ice shape geometries were scanned as a cloud of data points using a 3D laser scanner. The data point clouds were meshed using Geomagic software to create highly accurate models of the ice surface. The surface data was imported into Pointwise grid generation software to create the CFD surface and volume grids. It was determined that generating grids in Pointwise for complex 3D icing geometries was possible using various techniques that depended on the ice shape. Computations of the flow fields over these ice shapes were performed using the NASA National Combustion Code (NCC). Results for a rime ice shape for angle of attack conditions ranging from 0 to 10 degrees and for freestream Mach numbers of 0.10 and 0.18 are presented. For validation of the computational results, comparisons were made to test results from rapid-prototype models of the selected ice accretion shapes, obtained from a separate study in a subsonic wind tunnel at the University of Illinois at Urbana-Champaign. The computational and experimental results were compared for values of pressure coefficient and lift. Initial results show fairly good agreement for rime ice accretion simulations across the range of conditions examined. The glaze ice results are promising but require some further examination.

Aerodynamics↗

Computational Aerodynamic Analysis of Three-Dimensional Ice Shapes on a NACA 23012 Airfoil

The present study identifies a process for performing computational fluid dynamic calculations of the flow over full three-dimensional (3D) representations of complex ice shapes deposited on aircraft surfaces. Rime and glaze icing geometries formed on a NACA23012 airfoil were obtained during testing in the NASA Glenn Research Centers Icing Research Tunnel (IRT). The ice shape geometries were scanned as a cloud of data points using a 3D laser scanner. The data point clouds were meshed using Geomagic software to create highly accurate models of the ice surface. The surface data was imported into Pointwise grid generation software to create the CFD surface and volume grids. It was determined that generating grids in Pointwise for complex 3D icing geometries was possible using various techniques that depended on the ice shape. Computations of the flow fields over these ice shapes were performed using the NASA National Combustion Code (NCC). Results for a rime ice shape for angle of attack conditions ranging from 0 to 10 degrees and for freestream Mach numbers of 0.10 and 0.18 are presented. For validation of the computational results, comparisons were made to test results from rapid-prototype models of the selected ice accretion shapes, obtained from a separate study in a subsonic wind tunnel at the University of Illinois at Urbana-Champaign. The computational and experimental results were compared for values of pressure coefficient and lift. Initial results show fairly good agreement for rime ice accretion simulations across the range of conditions examined. The glaze ice results are promising but require some further examination.

Aerodynamics↗

User's manual for the NASA Lewis ice accretion/heat transfer prediction code with electrothermal deicer input

A version of LEWICE has been developed that incorporates a recently developed electrothermal deicer code, developed at the University of Toledo by William B. Wright. This was accomplished, in essence, by replacing a subroutine in LEWICE, called EBAL, which balanced the energies at the ice surface, with a subroutine called UTICE. UTICE performs this same energy balance, as well as handles all the time-timperature transients below the ice surface, for all of the layers of a composite blade as well as the ice layer itself. This new addition is set up in such a fashion that a user may specify any number of heaters, any heater chordwise length, and any heater gap desired. The heaters may be fired in unison, or they may be cycled with periods independent of each other. The heater intensity may also be varied. In addition, the user may specify any number of layers and thicknesses depthwise into the blade. Thus, the new addition has maximum flexibility in modeling virtually any electrothermal deicer installed into any airfoil. It should be noted that the model simulates both shedding and runback. With the runback capability, it can simulate the anti-icing mode of heater performance, as well as detect icing downstream of the heaters due to runback in unprotected portions of the airfoil. This version of LEWICE can be run in three modes. In mode 1, no conduction heat transfer is modeled (which would be equivalent to the original version of LEWICE). In mode 2, all heat transfer is considered due to conduction but no heaters are firing. In mode 3, conduction heat transfer where the heaters are engaged is modeled, with subsequent ice shedding. When run in the first mode, there is virtually identical agreement with the original version of LEWICE in the prediction of accreted ice shapes. The code may be run in the second mode to determine the effects of conduction on the ice accretion process.

Masiulaniec, Konstanty C.↗

Academic Airframe Icing Perspective

2-D ice accretion and aerodynamics reasonably well understood for engineering applications To significantly improve our current capabilities we need to understand 3-D: a) Important ice accretion physics and modeling not well understood in 3-D; and b) Aerodynamics unsteady and 3-D especially near stall. Larger systems issues important and require multidisciplinary team approach

Bragg, Mike↗

Experimental ice shape and performance characteristics for a multi-element airfoil in the NASA Lewis Icing Research Tunnel

A study of the ice accretion patterns and performance of characteristics of a multi-element airfoil was undertaken at the NASA-Lewis Icing Research Tunnel. Several configurations were examined to determine the ice shape and performance characteristics. The testing included glaze, rime, and mixed icing regimes. Tunnel cloud conditions were set to correspond to those typical of the operating environment for commercial transport aircraft. Measurements acquired included ice profile tracings and aerodynamic forces both during the accretion process and in a post-accretion evaluation over a range of angle of attack. Substantial ice accretions developed on the main wing, flaps, and slat surfaces. Force measurements indicate severe performance degradation, especially near CL max, for both light and heavy ice accretion. Frost was seen on the lower surface of the airfoil which was found to contribute significantly to the force components.

Berkowitz, Brian M.↗

LEWICE droplet trajectory calculations on a parallel computer

A parallel computer implementation (128 processors) of LEWICE, a NASA Lewis code used to predict the time-dependent ice accretion process for two-dimensional aerodynamic bodies of simple geometries, is described. Two-dimensional parallel droplet trajectory calculations are performed to demonstrate the potential benefits of applying parallel processing to ice accretion analysis. Parallel performance is evaluated as a function of the number of trajectories and the number of processors. For comparison, similar trajectory calculations are performed on single-processor Cray computers, and the best parallel results are found to be 33 and 23 times faster, respectively, than those of the Cray XMP and YMP.

Caruso, Steven C.↗

Aerodynamics of a finite wing with simulated ice

The effect of a simulated glaze ice accretion on the aerodynamic performance of a three-dimensional wing is studied experimentally. Results are reviewed from earlier two-dimensional tests which show the character of the large leading-edge separation bubbles caused by the simulated ice accretion. The 2-D bubbles are found to closely resemble well known airfoil laminar separation bubbles. For the 3-D experiments a semispan wing of effective aspect ratio five was mounted from the sidewall of the UIUC subsonic wind tunnel. The model uses a NACA 0012 airfoil section on a rectangular planform with interchangeable tip and root sections to allow for 0- and 30-degree sweep. A three-component sidewall balance was used to measure lift, drag and pitching moment on the clean and iced model. Fluorescent oil flow visualization has been performed on the iced model and reveals extensive spanwise and vortical flow in the separation bubble aft of the upper surface horn. Sidewall interaction and spanwise nonuniformity are also seen on the unswept model. Comparisons to the computed flow fields are shown. Results are also shown for roughness effects on the straight wing. Sand grain roughness on the ice shape is seen to have a different effect than isolated 3-D roughness elements.

Bragg, M. B.↗

The Effect of Ice Formations on Propeller Performance

Measurements of propeller efficiency loss due to ice formation are supplemented by an analysis to establish the magnitude of efficiency losses to be anticipated during flight in icing conditions. The measurements were made during flight in natural icing conditions; whereas the analysis consisted of an investIgation of changes in blade-section aerodynamic characteristics caused by ice formation and the resulting propeller efficiency changes. Agreement in the order of magnitude of eff 1- ciency losses to be expected is obtained between measured and analytical results. The results indicate that, in general, efficiency losses can be expected to be less than 10 percent; whereas maximum losses, which will be encountered only rarely, may be as high as 15 or 20 percent. Reported. losses larger than 15 or 20 percent, based on reductions in airplane performance, probably are due to ice accretions on other parts of the airplane. Blade-element theory is used in the analytical treatment, and calculations are made to show the degree to which the aerodynamic characteristics of a blade section. must be altered to produce various propeller efficiency losses. The effects of ice accretions on airfoil-section characteristics at subcritical speeds and their influence on drag-divergence Mach number are examined, and. the attendant maximum efficiency losses are computed. The effect of kinetic heating on the radial extent of ice formation is considered, and its influence on required length of blade heating shoes is discussed. It is demonstrated how the efficiency loss resulting from an icing encounter is influenced by the decisions of the pilot in adjusting the engine and propeller controls.

Neel, C. B., Jr.↗

An experimental study of airfoil icing characteristics

The present investigation is concerned with the icing characteristics of a general aviation wing section, taking into account selected results from icing wind tunnel tests of the wing section. Drag coefficient increases of about 130 percent were measured for a 20 minute glaze icing encounter and cruise aerodynamic conditions. The drag increased by about 40 percent for a 15 minute rime icing encounter and the same aerodynamic conditions. The considered wing section was found to be less sensitive to rime and glaze ice accretions occurring for the climb aerodynamic conditions. The observed aft frost layer growths contributed significantly to the measured increase in section drag coefficient level. Comparisons of measured increases in section drag coefficient due to the primary ice accretion with predictions of NACA drag correlation indicated agreement was as good as that for the original data upon which the correlation was based.

Shaw, R. J.↗

Modeling of a Turbofan Engine with Ice Crystal Ingestion in the NASA Propulsion System Laboratory

The main focus of this study is to apply a computational tool for the flow analysis of the turbine engine that has been tested with ice crystal ingestion in the Propulsion Systems Laboratory (PSL) at NASA Glenn Research Center. The PSL has been used to test a highly instrumented Honeywell ALF502R-5A (LF11) turbofan engine at simulated altitude operating conditions. Test data analysis with an engine cycle code and a compressor flow code was conducted to determine the values of key icing parameters, that can indicate the risk of ice accretion, which can lead to engine rollback (un-commanded loss of engine thrust). The full engine aerothermodynamic performance was modeled with the Honeywell Customer Deck specifically created for the ALF502R-5A engine. The mean-line compressor flow analysis code, which includes a code that models the state of the ice crystal, was used to model the air flow through the fan-core and low pressure compressor. The results of the compressor flow analyses included calculations of the ice-water flow rate to air flow rate ratio (IWAR), the local static wet bulb temperature, and the particle melt ratio throughout the flow field. It was found that the assumed particle size had a large effect on the particle melt ratio, and on the local wet bulb temperature. In this study the particle size was varied parametrically to produce a non-zero calculated melt ratio in the exit guide vane (EGV) region of the low pressure compressor (LPC) for the data points that experienced a growth of blockage there, and a subsequent engine called rollback (CRB). At data points where the engine experienced a CRB having the lowest wet bulb temperature of 492 degrees Rankine at the EGV trailing edge, the smallest particle size that produced a non-zero melt ratio (between 3 percent - 4 percent) was on the order of 1 micron. This value of melt ratio was utilized as the target for all other subsequent data points analyzed, while the particle size was varied from 1 micron - 9.5 microns to achieve the target melt ratio. For data points that did not experience a CRB which had static wet bulb temperatures in the EGV region below 492 degrees Rankine, a non-zero melt ratio could not be achieved even with a 1 micron ice particle size. The highest value of static wet bulb temperature for data points that experienced engine CRB was 498 degrees Rankine with a particle size of 9.5 microns. Based on this study of the LF11 engine test data, the range of static wet bulb temperature at the EGV exit for engine CRB was in the narrow range of 492 degrees Rankine - 498 degrees Rankine , while the minimum value of IWAR was 0.002. The rate of blockage growth due to ice accretion and boundary layer growth was estimated by scaling from a known blockage growth rate that was determined in a previous study. These results obtained from the LF11 engine analysis formed the basis of a unique “icing wedge.”

Turbo engine↗

Turbulent dispersion of the icing cloud from spray nozzles used in icing tunnels

To correctly simulate flight in natural icing conditions, the turbulence in an icing simulator must be as low as possible. But some turbulence is required to mix the droplets from the spray nozzles and achieve an icing cloud of uniform liquid water content. The goal for any spray system is to obtain the widest possible spray cloud with the lowest possible turbulence in the test section of a icing tunnel. This investigation reports the measurement of turbulence and the three-dimensional spread of the cloud from a single spray nozzle. The task was to determine how the air turbulence and cloud width are affected by spray bars of quite different drag coefficients, by changes in the turbulence upstream of the spray, the droplet size, and the atomizing air. An ice accretion grid, located 6.3 m downstream of the single spray nozzle, was used to measure cloud spread. Both the spray bar and the grid were located in the constant velocity test section. Three spray bar shapes were tested: the short blunt spray bar used in the NASA Lewis Icing Research Tunnel, a thin 14.6 cm chord airfoil, and a 53 cm chord NACA 0012 airfoil. At the low airspeed (56 km/hr) the ice accretion pattern was axisymmetric and was not affected by the shape of the spray bar. At the high airspeed (169 km/hr) the spread was 30 percent smaller than at the low airspeed. For the widest cloud the spray bars should be located as far upstream in the low velocity plenum of the icing tunnel. Good comparison is obtained between the cloud spread data and predicitons from a two-dimensional cloud mixing computer code using the two equation turbulence (k epsilon g) model.

Marek, C. J.↗

Status of 3D Ice Shape Measurement Effort

(1) Main goal of the Airframe Icing Technical Challenge is to achieve acceptance of experimental and computational icing simulation tools -SupercooledLarge Droplet Icing (SLD) conditions -3D airframe components including swept wings; (2) It is necessary to develop suitable means of recording and archiving fully 3D descriptions of experimental ice accretion geometry; (3) Past research has shown that commercial laser scanners have the potential to be adapted to this task; and (4) A research plan has been developed to implement and validate the use of this technology for experimental ice accretions.

Lee, Sam↗

Numerical Study on the Aerodynamics of an Iced Airfoil With Scale-Resolving Simulations

Scale-resolving simulations of the NACA 23012 airfoil with horn ice accretion on the leading edge are conducted using the hybrid Reynolds-averaged Navier–Stokes/large-eddy simulation (hybrid RANS/LES) and wall-modeled large-eddy simulation (WMLES) approaches implemented in the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework. Aerodynamic results at the Reynolds number of 1.8 million show good comparison with the experimental measurements at different angles of attack from pre-stall to post-stall regimes. The pressure plateaus caused by the flow separation and the recovery of pressure inside the separation bubble around the iced leading edge are well predicted with the scale-resolving simulations when sufficient grid resolution is used around the accreted ice. The unsteadiness of the turbulent flows around the iced airfoil is also examined through the turbulent kinetic energy with the Reynolds normal stress anisotropy. Kelvin–Helmholtz instability (KHI) arises at the shear layer triggered by the upper ice horn and leads to rapid laminar-to-turbulent transition over a large range of angle of attack. With the increase of the angle of attack, the region with high turbulence intensity induced by the unstable shear layer spreads quickly over the entire upper surface of the airfoil. The coherent KHI modes from the upper and lower ice horns are extracted using the spectral proper orthogonal decomposition (SPOD) technique. The SPOD modes extracted from the upper shear layer have large-scale variations in the spanwise direction and low-rank behavior where the energy of the leading SPOD mode at each Strouhal number of the KHI largely represents the total energy when the mode number in the spanwise direction is small.

ARMD↗

Numerical Study on the Aerodynamics of an Iced Airfoil With Scale-Resolving Simulations

Scale-resolving simulations of the NACA 23012 airfoil with horn ice accretion on the leading edge are conducted using the hybrid Reynolds-averaged Navier–Stokes/large-eddy simulation (hybrid RANS/LES) and wall-modeled large-eddy simulation (WMLES) approaches implemented in the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework. Aerodynamic results at the Reynolds number of 1.8 million show good comparison with the experimental measurements at different angles of attack from pre-stall to post-stall regimes. The pressure plateaus caused by the flow separation and the recovery of pressure inside the separation bubble around the iced leading edge are well predicted with the scale-resolving simulations when sufficient grid resolution is used around the accreted ice. The unsteadiness of the turbulent flows around the iced airfoil is also examined through the turbulent kinetic energy with the Reynolds normal stress anisotropy. Kelvin–Helmholtz instability (KHI) arises at the shear layer triggered by the upper ice horn and leads to rapid laminar-to-turbulent transition over a large range of angle of attack. With the increase of the angle of attack, the region with high turbulence intensity induced by the unstable shear layer spreads quickly over the entire upper surface of the airfoil. The coherent KHI modes from the upper and lower ice horns are extracted using the spectral proper orthogonal decomposition (SPOD) technique. The SPOD modes extracted from the upper shear layer have large-scale variations in the spanwise direction and low-rank behavior where the energy of the leading SPOD mode at each Strouhal number of the KHI largely represents the total energy when the mode number in the spanwise direction is small.

ARMD↗

Development and Validation of Physics Based Models for Ice Shedding

Calculations for ice accretion and shedding are presented for a model scale rotor in hover. The aerodynamic characteristics of the rotor are first computed using a combined blade element-momentum theory. The effective angles of attack, and the local flow velocity are used within the NASA Glenn solver LEWICE to estimate the collection efficiency. The computed convection efficiency and the surface pressure distribution from a panel method within LEWICE are used to estimate the ice accretion over the rotor blades over an elapsed time interval. Finally, a force balance approach is used to establish shedding events where the centrifugal force over the ice mass exceeds the adhesive forces at the rotor surface and the cohesive forces between adjacent masses of ice. Preliminary comparisons with test data acquired at the Pennsylvania State Icing Research Tunnel are presented. Sensitivity of the ice shedding events to surface roughness, adhesive strength, cohesive strength, and ambient conditions is discussed.

Gupta, Avani↗

Icing testing in the large Modane wind-tunnel on full-scale and reduced scale models

Icing tests on full scale models of parts of aircraft (wings, tailplanes, radome) equipped with actual de-icing systems were carried out in the large Modane wind tunnel of ONERA. For studying icing on the Concorde, it was necessary to use a 1/6 scale half model. The equations governing the relevant parameter ratios to obtain reasonably good similitude water catching and ice accretion are recalled. Despite the inherent limitations of this particular kind of testing, i.e., the impossibility of duplicating both the Mach and Reynolds conditions for the main flow pattern, it is possible to obtain on a reduced scale model a reasonably good representation of icing cloud catching and of the shape of resulting ice accretion.

Charpin, F.↗