SEARCH · Search NASA

Results for “ice flow model”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS2017

Comparisons of Mixed-Phase Icing Cloud Simulations with Experiments Conducted at the NASA Propulsion Systems Laboratory

This paper builds on previous work that compares numerical simulations of mixed-phase icing clouds with experimental data. The model couples the thermal interaction between ice particles and water droplets of the icing cloud with the flowing air of an icing wind tunnel for simulation of NASA Glenn Research Centers (GRC) Propulsion Systems Laboratory (PSL). Measurements were taken during the Fundamentals of Ice Crystal Icing Physics Tests at the PSL tunnel in March 2016. The tests simulated ice-crystal and mixed-phase icing that relate to ice accretions within turbofan engines.

Bartkus, Tadas↗

DOE OSTI2026

Xanthos-Lake Model Source Code

This repository contains the source code for Xanthos-Lake, a lake-modeling extension of the Xanthos framework that introduces a coupled lake component comprising the Xanthos-Lake Snow and Ice Model (xLSIM) and the Xanthos-Lake Water Balance Model (xLWBM). xLSIM is a basin-aware machine-learning model for lake snow, ice, and thermal conditions. It predicts monthly lake ice thickness, snow depth, snow-cover fraction, mixing-layer temperature, and lake ice fraction from meteorological forcing and lake surface-area information. It uses sequence-based deep-learning architectures, including Transformer and hybrid Long Short-Term Memory–Transformer (LSTM–Transformer) models, together with seasonal encoding, multi-lake learning, physical masking, and basin-level cryospheric and non-cryospheric classification. The training workflow uses Ray for scalable execution and includes optional Ray Tune hyperparameter optimization. Model predictions, observations, diagnostics, and feature-importance outputs are written in NetCDF. xLWBM is the water-balance component of the new lake framework. It simulates monthly lake storage, surface area, evaporation, inflow, outflow, and lake–groundwater exchange. It combines physical water-balance equations with calibrated bathymetric relationships, weir-based outlet flow, modified Penman open-water evaporation, groundwater head relaxation, Penman–Monteith snow and ice sublimation, and snow, ice, and thermal conditions supplied by xLSIM. The model calibrates lake parameters against satellite-derived surface-area data, using evaporation-based calibration where surface-area data are unavailable, and supports small, medium, and large lake classes. For large lakes, xLWBM is integrated with the managed-routing workflow so that lake storage and outflow interact directly with downstream river routing and reservoir operations. Together, xLSIM and xLWBM provide Xanthos with a coupled lake-modeling capability. xLSIM supplies the snow, ice, and thermal conditions that affect lake evaporation and snow- and ice-related water exchanges, while xLWBM translates those conditions into dynamic lake storage, surface area, evaporation, and discharge. In return, xLWBM supplies evolving lake surface area to xLSIM. This coupling enables Xanthos to represent lakes as active hydrologic components within basin-scale water-availability and routing simulations.

Machine Learning↗

NASA NTRS2004

Lunar and Planetary Science XXXV: Mars Polar Science and Exploration

The session, "Mars Polar Science and Exploration" included the following reports: Evidence for Possible Exposed Water Ice Deposits in Martian Low Latitude Chasms and Chaos; Stability and Exchange of Subsurface Ice on Mars; Evaporation Rates for Liquid Water and Ice Under Current Martian Conditions; Seasonal Cycle of Carbon Dioxide and Atmospheric Circulation in Mars Southern Hemisphere as Observed by Neutron Spectroscopy; Imaging Polarimetry of Mars with Hubble Space Telescope in 2003 Opposition; GCM Simulations of Tropical Ice Accumulations: Implications for Cold-based Glaciers; Numerical Modeling of Glaciers in Martian Paleoclimates; Valley Glaciers on Mars: Calculation of Flow Rate and Thickness; Internal Structure of the Southern Polar Cap of Mars and Formation Implications; Sublimation at the Base of a Seasonal CO2 Slab on Mars; Impact Crater Abundance of the Martian South Polar Layered Deposits from THEMIS Visible Imaging; Recent Changes in South-Polar-Polygonal Terrain During One Martian Year: Implications for Subsurface Ice-Wedges; and An Historical Search for Unfrozen Water at the Phoenix Landing Site.

