Search NASA⌕ Search

SEARCH · Search NASA

Results for “ice growth”

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 55 records · Page 3

Radar-Based Bayesian Estimation of Ice Crystal Growth Parameters within a Microphysical Model

The potential for polarimetric Doppler radar measurements to improve predictions of ice microphysical processes within an idealized model–observational framework is examined. In an effort to more rigorously constrain ice growth processes (e.g., vapor deposition) with observations of natural clouds, a novel framework is developed to compare simulated and observed radar measurements, coupling a bulk adaptive-habit model of vapor growth to a polarimetric radar forward model. Bayesian inference on key microphysical model parameters is then used, via a Markov chain Monte Carlo sampler, to estimate the probability distribution of the model parameters. The statistical formalism of this method allows for robust estimates of the optimal parameter values, along with (non-Gaussian) estimates of their uncertainty. To demonstrate this framework, observations from Department of Energy radars in the Arctic during a case of pristine ice precipitation are used to constrain vapor deposition parameters in the adaptive habit model. The resulting parameter probability distributions provide physically plausible changes in ice particle density and aspect ratio during growth. A lack of direct constraint on the number concentration produces a range of possible mean particle sizes, with the mean size inversely correlated to number concentration. Consistency is found between the estimated inherent growth ratio and independent laboratory measurements, increasing confidence in the parameter PDFs and demonstrating the effectiveness of the radar measurements in constraining the parameters. The combined Doppler and polarimetric observations produce the highest-confidence estimates of the parameter PDFs, with the Doppler measurements providing a stronger constraint for this case.

Robert S. Schrom↗

Development and test of a Microwave Ice Accretion Measurement Instrument (MIAMI)

The development of an ice accretion measurement instrument that is a highly sensitive, accurate, rugged and reliable microprocessor controlled device using low level microwave energy for non-instrusive real time measurement and recording of ice growth history, including ice thickness and accretion rate is discussed. Data is displayed and recorded digitally. New experimental data is presented, obtained with the instrument, which demonstrates its ability to measure ice growth on a two-dimensional airfoil. The device is suitable for aircraft icing protection. It may be mounted flush, non-intrusively, on any part of an aircraft skin including rotor blades and engine inlets.

Magenheim, B.↗

Can Milankovitch orbital variations initiate the growth of ice sheets in a general circulation model?

The possibility of initiating the growth of ice sheets by solar insolation variations is examined. The study is conducted using a climate model with three different orbital configurations corresponding to 116,000 and 106,000 yr before the present and a modified insolation field with greater reductions in summer insolation at high northern latitudes. Despite the reduced summer and fall insolation, the model fails to maintain snow cover through the summer at locations of suspected ice sheet initiation. The results suggest that there is a discrepancy between the model's response to Milankovitch perturbations and the geophysical evidence of ice sheet initiation. If the model results are correct, the growth of ice shown by geophysical evidence would have occurred in an extremely ablative environment, demanding a complicated strategy.

Rind, D.↗

Analysis of Experimental Ice Accretion Data and Assessment of a Thermodynamic Model During Ice Crystal Icing

This paper evaluates a thermodynamic ice crystal icing model that has been previously presented to describe the possible mechanisms of icing within the core of a turbofan jet engine. The model functions between two distinct ice accretions based on a surface energy balance: freeze-dominated icing and melt-dominated icing. Freeze-dominated icing occurs when liquid water (from melted ice crystals) freezes and accretes on a surface along with the existing ice of the impinging water and ice mass. This freeze-dominated icing is characterized as having strong adhesion to the surface. The amount of ice accretion is partially dictated by a freeze fraction, which is the fraction of impinging liquid water that freezes. Melt-dominated icing occurs as unmelted ice on a surface accumulates. This melt-dominated icing is characterized by weakly bonded surface adhesion. The amount of ice accumulation is partially dictated by a melt fraction, which is the fraction of impinging ice crystals that melts. Experimentally observed ice growth rates suggest that only a small fraction of the impinging ice remains on the surface, implying a mass loss mechanism such as splash, runback, bounce, or erosion. The fraction of mass loss must be determined in conjunction with the fraction of freezing liquid water or fraction of melting ice on an icing surface for a given ice growth rate. This mass loss parameter, however, along with the freeze fraction and melt fraction, are the only experimental parameters that are currently not measured directly. Using icing growth rates from ice crystal icing experiments, a methodology that has been previously proposed is used to determine these unknown parameters. This work takes ice accretion data from tests conducted by the National Aeronautics and Space Administration (NASA) at the Glenn Research Center in 2018 that examined the fundamental physics of ice crystal icing. This paper continues evaluation of the thermodynamic model from a previous effort, with additions to the model that account for sub-freezing temperatures that have been observed at the leading edge of the airfoil during icing. The predicted temperatures were generally in good agreement with measured temperatures. Other key findings include the total wet-bulb temperature being a good first order indicator of whether icing is freeze-dominated (sub-freezing values) or melt-dominated (above freezing). Maximum sticking efficiency values, the fraction of impinging mass that adheres to a surface, was calculated to be about 0.2, and retained this maximum value for a range of melt ratios (0.3 to 0.65 and possibly higher), which is defined as the ratio of liquid water content to total water content. Higher air velocities reduced the maximum sticking efficiency and shifted the icing regime to higher melt ratio values. Finally, the leading edge ice accretion angle was found to be related to ice growth (lower growth rates for smaller angles) and melt ratio (smaller melt ratios resulted in smaller angles, likely due to erosion effects).

