Mesoscale Simulations Investigating the Dynamic Strength of Partially Transitioned Water-Ice Mixtures Under Pressure-Shear Impact Loading
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
No abstract provided
In the past year, we have successfully developed the techniques necessary to conduct impact experiments on ice at very low temperatures. We employ the method of embedding gauges within a target to measure the shock wave and material properties. This means that our data are not model dependent; we directly measure the essential parameters needed for numerical simulations of impact cratering. Since then we have developed a new method for temperature control of icy targets that ensures temperature equilibrium throughout a porous target. Graduate student, Sarah Stewart-Mukhopadhyay, is leading the work on ices and porous materials as the main thrust of her thesis research. Our previous work has focused on icy materials with no porosity, and we propose to extend our research to include porous ice and porous ice-silicate mixtures. There is little shockwave data for porous ice, and none of the data was acquired under conditions applicable to the outer solar system. The solid ice Hugoniot is only defined for initial temperatures above -20 C. Our program uniquely measures the properties of ice at temperatures directly applicable to the solar system. Previous experiments were conducted at ambient temperatures soon after removing the target from a cold environment, usually just below freezing, or in a room just below freezing. Since ice has an extremely complicated phase diagram, it is important to conduct experiments at lower temperatures to determine the true outcome of impacts in the outer solar system. This research is complementary to other programs on icy materials. Our work focuses on the inherent material properties by measuring the shock wave directly; this complements the macroscopic observations and immediately provides the parameters necessary to extend this research to the gravity regime. Our numerical simulations of impacts in porous ice under very low gravity conditions, such as found on comets, show that the final crater size and shape is very dependent on the dynamic strength of the material.
After the loss of Columbia, there was great concern in the Space Shuttle program for the impact of debris against the leading edges of the Orbiter wings. It was quickly recognized that, in addition to impacts by foam, ice that formed on the liquid-oxygen bellows running down the outside of the External Tank could break free during launch and hit this sensitive area. A Center Director s Discretionary Fund (CDDF) project would concentrate on novel ideas that were potentially applicable. The most successful of the new concepts for ice mitigation involved shape memory alloy materials. These materials can be bent into a given shape and, when heated, will return to their original shape.
Both Saturn's rings and planetesimal disks are made up of particles in Keplerian orbits. Inelastic collisions between these particles regulate their dynamical evolution and possible aggregation. We present an experiment to simulate glancing collisions in Saturn's rings and in planetesimal disks and thus measure contributions to the energy loss for both normal and tangential velocity components. In this experiment, a spherical iceball mounted on a long-period, two dimensional pendulum is made to impact a flat ice surface in a low-temperature environment. This paper describes the experimental apparatus in detail and presents results for smooth unfrosted surfaces. The energy loss for tangential motion is suprisingly low, indicating that very little friction is present at low impact speeds for relatively smooth ice surfaces and temperatures near 100 K. We have also investigated room-temperature collisions of a rubber ball on a rough surface to understand the energy loss in situations where the tangential friction force is not small. In this analogous case, the energy loss is maximum for impact angles in the range 45 deg-60 deg.
A flow and ice particle trajectory analysis was performed for the booster of the Honeywell ALF502 engine. The analysis focused on two closely related conditions one of which produced an icing event and another which did not during testing of the ALF502 engine in the Propulsion Systems Lab (PSL) at NASA Glenn Research Center. The flow analysis was generated using the NASA Glenn GlennHT flow solver and the particle analysis was generated using the NASA Glenn LEWICE3D v3.63 ice accretion software. The inflow conditions for the two conditions were similar with the main differences being that the condition that produced the icing event was 6.8 K colder than the non-icing event case and the inflow ice water content (IWC) for the non-icing event case was 50% less than for the icing event case. The particle analysis, which considered sublimation, evaporation and phase change, was generated for a 5 micron ice particle with a sticky impact model and for a 24 micron median volume diameter (MVD), 7 bin ice particle distribution with a supercooled large droplet (SLD) splash model used to simulate ice particle breakup. The particle analysis did not consider the effect of the runback and re-impingement of water resulting from the heated spinner and anti-icing system. The results from the analysis showed that the amount of impingement for the components were similar for the same particle size and impact model for the icing and non-icing event conditions. This was attributed to the similar aerodynamic conditions in the booster for the two cases. The particle temperature and melt fraction were higher at the same location and particle size for the non-icing event than for the icing event case due to the higher incoming inflow temperature for the non-event case. The 5 micron ice particle case produced higher impact temperatures and higher melt fractions on the components downstream of the fan than the 24 micron MVD case because the average particle size generated by the particle breakup was larger than 5 microns which yielded less warming and melting. The analysis also showed that the melt fraction and wet bulb temperature icing criterion developed during tests in the Research Altitude Test Facility (RATFac) at the National Research Council (NRC) of Canada were useful in predicting icing events in the ALF502 engine. The development of an ice particle impact model which includes the effects of particle breakup, phase change, and surface state is necessary to further improve the prediction of ice particle transport with phase change through turbomachinery.
For the period 1978-1982, this paper examines the nature of the sea ice which forms over the Kuril Basin of the Okhotsk Sea and describes the impact of this ice on the regional oceanography. The oceanographic behavior during the heavy ice season associated with the cold 1979 winter is compared with the behavior during the lighter ice years of 1980 and 1982. Examination of the oceanography in the Okhotsk and the adjacent Pacific shows that the early summer water column structure depends on the heat loss from the Okhotsk during the preceding ice season, the total amount of Okhotsk ice formation, and, specifically, the amount of the ice formation in the Kuril Basin. Following the 1979 ice season, the upper 200-300 m of the Kuril Basin waters were cooler, less saline, and richer in oxygen than for the other years. This modification appears to be a process local to the Kuril Basin, driven by eddy-induced mixing, local cooling, and ice melting.
We investigate the chemical (hydroxyl and ozone) and dynamical response to changing from present day to pre-industrial conditions in the Goddard Institute for Space Studies General Circulation Model (GISS GMC). We identify three main improvements not included by many other works. Firstly, our model includes interactive cloud calculations. Secondly we reduce sulfate aerosol which impacts NOx partitioning hence Ox distributions. Thirdly we reduce sea surface temperatures and increase ocean ice coverage which impact water vapor and ground albedo respectively. Changing the ocean data (hence water vapor and ozone) produces a potentially important feedback between the Hadley circulation and convective cloud cover. Our present day run (run 1, control run) global mean OH value was 9.8 x 10(exp 5) molecules/cc. For our best estimate of pre-industrial conditions run (run 2) which featured modified chemical emissions, sulfate aerosol and sea surface temperatures/ocean ice, this value changed to 10.2 x 10(exp 5) molecules/cc. Reducing only the chemical emissions to pre-industrial levels in run 1 (run 3) resulted in this value increasing to 10.6 x 10(exp 5) molecules/cc. Reducing the sulfate in run 3 to pre-industrial levels (run 4) resulted in a small increase in global mean OH (10.7 x 10(exp 5) molecules/cc). Changing the ocean data in run 4 to pre-industrial levels (run 5) led to a reduction in this value to 10.3 x 10(exp 5) molecules/cc. Mean tropospheric ozone burdens were 262, 181, 180, 180, and 182 Tg for runs 1-5 respectively.
We conduct sensitivity experiments using a general circulation model that has an explicit water source tagging capability forced by prescribed composites of pre-industrial sea-ice concentrations (SICs) and corresponding sea surface temperatures (SSTs) to understand the impact of sea-ice anomalies on regional evaporation, moisture transport and source–receptor relationships for Antarctic precipitation in the absence of anthropogenic forcing. Surface sensible heat fluxes, evaporation and column-integrated water vapor are larger over Southern Ocean (SO) areas with lower SICs. Changes in Antarctic precipitation and its source attribution with SICs have a strong spatial variability. Among the tagged source regions, the Southern Ocean (south of 50∘ S) contributes the most (40 %) to the Antarctic total precipitation, followed by more northerly ocean basins, most notably the South Pacific Ocean (27%), southern Indian Ocean (16 %) and South Atlantic Ocean (11 %). Comparing two experiments prescribed with high and low pre-industrial SICs, respectively, the annual mean Antarctic precipitation is about 150 Gt yr−1 (or 6 %) more in the lower SIC case than in the higher SIC case. This difference is larger than the model-simulated interannual variability in Antarctic precipitation (99 Gt yr−1). The contrast in contribution from the Southern Ocean, 102 Gt yr−1, is even more significant compared to the interannual variability of 35 Gt yr−1 in Antarctic precipitation that originates from the Southern Ocean. The horizontal transport pathways from individual vapor source regions to Antarctica are largely determined by large-scale atmospheric circulation patterns. Vapor from lower-latitude source regions takes elevated pathways to Antarctica. In contrast, vapor from the Southern Ocean moves southward within the lower troposphere to the Antarctic continent along moist isentropes that are largely shaped by local ambient conditions and coastal topography. This study also highlights the importance of atmospheric dynamics in affecting the thermodynamic impact of sea-ice anomalies associated with natural variability on Antarctic precipitation. Our analyses of the seasonal contrast in changes of basin-scale evaporation, moisture flux and precipitation suggest that the impact of SIC anomalies on regional Antarctic precipitation depends on dynamic changes that arise from SIC–SST perturbations along with internal variability. The latter appears to have a more significant effect on the moisture transport in austral winter than in summer.
Combinations of sea ice freeboard and snow depth measurements from satellite data have the potential to provide a means to derive global sea ice thickness values. However, large differences in spatial coverage and resolution between the measurements lead to uncertainties when combining the data. High resolution airborne laser altimeter retrievals of snow-ice freeboard and passive microwave retrievals of snow depth taken in March 2006 provide insight into the spatial variability of these quantities as well as optimal methods for combining high resolution satellite altimeter measurements with low resolution snow depth data. The aircraft measurements show a relationship between freeboard and snow depth for thin ice allowing the development of a method for estimating sea ice thickness from satellite laser altimetry data at their full spatial resolution. This method is used to estimate snow and ice thicknesses for the Arctic basin through the combination of freeboard data from ICESat, snow depth data over first-year ice from AMSR-E, and snow depth over multiyear ice from climatological data. Due to the non-linear dependence of heat flux on ice thickness, the impact on heat flux calculations when maintaining the full resolution of the ICESat data for ice thickness estimates is explored for typical winter conditions. Calculations of the basin-wide mean heat flux and ice growth rate using snow and ice thickness values at the 70 m spatial resolution of ICESat are found to be approximately one-third higher than those calculated from 25 km mean ice thickness values.
The study of sea ice using airborne remote sensing platforms provides unique capabilities to measure a wide variety of sea ice properties. These measurements are useful for a variety of topics including model evaluation and improvement, assessment of satellite retrievals, and incorporation into climate data records for analysis of interannual variability and long-term trends in sea ice properties. In this paper we describe methods for the retrieval of sea ice thickness, freeboard, and snow depth using data from a multisensor suite of instruments on NASA's Operation IceBridge airborne campaign. We assess the consistency of the results through comparison with independent data sets that demonstrate that the IceBridge products are capable of providing a reliable record of snow depth and sea ice thickness. We explore the impact of inter-campaign instrument changes and associated algorithm adaptations as well as the applicability of the adapted algorithms to the ongoing IceBridge mission. The uncertainties associated with the retrieval methods are determined and placed in the context of their impact on the retrieved sea ice thickness. Lastly, we present results for the 2009 and 2010 IceBridge campaigns, which are currently available in product form via the National Snow and Ice Data Center
Computational methods for ice accretion can simulate the impact of water drops and ice crystals on an aircraft surface in a Lagrangian reference frame or in the Eulerian reference frame. In the Eulerian reference frame, particles are considered a continuous fluid while in the Lagrangian frame individual particle trajectories are calculated. Methods that use the Eulerian reference frame are typically easier to develop as established modules used for continuum mechanics can be leveraged. The Eulerian systems can also be faster since the user does not have to simulate millions of particles in order to achieve good results. It is imperative therefore that a Lagrangian method optimize the release points of trajectories such that accurate solutions can be obtained while minimizing as much as possible the number of trajectories computed. This paper will present a methodology for this refinement process and demonstrate its effectiveness on sample three dimensional test cases.
Computational methods for ice accretion can simulate the impact of water drops and ice crystals on an aircraft surface in a Lagrangian reference frame or in the Eulerian reference frame. In the Eulerian reference frame, particles are considered a continuous fluid while in the Lagrangian frame individual particle trajectories are calculated. Methods that use the Eulerian reference frame are typically easier to develop as established modules used for continuum mechanics can be leveraged. The Eulerian systems can also be faster since the user does not have to simulate millions of particles in order to achieve good results. It is imperative therefore that a Lagrangian method optimize the release points of trajectories such that accurate solutions can be obtained while minimizing as much as possible the number of trajectories computed. This paper will present a methodology for this refinement process and demonstrate its effectiveness on sample three dimensional test cases.
Mixed-phase clouds at high latitudes contribute to the uncertainty in predicting cloud feedbacks and climate sensitivity, mainly due to the complexity of microphysical processes that influence the partitioning between the supercooled liquid and ice phases, and hence, cloud radiative effects on regional scales. Particularly in Arctic mixed-phase clouds, the activation of ice-nucleating particles (INPs) from various aerosol populations remains a leading source of uncertainty. We developed an aerosol-cloud one-dimensional (AC-1D) model, which provides a novel framework to prognostically treat INP and ice crystal budgets while explicitly accounting for polydisperse and multicomponent aerosol that activate INPs following different freezing parameterizations. The AC-1D model is informed by large-eddy simulations to probe the impact of INP representation on predicted ice crystal number concentrations (N i ) and ice crystal budgets in mixed-phase Arctic stratus. We apply three immersion freezing (IMF) parameterizations, two time-independent (singular) and one time-dependent (classical nucleation theory), to predict the evolution of the INP reservoir and resulting ice crystal budget from polydisperse mineral dust, organic (humic-like substances), and sea spray aerosol particle size distributions. Our analysis focuses on how variations in aerosol number concentration and cloud system parameters such as cloud cooling rate, cloud-top entrainment rate, and ice crystal fall speed influence the INP reservoir and ice crystal budgets. Furthermore, this study investigates the competitive ice nucleation dynamics in mixed aerosol environments and provides a process-level quantification of the INP budget terms, which directly controls ice crystal budgets. For all studied case scenarios, the aerosol types and associated particle size distributions significantly impact INP and N i , and the choice between a time-dependent and a singular freezing description yields orders-of-magnitude differences in the predicted INP and N i over the 10 h simulation time, reflecting typical cloud lifetimes. Our results show that the influence of cloud cooling, INP entrainment, and sedimentation varies significantly depending on the chosen freezing parameterization. These findings underscore the critical need for robust IMF parameterizations and precise cloud system observations to enhance the accuracy of models in predicting mixed-phase cloud structure and evolution.
In the present study a modified split Hopkinson pressure bar (SHPB) was employed to study the effect of strain rate on the dynamic material response of ice. Disk-shaped ice specimens with flat, parallel end faces were either provided by Dartmouth College (Hanover, NH) or grown at Case Western Reserve University (Cleveland, OH). The SHPB was adapted to perform tests at high strain rates in the range 60 to 1400/s at test temperatures of -10 and -30 C. Experimental results showed that the strength of ice increases with increasing strain rates and this occurs over a change in strain rate of five orders of magnitude. Under these strain rate conditions the ice microstructure has a slight influence on the strength, but it is much less than the influence it has under quasi-static loading conditions. End constraint and frictional effects do not influence the compression tests like they do at slower strain rates, and therefore the diameter/thickness ratio of the samples is not as critical. The strength of ice at high strain rates was found to increase with decreasing test temperatures. Ice has been identified as a potential source of debris to impact the shuttle; data presented in this report can be used to validate and/or develop material models for ice impact analyses for shuttle Return to Flight efforts.
The resistance of current-design Pratt and Whitney Aircraft low aspect ratio advanced fiber reinforced epoxy matrix composite fan blades to foreign object damage (FOD) at STOL operating conditions was investigated. Five graphite/epoxy and five boron/epoxy wide chord fan blades with nickel plated stainless steel leading edge sheath protection were fabricated and impact tested. The fan blades were individually tested in a vacuum whirlpit under FOD environments. The FOD environments were typical of those encountered in service operations. The impact objects were ice balls, gravel, stralings and gelatin simulated birds. Results of the damage sustained from each FOD impact are presented for both the graphite boron reinforced blades. Tests showed that the present design composite fan blades, with wrap around leading edge protection have inadequate FOD impact resistance at 244 m/sec (800 ft/sec) tip speed, a possible STOL operating condition.
Results are presented on the following:(A) Ab initio quantum chemical studies of reactions in astrophysical ices.Theoretical electronic structure calculations were used to investigate reactions between formaldehyde (H2CO) and both hydrogen cyanide (HCN) and isocyanide (HNC) in search of other favorable reactions such as ammonia-formaldehyde addition, which was found in a recent theoretical study to be strongly enhanced when it occurs within cold ices.The present study examines further reactions between this product and H2CO in ices.(B) Heterogeneous hydrogenation of CO and H2CO on icy grain mantles.Formaldehyde (H2CO) and methanol (CH30H) are thought to be produced in the interstellar medium by the successive hydrogenation of carbon monoxide (CO) on grain surfaces. In the gas phase, the steps in which H adds to CO and H2CO possess modest barriers and are too inefficient to account for the observed abundances. Recent laboratory work has confirmed that formaldehyde and methanol are formed when H atoms are deposited on CO ice at 12 K. The present study employed ab initio quantum chemical calculations to investigate the impact of water ice on the sequential hydrogenation of CO.(C) Ice-bound condensed-phase reactions involving formic acid (HCOOH), methylenimine (CH2NH), hydrogen cyanide (HCN), hydrogen isocyanide (HNC), and ammonia ( 3) were investigated in order to characterize possible pathways to larger organic species that are efficient at the cold temperatures prevalent in cometary nuclei and the interstellar medium. (D) Pathways to glycine and other amino acids in ultraviolet-irradiated ices determined via quantum chemical modeling.(E) Photoionization in ultraviolet processing of astrophysical ice analogs at cryogenic temperatures.