Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ICE”

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 181 records · Page 10

Millennial–Scale Vulnerability of the Antarctic Ice Sheet to Regional Ice Shelf Collapse

The Antarctic Ice Sheet (AIS) remains the largest uncertainty in projections of future sea level rise. A likely climate–driven vulnerability of the AIS is thinning of floating ice shelves resulting from surface–melt–driven hydrofracture or incursion of relatively warm water into subshelf ocean cavities. The resulting melting, weakening, and potential ice shelf collapse reduces shelf buttressing effects. Upstream ice flow accelerates, causing thinning, grounding–line retreat, and potential ice sheet collapse. While high–resolution projections have been performed for localized Antarctic regions, full–continent simulations have typically been limited to low–resolution models. Here we quantify the vulnerability of the entire present–day AIS to regional ice shelf collapse on millennial timescales treating relevant ice flow dynamics at the necessary ~1–km resolution. Collapse of any of the ice shelves dynamically connected to the West Antarctic Ice Sheet (WAIS) is sufficient to trigger ice sheet collapse in marine–grounded portions of the WAIS. As a result, vulnerability elsewhere appears limited to localized responses.

54 ENVIRONMENTAL SCIENCES↗

Tidal Modulation of Ice-shelf Flow: a Viscous Model of the Ross Ice Shelf

Three stations near the calving front of the Ross Ice Shelf, Antarctica, recorded GPS data through a full spring-neap tidal cycle in November 2005. The data revealed a diurnal horizontal motion that varied both along and transverse to the long-term average velocity direction, similar to tidal signals observed in other ice shelves and ice streams. Based on its periodicity, it was hypothesized that the signal represents a flow response of the Ross Ice Shelf to the diurnal tides of the Ross Sea. To assess the influence of the tide on the ice-shelf motion, two hypotheses were developed. The first addressed the direct response of the ice shelf to tidal forcing, such as forces due to sea-surface slopes or forces due to sub-ice-shelf currents. The second involved the indirect response of ice-shelf flow to the tidal signals observed in the ice streams that source the ice shelf. A finite-element model, based on viscous creep flow, was developed to test these hypotheses, but succeeded only in falsifying both hypotheses, i.e. showing that direct tidal effects produce too small a response, and indirect tidal effects produce a response that is not smooth in time. This nullification suggests that a combination of viscous and elastic deformation is required to explain the observations.

Tidal↗

An Initial Study of the Fundamentals of Ice Crystal Icing Physics in the NASA Propulsion Systems Laboratory

This presentation shows results from an initial study of the fundamental physics of ice-crystal ice accretion using the NASA Propulsion Systems Lab (PSL). Ice accretion due to the ingestion of ice-crystals is being attributed to numerous jet-engine power-loss events. The NASA PSL is an altitude jet-engine test facility which has recently added a capability to inject ice particles into the flow. NASA is evaluating whether this facility, in addition to full-engine and motor-driven-rig tests, can be used for more fundamental ice-accretion studies that simulate the different mixed-phase icing conditions along the core flow passage of a turbo-fan engine compressor. The data from such fundamental accretion tests will be used to help develop and validate models of the accretion process. The present study utilized a NACA0012 airfoil. The mixed-phase conditions were generated by partially freezing the liquid-water droplets ejected from the spray bars. This presentation shows data regarding (1) the freeze out characteristics of the cloud, (2) changes in aerothermal conditions due to the presence of the cloud, and (3) the ice accretion characteristics observed on the airfoil model. The primary variable in this test was the PSL plenum humidity which was systematically varied for two duct-exit-plane velocities (85 and 135 ms) as well as two particle size clouds (15 and 50 m MVDi). The observed clouds ranged from fully glaciated to fully liquid, where the liquid clouds were at least partially supercooled. The air total temperature decreased at the test section when the cloud was activated due to evaporation. The ice accretions observed ranged from sharp arrow-like accretions, characteristic of ice-crystal erosion, to cases with double-horn shapes, characteristic of supercooled water accretions.

upwind schemes↗

ON THE FORMATION OF INTERSTELLAR WATER ICE: CONSTRAINTS FROM A SEARCH FOR HYDROGEN PEROXIDE ICE IN MOLECULAR CLOUDS

