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 91 records · Page 5

DiffESM: Conditional Emulation of Temperature and Precipitation in Earth System Models With 3D Diffusion Models

Earth system models (ESMs) are essential for understanding the interaction between human activities and the Earth's climate. However, the computational demands of ESMs often limit the number of simulations that can be run, hindering the robust analysis of risks associated with extreme weather events. While low-cost climate emulators have emerged as an alternative to emulate ESMs and enable rapid analysis of future climate, many of these emulators only provide output on at most a monthly frequency. This temporal resolution is insufficient for analyzing events that require daily characterization, such as heat waves or heavy precipitation. We propose using diffusion models, a class of generative deep learning models, to effectively downscale ESM output from a monthly to a daily frequency. Trained on a handful of ESM realizations, reflecting a wide range of radiative forcings, our DiffESM model takes monthly mean precipitation or temperature as input, and is capable of producing daily values with statistical characteristics close to ESM output. Combined with a low-cost emulator providing monthly means, this approach requires only a small fraction of the computational resources needed to run a large ensemble. We evaluate model behavior using a number of extreme metrics, showing that DiffESM closely matches the spatio-temporal behavior of the ESM output it emulates in terms of the frequency and spatial characteristics of phenomena such as heat waves, dry spells, or rainfall intensity.

54 ENVIRONMENTAL SCIENCES↗

Reaction–drift–diffusion models from master equations: application to material defects

We present a general method to produce well-conditioned continuum reaction–drift–diffusion equations directly from master equations on a discrete, periodic state space. We assume the underlying data to be kinetic Monte Carlo models (i.e. continuous-time Markov chains) produced from atomic sampling of point defects in locally periodic environments, such as perfect lattices, ordered surface structures or dislocation cores, possibly under the influence of a slowly varying external field. Our approach also applies to any discrete, periodic Markov chain. Here, the analysis identifies a previously omitted non-equilibrium drift term, present even in the absence of external forces, which can compete in magnitude with the reaction rates, thus being essential to correctly capture the kinetics. To remove fast modes which hinder time integration, we use a generalized Bloch relation to efficiently calculate the eigenspectrum of the master equation. A well conditioned continuum equation then emerges by searching for spectral gaps in the long wavelength limit, using an established kinetic clustering algorithm to define a proper reduced, Markovian state space.

36 MATERIALS SCIENCE↗

Modelling the effect of diffuse light on canopy photosynthesis in controlled environments

A layered canopy model was used to analyze the effects of diffuse light on canopy gross photosynthesis in controlled environment plant growth chambers, where, in contrast to the field, highly diffuse light can occur at high irradiance. The model suggests that high diffuse light fractions (approximately 0.7) and irradiance (1400 micromoles m-2 s-1) may enhance crop life-cycle canopy gross photosynthesis for hydroponic wheat by about 20% compared to direct light at the same irradiance. Our simulations suggest that high accuracy is not needed in specifying diffuse light fractions in chambers between approximately 0.7 and 1, because simulated photosynthesis for closed canopies plateau in this range. We also examined the effect of leaf angle distribution on canopy photosynthesis under growth chamber conditions, as these distributions determine canopy extinction coefficients for direct and diffuse light. We show that the spherical leaf angle distribution is not suitable for modeling photosynthesis of planophile canopies (e.g., soybean and peanut) in growth chambers. Also, the absorption of the light reflected from the surface below the canopy should generally be included in model simulations, as the corresponding albedo values in the photosynthetically active range may be quite high in growth chambers (e.g., approximately 0.5). In addition to the modeling implications, our results suggest that diffuse light conditions should be considered when drawing conclusions from experiments in controlled environments.

Non-NASA Center↗

Adaptive Multirate Mass Transfer (aMMT) Model: A New Approach to Upscale Regional‐Scale Transport Under Transient Flow Conditions

