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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Kinetic Model for the Reduction of Cu II Sites by NO + NH 3 and Reoxidation of NH 3 -Solvated Cu I Sites by O 2 and NO in Cu-SSZ-13

In this work, a kinetic model is developed for the reduction of Cu II sites by NO + NH 3 and the reoxidation of NH 3 -solvated Cu I sites by O 2 and NO in Cu-SSZ-13. Fourier transform infrared (FTIR) spectroscopy and spatially resolved capillary inlet mass spectrometry (SpaciMS) measurements during transient reactor experiments are utilized to identify the rate parameters associated with NO + NH 3 RHC (reduction half-cycle), proposed to occur via two distinct pathways involving adsorbed NH 3 and gas-phase NH 3 . The resulting NO + NH 3 RHC model is validated using spatiotemporal N 2 measurements covering a wide range of temperatures (200–450 °C) and space velocities (53,000–640,000 h –1 ). N 2 O formation is observed and modeled during NO + NH 3 RHC, with quantitative validation under standard selective catalytic reduction (SCR) conditions. Experimentally measured enthalpic and entropic changes associated with O 2 adsorption on NH 3 -solvated Cu I (ZCu(NH 3 ) 2 ) complexes [ Kamasamudram, K. Catal. Today 2010 , 151 (3–4), 212-222], along with activation energies estimated computationally for the intercage diffusion of ZCu(NH 3 ) 2 complexes [ Paolucci, C. Science 2017 , 357 (6 354), 898-903], are incorporated into a mean field kinetic model for the low-temperature oxidation half-cycle (OHC). Significant NH 3 release is observed during the isothermal oxidation of Cu I sites, attributed to desorption of NH 3 ligands from NH 3 -solvated Cu II dimers (Z 2 Cu 2 (NH 3 ) 4 O 2 ). Reduction of these dimeric complexes leads to the consumption of one NO/Cu II , contradicting the expected reduction stoichiometry. Inclusion of a global Arrhenius rate for the NO titration of Z 2 Cu 2 (NH 3 ) 4 O 2 complexes provides accurate representations of standard SCR on reduced and oxidized catalysts, predicting transient NO and NH 3 consumption between 150 and 250 °C as a function of hydrothermal aging. Deactivation of low-temperature standard SCR by NH 3 is observed at high NH 3 pressures, modeled via the formation of superoxo amino (ZCu(NH 3 ) 3 OO*) complexes during NH 3 titration of Z 2 Cu 2 (NH 3 ) 4 O 2 complexes [ Negri, C. J. Am. Chem. Soc. 2020 , 142 (37), 15884-15896]. The redox kinetic model presented here provides a foundational description of active site redox during low-temperature standard SCR, combining the recent kinetic, spectroscopic, and computational findings on the mechanism of standard SCR over Cu-SSZ-13.

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

Photochemistry of Saturn's Atmosphere: Hydrocarbon Chemistry and Comparisons with ISO Observations - 1

