SEARCH · Search NASA

Results for “ice flow equations”

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.

SourceResearch connectionFrom yearThrough year

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

87 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS2000

Convective Melting of Particles in Flow Under Microgravity Conditions

Study of melting of dispersed or packed solid particles in a fluid under gravity and microgravity conditions provides benchmark information for many engineering applications such as material processing, environmental assessment and protection and space fire protection. During such processes, packed or dispersed solid particles are interacting with fluid flow at above-melting temperatures. By unmasking the buoyancy effects in coupled flow, phase change and heat transport phenomena, a better understanding of melting rate in non-thermal equilibrium, convective conditions can be studied. A series of flight experiments were conducted onboard the NASA KC-135 Microgravity Research Airplane. The Particle Melting in Plow (PMF) module was designed to allow flow through the initially packed ice particles at controlled temperature and velocity. To achieve this, a close-loop flow system was designed. Video images were taken to record the visualization of the melting process, from which a time variation of packed particle thickness distribution at different times can be obtained by the image analysis method. The fluid temperature distribution within the melting zone is measured by thermocouples. An infrared camera was mounted from the top of the test section to record the ice-water thermal images at a given location. The results from thermal images yield local temperature variation between melting solid and liquid and local Stephan number. Typical results for a number of cases are presented. The mathematical model, describing mass, energy, and momentum balance equations for the liquid and solid phases, is presented. It is found that melting rate is influenced mainly by the ratio of Reynolds number (based on the initial particle diameter) to the Froud number, and Stephan number. At the absence of gravity, Froud number approaches zero, Reynolds number and Stephan number become dominant factors governing the melting rate. The numerically determined results are compared with the experimental ones. It is found that the discrepancy between the predicted and measured melting rate is largely due to the inaccuracy in the constitutive equations for effective thermal physical properties, such as effective thermal conductivity and diffusivity, and transport properties, such as particle interaction coefficient, and local heat transfer coefficient of particles.

Stewart, C.↗

NASA NTRS1991

Thermal models for icy satellites

The solid state greenhouse models were constrained for water ice surfaces using thermal eclipse data. As a result, it was determined that at depths of 3 cm or more Europa can have about a 30 K temperature enhancement over the diurnal mean surface temperature of about 97 K at the equator. It was also determined that the unusual shape of Europa's thermal flux curves during eclipses can be explained by the solid state greenhouse effect. The solid state greenhouse models were applied to the problem of the origin and energy sources of Triton's geyser-like plumes. Two classes of models are examined which may explain the driving energy for Triton's plumes. Triton's internal heatflow is recognized as a major contributor to its global energy balance. Internal heat flow from radioactive decay in Triton's interior along with absorbed thermal energy from Neptune total 5 to 20 pct. of the insolation absorbed by Triton, thus comprising a significant fraction of Triton's surface energy balance.

Brown, Robert H.↗

NASA NTRS1994

NASA Goddard Earth Sciences Graduate Student Program

The evolution of synoptic-scale dynamics associated with a middle and upper tropospheric cloud event that occurred on 26 November 1991 is examined. The case under consideration occurred during the FIRE CIRRUS-II Intensive Field Observing Period held in Coffeyville, KS during Nov. and Dec., 1991. Using data from the wind profiler demonstration network and a temporally and spatially augmented radiosonde array, emphasis is given to explaining the evolution of the kinematically-derived ageostrophic vertical circulations and correlating the circulation with the forcing of an extensively sampled cloud field. This is facilitated by decomposing the horizontal divergence into its component parts through a natural coordinate representation of the flow. Ageostrophic vertical circulations are inferred and compared to the circulation forcing arising from geostrophic confluence and shearing deformation derived from the Sawyer-Eliassen Equation. It is found that a thermodynamically indirect vertical circulation existed in association with a jet streak exit region. The circulation was displaced to the cyclonic side of the jet axis due to the orientation of the jet exit between a deepening diffluent trough and building ridge. The cloud line formed in the ascending branch of the vertical circulation with the most concentrated cloud development occurring in conjunction with the maximum large-scale vertical motion. The relationship between the large scale dynamics and the parameterization of middle and upper tropospheric clouds in large-scale models is discussed and an example of ice water contents derived from a parameterization forced by the diagnosed vertical motions and observed water vapor contents is presented.

Ackerman, Thomas P.↗

NASA NTRS1994

Extended field observations of cirrus clouds using a ground-based cloud observing system

