Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ice crystals”

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

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.↗

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.↗

Role of small ice crystals in radiative properties of cirrus: A case study, FIRE II, November 22, 1991

Aircraft observations of cirrus cloud were made near Coffeyville, Kansas, during November 1991 as part of the First ISCCP Regional Experiment II (FIRE II) project. Cloud ice particle spectra measurements were made using both a particle measuring system (PMS) 2DC probe and an ice particle replicator. Particles larger than 200 micrometers were column rosettes. The replicator shows the presence of large numbers of ice crystals smaller than 66 micrometers (two PMS size bins) that are not recorded by the PMS 2DC probe. Calculations based on the replicator data of the geometrical blocked area and absorption cross section of the cloud per unit volume show that small particles can contribute significantly to and sometimes dominate both the solar extinction and the infrared emission. Intercomparison is made of the ice particle size, area, and mass distributions determined by these different instruments. Power law relationships for area occluded by a crystal as a function of crystal maximum dimension were computed from the PMS 2DC data. The wavelength-dependent infrared absorption cross section per volume was computed using a simple model based on anomalous diffraction and area and mass dimensional relationships for the ice crystals.

Arnott, W. Patrick↗

The Effect of Roughness Model on Scattering Properties of Ice Crystals.

We compare stochastic models of microscale surface roughness assuming uniform and Weibull distributions of crystal facet tilt angles to calculate scattering by roughened hexagonal ice crystals using the geometric optics (GO) approximation. Both distributions are determined by similar roughness parameters, while the Weibull model depends on the additional shape parameter. Calculations were performed for two visible wavelengths (864 nm and 410 nm) for roughness values between 0.2 and 0.7 and Weibull shape parameters between 0 and 1.0 for crystals with aspect ratios of 0.21, 1 and 4.8. For this range of parameters we find that, for a given roughness level, varying the Weibull shape parameter can change the asymmetry parameter by up to about 0.05. The largest effect of the shape parameter variation on the phase function is found in the backscattering region, while the degree of linear polarization is most affected at the side-scattering angles. For high roughness, scattering properties calculated using the uniform and Weibull models are in relatively close agreement for a given roughness parameter, especially when a Weibull shape parameter of 0.75 is used. For smaller roughness values, a shape parameter close to unity provides a better agreement. Notable differences are observed in the phase function over the scattering angle range from 5deg to 20deg, where the uniform roughness model produces a plateau while the Weibull model does not.

Weibull density functions↗

Derivation of Physical and Optical Properties of Midlatitude Cirrus Ice Crystals for a Size-Resolved Cloud Microphysics Model

Single-crystal images collected in mid-latitude cirrus are analyzed to provide internally consistent ice physical and optical properties for a size-resolved cloud microphysics model, including single-particle mass, projected area, fall speed, capacitance, single-scattering albedo, and asymmetry parameter. Using measurements gathered during two flights through a widespread synoptic cirrus shield, bullet rosettes are found to be the dominant identifiable habit among ice crystals with maximum dimension (Dmax) greater than 100μm. Properties are therefore first derived for bullet rosettes based on measurements of arm lengths and widths, then for aggregates of bullet rosettes and for unclassified (irregular) crystals. Derived bullet rosette masses are substantially greater than reported in existing literature, whereas measured projected areas are similar or lesser, resulting in factors of 1.5-2 greater fall speeds, and, in the limit of large Dmax, near-infrared single-scattering albedo and asymmetry parameter (g) greater by approx. 0.2 and 0.05, respectively. A model that includes commonly imaged side plane growth on bullet rosettes exhibits relatively little difference in microphysical and optical properties aside from approx. 0:05 increase in mid-visible g primarily attributable to plate aspect ratio. In parcel simulations, ice size distribution, and g are sensitive to assumed ice properties.

asymmetry↗

Formation of Large (Approximately 100 micrometers) Ice Crystals Near the Tropical Tropopause