To investigate the details of hydrocarbon photochemistry on Saturn, we have developed a one-dimensional diurnally averaged model that couples hydrocarbon and oxygen photochemistry, molecular and eddy diffusion, radiative transfer, and condensation. The model results are compared with observations from the Infrared Space Observatory (ISO) to place tighter constraints on molecular abundances, to better define Saturn's eddy diffusion coefficient profile, and to identify important chemical schemes that control the abundances of the observable hydrocarbons in Saturn's upper atmosphere. From the ISO observations, we determine that the column 12 densities of CH3, CH3C2H, and C4H2 above 10 mbar are 4 (sup +2) (sub -1.5) x 10 (exp 13) cm (sup -2), (1.1 plus or minus 0.3) x 10 (exp 15) cm (exp -2), and (1.2 plus or minus 0.3) x 10 (exp 14) cm (sup -2), respectively. The observed ISO emission features also indicate C2H2 mixing ratios of 1.2 (sup +0.9) (sub -0.6) x 10 (exp -6) at 0.3 mbar and (2.7 plus or minus 0.8) x 10 (exp -7) at 1.4 mbar, and a C2H6 mixing ratio of (9 plus or minus 2.5) x 10 (exp -6) at 0.5 mbar. Upper limits are provided for C2H4, CH2CCH2, C3H8, and C6H2 sensitivity of the model results to variations in the eddy diffusion coefficient profile, the solar flux, the CH4 photolysis branching ratios, the atomic hydrogen influx, and key reaction rates are discussed in detail. We find that C4H2 and CH3C2H are particularly good tracers of important chemical processes and physical conditions in Saturn's upper atmosphere, and C2H6 is a good tracer of the eddy diffusion coefficient in Saturn's lower stratosphere. The eddy diffusion coefficient must be smaller than approximately 3 x 10 (exp 4) sq cm s (sup -1) at pressures greater than 1 mbar in order to reproduce the C2H6 abundance inferred from ISO observations. The eddy diffusion coefficients in the upper stratosphere could be constrained by observations of CH3 radicals if the low-temperature chemistry of CH3 were better understood. We also discuss the implications of our modeling for aerosol formation in Saturn's lower stratosphere-diacetylene, butane, and water condense between approximately 1 and 300 mbar in our model and will dominate stratospheric haze formation at nonauroral latitudes. Our photochemical models will be useful for planning observational sequences and for analyzing data from the upcoming Cassini mission.

Moses, Julianne I.↗
NASA NTRS1994

Mantle devolatilization and rheology in the framework of planetary evolution

We investigate the thermal history of an Earth-like planet with the help of a parameterized mantle convection model including the volatile exchange between mantle and surface reservoirs. The weakening of mantle silicates by dissolved volatiles is described by a functional relationship between creep rate and water fugacity. We use flow law parameters of diffusion creep in olivine under dry and wet conditions. The mantle degassing rate is considered as directly proportional to the seafloor spreading rate, which is also dependent on the mantle heat flow. To calculate the spreading rate, we assume that the heat flow under the mid-ocean ridges is double the average mantle heat flow. The rate of regassing also depends on the seafloor spreading rate as well as on other factors like the efficiency of volatile recycling through island arc volcanism. Both mechanisms (de- and regassing) are coupled self-consistently with the help of the parameterized convection model under implementation of a temperature and volatile-content-dependent mantle viscosity.

Franck, S.↗
DOE OSTI2019

Interaction of Iron with Alumina Refractory Under Flash Ironmaking Conditions

A novel flash ironmaking technology (FIT) based on the direct reduction of iron ore concentrate by a reducing gas such as hydrogen, natural gas, or coal gas has been developed at the University of Utah. A proper choice of refractories is expected to play an important role in the development of FIT and its scale-up. In this work, the interaction of iron with alumina refractory under flash ironmaking conditions has been investigated. A thermodynamic basis for the interaction phenomena has been developed by considering the Fe-Al-O system, and a kinetic model based on solid-state diffusion is proposed to describe the growth of the hercynite (FeAl{sub 2}O{sub 4}) spinel formed as a result of the interaction between Fe, O (from H{sub 2}O and CO{sub 2}), and Al{sub 2}O{sub 3}. Experiments were conducted with Fe powders and pure Al{sub 2}O{sub 3} in the temperature range 1200 to 1500 °C under gas atmospheres relevant to flash ironmaking, and the reacted samples were analyzed using XRD, SEM-EDX, and EPMA techniques. The analyses of reacted samples confirmed the formation of the hercynite (FeAl{sub 2}O{sub 4}) spinel phase and the thickness of the FeAl{sub 2}O{sub 4} spinel as measured using elemental line scans in SEM-EDX and EPMA varied in the range of 0.5 to 1.3 mm. Results showed that the proposed model suitably describes the growth of the hercynite spinel layer, which obeyed the parabolic rate law at all temperatures. The parabolic rate constants were obtained and the diffusion-controlled growth of the FeAl{sub 2}O{sub 4} layer strongly depended on temperature. Furthermore, using the kinetic model, an expression for the effective diffusivity ( (D{sub eff})-bar ) was obtained and its values at the experimental temperatures were determined. The solid-state diffusion was an activated process, as expected, with an activation energy value of 231 kJ/mol.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Grain growth kinetics of the gamma phase metallic uranium