Recent surface chemistry experiments have shown that the hydrogenation of molecular oxygen on interstellar dust grains is a plausible formation mechanism, via hydrogen peroxide (H{sub 2}O{sub 2}), for the production of water (H{sub 2}O) ice mantles in the dense interstellar medium. Theoretical chemistry models also predict the formation of a significant abundance of H{sub 2}O{sub 2} ice in grain mantles by this route. At their upper limits, the predicted and experimental abundances are sufficiently high that H{sub 2}O{sub 2} should be detectable in molecular cloud ice spectra. To investigate this further, laboratory spectra have been obtained for H{sub 2}O{sub 2}/H{sub 2}O ice films between 2.5 and 200 {mu}m, from 10 to 180 K, containing 3%, 30%, and 97% H{sub 2}O{sub 2} ice. Integrated absorbances for all the absorption features in low-temperature H{sub 2}O{sub 2} ice have been derived from these spectra. For identifying H{sub 2}O{sub 2} ice, the key results are the presence of unique features near 3.5, 7.0, and 11.3 {mu}m. Comparing the laboratory spectra with the spectra of a group of 24 protostars and field stars, all of which have strong H{sub 2}O ice absorption bands, no absorption features are found that can definitely be identified with H{sub 2}O{sub 2} ice. In the absence of definite H{sub 2}O{sub 2} features, the H{sub 2}O{sub 2} abundance is constrained by its possible contribution to the weak absorption feature near 3.47 {mu}m found on the long-wavelength wing of the 3 {mu}m H{sub 2}O ice band. This gives an average upper limit for H{sub 2}O{sub 2}, as a percentage of H{sub 2}O, of 9% {+-} 4%. This is a strong constraint on parameters for surface chemistry experiments and dense cloud chemistry models.

79 ASTRONOMY AND ASTROPHYSICS↗

Further Examinations of Bimodal SLD Ice Accretion in the NASA Icing Research Tunnel

An ice shape database has been created to document ice accretions on a 72-inch chord NACA 23012 airfoil model resulting from an exposure to a Supercooled Large Drop (SLD) icing cloud with a bimodal drop size distribution. The ice shapes created were documented with photographs, laser scanned surface measurements over a section of the model span, and measurement of the ice mass over the same section of each accretion. The icing conditions used in the test matrix were based upon previously measured ice shapes on the same model to connect the current database to previously measured information. Ice shapes resulting from the bimodal distribution as well as from equivalent mono-modal drop size distributions were obtained and compared. Results indicate that the ice shapes resulting from the bimodal drop size distributions had higher mass and volume values than their mono-modal distribution equivalents as well as having icing limits that extended further back on the chord of the model.

Aircraft↗

Ice Shape Analysis of an eVTOL Propeller in Forward Flight at the NASA Glenn Icing Research Tunnel

Advanced Air Mobility (AAM) introduces many novel electric vertical takeoff and landing (eVTOL) aircraft configurations for which the effect of icing is not well understood. While icing computational tools have often aided in the design and certification of conventional aircraft, experimental data is needed to support the development and validation of such tools for eVTOL applications. To investigate the icing phenomenon relevant to eVTOL aircraft, NASA Glenn Research Center developed a general-purpose propeller test stand for conducting fundamental icing research on electrically driven propellers in the Icing Research Tunnel (IRT). A 10-day test entry in the IRT was completed in March, 2023, with the goal of generating ice accretions on a non-proprietary propeller geometry under well-characterized conditions. Three carbon fiber propellers of diameters 0.610, 0.711, and 0.914 m (24, 28, and 36 in) were used during the test. Various parametric sweeps of cloud and operating conditions were performed to evaluate the sensitivity of the parameters on the resulting ice accretion characteristics such as ice thickness, mass, location (impingement limits), and type (glaze vs rime). Ice shapes were reviewed for their repeatability, correctness of trends, and sensitivity to parameter variations. The data from the test is being used to support the development of the ice accretion solver, GlennICE.

Icing↗

Ice Shape Analysis of an eVTOL Propeller in Forward Flight at the NASA Glenn Icing Research Tunnel