Ackerman, Thomas P.↗

NASA NTRS1994

High-lift flow-physics flight experiments on a subsonic civil transport aircraft (B737-100)

As part of the subsonic transport high-lift program, flight experiments are being conducted using NASA Langley's B737-100 to measure the flow characteristics of the multi-element high-lift system at full-scale high-Reynolds-number conditions. The instrumentation consists of hot-film anemometers to measure boundary-layer states, an infra-red camera to detect transition from laminar to turbulent flow, Preston tubes to measure wall shear stress, boundary-layer rakes to measure off-surface velocity profiles, and pressure orifices to measure surface pressure distributions. The initial phase of this research project was recently concluded with two flights on July 14. This phase consisted of a total of twenty flights over a period of about ten weeks. In the coming months the data obtained in this initial set of flight experiments will be analyzed and the results will be used to finalize the instrumentation layout for the next set of flight experiments scheduled for Winter and Spring of 1995. The main goal of these upcoming flights will be: (1) to measure more detailed surface pressure distributions across the wing for a range of flight conditions and flap settings; (2) to visualize the surface flows across the multi-element wing at high-lift conditions using fluorescent mini tufts; and (3) to measure in more detail the changes in boundary-layer state on the various flap elements as a result of changes in flight condition and flap deflection. These flight measured results are being correlated with experimental data measured in ground-based facilities as well as with computational data calculated with methods based on the Navier-Stokes equations or a reduced set of these equations. Also these results provide insight into the extent of laminar flow that exists on actual multi-element lifting surfaces at full-scale high-life conditions. Preliminary results indicate that depending on the deflection angle, the slat and flap elements have significant regions of laminar flow over a wide range of angles of attack. Boundary-layer transition mechanisms that were observed include attachment-line contamination on the slat and inflectional instability on the slat and fore flap. Also, the results agree fairly well with the predictions reported in a paper presented at last year's AIAA Fluid Dynamics Conference. The fact that extended regions of laminar flow are shown to exist on the various elements of the high-lift system raises the question what the effect is of loss of laminar flow as a result of insect contamiantion, rain or ice accumulation on high-life performance.

Vandam, Cornelis P.↗

NASA NTRS2019

Extratropical Large-Scale Traveling Weather Systems in the Southern Hemisphere on Mars

From late-autumn through early-spring, the middle- and high-latitudes of both hemispheres of Mars and its predominantly carbon-dioxide atmosphere support mean equator-to-pole thermal contrasts, and then, support a strong mean westerly polar vortex. Observations from orbiting spacecraft indicate that this intense mean baroclinicity-barotropicity supports large-scale eastward traveling weather systems (i.e., transient, traveling synoptic-period waves, on the order of the Rossby deformation scale). On Earth, extratropical weather disturbances arise from wind-shear instabilities, and these are critical components of the terrestrial global circulation. So it is the case for Mars. Large-scale traveling weather systems on Mars serve as agents in the transport of heat, momentum and scalar and tracer quantities (e.g., atmospheric dust, watervapor, ice clouds, chemical species, etc). Such weather systems interact with other large-scale atmospheric circulation components, namely, quasi-stationary (i.e., forced Rossby) modes; global thermal tidal modes; and then, upon large-/continental- geographical scales, upslope/ down-slope flows amongst high relief, low relief, impact basins, and volcanic rises, and more. The character of Mars' traveling extratropical weather disturbances in its southern hemisphere during late winter through early spring is investigated using a high-resolution Mars global climate model (i.e., Mars GCM), and one from the Agency's Mars Climate Modeling Center (MCMC) based at the NASA Ames Research Center. The climate model includes several complex atmospheric physical packages. With such physics modules, our global climate simulations present comparatively well with observations of the planet's current water cycle (Haberle et al.,2019). The climate model is "forced" with an annual dust cycle (i.e., nudged based on MGS/TES observations). Compared to the northern-hemisphere counterparts, the southern synoptic-period weather disturbances and accompanying frontal waves have smaller meridional and zonal scales, and are less intense. Influences of the zonally asymmetric (i.e., east-west varying) topography on southern large-scale weather are investigated, in addition to large-scale up-slope/down-slope flows and the diurnal cycle. A southern storm zone in late winter and early spring presents in the western hemisphere via orographic influences from the Tharsis highlands, and the Argyre and Hellas impact basins. Geographically localized transient-wave activity diagnostics are constructed that illuminate dynamical differences amongst the simulations and these are presented.

