Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

The effects of the diffuse radiation fields due to multiple scattering and thermal reradiation by dust on the dynamics and thermodynamics of a dusty cometary atmosphere

A self-consistent model of a dirty, clathrate cometary nucleus is extended to account for diffuse radiation fields caused by multiple scattering and thermal reradiation of the solar continuum by dust. The model is configured to fit conditions expected for the various spacecraft which will encounter Halley's comet at 0.89 AU. The atmosphere is assumed a chemically reactive dust-gas mixture in quasi-steady spherically symmetric expansion. The effect of electron-neutral ion collisions is accounted for, along with rate constants of the various species of clathrate ice particles and radiative transfer of the solar input in the UV, visible and near-IR intervals. The opacity of the circumnuclear dust is projected to be an order of magnitude greater than previous estimates, which severely impacts the potential visibility of the nucleus to the cameras of the flyby satellites.

Marconi, M. L.↗
NASA NTRS1985

Numerical simulation of compressible, turbulent, two-phase flow

A computer program for numerically simulating compressible, turbulent, two-phase flows is described and applied. Special attention is given to flows in which dust is ingested into the turbulent boundary layer behind shock waves moving over the earth's surface. it is assumed that the two phases are interpenetrating continua which are coupled by drag forces and heat transfer. The particle phase is assumed to be dilute, and turbulent effects are modeled by zero- and two-equation eddy viscosity models. An important feature of the turbulence modeling is the treatment of surface boundary conditions which control the ingestion of particles into the boundary layer by turbulent friction and diffusion. The numerical method uses second-order implicit upwind differencing of the inviscid terms of the equations and second-order central differencing of the viscous terms. A diagonal form of the implicit algorithm is used to improve efficiency, and the transformation to a curvilinear coordinate system is accomplished by the finite volume techniques. Applications to a series of representative flows include a two-phase nozzle flow, the steady flow of air over a sand bed, and the air flow behind a normal shock wave in uniform motion over a sand bed. Results of the latter two applications are compared with experimental results.

Coakley, t. J.↗
NASA NTRS1986

Fluctuation and rotation in diffusion-influenced monomolecular reactions

Models of partially absorbing boundaries for calculating diffusion-influenced chemical reactions are examined analytically. The limitations of diffusion-equation and random-walk models are indicated; the presence of an ambiguous length parameter in a standard radiation-boundary-condition (RBC) model is demonstrated; the RBC model for a boundary with oscillation (at a rate much faster than the rate of thermal collisions) between absorbing and nonabsorbing states is considered; an alternative stochastically switching absorber (SSA) model is introduced and shown to be physically consistent; the equivalence of one-dimensional RBC and SSA models under some conditions is shown; one time passage is calculated for an axially symmetric three-dimensional SSA model with anisotropic absorption; and the analogy between rotation and stochastic switching is explored.

Bashford, Donald↗
NASA NTRS1989

Theoretical study of the seasonal behavior of the global ionosphere at solar maximum

The seasonal behavior of the global ionosphere was studied using a time-dependent three-dimensional physical model (developed by Shunk and his coworkers) of the ionosphere at altitudes between 120 and 800 km. This model accounts for field-aligned diffusion, cross-field electrodynamic drifts both the equatorial region and at high latitudes, interhemispheric flow, thermospheric winds, polar wind escape, energy-dependent chemical reactions, neutral composition changes, ion production due to solar EUV radiation and auroral precipitation, thermal conduction, diffusion-thermal heat flow, and local heating and cooling processes. The model studies were carried out for both June and December solstice conditions at solar maximum and for low geomagnetic activity. The ionospheric features predicted by the model agreed qualitatively with the available measurements.

Sojka, J. J.↗
NASA NTRS1996

The Implementation and Evaluation of the Emergency Response Dose Assessment System (ERDAS) at Cape Canaveral Air Station/Kennedy Space Center