We report metallic uranium is a leading fuel form for sodium cooled fast reactors as an enabling technology of future nuclear energy systems. Mechanistic understanding of fuel behaviors and kinetics under thermodynamic equilibrium and highly non-equilibrium conditions are essential for evaluating fuel performance. It is important to understand and predict the grain and pore evolutions of metallic fuels under thermal and irradiation conditions. However, very limited data are available on the grain growth kinetics and mechanisms of pure gamma phase uranium. In this paper, the pure gamma uranium pellets with different grain structures were fabricated by combining high-energy ball milling and spark plasma sintering. Isothermal annealing tests were performed to investigate the grain growth behavior of the pure gamma phase uranium with different initial grain sizes. A parabolic relationship in grain growth with time was identified for the submicron-sized (374 nm) sample. In contrast, for the nano-sized (137 nm) sample, the grain growth shows a linear relationship with time. The activation energies of grain growth were determined as 199.5 KJ/mol and 80.6 KJ/mol for nano-sized and submicron-sized grain structures, respectively. For the nano-sized sample, the rate-control step of grain growth is dominated by the triple-junction migration, in which the grain boundary triple junction drags the grain growth, leading to a higher activation energy than the bulk diffusion. The dominating mechanism for the submicron-sized sample is grain boundary diffusion. The mechanistic understanding and critical data obtained on the kinetics of pure uranium phases will be useful to evaluate fuel behavior under thermodynamic equilibrium conditions and develop a high fidelity model to predict fuel performance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

Self-similar and diffusive expansion of nonextensive plasmas

Exact analytical self-similar solution is presented for free collisionless expansion of a two-component plasma of inertial ions and nonextensive electrons into vacuum, using the generalized nonextensive velocity distribution for electrons. Furthermore, a hydrodynamic model of plasma expansion in the presence of the ambipolar diffusion caused by collisions among the plasma species, such as electrons and ions, is developed and a Fokker-Planck-like generalized diffusion equation for steady-state expansion of a nonextensive electron-ion plasma is derived. For the case of generalized statistics and in the absence of particle diffusion, the density, velocity, electric potential, and field of expansion profiles are exactly obtained and studied in terms of the self-similar parameter. It is found that superthermal electrons lead to an accelerated expansion of plasma compared to that of Maxwellian electrons. It is also revealed that the nonextensivity parameter plays a fundamental role on the density, velocity, electric potential, and field configuration of the expansion. Therefore, one is able to distinguish three different regimes q < 1, q = 1, and q > 1 for expansion corresponding to sub-nonextensive, extensive, and super-nonextensive statistical profiles for electrons, respectively. Current research can provide useful information and suggests techniques for investigation of the involved statistical mechanism on the role of the energetic electron fluid in the expansion of plasma in strong pulsed laser-matter interaction experiments. It is also shown that the particle diffusion expansion mechanism becomes more dominant for relatively large values of the nonextensivity parameter, q.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

A simple model of multiphoton micromachining in silk hydrogels

High resolution three-dimensional voids can be directly written into transparent silk fibroin hydrogels using ultrashort pulses of near-infrared (NIR) light. Here, we propose a simple finite-element model that can be used to predict the size and shape of individual features under various exposure conditions. We compare predicted and measured feature volumes for a wide range of parameters and use the model to determine optimum conditions for maximum material removal. The simplicity of the model implies that the mechanism of multiphoton induced void creation in silk is due to direct absorption of light energy rather than diffusion of heat or other photoproducts, and confirms that multiphoton absorption of NIR light in silk is purely a 3-photon process.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Two-temperature effects in Hall-MHD simulations of the HIT-SI experiment

