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 199 records · Page 11

Plans and Preliminary Results of Fundamental Studies of Ice Crystal Icing Physics in the NASA Propulsion Systems Laboratory

This paper describes plans and preliminary results for using the NASA Propulsion Systems Lab (PSL) to experimentally study the fundamental physics of ice-crystal ice accretion. NASA is evaluating whether this facility, in addition to full-engine and motor-driven-rig tests, can be used for more fundamental ice-accretion studies that simulate the different mixed-phase icing conditions along the core flow passage of a turbo-fan engine compressor. The data from such fundamental accretion tests will be used to help develop and validate models of the accretion process. This paper presents data from some preliminary testing performed in May 2015 which examined how a mixed-phase cloud could be generated at PSL using evaporative cooling in a warmer-than-freezing environment.

Ice-Crystal Ice Accretion↗

Plans and Preliminary Results of Fundamental Studies of Ice Crystal Icing Physics in the NASA Propulsion Systems Laboratory

This presentation accompanies the paper titled Plans and Preliminary Results of Fundamental Studies of Ice Crystal Icing Physics in the NASA Propulsion Systems Laboratory. NASA is evaluating whether PSL, in addition to full-engine and motor-driven-rig tests, can be used for more fundamental ice-accretion studies that simulate the different mixed-phase icing conditions along the core flow passage of a turbo-fan engine compressor. The data from such fundamental accretion tests will be used to help develop and validate models of the accretion process. This presentation (and accompanying paper) presents data from some preliminary testing performed in May 2015 which examined how a mixed-phase cloud could be generated at PSL using evaporative cooling in a warmer-than-freezing environment.

Ice-Crystal Ice Accretion↗

Low temperature simulation of subliming boundary layer flow in Jupiter atmosphere

A low-temperature approximate simulation for the sublimation of a graphite heat shield under Jovian entry conditions is studied. A set of algebraic equations is derived to approximate the governing equation and boundary conditions, based on order-of-magnitude analysis. Characteristic quantities such as the wall temperature and the subliming velocity are predicted. Similarity parameters that are needed to simulate the most dominant phenomena of the Jovian entry flow are also given. An approximate simulation of the sublimation of the graphite heat shield is performed with an air-dry-ice model. The simulation with the air-dry-ice model may be carried out experimentally at a lower temperature of 3000 to 6000 K instead of the entry temperature of 14,000 K. The rate of graphite sublimation predicted by the present algebraic approximation agrees to the order of magnitude with extrapolated data. The limitations of the simulation method and its utility are discussed.

Chen, C. J.↗

Controls on Greenland Moulin Geometry and Evolution From the Moulin Shape Model

Nearly all meltwater from glaciers and ice sheets is routed englacially through moulins, which collectively comprise approximately 10–14% of the efficient englacial–subglacial hydrologic system. Therefore, the geometry and evolution of moulins has the potential to influence subglacial water pressure variations, ice motion, and the runoff hydrograph delivered to the ocean. We develop the Moulin Shape (MouSh) model, a time-evolving model of moulin geometry. MouSh models ice deformation around a moulin using both viscous and elastic rheologies and models melting within the moulin through heat dissipation from turbulent water flow, both above and below the water line. We force MouSh with idealized and realistic surface melt inputs. Our results show that variations in surface melt change the geometry of a moulin by approximately 30% daily and by over 100% seasonally. These size variations cause observable differences in moulin water storage capacity, moulin water levels, and subglacial channel size compared to a static, cylindrical moulin. Our results suggest that moulins are significant storage reservoirs for meltwater, with storage capacity and water levels varying over multiple timescales. Representing moulin geometry within subglacial hydrologic models would therefore improve their accuracy, especially over seasonal periods or in regions where overburden pressures are high.

Moulin↗

Ice Accretion Roughness Measurements and Modeling

