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 451 records · Page 25

Large Eddy Simulation of Cryogenic Injection Processes at Supercritical Pressure

This paper highlights results from the first of a series of hierarchical simulations aimed at assessing the modeling requirements for application of the large eddy simulation technique to cryogenic injection and combustion processes in liquid rocket engines. The focus is on liquid-oxygen-hydrogen coaxial injectors at a condition where the liquid-oxygen is injected at a subcritical temperature into a supercritical environment. For this situation a diffusion dominated mode of combustion occurs in the presence of exceedingly large thermophysical property gradients. Though continuous, these gradients approach the behavior of a contact discontinuity. Significant real gas effects and transport anomalies coexist locally in colder regions of the flow, with ideal gas and transport characteristics occurring within the flame zone. The current focal point is on the interfacial region between the liquid-oxygen core and the coaxial hydrogen jet where the flame anchors itself.

Oefelein, Joseph C.↗

A diffusion model for use with directional samplers

The paper presents a mathematical model for describing dispersion processes of airborne particles in the atmosphere. The process is described as a superposition of independent Brownian motion processes with drifts and a boundary at zero. It is assumed that the terrain is flat and of a homogeneous roughness. All sources are assumed to be point sources. The time dependencies of emission rates, wind speed, wind direction, and atmospheric conditions are taken into account.

Anbar, D.↗

Photochemistry of the atmosphere of Uranus

A detailed one-dimensional photochemical model incorporating vertical transport by eddy and molecular diffusion is used to study the importance of chemical and physical processes that control the hydrocarbon abundances in the stratosphere of Uranus. The study includes updated photochemical reaction rates and photolysis branching ratios; explicit calculation of the condensation rate of condensing species; and appropriate boundary conditions at the tropopause. It is found that the basic photochemical process generally produces low net conversion efficiencies of methane to higher hydrocarbons. However, even with these low efficiencies, the total condensation rate, which represents the net rate of conversion of methane to hydrocarbons, is roughly 2 x 10 to the -16th g/cm sec, in agreement with the upper range of the low condensation rate for the production of the stratospheric aerosol layer as determined by Pollack et al. (1987). Furthermore, a range of vertical transport rates is required to separately match the observational constraints on hydrocarbon abundances. The solar reflection data refer to the sunlit polar region of Uranus and require an eddy diffusion coefficient of K approx. = 50 sq cm/sec for globally averaged insolation or K approx = 1500 (Nh/N) exp 1/2 for overhead sun. The Voyager 2 Utraviolet Spectrometer solar occultation data refer to the equatorial regions, thus implying a large latitudinal gradient in stratospheric hydrocarbon abundance.

Summers, Michael E.↗

Recent Developments in the Radiation Belt Environment Model

The fluxes of energetic particles in the radiation belts are found to be strongly controlled by the solar wind conditions. In order to understand and predict the radiation particle intensities, we have developed a physics-based Radiation Belt Environment (RBE) model that considers the influences from the solar wind, ring current and plasmasphere. Recently, an improved calculation of wave-particle interactions has been incorporated. In particular, the model now includes cross diffusion in energy and pitch-angle. We find that the exclusion of cross diffusion could cause significant overestimation of electron flux enhancement during storm recovery. The RBE model is also connected to MHD fields so that the response of the radiation belts to fast variations in the global magnetosphere can be studied.Weare able to reproduce the rapid flux increase during a substorm dipolarization on 4 September 2008. The timing is much shorter than the time scale of wave associated acceleration.

Fok, M.-C.↗

Graphite Ablation and Thermal Response Simulation Under Arc-Jet Flow Conditions

The Two-dimensional Implicit Thermal Response and Ablation program, TITAN, was developed and integrated with a Navier-Stokes solver, GIANTS, for multidimensional ablation and shape change simulation of thermal protection systems in hypersonic flow environments. The governing equations in both codes are demoralized using the same finite-volume approximation with a general body-fitted coordinate system. Time-dependent solutions are achieved by an implicit time marching technique using Gauess-Siedel line relaxation with alternating sweeps. As the first part of a code validation study, this paper compares TITAN-GIANTS predictions with thermal response and recession data obtained from arc-jet tests recently conducted in the Interaction Heating Facility (IHF) at NASA Ames Research Center. The test models are graphite sphere-cones. Graphite was selected as a test material to minimize the uncertainties from material properties. Recession and thermal response data were obtained from two separate arc-jet test series. The first series was at a heat flux where graphite ablation is mainly due to sublimation, and the second series was at a relatively low heat flux where recession is the result of diffusion-controlled oxidation. Ablation and thermal response solutions for both sets of conditions, as calculated by TITAN-GIANTS, are presented and discussed in detail. Predicted shape change and temperature histories generally agree well with the data obtained from the arc-jet tests.

