Search NASA⌕ Search

SEARCH · Search NASA

Results for “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 127 records · Page 7

Ice Crystal Icing Research at NASA

Ice crystals found at high altitude near convective clouds are known to cause jet engine power-loss events. These events occur due to ice crystals entering a propulsion systems core flowpath and accreting ice resulting in events such as uncommanded loss of thrust (rollback), engine stall, surge, and damage due to ice shedding. As part of a community with a growing need to understand the underlying physics of ice crystal icing, NASA has been performing experimental efforts aimed at providing datasets that can be used to generate models to predict the ice accretion inside current and future engine designs. Fundamental icing physics studies on particle impacts, accretion on a single airfoil, and ice accretions observed during a rollback event inside a full-scale engine in the Propulsion Systems Laboratory are summarized. Low fidelity code development using the results from the engine tests which identify key parameters for ice accretion risk and the development of high fidelity codes are described. These activities have been conducted internal to NASA and through collaboration efforts with industry, academia, and other government agencies. The details of the research activities and progress made to date in addressing ice crystal icing research challenges are discussed.

icing↗

Ice Crystal Icing Research at NASA

Ice crystals found at high altitude near convective clouds are known to cause jet engine power-loss events. These events occur due to ice crystals entering a propulsion systems core flowpath and accreting ice resulting in events such as uncommanded loss of thrust (rollback), engine stall, surge, and damage due to ice shedding. As part of a community with a growing need to understand the underlying physics of ice crystal icing, NASA has been performing experimental efforts aimed at providing datasets that can be used to generate models to predict the ice accretion inside current and future engine designs. Fundamental icing physics studies on particle impacts, accretion on a single airfoil, and ice accretions observed during a rollback event inside a full-scale engine in the Propulsion Systems Laboratory are summarized. Low fidelity code development using the results from the engine tests which identify key parameters for ice accretion risk and the development of high fidelity codes are described. These activities have been conducted internal to NASA and through collaboration efforts with industry, academia, and other government agencies. The details of the research activities and progress made to date in addressing ice crystal icing research challenges are discussed.

upwind schemes↗

Predictability of the Minimum Sea Ice Extent from Winter Fram Strait Ice Area Export: Model vs Observations

Observations show predictive skill of the minimum sea ice extent (Min SIE) from late winter anomalous off-shore ice drift along the Eurasian coastline – leading to local ice thickness anomalies at the onset of the melt season – a signal then amplified by the ice-albedo feedback. We assess whether the observed seasonal predictability of September sea ice extent (Sept SIE) from Fram Strait Ice Area Export (FSIAE, a proxy for Eurasian coastal divergence) is present in Global Climate Model (GCM) large ensembles, namely the CESM2-LE, GISS-E2.1-G, FLOR-LE, CNRM-CM6-1, and CanESM5. All models show distinct periods where winter FSIAE anomalies are negatively correlated with the May sea ice thickness (May SIT) anomalies along the Eurasian coastline, and the following Sept Arctic SIE, as in observations. Counter-intuitively, several models show occasional periods where winter FSIAE anomalies are positively correlated with the following Sept SIE anomalies when the mean ice thickness is large or late in the simulation when the sea ice is thin, and/or when internal variability increases. More importantly, periods with weak correlation between FSIAE and the Sept SIE dominate, suggesting that summer melt processes dominate over late winter preconditioning and May SIT anomalies. In general, we find that the coupling between the FSIAE and ice thickness anomalies along the Eurasian coastline is an ubiquitous feature of GCMs and that the relationship with the following Sept SIE is dependent on the mean Arctic sea ice thickness.

off-shore ice drift↗

ICE Calculator 2.0: Final Report for Phase 1 of the National Initiative to Update the Interruption Cost Estimate (ICE) Calculator

