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 307 records · Page 17

Characteristics of cometary picked-up ions in a global model of Giacobini-Zinner

Energetic ions observed during the International Cometary Explorer (ICE) spacecraft flyby of comet Giacobini-Zinner provide information about both the constitution of comets and the plasma physical processes associated with their interaction with the solar wind. In this investigation the details of ion 'pickup,' in the limit where small-scale fluctuations in the plasma and magnetic field are neglected, are modeled by following the motion of a large number of initially cold, heavy (mass 18) ions in a global magnetohydrodynamic model of the local plasma and magnetic field. The results indicate how the background or macroscopic velocity and magnetic field structure of the comet can affect the average spatial and spectral characteristics of the observed cometary ions. These effects, which occur by virtue of forces associated with the compression and the curvature of the magnetic field in the presence of the stagnating plasma flow, can explain the double maxima in the time series of the energetic ion flux observed along the ICE trajectory.

Kimmel, C. D.↗

Phase transitions and 2D spherical convection in a large icy satellite

Three ice phase transitions, based on experimental and theoretical results, were incorporated into a 2D spherical convection model. These phase transitions are the exothermic Ice 1-2, the endothermic Ice 2-6, and the exothermic Ice 6-8 phase transitions. The fluid is assumed isoviscous and the anelastic liquid formulation was used. The object is an icy satellite whose physical characteristics resemble those of Ganymede or Callisto. Structural models, similar to those of Mueller and McKinnon, with different core sizes, i.e. with different degree of differentiation were studied. The shells, i.e. the mantle of the satellite, is heated both internally and from below to account for the decaying radiogenic heating and the heat flow from the solid core. The lower boundary of the mantle is rigid and isothermal, the upper boundary is isothermal. Calculations with rigid and shear stress free upper boundary are carried out in order to assess the role of the different boundary conditions. Two different Rayleigh numbers, depending on the assumed value of the viscosity, were used in the calculations and the thermal evolution of the satellite was studied. The suit of calculations presented demonstrates that phase transition cannot be ignored when the thermal evolution of a large icy satellite has to be studied. The importance of the construction of the realistic phase diagram for ice was also documented. On the other hand the degree of differentiation, simulated with the structures with different core's radius does not seem to drastically influence the overall behavior of the convective activity. Within the limits of our simulations, the post-accretional evolution seems to be independent on the degree of primordial differentiation.

Forni, O.↗

Interpretation of satellite-observed thunderstorm anvil structure

The radiative characteristics of an anvil plume sensed by satellite radiometers are studied to discern the reason for an observed 'V' pattern of colder temperatures and a corresponding internal warm region. The measured black body temperature is affected by the ice water content and particle distributions in the cloud layer. Attention is given to 'V' anvil structure examined by Haynesfield et al. (1979). The GOES satellite IR data are treated with a plane parallel geometric radiative transfer model and a particle trajectory model is used to simulate ice particle concentrations downwind from a strongly divergent cloud top. Output from the two models indicate the 'V' pattern arises from divergence in the storm anvil. The particle concentration was elevated within the V-shaped region, causing lower temperatures in the V arms. The V is produced by strong upper level flow as well as strong mass and moisture divergence in the storm outflow.

Heymsfield, G. M.↗

Background Error Covariance Estimation Using Information from a Single Model Trajectory with Application to Ocean Data Assimilation

An attractive property of ensemble data assimilation methods is that they provide flow dependent background error covariance estimates which can be used to update fields of observed variables as well as fields of unobserved model variables. Two methods to estimate background error covariances are introduced which share the above property with ensemble data assimilation methods but do not involve the integration of multiple model trajectories. Instead, all the necessary covariance information is obtained from a single model integration. The Space Adaptive Forecast error Estimation (SAFE) algorithm estimates error covariances from the spatial distribution of model variables within a single state vector. The Flow Adaptive error Statistics from a Time series (FAST) method constructs an ensemble sampled from a moving window along a model trajectory.SAFE and FAST are applied to the assimilation of Argo temperature profiles into version 4.1 of the Modular Ocean Model (MOM4.1) coupled to the GEOS-5 atmospheric model and to the CICE sea ice model. The results are validated against unassimilated Argo salinity data. They show that SAFE and FAST are competitive with the ensemble optimal interpolation (EnOI) used by the Global Modeling and Assimilation Office (GMAO) to produce its ocean analysis. Because of their reduced cost, SAFE and FAST hold promise for high-resolution data assimilation applications.

