Search NASA⌕ Search

SEARCH · Search NASA

Results for “Conditional diffusion models”

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 649 records · Page 36

Final Summary of Research Report to the National Aeronautics and Space Administration Origins of Solar Systems Program

The chondrites are aggregates of components (e.g. chondrules, chondrule rims and matrix) that formed in the nebula but, at present, there is no consensus on how any of these components formed or whether their formation produced or post dated the chemical fractionations between the chondrites. Chondrites are, at present, the most primitive Solar System objects available for laboratory study and the conditions under which their principle components formed would provide the most direct constraints for models of nebula formation and evolution. The conditions under which chondrules formed is of particular importance because, if their relative abundance in chondrites approximates that in the nebula, they are the products of one of the most energetic and pervasive processes that operated in the early Solar System. The goal of this proposal was to combine theoretical modeling with a comprehensive study of the elemental and isotopic compositions of the major components in unequilibrated ordinary chondrites (UOCs), with the aim of determining the conditions in the nebula at the time of their formation. The isotopes of volatile and moderately volatile elements should be particularly revealing of conditions during chondrule formation, as evaporation under most conditions would lead to isotopic mass fractionation. Modeling of chondrule and matrix formation requires the development of a kinetic model of evaporation and condensation, and calibration of this model against experiments. Cosmic spherules present an opportunity to test our evaporation models under flash heating conditions that would be difficult to simulate experimentally. However, there is surprisingly little known about the isotopic compositions of silicate cosmic spherules, and a number of questions need to be addressed. Is the range of compositions they exhibit due to evaporation? If they are, are the relative volatilities consistent with the models/experiments and are the isotopic fractionations consistent with Rayleigh conditions? For instance, do the alkalis and S evaporate prior to significant melting so that conditions did not meet the Rayleigh criteria of rapid diffusion? If so, their isotopic fractionation might be considerably suppressed. Could this mechanism of K loss apply to chondrule formation? The Fe isotopic fractionation during evaporation of silicates has not been measured, so cosmic spherules might provide a clue to whether FeO diffusion is fast enough to maintain Rayleigh conditions during evaporation. And so on.

O'D. Alexander, Conel↗

Evaporation and Accompanying Isotopic Fractionation of Sulfur from FE-S Melt During Shock Wave Heating

Chondrules probably formed by melting and subsequent cooling of solid precursors. Evaporation during chondrule melting may have resulted in depletion of volatile elements in chondrules. It is known that kinetic evaporation, especially evaporation from a melt, often leads to enrichment of heavy isotopes in an evaporation residue. However, no evidence for a large degree of heavy-isotope enrichment has been reported in chondrules for K, Mg, Si, and Fe (as FeO). The lack of isotopic fractionation has also been found for sulfur in troilites (FeS) within Bishunpur (LL3.1) and Semarkona (LL3.0) chondrules by an ion microprobe study. The largest fractionation, found in only one grain, was 2.7 +/- 1.4 %/amu, while all other troilite grains showed isotopic fractionations of <1 %/amu. The suppressed isotopic fractionation has been interpreted as results of (i) rapid heating of precursors at temperatures below the silicate solidus and (ii) diffusion-controlled evaporation through a surrounding silicate melt at temperatures above the silicate solidus. The kinetic evaporation model suggests that a rapid heating rate of >10(exp 4)-10(exp 6) K/h for a temperature range of 1000-1300 C is required to explain observed isotopic fractionations. Such a rapid heating rate seems to be difficult to be achieved in the X-wind model, but can be achieved in shock wave heating models. In this study, we have applied the sulfur evaporation model to the shock wave heating conditions of to evaluate evaporation of sulfur and accompanying isotopic fractionation during shock wave heating at temperatures below the silicate solidus.

Tachibana, S.↗

Permeation Rate Equations for Hydrogen and Deuterium in a Palladium-Silver Alloy