Source record↗

NASA NTRS2002

Heat Flow, Thermal Conductivity, and the Plausibility of the White Mars Hypothesis

Due to the low thermal conductivity of CO2 ice and clathrate vs. water ice, we find that liquid water reservoirs would not be confined to the deep subsurface as predicted by the controversial White Mars model, even assuming low global heat flow. Additional information is contained in the original extended abstract.

Urquhart, M. L.↗

NASA NTRS1980

Structure and evolution of Uranus and Neptune

Three-layer interior models of Uranus and Neptune with central rocky cores, mantles of water, methane, and ammonia (the 'ices'), and outer envelopes primarily composed of hydrogen and helium are presented. The models incorporate a new H2O equation of state based on experimental data which is considerably 'softer' than previous H2O equations of state. Corrections for interior temperatures approximately 5000 K are included in the models, and the thermal evolution of both planets is investigated using recent heat flow measurements. It is found that the evolutionary considerations are consistent with gravitational field data in supporting models with approximately solar abundances of 'ice' and 'rock'. Evolutionary considerations indicate that initial temperatures and luminosities for Uranus and Neptune were not substantially higher than the present value. Both planets apparently have relatively small approximately 1-2 earth masses) hydrogen-helium envelopes, with Neptune's envelope smaller than Uranus'. A monotonic trend is evident among the Jovian planets: all have central rock-ice cores of approximately 15 earth masses, but with hydrogen-helium envelopes which decrease in mass from Jupiter to Saturn to Uranus to Neptune.

Hubbard, W. B.↗

NASA NTRS1994

Ocean-ice interaction in the marginal ice zone using synthetic aperture radar imagery

Ocean-ice interaction processes in the marginal ice zone (MIZ) by wind, waves, and mesoscale features, such as up/downwelling and eddies are studied using Earth Remote-Sensing Satellite (ERS) 1 synthetic aperture radar (SAR) images and an ocean-ice interaction model. A sequence of seven SAR images of the MIZ in the Chukchi Sea with 3 or 6 days interval are investigated for ice edge advance/retreat. Simultaneous current measurements from the northeast Chukchi Sea, as well as the Barrow wind record, are used to interpret the MIZ dynamics. SAR spectra of waves in ice and ocean waves in the Bering and Chukchi Sea are compared for the study of wave propagation and dominant SAR imaging mechanism. By using the SAR-observed ice edge configuration and wind and wave field in the Chukchi Sea as inputs, a numerical simulation has been performed with the ocean-ice interaction model. After 3 days of wind and wave forcing the resulting ice edge configuration, eddy formation, and flow velocity field are shown to be consistent with SAR observations.

Liu, Antony K.↗

NASA NTRS2015

Convective Enhancement of Icing Roughness Elements in Stagnation Region Flows

To improve existing ice accretion simulation codes, more data regarding ice roughness and its effects on convective heat transfer are required. To build on existing research on this topic, this study used the Vertical Icing Studies Tunnel (VIST) at NASA Glenn Research to model realistic ice roughness in the stagnation region of a NACA 0012 airfoil. Using 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 3 surfaces, each with a different representation of ice roughness: 1) a control surface with no ice roughness, 2) a surface with ice roughness with element height scaled by 10x and streamwise rough zone width from the stagnation point scaled by 10x, and 3) a surface with ice roughness with element height scaled by 10x and streamwise rough zone width from the stagnation point scaled by 25x. Temperature data from the tests were recorded using an infrared camera and thermocouples imbedded in the test plate. From the temperature data, a convective heat transfer coefficient map was created for each case. Additional testing was also performed to validate the VIST's flow quality. These tests included five-hole probe and hot-wire probe velocity traces to provide flow visualization and to study boundary layer formation on the various test surfaces. The knowledge gained during the experiments will help improve ice accretion codes by providing heat transfer coefficient validation data and by providing flow visualization data helping understand current and future experiments performed in the VIST.

