Results for “ice flow model” — Page 15

NASA NTRS2013

Evidence for Water Ice 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.

Mini-RF radar↗

NASA NTRS2020

Application of Extended Messinger Models to Complex Geometries

Since, ice accretion can significantly degrade the performance and the stability of an airborne vehicle, it is imperative to be able to model it accurately. While ice accretion studies have been performed on airplane wings and helicopter blades in abundance, there are few that attempt to model the process on more complex geometries such as fuselages. This paper proposes a methodology that extends an existing in-house Extended Messinger solver to complex geometries by introducing the capability to work with unstructured grids and carry out spatial surface streamwise marching. For the work presented here commercial solvers such as STAR-CCM+ and ANSYS Fluent are used for the flow field and droplet dispersed phase computations. The ice accretion is carried out using an in-house icing solver called GT-ICE. The predictions by GT-ICE are compared to available experimental data, or to predictions by other solvers such as LEWICE and STAR-CCM+. Three different cases with varying levels of complexity are presented. The first case considered is a commercial transport airfoil, followed by a three-dimensional MS(1)-317 swept wing. Finally, ice accretion calculations performed on a Robin fuselage have been discussed. Good agreement with experimental data, where applicable, is observed. Differences between the ice accretion predictions by different solvers have been discussed.

Gupta, Avani↗

NASA NTRS1978

Documentation of a ground hydrology parameterization for use in the GISS atmospheric general circulation model

The moisture transport processes related to the earth's surface relevant to the ground circulation model GCM are presented. The GHM parametrizations considered are: (1) ground wetness and soil parameters; (2) precipitation; (3) evapotranspiration; (4) surface storage of snow and ice; and (5) runout. The computational aspects of the GHM using computer programs and flow charts are described.

Lin, J. D.↗

NASA NTRS2018

Total Temperature Measurements Using a Rearward Facing Probe in Supercooled Liquid Droplet and Ice Crystal Clouds

Engine Icing Performance loss: rollback, surge, flameout, and even internal engine damagePartial melting and refreeze of ice inside engine core (Mason et al., 2006). Ingestion of ice crystals and aggregates, mixed-phase droplets, or supercooled liquid dropletsNeed to better understand the conditions and properties that lead to engine icing.Simulation and analysis (physical and computational, and modeling)Test facilities (PSL, NRC, ...). Thermal and computational models and analysisProbesMultiple probes (aerothermal probes and ice cloud characterization probes and techniques). Total temperatureTraditional total temperature probes (vented forward facing)Heated total temperature probes (Goodyear). Rearward facing (developmental). Total temperature relevance. Thermal interaction between the icing cloud and air flow impinging particles contribute to kinetic heating effect (Gent et al., 2000). Measurement considerations Temperature sensor accuracy. Incomplete recovery of total temperature. Thermal surfaces (sources and sinks). Flow effects (viscous losses). Debris contamination, including icing and ice ingestion.

Agui, Juan H.↗

NASA NTRSDate not supplied

Comparisons of CFD Simulations of Icing Wind Tunnel Clouds with Experiments Conducted at the NASA Propulsions Systems Laboratory

This paper evaluates simulation predictions against experimental test data of icing clouds that were produced during 2018 ice crystal icing physics tests conducted at the NASA Propulsion Systems Laboratory icing wind tunnel. Aero-thermal and cloud parameters are set and known upstream at the tunnel inlet and spray system, but change as the cloud and air thermodynamically interact as the flowing masses reach the tunnel test section. Utilizing the ANSYS Fluent Discrete Phase Model function, 3D computational fluid dynamics (CFD) simulations were performed, capturing the thermodynamic interactions between the test parameters, and providing predictions of the aero thermal and cloud conditions at the tunnel test section. Simulation predictions were compared with test data that were measured at the tunnel exit plane. Evaluations focused on the cloud concentration (total water content), humidity content, and air temperature. CFD simulation predictions showed areas of agreement and disagreement. Simulations showed that cloud concentration profiles at the test section are strongly related to the initial spray nozzle pattern used at the tunnel inlet. Experimental data suggest that greater dispersion of the cloud occurred as the simulated cloud predicted areas of high and low cloud concentration compared to test data profiles. Simulations, however, captured the magnitude and location of the change in humidity content and the change in air temperature due to the presence of the cloud reasonably well, when compared to test data. This result would suggest that while the ANSYS simulation did not fully predict the spreading of the cloud as measured during experiment, it did capture evaporation and the molecular movements of air and vapor relatively well.

