Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mixture modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,171 records · Page 65

Experimental and Computational Investigation of Flow Boiling in Microgravity

This study explores use of Computational Fluid Dynamics (CFD) to predict near-saturated flow boiling of FC-72 in microgravity. The computational method employs transient analysis to predict detailed interfacial behavior and heat transfer characteristics along a rectangular channel heated along two opposite walls. Predicted results are validated against experimental temperature measurements and high-speed video images captured during a series of parabolic aircraft maneuvers for three sets of operating conditions which include variations of both mass velocity and wall heat flux. The computational method is based on the multi-phase volume of fluid (VOF) model, which is combined with appropriate phase change and turbulence models, and accounts for both shear-lift force on bubbles and conjugate heat transfer along the heating walls. A key advantage of the CFD method is ability to capture details that are very difficult to measure experimentally, including detailed spatial variations of bubble shape, void fraction, mixture fluid temperature, liquid velocity, and vapor velocity, results for which are presented for each of the three test cases. Different flow regimes predicted along the heated length exhibit a number of dominant mechanisms including bubble nucleation, bubble growth, coalescence, vapor blankets, interfacial waviness, and residual liquid sub-layer, all of which agree well with experiment. Vapor velocity is shown to increase appreciably along the heated length because of increased void fraction, while liquid velocity experiences large fluctuations. Non-equilibrium effects are accentuated with increasing mass velocity, contributing minor deviations of fluid temperature from simulations compared to those predicted by the analytical method. Predicted wall temperature is fairly uniform in the middle of the heated length but increases in the entrance region, due to sensible heat transfer in the subcooled liquid, and decreases toward the exit, mostly because of flow acceleration resulting from increased void fraction.

two-phase cooling; flow boiling; CFD; microgravity↗

Modeling of Commercial Turbofan Engine With Ice Crystal Ingestion: Follow-On

The occurrence of ice accretion within commercial high bypass aircraft turbine engines has been reported under certain atmospheric conditions. Engine anomalies have taken place at high altitudes that have been attributed to ice crystal ingestion, partially melting, and ice accretion on the compression system components. The result was degraded engine performance, and one or more of the following: loss of thrust control (roll back), compressor surge or stall, and flameout of the combustor. As ice crystals are ingested into the fan and low pressure compression system, the increase in air temperature causes a portion of the ice crystals to melt. It is hypothesized that this allows the ice-water mixture to cover the metal surfaces of the compressor stationary components which leads to ice accretion through evaporative cooling. Ice accretion causes a blockage which subsequently results in the deterioration in performance of the compressor and engine. The focus of this research is to apply an engine icing computational tool to simulate the flow through a turbofan engine and assess the risk of ice accretion. The tool is comprised of an engine system thermodynamic cycle code, a compressor flow analysis code, and an ice particle melt code that has the capability of determining the rate of sublimation, melting, and evaporation through the compressor flow path, without modeling the actual ice accretion. A commercial turbofan engine which has previously experienced icing events during operation in a high altitude ice crystal environment has been tested in the Propulsion Systems Laboratory (PSL) altitude test facility at NASA Glenn Research Center. The PSL has the capability to produce a continuous ice cloud which is ingested by the engine during operation over a range of altitude conditions. The PSL test results confirmed that there was ice accretion in the engine due to ice crystal ingestion, at the same simulated altitude operating conditions as experienced previously in flight. The computational tool was utilized to help guide a portion of the PSL testing, and was used to predict ice accretion could also occur at significantly lower altitudes. The predictions were qualitatively verified by subsequent testing of the engine in the PSL. In a previous study, analysis of select PSL test data points helped to calibrate the engine icing computational tool to assess the risk of ice accretion. This current study is a continuation of that data analysis effort. The study focused on tracking the variations in wet bulb temperature and ice particle melt ratio through the engine core flow path. The results from this study have identified trends, while also identifying gaps in understanding as to how the local wet bulb temperature and melt ratio affects the risk of ice accretion and subsequent engine behavior.

Jorgenson, Philip C. E.↗

Modeling of Commercial Turbofan Engine with Ice Crystal Ingestion; Follow-On