Advanced Air Mobility (AAM) introduces many novel electric vertical takeoff and landing (eVTOL) aircraft configurations for which the effect of icing is not well understood. While icing computational tools have often aided in the design and certification of conventional aircraft, experimental data is needed to support the development and validation of such tools for eVTOL applications. To investigate the icing phenomenon relevant to eVTOL aircraft, NASA Glenn Research Center developed a general-purpose propeller test stand for conducting fundamental icing research on electrically driven propellers in the Icing Research Tunnel (IRT). A 10-day test entry in the IRT was completed in March, 2023, with the goal of generating ice accretions on a non-proprietary propeller geometry under well-characterized conditions. Three carbon fiber propellers of diameters 0.610, 0.711, and 0.914 m (24, 28, and 36 in) were used during the test. Various parametric sweeps of cloud and operating conditions were performed to evaluate the sensitivity of the parameters on the resulting ice accretion characteristics such as ice thickness, mass, location (impingement limits), and type (glaze vs rime). Ice shapes were reviewed for their repeatability, correctness of trends, and sensitivity to parameter variations. The data from the test is being used to support the development of the ice accretion solver, GlennICE.

Icing↗

SPH non-Newtonian Model for Ice Sheet and Ice Shelf Dynamics

We propose a new three-dimensional smoothed particle hydrodynamics (SPH) non-Newtonian model to study coupled ice sheet and ice shelf dynamics. Most existing ice sheet numerical models use a grid-based Eulerian approach, and are usually restricted to shallow ice sheet and ice shelf approximations of the momentum conservation equation. SPH, a fully Lagrangian particle method, solves the full momentum conservation equation. SPH method also allows modeling of free-surface flows, large material deformation, and material fragmentation without employing complex front-tracking schemes, and does not require re-meshing. As a result, SPH codes are highly scalable. Numerical accuracy of the proposed SPH model is first verified by simulating a plane shear flow with a free surface and the propagation of a blob of ice along a horizontal surface. Next, the SPH model is used to investigate the grounding line dynamics of ice sheet/shelf. The steady position of the grounding line, obtained from our SPH simulations, is in good agreement with laboratory observations for a wide range of bedrock slopes, ice-to-fluid density ratios, and flux. We examine the effect of non-Newtonian behavior of ice on the grounding line dynamics. The non-Newtonian constitutive model is based on Glen's law for a creeping flow of a polycrystalline ice. Finally, we investigate the effect of a bedrock geometry on a steady-state position of the grounding line.

ice sheet, grounding line, smoothed particle hydro↗

Dynamic Inland Propagation of Thinning Due to Ice Loss at the Margins of the Greenland Ice Sheet

Mass-balance analysis of the Greenland ice sheet based on surface elevation changes observed by the European Remote-sensing Satellite (ERS) (1992-2002) and Ice, Cloud and land Elevation Satellite (ICESat) (2003-07) indicates that the strongly increased mass loss at lower elevations (<2000 m) of the ice sheet, as observed during 2003-07, appears to induce interior ice thinning at higher elevations. In this paper, we perform a perturbation experiment with a three-dimensional anisotropic ice-flow model (AIF model) to investigate this upstream propagation. Observed thinning rates in the regions below 2000m elevation are used as perturbation inputs. The model runs with perturbation for 10 years show that the extensive mass loss at the ice-sheet margins does in fact cause interior thinning on short timescales (i.e. decadal). The modeled pattern of thinning over the ice sheet agrees with the observations, which implies that the strong mass loss since the early 2000s at low elevations has had a dynamic impact on the entire ice sheet. The modeling results also suggest that even if the large mass loss at the margins stopped, the interior ice sheet would continue thinning for 300 years and would take thousands of years for full dynamic recovery.

ice dynamics↗

Impact of a New Sea Ice Thermodynamic Formulation in the CESM2 Sea Ice Component