CFD↗

NASA NTRSDate not supplied

Comparisons of CFD Simulations of Icing Wind Tunnel Clouds with Experiments Conducted at the NASA Propulsions Systems Laboratory

This paper evaluates simulation predictions against experimental test data of icing clouds that were produced during 2018 ice crystal icing physics tests conducted at the NASA Propulsion Systems Laboratory icing wind tunnel. Aero-thermal and cloud parameters are set and known upstream at the tunnel inlet and spray system, but change as the cloud and air thermodynamically interact as the flowing masses reach the tunnel test section. Utilizing the ANSYS Fluent Discrete Phase Model function, 3D computational fluid dynamics (CFD) simulations were performed, capturing the thermodynamic interactions between the test parameters, and providing predictions of the aero thermal and cloud conditions at the tunnel test section. Simulation predictions were compared with test data that were measured at the tunnel exit plane. Evaluations focused on the cloud concentration (total water content), humidity content, and air temperature. CFD simulation predictions showed areas of agreement and disagreement. Simulations showed that cloud concentration profiles at the test section are strongly related to the initial spray nozzle pattern used at the tunnel inlet. Experimental data suggest that greater dispersion of the cloud occurred as the simulated cloud predicted areas of high and low cloud concentration compared to test data profiles. Simulations, however, captured the magnitude and location of the change in humidity content and the change in air temperature due to the presence of the cloud reasonably well, when compared to test data. This result would suggest that while the ANSYS simulation did not fully predict the spreading of the cloud as measured during experiment, it did capture evaporation and the molecular movements of air and vapor relatively well.

CFD↗

NASA NTRS1994

An Experimental Study of the Flowfield on a Semispan Rectangular Wing with a Simulated Glaze Ice Accretion

Wind tunnel experiments were conducted in order to study the effect of a simulated glaze ice accretion on the flowfield of a semispan, reflection-plane, rectangular wing at Re = 1.5 million and M = 0.12. A laser Doppler velocimeter was used to map the flowfield on the upper surface of the model in both the clean and iced configurations at alpha = 0, 4, and 8 degrees angle of attack. At low angles of attack, the massive separation bubble aft of the leading edge ice horn was found to behave in a manner similar to laminar separation bubbles. At alpha = 0 and 4 degrees, the locations of transition and reattachment, as deduced from momentum thickness distributions, were found to be in good agreement with transition and reattachment locations in laminar separation bubbles. These values at y/b = 0.470, the centerline measurement location, matched well with data obtained on a similar but two dimensional model. The measured velocity profiles on the iced wing compared reasonably with the predicted profiles from Navier-Stokes computations. The iced-induced separation bubble was also found to have features similar to the recirculating region aft of rearward-facing steps. At alpha = 0 degrees and 4 degrees, reverse flow magnitudes and turbulence intensity levels were typical of those found in the recirculating region aft of rearward-facing steps. The calculated separation streamline aft of the ice horn at alpha = 4 degrees, y/b = 0.470 coincided with the locus of the maximum Reynolds normal stress. The maximum Reynolds normal stress peaked at two locations along the separation streamline. The location of the first peak-value coincided with the transition location, as deduced from the momentum thickness distributions. The location of the second peak was just upstream of reattachment, in good agreement with measurements of flows over similar obstacles. The intermittency factor in the vicinity of reattachment at alpha = 4 degrees, y/b = 0.470, revealed the time-dependent nature of the reattachment process. The size and extent of the separation bubble were found to be a function of angle of attack and the spanwise location. Three dimensional effects were found to be strongest at alpha = 8 degrees. The calculated separation and stagnation streamlines were found to vary little with spanwise location at alpha = 0 degrees. The calculated separation streamlines at alpha = 4 degrees revealed that the bubble was largest near the centerline measurement plane, whereas the tip-induced vortex flow and the model root-tunnel wall boundary-layer interaction reduced the size of the bubble. These effects were found to be most dramatic at alpha = 8 degrees.

Khodadoust, Abdollah↗