The occurrence of ice accretion within commercial high bypass aircraft turbine engines has been reported under certain atmospheric conditions. Engine anomalies have taken place at high altitudes that have been attributed to ice crystal ingestion, partially melting, and ice accretion on the compression system components. The result was degraded engine performance, and one or more of the following: loss of thrust control (roll back), compressor surge or stall, and flameout of the combustor. As ice crystals are ingested into the fan and low pressure compression system, the increase in air temperature causes a portion of the ice crystals to melt. It is hypothesized that this allows the ice-water mixture to cover the metal surfaces of the compressor stationary components which leads to ice accretion through evaporative cooling. Ice accretion causes a blockage which subsequently results in the deterioration in performance of the compressor and engine. The focus of this research is to apply an engine icing computational tool to simulate the flow through a turbofan engine and assess the risk of ice accretion. The tool is comprised of an engine system thermodynamic cycle code, a compressor flow analysis code, and an ice particle melt code that has the capability of determining the rate of sublimation, melting, and evaporation through the compressor flow path, without modeling the actual ice accretion. A commercial turbofan engine which has previously experienced icing events during operation in a high altitude ice crystal environment has been tested in the Propulsion Systems Laboratory (PSL) altitude test facility at NASA Glenn Research Center. The PSL has the capability to produce a continuous ice cloud which is ingested by the engine during operation over a range of altitude conditions. The PSL test results confirmed that there was ice accretion in the engine due to ice crystal ingestion, at the same simulated altitude operating conditions as experienced previously in flight. The computational tool was utilized to help guide a portion of the PSL testing, and was used to predict ice accretion could also occur at significantly lower altitudes. The predictions were qualitatively verified by subsequent testing of the engine in the PSL. In a previous study, analysis of select PSL test data points helped to calibrate the engine icing computational tool to assess the risk of ice accretion. This current study is a continuation of that data analysis effort. The study focused on tracking the variations in wet bulb temperature and ice particle melt ratio through the engine core flow path. The results from this study have identified trends, while also identifying gaps in understanding as to how the local wet bulb temperature and melt ratio affects the risk of ice accretion and subsequent engine behavior.

Jorgenson, Philip C. E.↗

Numerical Modeling of No Vent Filling of a Cryogenic Tank with Thermo-dynamic Vent System Assisted Injector

This paper presents a multi-node finite volume model of No Vent Filling (NVF) of a cryogenic tank with Thermo-dynamic Vent System (TVS) assisted injector, using the Generalized Fluid System Simulation Program (GFSSP), a general purpose flow network code. NVF tests were conducted in a CRYOgenic Orbital Testbed (CRYOTE) tank which was filled by liquid nitrogen using a TVS assisted injector with vent valve closed during the entire filling process. In TVS assisted injector, the liquid flow splits into two streams: one stream is routed to a Joule-Thomson (J-T) orifice where the flow immediately flashes from liquid to vapor or two phase mixture because the downstream to J-T leg is maintained at vacuum level.; the other stream is injected into the tank after being cooled by cold vapor of J-T leg in a heat exchanger. The flow through J-T leg is also used to cool the outer metal matrix of the injector which in turn cools the vapor in the tank ullage. The cooling of vapor in the ullage by cold injector surface reduces ullage pressure that allows liquid to enter and fill the tank. An integrated numerical model of the test set up was developed. The model included the dual lines of TVS assisted injector, tank and tank wall. The tank was discretized into multiple fluid nodes and branches to represent the ullage and liquid nitrogen and multiple solid nodes to represent the tank wall and structure. The heat transfer between solid to fluid was calculated from pool boiling correlations which include film, transition, and nucleate boiling, as well as natural convection during pre- and post-boiling. The model also accounts for the condensation of vapor at ullage-injector interface and when it comes in contact with the liquid spray. The predicted pressure, resident mass, wall temperature in the tank were compared with the test data.

Alok Majumdar↗

Composition dependence of ion diffusion coefficients in gas mixtures at arbitrary field strengths

