Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ice crystal 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

A Classification of Ice Crystal Habits Using Combined Lidar and Scanning Polarimeter Observations during the SEAC4RS Campaign

Using collocated NASA Cloud Physics Lidar (CPL) and Research Scanning Polarimeter (RSP) data from the Studies of Emissions and Atmospheric Composition, Clouds and Climate Coupling by Regional Surveys (SEAC4RS) campaign, a new observational-based method was developed which uses a K-means clustering technique to classify ice crystal habit types into seven categories: column, plates, rosettes, spheroids, and three different type of irregulars. Intercompared with the collocated SPEC, Inc., Cloud Particle Imager (CPI) data, the frequency of the detected ice crystal habits from the proposed method presented in the study agrees within 5% with the CPI-reported values for columns, irregulars, rosettes, and spheroids, with more disagreement for plates. This study suggests that a detailed ice crystal habit retrieval could be applied to combined space-based lidar and polarimeter observations such as CALIPSO and POLDER in addition to future missions such as the Aerosols, Clouds, Convection, and Precipitation (A-CCP).

Natalie Midzak↗

Update on the NASA Glenn Propulsion Systems Lab Icing and Ice Crystal Cloud Characterization (2017)

NASA Glenn's Propulsion Systems Lab, an altitude engine test facility, generates icing clouds with a spray system. While the spray system is used mostly to create ice crystal clouds (Appendix D/P), the 2017 cloud characterization effort added the requirement to produce exactly supercooled liquid clouds in Appendix C and Appendix O. Success was demonstrated to supercool the largest drops at the warmest conditions, but not freeze out the smallest drops at the coldest conditions. This paper documents primarily the total water content characterization methodology and results from an Iso-Kinetic Probe in ice crystals and Multi-Wire sensor in supercooled liquid, along with the cloud uniformity provided by light extinction tomography. Particle size distribution results from High Speed Imaging probes and a Phase Doppler Interferometer are discussed. Also, a new numerical model for tracking the thermodynamics of the air-drop interactions in PSL from the plenum toward the cloud characterization plane are noted. Both of these latter topic are more fully documented in companion papers.

Calibration↗

Fundamental Ice Crystal Accretion Physics Studies

Due to numerous engine power-loss events associated with high-altitude convective weather, ice accretion within an engine due to ice-crystal ingestion is being investigated. The National Aeronautics and Space Administration (NASA) and the National Research Council (NRC) of Canada are starting to examine the physical mechanisms of ice accretion on surfaces exposed to ice-crystal and mixed-phase conditions. In November 2010, two weeks of testing occurred at the NRC Research Altitude Facility utilizing a single wedge-type airfoil designed to facilitate fundamental studies while retaining critical features of a compressor stator blade or guide vane. The airfoil was placed in the NRC cascade wind tunnel for both aerodynamic and icing tests. Aerodynamic testing showed excellent agreement compared with CFD data on the icing pressure surface and allowed calculation of heat transfer coefficients at various airfoil locations. Icing tests were performed at Mach numbers of 0.2 to 0.3, total pressures from 93 to 45 kPa, and total temperatures from 5 to 15 C. Ice and liquid water contents ranged up to 20 and 3 grams per cubic meter, respectively. The ice appeared well adhered to the surface in the lowest pressure tests (45 kPa) and, in a particular case, showed continuous leading-edge ice growth to a thickness greater than 15 millimeters in 3 minutes. Such widespread deposits were not observed in the highest pressure tests, where the accretions were limited to a small area around the leading edge. The suction surface was typically ice-free in the tests at high pressure, but not at low pressure. The icing behavior at high and low pressure appeared to be correlated with the wet-bulb temperature, which was estimated to be above 0 C in tests at 93 kPa and below 0 C in tests at lower pressure, the latter enhanced by more evaporative cooling of water. The authors believe that the large ice accretions observed in the low pressure tests would undoubtedly cause the aerodynamic performance of a compressor component such as a stator blade to degrade significantly, and could damage downstream components if shed.

Currie, Tom↗

Fundamental Ice Crystal Accretion Physics Studies