Hollingsworth, Jeffery L.↗

NASA NTRS1996

Development of Three-Dimensional Flow Code Package to Predict Performance and Stability of Aircraft with Leading Edge Ice Contamination

In the work reported herein, a simplified, uncoupled, zonal procedure is 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 NPARC for the incremental forces and moments due to iced components. These are compared with wind tunnel force and moment data, supplied by the Boeing Company, examining longitudinal flight characteristics. Grid refinement improved the local flow features over previously reported work with no appreciable difference in the incremental ice effect. The computed lift curve slope with and without empennage ice matches the experimental value to within 1%, and the zero lift angle agrees to within 0.2 of a degree. The computed slope of the un-iced and iced aircraft longitudinal stability curve is within about 2% of the test data. This work demonstrates the feasibility of a zonal method for the icing analysis of complete aircraft or isolated components within the linear angle of attack range. In fact, this zonal technique has allowed for the viscous analysis of a complete aircraft with ice which is currently not otherwise considered tractable.

Aircraft icing↗

NASA NTRS2018

Parameterization of Mixed Layer and Deep-Ocean Mesoscales Including Nonlinearity

In 2011, Chelton et al. carried out a comprehensive census of mesoscales using altimetry data and reached the following conclusions: "essentially all of the observed mesoscale features are nonlinear" and "mesoscales do not move with the mean velocity but with their own drift velocity," which is "the most germane of all the nonlinear metrics." Accounting for these results in a mesoscale parameterization presents conceptual and practical challenges since linear analysis is no longer usable and one needs a model of nonlinearity. A mesoscale parameterization is presented that has the following features: 1) it is based on the solutions of the nonlinear mesoscale dynamical equations, 2) it describes arbitrary tracers, 3) it includes adiabatic (A) and diabatic (D) regimes, 4) the eddy-induced velocity is the sum of a Gent and McWilliams (GM) term plus a new term representing the difference between drift and mean velocities, 5) the new term lowers the transfer of mean potential energy to mesoscales, 6) the isopycnal slopes are not as flat as in the GM case, 7) deep-ocean stratification is enhanced compared to previous parameterizations where being more weakly stratified allowed a large heat uptake that is not observed, 8) the strength of the Deacon cell is reduced. The numerical results are from a stand-alone ocean code with Coordinated Ocean-Ice Reference Experiment I (CORE-I) normal-year forcing.

Baroclinic flows; Mesoscale processes; Mixing; Sub↗

NASA NTRS1953

Impingement of water droplets on wedges and diamond airfoils at supersonic speeds

An analytical solution has been obtained for the equations of motion of water droplets impinging on a wedge in a two-dimensional supersonic flow field with a shock wave attached to the wedge. The closed-form solution yields analytical expressions for the equation of the droplet trajectory, the local rate of impingement and the impingement velocity at any point on the wedge surface, and the total rate of impingement. The analytical expressions are utilized to determine the impingement on the forward surfaces of diamond airfoils in supersonic flow fields with attached shock waves. The results presented include the following conditions: droplet diameters from 2 to 100 microns, pressure altitudes from sea level to 30,000 feet, free-stream static temperatures from 420 degrees to 460 degrees R. Also, free-stream Mach numbers from 1.1 to 2.0, semi-apex angles for the wedge from 1.14 degrees to 7.97 degrees, thickness-to-chord ratios for the diamond airfoil from 0.02 to 0.14, chord lengths from 1 to 20 feet, and angles of attack from zero to the inverse tangent of the airfoil thickness-to-chord ratio.

SURFACE CONDITIONS - WING SECTIONS↗

NASA NTRS1993

(abstract) Multispecies Gas Flows in the Interior of Comets

A theoretical comet nucleus model has been written which employs 3 volatiles: CO, CO 2 , and H 2 O. It allows variations for equations of state: the Clausius/Clapeyron equation as well as other equations based on the latest thermodynamic data: the Brown/Ziegler relationship for CO 2 and CO at low pressure and the Lowe expression for water vapor at low pressure. It allows for variations in the physics of the outflow, from an essentially throttling mechanism (pure sublimation to space), to the Knudsen regime, to the Dusty-Gas Dynamic model. Variations in several parameters are explored. Results are compared with data from actual comets. None of the variations in parameter space produce a result that compares well with the comet data. All models yield a total gas production between 10&sup27; and 10&sup28; per second (low). This suggests that secondary mechanisms peripheral to this modeling effort, such as the opening and closing of cracks, a conductivity expression which depends upon porosity as well as temperature, etc., may play a more important role in the cometary efflux process than previously thought.