Icing↗

NASA NTRS2010

Current Methods Modeling and Simulating Icing Effects on Aircraft Performance, Stability, Control

Icing alters the shape and surface characteristics of aircraft components, which results in altered aerodynamic forces and moments caused by air flow over those iced components. The typical effects of icing are increased drag, reduced stall angle of attack, and reduced maximum lift. In addition to the performance changes, icing can also affect control surface effectiveness, hinge moments, and damping. These effects result in altered aircraft stability and control and flying qualities. Over the past 80 years, methods have been developed to understand how icing affects performance, stability, and control. Emphasis has been on wind-tunnel testing of two-dimensional subscale airfoils with various ice shapes to understand their effect on the flowfield and ultimately the aerodynamics. This research has led to wind-tunnel testing of subscale complete aircraft models to identify the integrated effects of icing on the aircraft system in terms of performance, stability, and control. Data sets of this nature enable pilot-in-the-loop simulations to be performed for pilot training or engineering evaluation of system failure impacts or control system design.

Ratvasky, Thomas P.↗

NASA NTRS2013

Aerodynamic Characterization of a Thin, High-Performance Airfoil for Use in Ground Fluids Testing

The FAA has worked with Transport Canada and others to develop allowance times for aircraft operating in ice-pellet precipitation. Wind-tunnel testing has been carried out to better understand the flowoff characteristics and resulting aerodynamic effects of anti-icing fluids contaminated with ice pellets using a thin, high-performance wing section at the National Research Council of Canada Propulsion and Icing Wind Tunnel. The objective of this paper is to characterize the aerodynamic behavior of this wing section in order to better understand the adverse aerodynamic effects of anti-icing fluids and ice-pellet contamination. Aerodynamic performance data, boundary-layer surveys and flow visualization were conducted at a Reynolds number of approximately 6.0×10(exp 6) and a Mach number of 0.12. The clean, baseline model exhibited leading-edge stall characteristics including a leading-edge laminar separation bubble and minimal or no separation on the trailing edge of the main element or flap. These results were consistent with expected 2-D aerodynamics and showed no anomalies that could adversely affect the evaluation of anti-icing fluids and ice-pellet contamination on the wing. Tests conducted with roughness and leading-edge flow disturbances helped to explain the aerodynamic impact of the anti-icing fluids and contamination. The stalling characteristics of the wing section with fluid and contamination appear to be driven at least partially by the effects of a secondary wave of fluid that forms near the leading edge as the wing is rotated in the simulated takeoff profile. These results have provided a much more complete understanding of the adverse aerodynamic effects of anti-icing fluids and ice-pellet contamination on this wing section. This is important since these results are used, in part, to develop the ice-pellet allowance times that are applicable to many different airplanes.

Broeren, Andy P.↗

NASA NTRS2013

Aerodynamic Characterization of a Thin, High-Performance Airfoil for Use in Ground Fluids Testing

The FAA has worked with Transport Canada and others to develop allowance times for aircraft operating in ice-pellet precipitation. Wind-tunnel testing has been carried out to better understand the flowoff characteristics and resulting aerodynamic effects of anti-icing fluids contaminated with ice pellets using a thin, high-performance wing section at the National Research Council of Canada Propulsion and Icing Wind Tunnel. The objective of this paper is to characterize the aerodynamic behavior of this wing section in order to better understand the adverse aerodynamic effects of anti-icing fluids and ice-pellet contamination. Aerodynamic performance data, boundary-layer surveys and flow visualization were conducted at a Reynolds number of approximately 6.0 x 10(exp 6) and a Mach number of 0.12. The clean, baseline model exhibited leading-edge stall characteristics including a leading-edge laminar separation bubble and minimal or no separation on the trailing edge of the main element or flap. These results were consistent with expected 2-D aerodynamics and showed no anomalies that could adversely affect the evaluation of anti-icing fluids and ice-pellet contamination on the wing. Tests conducted with roughness and leading-edge flow disturbances helped to explain the aerodynamic impact of the anti-icing fluids and contamination. The stalling characteristics of the wing section with fluid and contamination appear to be driven at least partially by the effects of a secondary wave of fluid that forms near the leading edge as the wing is rotated in the simulated takeoff profile. These results have provided a much more complete understanding of the adverse aerodynamic effects of anti-icing fluids and ice-pellet contamination on this wing section. This is important since these results are used, in part, to develop the ice-pellet allowance times that are applicable to many different airplanes.

