Search NASA⌕ Search

SEARCH · Search NASA

Results for “ice models”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Evaluation of a Thermodynamic Ice Crystal Icing Model Using Experimental Ice Accretion Data

This paper presents the evaluation of a thermodynamic ice crystal icing model, previously presented to describe the possible mechanisms of icing within the core of a turbofan jet engine. It has been proposed that there are two types of distinct ice accretions based on a surface energy balance: freeze-dominated icing and melt-dominated icing. In the former, ice accretion occurs where a freeze fraction (0 to 1) of melted ice crystals freezes on a surface, along with the existing ice of the impinging water and ice mass. This freeze-dominated icing is characterized by having strong adhesion to the surface. In the latter, icing occurs from accumulated unmelted ice on a surface, where a melt fraction (0 to 1) dictates the amount of unmelted impinged ice. This melt-dominated icing is characterized by weakly bonded surface adhesion. The 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. This mass loss parameter must be determined in conjunction with the fraction of freezing liquid water or fraction of melting ice on an icing surface. This loss parameter, however, along with the freeze and melt fraction, are the only experimental parameters that are currently not measured directly. Using reported icing growth rates from published ice crystal icing experiments, a methodology is proposed to determine these unknown parameters. This work takes reported ice accretion data from tests conducted by the National Aeronautics and Space Administration (NASA) in 2016 and tests NASA collaborated on with the National Research Council (NRC) of Canada in 2012 that examined the fundamental physics of ice crystal icing. Those research efforts sought to generate icing conditions representative of those that occur inside a jet engine when ingesting ice crystals. This paper presents the fundamental equations of the thermodynamic model, the methodology used to determine the aforementioned unknown icing parameters, and results from model evaluation using experimental data. In addition, this paper builds on the previously proposed model by adding a transient conduction term to explain ice growth behavior at the onset of experimental tests that was observed to be different from steady-state ice growth that occurred later in the test run.With the addition of this energy term, this becomes a quasi-steady model. A key finding from this work suggests that mass loss fractions can exceed 0.90 for steady ice growth periods. In addition, due to conductive heat fluxes when using a warmer-than-freezing airfoil, lower mass loss fraction values were calculated during the initial transient period.

Ice-Crystal Ice Accretion↗
NASA NTRS1984

The sensitivity of a one-dimensional thermodynamic sea ice model to changes in cloudiness

A thermodynamic sea ice-lead model is used to assess the importance of cloud cover changes to modeled ice thickness. For regions of either permanent multiyear ice or seasonal sea ice, the cloud amount variations have relatively little impact. However, for regions where the presence of summer ice is variable from year to year, the predicted ice thickness is strongly dependent on cloud cover. In general, with a snow covered surface, decreased cloud leads to surface cooling while increased cloud gives rise to a surface warming. For a melting bare ice surface, the reverse occurs. The ice model response time is too long for interannual variations in cloud amount to explain interannual variations in ice thickness and extent. Nevertheless, the implication of the results is that numerical modeling of sea ice distribution requires accurate cloud data or cloud prediction and that trends in cloud cover may lead to significant perturbations in sea ice extent and thickness.

Shine, K. P.↗
NASA NTRS2012

The Annual Glaciohydrology Cycle in the Ablation Zone of the Greenland Ice Sheet: Part 2. Observed and Modeled Ice Flow

Ice velocities observed in 2005/06 at three GPS stations along the Sermeq Avannarleq flowline, West Greenland, are used to characterize an observed annual velocity cycle. We attempt to reproduce this annual ice velocity cycle using a 1-D ice-flow model with longitudinal stresses coupled to a 1-D hydrology model that governs an empirical basal sliding rule. Seasonal basal sliding velocity is parameterized as a perturbation of prescribed winter sliding velocity that is proportional to the rate of change of glacier water storage. The coupled model reproduces the broad features of the annual basal sliding cycle observed along this flowline, namely a summer speed-up event followed by a fall slowdown event. We also evaluate the hypothesis that the observed annual velocity cycle is due to the annual calving cycle at the terminus. We demonstrate that the ice acceleration due to a catastrophic calving event takes an order of magnitude longer to reach CU/ETH ('Swiss') Camp (46km upstream of the terminus) than is observed. The seasonal acceleration observed at Swiss Camp is therefore unlikely to be the result of velocity perturbations propagated upstream via longitudinal coupling. Instead we interpret this velocity cycle to reflect the local history of glacier water balance.