Mass transfer of a gas through a selective, solid membrane is an effective method for separation of desired species. This selective permeability is evident in the flow of hydrogen and the isotope deuterium through palladium-silver metal alloy media. In this study, based upon Sieverts’s law and the Arrhenius diffusion equation, an empirical correlation was developed to determine the steady-state permeation rate R, dependence on media temperature T, gas supply pressure p(sub S), and gas backpressure p(sub B) on the lower pressure side. Because of an extensive range of experimental conditions and complete reporting of raw data, the research by Ackerman and Koskinas was used as a source for data allowing empirical equation fitting. Unfortunately, those authors reported best-fit equations that poorly represented their own results. To improve the modeling of the original data and demonstrate the quality of the measurements, the current study develops improved hydrogen and deuterium permeation rate equations: P = 4.22×10(exp −6) A[exp(−704/T)](sq. root p(sub S) − sq. root p(sub B))/t for hydrogen P = 2.12×10(exp −6) A[exp(−468/T )](sq. root p(sub S) − sq. root p(sub B))/t for deuterium for values of cross-sectional area A (sq.cm), medium thickness t (cm), pressure p (psia), and temperature T (Kelvin), giving a permeation rate in mole/minute. These equations model permeation rate data more closely than do several other existing literature sources.

Smith, Phillip J.↗

Velocity and pressure characteristics of a model SSME high pressure fuel turbopump

Under the present effort an experiment rig has been constructed, an instrumentation package developed and a series of mean and rms velocity and pressure measurements made in a turbopump which modelled the first stage of the Space Shuttle Main Engine (SSME) High Pressure Fuel Turbopump. The rig was designed so as to allow initial experiments with a single configuration consisting of a bell-mouth inlet, a flight impeller, a vaneless diffuser and a volute. Allowance was made for components such as inlet guide vanes, exit guide vanes, downstream pumps, etc. to be added in future experiments. This flexibility will provide a clear baseline set of experiments and allow evaluation in later experiments of the effect of adding specific components upon the pump performance properties. The rotational speed of the impeller was varied between 4260 and 7680 rpm which covered the range of scaled SSME rotation speeds when due allowance is made for the differing stagnation temperature, model to full scale. The results at the inlet obtained with rotational speeds of 4260, 6084 and 7680 rpm showed that the axial velocity at the bell-mouth inlet remained roughly constant at 2.2 of the bulk velocity at the exit of the turbopump near the center of the inlet, but it decreased rapidly with increasing radius at all three speeds. Reverse flow occurred at a radius greater than 0.9 R for all three speeds and the maximum negative velocity reduced from 1.3 of the bulk velocity at the exit of the turbopump at 4260 rpm to 0.35 at 7680 rpm, suggesting that operating at a speed closer to the design condition of 8700 rpm improved the inlet characteristics. The reverse flow caused positive prerotation at the impeller inlet which was negligibly small near the center but reached 0.7 of the impeller speed at the outer annulus. The results in the diffuser and the volute obtained at 7680 rpm show that the hub and shroud walls of the diffuser were characterized by regions of transient reverse flow with negative revolution-averaged velocity of 8 percent of the maximum forward revolution-averaged velocity at the center of the diffuser passage near the shroud wall.

Tse, D. G-N.↗

Stability of a diffuse linear pinch with axial boundaries

A formulation of the stability behavior of a finite-length pinch is presented. A general initial perturbation is expressed as a uniformly convergent sum over a complete discrete k set. A variational calculation is then performed, based on the energy principle, in which the end-boundary conditions appear as constraints. The requisite Lagrange multipliers mutually couple the elemental periodic excitations. The resulting extended form of delta-W still admits a proper second-variation treatment so that the minimization and stability considerations of Newcomb remain applicable. Comparison theorems are discussed as is the relevance of this end-effect model to the stability of solar coronal loops.

Einaudi, G.↗

Effects of confinement on partially premixed flames