In 2021, Berkeley Lab and Resource Innovations, Inc. launched the “ICE 2.0 Initiative” – a national study to refresh the underlying data and enhance the functionality of the ICE Calculator. The Initiative involves Berkeley Lab contracting with sponsoring utilities to administer identical, updated and comprehensive interruption cost surveys to statistically representative samples of each utility’s customers. Berkeley Lab and Resource Innovations then pool the survey results across the utilities and use them to update the analytical engines that drive the ICE Calculator. The ICE 2.0 Initiative is being conducted in phases. Each phase involves the administration of interruption cost surveys to the customers of sponsoring utilities, followed by an update to the ICE Calculator based on analysis of the pooled survey results. This report describes the activities and findings from Phase 1 of the ICE 2.0 Initiative. Phase 1 was sponsored by eight utilities: American Electric Power, Commonwealth Edison, Dominion Energy, Duke Energy, DTE Electric, Exelon, National Grid, and Puget Sound Energy. Phase 1 involved 11 customer interruption cost survey activities representing a total of 24 electricity distribution service territories, 23 of them located in the Eastern and Midwestern regions of the U.S. and one located in the Pacific Northwest. ICE 2.0 vs. 1.0 Comparison This memorandum compares customer power interruption costs estimated using the recently updated Interruption Cost Estimate (ICE) Calculator (“ICE 2.0”) to the original ICE Calculator (“ICE 1.0”). ICE 1.0 was developed in 2009 based on 15 independent power interruption cost surveys conducted by 10 electric utilities between 1989 and 2012. ICE 2.0 was developed in 2025 through a national initiative based on a consistent set of power interruption cost surveys and 11 surveying efforts conducted across 24 electric utility service territories between 2022 and 2024.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

The role of sea ice in 2 x CO2 climate model sensitivity. Part 1: The total influence of sea ice thickness and extent

As a first step in investigating the effects of sea ice changes on the climate sensitivity to doubled atmospheric CO2, the authors use a standard simple sea ice model while varying the sea ice distributions and thicknesses in the control run. Thinner ice amplifies the atmospheric temperature senstivity in these experiments by about 15% (to a warming of 4.8 C), because it is easier for the thinner ice to be removed as the climate warms. Thus, its impact on sensitivity is similar to that of greater sea ice extent in the control run, which provides more opportunity for sea ice reduction. An experiment with sea ice not allowed to change between the control and doubled CO2 simulations illustrates that the total effect of sea ice on surface air temperature changes, including cloud cover and water vapor feedbacks that arise in response to sea ice variations, amounts to 37% of the temperature sensitivity to the CO2 doubling, accounting for 1.56 C of the 4.17 C global warming. This is about four times larger than the sea ice impact when no feedbacks are allowed. The different experiments produce a range of results for southern high latitudes with the hydrologic budget over Antarctica implying sea level increases of varying magnitude or no change. These results highlight the importance of properly constraining the sea ice response to climate perturbations, necessitating the use of more realistic sea ice and ocean models.

Rind, D.↗

Quantification of Ice Accretions for Icing Scaling Evaluations

The comparison of ice accretion characteristics is an integral part of aircraft icing research. It is often necessary to compare an ice accretion obtained from a flight test or numerical simulation to one produced in an icing wind tunnel or for validation of an icing scaling method. Traditionally, this has been accomplished by overlaying two-dimensional tracings of ice accretion shapes. This paper addresses the basic question of how to compare ice accretions using more quantitative methods. For simplicity, geometric characteristics of the ice accretions are used for the comparison. One method evaluated is a direct comparison of the percent differences of the geometric measurements. The second method inputs these measurements into a fuzzy inference system to obtain a single measure of the goodness of the comparison. The procedures are demonstrated by comparing ice shapes obtained in the Icing Research Tunnel at NASA Glenn Research Center during recent icing scaling tests. The results demonstrate that this type of analysis is useful in quantifying the similarity of ice accretion shapes and that the procedures should be further developed by expanding the analysis to additional icing data sets.

Ruff, Gary A.↗

Advances in Measuring Antarctic Sea-Ice Thickness and Ice-Sheet Elevations with ICESat Laser Altimetry