The Emergency Response Dose Assessment System (ERDAS) is a system which combines the mesoscale meteorological prediction model RAMS with the diffusion models REEDM and HYPACT. Operators use a graphical user interface to run the models for emergency response and toxic hazard planning at CCAS/KCS. The Applied Meteorology Unit has been evaluating the ERDAS meteorological and diffusion models and obtained the following results: (1) RAMS adequately predicts the occurrence of the daily sea breeze during non-cloudy conditions for several cases. (2) RAMS shows a tendency to predict the sea breeze to occur slightly earlier and to move it further inland than observed. The sea breeze predictions could most likely be improved by better parameterizing the soil moisture and/or sea surface temperatures. (3) The HYPACT/REEDM/RAMS models accurately predict launch plume locations when RAMS winds are accurate and when the correct plume layer is modeled. (4) HYPACT does not adequately handle plume buoyancy for heated plumes since all plumes are presently treated as passive tracers. Enhancements should be incorporated into the ERDAS as it moves toward being a fully operational system and as computer workstations continue to increase in power and decrease in cost. These enhancements include the following: activate RAMS moisture physics; use finer RAMS grid resolution; add RAMS input parameters (e.g. soil moisture, radar, and/or satellite data); automate data quality control; implement four-dimensional data assimilation; modify HYPACT plume rise and deposition physics; and add cumulative dosage calculations in HYPACT.

Evans, Randolph J.↗
NASA NTRS1993

Numerical modeling of physical vapor transport under microgravity conditions: Effect of thermal creep and stress

One of the most promising applications of microgravity (micro-g) environments is the manufacture of exotic and high-quality crystals in closed cylindrical ampoules using physical vapor transport (PVT) processes. The quality enhancements are believed to be due to the absence of buoyant convection in the weightless environment - resulting in diffusion-limited transport of the vapor. In a typical experiment, solid-phase sample material is initially contained at one end of the ampoule. The sample is made to sublime into the vapor phase and deposit onto the opposite end by maintaining the source at an elevated temperature with respect to the deposit. Identification of the physical factors governing both the rates and uniformity of crystal growth, and the optimization of the micro-g technology, will require an accurate modeling of the vapor transport within the ampoule. Previous micro-g modeling efforts have approached the problem from a 'classical' convective/diffusion formulation, in which convection is driven by the action of buoyancy on thermal and solutal density differences. The general conclusion of these works have been that in low gravity environments the effect of buoyancy on vapor transport is negligible, and vapor transport occurs in a diffusion-limited mode. However, it has been recently recognized than in the non-isothermal (and often low total pressure) conditions encountered in ampoules, the commonly-assumed no-slip boundary condition to the differential equations governing fluid motion can be grossly unrepresentative of the actual situation. Specifically, the temperature gradients can give rise to thermal creep flows at the ampoule side walls. In addition, temperature gradients in the vapor itself can, through the action of thermal stress, lead to bulk fluid convection.

Mackowski, Daniel W.↗
NASA NTRS1990

The baroeffect and an appropriate momentum boundary condition

As disparate molecular weight gases isothermally diffuse between two ends of a capillary tube, they can support a pressure gradient. The magnitude of this pressure gradient depends critically on viscous wall stress and becomes a measure of the boundary condition. Here, a one-dimensional analytical model is proposed that allows for specular gas reflection from the wall. It explains anomalous (4/3) correction factors required previously to match experiment to baroeffect models. It predicts a new physical phenomenon, a surface-driven baroeffect for equal molecular weight gases. Diffusive slip contributions exceed the order of convective diffusion for Peclet number Pe less than 1 and approximately equal convective diffusion for Pe values between 1 and 4. For binary gases, this model further extends the baroeffect experiment as a means to find momentum accommodation coefficients, slip lengths, and external friction coefficients. A corrected binary diffusion coefficient is defined.

Noever, David A.↗
DOE OSTI2014

The exit-time problem for a Markov jump process

The purpose of this paper is to consider the exit-time problem for a finite-range Markov jump process, i.e, the distance the particle can jump is bounded independent of its location. Additionally, such jump diffusions are expedient models for anomalous transport exhibiting super-diffusion or nonstandard normal diffusion. We refer to the associated deterministic equation as a volume-constrained nonlocal diffusion equation. The volume constraint is the nonlocal analogue of a boundary condition necessary to demonstrate that the nonlocal diffusion equation is well-posed and is consistent with the jump process. A critical aspect of the analysis is a variational formulation and a recently developed nonlocal vector calculus. Finally, this calculus allows us to pose nonlocal backward and forward Kolmogorov equations, the former equation granting the various moments of the exit-time distribution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Distinguishing Surface and Bulk Reactivity: Concentration-Dependent Kinetics of Iodide Oxidation by Ozone in Microdroplets