Abstract The long‐term evaluation of regional‐scale groundwater quality needs efficient upscaling methods for transient flow. Upscaling techniques, such as the Multirate Mass Transfer (MRMT) method with constant upscaling parameters, have been used for transport with steady‐state flow, yet the upscaling parameters (i.e., rate coefficients) may be time dependent. This study proposed and validated an adaptive MRMT (aMMT) method by allowing the mass transfer coefficients in MRMT to change with the flow field. Advective‐dispersive contaminant transport simulated in a 3‐D heterogeneous medium was used as a reference solution. Equivalent transport under homogeneous flow conditions was evaluated by applying the MRMT and aMMT models for upscaling. The relationship between mass transfer coefficients and flow rates was fitted under steady‐state flow driven by various hydraulic gradients. A power law relationship was obtained, which was then used to update the mass transfer coefficients in each stress period under transient flow conditions in the aMMT method. Results indicated that for advection‐dominated transport, both the MRMT and aMMT methods can upscale the anomalous transport dynamics affected by subgrid heterogeneity under transient flow conditions. Whereas for diffusion‐dominated systems, the MRMT model failed to capture the tails of tracer breakthrough curves after the boundary condition changed, but the results from the aMMT model were significantly improved. However, if the overall flow direction changed, both MRMT and aMMT failed to represent the breakthrough curve tail generated by the heterogeneous system. The results point toward a promising path for upscaling transport in complex aquifers with transient flow.

Guo, Zhilin↗

Atomistic and cluster dynamics modeling of fission gas (Xe) diffusivity in TRISO fuel kernels

TRISO fuel particles are candidates for use in next generation reactors including gas reactors, fluoride salt-cooled high temperature reactors, and micro-reactors. The UCO fuel kernel consists of a uranium dioxide (UO) and uranium carbide mixture. The addition of UC helps suppress the formation of carbon monoxide gas, which led to failures during initial TRISO development. The addition of uranium carbide alters the chemistry of the UO kernel, which is known to influence performance parameters such as fission gas diffusivity, although the impact has not been quantified and no models exist that take the change in chemistry into account. Therefore, better understanding and more accurate models of the impact of chemistry on fuel performance are of high priority. In this paper, a first-principles density functional theory (DFT) and empirical potential based multi-scale study has been carried out to model the diffusivity of fission gas xenon (Xe) in UCO TRISO fuel kernels. The focus is on the UO component in the UCO fuel kernels, as that represents the largest volume fraction of the fuel kernels. The study relies on DFT and empirical potential calculations to determine Xe and point defect properties, which are then used in thermodynamic and kinetic models to predict diffusion for intrinsic conditions. In addition, the information is utilized in cluster dynamics simulations using the Centipede code to estimate the impact of irradiation on defect transport. Additionally, the presence of UC or UC in the UCO fuel kernels is shown to have a substantial impact on the UO non-stoichiometry by inducing oxygen vacancies and driving UO sub-stoichiometric, which causes much slower Xe diffusion in UCO compared to light water reactor UO fuel. The application of this model in fuel performance simulations using the Bison code is also demonstrated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A consistent and conservative model and its scheme for N -phase- M -component incompressible flows

Here, we propose a consistent and conservative model for multiphase and multicomponent incompressible flows, where there can be arbitrary numbers of phases and components. Each phase has a background fluid called the pure phase, each pair of phases is immiscible, and components are dissolvable in some specific phases. The model is developed based on the multiphase Phase-Field model including the contact angle boundary condition, the diffuse domain approach, and the analyses on the proposed consistency conditions for multiphase and multicomponent flows. The model conserves the mass of individual pure phases, the amount of each component in its dissolvable region, and thus the mass of the fluid mixture, and the momentum of the flow. It ensures that no fictitious phases or components can be generated and that the summation of the volume fractions from the Phase-Field model is unity everywhere so that there is no local void or overfilling. It satisfies a physical energy law and it is Galilean invariant. A corresponding numerical scheme is developed for the proposed model, whose formal accuracy is 2nd-order in both time and space. It is shown to be consistent and conservative and its solution is demonstrated to preserve the Galilean invariance and energy law. Numerical tests indicate that the proposed model and scheme are effective and robust to study various challenging multiphase and multicomponent flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Geometry-complete diffusion for 3D molecule generation and optimization

