Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational Icing”

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 109 records · Page 6

Ice Crystal Icing Engine Testing in the NASA Glenn Research Center's Propulsion Systems Laboratory: Altitude Investigation

The National Aeronautics and Space Administration (NASA) conducted a full scale ice crystal icing turbofan engine test using an obsolete Allied Signal ALF502-R5 engine in the Propulsion Systems Laboratory (PSL) at NASA Glenn Research Center. The test article used was the exact engine that experienced a loss of power event after the ingestion of ice crystals while operating at high altitude during a 1997 Honeywell flight test campaign investigating the turbofan engine ice crystal icing phenomena. The test plan included test points conducted at the known flight test campaign field event pressure altitude and at various pressure altitudes ranging from low to high throughout the engine operating envelope. The test article experienced a loss of power event at each of the altitudes tested. For each pressure altitude test point conducted the ambient static temperature was predicted using a NASA engine icing risk computer model for the given ambient static pressure while maintaining the engine speed.

Altitude Testing↗

Preliminary Results From a Heavily Instrumented Engine Ice Crystal Icing Test in a Ground Based Altitude Test Facility

Preliminary results from the heavily instrumented ALF502R-5 engine test conducted in the NASA Glenn Research Center Propulsion Systems Laboratory are discussed. The effects of ice crystal icing on a full scale engine is examined and documented. This same model engine, serial number LF01, was used during the inaugural icing test in the Propulsion Systems Laboratory facility. The uncommanded reduction of thrust (rollback) events experienced by this engine in flight were simulated in the facility. Limited instrumentation was used to detect icing on the LF01 engine. Metal temperatures on the exit guide vanes and outer shroud and the load measurement were the only indicators of ice formation. The current study features a similar engine, serial number LF11, which is instrumented to characterize the cloud entering the engine, detect/characterize ice accretion, and visualize the ice accretion in the region of interest. Data were acquired at key LF01 test points and additional points that explored: icing threshold regions, low altitude, high altitude, spinner heat effects, and the influence of varying the facility and engine parameters. For each condition of interest, data were obtained from some selected variations of ice particle median volumetric diameter, total water content, fan speed, and ambient temperature. For several cases the NASA in-house engine icing risk assessment code was used to find conditions that would lead to a rollback event. This study further helped NASA develop necessary icing diagnostic instrumentation, expand the capabilities of the Propulsion Systems Laboratory, and generate a dataset that will be used to develop and validate in-house icing prediction and risk mitigation computational tools. The ice accretion on the outer shroud region was acquired by internal video cameras. The heavily instrumented engine showed good repeatability of icing responses when compared to the key LF01 test points and during day-to-day operation. Other noticeable observations are presented.

Enigine Icing↗

Preliminary Results From a Heavily Instrumented Engine Ice Crystal Icing Test in a Ground Based Altitude Test Facility

Preliminary results from the heavily instrumented ALF502R-5 engine test conducted in the NASA Glenn Research Center Propulsion Systems Laboratory are discussed. The effects of ice crystal icing on a full scale engine is examined and documented. This same model engine, serial number LF01, was used during the inaugural icing test in the Propulsion Systems Laboratory facility. The uncommanded reduction of thrust (rollback) events experienced by this engine in flight were simulated in the facility. Limited instrumentation was used to detect icing on the LF01 engine. Metal temperatures on the exit guide vanes and outer shroud and the load measurement were the only indicators of ice formation. The current study features a similar engine, serial number LF11, which is instrumented to characterize the cloud entering the engine, detect/ characterize ice accretion, and visualize the ice accretion in the region of interest. Data were acquired at key LF01 test points and additional points that explored: icing threshold regions, low altitude, high altitude, spinner heat effects, and the influence of varying the facility and engine parameters. For each condition of interest, data were obtained from some selected variations of ice particle median volumetric diameter, total water content, fan speed, and ambient temperature. For several cases the NASA in-house engine icing risk assessment code was used to find conditions that would lead to a rollback event. This study further helped NASA develop necessary icing diagnostic instrumentation, expand the capabilities of the Propulsion Systems Laboratory, and generate a dataset that will be used to develop and validate in-house icing prediction and risk mitigation computational tools. The ice accretion on the outer shroud region was acquired by internal video cameras. The heavily instrumented engine showed good repeatability of icing responses when compared to the key LF01 test points and during day-to-day operation. Other noticeable observations are presented.