Expressions for the diffusion coefficient of ions in gas mixtures are obtained from momentum transfer theory, and are given in terms of the diffusion coefficients and drift velocities of the ions in the pure component gases. Blanc's law holds exactly at all field strengths if the mean free time between collisions is independent of velocity (Maxwell model), but otherwise there may be either positive or negative deviations from Blanc's law at high fields. Such deviations are of comparable magnitude for the diffusion coefficients and the mobility, but are not identical. Specific cases of inverse-power potentials are treated in further detail, and some numerical examples are given for rigid-sphere interactions.

Whealton, J. H.↗

Lean Limit Phenomena

The influence of stretch and preferential diffusion on premixed flame extinction and stability was investigated via two model flame configurations, namely the stagnation flame and the bunsen flame. Using a counterflow burner and a stagnation flow burner with a water-cooled wall, the effect of downstream heat loss on the extinction of a stretched premixed flame investigated for lean and rich propane/air and methane/air mixtures. It was demonstrated that extinction by stretch alone is possible only when the deficient reactant is the less mobile one. When it is the more mobile one, downstream heat loss or incomplete reaction is also needed to achieve extinction. The local extinction of bunsen flame tips and edges of hydrocarbon/air premixtures was investigated using a variety of burners. Results show that, while for both rich propane/air and butane/air mixtures tip opening occurs at a constant fuel equivalence ratio of 1.44 and is therefore independent of the intensity, uniformity, and configuration of the approach flow, for rich methane/air flames burning is intensified at the tip and therefore opening is not possible.

Law, C. K.↗

Measurements of physical properties of model Titan atmospheres

One aspect of the study of Titan's atmosphere is the elucidation of the chemical and physical nature of the aerosols. In order to facilitate this, a program to produce laboratory synthesized model materials for Titan's aerosol and to study their chemical and physical properties is now in progress. Various processes, including electric discharge, photolysis by ultraviolet light, and irradiation by energetic particles, will be used to produce the materials. A first set of experiments where a nominal Titan mixture (97%N2, 3% CH4, 0.2% H2) was subjected to pulsed high temperature shocks yielded a reddish brown waxy solid. This material was subjected to pyrolysis/gas chromatography, a technique that has been proposed as a method for analysis of the Titan aerosols. Preliminary results show the material to consist of simple hydrocarbons but little else, at least up to temperatures of 600 C. Since the material was colored, compounds other than those mentioned above must be present.

Scattergood, T. W.↗

Evaluation of the importance of the relative velocity during evaporation of drops in sprays

The importance of relative velocity past individual drops composing a burning spray with an envelope flame is estimated by applying the Bellan and Cuffel (1983) spray evaporation model to an ambient flow going around the spray without penetrating it. The model assumes a spherical spray composed of monodisperse uniformly distributed droplets which all move at the same speed, and relative gas/spray velocities, and spray radii used were representative of boiler and furnace operations. Numerical calculations indicate that the spray is more prone to ambient flow penetration for leaner mixtures (particularly for n = 1/cu cm), lower initial temperatures of gas phase and droplet, and more dilute spray configurations.

Bellan, J.↗

Nonlinear twist-kink instability of a coronal loop

Three-dimensional magnetoinductive particle simulations are used to demonstrate that the mechanical twisting motion applied to a magnetized plasma column induces a current aligned to the external magnetic field direction, pinches the plasma and magnetic fields, and stores the energy in poloidal magnetic fields. As the twist motion continues, the field lines locally begin to wrap around the plasma more than one revolution. A strong MHD instability sets in that is a mixture of kink and ballooning modes, releasing the magnetic energy and causing destruction of coherent column structure and flows of turbulent plasma. A similar episode ensues, exhibiting relaxation oscillations. The buildup of poloidal fields and structure and its sudden release driven by the twist motion may be a model for the solar coronal loop dynamics which exhibits a slow energy buildup with some photospheric motion and a sudden energy release by flares.

Zaidman, E. G.↗

Combustion instability analysis for liquid propellant rocket engines