Iodine oxidation reactions play an important role in environmental, biological, and industrial contexts. The multiphase reaction between aqueous iodide and ozone is of particular interest due to its prevalence in the marine atmosphere and unique reactivity at the air–water interface. Here, we explore the concentration dependence of the I – + O 3 reaction in levitated microdroplets under both acidic and basic conditions. To interpret the experimental kinetics, molecular simulations are used to benchmark a kinetic model, which enables insight into the reactivity of the interface, the nanometer-scale subsurface region, and the bulk interior of the droplet. For all experiments, a kinetic description of gas- and liquid-phase diffusion is critical to interpreting the results. We find that the surface dominates the iodide oxidation kinetics under concentrated and acidic conditions, with the reactive uptake coefficient approaching an upper limit of 10 –2 at pH 3. In contrast, reactions in the subsurface dominate under more dilute and alkaline conditions, with inhibition of the surface reaction at pH 12 and an uptake coefficient that is 10× smaller. In conclusion, the origin of a changing surface mechanism with pH is explored and compared to previous ozone-dependent measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2001

Rheology of Diabase: Implications for Tectonics on Venus and Mars

Two important goals of our experimental investigation of the rheological behavior of diabase rocks were: (1) to determine flow laws describing their creep behavior over wide ranges of temperature, stress and strain rate and (2) to develop an understanding of the physical mechanisms by which these rocks flow under laboratory conditions. With this basis, a primary objective then was to construct constitutive equations that can be used to extrapolate from laboratory to planetary conditions. We specifically studied the rheological properties of both natural rock samples and synthetic aggregates. The former provided constraints for geologic systems, while the latter defined the relative contributions of the constituent mineral phases and avoided the influence of glass/melt found in natural samples. In addition, partially molten samples of crustal rock composition were deformed in shear to large strains (greater than 200%) important in crustal environments. The results of this research yielded essential rheological properties essential for models of crustal deformation on terrestrial planets, specifically Venus and Mars, as well as on the geodynamical evolution of these planets. Over the past three years, we also completed our investigation of the creep behavior of water ice with applications to the glaciers, ice sheets and icy satellites. Constitutive equations were determined that describe flow over a wide ranged of stress, strain rate, grain size and temperature. In the case of ice, three creep regimes were delineate. Extrapolation demonstrates that dislocation glide and grain boundary sliding processes dominate flow in ice I under planetary conditions and that diffusion creep is not an important deformation mechanism either in the laboratory or on icy satellites. These results have already been incorporated by other investigators into models describing, for example, the thickness and stability of the ice shell on Europa and to unravel long-standing discrepancies between field observations on glaciers and laboratory results.

Kohlstedt, David L.↗
DOE OSTI2018

Quantitative Evaluation of Radial Diffusion and Local Acceleration Processes During GEM Challenge Events

We simulate the radiation belt electron flux enhancements during selected Geospace Environment Modeling (GEM) challenge events to quantitatively compare the major processes involved in relativistic electron acceleration under different conditions. Van Allen Probes observed significant electron flux enhancement during both the storm time of 17–18 March 2013 and non–storm time of 19–20 September 2013, but the distributions of plasma waves and energetic electrons for the two events were dramatically different. During 17–18 March 2013, the SYM-H minimum reached -130 nT, intense chorus waves (peak B w ~140 pT) occurred at 3.5 < L < 5.5, and several hundred keV to several MeV electron fluxes increased by ~2 orders of magnitude mostly at 3.5 < L < 5.5. During 19–20 September 2013, the SYM-H remained higher than -30 nT, modestly intense chorus waves (peak B w ~80 pT) occurred at L > 5.5, and electron fluxes at energies up to 3 MeV increased by a factor of ~5 at L > 5.5. The two electron flux enhancement events were simulated using the available wave distribution and diffusion coefficients from the GEM focus group Quantitative Assessment of Radiation Belt Modeling. By comparing the individual roles of local electron heating and radial transport, our simulation indicates that resonant interaction with chorus waves is the dominant process that accounts for the electron flux enhancement during the storm time event particularly near the flux peak locations, while radial diffusion by ultralow-frequency waves plays a dominant role in the enhancement during the non–storm time event. Incorporation of both processes reasonably reproduces the observed location and magnitude of electron flux enhancement.

58 GEOSCIENCES↗
DOE OSTI2014

Distributed Wind Diffusion Model Overview (Presentation)