Roughness on aircraft ice accretions is very important to the overall ice accretion process and to the resulting degradation in aircraft aerodynamic performance. Roughness enhances the local convection leading to more rapid ice accumulation rates, and roughness generates local flow perturbations that lead to higher skin friction. This paper presents 1) a review of the developments in ice shape three-dimensional laser scanning developed at NASA Glenn, 2) a review of the approach of McClain and Kreeger employed to characterize ice roughness evolution on an airfoil surface, and 3) a review of the experimental efforts that have been performed over the last five years to characterize, scale, and model ice roughness evolution physics.

SLD↗

Scaling of Lift Degradation Due to Anti-Icing Fluids Based Upon the Aerodynamic Acceptance Test

In recent years, the FAA has worked with Transport Canada, National Research Council Canada (NRC) and APS Aviation, Inc. to develop allowance times for aircraft operations in ice-pellet precipitation. These allowance times are critical to ensure safety and efficient operation of commercial and cargo flights. Wind-tunnel testing with uncontaminated anti-icing fluids and fluids contaminated with simulated ice pellets had been carried out at the NRC Propulsion and Icing Wind Tunnel (PIWT) to better understand the flow-off characteristics and resulting aerodynamic effects. The percent lift loss on the thin, high-performance wing model tested in the PIWT was determined at 8 deg. angle of attack and used as one of the evaluation criteria in determining the allowance times. Because it was unclear as to how performance degradations measured on this model were relevant to an actual airplane configuration, some means of interpreting the wing model lift loss was deemed necessary. This paper describes how the lift loss was related to the loss in maximum lift of a Boeing 737-200ADV airplane through the Aerodynamic Acceptance Test (AAT) performed for fluids qualification. A loss in maximum lift coefficient of 5.24% on the B737-200ADV airplane (which was adopted as the threshold in the AAT) corresponds to a lift loss of 7.3% on the PIWT model at 8 deg. angle of attack. There is significant scatter in the data used to develop the correlation related to varying effects of the anti-icing fluids that were tested and other factors. A statistical analysis indicated the upper limit of lift loss on the PIWT model was 9.2%. Therefore, for cases resulting in PIWT model lift loss from 7.3% to 9.2%, extra scrutiny of the visual observations is required in evaluating fluid performance with contamination.

Broeren, Andy↗

Martian glaciation and the flow of solid CO2

Results are discussed for a series of laboratory deformation experiments that were conducted to determine a flow law for solid CO2 so that a physical basis could be established for considering the possibility of CO2 glaciation in the Martian polar regions. The experimentally determined flow law indicates that (1) the concept of CO2 glacial flow at the Martian poles is not physically unrealistic, (2) the strength of CO2 near its sublimation point is considerably less than the strength of water ice near its melting point, and (3) CO2 is likely to flow over a broad range of temperatures below its melting point. A glacial model for the Martian polar caps is suggested in which CO2-rich layers deposited near the poles flowed outward laterally to relieve high internal shear stresses. It is shown that this model can explain the topography of the polar caps, the uniform layering of the layered deposits, and the general extent of the polar 'sediments'. Some problematic assumptions of the model are noted, including the large amounts of CO2 necessary to move the layered deposits and the need to increase vapor pressure and temperature with depth in the deposits.

Clark, B. R.↗

Long‐Term Support of an Active Subglacial Hydrologic System in Southeast Greenland by Firn Aquifers

The state of the subglacial hydrologic system, which can modify ice motion, is sensitive to the volume and rate of meltwater reaching it. Bare-ice regions rapidly transport meltwater to the bed via moulins, while in certain accumulation-zone regions, meltwater first flows through firn aquifers, which can introduce a substantial delay. We use a subglacial hydrological model forced with idealized meltwater input scenarios to test the effect of this delay on subglacial hydrology. We find that addition of firn-aquifer water to the subglacial system elevates the inland subglacial water pressure while reducing water pressure and enhancing subglacial channelization near the terminus. This effect dampens seasonal variations in subglacial water pressure and may explain regionally anomalous ice-velocity patterns observed in Southeast Greenland. As surface melt rates increase and firn aquifers expand inland, it is crucial to understand how inland drainage of meltwater affects the evolution of the subglacial hydrologic system.

Kristin Poinar↗