Ice-Crystal Icing↗

Ice crystal growth in a dynamic thermal diffusion chamber

Ice crystals were grown in a supersaturated environment produced by a dynamic thermal diffusion chamber, which employed two horizontal plates separated by a distance of 2.5 cm. Air was circulated between and along the 1.2 m length of the plates past ice crystals which nucleated and grew from a fiber suspended vertically between the two plates. A zoom stereo microscope with a magnification which ranged from 3X to 80X and both 35 mm still photographs and 16 mm time lapse cine films taken through the microscope were used to study the variation of the shape and linear growth rate of ice crystals as a function of the ambient temperature, the ambient supersaturation, and the forced ventilation velocity. The ambient growth conditions were varied over the range of temperature 0 to -40 C, over the range of supersaturation 4% to 50% with respect to ice, and over the range of forced ventilation velocities 0 cm/s to 20 cm/s.

Keller, V. W.↗

Cirrus cloud model parameterizations: Incorporating realistic ice particle generation

Recent cirrus cloud modeling studies have involved the application of a time-dependent, two dimensional Eulerian model, with generalized cloud microphysical parameterizations drawn from experimental findings. For computing the ice versus vapor phase changes, the ice mass content is linked to the maintenance of a relative humidity with respect to ice (RHI) of 105 percent; ice growth occurs both with regard to the introduction of new particles and the growth of existing particles. In a simplified cloud model designed to investigate the basic role of various physical processes in the growth and maintenance of cirrus clouds, these parametric relations are justifiable. In comparison, the one dimensional cloud microphysical model recently applied to evaluating the nucleation and growth of ice crystals in cirrus clouds explicitly treated populations of haze and cloud droplets, and ice crystals. Although these two modeling approaches are clearly incompatible, the goal of the present numerical study is to develop a parametric treatment of new ice particle generation, on the basis of detailed microphysical model findings, for incorporation into improved cirrus growth models. For example, the relation between temperature and the relative humidity required to generate ice crystals from ammonium sulfate haze droplets, whose probability of freezing through the homogeneous nucleation mode are a combined function of time and droplet molality, volume, and temperature. As an example of this approach, the results of cloud microphysical simulations are presented showing the rather narrow domain in the temperature/humidity field where new ice crystals can be generated. The microphysical simulations point out the need for detailed CCN studies at cirrus altitudes and haze droplet measurements within cirrus clouds, but also suggest that a relatively simple treatment of ice particle generation, which includes cloud chemistry, can be incorporated into cirrus cloud growth.

Sassen, Kenneth↗

Contrasts in Sea Ice Formation and Production in the Arctic Seasonal and Perennial Ice Zones

Four years (1997-2000) of RADARSAT Geophysical Processor System (RGPS) data are used to contrast the sea ice deformation and production regionally, and in the seasonal (SIZ) and perennial (PIZ) ice zones. Ice production is of seasonal ice in openings during the winter. 3-day estimates of these quantities are provided within Lagrangian elements initially 10 km on a side. A distinct seasonal cycle is seen in both zones with these estimates highest in the late fall and with seasonal minimums in the mid-winter. Regional divergence over the winter could be up to 30%. Spatially, the highest deformation is in the SIZ north of coastal Alaska. Both ice deformation and production are higher in the SIZ: deformation-related ice production in the SIZ (approx.0.5 m) is 1.5-2.3 times that of the PIZ (approx.0.3 m) - this is connected to ice strength and thickness. Atmospheric forcing and boundary layer structure contribute to only the seasonal and interannual variability. Seasonal ice growth in ice fractures accounts for approx.25-40% of the total ice production of the Arctic Ocean. By itself, this deformation-ice production relationship could be considered a negative feedback when thickness is perturbed. However, the overall effect on ice production in the face of increasing seasonal and thinner/weaker ice coverage could be modified by: local destabilization of the water column promoting overturning of warmer water due to increased brine rejection; and, the upwelling of the pynocline associated with increased occurrence of large shear motion in sea ice.