Chen, Y.-K.↗

Radiative transfer effects on reflected shock waves. II - Absorbing gas.

Radiative cooling effects behind a reflected shock wave are calculated for an absorbing-emitting gas by means of an expansion procedure in the small density ratio across the shock front. For a gray gas shock layer with an optical thickness of order unity or less the absorption integral is simplified by use of the local temperature approximation, whereas for larger optical thicknesses a Rosseland diffusion type of solution is matched with the local temperature approximation solution. The calculations show that the shock wave will attenuate at first and then accelerate to a constant velocity. Under appropriate conditions the gas enthalpy near the wall may increase at intermediate times before ultimately decreasing to zero. A two-band absorption model yields end-wall radiant-heat fluxes which agree well with available shock-tube measurements.

Su, F. Y.↗

The diffuse interstellar cloud toward HD 179406 (20 Aquilae)

An analysis of the diffuse interstellar cloud complex in front of HD 179406 (20 Aql) is presented. Along this sight line, multispectral absortion- and emission-line studies have uncovered at least three distinct velocity components due to individual clouds. A dominant velocity component is seen in both the absorption and emission-line data sets at 3 +/- 1 km/s. It is argued that the cloud associated with this velocity component is responsible for most of the atomic and all of the molecular gas in front of 20 Aql. The present chemical and physical analysis of the cloud combines the diagnostic tools of radio emission-line data with those of UV and optical absorption data. Using non-LTE models to synthesize the observed absorption profiles, (C-12)O and (C-13)O column densities along this line of sight are determined. The (C-12)O/(C-13)O abundance ratio was found to be 50 +/- 15, similar to that found by Wannier et al. toward Zeta Oph. The physical conditions of the cloud have been investigated using ultraviolet absorption lines. Measurements indicate that the dominant absorption cloud has a gas pressure similar to that found in the local diffuse molecular cloud in Ophiuchus with nT = 20,000/cu cm K.

Hanson, Margaret M.↗

A Validated All-Pressure Fluid Drop Model and Lewis Number Effects for a Binary Mixture

The differences between subcritical liquid drop and supercritical fluid drop behavior are discussed. Under subcritical, evaporative high emission rate conditions, a film layer is present in the inner part of the drop surface which contributes to the unique determination of the boundary conditions; it is this film layer which contributes to the solution's convective-diffusive character. In contrast, under supercritical condition as the boundary conditions contain a degree of arbitrariness due to the absence of a surface, and the solution has then a purely diffusive character. Results from simulations of a free fluid drop under no-gravity conditions are compared to microgravity experimental data from suspended, large drop experiments at high, low and intermediary temperatures and in a range of pressures encompassing the sub-and supercritical regime. Despite the difference between the conditions of the simulations and experiments (suspension vs. free floating), the time rate of variation of the drop diameter square is remarkably well predicted in the linear curve regime. The drop diameter is determined in the simulations from the location of the maximum density gradient, and agrees well with the data. It is also shown that the classical calculation of the Lewis number gives qualitatively erroneous results at supercritical conditions, but that an effective Lewis number previously defined gives qualitatively correct estimates of the length scales for heat and mass transfer at all pressures.

Harstad, K.↗

Numerical Modeling of Solidification in Space With MEPHISTO-4

A pre-flight analysis of the directional solidification of BiSn with MEPHISTO-4 is presented. Simplified Bridgman growth under microgravity conditions is simulated using a two dimensional finite element model. This numerical model is a single domain, pseudo-steady state model, and includes the effects of both thermal and solutal convection. The results show that for all orientations of the applied steady state gravity vector, of magnitude 1 micro-g, the directional solidification process remains diffusion controlled. The maximum convective velocity was found to be 4.424 x 10(exp -5) cm/s for the horizontal Bridgman growth configuration. This value is an order of magnitude lower than the growth velocity. The maximum and minimum values or solute concentration in the liquid at the crystal-melt interface were 13.867 at.% and 13.722 at.%, respectively. This gives a radial segregation value of xi = 1.046% at the interface. A secondary objective of this work was to compare the results obtained to those that consider thermal convection only (no solutal convection). It was found that the convective flow patterns in simulations which included solutal convection were significantly different from those which ignored solutal convection. The level of radial segregation predicted by the current simulations is an order of magnitude lower than that found in simulations which ignore solutal convection. The final aim was to investigate the effect of g-jitter on the crystal growth process. A simulation was performed to calculate the system response to a 1 second, 100 micro-g gravity impulse acting normal to the direction of growth. This pulse is consistent with that induced by Orbiter thruster firings. The results obtained indicate that such a gravity pulse causes an increase in the level of radial solute segregation at the interface from the steady state values. The maximum value of solute concentration in the liquid was found to be 13.888 at.%, the minimum value calculated was 13.706 at.%, yielding a radial segregation value of xi = 1.31% at the interface. These values occurred 126 seconds after the pulse terminated. Thus it is anticipated that the process will remain diffusion controlled even when subjected to such g-jitter.