NASA NTRS1991

Three-dimensional calculation of the mixing of radial jets from slanted slots with a reactive cylindrical crossflow

A numerical model that calculates the 3D chemically reacting flowfield in an experimental low emission combustor is described. The ICEd-ALE finite volume computational methodology is employed in this study. Radial jets issuing from slanted slots interact with a cylindrical axially flowing mainstream to produce a bulk swirl velocity downstream of the slot region. The swirl pattern at a given axial station is composed of a clockwise rotating region near the wall and a counterclockwise rotating region extending from the combustor centerline. The jet radial penetration and downstream swirl velocity axial development are shown to depend principally on the jet-to-mainstream momentum flux ratio.

Winowich, N. S.↗

NASA NTRS1989

Thermal histories of the samples of two KOSI comet nucleus simulation experiments

Temperatures recorded during two KOSI comet nucleus simulation experiments strongly suggest that heat transport by vapor flow into the interior of the sample is very important. Two comet nucleus simulation experiments have been done by the KOSI team in a big space simulator. The thermal evolution of the sample during insolation and the results of simplified thermal evolution calculations are discussed. The observed thermal histories cannot be explained by a simple model with heat transferred by heat conduction at a constant conductivity, so a coupled heat and mass transfer problem was considered. The porous ice matrix was assumed to have a constant thermal conductivity and to be in thermal equilibrium with vapor in the pores, the internal pressure being the vapor pressure. The vapor was modelled as an ideal gas because, at the temperatures relevant to the problem, the mean free path length of the vapor molecules is large in comparison with the pore dimensions. The heat capacity at constant volume per unit mass of the two phase mixture was also assumed constant. The vapor was allowed to flow and transfer heat in response to an internal pressure gradient.

Spohn, T.↗

NASA NTRS1986

The Comet Giacobini-Zinner magnetotail: Axial stresses and inferred near-nucleus properties

Utilizing the electron and magnetic field data from the ICE tail traversal of comet Giacobini-Zinner along with the MHD equations, a steady state, stress balance model of the cometary magnetotail was developed, and used to infer important but unmeasured ion properties within the magnetotail at ICE and upstream at the average point along each streamline where cometary ions are picked-up. The derived tailward ion flow speed at ICE is quite constant at approx. -20 to -30 km/sec across the entire tail. The flow velocity, ion temperature, density, and ion source rates upstream from the lobes (current sheet) at the average pick-up locations are approx. -75 km/sec (approx. -12), approx. 4 million K (approx. 100,000), approx. 20 cc (approx. 400), and approx. 15 cu cm/sec. Gradients in the plasma properties between the two regions are quite strong. Implications of inferred plasma properties for the near-nucleus region and for cometary magnetotail formation are examined.

Mccomas, D. J.↗

NASA NTRS2017

Mathematical Modeling of Electrodynamics Near the Surface of Earth and Planetary Water Worlds

An interesting feature of planetary bodies with hydrospheres is the presence of an electrically conducting shell near the global surface. This conducting shell may typically lie between relatively insulating rock, ice, or atmosphere, creating a strong constraint on the flow of large-scale electric currents. All or parts of the shell may be in fluid motion relative to main components of the rotating planetary magnetic field (as well as the magnetic fields due to external bodies), creating motionally-induced electric currents that would not otherwise be present. As such, one may expect distinguishing features in the types of electrodynamic processes that occur, as well as an opportunity for imposing specialized mathematical methods that efficiently address this class of application. The purpose of this paper is to present and discuss such specialized methods. Specifically, thin-shell approximations for both the electrodynamics and fluid dynamics are combined to derive simplified mathematical formulations describing the behavior of these electric currents as well as their associated electric and magnetic fields. These simplified formulae allow analytical solutions featuring distinct aspects of the thin-shell electrodynamics in idealized cases. A highly efficient numerical method is also presented that is useful for calculations under inhomogeneous parameter distributions. Finally, the advantages as well as limitations in using this mathematical approach are evaluated. This evaluation is presented primarily for the generic case of bodies with water worlds or other thin spherical conducting shells. More specific discussion is given for the case of Earth, but also Europa and other satellites with suspected oceans.

Tyler, Robert H.↗

NASA NTRS1985

