Search NASA⌕ Search

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.

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

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Flow Quality Measurements in an Aerodynamic Model of NASA Lewis' Icing Research Tunnel

As part of an ongoing effort to improve the aerodynamic flow characteristics of the Icing Research Tunnel (IRT), a modular scale model of the facility was fabricated. This 1/10th-scale model was used to gain further understanding of the flow characteristics in the IRT. The model was outfitted with instrumentation and data acquisition systems to determine pressures, velocities, and flow angles in the settling chamber and test section. Parametric flow quality studies involving the insertion and removal of a model of the IRT's distinctive heat exchanger (cooler) and/or of a honeycomb in the settling chamber were performed. These experiments illustrate the resulting improvement or degradation in flow quality.

Canacci, Victor A.↗
NASA NTRS2002

Ice Shelves and Landfast Ice on the Antarctic Perimeter: Revised Scope of Work

Ice shelves respond quickly and profoundly to a warming climate. Within a decade after mean summertime temperature reaches approx. O C and persistent melt pending is observed, a rapid retreat and disintegration occurs. This link was documented for ice shelves in the Antarctic Peninsula region (the Larsen 'A', 'B' and Wilkins Ice shelves) by the results of a previous grant under ADRO-1. Modeling of ice flow and the effects of meltwater indicated that melt pending accelerates shelf breakup by increasing fracture penetration. SAR data supplemented an AVHRR- and SSM/I-based image analysis of extent and surface characteristic changes. This funded grant is a revised, scaled-down version of an earlier proposal under the ADRO-2 NRA. The overall objective remains the same: we propose to build on the previous study by examining other ice shelves of the Antarctic and incorporate an examination of the climate-related characteristics of landfast ice. The study now considers just a few shelf and fast ice areas for study, and is funded for two years. The study regions are the northeastern Ross Ice Shelf, the Larsen 'B' and 'C' shelves, fast ice and floating shelf ice in the Pine Island Glacier area, and fast ice along the Wilkes Land coast. Further, rather than investigating a host of shelf and fast ice processes, we will home in on developing a series of characteristics associated with climate change over shelf and fast ice areas. Melt pending and break-up are the end stages of a response to a warming climate that may begin with increased melt event frequency (which changes both albedo and emissivity temporarily), changing firn backscatter (due to percolation features), and possibly increased rifting of the shelf surface. Fast ice may show some of these same processes on a seasonal timescale, providing insight into shelf evolution.

Scambos, Ted↗
NASA NTRS2019

Physically Based and Stochastic Models for Greenland Moulin Formation, Longevity, and Spatial Distribution

Nearly all proglacial water discharge from the Greenland Ice Sheet is routed englacially, from the surface to the bed, via moulins. Identification of moulins in high-resolution imagery is a frequent topic of study, but the processes controlling how and where moulins form remain poorly understood. We seek to leverage information gained from the development of a physical model of moulin formation, remotely sensed ice-sheet data products, and an analytic model of ice-flow perturbations to develop a predictive stochastic model of moulin distribution across Greenland. Here we present initial results from the physical model of moulin formation and characterize the sensitivity of moulin geometry to a range of model parameters. This parameterization of moulin formation is the first step in developing a stochastic model that will be a predictive, computationally efficient representation of the englacial hydrologic system.

Andrews, Lauren C.↗
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↗
DOE OSTI2025

Computationally efficient subglacial drainage modelling using Gaussian process emulators: GlaDS-GP v1.0

Subglacial drainage models represent water flow at the ice–bed interface through coupled distributed and channelized systems to determine water pressure, discharge, and drainage system geometry. While they are used to understand processes such as the relationship between surface melt and ice flow, the number of uncertain model parameters and the computational cost of running models makes it difficult to adequately explore the high-dimensional parameter space and evaluate uncertainty in model predictions. Here, we develop Gaussian process (GP) emulators that make fast predictions with associated uncertainty of subglacial drainage model outputs. Using a truncated principal component (PC) basis representation, we construct a GP emulator for diurnally averaged subglacial water pressure. We also explore emulation of scalar variables describing drainage efficiency and configuration. We train the emulators using ensembles of up to 512 simulations varying eight parameters of the Glacier Drainage System (GlaDS) model on a synthetic domain intended to represent an ice-sheet margin. The emulators make predictions ∼ 1000 times faster than GlaDS simulations, with errors <3 % for the water pressure field and ∼ 5 %–9 % for drainage efficiency and configuration. We apply the emulators to explore the eight-dimensional parameter space by computing variance-based parameter sensitivity indices, finding that three parameters (ice flow coefficient, bed bump aspect ratio, and the subglacial cavity system conductivity) explain 90 % of the variance in modelled water pressure in response to parameter changes. The GP emulator approach described here is well suited to integrating observational data with models to make calibrated, credible predictions of subglacial drainage.