Due to numerous engine power-loss events associated with high-altitude convective weather, ice accretion within an engine due to ice crystal ingestion is being investigated. The National Aeronautics and Space Administration (NASA) and the National Research Council (NRC) of Canada are starting to examine the physical mechanisms of ice accretion on surfaces exposed to ice-crystal and mixed-phase conditions. In November 2010, two weeks of testing occurred at the NRC Research Altitude Facility utilizing a single wedge-type airfoil designed to facilitate fundamental studies while retaining critical features of a compressor stator blade or guide vane. The airfoil was placed in the NRC cascade wind tunnel for both aerodynamic and icing tests. Aerodynamic testing showed excellent agreement compared with CFD data on the icing pressure surface and allowed calculation of heat transfer coefficients at various airfoil locations. Icing tests were performed at Mach numbers of 0.2 to 0.3, total pressures from 93 to 45 kPa, and total temperatures from 5 to 15 C. Ice and liquid water contents ranged up to 20 and 3 g/m3, respectively. The ice appeared well adhered to the surface in the lowest pressure tests (45 kPa) and, in a particular case, showed continuous leading-edge ice growth to a thickness greater than 15 mm in 3 min. Such widespread deposits were not observed in the highest pressure tests, where the accretions were limited to a small area around the leading edge. The suction surface was typically ice-free in the tests at high pressure, but not at low pressure. The icing behavior at high and low pressure appeared to be correlated with the wet-bulb temperature, which was estimated to be above 0 C in tests at 93 kPa and below 0 C in tests at lower pressure, the latter enhanced by more evaporative cooling of water. The authors believe that the large ice accretions observed in the low pressure tests would undoubtedly cause the aerodynamic performance of a compressor component such as a stator blade to degrade significantly, and could damage downstream components if shed.

Struk, Peter M.↗

A Flexible Parameterization for Shortwave Optical Properties of Ice Crystals

A parameterization is presented that provides extinction cross section sigma (sub e), single-scattering albedo omega, and asymmetry parameter (g) of ice crystals for any combination of volume, projected area, aspect ratio, and crystal distortion at any wavelength in the shortwave. Similar to previous parameterizations, the scheme makes use of geometric optics approximations and the observation that optical properties of complex, aggregated ice crystals can be well approximated by those of single hexagonal crystals with varying size, aspect ratio, and distortion levels. In the standard geometric optics implementation used here, sigma (sub e) is always twice the particle projected area. It is shown that omega is largely determined by the newly defined absorption size parameter and the particle aspect ratio. These dependences are parameterized using a combination of exponential, lognormal, and polynomial functions. The variation of (g) with aspect ratio and crystal distortion is parameterized for one reference wavelength using a combination of several polynomials. The dependences of g on refractive index and omega are investigated and factors are determined to scale the parameterized (g) to provide values appropriate for other wavelengths. The parameterization scheme consists of only 88 coefficients. The scheme is tested for a large variety of hexagonal crystals in several wavelength bands from 0.2 to 4 micron, revealing absolute differences with reference calculations of omega and (g) that are both generally below 0.015. Over a large variety of cloud conditions, the resulting root-mean-squared differences with reference calculations of cloud reflectance, transmittance, and absorptance are 1.4%, 1.1%, and 3.4%, respectively. Some practical applications of the parameterization in atmospheric models are highlighted.

parameterization↗

Charge Transfer Process During Collision of Riming Graupel Pellet with Small Ice Crystals within a Thundercloud

A charge transfer process during the collision of a riming graupel pellet and an ice-crystal at low temperature is proposed. During riming, the surface structure of graupel deviates from perfect crystalline structure. A concept of quasi-solid layer (QSL) formation on the surface is introduced. This QSL contains defects formed during riming. In absence of impurities, positively charged X-defect abundance is considered in the outer layer. These defects are assumed to be the charge carriers during the charge transfer process. Some part of the QSL is stripped off by the colliding ice crystals, which thereby gain some positive charge, leaving the graupel pellet negatively charged. With the proposed model, fC to pC of charge transfer is observed per collision. A transition temperature between -10 C to -15 C is also noted beyond which the QSL concept does not hold. This transition temperature is dependent on the bulk liquid water content of the cloud.

Datta, Saswati↗

Update on the NASA Glenn PSL Ice Crystal Cloud Characterization (2016)

NASA Glenn's Propulsion Systems Laboratory (PSL) is an altitude engine research test facility capable of producing ice-crystal and supercooled liquid clouds. The cloud characterization parameter space is fairly large and complex, but the phase of the cloud seems primarily governed by wet bulb temperature. The presentation will discuss some of the issues uncovered through four cloud characterization efforts to date, as well as some of instrumentation that has been used to characterize cloud parameters including cloud uniformity, bulk total water content, median volumetric diameter and max-diameter, percent freeze-out, relative humidity, and an update on the NASA Glenn PSL Ice Crystal Cloud Characterization (2016).

Calibration↗

General equations for the motions of ice crystals and water drops in gravitational and electric fields