Partially premixed combustion is an intermediate regime between the limiting cases of premixed and nonpremixed combustion. Although combustion problems are generally approached from one of these two limiting cases, there are many practical situations where flames cannot be considered as purely premixed or nonpremixed, and thus the partially premixed approach must be used. In partially premixed combustion, mechanisms from the premixed and nonpremixed regimes can coexist, and as a result some interesting new phenomena can arise. One such phenomenon is the flame stabilization in laminar mixing layers by triple flames. One of the first observations of triple flames was made by Phillips (1965), who investigated a triple flame propagating in a methane mixing layer. Kioni et al. (1993) also examined triple flames both experimentally and numerically. There have also been numerous analytical studies on the shape and propagation of triple flames under various assumptions by Dold (1989), Dold et al. (1991), and Hartley and Dold (1991). In terms of modeling, Muller et al. (1994) have combined the flamelet formulations for premixed and nonpremixed combustion in order to treat lifted diffusion flames. One common feature in the analytical and numerical studies mentioned above is the assumption of zero heat release, which is necessary to make the problem tractable. The effect of heat release on triple flames was investigated by Ruetsch et al. (1995), where for the unconfined case, flame speeds larger than their premixed counterparts were found. One of the most important practical situations in which these conditions arise is in lifted turbulent jet diffusion flames. At a critical velocity the burning zone of a fuel jet lifts off from the nozzle, moves to increasing distances as the jet velocity increases, and finally blows off. The mechanisms that control these phenomena, i.e. that determine the stability of these flames, are still not understood. In addition to regions where diffusion flame stabilization takes place, partially premixed conditions also exist during the ignition process in nonpremixed systems. Numerical simulations by Reveillon et al. (1994) of the ignition process in a weakly stirred mixture of fuel and oxidizer show that triple flames propagate along lines of stoichiometric mixture fraction throughout the fluid. In addition, Peters (1994) notes that NO(x) emissions are likely to be large in such transient cases, and therefore an understanding of triple flames can provide information concerning pollutant formation. This study extends the work previously done and examines the effects of lateral confinement on partially premixed flames. Once again, we study both the flame structure and propagation.

Ruetsch, G. R.↗

Effects of Fuel Preheat on Soot Formation in Microgravity Laminar Diffusion Flames

Nonbuoyant flames offer themselves as an attractive and promising platform to gain a better understanding of soot mechanisms. The effects of buoyancy can be eliminated temporarily in drop towers which sustain brief intervals of reduced gravity-typically lower than 10(exp -3)g- extending up to several seconds at a time. Microgravity facilities have been employed to show that nonbuoyant flames are longer, wider and sootier than their normal-gravity counterparts. Sunderland et al. recently verified the existence of smoke point in laminar nonbuoyant flames. As reported, microgravity flames operating above their smoke point displayed a blunt tip and much broader soot-containing regimes in comparison to their buoyant counterparts. Mortazavi et al. established that residence times in microgravity laminar jet diffusion flames with Re=0(100) tend to be proportional to burner diameter and inversely proportional to burner exit velocity. This offers the capability to alter residence times in nonbuoyant laminar jet diffusion flames when varying the burner exit diameters and velocities. Megaridis et al. presented a quantitative definition of the soot-field structure within laminar microgravity jet diffusion flames which operated well above their smoke point. The experimental methodology involved a full-field laser-light extinction technique and jet diffusion flames of nitrogen-diluted (50% vol.) acetylene fuel burning in quiescent air at atmospheric pressure. The work was conducted at the 2.2s drop tower of the NASA Lewis Research Center (NASA-LeRC). Parallel work on 1-g flames was also presented in (6) to facilitate comparisons on the effect of gravity on the soot fields. As reported, the soot spatial distributions in 0-g flames did not change in a detectable manner after 1s within a typical 2.2s experiment. During that period, the soot field was shown to sustain a pronounced annular structure throughout the luminous nonbuoyant-flame zone. The maximum soot volume fraction measured at 0-g was nearly a factor of two higher than that at 1-g, thus confirming the enhanced sooting tendency of nonbuoyant flames. Greenberg and Ku presented a similar study and reported trends that matched those of for the 50% (vol.) nitrogen-diluted acetylene fuel. Furthermore, they examined pure acetylene flames and reported similar trends with respect to the influence of gravity on maximum soot volume fractions and flame cross-section-averaged soot loadings. Both studies clearly demonstrated the improved spatial resolution of microgravity flames compared to their normal-gravity counterparts. The current study evaluates the influence of moderate fuel preheat on soot formation within 0-g laminar gas jet diffusion flames. While fuel temperature variations have little influence on residence times in 1-g, they have a much more significant effect in 0-g. The primary objective of this program is to quantify this effect and its consequences on sooting by comparing soot volume fraction distributions under preheated and unpreheated-fuel conditions. Furthermore, the current work aims at expanding the limited soot database available for nonbuoyant flames. Soot fields in such flames can be used to perform additional tests of recently developed soot sub-models which have the potential to become powerful predicting tools in combustion design.