A two-temperature Hall-magnetohydrodynamic (Hall-MHD) model, which evolves the electron and ion temperatures separately, is implemented in the PSI-Tet 3D MHD code and used to model plasma dynamics in the Helicity Injected Torus–Steady Inductive (HIT-SI) experiment. The two-temperature model is utilized for HIT-SI simulations in both the PSI-Tet and NIMROD codes at a number of different injector frequencies in the 14.5–68.5 kHz range. At all frequencies, the NIMROD two-temperature model results in increased toroidal current, lower chord-averaged density, higher average temperatures, outward radial shift of the current centroid, and axial symmetrization of the current centroid, relative to the single-temperature NIMROD simulations. The two-temperature PSI-Tet model illustrates similar trends, but at high frequency operation, it exhibits lower electron temperature, smaller toroidal current, and decreased axial symmetrization with respect to the single-temperature PSI-Tet model. With all models, average temperatures and toroidal currents increase with the injector frequency. Power balance and heat fluxes to the wall are calculated for the two-temperature PSI-Tet model and illustrate considerable viscous and compressive heating, particularly at high injector frequency. Parameter scans are also presented for artificial diffusivity, wall temperature, and density. Both artificial diffusivity and the density boundary condition significantly modify the plasma density profiles, leading to larger average temperatures, toroidal current, and relative density fluctuations at low densities. As a result, a low density simulation achieves sufficiently high current gain (G > 5) to generate significant volumes of closed flux lasting 1–2 injector periods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Surrogate modeling of Cellular-Potts agent-based models as a segmentation task using the U-Net neural network architecture

The Cellular-Potts model is a powerful and ubiquitous framework for developing computational models for simulating complex multicellular biological systems. Cellular-Potts models (CPMs) are often computationally expensive due to the explicit modeling of interactions among large numbers of individual model agents and diffusive fields described by partial differential equations (PDEs). In this work, we develop a convolutional neural network (CNN) surrogate model using a U-Net architecture that accounts for periodic boundary conditions. We use this model to accelerate the evaluation of a mechanistic CPM previously used to investigate in vitro vasculogenesis. The surrogate model was trained to predict 100 computational steps ahead (Monte-Carlo steps, MCS), accelerating simulation evaluations by a factor of 562 times compared to single-core CPM code execution on CPU. Over short timescales of up to 3 recursive evaluations, or 300 MCS, our model captures the emergent behaviors demonstrated by the original Cellular-Potts model such as vessel sprouting, extension and anastomosis, and contraction of vascular lacunae. This approach demonstrates the potential for deep learning to serve as a step toward efficient surrogate models for CPM simulations, enabling faster evaluation of computationally expensive CPM simulations of biological processes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

First-Principles Studies of Tritium Species Diffusivity Across the Interfaces of Ni-NiZr Alloy-Zircaloy-4

A tritium producing assembly, also known as tritium producing burnable absorber rods (TPBARs), consists of a metal getter tube located between the cladding and γ-LiAlO2 pellets. The metal getter tube is composed of a nickel (Ni) layer coated on Zircaloy-4. The getter assembly is used to capture tritium (3H) species (mainly 3H2 and 3H2O) generated from γ-LiAlO2 pellets during irradiation. Exploring 3H species (3H2, 3H2O) dissociation on the Ni surface and diffusion in the Ni layer and across the interface of Ni-plated Zircaloy-4 getters can provide insights on tritium transport and retention in the pellets and the getter materials. . In FY25, based on the ideal interface model of Ni-Zircaloy-4 generated in FY24, we will further explore the diffusion pathways of 3H across such Ni-ZrNi-alloy layer-Zircalory-4 interfaces under different conditions, including with oxide/hydroxide clusters on the Ni layer, and impurities located in ZrNi-alloy and Zircalory-4 layers. Due to the size limitation of DFT simulation, we will separate the 3-layer system into 3 sub-systems: Ni-ZrNi interface, ZrNi alloy layer, ZrNi-alloy-Zircaloy-4 interface, and simulate the diffusion barriers for tritium.

diffusion barrier↗
DOE OSTI2022

A PDE Model of Breast Tumor Progression in MMTV-PyMT Mice