Abstract Generative deep learning methods have recently been proposed for generating 3D molecules using equivariant graph neural networks (GNNs) within a denoising diffusion framework. However, such methods are unable to learn important geometric properties of 3D molecules, as they adopt molecule-agnostic and non-geometric GNNs as their 3D graph denoising networks, which notably hinders their ability to generate valid large 3D molecules. In this work, we address these gaps by introducing the Geometry-Complete Diffusion Model (GCDM) for 3D molecule generation, which outperforms existing 3D molecular diffusion models by significant margins across conditional and unconditional settings for the QM9 dataset and the larger GEOM-Drugs dataset, respectively. Importantly, we demonstrate that GCDM’s generative denoising process enables the model to generate a significant proportion of valid and energetically-stable large molecules at the scale of GEOM-Drugs, whereas previous methods fail to do so with the features they learn. Additionally, we show that extensions of GCDM can not only effectively design 3D molecules for specific protein pockets but can be repurposed to consistently optimize the geometry and chemical composition of existing 3D molecules for molecular stability and property specificity, demonstrating new versatility of molecular diffusion models. Code and data are freely available on GitHub .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Turbulence models and their applications to the prediction of internal flows - A review

The paper presents a brief account of various turbulence models employed in the computation of turbulent flows, and evaluates the application of these models to selected internal flow configurations. The main conclusions of this analysis are: (1) the k-epsilon model is used in a majority of all the two-dimensional flow calculations reported in the literature; (2) modified forms of the k-epsilon model improve the performance for flows with streamline curvature and heat transfer; (3) for flows with swirl, the k-epsilon model performs rather poorly (the algebraic stress model performs better in this case); and (4) for flows with regions of secondary flow (noncircular duct flows), the algebraic stress model performs fairly well. Two important factors in the numerical solution of the model equations, namely false diffusion and inlet boundary conditions, are discussed. The existence of countergradient transport and its implications in turbulence modeling are examined. Finally, some recommendations for improving the model performance are made. The need for detailed experimental data in flows with strong curvature is emphasized.

Nallasamy, M.↗

A critical evaluation of various turbulence models as applied to internal fluid flows

Models employed in the computation of turbulent flows are described and their application to internal flows is evaluated by examining the predictions of various turbulence models in selected flow configurations. The main conclusions are: (1) the k-epsilon model is used in a majority of all the two-dimensional flow calculations reported in the literature; (2) modified forms of the k-epsilon model improve the performance for flows with streamline curvature and heat transfer; (3) for flows with swirl, the k-epsilon model performs rather poorly; the algebraic stress model performs better in this case; and (4) for flows with regions of secondary flow (noncircular duct flows), the algebraic stress model performs fairly well for fully developed flow, for developing flow, the algebraic stress model performance is not good; a Reynolds stress model should be used. False diffusion and inlet boundary conditions are discussed. Countergradient transport and its implications in turbulence modeling is mentioned. Two examples of recirculating flow predictions obtained using PHOENICS code are discussed. The vortex method, large eddy simulation (modeling of subgrid scale Reynolds stresses), and direct simulation, are considered. Some recommendations for improving the model performance are made. The need for detailed experimental data in flows with strong curvature is emphasized.

Nallasamy, M.↗

Pore-scale simulation of microbial growth using a genome-scale metabolic model: Implications for Darcy-scale reactive transport