Konsur, Bogdan↗

Structure of confined laminar spray diffusion flames: Numerical investigation

The structure of confined laminar spray diffusion flames is investigated numerically by solving the gas-phase conservation equations for mass species, continuity, momentum, and energy and the liquid-phase equations for droplet position, velocity, size, and temperature. A one-step global reaction scheme along with six equilibrium reactions are employed to model the flame chemistry. Monodisperse as well as polydisperse sprays are considered. The numerical results demonstrate that liquid spray flames substantially differ from gaseous flames in their structure, i.e., temperature, concentration, and velocity fields, shape, and dimensions under the same conditions. Spray flames are predicted to be taller and narrower than their counterpart gaseous ones and their shapes are almost cylindrical. This is in agreement with experimental observations. The numerical computations also show that the use of the equilibrium reactions with the one-step reaction scheme decreases the flame temperature compared to the one-step reaction scheme without the equilibrium reactions and more importantly increases the surface area of the flame zone due to a phenomenon termed 'equilibrium broadening.' The spray flames also possess a finite thickness with minimal overlap of the fuel and oxygen species. A case for which a fuel-mixture consisting of 20 to 80 percent gas-liquid by mass is introduced into the combustor is also investigated and compared with predictions using only gaseous or liquid fuel.

Mawid, M. A.↗

The drift field model applied to the lithium-containing silicon solar cell

The drift field model used by Wolf (1963) to calculate the short-circuit current was extended to permit calculations of the open-circuit voltage and the maximum power under conditions of illumination of either tungsten (2800 C) source or air mass zero sunlight. Voltages were calculated using an expression for the drift field diode saturation current. The model, applied to the oxygen-rich (C-13 group) lithium solar cells, was used to calculate the pre- and post-electron bombardment trends of the open-circuit voltage, maximum power, and short-circuit current for lithium gradients in the range from 10 to the 18th power to 10 to the 19th power Li/cm to the 4th power. Diffusion length degradation and carrier removal effects were sufficient to predict the cell performance up to 3 x 10 to the 14th electrons per sq cm. Beyond this fluence it was necessary to include drift field removal effects.

Godlewski, M. P.↗

Numerical modeling of Bridgman growth of PbSnTe in a magnetic field

In this work we study heat and mass transport, fluid motion, and solid/liquid phase change in the process of steady Bridgman growth of Pb(.8)Sn(.2)Te (LTT) in an axially-imposed uniform magnetic field under terrestrial and microgravity conditions. In particular, this research is concerned with the interrelationships among segregation, buoyancy-driven convection, and magnetic damping in the LTT melt. The main objectives are to provide a quantitative understanding of the complex transport phenomena during solidification of the nondilute binary of LTT, to provide estimates of the strength of magnetic field required to achieve the desired diffusion-dominated growth, and to assess the role of magnetic damping for space and earth based control of the buoyancy-induced convection. The problem was solved by using FIDAP and numerical results for both vertical and horizontal growth configurations with respect to the acceleration of gravity vector are presented.

Yao, Minwu↗

Effects of G-Jitter on Interfacial Dynamics of Two Miscible Liquids: Application of MIM