The evolution of breast tumors greatly depends on the interaction network among different cell types, including immune cells and cancer cells in the tumor. This study takes advantage of newly collected rich spatio-temporal mouse data to develop a data-driven mathematical model of breast tumors that considers cells’ location and key interactions in the tumor. The results show that cancer cells have a minor presence in the area with the most overall immune cells, and the number of activated immune cells in the tumor is depleted over time when there is no influx of immune cells. Interestingly, in the case of the influx of immune cells, the highest concentrations of both T cells and cancer cells are in the boundary of the tumor, as we use the Robin boundary condition to model the influx of immune cells. In other words, the influx of immune cells causes a dominant outward advection for cancer cells. We also investigate the effect of cells’ diffusion and immune cells’ influx rates in the dynamics of cells in the tumor micro-environment. Sensitivity analyses indicate that cancer cells and adipocytes’ diffusion rates are the most sensitive parameters, followed by influx and diffusion rates of cytotoxic T cells, implying that targeting them is a possible treatment strategy for breast cancer.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

A Model to Account for the Effects of Load Ratio and Hydrogen Pressure on the Fatigue Crack Growth Behavior of Pressure Vessel Steels

A phenomenological model for estimating the effects of load ratio R and hydrogen pressure P H2 on the hydrogen-assisted fatigue crack growth rate (HA–FCGR) behavior in the transient and steady-state regimes of pressure vessel steels is described. The “transient regime” is identified with crack growth within a severely embrittled zone of intense plasticity at the crack tip. The “steady-state” behavior is associated with the crack growing into a region of comparatively lower hydrogen concentration located further away from the crack tip. The model treats the effects of R and P H2 as being functionally separable. In the transient regime, the effects of the hydrogen pressure on the HA–FCGR behavior were negligible but were significant in the steady-state regime. The hydrogen concentration in the steady-state region is modeled as being dependent on the kinetics of lattice diffusion, which is sensitive to pressure. Experimental HA–FCGR data from the literature were used to validate the model. The new model was shown to be valid over a wide range of conditions that ranged between –1 ≤ R ≤ 0.8 and 0.02 ≤ P H2 ≤ 103 MPa for pressure vessel steels.

36 MATERIALS SCIENCE↗
NASA NTRS1977

Prediction of laminar and turbulent boundary layer flow separation in V/STOL engine inlets

The paper provides a theoretical description of the development of the boundary layer on the lip and diffuser surface of a subsonic inlet at arbitrary operating conditions of mass flow rate, freestream velocity and incidence angle. Both laminar separation on the lip and turbulent separation in the diffuser are discussed. The agreement of the theoretical results with model experimental data illustrates the capability of the theory to predict separation. The effects of throat Mach number, inlet size, and surface roughness on boundary-layer development and separation are illustrated.

Chou, D. C.↗
NASA NTRS1987

Current-driven magnetohydrodynamic thermal instabilities in sheared fields

Approximate analytic solutions are sought for the dispersion relation for the MHD stability of magnetized medium in current-driven filamentation modes such as those observed in the solar atmosphere. The magnetic field is assumed to have a self-consistent sheared equilibrium structure. The analysis is carried out in the small wavenumber regime, where shear length is similar to the mode wavelength. Instability is found to depend on the ratio between the thermal and magnetic diffusivities, i.e., the Prandtl number, which identifies the unstable transverse wavenumbers. The instability conditions are expressed in an algebraic equation amenable to numerical solution. Results are provided from use of the model to determine the maximum growth rate and typical scale lengths of instabilities in a precoronal atmosphere and the lower transition region.

Bodo, G.↗
NASA NTRS1993

Turbulent structure of a wingtip vortex in the near field

The turbulent rollup of a vortex generated by a rectangular wing has been investigated. Extensive mean and turbulence measurements of the flowfield on a wingtip and in the near field have been completed. Velocity fluctuation measurements show that the near-field core is not laminar. A large axial velocity excess was found to exist in the core of the vortex. A momentum balance in the near-field of the wingtip showed that the magnitude of the core Reynolds-stress gradient terms are the same order as the largest terms in the governing equations. Navier-Stokes computations of the identical configuration, including wind tunnel walls and using measured inflow and outflow boundary conditions, reproduced many of the features of the experiment. Inherent limitations of the Baldwin-Barth turbulence model combined with limited grid resolution caused the computed vortex core to be more diffuse than desired. The momentum balance also demonstrated that the level of numerically generated false diffusion in the vortex core is relatively high.