turbomachinery↗

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.

Potapczuk, Mark G.↗

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.

Potapczuk, Mark G.↗

Role of Wind Tunnels and Computer Codes in the Certification and Qualification of Rotorcraft for Flight in Forecast Icing

The cost and time to certify or qualify a rotorcraft for flight in forecast icing has been a major impediment to the development of ice protection systems for helicopter rotors. Development and flight test programs for those aircraft that have achieved certification or qualification for flight in icing conditions have taken many years, and the costs have been very high. NASA, Sikorsky, and others have been conducting research into alternative means for providing information for the development of ice protection systems, and subsequent flight testing to substantiate the air-worthiness of a rotor ice protection system. Model rotor icing tests conducted in 1989 and 1993 have provided a data base for correlation of codes, and for the validation of wind tunnel icing test techniques. This paper summarizes this research, showing test and correlation trends as functions of cloud liquid water content, rotor lift, flight speed, and ambient temperature. Molds were made of several of the ice formations on the rotor blades. These molds were used to form simulated ice on the rotor blades, and the blades were then tested in a wind tunnel to determine flight performance characteristics. These simulated-ice rotor performance tests are discussed in the paper. The levels of correlation achieved and the role of these tools (codes and wind tunnel tests) in flight test planning, testing, and extension of flight data to the limits of the icing envelope are discussed. The potential application of simulated ice, the NASA LEWICE computer, the Sikorsky Generalized Rotor Performance aerodynamic computer code, and NASA Icing Research Tunnel rotor tests in a rotorcraft certification or qualification program are also discussed. The correlation of these computer codes with tunnel test data is presented, and a procedure or process to use these methods as part of a certification or qualification program is introduced.

Flemming, Robert J.↗

A theory of wave scatter from an inhomogeneous medium with a slightly rough boundary and its application to sea ice

An analytical theory of electromagnetic wave scattering from an inhomogeneous medium with a slightly rough boundary surface is formulated. The inhomogeneity in the medium is assumed to vary continuously in the vertical direction and to have a small random variation in the horizontal direction. The medium is assumed to consist of two layers. Maxwell's equations are solved by using the small perturbation method together with Fourier transform technique. The resulting differential equations are solved by using WKB and variation of parameter methods. Field amplitudes in each medium are determined by taking boundary conditions into account. The expressions for first order polarized radar backscatter cross-section are obtained. An attempt is made to apply the developed theory to compute sea ice scatter. Numerical calculations are performed for polarized radar backscatter cross-section at two frequencies, 13.3 GHz and 400 MHz. It is shown that WKB method is applicable at both of these frequencies. Theoretical results are compared with the experimental results obtained from NASA Earth Resources Program mission 126. Theoretical results and experimental results are in good agreement.

Parashar, S. K.↗

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.

Gupta, Avani↗

Numerical simulation of the flowfield over ice accretion shapes

The primary goals are directed toward the development of a numerical method for computing flow about ice accretion shapes and determining the influence of these shapes on flow degradation. It is expedient to investigate various aspects of icing independently in order to assess their contribution to the overall icing phenomena. The specific aspects to be examined include the water droplet trajectories with collection efficiencies and phase change on the surface, the flowfield about specified shapes including lift, drag, and heat transfer distribution, and surface roughness effects. The configurations computed were models of ice accretion shapes formed on a circular cylinder in the NASA Lewis Icing Research Tunnel. An existing Navier-Stokes program was modified to compute the flowfield over four shapes (2, 5, and 15 minute models of glaze ice, and a 15 minute accumulation of rime ice).

Scott, J. N.↗

The NASA aircraft icing research program

The objective of the NASA aircraft icing research program is to develop and make available to industry icing technology to support the needs and requirements for all weather aircraft designs. Research is being done for both fixed and rotary wing applications. The NASA program emphasizes technology development in two key areas: advanced ice protection concepts and icing simulation (analytical and experimental). The computer code development/validation, icing wind tunnel testing, and icing flight testing efforts which were conducted to support the icing technology development are reviewed.