General equations for the Reynolds number of a variety of types of ice crystals and water drops are given in terms of the Davies, Bond, and Knudsen numbers. The equations are in terms of the basic physical parameters of the system and are valid for calculating velocities in gravitational and electric fields over a very wide range of sizes and atmospheric conditions. The equations are asymptotically matched at the bottom and top of the size spectrum, useful when checking large computer codes. A numerical system for specifying the dimensional properties of ice crystals is introduced. Within the limits imposed by such variables as particle density, which have large deviations, the accuracy of velocities appears to be within 10 percent over the entire range of sizes of interest.

Nisbet, John S.↗

General equations for the motions of ice crystals and water drops in gravitational and electric fields

General equations for the Reynolds number of a variety of types of ice crystals and water drops are given in terms of the Davies, Bond, and Knudsen numbers. The equations are in terms of the basic physical parameters of the system and are valid for calculating velocities in gravitational and electric fields over a very wide range of sizes and atmospheric conditions. The equations are asymptotically matched at the bottom and top of the size spectrum, useful when checking large computer codes. A numerical system for specifying the dimensional properties of ice crystals is introduced. Within the limits imposed by such variables as particle density, which have large deviations, the accuracy of velocities appears to be within 10 percent over the entire range of sizes of interest.

Nisbet, John S.↗

Molecular-resolution imaging of ice crystallized from liquid water by cryogenic liquid-cell TEM

Despite the ubiquity of ice, a molecular-resolution image of nanoscopic defects or microstructures in ice crystallized from liquid water has never been obtained. This is mainly due to the difficulties in preparing and preserving crystalline ice samples that can survive under high-resolution imaging conditions. Here, we report the stabilization and Å-resolution electron imaging of ice Ih crystallized from liquid water by developing cryogenic liquid-cell transmission electron microscopy (CRYOLIC-TEM). We combine lattice mapping with molecular dynamics simulations to reveal that ice formation is highly tolerant to nanoscale defects such as misoriented subdomains and trapped gas bubbles, which are stabilized by molecular-scale structural motifs. Importantly, bubble surfaces adopt low-energy nanofacets and create negligible strain fields in the surrounding crystal. These bubbles can dynamically nucleate, grow, migrate, dissolve, and coalesce under electron irradiation and be monitored in situ near a steady state. This work improves our understanding of water crystallization behaviors at a molecular spatial resolution.

Materials science↗

Simulation of Fluid Flow and Collection Efficiency for a SEA Inc. Multi-Element Probe and Ice Crystal Detector Using GlennICE

Numerical simulation results of fluid flow and collection efficiency of the Science Engineering Associates Inc. Multi-Element Probe (Multiwire) and Ice Crystal Detector (ICD) are presented. Fluid flow simulations were conducted using NASA's FUN3D while collection efficiency simulations were conducted using NASA's LEWICE3D software and GlennICE software. For both probes, 3D unsteady flow results were time averaged. Simulations were computed for freestream velocities ranging from 85 to 185 m/s and freestream total pressures of 44.8 and 93.1 kPa. Collection efficiency results were computed for four spherical particle diameter sizes of 5, 20, 50, and 100 µm. GlennICE collection efficiency results for the multiwire were compared with previously published collection efficiency values calculated using LEWICE3D. Numerical collection efficiency results for the Ice Crystal Detector are presented for the first time.

aircraft, icing, Computational Fluid Dynamics↗

Simulation of Fluid Flow and Collection Efficiency for a SEA Inc. Multi-Element Probe and Ice Crystal Detector Using GlennICE

Numerical simulation results of fluid flow and collection efficiency of the Science Engineering Associates Inc. Multi-Element Probe (Multiwire) and Ice Crystal Detector (ICD) are presented. Fluid flow simulations were conducted using NASA's FUN3D while collection efficiency simulations were conducted using NASA's LEWICE3D software and GlennICE software. For both probes, 3D unsteady flow results were time-averaged. Simulations were computed for free steam velocities ranging from 85 to 185 m/s and freestream total pressures of 44.8 and 93.1 kPa. Collection efficiency results were computed for four spherical particle diameter sizes of 5, 20, 50, and 100 µm. GlennICE collection efficiency results for the multiwire were compared with previously published collection efficiency values calculated using LEWICE3D. Numerical collection efficiency results for the Ice Crystal Detector are presented for the first time.

aircraft, icing, Computational Fluid Dynamics↗

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↗

Analysis of the Honeywell Uncertified Research Engine (HURE) with Ice Crystal Cloud Ingestion at Simulated Altitudes