Simpson, James E.↗

Melt Stabilization of PbSnTe in a Magnetic Field

Both the experimental observation and numerical simulation indicate that the Bridgman growth of PbSnTe under the microgravity environment in space is still greatly influenced by buoyancy-induced convection. The application of a magnetic field during the semiconductor growth can dampen the convective flow in the metal-like melt. However, for Bridgman growth of PbSnTe on earth (with either vertical or horizontal configuration), both experimental observation and numerical modeling suggest that even with a strong magnetic furnace (5-Tesla constant axial magnetic field), the convective flow in the melt still cannot be sufficiently suppressed to reach the diffusion-controlled level. In order to completely dampen the buoyancy-induced convection on earth, estimates based on scaling analysis indicate that for common experimental conditions, an extremely high magnetic field is required, far beyond the capacity of the experimental apparatus currently available. Therefore, it is proposed that only the combination of microgravity environment and magnetic damping will produce the desired diffusion-controlled growth state for this particular material. The primary objectives of this study are to provide a quantitative understanding of the complex transport phenomena during solidification of non-dilute binarys, to furnish a numerical tool for furnace design and growth condition optimization, to provide estimates of the required magnetic field strength for low gravity growth, and to assess the role of magnetic damping for space and earth control of the double-diffusive convection. As an integral part of a NASA research program, our numerical simulation supports both the flight and ground-based experiments in an effort to bring together a complete picture of the complex physical phenomena involved in the crystal growth process. For Bridgman growth of PbSnTe under microgravity (with both vertical and horizontal configurations), the simulations suggest that a moderate axial magnetic field of only a few kilo-Gauss in strength could effectively eliminate buoyancy-induced convection in the melt and control solute segregation. Therefore, this work confirms the idea that the combination of microgravity environment and the magnetic damping will indeed be sufficient to produce the desired diffusion-controlled growth state for PbSnTe.

Fripp, Archibald L.↗

Soil Respiration Phenology Improves Modeled Phase of Terrestrial Net Ecosystem Exchange in Northern Hemisphere

In the northern hemisphere, terrestrial ecosystems transition from net sources of CO2 to the atmosphere in winter to net ecosystem carbon sinks during spring. The timing (or phase) of this transition, determined by the balance between ecosystem respiration (RECO) and primary production, is key to estimating the amplitude of the terrestrial carbon sink. We diagnose an apparent phase bias in the RECO and net ecosystem exchange (NEE) seasonal cycles estimated by the Terrestrial Carbon Flux (TCF) model framework and investigate its link to soil respiration mechanisms. Satellite observations of vegetation canopy conditions, surface meteorology, and soil moisture from the NASA SMAP Level 4 Soil Moisture product are used to model a daily carbon budget for a global network of eddy covariance flux towers. Proposed modifications to TCF include: the inhibition of foliar respiration in the light (the Kok effect); a seasonally varying litterfall phenology; an O2 diffusion limitation on heterotrophic respiration (RH); and a vertically resolved soil decomposition model. We find that RECO phase bias can result from bias in RECO magnitude and that mechanisms which reduce northern spring RECO, like substrate and O2 diffusion limitations, can mitigate the phase bias. A vertically resolved soil decomposition model mitigates this bias by temporally segmenting and lagging RH. Applying these model enhancements at Continuous Soil Respiration (COSORE) sites verifies their improvement of RECO and NEE skill compared to in situ observations (up to ∆RMSE = −0.76 g C m−2 d −1 35 ). Ultimately, these mechanisms can improve prior estimates of NEE for atmospheric inversion studies.

Soil respiration↗

Binary Oscillatory Crossflow Electrophoresis