ice flow↗
NASA NTRS2015

Recent Advances in the LEWICE Icing Model

This paper will describe two recent modifications to the Glenn ICE software. First, a capability for modeling ice crystals and mixed phase icing has been modified based on recent experimental data. Modifications have been made to the ice particle bouncing and erosion model. This capability has been added as part of a larger effort to model ice crystal ingestion in aircraft engines. Comparisons have been made to ice crystal ice accretions performed in the NRC Research Altitude Test Facility (RATFac). Second, modifications were made to the run back model based on data and observations from thermal scaling tests performed in the NRC Altitude Icing Tunnel.

fluid dynamics↗
NASA NTRS2018

The NASA Eulerian Snow on Sea Ice Model (NESOSIM) v1.0: initial model development and analysis

The NASA Eulerian Snow On Sea Ice Model (NESOSIM) is a new, open-source snow budget model that is currentlyconfigured to produce daily estimates of the depth and density of snow on sea ice across the Arctic Ocean throughthe accumulation season. NESOSIM has been developed in a three-dimensional Eulerian framework and includes two(vertical) snow layers and several simple parameterizations (accumulation, wind packing, advection_divergence, blowingsnow lost to leads) to represent key sources and sinks of snow on sea ice. The model is forced with daily inputs of snowfall and near-surface winds (from reanalyses), sea ice concentration (from satellite passive microwave data) and sea ice drift (from satellite feature tracking) during the accumulation season (August through April). In this study, we present the NESOSIM formulation, calibration efforts, sensitivity studies and validation efforts across an Arctic Ocean domain (100 km horizontal resolution). The simulated snow depth and density are calibrated with in situ data collected on drifting ice stations during the 1980s. NESOSIM shows strong agreement with the in situ seasonal cycles of snow depth and density, and shows good (moderate) agreement with the regional snow depth (density) distributions. NESOSIM is run for a contemporary period (2000 to 2015), with the results showing strong sensitivity to the reanalysis derived snowfall forcing data, with the Modern-Era Retrospective analysis for Research and Applications (MERRA) and the Japanese Meteorological Agency 55-year reanalysis (JRA-55) forced snow depths generally higher than ERA-Interim,and the Arctic System Reanalysis (ASR) generally lower. We also generate and force NESOSIM with a consensus "median" daily snowfall dataset from these reanalyses.The results are compared against snow depth estimates derived from NASA's Operation IceBridge (OIB) snow radar data from 2009 to 2015, showing moderate_strong correlations and root mean squared errors of 10 cm depending on the OIB snow depth product analyzed, similar to the comparisons between OIB snow depths and the commonly used modified Warren snow depth climatology. Potential improvements to this initial NESOSIM formulation are discussed in the hopes of improving the accuracy and reliability of these simulated snow depths and densities.

Petty, Alek A.↗
NASA NTRS2007

A New Sea Ice Model: Pancakes to Metrics

We are developing, testing, and validating a new sea ice dynamics model that treats the ice cover as an elastic/decohesive material in the permanent pack and includes the correct frazil/pancake behavior in the marginal zone. Two salient features of present ice dynamics models are that they do not: 1) reproduce the oriented fracture patterns of openings and closings in the pack ice, and 2) accurately model the effects of frazil/pancake ice formation in the ice margin. These poorly modeled areas account for a substantial portion of the ice growth, turbulent heat flux to the atmosphere, salt flux to the ocean, and energy dissipation due to slippage, ridging, and rafting, in the Arctic. Existing sea ice models have shown limited success in predicting the degree to which a lead will open for prescribed or observed forcing conditions. An important aspect of the new model we are developing is that the existence of cracks, along with their orientation, opening, and closing, is predicted. To put this effort in perspective a short history of ice dynamics modeling and data collection is presented. The RGPS data set is used to validate the model. As part of the testing and validation of the model, we are working on a new metric for comparing linear features (leads and ridges) in the data and model to be used in data assimilation for this model. The model framework is presented as well as some results showing the creation and development of leads in a simulation of ice dynamics in the Beaufort Sea. Other presentations by the authors will show other results from this effort.

