Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ice flow model 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.

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

52 records · Page 3

A model of Comet P/Giacobini-Zinner

A computer model of Comet P/Giacobini-Zinner is presented which contains photoprocesses, gas-phase chemical kinetics, energy balance, multifluid hydrodynamics with a transition to free molecular flow, and solar wind interaction. Recently the physics for electrons in the model has been improved by including electron impact ionization and dissociation and separately accounting for electron energetics. Electron heating and cooling mechanisms include photoprocesses, recombination processes, inelastic and elastic collisions with heavy molecules, and expansion cooling. The model incorporates an internally consistent interaction of the solar wind with the coma gas using the axisymmetric ideal fluid dynamic equations. The nuclear size and composition have been chosen to make the calculations relevant to the September 11, 1985 International Cometary Explorer (ICE) encounter with Comet P/Giacobini-Zinner. Model profiles of the temperature, velocity, and number density of the electrons are in good agreement with measurements along ICE's trajectory. The results indicate the probe passed through a region of the coma at the onset of the plasma tail.

Boice, D. C.↗

Revision of Paschen's Law Relating to the ESD of Aerospace Vehicle Surfaces

The purpose of this work is to develop a version of Paschens law that takes into account the flow of ambient gas past electrode surfaces. Paschens law does not consider the flow of gas past an aerospace vehicle whose surfaces may be triboelectrically charged by dust or ice crystal impingement while traversing the atmosphere. The basic hypothesis of this work is that the number of electron-ion pairs created per unit distance between electrode surfaces is mitigated by the electron-ion pairs removed per unit distance by the flow of gas. The revised theoretical model must be a function of the mean velocity vxm of the ambient gas and reduce to Paschens law when the mean velocity is zero. A new theoretical formulation of Paschens law, taking into account the Mach number and compressible dynamic pressure, derived by the authors, will be discussed. This equation has been evaluated by wind tunnel experimentation. Initial data of the baseline wind tunnel experiments show results consistent with the hypothesis. This work may enhance the safety of aerospace vehicles through a redefinition of electrostatic launch commit criteria. It is also possible for new products, such as antistatic coatings, to be formulated based on this data.

discharge↗

Revision of Paschen's Law Relating to the ESD of Aerospace Vehicle Surfaces

The purpose of this work is to develop a version of Paschen's law that takes into account the flow of ambient gas past electrode surfaces. Paschen's law does not consider the flow of gas past an aerospace vehicle whose surfaces may be triboelectrically charged by dust or ice crystal impingement while traversing the atmosphere. The basic hypothesis of this work is that the number of electron-ion pairs created per unit distance between electrode surfaces is mitigated by the electron-ion pairs removed per unit distance by the flow of gas. The revised theoretical model must be a function of the mean velocity vxm of the ambient gas and reduce to Paschen's law when the mean velocity is zero. A new theoretical formulation of Paschen's law, taking into account the Mach number and compressible dynamic pressure, derived by the authors, will be discussed. This equation has been evaluated by wind tunnel experimentation. Initial data of the baseline wind tunnel experiments show results consistent with the hypothesis. This work may enhance the safety of aerospace vehicles through a redefinition of electrostatic launch commit criteria. It is also possible for new products, such as antistatic coatings, to be formulated based on this data.

electrostatic↗

Dynamic Gas Flow Effects on the ESD of Aerospace Vehicle Surfaces

The purpose of this work is to develop a version of Paschen's Law that takes into account the flow of ambient gas past electrode surfaces. Paschen's Law does not consider the flow of gas past an aerospace vehicle, whose surfaces may be triboelectrically charged by dust or ice crystal impingement while traversing the atmosphere. The basic hypothesis of this work is that the number of electron-ion pairs created per unit distance between electrode surfaces is mitigated by the electron-ion pairs removed per unit distance by the flow of gas. The revised theoretical model must be a function of the mean velocity, v (sub xm), of the ambient gas and reduce to Paschen's law when the gas mean velocity, v (sub xm) equals 0. A new theoretical formulation of Paschen's Law, taking into account the Mach number and dynamic pressure, derived by the authors, will be discussed. This equation was evaluated by wind tunnel experimentation whose results were consistent with the model hypothesis.

Electrostatic↗

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.↗

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 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.↗

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

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.↗

Roughness Parameter Optimization of the McClain Model in GlennICE

Aircraft icing remains a significant threat to aviation safety. Software that predicts the impingement and ice accretion on full aircraft geometries and aircraft components are in demand and NASA Glenn is committed to produce software that meets this need. One of the key parameters affecting an accurate prediction of iced geometry is the effect of ice roughness on the heat transfer coefficient. While many efforts have been made to implement the roughness in the flow solver, this report takes a correlation for roughness height distribution that is based on experimental measurements and demonstrates how to relate those measurements to an augmentation to the heat transfer coefficient provided by the flow solution. The outcome of this effort was the calibration of defaults for user supplied parameters to this correlation through comparison with 95 large glaze conditions from experiment by adjusting user-supplied parameters in the roughness augmentation equation.

Aircraft Icing↗

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.↗

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↗

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.↗

(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↗

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

This paper presents 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 paper presents 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.

Evaporative Cooling↗

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↗

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↗