58 GEOSCIENCES↗
NASA NTRS2022

GBaTSv2: a revised synthesis of the likely basal thermal state of the Greenland Ice Sheet

The basal thermal state (frozen or thawed) of the Greenland Ice Sheet is under-constrained due to few direct measurements, yet knowledge of this state is becoming increasingly important to interpret modern changes in ice flow. The first synthesis of this state relied on inferences from widespread airborne and satellite observations and numerical models, for which most of the underlying datasets have since been updated. Further, new and independent constraints on the basal thermal state have been developed from analysis of basal and englacial reflections observed by airborne radar sounding. Here we synthesize constraints on the Greenland Ice Sheet's basal thermal state from boreholes, thermomechanical ice-flow models that participated in the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6; Coupled Model Intercomparison Project Phase 6), IceBridge BedMachine Greenland v4 bed topography, Making Earth Science Data Records for Use in Research Environments (MEaSUREs) Multi-Year Greenland Ice Sheet Velocity Mosaic v1 and multiple inferences of a thawed bed from airborne radar sounding. Most constraints can only identify where the bed is likely thawed rather than where it is frozen. This revised synthesis of the Greenland likely Basal Thermal State version 2 (GBaTSv2) indicates that 33 % of the ice sheet's bed is likely thawed, 40 % is likely frozen and the remainder (28 %) is too uncertain to specify. The spatial pattern of GBaTSv2 is broadly similar to the previous synthesis, including a scalloped frozen core and thawed outlet-glacier systems. Although the likely basal thermal state of nearly half (46 %) of the ice sheet changed designation, the assigned state changed from likely frozen to likely thawed (or vice versa) for less than 6 % of the ice sheet. This revised synthesis suggests that more of northern Greenland is likely thawed at its bed and conversely that more of southern Greenland is likely frozen, both of which influence interpretation of the ice sheet's present subglacial hydrology and models of its future evolution. The GBaTSv2 dataset, including both code that performed the analysis and the resulting datasets, is freely available at https://doi.org/10.5281/zenodo.6759384 (MacGregor, 2022).

Joseph A. MacGregor↗
NASA NTRS1992

LDV measurements on a rectangular wing with a simulated glaze ice accretion

LDV measurement results are presented for the upper surface of a rectangular semispan wing with and without simulated glaze ice accretion. Inspection of the model centerline flow field indicates that a large region of reverse flow exists aft of the ice horn on the iced model. At alpha = 0 deg, this region extends to 7 percent chord, while at alpha = 4 deg the bubble grows to more than 12 percent chord. At alpha = 8 deg, the time-averaged separation bubble is measured well beyond 50 percent chord. Experimental and computational flow visualization support these findings. The flow in the vicinity of the ice shape contains many of the features of flow over a backward-facing step.

Khodadoust, A.↗
NASA NTRS1993

Ice Accretion Prediction for a Typical Commercial Transport Aircraft

Ice accretion calculations were made for a modern commercial transport using the NASA Lewis LEWICE3D ice accretion code. The ice accretion calculations were made for the wing and horizontal tail using both isolated flow models and flow models incorporating the entire airplane. The isolated flow model calculations were made to assess the validity of using these simplified models in lieu of the entire model in the ice accretion analysis of full aircraft. Ice shapes typifying a rime and a mixed ice shape were generated for a 30 minute hold condition. In general, the calculated ice shapes looked reasonable and appeared representative of a rime and a mixed ice conditions. The isolated flow model simplification was good for the main wing except at the root where it overpredicted the amount of accreted ice relative to the full aircraft flow model. For the horizontal tail the size and amount of predicted ice compared well for the two flow models, but the position of the accretions were more towards the upper surface for the aircraft flow model relative to the isolated flow model. This was attributed to downwash from the main wing which resulted in a lower effective angle of attack for the aircraft tail.

Bidwell, C. S.↗
NASA NTRS1993

Ice accretion prediction for a typical commercial transport aircraft

Ice accretion calculations were made for a modern commercial transport using the NASA Lewis LEWICE3D ice accretion code. The ice accretion calculations were made for the wing and horizonal tail using both isolated flow models and flow models incorporating the entire airplane. The isolated flow model calculations were made to assess the validity of using these simplified models in lieu of the entire model in the ice accretion analysis for full aircraft. Ice shapes typifying a rime and a mixed ice shape were generated for a 30 minute hold condition. In general, the calculated ice shapes looked reasonable and appeared representative of a rime and a mixed ice conditions. The isolated flow model simplification was good for the main wing except at the root where it overpredicted the amount of accreted ice relative to the full aircraft flow model. For the horizontal tail the size and amount of predicted ice compared well for the two flow models, but the position of the accretions were more towards the upper surface for the aircraft flow model relative to the isolated flow model. This was attributed to downwash from the main wing which resulted in a lower effective angle of attack for the aircraft tail.