Development and Validation of Physics Based Models for Ice Shedding

Calculations for ice accretion and shedding are presented for a model scale rotor in hover. The aerodynamic characteristics of the rotor are first computed using a combined blade element-momentum theory. The effective angles of attack, and the local flow velocity are used within the NASA Glenn solver LEWICE to estimate the collection efficiency. The computed convection efficiency and the surface pressure distribution from a panel method within LEWICE are used to estimate the ice accretion over the rotor blades over an elapsed time interval. Finally, a force balance approach is used to establish shedding events where the centrifugal force over the ice mass exceeds the adhesive forces at the rotor surface and the cohesive forces between adjacent masses of ice. Preliminary comparisons with test data acquired at the Pennsylvania State Icing Research Tunnel are presented. Sensitivity of the ice shedding events to surface roughness, adhesive strength, cohesive strength, and ambient conditions is discussed.

Gupta, Avani↗

Atmospheric boundary layer modification in the marginal ice zone

A case study of the Andreas et al. (1984) data on atmospheric boundary layer modification in the marginal ice zone is made. The model is a two-dimensional, multilevel, linear model with turbulence, lateral and vertical advection, and radiation. Good agreement between observed and modeled temperature cross sections is obtained. In contrast to the hypothesis of Andreas et al., the air flow is found to be stable to secondary circulations. Adiabatic lifting and, at long fetches, cloud top longwave cooling, not an air-to-surface heat flux, dominate the cooling of the boundary layer. The accumulation with fetch over the ice of changes in the surface wind field is shown to have a large effect on estimates of the surface wind stress. It is speculated that the Andreas et al. estimates of the drag coefficient over the compact sea ice are too high.

Bennett, Theodore J., Jr.↗

Numerical Investigation of Heat Transfer on the NASA SIDRM Model

NASA has developed a Simulated Inter-Duct Research Model (SIDRM) to study ice crystal icing in aircraft engines. A numerical investigation has been completed using the NASA in-house code, GlennHT. For each flow condition, many different wall thermal boundary conditions were imposed. Six uniform wall temperature cases were generated, with wall temperature relative to the freestream total temperature of 0.85, 0.90, 0.95, 1.05, 1.10, and 1.15. In addition, five cases with uniform heat flux were generated. Heat fluxes were chosen to generate average wall temperatures in the range explored by the uniform wall temperature cases. Included in the uniform heat flux cases is the adiabatic case where heat flux equals zero. The results of these simulations are then used to quantify the degree to which the heat transfer coefficient depends on wall thermal boundary condition. It is seen that the heat transfer coefficient is fairly independent of wall thermal boundary condition, as expected. However, it would appear that, at least in some locations, it may be prudent to attempt to correlate a variation of heat transfer coefficient with wall thermal boundary condition.

Ice Crystal Icing↗

Simulation Modeling Requirements for Loss-of-Control Accident Prevention of Turboprop Transport Aircraft

In-flight loss of control remains the leading contributor to aviation accident fatalities, with stall upsets being the leading causal factor. The February 12, 2009. Colgan Air, Inc., Continental Express flight 3407 accident outside Buffalo, New York, brought this issue to the forefront of public consciousness and resulted in recommendations from the National Transportation Safety Board to conduct training that incorporates stalls that are fully developed and develop simulator standards to support such training. In 2010, Congress responded to this accident with Public Law 11-216 (Section 208), which mandates full stall training for Part 121 flight operations. Efforts are currently in progress to develop recommendations on implementation of stall training for airline pilots. The International Committee on Aviation Training in Extended Envelopes (ICATEE) is currently defining simulator fidelity standards that will be necessary for effective stall training. These recommendations will apply to all civil transport aircraft including straight-wing turboprop aircraft. Government-funded research over the previous decade provides a strong foundation for stall/post-stall simulation for swept-wing, conventional tail jets to respond to this mandate, but turboprops present additional and unique modeling challenges. First among these challenges is the effect of power, which can provide enhanced flow attachment behind the propellers. Furthermore, turboprops tend to operate for longer periods in an environment more susceptible to ice. As a result, there have been a significant number of turboprop accidents as a result of the early (lower angle of attack) stalls in icing. The vulnerability of turboprop configurations to icing has led to studies on ice accumulation and the resulting effects on flight behavior. Piloted simulations of these effects have highlighted the important training needs for recognition and mitigation of icing effects, including the reduction of stall margins. This paper addresses simulation modeling requirements that are unique to turboprop transport aircraft and highlights the growing need for aerodynamic models suitable for stall training for these configurations. A review of prominent accidents that involved aerodynamic stall is used to illustrate various modeling features unique to turboprop configurations and the impact of stall behavior on susceptibility to loss of control that has led to new training requirements. This is followed by an overview of stability and control behavior of straight-wing turboprops, the related aerodynamic characteristics, and a summary of recent experimental studies on icing effects. In addition, differences in flight dynamics behavior between swept-wing jets and straight-wing turboprop configurations are discussed to compare and contrast modeling requirements. Specific recommendations for aerodynamic models along with further research needs and data measurements are also provided. 1