Distributed wind market demand is driven by current and future wind price and performance, along with several non-price market factors like financing terms, retail electricity rates and rate structures, future wind incentives, and others. We developed a new distributed wind technology diffusion model for the contiguous United States that combines hourly wind speed data at 200m resolution with high resolution electricity load data for various consumer segments (e.g., residential, commercial, industrial), electricity rates and rate structures for utility service territories, incentive data, and high resolution tree cover. The model first calculates the economics of distributed wind at high spatial resolution for each market segment, and then uses a Bass diffusion framework to estimate the evolution of market demand over time. The model provides a fundamental new tool for characterizing how distributed wind market potential could be impacted by a range of future conditions, such as electricity price escalations, improvements in wind generator performance and installed cost, and new financing structures. This paper describes model methodology and presents sample results for distributed wind market potential in the contiguous U.S. through 2050.

17 WIND ENERGY↗
DOE OSTI2024

Modeling Ni Coarsening and Migration in the Hydrogen Electrode of Solid oxide Cells under Operating Conditions

Ni coarsening and migration is the most important degradation in the hydrogen electrode of solid oxide cells. This presentation reviews our works on Ni coarsening and migration in NETL over the past years. We used simulation techniques including phase-field modeling and density function theory to investigate possible mechanisms of Ni coarsening and migration, including self-diffusion of Ni, diffusion of gaseous Ni(OH)2, surface diffusion of Ni(OH)x and Ni-YSZ wettability change, and examined the effect of operating conditions on these mechanisms. So far, none of the mechanisms can fully explain the experiments. Remaining questions and possible paths forward are summarized and discussed.

Lei, Yinkai↗
DOE OSTI2015

Monte Carlo Simulations of Coupled Diffusion and Surface Reactions during the Aqueous Corrosion of Borosilicate Glasses

Borosilicate nuclear waste glasses develop complex altered layers as a result of coupled processes such as hydrolysis of network species, condensation of Si species, and diffusion. However, diffusion has often been overlooked in Monte Carlo models of the aqueous corrosion of borosilicate glasses. Therefore, three different models for dissolved Si diffusion in the altered layer were implemented in a Monte Carlo model and evaluated for glasses in the compositional range (75-x) mol% SiO2 (12.5+x/2) mol% B2O3 and (12.5+x/2) mol% Na2O, where 0 ≤ x ≤ 20%, and corroded in static conditions at a surface-to-volume ratio of 1000 m-1. The three models considered instantaneous homogenization (M1), linear concentration gradients (M2), and concentration profiles determined by solving Fick’s 2nd law using a finite difference method (M3). Model M3 revealed that concentration profiles in the altered layer are not linear and show changes in shape and magnitude as corrosion progresses, unlike those assumed in model M2. Furthermore, model M3 showed that, for borosilicate glasses with a high forward dissolution rate compared to the diffusion rate, the gradual polymerization and densification of the altered layer is significantly delayed compared to models M1 and M2. Models M1 and M2 were found to be appropriate models only for glasses with high release rates such as simple borosilicate glasses with low ZrO2 content.

borosilicates↗
DOE OSTI2017

Boundary Conditions for Diffusion-Mediated Processes within Linear Nanopores: Exact Treatment of Coupling to an Equilibrated External Fluid

In this paper, we consider a variety of diffusion-mediated processes occurring within linear nanopores, but which involve coupling to an equilibrated external fluid through adsorption and desorption. By determining adsorption and desorption rates through a set of tailored simulations, and by exploiting a spatial Markov property of the models, we develop a formulation for performing efficient pore-only simulations of these processes. Coupling to the external fluid is described exactly through appropriate nontrivial boundary conditions at the pore openings. This formalism is applied to analyze the following: (i) tracer counter permeation (TCP) where different labeled particles adsorb into opposite ends of the pore and establish a nonequilibrium steady state; (ii) tracer exchange (TE) with exchange of differently labeled particles within and outside the pore; (iii) catalytic conversion reactions where a reactant in the external fluid adsorbs into the pore and converts to a product which may desorb. The TCP analysis also generates a position-dependent generalized tracer diffusion coefficient, the form of which controls behavior in the TE and catalytic conversion processes. We focus on the regime of single-file diffusion within the pore which produces the strongest correlations and largest deviations from mean-field type behavior. Finally, behavior is quantified precisely via kinetic Monte Carlo simulations but is also captured with appropriate analytic treatments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1990

In situ analysis of CO during chemisorption and oxidation on graphite: Supported Pt by FTIR-microspectrometry