The multi-dimensional numerical model has been developed to analyze the nonlinear combustion instabilities in liquid-fueled engines. The present pressure-based approach can handle the implicit pressure-velocity coupling in a non-iterative way. The additional scalar conservation equations for the chemical species, the energy, and the turbulent transport quantities can be handled by the same predictor-corrector sequences. This method is time-accurate and it can be applicable to the all-speed, transient, multi-phase, and reacting flows. Special emphasis is given to the acoustic/vaporization interaction which may act as the crucial rate-controlling mechanism in the liquid-fueled rocket engines. The subcritical vaporization is modeled to account for the effects of variable thermophysical properties, non-unitary Lewis number in the gas-film, the Stefan flow effect, and the effect of transient liquid heating. The test cases include the one-dimenisonal fast transient non-reacting and reacting flows, and the multi-dimensional combustion instabilities encountered in the liquid-fueled rocket thrust chamber. The present numerical model successfully demonstrated the capability to simulate the fast transient spray-combusting flows in terms of the limiting-cycle amplitude phenomena, correspondence between combustion and acoustics, and the steep-fronted wave and flame propagation. The investigated parameters include the spray initial conditions, air-fuel mixture ratios, and the engine geometry. Stable and unstable operating conditions are found for the liquid-fueled combustors. Under certain conditions, the limiting cycle behavior of the combusting flowfields is obtained. The numerical results indicate that the spray vaporization processes play an important role in releasing thermal energy and driving the combustion instability.

Kim, Y. M.↗

Analysis of the Westland Data Set

The "Westland" set of empirical accelerometer helicopter data with seeded and labeled faults is analyzed with the aim of condition monitoring. The autoregressive (AR) coefficients from a simple linear model encapsulate a great deal of information in a relatively few measurements; and it has also been found that augmentation of these by harmonic and other parameters call improve classification significantly. Several techniques have been explored, among these restricted Coulomb energy (RCE) networks, learning vector quantization (LVQ), Gaussian mixture classifiers and decision trees. A problem with these approaches, and in common with many classification paradigms, is that augmentation of the feature dimension can degrade classification ability. Thus, we also introduce the Bayesian data reduction algorithm (BDRA), which imposes a Dirichlet prior oil training data and is thus able to quantify probability of error in all exact manner, such that features may be discarded or coarsened appropriately.

Wen, Fang↗

Effective Forces Between Colloidal Particles

Colloidal suspensions have proven to be excellent model systems for the study of condensed matter and its phase behavior. Many of the properties of colloidal suspensions can be investigated with a systematic variation of the characteristics of the systems and, in addition, the energy, length and time scales associated with them allow for experimental probing of otherwise inaccessible regimes. The latter property also makes colloidal systems vulnerable to external influences such as gravity. Experiments performed in micro-ravity by Chaikin and Russell have been invaluable in extracting the true behavior of the systems without an external field. Weitz and Pusey intend to use mixtures of colloidal particles with additives such as polymers to induce aggregation and form weak, tenuous, highly disordered fractal structures that would be stable in the absence of gravitational forces. When dispersed in a polarizable medium, colloidal particles can ionize, emitting counterions into the solution. The standard interaction potential in these charged colloidal suspensions was first obtained by Derjaguin, Landau, Verwey and Overbeek. The DLVO potential is obtained in the mean-field linearized Poisson-Boltzmann approximation and thus has limited applicability. For more precise calculations, we have used ab initio density functional theory. In our model, colloidal particles are charged hard spheres, the counterions are described by a continuum density field and the solvent is treated as a homogeneous medium with a specified dielectric constant. We calculate the effective forces between charged colloidal particles by integrating over the solvent and counterion degrees of freedom, taking into account the direct interactions between the particles as well as particle-counterion, counterion-counterion Coulomb, counterion entropic and correlation contributions. We obtain the effective interaction potential between charged colloidal particles in different configurations. We evaluate two- and three-body forces in the bulk as well as study the influence of soft walls. We qualitatively explain the effects of the walls on the forces and demonstrate that many-body effects are negligible in our system. With adjustments in the parameters, the DLVO pair-potential can describe the results quantitatively. Besides electrostatic interactions, entropic depletion effects that arise from (hard-core) exclusion play an important role in determining the behavior of multi-component colloidal suspensions. A standard theory for depletion forces is due to Asakura and Oosawa and is based on the ideal gas approximation. To go beyond this approximation, we have studied entropic forces in molecular dynamics simulations of systems of hard spheres (the effects of the solvent have been ignored). The effective depletion forces for these systems can be found either from equilibrium distribution functions or from direct momentum transfer calculations. Our results obtained by either method show qualitative differences from the Asakura-Oosawa forces, indicating a longer range, higher value at contact and most importantly a more complicated structure, comprising of several maxima and minima. Our calculations include the determination of effective forces between two spheres, a hard sphere and a wall, and the behavior of a hard sphere near a step-edge and a corner. We also demonstrate that such entropic forces do not necessarily satisfy pairwise additivity.