Bidwell, C. S.↗
NASA NTRS1986

Experimental evidence for modifying the current physical model for ice accretion on aircraft surfaces

Closeup movies, still photographs, and other experimental data suggest that the current physical model for ice accretion needs significant modification. At aircraft airspeeds there was no flow of liquid over the surface of the ice after a short initial flow, even at barely subfreezing temperatures. Instead, there were very large stationary drops on the ice surface that lose water from their bottoms by freezing and replenish their liquid by catching the microscopic cloud droplets. This observation disagrees with the existing physical model, which assumes there is a thin liquid film continuously flowing over the ice surface. With no such flow, the freezing-fraction concept of the model fails when a mass balance is performed on the surface water. Rime ice does, as the model predicts, form when the air temperature is low enough to cause the cloud droplets to freeze almost immediately on impact. However, the characteristic shapes of horn-glaze ice or rime ice are primarily caused by the ice shape affecting the airflow locally and consequently the droplet catch and the resulting ice shape. Ice roughness greatly increases the heat transfer coefficient, stops the movement of drops along the surface, and may also affect the airflow initially and thereby the droplet catch. At high subreezing temperatures the initial flow and shedding of surface drops have a large effect on the ice shape. At the incipient freezing limit, no ice forms.

Olsen, W.↗
NASA NTRS1999

Ice Flow in the Northeast Greenland Ice Stream

As the use of satellite data for icesheet research has accelerated, new discoveries have changed the way that we think about ice discharge, and new tools have changed the way that we measure it. We have begun to apply a combination of remote sensing and ice sheet modeling techniques to study ice flow in the recently identified northeast Greenland ice stream. This ice stream shows evidence of organized flow far into the interior and has low-slope areas of rapid flow and regions of enhanced shear that resemble the ice streams of West Antarctica. The onset of rapid flow close to the ice divide and the evolution of its flow pattern, however, make this ice stream unique. Our investigation is aimed at increasing our understanding of this ice stream and its origins and role in the discharge of ice from the ice sheet and the effects that discharge may have on the history and evolution of the ice in northern Greenland. The ice flow data comes from satellite radar interferometry tied to reference points. We have applied this technique to measure the velocity over the entire ice stream, which is more than 700 km in length. The velocity contours reveal distinct, well-developed shear margins, which are coincident with linear features in the Synthetic Aperture Radar (SAR) amplitude imagery. It is surprising to see such defined, enhanced flow so far inland and so close to the ice divide. The velocity data and topography data are being used to feed a modeling effort that allows us to study the dynamics and character rapid flow. Additional information is contained in the original.

Joughin, Ian↗
NASA NTRS2005

Rheology of the Ronne Ice Shelf, Antarctica, Inferred from Satellite Radar Interferometry Data using an Inverse Control Method

The Antarctic Ice Sheet is surrounded by large floating ice shelves that spread under their own weight into the ocean. Ice shelf rigidity depends on ice temperature and fabrics, and is influenced by ice flow and the delicate balance between bottom and surface accumulation. Here, we use an inverse control method to infer the rigidity of the Ronne Ice Shelf that best matches observations of ice velocity from satellite radar interferometry. Ice rigidity, or flow law parameter B, is shown to vary between 300 and 900 kPa a(sup 1/3). Ice is softer along the side margins due to frictional heating, and harder along the outflow of large glaciers, which advect cold continental ice. Melting at the bottom surface of the ice shelf increases its rigidity, while freezing decreases it. Accurate numerical modelling of ice shelf flow must account for this spatial variability in mechanical characteristics.

ice velocity↗
DOE OSTI2026

An Approach for Simulating Transport of Nonspherical, Sublimating Particles in Rarefied Gas Flows via the Direct Simulation Monte Carlo Method

This paper extends the particle Green’s function approach of Gallis et al. [Phys. Fluids 13, 3482–3492 (2001)] to compute the force and heat transfer on nonspherical, sublimating particles in rarefied gas flows. Green’s functions are derived for rods and plates, which are representative of ice particle shapes, and a numerical procedure is developed for computing the Green’s functions for a general discretized convex shape. Additionally, a sublimation model is developed to simulate the evolution of ice particles in rarefied gas flows. These models are implemented in a direct simulation Monte Carlo code and demonstrated for the case of ice-particle-laden hypersonic flow over a ramp.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

HdsaLib

HdsaLib performs hyper-differential sensitivity analysis in optimization problems that are constrained by partial differential equation problems. The software aids researchers in understanding the impact of uncertainties on their optimization models. The tool has been applied for use in ice sheet modeling, fluid flow modeling, atmospheric tracer transport modeling, and subsurface flow modeling. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

van Bloemen Waanders, Bart [Sandia National Lab. (↗