NASA's Ice, Cloud and Land Elevation Satellite (ICESat) has been measuring elevations of the Antarctic ice sheet and sea-ice freeboard elevations with unprecedented accuracy. Since February 20,2003, data has been acquired during three periods of laser operation varying from 36 to 54 days, which is less than the continuous operation of 3 to 5 years planned for the mission. The primary purpose of ICESat is to measure time-series of ice-sheet elevation changes for determination of the present-day mass balance of the ice sheets, study of associations between observed ice changes and polar climate, and estimation of the present and future contributions of the ice sheets to global sea level rise. ICESat data will continue to be acquired for approximately 33 days periods at 3 to 6 month intervals with the second of ICESat's three lasers, and eventually with the third laser. The laser footprints are about 70 m on the surface and are spaced at 172 m along-track. The on-board GPS receiver enables radial orbit determinations to an accuracy better than 5 cm. The orbital altitude is around 600 km at an inclination of 94 degrees with a 8-day repeat pattern for the calibration and validation period, followed by a 91 -day repeat period for the rest of the mission. The expected range precision of single footprint measurements was 10 cm, but the actual range precision of the data has been shown to be much better at 2 to 3 cm. The star-tracking attitude-determination system should enable footprints to be located to 6 m horizontally when attitude calibrations are completed. With the present attitude calibration, the elevation accuracy over the ice sheets ranges from about 30 cm over the low-slope areas to about 80 cm over areas with slopes of 1 to 2 degrees, which is much better than radar altimetry. After the first period of data collection, the spacecraft attitude was controlled to point the laser beam to within 50 m of reference surface tracks over the ice sheets. Detection of ice elevation changes over select areas of the ice sheet is demonstrated with using both crossover analysis and precise-repeat track analysis. Sea ice freeboard-height distributions over the Antarctic sea pack are derived over distances of 50 km and converted into maps of average freeboard thickness and sea-ice thickness.

Zwally, H. Jay↗

Effect of the Inhomogeneity of Ice Crystals on Retrieving Ice Cloud Optical Thickness and Effective Particle Size

Spherical or spheroidal air bubbles are generally trapped in the formation of rapidly growing ice crystals. In this study the single-scattering properties of inhomogeneous ice crystals containing air bubbles are investigated. Specifically, a computational model based on an improved geometric-optics method (IGOM) has been developed to simulate the scattering of light by randomly oriented hexagonal ice crystals containing spherical or spheroidal air bubbles. A combination of the ray-tracing technique and the Monte Carlo method is used. The effect of the air bubbles within ice crystals is to smooth the phase functions, diminish the 22deg and 46deg halo peaks, and substantially reduce the backscatter relative to bubble-free particles. These features vary with the number, sizes, locations and shapes of the air bubbles within ice crystals. Moreover, the asymmetry factors of inhomogeneous ice crystals decrease as the volume of air bubbles increases. Cloud reflectance lookup tables were generated at wavelengths 0.65 m and 2.13 m with different air-bubble conditions to examine the impact of the bubbles on retrieving ice cloud optical thickness and effective particle size. The reflectances simulated for inhomogeneous ice crystals are slightly larger than those computed for homogenous ice crystals at a wavelength of 0.65 microns. Thus, the retrieved cloud optical thicknesses are reduced by employing inhomogeneous ice cloud models. At a wavelength of 2.13 microns, including air bubbles in ice cloud models may also increase the reflectance. This effect implies that the retrieved effective particle sizes for inhomogeneous ice crystals are larger than those retrieved for homogeneous ice crystals, particularly, in the case of large air bubbles.

Xie, Yu↗

Spatially Averaged Ice Contents of Ice-Wedge Polygon Cross-Sections to 3-m Depth, July 2013, Utqiagvik, Alaska