Tehver, Riina↗

A Satellite-based Assessment of Trans-Pacific Transport of Pollution Aerosol

It has been well documented that pollution aerosol and dust from East Asia can transport across the North Pacific basin, reaching North America and beyond. Such intercontinental transport extends the impact of aerosols for climate change, air quality, atmospheric chemistry, and ocean biology from local and regional scales to hemispheric and global scales. Long term, measurement-based studies are necessary to adequately assess the implications of these wider impacts. A satellite-based assessment can augment intensive field campaigns by expanding temporal and spatial scales and also serve as constraints for model simulations. Satellite imagers have been providing a wealth of evidence for the intercontinental transport of aerosols for more than two decades. Quantitative assessments, however, became feasible only recently as a result of the much improved measurement accuracy and enhanced new capabilities of satellite sensors. In this study, we generated a 4-year (2002 to 2005) climatology of optical depth for pollution aerosol (defined as a mixture of aerosols from urbanlindustrial pollution and biomass burning in this study) over the North Pacific from MODerate resolution Imaging Spectro-radiometer (MODIS) observations of fine- and coarse-mode aerosol optical depths. The pollution aerosol mass loading and fluxes were then calculated using measurements of the dependence of aerosol mass extinction efficiency on relative humidity and of aerosol vertical distributions from field campaigns and available satellite observations in the region. We estimated that about 18 Tg/year pollution aerosol is exported from East Asia to the northwestern Pacific Ocean, of which about 25% reaches the west coast of North America. The pollution fluxes are largest in spring and smallest in summer. For the period we have examined the strongest export and import of pollution particulates occurred in 2003, due largely to record intense Eurasia wildfires in spring and summer. The overall uncertainty of pollution fluxes is estimated at about 80%. A reduction of uncertainty can be achieved with a better characterization of pollution aerosol through integrating emerging A-Train measurements. Simulations by the Goddard Chemistry Aerosol Radiation and Transport (GOCART) and Global Modeling Initiative (GMI) models agree quite well with the satellite-based estimates of annual and latitudeintegrated fluxes, with larger model-satellite differences in latitudinal variations of fluxes.

Yu, Hongbin↗

Thermal model for impact breccia lithification - Manicouagan and the moon

The thermal model of Simonds (1975) is extended to the full spectrum of impact-produced rocks ranging from fragmental breccias to impact melts, with reference to the Manicouagan impact structure in Quebec. This is done by relating the basic textural features of impact-lithified rocks to variations in the mixture of superheated impact-fused material originating near the point of impact and much cooler fragmented debris originating farther from the point of impact.

Simonds, C. H.↗

Composition of the Cayley Formation at Apollo 16 as inferred from impact melt splashes

Abundances of major and trace elements and magnetic properties of 50 impact melt splashes (IMSs) from the Apollo 16 landing site are analzyed to determine the composition of their meteoritic component. MgO-Sc and Ca-Sc variation diagrams and least-squares mixing models are utilized to analyze the IMS, soil, and rock data. Consideration is given to progenitor lithologies of the IMS, the number of impact events represented by the IMS, and the heterogeneity of impact melts from single events. It is observed that the IMSs are composed of either a mixture of anorthosite and low-Sc impact melt rocks or anorthositic norite. It is determined that the surface Cayley layer is composed of TiO2, MgO, Sc, and La concentrations of 0.69, and 7.1 wt pct and 10.5 and 21.2 microg/g, respectively and 0.38 and 5.9 wt pct and 6.1 and 11.8 microg/g, respectively, for the subsurface Cayley layer. The Descartes Formation composition is estimated as TiO2, MgO, Sc, and La concentrations of 0.25, and 3.5 wt pct and 7.7 and 2.2 microg/g, respectively.