Error Covariance↗

Evidence for Water Ie on the Moon: Results for Anomalous Polar Craters from the LRO Mini-RF Imaging Radar

The Mini-RF radar instrument on the Lunar Reconnaissance Orbiter spacecraft mapped both lunar poles in two different RF wavelengths (complete mapping at 12.6 cm S-band and partial mapping at 4.2 cm X-band) in two look directions, removing much of the ambiguity of previous Earth- and spacecraft-based radar mapping of the Moon's polar regions. The poles are typical highland terrain, showing expected values of radar cross section (albedo) and circular polarization ratio (CPR). Most fresh craters display high values of CPR in and outside the crater rim; the pattern of these CPR distributions is consistent with high levels of wavelength-scale surface roughness associated with the presence of block fields, impact melt flows, and fallback breccia. A different class of polar crater exhibits high CPR only in their interiors, interiors that are both permanently dark and very cold (less than 100 K). Application of scattering models developed previously suggests that these anomalously high-CPR deposits exhibit behavior consistent with the presence of water ice. If this interpretation is correct, then both poles may contain several hundred million tons of water in the form of relatively "clean" ice, all within the upper couple of meters of the lunar surface. The existence of significant water ice deposits enables both long-term human habitation of the Moon and the creation of a permanent cislunar space transportation system based upon the harvest and use of lunar propellant.

polar↗

Warm-Based Glacial Landforms and Past Climate Signals: Icelandic Eskers as Proxies for Small Sinuous Ridges on Mars

The surface of Mars exhibits plentiful geomorphic evidence for glacial processes, but the history of the cryosphere on Mars remains debated [e.g., 1, 2]. Recent hydrological modeling predicts that warm-based ice sheets on Mars may have resulted in subglacial incised channels (e.g., valley networks) and eskers [3]. Evidence on Mars for warm-based glaciation includes candidate moraines, bedrock channels, and candidate eskers [e.g., 4-6]. Eskers are sinuous ridges comprising glaciofluvial sediments deposited by meltwater flowing through tunnels at the bed of warm-based large glaciers [e.g., 7, 8]. They provide direct evidence of basal melting of ice and preserve glaciofluvial sediments. Warm-based glacial transport is primarily driven by meltwater and can chemically weather the underlying bedrock. An alternative signature of past climate on Mars may thus be the geochemical record, due to alteration by interactions with glacial meltwater, preserved in esker sediments and strata. At present, the sedimentology and stratigraphy of terrestrial eskers is poorly understood. This study will combine sedimentology, geochemistry, and geomorphology to better constrain the unique properties of eskers resulting from warm-based glaciation on Mars-like substrates by leveraging a natural laboratory in Iceland.

A. M. Rutledge↗

Correlations Among Ice Measurements, Impingement Rates Icing Conditions, and Drag Coefficients for Unswept NACA 65A004 Airfoil

An empirical relation has been obtained by which the change in drag coefficient caused by ice formations on an unswept NACA 65AO04 airfoil section can be determined from the following icing and operating conditions: icing time, airspeed, air total temperature, liquid-water content, cloud droplet impingement efficiencies, airfoil chord length, and angles of attack. The correlation was obtained by use of measured ice heights and ice angles. These measurements were obtained from a variety of ice formations, which were carefully photographed, cross-sectioned, and weighed. Ice weights increased at a constant rate with icing time in a rime icing condition and at progressively increasing rates in glaze icing conditions. Initial rates of ice collection agreed reasonably well with values predicted from droplet impingement data. Experimental droplet impingement rates obtained on this airfoil section agreed with previous theoretical calculations for angles of attack of 40 or less. Disagreement at higher angles of attack was attributed to flow separation from the upper surface of the experimental airfoil model.

Gray, Vernon H.↗

Background Error Covariance Estimation using Information from a Single Model Trajectory with Application to Ocean Data Assimilation into the GEOS-5 Coupled Model