This dataset contains spatially averaged estimates of ice content (volume percent) for two-dimensional cross-sectional profiles (from trough center to trough center to a depth of 3 meters) of low- and flat-centered ice-wedge polygons (three of each type) located near Utqiagvik, Alaska. A combination of soil pits, trenches, and cores were used to describe, sample, and map the cross-section stratigraphy of soil horizons and ice wedges for each polygon at 6 depth intervals. Observed soil horizons were assigned to four types with increasing amounts of organic components (mineral, mineral/organic, organic/mineral, and organic). The average ice contents of each soil horizon type below the permafrost boundary and wedge ice were weighted by their cross-sectional area fractions to calculate spatially averaged estimates of ice content for each polygon. In the active layer, spatially averaged estimates of volumetric water contents were similarly determined and reported here as “ice” content to enable estimates of the soil’s structurally competent porosity and excess ice in permafrost layers. In this dataset, the file AK13_ice_contents.csv includes cross-sectional area fractions and spatially averaged ice contents for soil layers and horizon types at the six depth intervals. The file AK13_permafrost_ice_fractions.csv contains the calculated partitioning of ice contents into pore ice and excess ice fractions for the three permafrost-dominated depth intervals. In addition, there is a data dictionary file for each of these data files and a file-level metadata file. These data were generated by the Department of Energy’s Soil Carbon Response to Environmental Change Scientific Focus Area and were used as inputs to model simulations examining the consequences of thaw-affected subsidence and microtopography change on active layer thickness of low-relief polygonal tundra landscapes in a warming Arctic.

54 ENVIRONMENTAL SCIENCES↗

A Method for Calculating the Heat Required for Windshield Thermal Ice Prevention Based on Extensive Flight Tests in Natural Icing Conditions

An equation is presented for calculating the heat flow required from the surface of an internally heated windshield in order to prevent the formation of ice accretions during flight in specified icing conditions. To ascertain the validity of the equation, comparison is made between calculated values of the heat required and measured values obtained for test windshields in actual flights in icing conditions. The test windshields were internally heated and provided data applicable to two common types of windshield configurations; namely the V-type and the type installed flush with the fuselage contours. These windshields were installed on a twin-engine cargo airplane and the icing flights were conducted over a large area of the United States during the winters of 1945-46 and 1946-47. In addition to the internally heated windshield investigation, some test data were obtained for a windshield ice-prevention system in which heated air was discharged into the windshield boundary layer. The general conclusions resulting from this investigation are as follows: 1) The amount of heat required for the prevention of ice accretions on both flush- and V-type windshields during flight in specified icing conditions can be calculated with a degree of accuracy suitable for design purposes. 2) A heat flow of 2000 to 2500 Btu per hour per square foot is required for complete and continuous protection of a V-type windshield in fight at speeds up to 300 miles per hour in a moderate cumulus icing condition. For the same degree of protection and the same speed range, a value of 1000 Btu per hour per square foot suffices in a moderate stratus icing condition. 3) A heat supply of 1000 Btu per hour per square foot is adequate for a flush windshield located well aft of the fuselage stagnation region, at speeds up to 300 miles per hour, for flight in both stratus and moderate cumulus icing conditions. 4) The external air discharge system of windshield thermal ice prevention is thermally inefficient and requires a heat supply approximately 20 times that required for an internal system having the same performance.

ICE FORMATION↗

Basal melt beneath whillans ice stream and ice streams A and C

We have used a recently derived map of the velocity of Whillans Ice Stream and Ice Streams A and C to help estimate basal melt. Temperature was modeled with a simple vertical advection-diffusion equation, 'tuned' to match temperature profiles. We find that most of the melt occurs beneath the tributaries where larger basal shear stresses and thicker ice favors greater melt (e.g., 10-20 mm/yr). The occurrence of basal freezing is predicted beneath much of the ice plains of Ice Stream C and Whillans Ice Stream. Modelled melt rates for when Ice Stream C was active suggest there was just enough melt water generated in its tributaries to balance basal freezing on its ice plain. Net basal melt for Whillans Ice Stream is positive due to smaller basal temperature gradients. Modelled temperatures on Whillans Ice Stream, however, were constrained by a single temperature profile at UpB. Basal temperature gradients for Whillans B1 and Ice Stream A may have conditions more similar to those beneath Ice Streams C and D, in which case, there may not be sufficient melt to sustain motion. This would be consistent with the steady deceleration of Whillans stream over the last few decades.