Kwok, R.↗

Contrasts in Sea Ice Deformation and Production in the Arctic Seasonal and Perennial Ice Zones

Four years (1997-2000) of RADARSAT Geophysical Processor System (RGPS) data are used to contrast the sea ice deformation and production regionally, and in the seasonal (SIZ) and perennial (PIZ) ice zones. Ice production is of seasonal ice in openings during the winter. Three-day estimates of these quantities are provided within Lagrangian elements initially 10 km on a side. A distinct seasonal cycle is seen in both zones with these estimates highest in the late fall and with seasonal minimums in the midwinter. Regional divergence over the winter could be up to 30%. Spatially, the highest deformation is seen in the SIZ north of coastal Alaska. Both ice deformation and production are higher in the SIZ: deformation-related ice production in the SIZ (approx.0.5 m) is 1.5-2.3 times that of the PIZ (approx.0.3 m): this is connected to ice strength and thickness. Atmospheric forcing and boundary layer structure contribute to only the seasonal and interannual variability. Seasonal ice growth in ice fractures accounts for approx.25-40% of the total ice production of the Arctic Ocean. Uncertainties in these estimates are discussed. By itself, this deformation-ice production relationship could be considered a negative feedback when thickness is perturbed. However, the overall effect on ice production in the face of increasing seasonal and thinner/weaker ice coverage could be modified by local destabilization of the water column promoting overturning of warmer water due to increased brine rejection; and the upwelling of the pynocline associated with increased occurrence of large shear motion in sea ice. Divergence is shown to be negligibly correlated to cyclonic motion in summer and winter in both ice zones.

Arctic Ocean↗

Thin Sea-Ice Thickness as Inferred from Passive Microwave and In Situ Observations

Since microwave radiometric signals from sea-ice strongly reflect physical conditions of a layer near the ice surface, a relationship of brightness temperature with thickness is possible especially during the early stages of ice growth. Sea ice is most saline during formation stage and as the salinity decreases with time while at the same time the thickness of the sea ice increases, a corresponding change in the dielectric properties and hence the brightness temperature may occur. This study examines the extent to which the relationships of thickness with brightness temperature (and with emissivity) hold for thin sea-ice, approximately less than 0.2 -0.3 m, using near concurrent measurements of sea-ice thickness in the Sea of Okhotsk from a ship and passive microwave brightness temperature data from an over-flying aircraft. The results show that the brightness temperature and emissivity increase with ice thickness for the frequency range of 10-37 GHz. The relationship is more pronounced at lower frequencies and at the horizontal polarization. We also established an empirical relationship between ice thickness and salinity in the layer near the ice surface from a field experiment, which qualitatively support the idea that changes in the near-surface brine characteristics contribute to the observed thickness-brightness temperature/emissivity relationship. Our results suggest that for thin ice, passive microwave radiometric signals contain, ice thickness information which can be utilized in polar process studies.

Naoki, Kazuhiro↗

Ultrasonic techniques for aircraft ice accretion measurement

Ultrasonic pulse-echo measurements of ice growth on cylinders and airfoils exposed to both artificial (icing wind tunnel) and natural (flight) icing conditions are presented. An accuracy of + or - 0.5 mm is achieved with the present method. The ultrasonic signal characteristics associated with each of the two types of icing regimes identified, wet and dry ice growth, are discussed. Heat transfer coefficients are found to be higher in the wind tunnel environment than in flight. Results for ice growth on airfoils have also been obtained using an array of ultrasonic transducers. Icing profiles obtained during flight are compared with mechanical and stereo image measurements.

Hansman, R. John, Jr.↗

Importance of aggregation and small ice crystals in cirrus clouds, based on observations and an ice particle growth model