Recent advances in microbiology have enabled the quantitative simulation of microbial metabolism and growth based on genome-scale characterization of metabolic pathways and fluxes. We have incorporated a genome-scale metabolic model of the iron-reducing bacteria Geobacter sulfurreducens into a pore-scale simulation of microbial growth based on coupling of iron reduction to oxidation of a soluble electron donor (acetate). In our model, fluid flow and solute transport is governed by a combination of the Navier-Stokes and advection-diffusion-reaction equations. Microbial growth occurs only on the surface of soil grains where solid-phase mineral iron oxides are available. Mass fluxes of chemical species associated with microbial growth are described by the genome-scale microbial model, implemented using a constraint-based metabolic model, and provide the Robin-type boundary condition for the advection-diffusion equation at soil grain surfaces. Conventional models of microbially-mediated subsurface reactions use a lumped reaction model that does not consider individual microbial reaction pathways, and describe reactions rates using empirically-derived rate formulations such as the Monod-type kinetics. We have used our pore-scale model to explore the relationship between genome-scale metabolic models and Monod-type formulations, and to assess the manifestation of pore-scale variability (microenvironments) in terms of apparent Darcy-scale microbial reaction rates. The genome-scale model predicted lower biomass yield, and different stoichiometry for iron consumption, in comparisonto prior Monod formulations based on energetics considerations. We were able to fit an equivalent Monod model, by modifying the reaction stoichiometry and biomass yield coefficient, that could effectively match results of the genome-scale simulation of microbial behaviors under excess nutrient conditions, but predictions of the fitted Monod model deviated from those of the genome-scale model under conditions in which one or more nutrients were limiting. The fitted Monod kinetic model was also applied at the Darcy scale; that is, to simulate average reaction processes at the scale of the entire pore-scale model domain. As we expected, even under excess nutrient conditions for which the Monod and genome-scale models predicted equal reaction rates at the pore scale, the Monod model over-predicted the rates of biomass growth and iron and acetate utilization when applied at the Darcy scale. This discrepancy is caused by an inherent assumption of perfect mixing over the Darcy-scale domain, which is clearly violated in the pore-scale models. These results help to explain the need to modify the flux constraint parameters in order to match observations in previous applications of the genome-scale model at larger scales. These results also motivate further investigation of quantitative multi-scale relationships between fundamental behavior at the pore scale (where genome-scale models are appropriately applied) and observed behavior at larger scales (where predictions of reactive transport phenomena are needed).

bioremediation↗

Numerical simulation of temperature field, microstructure evolution and mechanical properties of HSS during hot stamping

The hot stamping of boron steels is widely used to produce ultra high strength automobile components without any spring back. The ultra high strength of final products is attributed to the fully martensitic microstructure that is obtained through the simultaneous forming and quenching of the hot blanks after austenization. In the present study, a mathematical model incorporating both heat transfer and the transformation of austenite is presented. A FORTRAN program based on finite element technique has been developed which permits the temperature distribution and microstructure evolution of high strength steel during hot stamping process. Two empirical diffusion-dependent transformation models under isothermal conditions were employed respectively, and the prediction capability on mechanical properties of the models were compared with the hot stamping experiment of an automobile B-pillar part.

36 MATERIALS SCIENCE↗

Mathematical simulation of heat and mass transfer in a cylindrical channel versus the tangential momentum accommodation coefficient

The process of heat and mass transfer in a long cylindrical channel has been considered in terms of the mirror–diffuse model of the Maxwell boundary condition. The Hazen–Williams equation is used as a basic equation of the process kinetics. A constant temperature gradient is maintained in the channel. The heat and mass fluxes through the cross section of the channel versus the tangential momentum accommodation coefficient have been calculated in a wide range of the Knudsen number. The heat flux profiles have been constructed. A comparison with relevant published data has been carried out.

42 ENGINEERING↗

Purging of multilayer insulation by gas diffusion

An experimental investigation was conducted to determine the time required to purge a multilayer insulation (MLI) panel with gaseous helium by means of gas diffusion to obtain a condensable (nitrogen) gas concentration of less than 1 percent within the panel. Two flat, rectangular MLI panel configurations, one incorporating a butt joint, were tested. The insulation panels consisted of 15 double-aluminized Mylar radiation shields separated by double silk net spacers. The test results indicated that the rate which the condensable gas concentration at the edge or at the butt joint of an MLI panel was reduced was a significant factor in the total time required to reduce the condensable gas concentration within the panel to less than 1 percent. The experimental data agreed well with analytical predictions made by using a simple, one-dimensional gas diffusion model in which the boundary conditions at the edge of the MLI panel were time dependent.