Coon, M.↗
NASA NTRS1993

Modeling ice streams: Derived quantities

The model addressed is a finite-element, map-plane, time-dependent, column-averaged continuity equation solver. The key to the fitting process involves the balance between ice motion dominated by flow in internal layers, and ice motion dominated by sliding at the bed. The fitting process involves an iterative process carried out in the time domain. Beginning with the portion of the ice sheet being modeled identical to the present ice sheet with uniform flow, sliding, and fraction specified at nominal values, the model monitors each nodal point surface elevation. As the calculated surface elevation deviates from the present surface, a correction proportional to the difference is applied to selected parameter sets. This correction is in a sense that would tend to improve the fit at the particular nodal point. A calculated surface elevation that was higher than the present surface would result in an increased fraction, which would tend to lower the calculated surface (if the flow or sliding constant were being used as the fitting parameter, they would be lowered to improve the fit). This process is allowed to proceed as long as is necessary for the situation to stabilize. Typically, this takes tens of thousands of model years, but the rate is dependent on other external forcings such as the accumulation rate. The primary result is that while a typical sample of ice streams from around Antarctica can be fitted quite reasonably using only the fraction of the velocity due to sliding, a different mechanism seems to be in play along the Siple Coast, where reduced sliding constants are required to attain a reasonable fit. Flow is more strongly channelized in this region, and velocities are, in general, higher than are observed in other regions. It is unlikely that the mechanism that controls the ice movement along the Siple Coast is exactly similar to the mechanisms in the other ice streams. The concept of deformable sediments and their contribution to the fast flow along the Siple Coast may have limited applicability to other Antarctic ice streams.

Fastook, James↗
NASA NTRS1985

Sensitivity of a thermodynamic sea ice model with leads to time step size

The characteristics of sea ice models, developed to study the physics of the growth and melt of ice at the ocean surface and the variations in ice extent, depend on the size of the time step. Thus, to study longer-term variations within a reasonable computer budget, a model with a scheme allowing longer time steps has been constructed. However, the results produced by the model can definitely depend on the length of the time step. The sensitivity of a model to time-step size can be reduced by appropriate approaches. The present investigation is concerned with experiments which use a formulation of a lead parameterization that can be considered as a first step toward the development of a lead parameterization suitable for a use in long-term climate studies.

Ledley, T. S.↗
NASA NTRS1982

Sensitivity of a climatologically-driven sea ice model to the ocean heat flux

Ocean heat flux sensitivity was studied on a numerical model of sea ice covering the Weddell Sea region of the southern ocean. The model is driven by mean monthly climatological atmospheric variables. For each model run, the ocean heat flux is uniform in both space and time. Ocean heat fluxes below 20 W m to the minus 2 power do not provide sufficient energy to allow the ice to melt to its summertime thicknesses and concentrations by the end of the 14 month simulation, whereas ocean heat fluxes of 30 W m to the minus 2 power and above result in too much ice melt, producing the almost total disappearance of ice in the Weddell Sea by the end of the 14 months. These results are dependent on the atmospheric forcing fields.

Parkinson, C. L.↗
NASA NTRS2005

Current Experimental Basis for Modeling Ice Accretions on Swept Wings

This work presents a review of the experimental basis for modeling ice accretions on swept wings. Experimental work related to ice accretion physics on swept wings conducted between 1954 and 2004 is reviewed. Proposed models or explanations of scallop formations are singled out and discussed. Special emphasis is placed on reviewing the work done to determine the basic macroscopic mechanisms of scallop formation. The role of feather growth and its connection to scallop growth is discussed. Conceptual steps in modeling scallop formations are presented. Research elements needed for modeling are discussed.

Vargas, Mario↗
DOE OSTI2024

Refactoring the elastic–viscous–plastic solver from the sea ice model CICE v6.5.1 for improved performance