We designed an experiment to examine the effects of g-jitter on mixing of two miscible liquids using the Microgravity Vibration Isolation Mount (MIM). The global bifurcation of the interface was observed with the MIM operating alternatively to either transmit the g-jitter, isolate from the g-jitter or to provide controlled vibration levels with well defined amplitude and frequency content. With the MIM in isolation mode, the interface remains stationary indicating buoyancy induced convection is negligibly small such that mixing occurs via intrinsic mass diffusion without the masking effect of vibration driven convection. Analytical and computational results are in agreement with the experimental findings. Operation of the MIM in forced mode with conditions typical of g-jitter shows that vibration induced convective flows can excite instability mechanisms such as Kelvin-Helmholtz to generate large amplitude quasi-stationary waves oriented vertically for various cases with Stokes-Reynolds number in the range of 0.003 to 0.5. The two and four mode quasi-stationary waves are also predicted with a mathematical model. Though unplanned, the effect of a primary thruster filing was captured and shown to cause a catastrophic bifurcation, enhancing local mass transport. In light of the findings, experiments planned for the International Space Station should consider the potential effects of g-jitter.

Duval, Walter M. B.↗

Analytical model for the density distribution in the Io plasma torus

An analytical model is developed for the diffusive equilibrium plasma density distribution in the Io plasma torus. The model has been employed successfully to follow the ray path of plasma waves in the multi-ion Jovian magnetosphere; it would also be valuable for other studies of the Io torus that require a smooth and continuous description of the plasma density and its gradients. Validity of the analytical treatment requires that the temperature of thermal electrons be much lower than the ion temperature and that superthermal electrons be much less abundant than the thermal electrons; these two conditions are satisfied in the warm outer region of the Io torus from L = 6 to L = 10. The analytical solutions agree well with exact numerical calculations for the most dense portion of the Io torus within 30 deg of the equator.

Mei, YI↗

Unsteady Crystal Growth Due to Step-Bunch Cascading

Based on our experimental findings of growth rate fluctuations during the crystallization of the protein lysozym, we have developed a numerical model that combines diffusion in the bulk of a solution with diffusive transport to microscopic growth steps that propagate on a finite crystal facet. Nonlinearities in layer growth kinetics arising from step interaction by bulk and surface diffusion, and from step generation by surface nucleation, are taken into account. On evaluation of the model with properties characteristic for the solute transport, and the generation and propagation of steps in the lysozyme system, growth rate fluctuations of the same magnitude and characteristic time, as in the experiments, are obtained. The fluctuation time scale is large compared to that of step generation. Variations of the governing parameters of the model reveal that both the nonlinearity in step kinetics and mixed transport-kinetics control of the crystallization process are necessary conditions for the fluctuations. On a microscopic scale, the fluctuations are associated with a morphological instability of the vicinal face, in which a step bunch triggers a cascade of new step bunches through the microscopic interfacial supersaturation distribution.

Vekilov, Peter G.↗

Navier-Stokes Aerodynamic Simulation of the V-22 Osprey on the Intel Paragon MPP

The paper will describe the Development of a general three-dimensional multiple grid zone Navier-Stokes flowfield simulation program (ENS3D-MPP) designed for efficient execution on the Intel Paragon Massively Parallel Processor (MPP) supercomputer, and the subsequent application of this method to the prediction of the viscous flowfield about the V-22 Osprey tiltrotor vehicle. The flowfield simulation code solves the thin Layer or full Navier-Stoke's equation - for viscous flow modeling, or the Euler equations for inviscid flow modeling on a structured multi-zone mesh. In the present paper only viscous simulations will be shown. The governing difference equations are solved using a time marching implicit approximate factorization method with either TVD upwind or central differencing used for the convective terms and central differencing used for the viscous diffusion terms. Steady state or Lime accurate solutions can be calculated. The present paper will focus on steady state applications, although time accurate solution analysis is the ultimate goal of this effort. Laminar viscosity is calculated using Sutherland's law and the Baldwin-Lomax two layer algebraic turbulence model is used to compute the eddy viscosity. The Simulation method uses an arbitrary block, curvilinear grid topology. An automatic grid adaption scheme is incorporated which concentrates grid points in high density gradient regions. A variety of user-specified boundary conditions are available. This paper will present the application of the scalable and superscalable versions to the steady state viscous flow analysis of the V-22 Osprey using a multiple zone global mesh. The mesh consists of a series of sheared cartesian grid blocks with polar grids embedded within to better simulate the wing tip mounted nacelle. MPP solutions will be shown in comparison to equivalent Cray C-90 results and also in comparison to experimental data. Discussions on meshing considerations, wall clock execution time, load balancing, and scalability will be provided.