comets↗

NASA NTRS2017

Large-Scale Traveling Weather Systems in Mars Southern Extratropics

Between late fall and early spring, Mars' middle- and high-latitude atmosphere supports strong mean equator-to-pole temperature contrasts and an accompanying mean westerly polar vortex. Observations from both the MGS Thermal Emission Spectrometer (TES) and the MRO Mars Climate Sounder (MCS) indicate that a mean baroclinicity-barotropicity supports intense, large-scale eastward traveling weather systems (i.e., transient synoptic-period waves). Such extratropical weather disturbances are critical components of the global circulation as they serve as agents in the transport of heat and momentum, and generalized scalar/tracer quantities (e.g., atmospheric dust, water-vapor and ice clouds). The character of such traveling extratropical synoptic disturbances in Mars' southern hemisphere during late winter through early spring is investigated using a moderately high-resolution Mars global climate model (Mars GCM). This Mars GCM imposes interactively-lifted and radiatively-active dust based on a threshold value of the surface stress. The model exhibits a reasonable "dust cycle" (i.e., globally averaged, a dustier atmosphere during southern spring and summer occurs). Compared to the northern-hemisphere counterparts, the southern synoptic-period weather disturbances and accompanying frontal waves have smaller meridional and zonal scales, and are far less intense. Influences of the zonally asymmetric (i.e., east-west varying) topography on southern large-scale weather are investigated, in addition to large-scale up-slope/down-slope flows and the diurnal cycle. A southern storm zone in late winter and early spring presents in the western hemisphere via orographic influences from the Tharsis highlands, and the Argyre and Hellas impact basins. Geographically localized transient-wave activity diagnostics are constructed that illuminate dynamical differences amongst the simulations and these are presented.

Hollingsworth, Jeffery L.↗

NASA NTRS1953

Maximum Evaporation Rates of Water Droplets Approaching Obstacles in the Atmosphere Under Icing Conditions

When a closed body or a duct envelope moves through the atmosphere, air pressure and temperature rises occur ahead of the body or, under ram conditions, within the duct. If cloud water droplets are encountered, droplet evaporation will result because of the air-temperature rise and the relative velocity between the droplet and stagnating air. It is shown that the solution of the steady-state psychrometric equation provides evaporation rates which are the maximum possible when droplets are entrained in air moving along stagnation lines under such conditions. Calculations are made for a wide variety of water droplet diameters, ambient conditions, and flight Mach numbers. Droplet diameter, body size, and Mach number effects are found to predominate, whereas wide variation in ambient conditions are of relatively small significance in the determination of evaporation rates. The results are essentially exact for the case of movement of droplets having diameters smaller than about 30 microns along relatively long ducts (length at least several feet) or toward large obstacles (wings), since disequilibrium effects are then of little significance. Mass losses in the case of movement within ducts will often be significant fractions (one-fifth to one-half) of original droplet masses, while very small droplets within ducts will often disappear even though the entraining air is not fully stagnated. Wing-approach evaporation losses will usually be of the order of several percent of original droplet masses. Two numerical examples are given of the determination of local evaporation rates and total mass losses in cases involving cloud droplets approaching circular cylinders along stagnation lines. The cylinders chosen were of 3.95-inch (10.0+ cm) diameter and 39.5-inch 100+ cm) diameter. The smaller is representative of icing-rate measurement cylinders, while with the larger will be associated an air-flow field similar to that ahead of an airfoil having a leading-edge radius comparable with that of the cylinder. It is found that the losses are less than 5 percent. It is concluded that such losses are, in general, very small (less than 1 percent) in the case of smaller obstacles (of icing-rate measurement- cylinder size); the motional dynamics are such, however, that exceptions will occur by reason of failure of very small droplets (moving along stagnation lines) to impinge upon obstacle surfaces. In such cases, the droplets will evaporate completely.

Lowell, H. H.↗

NASA NTRS2005

Capture of Hydrogen Using ZrNi