Electrophoresis has long been recognized as an effective analytic technique for the separation of proteins and other charged species, however attempts at scaling up to accommodate commercial volumes have met with limited success. In this report we describe a novel electrophoretic separation technique - Binary Oscillatory Crossflow Electrophoresis (BOCE). Numerical simulations indicate that the technique has the potential for preparative scale throughputs with high resolution, while simultaneously avoiding many problems common to conventional electrophoresis. The technique utilizes the interaction of an oscillatory electric field and a transverse oscillatory shear flow to create an active binary filter for the separation of charged protein species. An oscillatory electric field is applied across the narrow gap of a rectangular channel inducing a periodic motion of charged protein species. The amplitude of this motion depends on the dimensionless electrophoretic mobility, alpha = E(sub o)mu/(omega)d, where E(sub o) is the amplitude of the electric field oscillations, mu is the dimensional mobility, omega is the angular frequency of oscillation and d is the channel gap width. An oscillatory shear flow is induced along the length of the channel resulting in the separation of species with different mobilities. We present a model that predicts the oscillatory behavior of charged species and allows estimation of both the magnitude of the induced convective velocity and the effective diffusivity as a function of a in infinitely long channels. Numerical results indicate that in addition to the mobility dependence, the steady state behavior of solute species may be strongly affected by oscillating fluid into and out of the active electric field region at the ends of the cell. The effect is most pronounced using time dependent shear flows of the same frequency (cos((omega)t)) flow mode) as the electric field oscillations. Under such conditions, experiments indicate that solute is drawn into the cell from reservoirs at both ends of the cell leading to a large mass build up. As a consequence, any initially induced mass flux will vanish after short times. This effect was not captured by the infinite channel model and hence numerical and experimental results deviated significantly. The revised model including finite cell lengths and reservoir volumes allowed quantitative predictions of the time history of the concentration profile throughout the system. This latter model accurately describes the fluxes observed for both oscillatory flow modes in experiments using single protein species. Based on the results obtained from research funded under NASA grant NAG-8-1080.S, we conclude that binary separations are not possible using purely oscillatory flow modes because of end effects associated with the cos((omega)t) mode. Our research shows, however, that a combination of cos(2(omega)t) and steady flow should lead to efficient separation free of end effects. This possibility is currently under investigation.

Molloy, Richard F.↗

Kinetic Monte Carlo Simulations of Rod Eutectics and the Surface Roughening Transition in Binary Alloys

In this paper we report on two related topics: Kinetic Monte Carlo simulations of the steady state growth of rod eutectics from the melt, and a study of the surface roughness of binary alloys. We have implemented a three dimensional kinetic Monte Carlo (kMC) simulation with diffusion by pair exchange only in the liquid phase. Entropies of fusion are first chosen to fit the surface roughness of the pure materials, and the bond energies are derived from the equilibrium phase diagram, by treating the solid and liquid as regular and ideal solutions respectively. A simple cubic lattice oriented in the {100} direction is used. Growth of the rods is initiated from columns of pure B material embedded in an A matrix, arranged in a close packed array with semi-periodic boundary conditions. The simulation cells typically have dimensions of 50 by 87 by 200 unit cells. Steady state growth is compliant with the Jackson-Hunt model. In the kMC simulations, using the spin-one Ising model, growth of each phase is faceted or nonfaceted phases depending on the entropy of fusion. There have been many studies of the surface roughening transition in single component systems, but none for binary alloy systems. The location of the surface roughening transition for the phases of a eutectic alloy determines whether the eutectic morphology will be regular or irregular. We have conducted a study of surface roughness on the spin-one Ising Model with diffusion using kMC. The surface roughness was found to scale with the melting temperature of the alloy as given by the liquidus line on the equilibrium phase diagram. The density of missing lateral bonds at the surface was used as a measure of surface roughness.

Bentz, Daniel N.↗

Optical observations related to the molecular chemistry in diffuse interstellar clouds

Observations, which have been published since 1979, of molecular species in diffuse clouds are discussed. Particular attention is given to the ultraviolet measurements of CO with the Copernicus and IUE satellites and to ground-based optical measurements of CH, CH(+), CN, and 02. These data encompass large enough samples to test the chemical schemes expected to occur in diffuse clouds. Upper limits for other species (e.g., H2O, H2O(+), and C3) place restrictions on the pathways for molecular production. Moreover, analysis of the rotational distribution of the C2 molecule results in the determination of the physical conditions of the cloud. These parameters, including density, temperature, and the intensity of the radiation field, are necessary for modeling the chemistry.

Federman, S. R.↗

Evaluation of transport conditions during physical vapor transport growth of opto-electronic crystals