Shaw, Robert J.↗

The NASA aircraft icing research program

The objective of the NASA aircraft icing research program is to develop and make available to industry icing technology to support the needs and requirements for all-weather aircraft designs. Research is being done for both fixed wing and rotary wing applications. The NASA program emphasizes technology development in two areas, advanced ice protection concepts and icing simulation. Reviewed here are the computer code development/validation, icing wind tunnel testing, and icing flight testing efforts.

Shaw, Robert J.↗

Simulation of Fluid Flow and Collection Efficiency for an SEA Multi-element Probe

Numerical simulations of fluid flow and collection efficiency for a Science Engineering Associates (SEA) multi-element probe are presented. Simulation of the flow field was produced using the Glenn-HT Navier-Stokes solver. Three-dimensional unsteady results were produced and then time averaged for the heat transfer and collection efficiency results. Three grid densities were investigated to enable an assessment of grid dependence. Simulations were completed for free stream velocities ranging from 85-135 meters per second, and free stream total pressure of 44.8 and 93.1 kilopascals (6.5 and 13.5 pounds per square inch absolute). In addition, the effect of angle of attack and yaw were investigated by including 5 degree deviations from straight for one of the flow conditions. All but one of the cases simulated a probe in isolation (i.e. in a very large domain without any support strut). One case is included which represents a probe mounted on a support strut within a finite sized wind tunnel. Collection efficiencies were generated, using the LEWICE3D code, for four spherical particle sizes, 100, 50, 20, and 5 micron in diameter. It was observed that a reduction in velocity of about 20% occurred, for all cases, as the flow entered the shroud of the probe. The reduction in velocity within the shroud is not indicative of any error in the probe measurement accuracy. Heat transfer results are presented which agree quite well with a correlation for the circular cross section heated elements. Collection efficiency results indicate a reduction in collection efficiency as particle size is reduced. The reduction with particle size is expected, however, the results tended to be lower than the previous results generated for isolated two-dimensional elements. The deviation from the two-dimensional results is more pronounced for the smaller particles and is likely due to the reduced flow within the protective shroud. As particle size increases differences between the two-dimensional and three dimensional results become negligible. Taken as a group, the total collection efficiency of the elements including the effects of the shroud has been shown to be in the range of 0.93 to 0.99 for particles above 20 microns. The 3D model has improved the estimated collection efficiency for smaller particles where errors in previous estimates were more significant.

Water Content↗

Viscous Three-Dimensional Simulation of Flow in an Axial Low Pressure Compressor at Engine Icing Operating Points

Viscous three-dimensional simulations of the Honeywell ALF502R-5 low pressure compressor (sometimes called a booster) using the NASA Glenn code GlennHT have been carried out. A total of ten simulations were produced. Five operating points are investigated, with each point run with two different wall thermal conditions. These operating points are at, or near, points where engine icing has been determined to be likely. In the future, the results of this study will be used for further analysis such as predicting collection efficiency of ice particles and ice growth rates at various locations in the compressor. A mixing plane boundary condition is used between each blade row, resulting in convergence to steady state within each blade row. The k-omega turbulence model of Wilcox, combined with viscous grid spacing near the wall on the order of one, is used to resolve the turbulent boundary layers. For each of the operating points, heat transfer coefficients are generated on the blades and walls. The heat transfer coefficients are produced by running the operating point with two different wall thermal conditions and then solving simultaneously for the heat transfer coefficient and adiabatic wall temperature at each point. Average Nusselt numbers are calculated for the most relevant surfaces. The values are seen to scale with Reynolds number to approximately a power of 0.7. Additionally, images of surface distribution of Nusselt number are presented. Qualitative comparison between the five operating points show that there is relatively little change in the character of the distribution. The dominant observed effect is that of an overall scaling, which is expected due to Reynolds number differences. One interesting aspect about the Nusselt number distribution is observed on the casing (outer diameter) downstream of the exit guide vanes (EGVs). The Nusselt number is relatively high between the pairs of EGVs, with two lower troughs downstream of each EGV trailing edge. This is of particular interest since rather complex ice shapes have been observed in that region.

compressor↗