Zilliac, Gregory G.↗
NASA NTRS1998

Experimental Studies of Hydrogenation and Other Reactions on Surfaces Under Astrophysically Relevant Conditions

The goal of our project is to study hydrogen recombination reactions on solid surfaces under conditions that are relevant in astrophysics. Laboratory experiments were conducted using low-flux, cold atomic H and D beams impinging on a sample kept under ultra high vacuum conditions. Realistic analogues of interstellar dust grains were used. Our results show that current models for hydrogen recombination reactions have to be modified to take into account the role of activated diffusion of H on surfaces even at low temperature.

Vidali, Gianfranco↗
NASA NTRS2001

Simulation of the Effect of Realistic Space Vehicle Environments on Binary Metal Alloys

Simulations that assess the effect of space vehicle acceleration environments on the solidification of Pb-Sb alloys are reported. Space microgravity missions are designed to provide a near zero-g acceleration environment for various types of scientific experiments. Realistically. these space missions cannot provide a perfect environment. Vibrations caused by crew activity, on-board experiments, support systems stems (pumps, fans, etc.), periodic orbital maneuvers, and water dumps can all cause perturbations to the microgravity environment. In addition, the drag on the space vehicle is a source of acceleration. Therefore, it is necessary to predict the impact of these vibration-perturbations and the steady-state drag acceleration on the experiments. These predictions can be used to design mission timelines. so that the experiment is run during times that the impact of the acceleration environment is acceptable for the experiment of interest. The simulations reported herein were conducted using a finite element model that includes mass, species, momentum, and energy conservation. This model predicts the existence of "channels" within the processing mushy zone and subsequently "freckles" within the fully processed solid, which are the effects of thermosolutal convection. It is necessary to mitigate thermosolutal convection during space experiments of metal alloys, in order to study and characterize diffusion-controlled transport phenomena (microsegregation) that are normally coupled with macrosegregation. The model allows simulation of steady-state and transient acceleration values ranging from no acceleration (0 g). to microgravity conditions (10(exp -6) to 10(exp -3) g), to terrestrial gravity conditions (1 g). The transient acceleration environments simulated were from the STS-89 SpaceHAB mission and from the STS-94 SpaceLAB mission. with on-orbit accelerometer data during different mission periods used as inputs for the simulation model. Periods of crew exercise, quiet (no crew activity), and nominal conditions from STS-89 were used as simulation inputs as were periods of nominal. overboard water-dump, and free-drift (no orbit maneuvering operations) from STS-94. Steady-state acceleration environments of 0.0 and 10(exp -6) to 10(exp -1) g were also simulated, to serve as a comparison to the transient data and to assess an acceptable magnitude for the steady-state vehicle drag

Westra, Douglas G.↗
DOE OSTI2021

Carbonation of model cement pastes: The mineralogical origin of microstructural changes and shrinkage

This study explored the reactive processes of atmospheric carbonation and the consequences with respect to cementitious materials. Two model pastes were used: hydrated C{sub 3}S (including C-S-H and portlandite) and a paste prepared by hydrating a blend of C{sub 3}S and nanosilica (including C-S-H only). The two pastes were carbonated under accelerated conditions in the laboratory. The resulting mineralogical assemblage was examined using X-ray diffraction, thermogravimetric analysis and nuclear magnetic resonance. The microstructural changes were studied by X-ray tomography and porosimetry, and their macroscopic impacts were evaluated through gas diffusion and shrinkage measurements. The use of model pastes allowed for the evaluation of the change in solid volume induced by the carbonation of C-S-H. C-S-H decalcification and subsequent silica chain polymerisation were found responsible for carbonation shrinkage (and potentially cracking). Finally, the results highlight the protective role of portlandite: portlandite helped in limiting C-S-H decalcification and then reducing carbonation shrinkage and cracking.

36 MATERIALS SCIENCE↗