Morris, Richard V.↗

The oscillations of superfluid neutron stars

The effects of superfluid hydrodynamics on the oscillations of neutron stars are investigated. The equations describing the small-amplitude pulsations of a neutron star's neutron-proton superfluid mixture are reduced to a system of three second-order equations for three scalar potentials. A variational principle is developed from which the frequencies of the modes of these superfluid oscillations may be estimated. These pulsations are studied by finding analytical solutions to the equations for simple uniform models of nonrotating neutron stars and numerical solutions for realistic models (including rotation). These solutions reveal that the lowest frequency modes are almost indistinguishable from those based on ordinary-fluid hydrodynamics. The analytical solutions also reveal the existence of a new set of modes having no ordinary-fluid counterpart.

Lindblom, Lee↗

Evidence for Methane Segregation at the Surface of Pluto

In May 1995, a set of spectrophotometric curves of the system Pluto-Charon was recorded with the UKIRT telescope equipped with the spectrometer CGS4. As for the previous observations, the spectra cover a part of the near infrared range, between 1.4 and 2.55 micrometers, but with a higher resolution of approximately 700. In both the 1992 and 1995 data, the existence of solid methane is confirmed by numerous absorption bands, and the carbon monoxide and the nitrogen ices are identified by their respective signatures at 2.35 and 2.15 um. The solid nitrogen seems to be the principal icy component and forms a matrix in which the CH4 and CO molecules are diluted. However a spectroscopic analysis of the 1995 observations indicates that pure methane may coexist with its diluted phase in N2. In order to derive the horizontal and vertical distribution of these different species and to obtain some quantitative information about their characteristics, we have modeled the spectrum of May 15 that corresponds to the maximum of Pluto's visible light curve. This was achieved by means of a radiative transfer algorithm dealing with compact and stratified media. Among the various representations we have tested to describe the surface of Pluto, only a geographical mixture of three distinct units explains all the significant structures of the analyzed spectrum. The first unit is a thin granular layer of pure CH4 covering a compact polycrystalline substratum of N2-CH4-CO, which are in a molecular mixture (concentrations of and CO of the order of 0.45%, 0.1-0.2% respectively). It covers about 70% of the observed area and corresponds to volatile deposits that are sublimating under solar illumination. The second unit is either (a) a single thick layer of pure granular methane or (b) a unit similar to the first unit but with the two components inverted (i.e. with CH4 forming a substratum and the N2-CH4-CO mixture a superficial layer of fine grains). Covering 20% of the surface, it represents some old surfaces that have been sublimated for a long time, and eventually recovered later by very small amounts of fresh deposits of the molecular mixture N2-CH4-CO. Finally, the third unit may result from the condensation of very fine grains of nearly pure N2. It covers the remainder of the surface (about 10%). All these results allow a better understanding of the processes of deposition, metamorphism, sublimation and transport affecting the different ices detected on Pluto during its climatic cycles.

Doute, S.↗

Numerical Study of Ammonia Leak and Dispersion in the International Space Station

Release of ammonia into the International Space Station (ISS) cabin atmosphere can occur if the water/ammonia barrier breach of the active thermal control system (ATCS) interface heat exchanger (IFHX) happens. After IFHX breach liquid ammonia is introduced into the water-filled internal thermal control system (ITCS) and then to the cabin environment through a ruptured gas trap. Once the liquid water/ammonia mixture exits ITCS, it instantly vaporizes and mixes with the U.S. Laboratory cabin air that results in rapid deterioration of the cabin conditions. The goal of the study is to assess ammonia propagation in the Station after IFHX breach to plan the operation procedure. A Computational Fluid Dynamics (CFD) model for accurate prediction of airflow and ammonia transport within each of the modules in the ISS cabin was developed. CFD data on ammonia content in the cabin aisle way of the ISS and, in particular, in the Russian On- Orbit Segment during the period of 15 minutes after gas trap rupture are presented for four scenarios of rupture response. Localized effects of ammonia dispersion and risk mitigation are discussed.

Son, Chang H.↗