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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Ice Load Project Final Technical Report

As interest and investment in offshore wind projects increase worldwide, some turbines will be installed in locations where ice of significant thickness forms on the water surface. This ice moves under the driving forces of wind, current, and thermal effects and may result in substantial forces on bottom-fixed support structures. The North and Baltic Seas in Europe have begun to see significant wind energy development and the Great Lakes of the United States and Canada may host wind energy development in the near future. Design of the support structures for these projects is best performed through the use of an integrated tool that can calculate the cumulative effects of forces due to turbine operations, wind, waves, and floating ice. The dynamic nature of ice forces requires that these forces be included in the design simulations, rather than added as static forces to simulation results. The International Electrotechnical Commission (IEC) standard[2] for offshore wind turbine design and the International Organization for Standardization (ISO) standard[3] for offshore structures provide requirements and algorithms for the calculation of forces induced by surface ice; however, currently none of the major wind turbine dynamic simulation codes provides the ability to model ice loads. The scope of work of the project described in this report includes the development of a suite of subroutines, collectively named IceFloe, that meet the requirements of the IEC and ISO standards and couples with four of the major wind turbine dynamic simulation codes. The mechanisms by which ice forces impinge on offshore structures generally include the forces required for crushing of the ice against vertical-sided structures and the forces required to fracture the ice as it rides up on conical-sided structures. Within these two broad categories, the dynamic character of the forces with respect to time is also dependent on other factors such as the velocity and thickness of the moving ice and the response of the structure. In some cases, the dynamic effects are random and in other cases they are deterministic, such as the effect of structural resonance and coupling of the ice forces with the defection of the support structure. The initial versions of the IceFloe routines incorporate modules that address these varied force and dynamic phenomena with seven alternative algorithms that can be specified by the user. The IceFloe routines have been linked and tested with four major wind turbine aeroelastic simulation codes: FAST, a tool developed under the management of the National Renewable Energy Laboratory (NREL) and available free of charge from its web site; Bladed[4], a widely-used commercial package available from DNV GL; ADAMS[5], a general purpose multi-body simulation code used in the wind industry and available from MSC Software; and HAWC2[6], a code developed by and available for purchase from Danmarks Tekniske Universitet (DTU). Interface routines have been developed and tested with full wind turbine simulations for each of these codes and the source code and example inputs and outputs are available from the NREL website.

17 WIND ENERGY↗
NASA NTRS2016

Generation of Fullspan Leading-Edge 3D Ice Shapes for Swept-Wing Aerodynamic Testing

The deleterious effect of ice accretion on aircraft is often assessed through dry-air flight and wind tunnel testing with artificial ice shapes. This paper describes a method to create fullspan swept-wing artificial ice shapes from partial span ice segments acquired in the NASA Glenn Icing Reserch Tunnel for aerodynamic wind-tunnel testing. Full-scale ice accretion segments were laser scanned from the Inboard, Midspan, and Outboard wing station models of the 65% scale Common Research Model (CRM65) aircraft configuration. These were interpolated and extrapolated using a weighted averaging method to generate fullspan ice shapes from the root to the tip of the CRM65 wing. The results showed that this interpolation method was able to preserve many of the highly three dimensional features typically found on swept-wing ice accretions. The interpolated fullspan ice shapes were then scaled to fit the leading edge of a 8.9% scale version of the CRM65 wing for aerodynamic wind-tunnel testing. Reduced fidelity versions of the fullspan ice shapes were also created where most of the local three-dimensional features were removed. The fullspan artificial ice shapes and the reduced fidelity versions were manufactured using stereolithography.

Icing↗
NASA NTRS2023

Modeling Sea Ice Albedo and Transmittance Measurements With A Fully-Coupled Radiative Transfer Model