Water, as ice, is thought to reside in craters at the lunar poles along with CH4 and H2 . A proposed robotic mission for 2012 will utilize metal/metal hydrides for H2 recovery. Specifications are 99% capture of H2 initially at 5 bar and 100C (or greater), and degassing completely at 300C. Of 47-systems examined using the van't Hoff equation, 4 systems, Mg/MgH2, Mg2Ni/Mg2NiH4, ZrNi/ZrNiH2.8, and Pd/PdH0.77, were considered likely candidates for further examination. It is essential, when selecting a system, to also examine questions regarding activation, kinetics, cyclic stability, and gas impurity effects. After considering those issues, ZrN1 was selected as the most promising candidate, as it is easily activated and rapidly forms ZrNiH 2.8 . In addition, it resists oxide poisoning by CO2, and H2O, while some oxidation by O2 is recommended for improved activation . The presence of hydrogen in the as received Zr-Ni alloy from Alfa Aesar posed additional technical problems. X-ray diffraction of the Zr-Ni powder (-325 mesh), with a Zr:Ni wt% ratio of 70:30, was found to consist of ZrH2, ZrNiH2.8, and ZrNi. ZrH2 in the alloy presented the risk that after degassing that both Zr and ZrNi would be present, and thus lead to erroneous results regarding the reactivity of ZrNi with H2 . Fortunately, ZrH2 is a highly stable hydride that does not degas H2 to any significant extent at temperatures below 300C. Based on equilibrium calculations for the decomposition of ZrH2, only 1 millionth of the hydride decomposed at 300C under a N2 atmosphere flowing at 25 ccm for 64 hours, the longest time for pretreatment employed in the investigation. It was possible, from the X-ray results and knowledge of the Zr:Ni ratio, to compute the composition of a pretreated specimen as being 76 wt% ZrNi and the balance ZrH2.

Patton, Lisa↗

NASA NTRS2017

An Initial Study of the Fundamentals of Ice Crystal Icing Physics in the NASA Propulsion Systems Laboratory

This presentation shows results from an initial study of the fundamental physics of ice-crystal ice accretion using the NASA Propulsion Systems Lab (PSL). Ice accretion due to the ingestion of ice-crystals is being attributed to numerous jet-engine power-loss events. The NASA PSL is an altitude jet-engine test facility which has recently added a capability to inject ice particles into the flow. 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. The present study utilized a NACA0012 airfoil. The mixed-phase conditions were generated by partially freezing the liquid-water droplets ejected from the spray bars. This presentation shows data regarding (1) the freeze out characteristics of the cloud, (2) changes in aerothermal conditions due to the presence of the cloud, and (3) the ice accretion characteristics observed on the airfoil model. The primary variable in this test was the PSL plenum humidity which was systematically varied for two duct-exit-plane velocities (85 and 135 ms) as well as two particle size clouds (15 and 50 m MVDi). The observed clouds ranged from fully glaciated to fully liquid, where the liquid clouds were at least partially supercooled. The air total temperature decreased at the test section when the cloud was activated due to evaporation. The ice accretions observed ranged from sharp arrow-like accretions, characteristic of ice-crystal erosion, to cases with double-horn shapes, characteristic of supercooled water accretions.

upwind schemes↗

NASA NTRS2017

Comparisons of Mixed-Phase Icing Cloud Simulations with Experiments Conducted at the NASA Propulsion Systems Laboratory

This paper builds on previous work that compares numerical simulations of mixed-phase icing clouds with experimental data. The model couples the thermal interaction between ice particles and water droplets of the icing cloud with the flowing air of an icing wind tunnel for simulation of NASA Glenn Research Centers (GRC) Propulsion Systems Laboratory (PSL). Measurements were taken during the Fundamentals of Ice Crystal Icing Physics Tests at the PSL tunnel in March 2016. The tests simulated ice-crystal and mixed-phase icing that relate to ice accretions within turbofan engines. Experimentally measured air temperature, humidity, total water content, liquid and ice water content, as well as cloud particle size, are compared with model predictions. The model showed good trend agreement with experimentally measured values, but often over-predicted aero-thermodynamic changes. This discrepancy is likely attributed to radial variations that this one-dimensional model does not address. One of the key findings of this work is that greater aero-thermodynamic changes occur when humidity conditions are low. In addition a range of mixed-phase clouds can be achieved by varying only the tunnel humidity conditions, but the range of humidities to generate a mixed-phase cloud becomes smaller when clouds are composed of smaller particles. In general, the model predicted melt fraction well, in particular with clouds composed of larger particle sizes.

Thermal model↗