The Honeywell Uncertified Research Engine (HURE), a research version of a turbofan engine that never entered production, was tested in the NASA Propulsion System Laboratory (PSL), an altitude test facility at the NASA Glenn Research Center. The PSL is a facility that is equipped with water spray bars capable of producing an ice cloud consisting of ice particles, having a controlled particle diameter and concentration in the air flow. To develop the test matrix of the HURE, numerical analysis of flow and ice particle thermodynamics was performed on the compression system of the turbofan engine to predict operating conditions that could potentially result in a risk of ice accretion due to ice crystal ingestion. The goal of the test matrix was to have ice accrete in two regions of the compression system: region one, which consists of the fan-stator through the inlet guide vane (IGV), and region two which is the first stator within the high pressure compressor. The predictive analyses were performed with the mean line compressor flow modeling code (COMDES-MELT) which includes an ice particle model. The HURE engine was tested in PSL with the ice cloud over the range of operating conditions of altitude, ambient temperature, simulated flight Mach number, and fan speed with guidance from the analytical predictions. The engine was fitted with video cameras at strategic locations within the engine compression system flow path where ice was predicted to accrete, in order to visually confirm ice accretion when it occurred. In addition, traditional compressor instrumentation such as total pressure and temperature probes, static pressure taps, and metal temperature thermocouples were installed in targeted areas where the risk of ice accretion was expected. The current research focuses on the analysis of the data that was obtained after testing the HURE engine in PSL with ice crystal ingestion. The computational method (COMDES-MELT) was enhanced by computing key parameters through the fan- stator at multiple span wise locations, in order to increase the fidelity with the current mean-line method. The Icing Wedge static wet bulb temperature thresholds were applicable for determining the risk of ice accretion in the fan-stator, which is thought to be an adiabatic region. At some operating conditions near the splitter-lip region, other sources of heat (non-adiabatic walls) were suspected to be the cause of accretion, and the Icing Wedge was not applicable to predict accretion at that location. A simple order-of-magnitude heat transfer model was implemented into the COMDES-MELT code to estimate the wall temperature minimum and maximum thresholds that support ice accretion, as observed by video confirmation. The results from this model spanned the range of wall temperatures measured on a previous engine that experienced ice accretion at certain operating conditions.

Veres, Joseph P.↗

Analysis of the Honeywell Uncertified Research Engine (HURE) with Ice Crystal Cloud Ingestion at Simulated Altitudes

The Honeywell Uncertified Research Engine (HURE), a research version of a turbofan engine that never entered production, was tested in the NASA Propulsion System Laboratory (PSL), an altitude test facility at the NASA Glenn Research Center. The PSL is a facility that is equipped with water spray bars capable of producing an ice cloud consisting of ice particles, having a controlled particle diameter and concentration in the air flow. To develop the test matrix of the HURE, numerical analysis of flow and ice particle thermodynamics was performed on the compression system of the turbofan engine to predict operating conditions that could potentially result in a risk of ice accretion due to ice crystal ingestion. The goal of the test matrix was to have ice accrete in two regions of the compression system: region one, which consists of the fan-stator through the inlet guide vane (IGV), and region two which is the first stator within the high pressure compressor. The predictive analyses were performed with the mean line compressor flow modeling code (COMDES-MELT) which includes an ice particle model. The HURE engine was tested in PSL with the ice cloud over the range of operating conditions of altitude, ambient temperature, simulated flight Mach number, and fan speed with guidance from the analytical predictions. The engine was fitted with video cameras at strategic locations within the engine compression system flow path where ice was predicted to accrete, in order to visually confirm ice accretion when it occurred. In addition, traditional compressor instrumentation such as total pressure and temperature probes, static pressure taps, and metal temperature thermocouples were installed in targeted areas where the risk of ice accretion was expected. The current research focuses on the analysis of the data that was obtained after testing the HURE engine in PSL with ice crystal ingestion. The computational method (COMDES-MELT) was enhanced by computing key parameters through the fan-stator at multiple span wise locations, in order to increase the fidelity with the current mean-line method. The Icing Wedge static wet bulb temperature thresholds were applicable for determining the risk of ice accretion in the fan-stator, which is thought to be an adiabatic region. At some operating conditions near the splitter-lip region, other sources of heat (non-adiabatic walls) were suspected to be the cause of accretion, and the Icing Wedge was not applicable to predict accretion at that location. A simple order-of-magnitude heat transfer model was implemented into the COMDES-MELT code to estimate the wall temperature minimum and maximum thresholds that support ice accretion, as observed by video confirmation. The results from this model spanned the range of wall temperatures measured on a previous engine that experienced ice accretion at certain operating conditions.