Abstract The sea ice component of the Community Earth System Model version 2 (CESM2) contains new “mushy‐layer” physics that simulates prognostic salinity in the sea ice, with consequent modifications to sea ice thermodynamics and the treatment of melt ponds. The changes to the sea ice model and their influence on coupled model simulations are described here. Two simulations were performed to assess the changes in the vertical thermodynamics formulation with prognostic salinity compared to a constant salinity profile. Inclusion of the mushy layer thermodynamics of Turner et al. (2013, https://doi.org/10.1002/jgrc.20171 ) in a fully coupled Earth system model produces thicker and more extensive sea ice in the Arctic, with relatively unchanged sea ice in the Antarctic compared to simulations using a constant salinity profile. While this is consistent with the findings of uncoupled ice‐ocean model studies, the role of the frazil and congelation growth is more important in fully coupled simulations. Melt pond drainage is also an important contribution to simulated ice thickness differences as also found in the uncoupled simulations of Turner and Hunke (2015; https://doi.org/10.1002/2014JC010358 ). However, it is an interaction of the ponds and the snow fraction that impacts the surface albedo and hence the top melt. The changes in the thermodynamics and resulting ice state modify the ice‐ocean‐atmosphere fluxes with impacts on the atmosphere and ocean states, particularly temperature.

Bailey, David A.↗

Subglacial Discharge Effects on Antarctic Ice‐Shelf Basal Melt and the Southern Ocean in a Global, Coupled Ocean—Sea‐Ice Model

Subglacial freshwater from beneath Antarctica enters the ocean at depth, enhancing ice-shelf melting and affecting Southern Ocean properties. To study these effects in an Antarctic-wide context, we use a continental-scale subglacial hydrology model that calculates grounding line freshwater flux for a global, coupled ocean—sea-ice model. We find that subglacial discharge impacts melt rates primarily through continental shelf temperature modification, contrasting with findings from regional studies that do not permit large-scale adjustments. The consequence is that Antarctic melt rates scale with subglacial discharge more strongly than inferred from regional studies. We also find that the addition of buoyancy at depth facilitates heat upwelling to the surface, resulting in higher sea ice volume downstream of cold ice shelves and lower sea ice volume downstream of warm ice shelves. This highlights the drawbacks of simplifications in previous global studies that deposit Antarctic meltwater at the ocean surface and find uniform ocean surface cooling and sea-ice growth. While the patterns we find are robust, we conclude that the addition of subglacial discharge at present-day rates has a small effect on basal melt rates, hydrography, and sea ice. However, stronger discharge can have significant effects and can even accelerate a shift from low to high melting for ice shelves close to such a tipping point. Finally, the importance of feedbacks between enhanced cavity overturning and continental shelf conditions poses a complication for parameterizing subglacial discharge effects on melting for ice-sheet projections that do not include a coupled ocean component.

54 ENVIRONMENTAL SCIENCES↗

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

ICE 2.0 Phase 2 Final Report In 2021, Berkeley Lab and Resource Innovations, Inc. launched the “ICE Calculator 2 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 Calculator 2 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 and 2 of the ICE Calculator 2 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 2 was sponsored by six utilities: Empire District Electric Company, Evergy Missouri, Pacific Gas & Electric, San Diego Gas & Electric, Southern California Edison, and Union Electric. Phase 1 and 2 involved 15 customer interruption cost survey activities representing a total of 30 electricity distribution service territories. ICE Calculator Version 2.0 and 2.2 Comparison This memorandum describes–at a high-level–the improvements in interruption cost estimates for the version 2.2 of the ICE Calculator (released February 2026) compared to version 2.0 (released in April 2025). Version 2.2 of the ICE Calculator corresponds to Phase 2 of the initiative, while version 2.0 corresponds to Phase 1. The improvements in version 2.2 result from both a significant increase in the number of customer responses that have been collected and the identification of seven additional factors that help estimate customer interruption costs.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Mid-Holocene Antarctic sea-ice increase driven by marine ice sheet retreat

Abstract. Over recent decades Antarctic sea-ice extent has increased, alongside widespread ice shelf thinning and freshening of waters along the Antarctic margin. In contrast, Earth system models generally simulate a decrease in sea ice. Circulation of water masses beneath large-cavity ice shelves is not included in current Earth System models and may be a driver of this phenomena. We examine a Holocene sediment core off East Antarctica that records the Neoglacial transition, the last major baseline shift of Antarctic sea ice, and part of a late-Holocene global cooling trend. We provide a multi-proxy record of Holocene glacial meltwater input, sediment transport, and sea-ice variability. Our record, supported by high-resolution ocean modelling, shows that a rapid Antarctic sea-ice increase during the mid-Holocene (∼ 4.5 ka) occurred against a backdrop of increasing glacial meltwater input and gradual climate warming. We suggest that mid-Holocene ice shelf cavity expansion led to cooling of surface waters and sea-ice growth that slowed basal ice shelf melting. Incorporating this feedback mechanism into global climate models will be important for future projections of Antarctic changes.

58 GEOSCIENCES↗

ISMIP6-based projections of ocean-forced Antarctic Ice Sheet evolution using the Community Ice Sheet Model

Abstract. The future retreat rate for marine-based regions of the Antarctic Ice Sheet is one of the largest uncertainties in sea-level projections. The Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6) aims to improve projections and quantify uncertainties by running an ensemble of ice sheet models with atmosphere and ocean forcing derived from global climate models. Here, the Community Ice Sheet Model (CISM) is used to run ISMIP6-based projections of ocean-forced Antarctic Ice Sheet evolution. Using multiple combinations of sub-ice-shelf melt parameterizations and calibrations, CISM is spun up to steady state over many millennia. During the spin-up, basal friction parameters and basin-scale thermal forcing corrections are adjusted to optimize agreement with the observed ice thickness. The model is then run forward for 550 years, from 1950–2500, applying ocean thermal forcing anomalies from six climate models. In all simulations, the ocean forcing triggers long-term retreat of the West Antarctic Ice Sheet, especially in the Filchner–Ronne and Ross sectors. Mass loss accelerates late in the 21st century and then rises steadily for several centuries without leveling off. The resulting ocean-forced sea-level rise at year 2500 varies from about 150 to 1300 mm, depending on the melt scheme and ocean forcing. Further experiments show relatively high sensitivity to the basal friction law, moderate sensitivity to grid resolution and the prescribed collapse of small ice shelves, and low sensitivity to the stress-balance approximation. The Amundsen sector exhibits threshold behavior, with modest retreat under many parameter settings but complete collapse under some combinations of low basal friction and high thermal forcing anomalies. Large uncertainties remain, as a result of parameterized sub-shelf melt rates, simplified treatments of calving and basal friction, and the lack of ice–ocean coupling.