For chemisorption and oxidation on Pt/HOPG (highly-orientated pyrolytic) graphite, reflectance Fourier Transform Infrared (FTIR)-microspectrometry reveals a variable state and reactivity for CO. Even for model surface science systems, where surface heterogeneity is minimal, surface diffusion may be too slow relative to the reaction rate to avoid segregation of reactants into surface islands under steady-state conditions. Thus in CO oxidation on Pt (where the relevant surface diffusion coefficients are such that D sub O less than D sub CO) then reactant CO islands exists at the perimeters of which the surface reaction is thought to occur. Furthermore CO can chemisorb on metals in linear and bridge forms to extents which vary with precise faces predominantly exposed coverage, etc. Infrared has long been used to probe the nature of adsorbed CO on model film and heterogeneous surfaces, but it may now be that FTIR-microspectrometry will allow the state of this adsorbate and reactant to be investigated with a spatial resolution of 4.4 microns on model (and real) catalytic surfaces.

Self, Valerie A.↗
NASA NTRS1999

Turbulent Diffusivities for Momentum, Heat, Salt and Passive Scalars

A program began ten years ago to build a turbulence model to describe high Reynolds numbers flows. Specifically, the aim was to devise a model that would satisfy two basic conditions: reproduce available turbulence data (laboratory, numerical simulations, etc.) concerning shear driven flows, buoyancy driven flows, 2D turbulence, freely decaying turbulence etc., and yet be manageable so as to be used, for example, in GCM'S. The model was presented in a series of papers that have appeared in Physics of Fluids since 1996. A total of about 80 turbulence statistics were reproduced. The model has no adjustable parameters. The next step was to apply the model to construct the vertical diffusities K for moment K(sub m), temperature K(sub h), salt K(sub s), and passive scalars K(sub c). First, we took K(sub c) = K(sub h) and tested the model using the GFDL ocean model. The results for the profiles of T and S vs. depth are indistinguishable from those derided using the latest model, the KPP model by the NCAR group. Presently, we are running the same GFDL code relaxing the assumptive K(sub c) = K(sub h). Indeed, the turbulence model yelds a salt diffusivity that depends on R(sub i) and R(sub rho) = ((beta)(delta)S/(alpha)(delta)T/(delta)z) in such a way that K(sub s) may be quite different from K(sub h). Salt fingers and double diffusivity laboratory data are reproduced. Results from the ocean model will be available shortly. Finally, we are trying to derive the horizontal diffusivities with the goal of providing a physically acceptable representation of mesoscale eddies. The recently suggested GMW parameterization has improved several O-GCM results and the goal here is to try to derive/justify it from a turbulence model and/or propose improvements/modifications. Theoretical work is in progress.

Canuto, V. M.↗
NASA NTRS1999

Turbulent Diffusivities for Momentum, Heat, Salt and Passive Scalars

A program began ten years ago to build a turbulence model to describe high Reynolds numbers flows. Specifically, the aim was to devise a model that would satisfy two basic conditions: reproduce available turbulence data (laboratory, numerical simulations, etc.) concerning shear driven flows, buoyancy driven flows, 2D turbulence, freely decaying turbulence etc., and yet be manageable so as to be used, for example, in GCM'S. The model was presented in a series of papers that have appeared in Physics of FLuids since 1996. A total of about 80 turbulence statistics were reproduced. The model has no adjustable parameters. The next step was to apply the model to construct the vertical diffusivities K for moment K(sub m), temperature K(sub h), salt K(sub g) and passive scalars K(sub c). First, we took K(sub g) = K(sub h) and tested the model using the Geophysical Fluid Dynamics Laboratory (GFDL) ocean model. The results for the profiles of T and S vs. depth are indistinguishable from those derided using the latest model, the KPP model by the NCAR group. Presently, we are running the same GFDL code relaxing the assumptive K(sub g) = K(sub h). Indeed, the turbulence model yields a salt diffusivity that depends on Ri and R Rho (= Beta Sigma S/Sigma z/Alpha Sigma T/Sigma z) in such a way that K(sub g) may be quite different from K(sub h). Salt fingers and double diffusivity laboratory data are reproduced. Results from the ocean model will be available shortly. Finally, we are trying to derive the horizontal diffusivities with the goal of providing a physically acceptable representation of mesoscale eddies. The recently suggested GMW parameterization has improved several O-GCM results and the goal here is to try to derive/justify it from a turbulence model and/or propose improvements/modifications. Theoretical work is in progress.

Canuto, Vittorio↗