glaciers ice sheets↗

Method to Generate Full-Span Ice Shape on Swept Wing Using Icing Tunnel Data

There is a collaborative research program by NASA, FAA, ONERA, and university partners to improve the fidelity of experimental and computational simulation methods for swept-wing ice accretion formulations and resultant aerodynamic effects on large transport aircraft. This research utilizes a 65 scale Common Research Model as the baseline configuration. In order to generate the ice shapes for the aerodynamic testing, ice-accretion testing will be conducted in the NASA Icing Research Tunnel utilizing hybrid model from the 20, 64, and 83 spanwise locations. The models will have full-scale leading edges with truncated chord in order to fit the IRT test section. The ice shapes from the IRT tests will be digitized using a commercially available articulated-arm 3D laser scanning system. The methodology to acquire 3D ice shapes using a laser scanner was developed and validated in a previous research effort. Each of these models will yield a 1.5ft span of ice than can be used. However, a full-span ice accretion will require 75 ft span of ice. This means there will be large gaps between these spanwise ice sections that must be filled, while maintaining all of the important aerodynamic features. A method was developed to generate a full-span ice shape from the three 1.5 ft span ice shapes from the three models.

swept↗

Characterization of Ice Roughness Variations in Scaled Glaze Icing Conditions

Because of the significant influence of surface tension in governing the stability and breakdown of the liquid film in flooded stagnation regions of airfoils exposed to glaze icing conditions, the Weber number is expected to be a significant parameter governing the formation and evolution of ice roughness. To investigate the influence of the Weber number on roughness formation, 53.3-cm (21-in.) and 182.9-cm (72-in.) NACA 0012 airfoils were exposed to flow conditions with essentially the same Weber number and varying stagnation collection efficiency to illuminate similarities of the ice roughness created on the different airfoils. The airfoils were exposed to icing conditions in the Icing Research Tunnel (IRT) at the NASA Glenn Research Center. Following exposure to the icing event, the airfoils were then scanned using a ROMER Absolute Arm scanning system. The resulting point clouds were then analyzed using the self-organizing map approach of McClain and Kreeger (2013) to determine the spatial roughness variations along the surfaces of the iced airfoils. The roughness characteristics on each airfoil were then compared using the relative geometries of the airfoil. The results indicate that features of the ice shape and roughness such as glaze-ice plateau limits and maximum airfoil roughness were captured well by Weber number and collection efficiency scaling of glaze icing conditions. However, secondary ice roughness features relating the instability and waviness of the liquid film on the glaze-ice plateau surface are scaled based on physics that were not captured by the local collection efficiency variations.

Roughness Formation↗

Ice Crystal Icing Research at NASA

Ice crystals found at high altitude near convective clouds are known to cause jet engine power-loss events. These events occur due to ice crystals entering a propulsion system's core flowpath and accreting ice resulting in events such as uncommanded loss of thrust (rollback), engine stall, surge, and damage due to ice shedding. As part of a community with a growing need to understand the underlying physics of ice crystal icing, NASA has been performing experimental efforts aimed at providing datasets that can be used to generate models to predict the ice accretion inside current and future engine designs. Fundamental icing physics studies on particle impacts, accretion on a single airfoil, and ice accretions observed during a rollback event inside a full-scale engine in the Propulsion Systems Laboratory are summarized. Low fidelity code development using the results from the engine tests which identify key parameters for ice accretion risk and the development of high fidelity codes are described. These activities have been conducted internal to NASA and through collaboration efforts with industry, academia, and other government agencies. The details of the research activities and progress made to date in addressing ice crystal icing research challenges are discussed.

crysta↗

Impacts of a mushy‐layer thermodynamic approach in global sea‐ice simulations using the CICE sea‐ice model