54 ENVIRONMENTAL SCIENCES↗

Sea ice - Multiyear cycles and white ice

The multiyear thickness cycles represent one of the interesting features of the sea ice studies performed by Semtner (1976) and Washington et al. (1976) with simple thermodynamic models of sea ice. In the present article, a description is given of results which show that the insulating effect of snow on the surface of the sea ice is important in producing these multiyear cycles given the physics included in the model. However, when the formation of white ice is included, the cycles almost disappear. White ice is the ice which forms at the snow-ice interface when the snow layer becomes thick enough to depress the ice below the water level. Water infiltrates the snow by coming through the ice at leads and generally freezes there, forming white ice.

Ledley, T. S.↗

Satellite-derived ice data sets no. 2: Arctic monthly average microwave brightness temperatures and sea ice concentrations, 1973-1976

A summary data set for four years (mid 70's) of Arctic sea ice conditions is available on magnetic tape. The data include monthly and yearly averaged Nimbus 5 electrically scanning microwave radiometer (ESMR) brightness temperatures, an ice concentration parameter derived from the brightness temperatures, monthly climatological surface air temperatures, and monthly climatological sea level pressures. All data matrices are applied to 293 by 293 grids that cover a polar stereographic map enclosing the 50 deg N latitude circle. The grid size varies from about 32 X 32 km at the poles to about 28 X 28 km at 50 deg N. The ice concentration parameter is calculated assuming that the field of view contains only open water and first-year ice with an ice emissivity of 0.92. To account for the presence of multiyear ice, a nomogram is provided relating the ice concentration parameter, the total ice concentration, and the fraction of the ice cover which is multiyear ice.

Parkinson, C. L.↗

Analytical ice shape predictions for flight in natural icing conditions

LEWICE is an analytical ice prediction code that has been evaluated against icing tunnel data, but on a more limited basis against flight data. Ice shapes predicted by LEWICE is compared with experimental ice shapes accreted on the NASA Lewis Icing Research Aircraft. The flight data selected for comparison includes liquid water content recorded using a hot wire device and droplet distribution data from a laser spectrometer; the ice shape is recorded using stereo photography. The main findings are as follows: (1) An equivalent sand grain roughness correlation different from that used for LEWICE tunnel comparisons must be employed to obtain satisfactory results for flight; (2) Using this correlation and making no other changes in the code, the comparisons to ice shapes accreted in flight are in general as good as the comparisons to ice shapes accreted in the tunnel (as in the case of tunnel ice shapes, agreement is least reliable for large glaze ice shapes at high angles of attack); (3) In some cases comparisons can be somewhat improved by utilizing the code so as to take account of the variation of parameters such as liquid water content, which may vary significantly in flight.

Berkowitz, Brian M.↗