Crider, Dennis↗

Evolution and Transport of Water in the Upper Regolith of Mars

Long standing theoretical predictions [1-3], as well as recent spacecraft observations [4] indicate that large quantities of ice is present in the high latitudes upper decimeters to meters of the Martian regolith. At shallower depths and warmer locations small amounts of H2O, either adsorbed or free, may be present transiently. An understanding of the evolution of water based on theoretical and experimental considerations of the processes operating at the Martian environment is required. In particular, the porosity, diffusivity, and permeability of soils and their effect on water vapor transport under Mars-like conditions have been estimated, but experimental validation of such models is lacking. Goal: Three related mechanisms may affect water transport in the upper Martian regolith. 1) diffusion along a concentration gradient under isobaric conditions, 2) diffusion along a thermal gradient, which may give rise to a concentration gradient as ice sublimes or molecules desorb from the regolith, and 3) hydraulic flow, or mass motion in response to a pressure gradient. Our combined theoretical and experimental investigation seeks to disentangle these mechanisms and determine which process(es) are dominant in the upper regolith over various timescales. A detailed one-dimensional model of the upper regolith is being created which incorporates water adsorption/ desorption, condensation, porosity, diffusivity, and permeability effects. Certain factors such as diffusivity are difficult to determine theoretically due to the wide range of intrinsic grain properties such as particle sizes, shapes, packing densities, and emergent properties such as tortuosity. An experiment is being designed which will allow us to more accurately determine diffusivity, permeability, and water desorption isotherms for regolith simulants.

Hudson, T. L.↗

Hybrid Glacier Inventory, Gravimetry and Altimetry (HIGA) Mass Balance Product for Greenland and the Canadian Arctic

We present a novel inversion algorithm that generates a mass balance field that is simultaneously consistent with independent observations of glacier inventory derived from optical imagery, cryosphere-attributed mass trends derived from satellite gravimetry, and ice surface elevation trends derived from airborne and satellite altimetry. We use this algorithm to assess mass balance across Greenland and the Canadian Arctic over the Sep-2003 to Oct- 2009 period at 26 kilometers resolution. We evaluate local algorithm-inferred mass balance against forty in situ point observations. This evaluation yields a root mean squared error (RMSE) of 0.15 mWE/a ( 0.15 meters (water equivalent) per annum), and highlights a paucity of in situ observations from regions of high dynamic mass loss and peripheral glaciers. We assess mass losses of 212 plus or minus 67 Gigatons per annum to the Greenland ice sheet proper, 38 Gigatons per annum to peripheral glaciers in Greenland, and 42 Gigatons per annum to glaciers in the Canadian Arctic. These magnitudes of mass loss are dependent on the gravimetry-derived spherical harmonic mass trend we invert. We spatially partition the transient glacier continuity equation by differencing algorithm inferred mass balance from modeled surface mass balance, in order to solve the horizontal divergence of ice flux as a residual. This residual ice dynamic field infers flux divergence (or submergent flow) in the ice sheet accumulation area and at tidewater margins, and flux convergence (or emergent flow) in land-terminating ablation areas, which is consistent with continuum mechanics theory.