Comprehensive modeling of a liquid rocket combustion chamber

An analytical model for the simulation of detailed three-phase combustion flows inside a liquid rocket combustion chamber is presented. The three phases involved are: a multispecies gaseous phase, an incompressible liquid phase, and a particulate droplet phase. The gas and liquid phases are continuum described in an Eulerian fashion. A two-phase solution capability for these continuum media is obtained through a marriage of the Implicit Continuous Eulerian (ICE) technique and the fractional Volume of Fluid (VOF) free surface description method. On the other hand, the particulate phase is given a discrete treatment and described in a Lagrangian fashion. All three phases are hence treated rigorously. Semi-empirical physical models are used to describe all interphase coupling terms as well as the chemistry among gaseous components. Sample calculations using the model are given. The results show promising application to truly comprehensive modeling of complex liquid-fueled engine systems.

Liang, P.-Y.↗

NASA NTRSDate not supplied

Utilization of Streamtubes to Analyze the Physical Interaction of a Dispersed Cloud with the CRM65 Hybrid Midspan Model

Recently there have been numerous efforts to identify the relevant icing physics related to the formation of complex three dimensional features, sometimes referred to as ’scallops’, on swept wings. However, much of the physics is still not well understood. This paper computationally investigates the interaction of the icing cloud with the 65 percent Common Research Model (CRM65). Both the interaction with an uniced model, and the interaction with a representative simulated three dimensional ice accretion are analyzed. Preliminary results suggest that the liquid water content increases near the aerodynamic body. For small droplets on the uniced geometry, the particle velocity vector becomes nearly parallel with the aerodynamic model resulting in a low value of collection efficiency. When a three dimensional feature with a length scale much smaller than the leading edge of the airfoil is introduced into the flow, the impingement of small particles can become significantly more perpendicular to the particle velocity vector in the region of these features. Since the liquid water content near the body has increased due to the interaction with the larger scale features of the aerodynamic model, i.e. the leading edge of the swept wing, extremely high collection efficiency is observed. These results suggest that three dimensional features are likely a significant physical driver that should be modeled in some capacity when simulating the impingement of a cloud on an aerodynamic model.

Icing↗

NASA NTRSDate not supplied

Utilization of Streamtubes to Analyze the Physical Interaction of a Dispersed Cloud with the CRM65 Hybrid Midspan Model

Recently there have been numerous efforts to identify the relevant icing physics related to the formation of complex three dimensional features, sometimes referred to as ’scallops’, on swept wings. However, much of the physics is still not well understood. This paper computationally investigates the interaction of the icing cloud with the 65 percent Common Research Model (CRM65). Both the interaction with an uniced model, and the interaction with a representative simulated three dimensional ice accretion are analyzed. Preliminary results suggest that the liquid water content increases near the aerodynamic body. For small droplets on the uniced geometry, the particle velocity vector becomes nearly parallel with the aerodynamic model resulting in a low value of collection efficiency. When a three dimensional feature with a length scale much smaller than the leading edge of the airfoil is introduced into the flow, the impingement of small particles can become significantly more perpendicular to the particle velocity vector in the region of these features. Since the liquid water content near the body has increased due to the interaction with the larger scale features of the aerodynamic model, i.e. the leading edge of the swept wing, extremely high collection efficiency is observed. These results suggest that three dimensional features are likely a significant physical driver that should be modeled in some capacity when simulating the impingement of a cloud on an aerodynamic model.

Icing↗

NASA NTRS1991

Surface reaction and pore diffusion in flow-tube reactors

The interaction of gas diffusion with surface reaction in porous solids is discussed and applied specifically to heterogeneous rate measurements in flow-tube reactors. External diffusion to the outer surface of a reactive solid, internal diffusion within the pores, surface reaction, and laminar flow are considered. A procedure is developed to correct observed surface rate constants for the interaction of these processes. Measured surface areas and bulk densities are used to construct a semiempirical model for porous diffusion in vapor-formed HNO3-H2O ices which are used to simulate polar stratospheric cloud surfaces. The model is tested experimentally by varying the thickness of these ices from about 15 to 120 microns. The results are consistent with the model predictions and show that the HNO3-H2O ices used are highly porous, and the internal surface must be considered in calculating kinetic parameters from observed loss rates. The best fit of the data yields a tortuosity factor of 3.3 +/-1.1 for the ice substrates.