Veres, Joseph P.↗

Analysis of the Honeywell Uncertified Research Engine (HURE) with Ice Crystal Cloud Ingestion at Simulated Altitudes: Public Version

The Honeywell Uncertified Research Engine (HURE), a research version of a turbofan engine that never entered production, was tested in the NASA Propulsion System Laboratory (PSL), an altitude test facility at the NASA Glenn Research Center. The PSL is a facility that is equipped with water spray bars capable of producing an ice cloud consisting of ice particles, having a controlled particle diameter and concentration in the air flow. In preparation for testing of the HURE, numerical analysis of flow and ice particle thermodynamics was performed on the compression system of the turbofan engine to predict operating conditions that could potentially result in a risk of ice accretion due to ice crystal ingestion. The results of those analyses formed the basis of the test matrix. The goal of the test matrix was to have ice accrete in two regions of the compression system: region one, which consists of the fan-stator through the inlet guide vane (IGV), and region two which is the first stator within the high pressure compressor. The predictive analyses were performed with the mean line compressor flow modeling code (COMDES-MELT) which includes an ice particle model. Together these comprise a one-dimensional icing tool. The HURE engine was tested in PSL with the ice cloud over the range of operating conditions of altitude, ambient temperature, simulated flight Mach number, and fan speed with guidance from the analytical predictions. The engine was fitted with video cameras at strategic locations within the engine compression system flow path where ice was predicted to accrete, in order to visually confirm ice accretion when it occurred. In addition, traditional compressor instrumentation such as total pressure and temperature probes, static pressure taps, and metal temperature thermocouples were installed in targeted areas where the risk of ice accretion was expected. The current research focuses on the analysis of the data that was obtained after testing the HURE engine in PSL with ice crystal ingestion. The computational method was enhanced by computing key parameters through the fan-stator at multiple spanwise locations, in order to increase the fidelity with the current mean-line method. In addition, other sources of heat (non-adiabatic walls) were suspected to be the cause of accretion near the splitter-lip and shroud. Since there were no thermocouples near the splitter, a simple order of magnitude heat transfer model was implemented to estimate the wall temperature. Future analyses will require a higher fidelity thermal analysis of the compression system metal walls to accurately calculate the total heat flux to the ice particle. For many data points analyzed, there were differences between the thermodynamic system model and the measured test data that may partially be responsible for uncertainties with the results of the current analyses.

Turbomachinery↗

Effects of Ice-Crystal Structure on Halo Formation: Cirrus Cloud Experimental and Ray-Tracing Modeling Studies

During the 1986 Project FIRE (First International Satellite Cloud Climatology Project Regional Experiment) field campaign, four 22 deg halo-producing cirrus clouds were studied jointly from a ground-based polarization lidar and an instrumented aircraft. The lidar data show the vertical cloud structure and the relative position of the aircraft, which collected a total of 84 slides by impaction, preserving the ice crystals for later microscopic examination. Although many particles were too fragile to survive impaction intact, a large fraction of the identifiable crystals were columns and radial bullet rosettes, with both displaying internal cavitations and radial plate-column combinations. Particles that were solid or displayed only a slight amount of internal structure were relatively rare, which shows that the usual model postulated by halo theorists, i.e., the randomly oriented, solid hexagonal crystal, is inappropriate for typical cirrus clouds. With the aid of new ray-tracing simulations for hexagonal hollow-ended column and bullet-rosette models, we evaluate the effects of more realistic ice-crystal structures on halo formation and lidar depolarization and consider why the common halo is not more common in cirrus clouds.

Sassen, Kenneth↗

Measurements of the aerosol and ice crystal populations in tropical stratospheric cumulonimbus anvils

A pair of particle size spectrometers was flown aboard the NASA U-2 operating from the Canal Zone to make measurements of the aerosol and ice crystal budgets in cirrus produced by thunderstorms in the tropics. Measurements indicate that fairly large crystals up to 1 mm in size are injected into the stratosphere during cumulonimbus activity. Ice water contents range from a few thousandths to a few hundredths of a gram per cubic meter. Because the ambient temperature is typically around -80 C the mass of the larger crystals largely returns to lower altitudes before evaporating. Aerosol size distributions indicate a curious narrow growth mode between 0.15 and 0.2 micron in cirrus anvils which is absent outside these clouds. The presence of this narrow mode is attributed to near water saturation produced at cloud top. Nucleation of new ice crystals as well as aerosol is hypothesized

Knollenberg, R. G.↗