A rigorous treatment of the sea ice medium has been incorporated in the advanced Coupled Ocean-Atmosphere Radiative Transfer (COART) model. The inherent optical properties (IOPs) of brine pockets and air bubbles over the 0.25-4.0 μm spectral region are parameterized as a function of the vertical profile of the sea ice physical properties (temperature, salinity and density). We test the model performance using available albedo and transmittance measurements collected during the Impacts of Climate on the Ecosystems and Chemistry of the Arctic Pacific Environment (ICESCAPE) and the Surface Heat Budget of the Arctic Ocean (SHEBA) field campaigns. The observations are adequately simulated when at least three layers are used to represent bare (first-year and multi-year) ice, including a thin top layer characterized by low density and high scattering. Two layers can be sufficient to model isolated cases of multi-year ice, and apply well to ponded ice except for shallow ponds over thick ice. The albedo and transmittance of ponded ice in the visible are mainly determined by the optical properties of the ice underlying the water layer used to model the pond. Sensitivity results indicate that the air volume or ice density has the largest impact on the simulated fluxes. Possible contamination from light-absorbing impurities, such as black carbon or ice algae, is also implemented in the model and is able to effectively reduce the albedo and transmittance in the visible spectrum to further improve the model-observation agreement. The agreement between the observed and modeled spectra validates the parameterization of the sea ice IOPs, and endorses COART as an accurate tool for radiation studies in the cryosphere.