Recent high-altitude aircraft measurements with in situ imaging instruments indicated the presence of relatively large (approx.100 microns length), thin (aspect ratios of approx.6:1 or larger) hexagonal plate ice crystals near the tropical tropopause in very low concentrations (<0.01/L). These crystals were not produced by deep convection or aggregation. We use simple growth-sedimentation calculations as well as detailed cloud simulations to evaluate the conditions required to grow the large crystals. Uncertainties in crystal aspect ratio leave a range of possibilities, which could be constrained by knowledge of the water vapor concentration in the air where the crystal growth occurred. Unfortunately, water vapor measurements made in the cloud formation region near the tropopause with different instruments ranged from <2 ppmv to approx.3.5 ppmv. The higher water vapor concentrations correspond to very large ice supersaturations (relative humidities with respect to ice of about 200%). If the aspect ratios of the hexagonal plate crystals are as small as the image analysis suggests (6:1, see companion paper (Lawson et al., 2008)) then growth of the large crystals before they sediment out of the supersaturated layer would only be possible if the water vapor concentration were on the high end of the range indicated by the different measurements (>3 ppmv). On the other hand, if the crystal aspect ratios are quite a bit larger (approx.10:1), then H2O concentrations toward the low end of the measurement range (approx.2-2.5 ppmv) would suffice to grow the large crystals. Gravity-wave driven temperature and vertical wind perturbations only slightly modify the H2O concentrations needed to grow the crystals. We find that it would not be possible to grow the large crystals with water concentrations less than 2 ppmv, even with assumptions of a very high aspect ratio of 15 and steady upward motion of 2 cm/s to loft the crystals in the tropopause region. These calculations would seem to imply that the measurements indicating water vapor concentrations less than 2ppmv are implausible, but we cannot rule out the possibility that higher humidity prevailed upstream of the aircraft measurements and the air was dehydrated by the cloud formation. Simulations of the cloud formation with a detailed model indicate that homogeneous freezing should generate ice concentrations larger than the observed concentrations (20/L), and even concentrations as low as 20/L should have depleted the vapor in excess of saturation and prevented growth of large crystals. It seems likely that the large crystals resulted from ice nucleation on effective heterogeneous nuclei at low ice supersaturations. Improvements in our understanding of detailed cloud microphysical processes require resolution of the water vapor measurement discrepancies in these very cold, dry regions of the atmosphere.

Jensen, E. J.↗

Advanced Optical Diagnostics for Ice Crystal Cloud Measurements in the NASA Glenn Propulsion Systems Laboratory

A light extinction tomography technique has been developed to monitor ice water clouds upstream of a direct connected engine in the Propulsion Systems Laboratory (PSL) at NASA Glenn Research Center (GRC). The system consists of 60 laser diodes with sheet generating optics and 120 detectors mounted around a 36-inch diameter ring. The sources are pulsed sequentially while the detectors acquire line-of-sight extinction data for each laser pulse. Using computed tomography algorithms, the extinction data are analyzed to produce a plot of the relative water content in the measurement plane. To target the low-spatial-frequency nature of ice water clouds, unique tomography algorithms were developed using filtered back-projection methods and direct inversion methods that use Gaussian basis functions. With the availability of a priori knowledge of the mean droplet size and the total water content at some point in the measurement plane, the tomography system can provide near real-time in-situ quantitative full-field total water content data at a measurement plane approximately 5 feet upstream of the engine inlet. Results from ice crystal clouds in the PSL are presented. In addition to the optical tomography technique, laser sheet imaging has also been applied in the PSL to provide planar ice cloud uniformity and relative water content data during facility calibration before the tomography system was available and also as validation data for the tomography system. A comparison between the laser sheet system and light extinction tomography resulting data are also presented. Very good agreement of imaged intensity and water content is demonstrated for both techniques. Also, comparative studies between the two techniques show excellent agreement in calculation of bulk total water content averaged over the center of the pipe.

Bencic, Timothy J.↗