Sumner, I. E.↗

Purging of a multilayer insulation with dacron tuft spacer by gas diffusion

The time and purge gas usage required to purge a multilayer insulation (MLI) panel with gaseous helium by means of gas diffusion to obtain a condensable gas (nitrogen) concentration of less than 1 percent within the panel are stipulated. Two different, flat, rectangular MLI panels, one incorporating a butt joint, were constructed of of 11 double-aluminized Mylar (DAM) radiation shields separated by Dacron tuft spacers. The DAM/Dacron tuft concept is known commercially as Superfloc. The nitrogen gas concentration as a function of time within the MLI panel could be adequately predicted by using a simple, one dimensional gas diffusion model in which the boundary conditions at the edge of the MLI panel were time dependent. The time and purge gas usage required to achieve 1 percent nitrogen gas concentration within the MLI panel varied from 208 to 86 minutes and 34.1 to 56.5 MLI panel purge volumes, respectively, for gaseous helium purge rates from 10 to 40 MLI panel volumes per hour.

Sumner, I. E.↗

Reliability of degrading dynamic systems subject to dynamic random loads

Reliability was determined for two degrading dynamic systems subject to random load processes. Damage is caused by loss of components for Daniels systems and crack extension for plates with cracks. The analysis accounted for the coupling between response and current damage state of the system. It is based on mean crossing rates of conditional processes and properties of diffusion models. Simple systems are used to illustrate proposed methods for estimating reliability.

Grigoriu, Mircea↗

Improved performance of a near-field thermophotovoltaic device by a back gapped reflector

Various spectral control techniques can be applied to improve the performance of a thermophotovoltaic (TPV) device. For example, a back surface reflector (BSR) is a common structure to improve the performance of TPV devices. A conventional metal BSR structure enhances the photogeneration rate by increasing the absorption probability of photons via back surface reflections, affording a second chance for absorption. However, the effects of surface passivation and external luminescence introduced by BSR structures have been previously ignored, which potentially decreases the performance of TPV devices. Recently, a back gapped reflector (BGR) structure was proposed to greatly improve the performance of far-field TPV devices by reducing imperfect reflections at the semiconductor-metal interface. In the present work, the performance improvement on a thin-film, near-field InAs TPV device with a BGR is investigated, comparing its performance to that with a conventional metal BSR. Surface passivation conditions are also investigated to further improve the performance of TPV devices with back reflectors. Further, the output power and efficiency are calculated using an iterative model combining fluctuational electrodynamics and the full drift diffusion model. For the well-passivated condition, when the BSR is replaced by the BGR, the calculated conversion efficiency of the near-field TPV was improved from 16.4% to 21% and the output power was increased by 10%. Finally, the absorption of the back reflectors and external luminescence loss are analyzed to explain the performance improvement.

external luminescence↗

Modeling of the SiC chemical vapor deposition process and comparison with experimental results

A model of the NASA Lewis SiC CVD process is described, and the results presented. A key feature of the model is the direct coupling of gas-phase chemical kinetics with diffusion. A deposition chemistry model using reactive sticking coefficients for each gas-phase species provides the diffusion boundary conditions. SiC deposition rates predicted by the model agree reasonably well with experimentally observed rates.

Annen, K. D.↗

A probability distribution model for rain rate

A systematic approach is suggested for modeling the probability distribution of rain rate. Rain rate, conditional on rain and averaged over a region, is modeled as a temporally homogeneous diffusion process with appropiate boundary conditions. The approach requires a drift coefficient-conditional average instantaneous rate of change of rain intensity-as well as a diffusion coefficient-the conditional average magnitude of the rate of growth and decay of rain rate about its drift. Under certain assumptions on the drift and diffusion coefficients compatible with rain rate, a new parametric family-containing the lognormal distribution-is obtained for the continuous part of the stationary limit probability distribution. The family is fitted to tropical rainfall from Darwin and Florida, and it is found that the lognormal distribution provides adequate fits as compared with other members of the family and also with the gamma distribution.

Kedem, Benjamin↗