This study focuses on the performance of the elastic–viscous–plastic (EVP) dynamical solver within the sea ice model, CICE v6.5.1. The study has been conducted in two steps. First, the standard EVP solver was extracted from CICE for experiments with refactored versions, which are used for performance testing. Second, one refactored version was integrated and tested in the full CICE model to demonstrate that the new algorithms do not significantly impact the physical results. The study reveals two dominant bottlenecks, namely (1) the number of Message Parsing Interface (MPI) and Open Multi-Processing (OpenMP) synchronization points required for halo exchanges during each time step combined with the irregular domain of active sea ice points and (2) the lack of single-instruction, multiple-data (SIMD) code generation. The standard EVP solver has been refactored based on two generic patterns. The first pattern exposes how general finite differences on masked multi-dimensional arrays can be expressed in order to produce significantly better code generation by changing the memory access pattern from random access to direct access. The second pattern takes an alternative approach to handle static grid properties. The measured single-core performance improvement is more than a factor of 5 compared to the standard implementation. The refactored implementation of strong scales on the Intel® Xeon® Scalable Processors series node until the available bandwidth of the node is used. For the Intel® Xeon® CPU Max series, there is sufficient bandwidth to allow the strong scaling to continue for all the cores on the node, resulting in a single-node improvement factor of 35 over the standard implementation. This study also demonstrates improved performance on GPU processors.

58 GEOSCIENCES↗
DOE OSTI2024

Sensitivities of the West Greenland Current to Greenland Ice Sheet Meltwater in a Mesoscale Ocean/Sea Ice Model

Abstract Meltwater from the Greenland Ice Sheet can alter the continental shelf/slope circulation and cross-shelf freshwater fluxes and limit deep convection in adjacent basins through surface freshening. We explore the impacts on the West Greenland Current and eastern Labrador Sea with different vertical distributions of the meltwater forcing. In this study, we present the results from global coupled ocean/sea ice simulations, forced with atmospheric reanalysis, that are mesoscale eddy-active (∼2–3-km horizontal spacing) and eddy-permitting (∼6–7-km horizontal spacing) in the study region. We compare the West Greenland Current in mesoscale eddy-active and eddy-permitting without meltwater to highlight the role of small-scale features. The mesoscale eddy-active configuration is then used to assess the change in the eastern Labrador Sea when meltwater is added to the surface or vertically distributed to account for mixing within fjords. In both simulations with meltwater, the West Greenland and West Greenland Coastal Currents are faster than in the simulation with no meltwater; their mean surface speeds are the highest in the vertical distribution case. In the latter case, there is enhanced baroclinic conversion at the shelf break compared to the simulation with no meltwater. When meltwater is vertically distributed, there is an increase in baroclinic conversion at the shelf break associated with increased eddy kinetic energy. In addition, in the eastern Labrador Sea, the salinity is lower and the meltwater volume is greater when meltwater is vertically distributed. Therefore, the West Greenland Current is sensitive to how meltwater is added to the ocean with implications for the freshening of the Labrador Sea. Significance Statement Our goal is to understand how the flux of freshwater across the West Greenland continental slope into the Labrador Sea is modified by meltwater from the Greenland Ice Sheet. We compare the simulations of the ocean that capture key dynamics along the West Greenland continental slope that have no meltwater, meltwater added to the ocean surface, and meltwater distributed vertically to represent the mixing within fjords. When meltwater is added, the currents along the continental slope are faster, with the greatest increase when meltwater is vertically distributed. In that case, there is enhanced freshening of the Labrador Sea because modified density gradients generate more eddies. Proper representation of the vertical structure of meltwater is important for projecting the impact of freshwater on the subpolar North Atlantic.

58 GEOSCIENCES↗
NASA NTRS1985

A coupled ice-ocean model of ice breakup and banding in the marginal ice zone

A coupled ice-ocean numerical model for the marginal ice zone is considered. The model consists of a nonlinear sea ice model and a two-layer (reduced gravity) ocean model. The dependence of the upwelling response on wind stress direction is discussed. The results confirm earlier analytical work. It is shown that there exist directions for which there is no upwelling, while other directions give maximum upwelling in terms of the volume of uplifted water. The ice and ocean is coupled directly through the stress at the ice-ocean interface. An interesting consequence of the coupling is found in cases when the ice edge is almost stationary. In these cases the ice tends to break up a few tenths of kilometers inside of the ice edge.

Smedstad, O. M.↗
DOE OSTI2025

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↗
NASA NTRS1997

Scale Model Icing Research Tunnel