Keyser, Leon F.↗

NASA NTRS2009

Old Faithful Model for Radiolytic Gas-Driven Cryovolcanism at Enceladus

A new model is presented on how chemically driven cryovolcanism might contribute to episodic outgassing at the icy moon Enceladus and potentially elsewhere including Europa and Kuiper Belt Objects. Exposed water ices can become oxidized from radiolytic chemical alteration of near-surface water ice by space environment irradiation. In contact with primordially abundant reductants such as NH3, CH4, and other hydrocarbons, the product oxidants can react exothermically to produce volatile gases driving cryovolcanism via gas-piston forces on any subsurface liquid reservoirs. Radiolytic oxidants such as H2O2 and O2 can continuously accumulate deep in icy regoliths and be conveyed by rheological flows to subsurface chemical reaction zones over million-year time scales indicated by cratering ages for active regions of Enceladus and Europa. Surface blanketing with cryovolcanic plume ejecta would further accelerate regolith burial of radiolytic oxidants. Episodic heating from transient gravitational tides, radioisotope decay, impacts, or other geologic events might occasionally accelerate chemical reaction rates and ignite the exothermic release of cumulative radiolytic oxidant energy. The time history for the suggested "Old Faithful" model of radiolytic gas-driven cryovolcanism at Enceladus and elsewhere therefore consists of long periods of chemical energy accumulation punctuated by much briefer episodes of cryovolcanic activity. The most probable sequence for detection of activity in the current epoch is a long evolutionary phase of slow but continuous oxidant accumulation over billions of years followed by continuous but variable high activity over the past 10(exp 7)-10(exp 8) years. Detectable cryovolcanic activity could then later decline due to near-total oxidation of the rheologically accessible ice crust and depletion the accessible reductant abundances, as may have already occurred for Europa in the more intense radiation environment of Jupiter's magnetosphere. Astrobiological potential of Enceladus could correspondingly be higher than at Europa due to a less extreme state of oxidation and greater residual abundance of organics.

Cooper, John F.↗

NASA NTRS1991

Preparation of the Ice Certification of the Dornier 328 Regional Airliner By Numerical Simulation and By Ground Test

The Dornier 328, a new regional airliner, is to be qualified to FAR25/JAR25 requirements for operation into known icing conditions. All lifting surfaces are pneumatically deiced. Impingement limits were determined using a validated computational fluid dynamic (CFD) code. Wind tunnel tests with a model of the whole aircraft and the isolated empennage equipped with artificial ice shapes showed a degradation of handling characteristics due to ice accretion. Extensive two and three dimensional CFD calculations evaluated ice accretion at the deflection elevator horn. Icing tunnel tests were performed with a model of the horizontal tail with icing similitude fulfilled. The Dornier Do328 powerplant air induction system has successfully completed a program of icing tests in support of aircraft certification for operation into known icing conditions. Potential flight icing conditions were extensively analyzed using a CFD potential flow code and 3-D particle tracking routine to identify critical flight icing conditions. The tunnel test program verified the analytically predicted critical ice accretion surfaces of the air induction system, the adequacy of the ice protection provisions and demonstrated compliance with the applicable (JAR25) regulatory material.

D Welte↗

NASA NTRS2004

Mangala Valles, Mars: Investigations of the Source of Flood Water and Early Stages of Flooding

Mangala Valles, an approx. 900 km long north-south trending outflow channel located southwest of the Tharsis rise, extends northward from one of the Memnonia Fossae graben across the southern highlands, terminating at the dichotomy boundary. Previous Viking-based analyses suggest that the water that carved the channel was expelled from the graben, possibly during two distinct flood events, one in the Late Hesperian and one in the Latest Hesperian/Early Amazonian. The mechanism by which the water was transported to the graben, and ultimately to the surface remained ambiguous, although two general scenarios were proposed: melting of near surface ground ice via nearby Tharsis lava flows, and tapping of a near surface aquifer via faulting associated with the graben. Here we use MOLA altimetric data and MOC and THEMIS images to reexamine Mangala Valles and the surrounding region. Further, we develop a new model for the production and transport of the floodwater.

Ghatan, Gil J.↗