Altimetry↗

Ice-Crystal Icing Accretion Studies at the NASA Propulsion Systems Laboratory

This paper describes an ice-crystal icing experiment conducted at the NASA Propulsion System Laboratory during June 2018. This test produced ice shape data on an airfoil for different test conditions similar to those inside the compressor region of a turbo-fan jet engine. Mixed-phase icing conditions were generated by partially freezing out a water spray using the relative humidity of flow as the primary parameter to control freeze-out. The paper presents the ice shape data and associated conditions which include pressure, velocity, temperature, humidity, total water content, melt ratio, and particle size distribution. The test featured a new instrument traversing system which allowed surveys of the flow and cloud. The purpose of this work was to provide experimental ice shape data and associated conditions to help develop and validate ice-crystal icing accretion models. The results support previous experimental observations of a minimum melt-ratio threshold for accretion to occur as well as the existence of a plateau region where the icing severity is high for a range of melt ratios. However, a maximum limit for melt ratio, which is suggested in the ice crystal icing literature, was not observed perhaps complicated by the potential for some supercooling of the water at these conditions.

Engine Icing↗

Numerical Simulations of Non-Newtonian Convection in Ice: Application to Europa

Numerical simulations of solid state convection in Europa's ice shell have so far been limited to consideration of Newtonian flow laws, where the viscosity of ice is strongly dependent upon temperature, predicting that a stagnant lid should form at the top (10-40%) of a convecting ice shell. Such large thicknesses seem to contradict estimates of the effective elastic thickness of Europa s ice shell during its geologically active period. Recent laboratory experiments characterize the rheology of ice as the sum of contributions from several temperature and strain rate-dependent creep mechanisms. We present the results of numerical simulations of convection within Europa s ice shell using the finite-element model Citcom, applying the non-Newtonian rheology of grain boundry sliding. Our calculations suggest a shallower brittle/ductile transition and larger interior convective velocities compared to Newtonian rheology. The flow field is time-dependent, with small, localized upwellings and downwellings at the thermal boundary layers that have minimal topographic expression at the surface.

Barr, A. C.↗

SSME fuel preburner phase 1/phase 2+ comparison

An analytical comparison of the predicted combustion flowfields for the Phase 1 and Phase 2+ fuel preburner configurations using an advanced CFD combustion code is currently underway at Rocketdyne. The Phase 2+ injector has a modified element design. The redesigned element provides for improved atomization through a shift in the injector element/face bleed flow split while maintaining an acceptable pressure drop. Cold-flow tests of the new design have been encouraging. As a complement to those tests, the Advanced Rocket Injector Combustor Code (ARICC) was selected to conduct a detailed combustion flowfield simulation. ARICC models two-dimensional (or axisymmetric), transient, turbulent, two-phase mixing and combustion flowfields. It is derived from a member of the Los Alamos ICEd-ALE family of codes. Unique features of ARICC include explicit representation of the coaxial LOX jet and hydrogen gas flows, distributed atomization processes, droplet breakup, and supercritical vaporation processes. Predictions of temperature and OH concentration profiles for the Phase 1 and Phase 2+ injector flowfields indicate a trend toward more uniform temperature distributions and shorter flame lengths with the Phase 2+ design.

Jensen, R. J.↗

CFD zonal modeling of leading-edge ice effects for a complete aircraft

A simplified, uncoupled zonal procedure was utilized to assess the capability of numerically simulating icing effects on a Boeing 727-200 aircraft. The computational approach combines potential flow, plus boundary layer simulations by VSAERO for the un-iced aircraft forces and moments, with Navier-Stokes simulations by ARC3D for the incremental forces and moments due to iced components. These are compared with wind tunnel longitudinal force and moment data. Although the computational results compared favorably with the test data in the linear angle of attack range, it is clear that for general aircraft icing calculations, a multiblock Navier-Stokes code will be required for the viscous component of this zonal method.

Summa, J. M.↗