Coupled Ocean-9 Atmosphere Radiative Transfer (COA↗
NASA NTRS2024

Physics of the Seasonal Sea Ice Zone

The seasonal sea ice zone encompasses the region between the winter maximum and summer minimum sea ice extent. In both the Arctic and Antarctic, the majority of the ice cover can now be classified as seasonal. Here, we review the sea ice physics that governs the evolution of seasonal sea ice in the Arctic and Antarctic, spanning sea ice growth, melt, and dynamics and including interactions with ocean surface waves as well as other coupled processes. The advent of coupled wave–ice modeling and discrete-element modeling, together with improved and expanded satellite observations and field campaigns, has yielded advances in process understanding. Many topics remain in need of further investigation, including rheologies appropriate for seasonal sea ice, wave-induced sea ice fracture, welding for sea ice freeze-up, and the distribution of snow on seasonal sea ice. Future research should aim to redress biases (such as disparities in focus between the Arctic and Antarctic and between summer and winter processes) and connect observations to modeling across spatial scales.

sea ice↗
NASA NTRS2020

COMDES-MELT: A Turbofan Engine Icing Risk Analysis Tool, User’s Manual

The computational tool COMDES-MELT was developed to predict the susceptibility of turbofan engines to ice accretion when flying in an ice crystal environment. COMDES-MELT is a first-generation computational tool that can estimate the conditions at which the accretion of ice can occur, location within the turbofan engine compression system, and at which engine operating points within the flight trajectory. This is accomplished by tracking the parameters of static wet bulb temperature, ice particle melt ratio, and the ice-water to airflow rates ratio as ice crystals are ingested into the engine compression system. It has been used successfully to predict engine operating points for simulation in an altitude wind tunnel with ice crystal ingestion. These engine tests occur in a simulated ice crystal cloud environment. This predictive analysis has been used to generate the test matrix in preparation for tests in the Propulsion Systems Laboratory, an altitude test facility at the NASA Glenn Research Center. Furthermore, due to the speed of the code during an engine icing test, adjustments to the test matrix can be provided to the test engineer in real time to focus on a particular area of ice accretion susceptibility. The use of the COMDES-MELT code will be examined in this manual.

Turbofan engine↗
DOE OSTI2011

Effects of ice number concentration on dynamics of a shallow mixed-phase stratiform cloud

Previous modeling studies have shown a high sensitivity of simulated properties of mixed-phase clouds to ice number concentration, Ni, with many models losing their ability to maintain the liquid phase as Ni increases. Although models differ widely at what Ni the mixed-phase cloud becomes unstable, the transition from a mixed-phase to an ice only cloud in many cases occurs over a narrow range of ice concentration. To gain better understanding of this non-linear model behavior, in this study, we analyze simulations of a mixed-phase stratiform Artic cloud observed on 26 April 2008 during recent Indirect and Semi-Direct Aerosol Campaign (ISDAC). The BASE simulation, in which Ni is constrained to match the measured value, produces a long-lived cloud in a quasi steady state similar to that observed. The simulation without the ice (NO_ICE) produces a comparable but slightly thicker cloud because more moisture is kept in the mixed layer due to lack of precipitation. When Ni is quadrupled relative to BASE (HI_ICE), the cloud starts loosing liquid water almost immediately and the liquid water path is reduced by half in less than two hours. The changes in liquid water are accompanied by corresponding reduction in the radiative cooling of the layer and a slow down in the vertical mixing, confirming the important role of interactions among microphysics, radiation and dynamics in this type of clouds. Deviations of BASE and HI_ICE from NO_ICE are used to explore the linearity of the model response to variation in Ni. It is shown that at early stages, changes in liquid and ice water as well as in radiative cooling/heating rates are proportional to the Ni change, while changes in the vertical buoyancy flux are qualitatively different in HI_ICE compared to BASE. Thus, while the positive feedback between the liquid water path and radiative cooling of the cloud layer is essential for glaciation of the cloud at higher Ni, the non-linear (with respect to Ni) reduction in positive buoyancy flux within and below the liquid cloud layer emerges as the process that influences the initial response of the mixed layer dynamics to the appearance of ice and subsequently determines the sustainability of liquid water in the cloud in this case.

58 GEOSCIENCES↗
DOE OSTI2021

Unraveling the Combined Effects of Ice and Permafrost on Arctic Delta Morphodynamics

Arctic delta channel dynamics and coastal deposition patterns are uniquely affected by ice and permafrost. Building on previous experiments that documented the individual effects of ice and permafrost, we conducted new numerical modeling experiments to examine the combined influences of ice thickness and permafrost erodibility. We quantified the effects of ice, permafrost, and their combination on channel migration, avulsion, and persistence; onshore and nearshore deposition; and large-scale delta morphology. Both ice and permafrost limited channel mobility, enhanced overbank deposition, and roughened shorelines. Thick ice increased offshore channel incision and sediment delivery. Permafrost stabilized channels and focused deposition that led to enhanced channel levee development. The increased levee deposition limited widespread channel infilling leading to the preservation of abandoned channels. Levee growth was enhanced by overbank flooding due to ice-driven changes in water surface elevations. Furthermore, this flooding led to more extensive levees on deltas with both ice and permafrost compared to deltas with only permafrost. Ice also preserved small channels, leading to size-dependent channel persistence: the presence of ice in the channels reduced the fraction of time small channels remained active relative to larger channels. Resistant permafrost damped this effect, as smaller and more distal channels were active earlier in the flood on deltas with resistant permafrost likely due to the smaller number of channels on permafrost-laden deltas. Finally, our results suggest that the loss of both ice and permafrost will increase channel mobility, decrease overbank flooding, increase shoreline progradation rates, and decrease offshore transport of riverine fluxes.

58 GEOSCIENCES↗
DOE OSTI2023

Assessing the Robustness of Arctic Sea Ice Bi‐Stability in the Presence of Atmospheric Feedbacks

Abstract Arctic sea‐ice loss is influenced by multiple positive feedbacks, sparking concerns of accelerated loss in the coming years or even a tipping point, where a sea‐ice equilibrium disappears at a given CO 2 value and sea ice rapidly evolves to a new steady state. Such a tipping point would imply a bi‐stability of the Arctic climate—where multiple steady‐state Arctic climates are possible at the same CO 2 value. Previous works have sought to establish the existence of bi‐stability using a range of models, from zero‐dimensional sea ice thermodynamic models to fully coupled global climate models, with conflicting results. Here, we present a new model of the Arctic that includes both sea‐ice thermodynamics and key atmospheric feedbacks in a simple framework. We exploit the model's simplicity to identify physical mechanisms that control the timing and extent of sea‐ice bi‐stability, and the abruptness of ice loss. We show that longwave radiation feedbacks can have a strong influence on Arctic surface climate from atmospheric temperature increases alone, even without major contributions from clear‐sky moisture or convective clouds suggested previously. While winter sea‐ice bi‐stability is robust to changes in uncertain model parameters in this study, summer sea ice is more sensitive. Finally, our model indicates that positive feedbacks may modulate the CO 2 threshold of sea‐ice loss and the width of bi‐stability much more strongly than the abruptness of loss. These results lead to a comprehensive understanding of the conditions that favor Arctic sea‐ice bi‐stability, particularly the role of atmospheric feedbacks, in both future and past climates.

Meteorology & Atmospheric Sciences↗
DOE OSTI2023

Atmospheric ice nucleation

Atmospheric ice nucleation is crucial for global precipitation and affects the structure, lifetime and reflectivity of clouds, thereby impacting climate. Ice nucleates in various ways from aerosol particles, termed ice-nucleating particles, over an extensive temperature and humidity range. Quantifying the kinetic and thermodynamic regimes of nucleation is necessary to relate fundamental physics to theoretically based predictions of ice formation for implementation in cloud and climate models. We review how the molecular picture of ice nucleation has advanced in recent years and consequential impacts on the interpretation and parameterization of ice nucleation. Advances include the role of interfacial free energy and pressure on ice nucleation rates, mobility regions of water that generate the critical ice nucleus, classical and non-classical pathways of nucleation, the type of ice polymorph that forms, the impact of solutes on freezing and the role of nanopores as surface features promoting ice nucleation. We also introduce currently debated and evaluated freezing parameterizations for application in model environments. Finally, we outline what we believe are the current needs for improving predictive understanding of ice nucleation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Measured and modeled acoustic propagation underneath the rough Arctic sea-ice

A characteristic surface duct beneath the sea-ice in the Marginal Ice Zone causes acoustic waves to be trapped and continuously interact with the sea-ice. The reflectivity of the sea-ice depends on the thickness, the elastic properties, and its roughness. This work focuses on the influence of sea-ice roughness on long-range acoustic propagation, and on how well the arrival structure can be predicted by the full wave integration model OASES. In 2013, acoustic signals centered at 900 Hz were transmitted every hour for three days between ice-tethered buoys in a drifting network in the Fram Strait. The experiment was set up to study the signal stability in the surface channel below the sea-ice. Oceanographic profiles were collected during the experiment, while a statistical description of the rough sea-ice was established based on historical ice-draft measurements. This environmental description is used as input to the range independent version of OASES. The model simulations correspond fairly well with the observations, despite that a flat bathymetry is used and the sea-ice roughness cannot be fully approximated by the statistical representation used in OASES. Long-range transmissions around 900 Hz are found to be more sensitive to the sea-ice roughness than the elastic parameters.

Hope, Gaute↗
DOE OSTI2013

Ice-Crystallization Kinetics in the Catalyst Layer of a Proton-Exchange-Membrane Fuel Cell

Nucleation and growth of ice in the catalyst layer of a proton-exchange-membrane fuel cell (PEMFC) are investigated using isothermal differential scanning calorimetry and isothermal galvanostatic cold-starts. Isothermal ice-crystallization rates and ice-nucleation rates are obtained from heat-flow and induction-time measurements at temperatures between 240 and 273 K for four commercial carbon-support materials with varying ionomer fraction and platinum loading. Measured induction times follow expected trends from classical nucleation theory and reveal that the carbon-support material and ionomer fraction strongly impact the onset of ice crystallization. Conversely, dispersed platinum particles play little role in ice crystallization. Following our previous approach, a nonlinear ice-crystallization rate expression is obtained from Johnson-Mehl-Avrami-Kolmogorov (JMAK) theory. A validated rate expression is now available for predicting ice crystallization within water-saturated catalyst layers. Using a simplified PEMFC isothermal cold-start continuum model, we compare cell-failure time predicted using the newly obtained rate expression to that predicted using a traditional thermodynamic-based approach. From this comparison, we identify conditions under which including ice-crystallization kinetics is critical and elucidate the impact of freezing kinetics on low-temperature PEMFC operation. The numerical model illustrates that cell-failure time increases with increasing temperature due to a longer required time for ice nucleation. Hence, ice-crystallization kinetics is critical when induction times are long (i.e., in the "nucleation-limited" regime for T > 263 K). Cell-failure times predicted using ice-freezing kinetics are in good agreement with the isothermal cold-starts, which also exhibit long and distributed cell-failure times for T > 263 K. These findings demonstrate a significant departure from cell-failure times predicted using the thermodynamic-based approach.

Dursch, T. J.↗
DOE OSTI2023

Vertical Structure and Ice Production Processes of Shallow Convective Postfrontal Clouds over the Southern Ocean in MARCUS. Part I: Observational Study

A study of the vertical structure of postfrontal shallow clouds in the marine boundary layer over the Southern Ocean is presented here. The central question of this two-part study regards cloud phase (liquid/ice) of precipitation, and the associated growth mechanisms. In this first part, data from the Measurements of Aerosols, Radiation, and Clouds over the Southern Ocean (MARCUS) field campaign are analyzed, starting with a 75-h case with continuous sea surface-based thermal instability, modest surface heat fluxes, an open-cellular mesoscale organization, and very few ice nucleating particles (INPs). The clouds are mostly precipitating and shallow (tops mostly around 2 km above sea level), with weak up- and downdrafts, and with cloud-top temperatures generally around –18° to –10°C. Further, the case study is extended to three other periods of postfrontal shallow clouds in MARCUS. While abundant supercooled liquid water is commonly present, an experimental cloud-phase algorithm classifies nearly two-thirds of clouds in the 0° to –5°C layer as containing ice (cloud ice, snow, or mixed phase), implying that much of the precipitation grows through cold-cloud processes. The best predictors of ice presence are cloud-top temperature, cloud depth, and INP concentration. Measures of convective activity and turbulence are found to be poor indicators of ice presence in the studied environment. The water-phase distribution in this cloud regime is explored through numerical simulations in Part II. Significance Statement: Climate models generally predict a lower albedo than observed over the Southern Ocean, and this is largely attributed to a lack of cloudiness, especially in the postfrontal cold sector of midlatitude cyclones. This in turn may be due to an excess of ice in these simulated clouds, resulting in rapid precipitation fallout and an overly brief cloud lifespan. The objective of this study is to examine whether shallow postfrontal clouds over the Southern Ocean are dominated by supercooled drops, or by snow and ice, using data collected by a U.S. Department of Energy Atmospheric Radiation Measurement Mobile Facility deployed aboard an Australian Antarctic supply vessel. We find that these clouds contain much supercooled liquid, even though cloud-top temperatures generally are around –18° to –8°C, and that about two-thirds of the clouds just above the freezing level contain ice. Much of the precipitation appears to grow through cold-cloud processes above the freezing level, rather than drizzle/rain. Updrafts and/or turbulence in convection or in cloud-top generating cells do not initiate much ice, compared to observations elsewhere in a similar temperature range. This may be attributable to the extremely low concentration of ice nucleating particles in this environment. Ultimately, the deepest clouds with the coldest cloud tops are most likely to be ice dominated.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Development of a triple-moment ice-phase cloud microphysics scheme and its application to the Single Column Atmosphere Model

Parameterization of cloud microphysics is critical for accurate simulation of weather and climate, in which the characteristics of cloud particle spectrum inevitably further affect climate simulation by changing the cloud evolution and cloud radiation effects. The popular currently used double-moment cloud microphysics schemes in numerical models can predict the intercept (N 0 ) and slope (λ) parameters of cloud particle spectrum but cannot predict the spectral shape parameter (μ), which hinders accurate description of cloud physical processes in climate models. Therefore, in the present study, we built upon the ideas of previously developed triple-moment cloud microphysics scheme, considered radar reflectivity factor as the third predictor in addition to number and mass concentration, and introduced its related prediction equation to the Single Column Atmosphere Model Version 5.3 (SCAM5.3). Moreover, the relevant microphysical process formulas in the model were revised, and a triple-moment ice-cloud microphysics scheme was constructed to predict the μ of ice particle spectrum. Based on this model, a 29-day case of the Atmospheric Radiation Measurement Program in the summer of 1997 (ARM97) was simulated, and the differences in cloud fraction and radiation simulations between the double- and triple-moment schemes were analysed. The μ of ice particle spectrum predicted by the triple-moment scheme mainly ranged from 0 to 4; the peak value was around 2, and at least 85% of the μ values were greater than the default 0. Therefore, the μ=0 setting in the double-moment scheme is unreasonable. Moreover, the developed triple-moment ice-cloud microphysics scheme yielded a narrower ice particle spectrum, which was closer to the observation results of previous studies, than the double-moment cloud microphysics scheme (e.g., height=233 hPa). Furthermore, compared with the double-moment scheme, the triple-moment scheme achieved closer simulations of the cloud fraction observations, particularly for ice clouds in the upper levels. Considering the close association between cloud fraction and radiation, the error with the triple-moment scheme was smaller than that with the double-moment scheme regardless of the shortwave (surface downward shortwave radiation, net downward shortwave radiation at the top of the atmosphere, and net shortwave radiation at the surface) or longwave (surface downward longwave radiation and net longwave radiation at the surface) flux density. Finally, the improvement mechanism of the triple-moment scheme on cloud fraction and radiation simulations was explored. The major reason is that the triple-moment scheme weakens the autoconversion of ice crystal to form snow and enhances the growth process of ice crystal deposition, thereby increasing the mass concentration of ice crystals and ice cloud fraction. Overall, the developed triple-moment ice-cloud microphysics scheme can help improve the simulation ability of models for cloud-climate feedback and other processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

A New Discrete Element Sea-Ice Model for Earth System Modeling

Sea ice forms a frozen crust of sea water oating in high-latitude oceans. It is a critical component of the Earth system because its formation helps to drive the global thermohaline circulation, and its seasonal waxing and waning in the high north and Southern Ocean signi cantly affects planetary albedo. Usually 4{6% of Earth's marine surface is covered by sea ice at any one time, which limits the exchange of heat, momentum, and mass between the atmosphere and ocean in the polar realms. Snow accumulates on sea ice and inhibits its vertical growth, increases its albedo, and contributes to pooled water in melt ponds that darken the Arctic ice surface in the spring. Ice extent and volume are subject to strong seasonal, inter-annual and hemispheric variations, and climatic trends, which Earth System Models (ESMs) are challenged to simulate accurately (Stroeve et al., 2012; Stocker et al., 2013). This is because there are strong coupled feedbacks across the atmosphere-ice-ocean boundary layers, including the ice-albedo feedback, whereby a reduced ice cover leads to increased upper ocean heating, further enhancing sea-ice melt and reducing incident solar radiation re ected back into the atmosphere (Perovich et al., 2008). A reduction in perennial Arctic sea-ice during the satellite era has been implicated in mid-latitude weather changes, including over North America (Overland et al., 2015). Meanwhile, most ESMs have been unable to simulate observed inter-annual variability and trends in Antarctic sea-ice extent during the same period (Gagne et al., 2014).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Higher Efficiency, Demand Flexible Refrigerator with On-Demand Micro-Vibrational De-icing Technology

Refrigerator technology has advanced significantly over the last couple of decades. Today’s refrigerators use only about 25% of the energy that was required to power models built in 1975. Even as they continually improve efficiency to meet standards, refrigerators have increased in size by almost 20%, added energy-consuming features such as through-the-door ice, and provide more benefits than ever before. However, a few challenges and technology gaps are preventing further improvement of the demand responsiveness and efficiency of the refrigerators. One of the major technology gaps in existing refrigerators is their outdated de-icing process. When the evaporator generates frost, an old-fashioned resistive heating element melts the ice. Most refrigerators have a timed defrost cycle, rather than an active system that could monitor the state of the frost. In these systems, not only is the precious electricity used at its least efficient form of conversion (direct conversion of electricity to heat), but also all the latent heat associated with the ice is wasted during the melting process. On top of that, the refrigerator needs to work harder to pull the temperature down after defrosting, and, last but not least, the food quality is severely impacted by the temperature swings during the defrost cycle. According to a study, the EU alone wastes 89 million tons of food in the supply chain every year. Any temperature swing during defrosting (about 6F according to Emerson for low-temperature cases) can negatively impact the shelf life of meat and other products for multiple days. All these issues can happen during the peak demand time of the electric grid. Unlike the conventional systems, the proposed novel advanced micro-vibrational deicing process uses no heat for defrosting. Instead, it uses the micro vibrations generated by a piezoelectric or vibration-generating module to mechanically break ice from the heat exchanger almost instantaneously. The project titled “Higher Efficiency, Demand Flexible Refrigerator with On-Demand Micro-Vibrational De-icing Technology, performed by Ultrasonic Technology Solutions, LLC (UTS) of Knoxville, TN, in collaboration with Emerson (now Copeland), represents the final phase of a multi-year effort funded under the U.S. Department of Energy’s Building Technologies Office (BTO) BENEFIT FOA 2020. Initiated on October 1, 2021, and completed after a nine-month no-cost extension ending September 30, 2025, this project aimed to develop and validate a novel micro-vibrational mechanical defrosting system, achieving more than 25% improvement in defrosting energy efficiency over conventional baseline defrosting technologies. Over sixteen quarters, the project advanced from fundamental ice-mechanical characterization and prototype development to full-scale system integration and validation. Initial efforts established project management infrastructure and characterized ice adhesion properties, followed by the design and fabrication of early aluminum-based prototypes for resonance frequency testing. Subsequent quarters saw rapid technical progression, including the identification of optimal piezoelectric and motor-based vibration mechanisms, the demonstration of effective de-icing over 6x6-inch aluminum surfaces. The team achieved its Go/No-Go milestone by exceeding the 25% energy-efficiency improvement target—reaching up to 3,340% under optimized conditions—and later confirmed that motor-driven systems offered superior performance and energy efficiency compared to piezoelectric alternatives. Continued refinement led to the development of amplifier systems on printed circuit boards, improved control and instrumentation hardware, and integration into full-scale heat exchanger (HX) prototypes at both UTS and Copeland facilities. Multiple vibration-mounting studies and frost-growth experiments guided mechanical optimization and noise-mitigation strategies, achieving a 17.5 dB reduction in sound pressure level and verifying robust mechanical performance. Advanced analyses, including modal and harmonic simulations, established a quantitative understanding of vibrational behavior and de-icing efficiency across >1000 cm² systems. The final project phase successfully demonstrated scalable integration within reach-in and chest freezer prototypes, confirmed >25% efficiency improvements in large-area systems, and completed a comprehensive business model and scale-up strategy identifying electric defrost systems as the primary beachhead market. The culmination of this DOE-supported effort establishes micro-vibrational defrosting as a viable, high-efficiency, low-noise, and demand-flexible de-icing technology, paving the way for commercial deployment and broader application in next-generation refrigeration systems.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Re-evaluating cloud chamber constraints on depositional ice growth in cirrus clouds – Part 1: Model description and sensitivity tests

Abstract. Ice growth from vapor deposition is an important process for the evolution of cirrus clouds, but the physics of depositional ice growth at the low temperatures (<235 K) characteristic of the upper troposphere and lower stratosphere is not well understood. Surface attachment kinetics, generally parameterized as a deposition coefficient αD, control ice crystal habit and also may limit growth rates in certain cases, but significant discrepancies between experimental measurements have not been satisfactorily explained. Experiments on single ice crystals have previously indicated the deposition coefficient is a function of temperature and supersaturation, consistent with growth mechanisms controlled by the crystal's surface characteristics. Here we use observations from cloud chamber experiments in the Aerosol Interactions and Dynamics in the Atmosphere (AIDA) aerosol and cloud chamber to evaluate surface kinetic models in realistic cirrus conditions. These experiments have rapidly changing temperature, pressure, and ice supersaturation such that depositional ice growth may evolve from diffusion limited to surface kinetics limited over the course of a single experiment. In Part 1, we describe the adaptation of a Lagrangian parcel model with the Diffusion Surface Kinetics Ice Crystal Evolution (DiSKICE) model (Zhang and Harrington, 2014) to the AIDA chamber experiments. We compare the observed ice water content and saturation ratios to that derived under varying assumptions for ice surface growth mechanisms for experiments simulating ice clouds between 180 and 235 K and pressures between 150 and 300 hPa. We found that both heterogeneous and homogeneous nucleation experiments at higher temperatures (>205 K) could generally be modeled consistently with either a constant deposition coefficient or the DiSKICE model assuming growth on isometric crystals via abundant surface dislocations. Lower-temperature experiments showed more significant deviations from any depositional growth model, with different ice growth rates for heterogeneous and homogeneous nucleation experiments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Geometric remapping of particle distributions in the Discrete Element Model for Sea Ice (DEMSI v0.0)

Abstract. A new sea ice dynamical core, the Discrete Element Model for Sea Ice (DEMSI), is under development for use in coupled Earth system models. DEMSI is based on the discrete element method, which models collections of ice floes as interacting Lagrangian particles. In basin-scale sea ice simulations the Lagrangian motion results in significant convergence and ridging, which requires periodic remapping of sea ice variables from a deformed particle configuration back to an undeformed initial distribution. At the resolution required for Earth system models we cannot resolve individual sea ice floes, so we adopt the sub-grid-scale thickness distribution used in continuum sea ice models. This choice leads to a series of hierarchical tracers depending on ice fractional area or concentration that must be remapped consistently. The circular discrete elements employed in DEMSI help improve the computational efficiency at the cost of increased complexity in the effective element area definitions for sea ice cover that are required for the accurate enforcement of conservation. An additional challenge is the accurate remapping of element values along the ice edge, the location of which varies due to the Lagrangian motion of the particles. In this paper we describe a particle-to-particle remapping approach based on well-established geometric remapping ideas that enforces conservation, bounds preservation, and compatibility between associated tracer quantities, while also robustly managing remapping at the ice edge. One element of the remapping algorithm is a novel optimization-based flux correction that enforces concentration bounds in the case of nonuniform motion. We demonstrate the accuracy and utility of the algorithm in a series of numerical test cases.

58 GEOSCIENCES↗
DOE OSTI2018

Design and results of the ice sheet model initialisation experiments initMIP-Greenland: an ISMIP6 intercomparison

Abstract. Earlier large-scale Greenland ice sheet sea-level projections (e.g. those run during the ice2sea and SeaRISE initiatives) have shown that ice sheet initial conditions have a large effect on the projections and give rise to important uncertainties. The goal of this initMIP-Greenland intercomparison exercise is to compare, evaluate, and improve the initialisation techniques used in the ice sheet modelling community and to estimate the associated uncertainties in modelled mass changes. initMIP-Greenland is the first in a series of ice sheet model intercomparison activities within ISMIP6 (the Ice Sheet Model Intercomparison Project for CMIP6), which is the primary activity within the Coupled Model Intercomparison Project Phase 6 (CMIP6) focusing on the ice sheets. Two experiments for the large-scale Greenland ice sheet have been designed to allow intercomparison between participating models of (1) the initial present-day state of the ice sheet and (2) the response in two idealised forward experiments. The forward experiments serve to evaluate the initialisation in terms of model drift (forward run without additional forcing) and in response to a large perturbation (prescribed surface mass balance anomaly); they should not be interpreted as sea-level projections. We present and discuss results that highlight the diversity of data sets, boundary conditions, and initialisation techniques used in the community to generate initial states of the Greenland ice sheet. We find good agreement across the ensemble for the dynamic response to surface mass balance changes in areas where the simulated ice sheets overlap but differences arising from the initial size of the ice sheet. The model drift in the control experiment is reduced for models that participated in earlier intercomparison exercises.

54 ENVIRONMENTAL SCIENCES↗