The 1 November 1986 FIRE I case study was used to test an ice particle growth model which predicts bimodal size spectra in cirrus clouds. The model was developed from an analytically based model which predicts the height evolution of monomodal ice particle size spectra from the measured ice water content (IWC). Size spectra from the monomodal model are represented by a gamma distribution, N(D) = N(sub o)D(exp nu)exp(-lambda D), where D = ice particle maximum dimension. The slope parameter, lambda, and the parameter N(sub o) are predicted from the IWC through the growth processes of vapor diffusion and aggregation. The model formulation is analytical, computationally efficient, and well suited for incorporation into larger models. The monomodal model has been validated against two other cirrus cloud case studies. From the monomodal size spectra, the size distributions which determine concentrations of ice particles less than about 150 mu m are predicted.

Mitchell, David L.↗

Analysis of Ice Mass Growth Over Time on the CRM65 Midspan Hybrid Model

The Aeronautics Research Mission Directorate at NASA is developing and applying tools to enable future technologies towards sustainable flight. Aircraft icing has been identified as a potential barrier to entry into service for innovative designs necessitating improvements to computational ice accretion tools. NASA is developing the Glenn Icing Computational Environment (GlennICE) to address deficiencies in the computational modeling capabilities of previously developed ice accretion solvers. To benchmark and improve the ability to model highly three-dimensional ice accretion, high quality validation data against experimental data is required. The CRM65 Midspan Hybrid geometry was previously tested at the NASA Icing Research Tunnel to generate experimental data for swept wing geometries typical for commercial transport aircraft. As a part of a 2018 icing test campaign, experimental data characterizing the relationship between ice accretion time and ice mass growth was obtained and can be leveraged for use in validation of computational tools. The desire for computational ice accretion solvers to predict ice shapes profiles accreted experimentally has often overshadowed the comparison to the mass and bulk volume of ice accreted. To address this deficiency, an analysis is presented in which GlennICE is applied to simulations of the CRM65 Midspan Hybrid model tested in the NASA Icing Research Tunnel. Results from the computational fluid dynamics simulations compared favorably to the experimental pressure coefficient data, thus validating the modeling setup. The experimental data showed excellent repeatability for the 15.0 minute accretion time. The comparisons between the experimental and computational ice mass over time showed good agreement up to 10.0 minutes after which the ice mass was underpredicted. The experimental ice mass was largely linear with some nonlinear data. The bulk volume of ice accreted experimentally compared well to GlennICE for the scanned ice shapes and mean combined cross section ice shapes, but was underpredicted for the maximum combined cross section ice shapes at longer accretion times. The experimental minimum combined cross section, mean combined cross section, and maximum combined cross section profiles when compared to GlennICE show good agreement for the mean combined cross section up to 15.0 minutes. The analyses show that with a single-shot method, GlennICE currently underpredicts the ice mass for longer accretion times, is not able to match the bulk volume of the maximum combined cross section due to dominating scallop features, and future work is required to generate a more generalized ice bulk density model.

Icing↗

Analysis of Ice Mass Growth Over Time on the CRM65 Midspan Hybrid Model

The Aeronautics Research Mission Directorate at NASA is developing and applying tools to enable future technologies towards sustainable flight. Aircraft icing has been identified as a potential barrier to entry into service for innovative designs necessitating improvements to computational ice accretion tools. NASA is developing the Glenn Icing Computational Environment (GlennICE) to address deficiencies in the computational modeling capabilities of previously developed ice accretion solvers. To benchmark and improve the ability to model highly three-dimensional ice accretion, high quality validation data against experimental data is required. The CRM65 Midspan Hybrid geometry was previously tested at the NASA Icing Research Tunnel to generate experimental data for swept wing geometries typical for commercial transport aircraft. As a part of a 2018 icing test campaign, experimental data characterizing the relationship between ice accretion time and ice mass growth was obtained and can be leveraged for use in validation of computational tools. The desire for computational ice accretion solvers to predict ice shapes profiles accreted experimentally has often overshadowed the comparison to the mass and bulk volume of ice accreted. To address this deficiency, an analysis is presented in which GlennICE is applied to simulations of the CRM65 Midspan Hybrid model tested in the NASA Icing Research Tunnel. Results from the computational fluid dynamics simulations compared favorably to the experimental pressure coefficient data, thus validating the modeling setup. The experimental data showed excellent repeatability for the 15.0 minute accretion time. The comparisons between the experimental and computational ice mass over time showed good agreement up to 10.0 minutes after which the ice mass was underpredicted. The experimental ice mass was largely linear with some nonlinear data. The bulk volume of ice accreted experimentally compared well to GlennICE for the scanned ice shapes and mean combined cross section ice shapes, but was underpredicted for the maximum combined cross section ice shapes at longer accretion times. The experimental minimum combined cross section, mean combined cross section, and maximum combined cross section profiles when compared to GlennICE show good agreement for the mean combined cross section up to 15.0 minutes. The analyses show that with a single-shot method, GlennICE currently underpredicts the ice mass for longer accretion times, is not able to match the bulk volume of the maximum combined cross section due to dominating scallop features, and future work is required to generate a more generalized ice bulk density model.