Broeren, Andy P.↗

NASA NTRS2002

Implementation of the NCAR Community Land Model (CLM) in the NASA/NCAR finite-volume Global Climate Model (fvGCM)

In this study, the NCAR CLM version 2.0 land-surface model was integrated into the NASA/NCAR fvGCM. The CLM was developed collaboratively by an open interagency/university group of scientists and based on well-proven physical parameterizations and numerical schemes that combine the best features of BATS, NCAR-LSM, and IAP94. The CLM design is a one-dimensional point model with 1 vegetation layer, along with sub-grid scale tiles. The features of the CLM include 10-uneven soil layers with water, ice, and temperature states in each soil layer, and five snow layers, with water flow, refreezing, compaction, and aging allowed. In addition, the CLM utilizes two-stream canopy radiative transfer, the Bonan lake model and topographic enhanced streamflow based on TOPMODEL. The DAO fvGCM uses a genuinely conservative Flux-Form Semi-Lagrangian transport algorithm along with terrain- following Lagrangian control-volume vertical coordinates. The physical parameterizations are based on the NCAR Community Atmosphere Model (CAM-2). For our purposes, the fvGCM was run at 2 deg x 2.5 deg horizontal resolution with 55 vertical levels. The 10-year climate from the fvGCM with CLM2 was intercompared with the climate from fvGCM with LSM, ECMWF and NCEP. We concluded that the incorporation of CLM2 did not significantly impact the fvGCM climate from that of LSM. The most striking difference was the warm bias in the CLM2 surface skin temperature over desert regions. We determined that the warm bias can be partially attributed to the value of the drag coefficient for the soil under the canopy, which was too small resulting in a decoupling between the ground surface and the canopy. We also discovered that the canopy interception was high compared to observations in the Amazon region. A number of experiments were then performed focused on implementing model improvements. In order to correct the warm bias, the drag coefficient for the soil under the canopy was considered a function of LAI (Leaf Area Index). Analysis of the results revealed that there was a substantial impact, and the warm and dry bias in the CLM2 was significantly reduced. For the interception scheme, the canopy throughfall was increased to allow for more infiltration of precipitation into the soil, resulting in increased low-level moisture and a decrease in the interception loss ratio (canopy evaporation to precipitation).

Radakovich, Jon D.↗

NASA NTRS2022

Ambiguous Stability of Glaciers at Bed Peaks

Increasing ice flux from glaciers retreating over deepening (retrograde) bed topography has been implicated in the recent acceleration of mass loss from the Greenland and Antarctic ice sheets. We show in observations that some glaciers have remained at peaks in bed topography without retreating despite enduring significant changes in climate. Observations also indicate that some glaciers which persist at bed peaks undergo sudden retreat years or decades after the onset of local ocean or atmospheric warming. Using model simulations, we show that persistence of a glacier at a bed peak is caused by ice slowing as it flows up a reverse-sloping bed to the peak. Persistence at bed peaks may lead to two very different future behaviors for a glacier: one where it persists at a bed peak indefinitely, and another where it retreats from the bed peak after potentially long delays following climate forcing. However, it is nearly impossible to distinguish which of these two future behaviors will occur from current observations. We conclude that inferring glacier stability from observations of persistence obscures our true commitment to future sea-level rise under climate change. We recommend that further research is needed on seemingly stable glaciers to determine their likely future.

ice dynamics↗

DOE OSTI2024

Glaciation of mixed-phase clouds: insights from bulk model and bin-microphysics large-eddy simulation informed by laboratory experiment