NASA Lewis Research Center's Icing Research Tunnel (IRT) is the world's largest refrigerated wind tunnel and one of only three icing wind tunnel facilities in the United States. The IRT was constructed in the 1940's and has been operated continually since it was built. In this facility, natural icing conditions are duplicated to test the effects of inflight icing on actual aircraft components as well as on models of airplanes and helicopters. IRT tests have been used successfully to reduce flight test hours for the certification of ice-detection instrumentation and ice protection systems. To ensure that the IRT will remain the world's premier icing facility well into the next century, Lewis is making some renovations and is planning others. These improvements include modernizing the control room, replacing the fan blades with new ones to increase the test section maximum velocity to 430 mph, installing new spray bars to increase the size and uniformity of the artificial icing cloud, and replacing the facility heat exchanger. Most of the improvements will have a first-order effect on the IRT's airflow quality. To help us understand these effects and evaluate potential improvements to the flow characteristics of the IRT, we built a modular 1/10th-scale aerodynamic model of the facility. This closed-loop scale-model pilot tunnel was fabricated onsite in the various shops of Lewis' Fabrication Support Division. The tunnel's rectangular sections are composed of acrylic walls supported by an aluminum angle framework. Its turning vanes are made of tubing machined to the contour of the IRT turning vanes. The fan leg of the tunnel, which transitions from rectangular to circular and back to rectangular cross sections, is fabricated of fiberglass sections. The contraction section of the tunnel is constructed from sheet aluminum. A 12-bladed aluminum fan is coupled to a turbine powered by high-pressure air capable of driving the maximum test section velocity to 550 ft/sec (Mach 0.45). The air turbine and instrumentation are housed inside a fiberglass nacelle. Total and static pressure measurements can be taken around the loop, and velocity and flow angularity measurements can be taken with hot-wire and five-hole probes at specific locations. The Scale Model Icing Research Tunnel (SMIRT) is undergoing checkout tests to determine how its airflow characteristics compare with the IRT. Near-term uses for this scale-model tunnel include determining the aerodynamic effects of replacing the 52-yearold W-shaped heat exchanger with a flat-faced heat exchanger. SMIRT is an integral part of the improvements planned for the IRT because testing the proposed IRT improvements in a scale-model tunnel will lower costs and improve productivity.

Canacci, Victor A.↗
NASA NTRS2016

Convection from Hemispherical and Conical Model Ice Roughness Elements in Stagnation Region Flows

To improve ice accretion prediction codes, more data regarding ice roughness and its effects on convective heat transfer are required. The Vertical Icing Studies Tunnel (VIST) at NASA Glenn Research was used to model realistic ice roughness in the stagnation region of a NACA 0012 airfoil. In the VIST, a test plate representing the leading 2% chord of the airfoil was subjected to flows of 7.62 m/s (25 ft/s), 12.19 m/s (40 ft/s), and 16.76 m/s (55 ft/s). The test plate was fitted with multiple surfaces or sets of roughness panels, each with a different representation of ice roughness. The sets of roughness panels were constructed using two element distribution patterns that were created based on a laser scan of an iced airfoil acquired in the Icing Research Tunnel at NASA Glenn. For both roughness patterns, surfaces were constructed using plastic hemispherical elements, plastic conical elements, and aluminum conical elements. Infrared surface thermometry data from tests run in the VIST were used to calculate area averaged heat transfer coefficient values. The values from the roughness surfaces were compared to the smooth control surface, showing convective enhancement as high as 400% in some cases. The data gathered during this study will ultimately be used to improve the physical modeling in LEWICE or other ice accretion codes and produce predictions of in-flight ice accretion on aircraft surfaces with greater confidence.

Icing Physics↗
NASA NTRS1998

Process Developed for Forming Urethane Ice Models

A new process for forming ice shapes on an aircraft wing was developed at the NASA Lewis Research Center. The innovative concept was formed by Lewis' Icing Research Tunnel (IRT) team, and the hardware was manufactured by Lewis' Manufacturing Engineering Division. This work was completed to increase our understanding of the stability and control of aircraft during icing conditions. This project will also enhance our evaluation of true aerodynamic wind tunnel effects on aircraft. In addition, it can be used as a design tool for evaluating ice protection systems.

Vannuyen, Thomas↗