Icing↗

Fall Freeze-up of Sea Ice in the Beaufort-Chukchi Seas Using ERS-1 SAR and Buoy Data

The lowering of air temperatures below freezing in the fall indicates the end of summer melt and the onset of steady sea ice growth. The thickness and condition of ice that remains at the end of summer has ramifications for the thickness that that ice will attain at the end of the following winter. This period also designates a shifting of key fluxes from upper ocean freshening from ice melt to increased salinity from brine extraction during ice growth. This transitional period has been examined in the Beaufort and Chukchi Seas using ERS-1 SAR imagery and air temperatures from drifting buoys during 1991 and 1992. The SAR imagery is used to examine the condition and types of ice present in this period. Much of the surface melt water has drained off at this time. Air temperatures from drifting buoys coincident in time and within 100 km radius of the SAR imagery have been obtained...

Holt, B.↗

(abstract) Variations in Polarimetric Backscatter of Saline Ice Grown Under Diurnal Thermal Cycling Condition

An experiment was carried out in January 1994 at the Geophysical Research Facility in the Cold Regions Research and Engineering Laboratory. To investigate effects on polarimetric scattering signatures of sea ice growth under diurnal temperature variations, an ice sheet was grown for 2.5 days for the thickness of 10 cm and a polarimetric radar operating at C-band was used to obtain backscattering data in conjunction with ice-characterization measurements. The ice sheet was grown in the late morning of January 19, 1994. The initial growth rate was slow due to high insolation and temperature. As the air temperature dropped during the night, the growth rate increased significantly. The air temperature changed drastically from about -10(deg)C to -35(deg)C between day and night. The temperature cycle was repeated during the next day and the growth rate varied in the same manner. The surface of the ice was partially covered by frost flowers and the areal coverage increased as the ice became thicker. Throughout the ice growth duration of 2.5 days, polarimetric backscatter data were collected at roughly every centimeter of ice growth. For each set of radar measurements of saline ice, a set of calibration measurements was carried out with trihedrial corner reflectors and a metallic sphere. Measured polarimetric backscattering coefficients of the ice sheet reveal a strong correlation between radar data and temperature variations. As the temperature increased (decreased), the backscatter increased (decreased) correspondingly. From the ice-characterization data, temperatures of the air, at the ice-air interface, and in the ice layer had the same variation trend. Another interesting experimental observation is that the salinity measured as a function of ice depth from a sample of 10-cm thich ice indicated that the salinity variations had a similar cycle as the temperature; i.e., the salinity profile recorded the history of the temperature variations. Characterization data of the ice sheet are used in a theoretical model for scattering from saline ice with frost cover to explain the observed polarimetric signatures.

characterization polarimetric radar salinity tempe↗

Contrail Formation in Aircraft Wakes Using Large-Eddy Simulations

In this work we analyze the issue of the formation of condensation trails ("contrails") in the near-field of an aircraft wake. The basic configuration consists in an exhaust engine jet interacting with a wing-tip training vortex. The procedure adopted relies on a mixed Eulerian/Lagrangian two-phase flow approach; a simple micro-physics model for ice growth has been used to couple ice and vapor phases. Large eddy simulations have carried out at a realistic flight Reynolds number to evaluate the effects of turbulent mixing and wake vortex dynamics on ice-growth characteristics and vapor thermodynamic properties.

Paoli, R.↗

Dendritic Growth Velocity of Ice Crystals

The dendritic growth velocity of ice crystals is measured at undercooling levels between 5-8.5 K using a unique apparatus. The result shows that the measured velocities are smaller that the theoretical prediction based on the thermal diffusion theory with the scalling constant, Q= 0.025, which has been shown to agree with the experimental results at lower undercooling levels between 0.4 -2.5 K.

Dendritic Ice Growth velocity Undercooling Interfa↗