Vadyak, Joseph↗

Gravitational Influences on Flame Propagation through Non-Uniform, Premixed Gas Systems

Flame propagation through non-uniformly premixed (or layered) gases has importance both in useful combustion systems and in unintentional fires. As summarized previously, non-uniform premixed gas combustion receives scant attention compared to the more usual limiting cases of diffusion or uniformly premixed flames, especially regarding the role gravity plays. This paper summarizes our progress on furthering the knowledge of layered combustion, in which a fuel concentration gradient exists normal to the direction of flame spread. We present experimental and numerical results for flame spread through propanol-air layers formed near the flash point temperature (25 C) or near the stoichiometric temperature (33 C). Both the model and experimental results show that the removal of gravity results in a faster spreading flame, by as much as 80% depending on conditions. This is exactly the opposite effect as that predicted by an earlier model reported. We also found that having a gallery lid results in faster flame spread, an effect more pronounced at normal gravity, demonstrating the importance of enclosure geometry. Also reported here is the beginning of our spectroscopic measurements of fuel vapor.

Miller, Fletcher J.↗

The formation of spikes in the displacement of miscible fluids

We report on experiments in which a more viscous fluid displaces a less viscous one in a vertical cylindrical tube. These experiments were performed using silicone oils in a vertical pipette of small diameter. The more viscous fluid also had a slightly larger density than the less viscous fluid. In the initial configuration, the fluids were at rest, and the interface was nominally flat. A dye was added to the more viscous fluid for ease of observation of the interface between the fluids. The flow was initiated by pumping the more viscous fluid into the less viscous one. The displacement velocity was such that the Reynolds number was smaller than unity and the Peclet number for mass transfer between the fluids was large compared to unity. For upward displacement of the more viscous fluid from an initially stable configuration, an axisymmetric finger was observed under all conditions. However, a needle-shaped spike was seen to propagate from the main finger in many cases, similar to that observed by Petitjeans and Maxworthy for the displacement of a more viscous fluid by a less viscous one.

Rheology/methods↗

Numerical Simulation of the High-Speed Leg of the National Transonic Facility

Numerical simulations for the flow inside the high-speed leg of the National Transonic Facility was conducted. The NASA Tetrahedral Unstructured Software System (TetrUSS) with its USM3D_ME solver was used to perform the numerical simulations. USM3D_ME is developed and maintained by NASA Langley Research Center. Simulations were conducted for three configurations: empty tunnel, body of revolution installed, and NASA Common Research Model installed in the test section. Simulations were performed for a test section Mach number of 0.7 and 0.85 and a corresponding Reynolds number of 8 million per ft. A controller was developed that automated dynamic outflow boundary and streamlined the process of running multiple simulations. The use of a dynamic outflow boundary was the key parameter to drive simulation to the desired tunnel conditions. The numerical simulations captured expected flow features in NTF test section and in the tunnel plenum. The simulations depicted that separation of the flow inside the diffuser appears to be asymmetric and more extensive toward the top and bottom walls. Analysis of the flow field were conducted, and the computed drag coefficient compared to wind tunnel data.

NTF↗

The structure of dilute combusting sprays

An experimental and theoretical study of drop processes in a turbulent flame is described. The experiments involved a monodisperse (105 and 180 micro m initial diameter) stream of methanol drops injected at the base of a turbulent methane-fueled diffusion flame burning in still air. The following measurements were made: mean and fluctuating phase velocities, mean drop number flux, drop-size distributions and mean gas-phase temperatures. Measurements were compared with predictions of two separated flow models: (1) deterministic separated flow, where drop-turbulence interactions are ignored; and (2) stochastic separated flow, where drop-turbulence interactions are considered using random-walk computations. The stochastic separated flow analysis yielded best agreement with measurements, since it provides for turbulent dispersion of drops which was important for present test conditions (and probably for most combusting sprays as well). Distinguishing the presence or absence of envelope flames around the drops, however, was relatively unimportant for present test conditions, since the drops spent most of their lifetime in fuel-rich regions of the flow where this distinction is irrelevant.

Shuen, J. S.↗