Transport conditions were evaluated during the vapor phase growth of mercurous chloride crystals in a closed tube. Experimentally observed growth rates were much smaller than those calculated by the Hertz-Knudsen (H-K) equation. The Arrhenius behavior of growth rate with the temperature was used to derive the sticking coefficient. A one-dimensional diffusion model was used to calculate the total mass flux and was compared with the condensing flux. It was predicted that growth occurred in the convecto-diffusive range.

Singh, N. B.↗

Oxygenic Photosynthesis and the Oxidation State of Mars

The oxidation state of the Earth's surface is one of the most obvious indications of the effect of life on this planet. The surface of Mars is highly oxidized, as evidenced by its red color, but the connection to life is less apparent. Two possibilities can be considered. First, the oxidant may be photochemically produced in the atmosphere. In this case the fundamental source of O2 is the loss of H2 to space and the oxidant produced is H2O2. This oxidant would accumulate on the surface and thereby destroy any organic material and other reductants to some depth. Recent models suggest that diffusion limits this depth to a few meters. An alternative source of oxygen is biological oxygen production followed by sequestration of organic material in sediments - as on the Earth. In this case, the net oxidation of the surface was determined billions of years ago when Mars was a more habitable planet and oxidative conditions could persist to great depths, over 100 m. Below this must be a compensating layer of biogenic organic material. Insight into the nature of past sources of oxidation on Mars will require searching for organics in the martian subsurface and sediments.

Hartman, Hyman↗

Great (Flame) Balls of Fire! Structure of Flame Balls at Low Lewis-number-2 (SOFBALL-2)

Everyone knows that an automobile engine wastes fuel and energy when it runs with a fuel-rich mixture. 'Lean' burning, mixing in more air and less fuel, is better for the environment. But lean mixtures also lead to engine misfiring and rough operation. No one knows the ultimate limits for lean operation, for 'weak' combustion that is friendly to the environment while still moving us around. This is where the accidental verification of a decades-old prediction may have strong implications for designing and running low-emissions engines in the 21st century. In 1944, Soviet physicist Yakov Zeldovich predicted that stationary, spherical flames are possible under limited conditions in lean fuel-air mixtures. Dr. Paul Ronney of the University of Southern California accidentally discovered such 'flame balls' in experiments with lean hydrogen-air mixtures in 1984 during drop-tower experiments that provided just 2.2 seconds of near weightlessness. Experiments aboard NASA's low-g aircraft confirmed the results, but a thorough investigation was hampered by the aircraft's bumpy ride. And stable flame balls can only exist in microgravity. The potential for investigating combustion at the limits of flammability, and the implications for spacecraft fire safety, led to the Structure of Flame Balls at Low Lewis-number (SOFBALL) experiment flown twice aboard the Space Shuttle on the Microgravity Sciences Laboratory-1 (MSL-1) in 1997. Success there led to the planned reflight on STS-107. Flame balls are the weakest fires yet produced in space or on Earth. Typically each flame ball produced only 1 watt of thermal power. By comparison, a birthday candle produces 50 watts. The Lewis-number measures the rate of diffusion of fuel into the flame ball relative to the rate of diffusion of heat away from the flame ball. Lewis-number mixtures conduct heat poorly. Hydrogen and methane are the only fuels that provide low enough Lewis-numbers to produce stable flame balls, and even then only for very weak, barely flammable mixtures. Nevertheless, under these conditions flame balls give scientists the opportunity to test models in one of the simplest combustion experiments possible. SOFBALL-2 science objectives include: Improving our understanding of the flame ball phenomenon; Determining the conditions under which flame balls exist; Testing predictions of flame ball lifetimes; Acquiring more precise data for critical model comparison.

Ronney, Paul↗

Multiscale Computer Simulation of Failure in Aerogels

Aerogels have been of interest to the aerospace community primarily for their thermal properties, notably their low thermal conductivities. While such gels are typically fragile, recent advances in the application of conformal polymer layers to these gels has made them potentially useful as lightweight structural materials as well. We have previously performed computer simulations of aerogel thermal conductivity and tensile and compressive failure, with results that are in qualitative, and sometimes quantitative, agreement with experiment. However, recent experiments in our laboratory suggest that gels having similar densities may exhibit substantially different properties. In this work, we extend our original diffusion limited cluster aggregation (DLCA) model for gel structure to incorporate additional variation in DLCA simulation parameters, with the aim of producing DLCA clusters of similar densities that nevertheless have different fractal dimension and secondary particle coordination. We perform particle statics simulations of gel strain on these clusters, and consider the effects of differing DLCA simulation conditions, and the resultant differences in fractal dimension and coordination, on gel strain properties.

Good, Brian S.↗