Simulation of Fluid Flow and Collection Efficiency at Low Stokes Number

Numerical simulations of fluid flow and collection efficiency for a Science Engineering Associates (SEA) Robust probe are presented. The cases cover a range of freestream velocity from 10-135 m/s, and particle sizes from 1-200 µm. This combination produces results for Stokes numbers in the range of 0.07-1092. A blended modified Stokes (defined in the paper) is introduced which covers a range of 0.07-146. Of particular interest is the collection efficiency of the sensing element (a forward-facing half-pipe), Etot, in the low Stokes number regime. It was observed that the collection efficiency of the sensing area approached one for large Stokes numbers (>100) as expected. As blended modified Stokes number approached one, the Etot dropped to approximately 0.5. By blended modified Stokes numbers of 0.3, the Etot has dropped to essentially zero. Plotting collection versus blended modified Stokes number, the different freestream velocities collapsed quite well on to a single curve. The paper presents a discussion on the flowfield and local collection efficiency, beta, in an effort to explain the predictions. Of particular interest is what is essentially a shadow zone in plain sight. As the Stokes number is reduced, the local beta along the centerline of the half-pipe goes to zero. This occurs at Stokes numbers where the overall collection, Etot, is still appreciable. An empirical relation was developed to allow for corrected collection values at low Stokes numbers. The correlation remains within 2% for blended modified Stokes numbers above 1, and within 6% for all cases.

Instrumentation↗

Simulation of Fluid Flow and Collection Efficiency at Low Stokes Number

Numerical simulations of fluid flow and collection efficiency for a Science Engineering Associates (SEA) Robust probe are presented. The cases cover a range of freestream velocity from 10-135 m/s, and particle sizes from 1-200 µm. This combination produces results for Stokes numbers in the range of 0.07-1092. A blended modified Stokes (defined in the paper) is introduced which covers a range of 0.07-146. Of particular interest is the collection efficiency of the sensing element (a forward-facing half-pipe), Etot, in the low Stokes number regime. It was observed that the collection efficiency of the sensing area approached one for large Stokes numbers (>100) as expected. As blended modified Stokes number approached one, the Etot dropped to approximately 0.5. By blended modified Stokes numbers of 0.3, the Etot has dropped to essentially zero. Plotting collection versus blended modified Stokes number, the different freestream velocities collapsed quite well on to a single curve. The paper presents a discussion on the flowfield and local collection efficiency, beta, in an effort to explain the predictions. Of particular interest is what is essentially a shadow zone in plain sight. As the Stokes number is reduced, the local beta along the centerline of the half-pipe goes to zero. This occurs at Stokes numbers where the overall collection, Etot, is still appreciable. An empirical relation was developed to allow for corrected collection values at low Stokes numbers. The correlation remains within 2% for blended modified Stokes numbers above 1, and within 6% for all cases.

Instrumentation↗

Turbulence modeling in aircraft icing

The Icing and Cryogenic Technology Branch develops computational tools which predict ice growth on aircraft surfaces and uses existing CFD technology to evaluate the aerodynamic changes associated with such accretions. Surface roughness, transition location, and laminar, transition, or turbulent convective heat transfer all influence the ice growth process on aircraft surfaces. Turbulence modeling is a critical element within the computational tools used for both ice shape prediction and for performance degradation evaluation.

Potapczuk, Mark G.↗

Verification, Validation, and Solution Quality in Computational Physics: CFD Methods Applied to Ice Sheet Physics

Procedures and methods for veri.cation of coding algebra and for validations of models and calculations used in the aerospace computational fluid dynamics (CFD) community would be ef.cacious if used by the glacier dynamics modeling community. This paper presents some of those methods, and how they might be applied to uncertainty management supporting code veri.cation and model validation for glacier dynamics. The similarities and differences between their use in CFD analysis and the proposed application of these methods to glacier modeling are discussed. After establishing sources of uncertainty and methods for code veri.cation, the paper looks at a representative sampling of veri.cation and validation efforts that are underway in the glacier modeling community, and establishes a context for these within an overall solution quality assessment. Finally, a vision of a new information architecture and interactive scienti.c interface is introduced and advocated.

Thompson, David E.↗