Abstract. Mixed-phase clouds affect precipitation and radiation differently from liquid and ice clouds, posing greater challenges to their representation in numerical simulations. Recent laboratory experiments using the Pi Cloud Chamber explored cloud glaciation conditions based on increased injection of ice-nucleating particles. In this study, we use two approaches to reproduce the results of the laboratory experiments: a bulk scalar mixing model and large-eddy simulation (LES) with bin microphysics. The first approach assumes a well-mixed domain to provide an efficient assessment of the mean cloud properties for a wide range of conditions. The second approach resolves the energy-carrying turbulence, the particle size distribution, and their spatial distribution to provide more details. These modeling approaches enable a separate and detailed examination of liquid and ice properties, which is challenging in the laboratory. Both approaches demonstrate that, with an increased ice number concentration, the flow and microphysical properties exhibit the same changes in trends. Additionally, both approaches show that the ice integral radius reaches the theoretical glaciation threshold when the cloud is subsaturated with respect to liquid water. The main difference between the results of the two approaches is that the bulk model allows for the complete glaciation of the cloud. However, LES reveals that, in a dynamic system, the cloud is not completely glaciated as liquid water droplets are continuously produced near the warm lower boundary and subsequently mixed into the chamber interior. These results highlight the importance of the ice mass fraction in distinguishing the mixed-phase clouds and ice clouds.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2001

Rheology of Diabase: Implications for Tectonics on Venus and Mars

Two important goals of our experimental investigation of the rheological behavior of diabase rocks were: (1) to determine flow laws describing their creep behavior over wide ranges of temperature, stress and strain rate and (2) to develop an understanding of the physical mechanisms by which these rocks flow under laboratory conditions. With this basis, a primary objective then was to construct constitutive equations that can be used to extrapolate from laboratory to planetary conditions. We specifically studied the rheological properties of both natural rock samples and synthetic aggregates. The former provided constraints for geologic systems, while the latter defined the relative contributions of the constituent mineral phases and avoided the influence of glass/melt found in natural samples. In addition, partially molten samples of crustal rock composition were deformed in shear to large strains (greater than 200%) important in crustal environments. The results of this research yielded essential rheological properties essential for models of crustal deformation on terrestrial planets, specifically Venus and Mars, as well as on the geodynamical evolution of these planets. Over the past three years, we also completed our investigation of the creep behavior of water ice with applications to the glaciers, ice sheets and icy satellites. Constitutive equations were determined that describe flow over a wide ranged of stress, strain rate, grain size and temperature. In the case of ice, three creep regimes were delineate. Extrapolation demonstrates that dislocation glide and grain boundary sliding processes dominate flow in ice I under planetary conditions and that diffusion creep is not an important deformation mechanism either in the laboratory or on icy satellites. These results have already been incorporated by other investigators into models describing, for example, the thickness and stability of the ice shell on Europa and to unravel long-standing discrepancies between field observations on glaciers and laboratory results.

Kohlstedt, David L.↗

NASA NTRS2004

A Computational Icing Effects Study for a Three-Dimensional Wing: Comparison with Experimental Data and Investigation of Spanwise Variation

Steady-state solutions of the Reynolds-averaged Navier-Stokes (RANS) equations were computed using the Colbalt flow solver for a constant-section, rectangular wing based on an extruded two-dimensional glaze ice shape. The one equation Spalart-Allmaras turbulence model was used. The results were compared with data obtained from a recent wind tunnel test. Computed results indicate that the steady RANS solutions do not accurately capture the recirculating region downstream of the ice accretion, even after a mesh refinement. The resulting predicted reattachment is farther downstream than indicated by the experimental data. Additionally, the solutions computed on a relatively coarse baseline mesh had detailed flow characteristics that were different from those computed on the refined mesh or the experimental data. Steady RANS solutions were also computed to investigate the effects of spanwise variation in the ice shape. The spanwise variation was obtained via a bleeding function that merged the ice shape with the clean wing using a sinusoidal spanwise variation. For these configurations, the results predicted for the extruded shape provided conservative estimates for the performance degradation of the wing. Additionally, the spanwise variation in the ice shape and the resulting differences in the flow fields did not significantly change the location of the primary reattachment.

Thompson, D.↗