An attractive property of ensemble data assimilation methods is that they provide flow dependent background error covariance estimates which can be used to update fields of observed variables as well as fields of unobserved model variables. Two methods to estimate background error covariances are introduced which share the above property with ensemble data assimilation methods but do not involve the integration of multiple model trajectories. Instead, all the necessary covariance information is obtained from a single model integration. The Space Adaptive Forecast error Estimation (SAFE) algorithm estimates error covariances from the spatial distribution of model variables within a single state vector. The Flow Adaptive error Statistics from a Time series (FAST) method constructs an ensemble sampled from a moving window along a model trajectory. SAFE and FAST are applied to the assimilation of Argo temperature profiles into version 4.1 of the Modular Ocean Model (MOM4.1) coupled to the GEOS-5 atmospheric model and to the CICE sea ice model. The results are validated against unassimilated Argo salinity data. They show that SAFE and FAST are competitive with the ensemble optimal interpolation (EnOI) used by the Global Modeling and Assimilation Office (GMAO) to produce its ocean analysis. Because of their reduced cost, SAFE and FAST hold promise for high-resolution data assimilation applications.

Data Assimilation↗

Greenland Ice Sheet Mass Balance: Distribution of Increased Mass Loss with Climate Warming; 2003-07 Versus 1992-2002