Abstract We perform global simulations of the Los Alamos sea‐ice model, CICE, with a new thermodynamics component that has a fully prognostic, variable bulk salinity vertical profile based on mushy layer physics. The processes of gravity drainage, melt‐water flushing and snow‐ice formation are parameterized to allow the bulk salinity to evolve with time. We analyze the seasonal and spatial variation of sea‐ice bulk salinity, area, volume and thickness and compare these quantities to simulations using the previous thermodynamic component. Adjusting one of the gravity drainage parameters, we find good agreement between simulation results and fieldwork ice‐core observations of sea‐ice bulk salinity. As expected, bulk salinity is highest during periods of ice growth and lowest after periods of ice melt. In the northern hemisphere the new thermodynamics component produces thicker ice than the previous thermodynamics component. Of the nine major differences between the two models, differences in how salinities are calculated and how melt‐pond flushing is parameterized are the principal causes of this thickness difference. Thickness differences are smaller in the southern hemisphere than in the northern hemisphere since a greater fraction of ice melts, and differences cannot accumulate year‐on‐year. Model differences in how ice thickness changes and snow‐ice formation are calculated are the most important causes of the different thickness between the two thermodynamic components in the southern hemisphere. The melt‐pond area and volume are found to be highly sensitive to a parameter choice controlling drainage through macroscopic holes in the ice, in both hemispheres.

Turner, Adrian K.↗

Ice nucleation activity of diesel soot particles at Cirrus relevant conditions: Effects of hydration, secondary organics coating, hydration, soot morphology, and coagulation

The role of atmospheric relevant soot particles that are processed in the atmosphere toward ice nucleation at cirrus cloud condition is poorly understood. In this study, the ice nucleating properties of diesel soot particles subjected to various physical and chemical aging treatments were investigated at temperatures ranging from -40 to -50 °C. We show that bare soot particles nucleate ice in deposition mode, but coating with secondary organics suppresses the heterogeneous ice nucleation potential of soot particles requiring homogeneous freezing threshold conditions. However, the ice nucleation efficiency of soot particles coated with an aqueous organic layer was similar to bare soot particles. Hydration of bare soot particles slightly enhanced the ice nucleation efficiency, and the IN abilities of compact soot particles (roundness = ~ 0.6) were similar to bare lacey soot particles (roundness = ~ 0.4). These results indicate that ice nucleation properties are sensitive to the various aging treatments.

Environmental Molecular Sciences Laboratory↗

On the Role of Airborne Ice Nucleating Particles in Primary and Secondary Ice Formation Processes in Convective Midlatitude Clouds

Abstract Formation of ice in mixed-phase clouds may proceed via primary nucleation and secondary ice production (SIP). Primary nucleation involves the presence of ice nucleating particles (INPs), while SIP follows primary nucleation, resulting in ice crystal number concentrations N ice exceeding the number of INPs. The ice formation pathways in young congestus clouds are not well understood nor are the environmental conditions favorable for the onset of SIP. Coincident airborne measurements of INPs and N ice collected during the Secondary Production of Ice in Cumulus Experiment (SPICULE) campaign over the U.S. central Great Plains are reported for young congestus clouds. Number concentrations of INPs and N ice at temperatures between −10° and −20°C were used to categorize whether SIP was active in each analyzed cloud pass. Observational data suggested that fragmentation of freezing droplets (FFDs) was occurring and the most likely SIP mechanism for ice enhancement during the early stages of the sampled clouds, where higher cloud-base temperatures and stronger updrafts produced conditions favorable for the onset of FFD. Simple model simulations of a congestus cloud were used to compare predictions of N ice from newly implemented SIP mechanisms (including FFD) alongside those from the existing Hallett–Mossop (HM) mechanism. Although predicted N ice , when all SIP mechanisms were active, was in good agreement with observations, the dominant SIP mechanism in the model was predicted to be HM. The rate of heterogeneous freezing in raindrops was likely inhibiting realistic FFD rates in the simulated cloud. These discrepancies between the observations and simulations underscore the need for extended laboratory studies to formulate improved model representations of SIP in convective clouds.

Patnaude, Ryan J. [Department of Atmospheric Scien↗