NASA NTRS2016

Modeling Antarctic Subglacial Lake Filling and Drainage Cycles

The growth and drainage of active subglacial lakes in Antarctica has previously been inferred from analysis of ice surface altimetry data. We use a subglacial hydrology model applied to a synthetic Antarctic ice stream to examine internal controls on the filling and drainage of subglacial lakes. Our model outputs suggest that the highly constricted subglacial environment of our idealized ice stream, combined with relatively high rates of water flow funneled from a large catchment, can combine to create a system exhibiting slow-moving pressure waves. Over a period of years, the accumulation of water in the ice stream onset region results in a buildup of pressure creating temporary channels, which then evacuate the excess water. This increased flux of water beneath the ice stream drives lake growth. As the water body builds up, it steepens the hydraulic gradient out of the overdeepened lake basin and allows greater flux. Eventually this flux is large enough to melt channels that cause the lake to drain. Lake drainage also depends on the internal hydrological development in the wider system and therefore does not directly correspond to a particular water volume or depth. This creates a highly temporally and spatially variable system, which is of interest for assessing the importance of subglacial lakes in ice stream hydrology and dynamics.

drainage↗

NASA NTRS2014

Computational Aerodynamic Analysis of Three-Dimensional Ice Shapes on a NACA 23012 Airfoil

The present study identifies a process for performing computational fluid dynamic calculations of the flow over full three-dimensional (3D) representations of complex ice shapes deposited on aircraft surfaces. Rime and glaze icing geometries formed on a NACA23012 airfoil were obtained during testing in the NASA Glenn Research Centers Icing Research Tunnel (IRT). The ice shape geometries were scanned as a cloud of data points using a 3D laser scanner. The data point clouds were meshed using Geomagic software to create highly accurate models of the ice surface. The surface data was imported into Pointwise grid generation software to create the CFD surface and volume grids. It was determined that generating grids in Pointwise for complex 3D icing geometries was possible using various techniques that depended on the ice shape. Computations of the flow fields over these ice shapes were performed using the NASA National Combustion Code (NCC). Results for a rime ice shape for angle of attack conditions ranging from 0 to 10 degrees and for freestream Mach numbers of 0.10 and 0.18 are presented. For validation of the computational results, comparisons were made to test results from rapid-prototype models of the selected ice accretion shapes, obtained from a separate study in a subsonic wind tunnel at the University of Illinois at Urbana-Champaign. The computational and experimental results were compared for values of pressure coefficient and lift. Initial results show fairly good agreement for rime ice accretion simulations across the range of conditions examined. The glaze ice results are promising but require some further examination.

Aerodynamics↗

NASA NTRS2014

Computational Aerodynamic Analysis of Three-Dimensional Ice Shapes on a NACA 23012 Airfoil

The present study identifies a process for performing computational fluid dynamic calculations of the flow over full three-dimensional (3D) representations of complex ice shapes deposited on aircraft surfaces. Rime and glaze icing geometries formed on a NACA23012 airfoil were obtained during testing in the NASA Glenn Research Center's Icing Research Tunnel (IRT). The ice shape geometries were scanned as a cloud of data points using a 3D laser scanner. The data point clouds were meshed using Geomagic software to create highly accurate models of the ice surface. The surface data was imported into Pointwise grid generation software to create the CFD surface and volume grids. It was determined that generating grids in Pointwise for complex 3D icing geometries was possible using various techniques that depended on the ice shape. Computations of the flow fields over these ice shapes were performed using the NASA National Combustion Code (NCC). Results for a rime ice shape for angle of attack conditions ranging from 0 to 10 degrees and for freestream Mach numbers of 0.10 and 0.18 are presented. For validation of the computational results, comparisons were made to test results from rapid-prototype models of the selected ice accretion shapes, obtained from a separate study in a subsonic wind tunnel at the University of Illinois at Urbana-Champaign. The computational and experimental results were compared for values of pressure coefficient and lift. Initial results show fairly good agreement for rime ice accretion simulations across the range of conditions examined. The glaze ice results are promising but require some further examination.

Aerodynamics↗