We derive mass changes of the Greenland ice sheet (GIS) for 2003-07 from ICESat laser altimetry and compare them with results for 1992-2002 from ERS radar and airborne laser altimetry. The GIS continued to grow inland and thin at the margins during 2003 07, but surface melting and accelerated flow significantly increased the marginal thinning compared with the 1990s. The net balance changed from a small loss of 7 plus or minus 3 Gt a 1(sup -1) in the 1990s to 171 plus or minus 4 Gt a (sup -1) for 2003-07, contributing 0.5 mm a(sup -1) to recent global sea-level rise. We divide the derived mass changes into two components: (1) from changes in melting and ice dynamics and (2) from changes in precipitation and accumulation rate. We use our firn compaction model to calculate the elevation changes driven by changes in both temperature and accumulation rate and to calculate the appropriate density to convert the accumulation-driven changes to mass changes. Increased losses from melting and ice dynamics (17-206 Gt a(sup-1) are over seven times larger than increased gains from precipitation (10 35 Gt a(sup-1) during a warming period of approximately 2 K (10 a)(sup -1) over the GIS. Above 2000m elevation, the rate of gain decreased from 44 to 28 Gt a(sup-1), while below 2000m the rate of loss increased from 51 to 198 Gt a(sup-1). Enhanced thinning below the equilibrium line on outlet glaciers indicates that increased melting has a significant impact on outlet glaciers, as well as accelerating ice flow. Increased thinning at higher elevations appears to be induced by dynamic coupling to thinning at the margins on decadal timescales.

Zwally, H. Jay↗

The origin of fluvial valleys and early geologic history, Aeolis Quadrangle, Mars

A valley genesis model of the southern Aeolis Quadrangle, eastern Mars, is proposed. In response to widespread effusive volcanism in interstratified, ice-reach terrains, local subsidence occurred along fractures/faults, and produced scarps that intersected local aquifers. Heated spring discharges issuing from these scarps may have carved the valleys through a combination of downvalley water, ice, and debris flows, and headward erosion along stratal discontinuities and individual conduit faults and fractures. Fifty-six asymmetric scarps or ridges that are probable thrust faults are mapped. The faults exhibit an orientation vector mean of N 63 deg W +/- 11 deg and they transect the lava plains and the older plateau sequence units. The vector mean for the 264 valleys mapped is N 48 deg W +/- 12 deg, with a larger dispersion about the mean. The similar orientations displayed by thrust faults and valley axes suggest that valley locations are partly controlled by preexisting thrust faults and fracture systems.

Brakenridge, G. Robert↗

Insolation and the Precession Index

Simple nonlinear climate models yield a precession index-like term in the temperature. Despite its importance in the geologic record, the precession index e sin omega, where e is the Earth's orbital eccentricity and omega is the Sun's perigee in the geocentric frame, is not present in the insolation at the top of the atmosphere. Hence there is no one-for-one mapping of 23,000 and 19,000 year periodicities from the insolation to the paleoclimate record; a nonlinear climate model is needed to produce these periods. Two such models, a grey body and an energy balance climate model with an added quadratic term, produce e sin omega terms in temperature. These terms, which without feedback mechanisms achieve extreme values of about plus or minus 0.48 K for the grey body and plus or minus 0.64 K for the energy balance model, simultaneously cool one hemisphere while they warm the other. Moreover, they produce long-term cooling in the northern hemisphere when the Sun's perigee is near northern solstice and long-term warming in the northern hemisphere when the perigee is near southern solstice. Thus this seemingly paradoxical mechanism works against the standard model which requires cool northern summers (Sun far from Earth in northern summer) to build up northern ice sheets, so that if the standard model is correct it may be more efficient than previously thought. Alternatively, the new mechanism could possibly be dominant and indicate southern hemisphere control of the northern ice sheets, wherein the southern oceans undergo a long-term cooling when the Sun is close to the Earth during southern summer. The cold water eventually flows north, cooling the northern hemisphere. This might explain why the northern oceans lag the southern ones when it comes to orbital forcing.

Rubincam, David Parry↗

The Precession Index and a Nonlinear Energy Balance Climate Model

A simple nonlinear energy balance climate model yields a precession index-like term in the temperature. Despite its importance in the geologic record, the precession index e sin (Omega)S, where e is the Earth's orbital eccentricity and (Omega)S is the Sun's perigee in the geocentric frame, is not present in the insolation at the top of the atmosphere. Hence there is no one-for-one mapping of 23,000 and 19,000 year periodicities from the insolation to the paleoclimate record; a nonlinear climate model is needed to produce these long periods. A nonlinear energy balance climate model with radiative terms of form T n, where T is surface temperature and n less than 1, does produce e sin (omega)S terms in temperature; the e sin (omega)S terms are called Seversmith psychroterms. Without feedback mechanisms, the model achieves extreme values of 0.64 K at the maximum orbital eccentricity of 0.06, cooling one hemisphere while simultaneously warming the other; the hemisphere over which perihelion occurs is the cooler. In other words, the nonlinear energy balance model produces long-term cooling in the northern hemisphere when the Sun's perihelion is near northern summer solstice and long-term warming in the northern hemisphere when the aphelion is near northern summer solstice. (This behavior is similar to the inertialess gray body which radiates like T 4, but the amplitude is much lower for the energy balance model because of its thermal inertia.) This seemingly paradoxical behavior works against the standard Milankovitch model, which requires cool northern summers (Sun far from Earth in northern summer) to build up northern ice sheets, so that if the standard model is correct it must be more efficient than previously thought. Alternatively, the new mechanism could possibly be dominant and indicate southern hemisphere control of the northern ice sheets, wherein the southern oceans undergo a long-term cooling when the Sun is far from the Earth during northern summer. The cold water eventually flows north, cooling the northern hemisphere. This might explain why the northern oceans lag the southern ones when it comes to orbital forcing.

Rubincam, David↗

Weber Number Tests in the NASA Icing Research Tunnel

A study of the Weber Number effects on droplets in the NASA Icing Research Tunnel is described. The work focuses on examining the droplet Weber Number effects observed for droplets accelerated by air flow in the contraction section of the Icing Research Tunnel to the test section. These results will aid in Supercooled Large Drop facility design studies. Measurements acquired with the Phase Doppler Interferometer and High Speed Imaging Dual Range Flight Probes at a series of locations through the contraction are presented alongside a 1D numerical model developed during this study to aid interpretation of the experimental results. An estimate of the maximum Weber Number observed in the Icing Research Tunnel for varying drop sizes up to 1000 m is presented and provided for incorporation into future design studies. Finally, experimental results coupled with a numerical model indicate that breakup of drops up to 1000 m is not occurring in the NASA Icing Research Tunnel up to 129 m/s.

Icing Wind Tunnel↗

Weber Number Tests in the NASA Icing Research Tunnel

A study of the Weber Number effects on droplets in the NASA Icing Research Tunnel is described. The work focuses on examining the droplet Weber Number effects observed for droplets accelerated by air flow in the contraction section of the Icing Research Tunnel to the test section. These results will aid in Supercooled Large Drop facility design studies. Measurements acquired with the Phase Doppler Interferometer and High Speed Imaging Dual Range Flight Probes at a series of locations through the contraction are presented alongside a 1D numerical model developed during this study to aid interpretation of the experimental results. An estimate of the maximum Weber Number observed in the Icing Research Tunnel for varying drop sizes up to 1000 μm is presented and provided for incorporation into future design studies. Finally, experimental results coupled with a numerical model indicate that breakup of drops up to 1000 μm is not occurring in the NASA Icing Research Tunnel up to 129 m/s.

Icing Wind Tunnel↗

The physical mechanism of comet outbursts: An experiment

During a series of impact experiments into regolith-like powders at the NASA Ames Research Center Vertical Gun Facility in 1976, I observed and filmed a unique anomalous event that may illuminate outburst mechanisms in comets. During one test, a new batch of basalt powder (half the mass in particles less than 800 microns in diameter) retained some air pressure while the vacuum chamber was being evacuated. As a result, the projectile impacted into gas-charged regolith. Instead of ejecting the normal, relatively negligible amount of debris, the disturbance triggered a major eruption that lasted at least 18 seconds. The experimental results have been recently re-analyzed with reference to cometary phenomena. A series of frames from this eruption experiment are shown. The ejecta velocities of 150 to 300 cm/s would have been sufficient to drive debris into the coma of a comet nucleus smaller than a few kilometers diameter. The event suggests a mechanism for comet outbursts, discussed briefly by Hartmann et al.: the pore space in a layer of regolith, possibly with weak effective tensile strength, becomes gas charged as ice slowly sublimates. Once the effective tensile strength is exceeded by the gas pressure, the surface fails locally, triggering an eruption such as photographed here. This model is consistent with the emerging view of regolith materials on comets and is closest to the recent model of Rickman et al. The earlier models generally picture a more uniform flow of debris off the comet, not outbursts. Rickman et al. allow gas pressure to build until it matches the overburden pressure, followed by 'instantaneous blow-off'. They note that as soon as the mantle is found to be unstable, we consider it to be instantaneously swept away by the gas pressure. The main new points made here are that the experiment gives a more realistic view of the blow-off process after surface failure occurs, and the present model gives a recharge mechanism that can explain recurrent outbursts on comets such as P/Schwassmann-Wachmann 1 and 2060 Chiron. In fact, the resulting jets resemble distinct jet structures in high-resolution comet comae.

Hartmann, William K.↗

Effect of Geometric Fidelity on the Aerodynamics of a Swept Wing with Scalloped Ice Accretion

The effect of geometric fidelity on the aerodynamics of a swept wing with a “scallop” ice shape was studied. Three geometric fidelity versions of the ice shape were studied. The High Fidelity ice shape maintained all of the highly 3D features of the scallop shape. The 3D Smooth ice shape was smoothly lofted over the High Fidelity shape in order to eliminate all of the locally 3D feature. The third ice shape was the 3D Smooth ice shape with roughness grit attached to the surface. The two 3D Smooth versions of the ice shape exhibited a flowfield characterized by a leading edge separation bubble that rolled into spanwise running vortex. The surface pressure data showed classic leading edge separation bubbles that are observed on 2D airfoils with leading edge horn ice shapes. The High Fidelity ice shape exhibited a flowfield characterized by streamwise vortices that formed downstream of the ice shape. The streamwise vortices observed downstream of the High Fidelity ice shape appeared to reduce the size of the separation suction peak observed downstream of the 3D Smooth ice shapes. This reduced the lift for the High Fidelity ice shape, when compared to the 3D Smooth ice shapes. However, these streamwise vortices may have allowed the flow to remain attached longer and slightly increased the stall angle of attack, even though the maximum lift was lower.

Aircraft icing↗

Mathematical modeling of ice accretion on airfoils

The progress toward development of a computer model suitable for predicting icing behavior on airfoils over a wide range of environmental conditions and airfoils shapes is reported. The LEWICE program was formulated to solve a set of equations which describe the physical processes which occur during accretion of ice on an airfoil, including heat transfer in a time dependent mode, with the restriction that the flow must be describable by a two-dimensional flow code. Input data comprises the cloud liquid water content, mean droplet diameter, ambient air temperature, air velocity, and relative humidity. A potential flowfield around the airfoil is calculated, along with the droplet trajectories within the flowfield, followed by local values of water droplet collection efficiency at the impact points. Both glaze and rime ice conditions are reproduced, and comparisons with test results on icing of circular cylinders showed good agreement with the physical situation.

Macarthur, C. D.↗

SmaggIce 2.0: Additional Capabilities for Interactive Grid Generation of Iced Airfoils

The Surface Modeling and Grid Generation for Iced Airfoils (SmaggIce) software toolkit has been extended to allow interactive grid generation for multi-element iced airfoils. The essential phases of an icing effects study include geometry preparation, block creation and grid generation. SmaggIce Version 2.0 now includes these main capabilities for both single and multi-element airfoils, plus an improved flow solver interface and a variety of additional tools to enhance the efficiency and accuracy of icing effects studies. An overview of these features is given, especially the new multi-element blocking strategy using the multiple wakes method. Examples are given which illustrate the capabilities of SmaggIce for conducting an icing effects study for both single and multi-element airfoils.

Kreeger, Richard E.↗