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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Modeling gas–shell mixing in ICF with separated reactants

Mixing between fuel and shell materials in ICF implosions can affect implosion dynamics and even prevent ignition. We use data from a series of separated reactant experiments on the National Ignition Facility to calibrate and test the predictive power of gas–shell mix models. Two models are used to estimate fuel–shell mix: a Reynolds-averaged turbulence model and molecular diffusion. Minor uncertainties in capsule manufacture, experimental conditions, and values for mix model parameters produce significant variation in simulation results. Using input/output pairs from 1D simulations, we train Gaussian process surrogate models to predict experimental quantities of interest. The surrogates are used to construct posteriors for mix model parameters by marginalizing over uncertainties in capsule manufacture and experimental conditions. Mix models are calibrated with a subset of experimental data (neutron yields, ion temperature, and bang time) and tested using the remaining data. In general, both the diffusion and turbulence model correctly predict experimental DT and TT neutron yields. Despite having more free parameters, the turbulence model underpredicts ion temperature at high convergence ratio. Furthermore, the simpler diffusion model correctly predicts these temperatures, suggesting nonhydrodynamic gas–shell mix. The computational model consistently overpredicts DD neutron yield, indicating possible shortcomings outside of the mix model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Probabilistic Analysis of Long-Term Degradation of Microwave Cavity Flow Sensor

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Final Report on Development of Microwave Resonant Cavity Transducer for Fluid Flow Sensing: Development of Microwave Cavity Flow Meter for Advanced Reactor High Temperature Fluids

We are investigating a microwave resonant cavity transducer for flow sensing in the vessel of a high temperature fluid advanced reactor (AR), such as a molten salt cooled reactor (MSCR) or a sodium fast reactor (SFR). This transducer is a hollow metallic cylindrical cavity, with one of the flat walls of the cylinder flexible enough to undergo microscopic deflection due to dynamic fluid pressure. Membrane deflection leads to a shift in the microwave resonant frequency, which can be detected with a spectrum analyzer. We have developed a continuum electromechanics model of the microwave resonant cavity transducer performance. The deflection of the membrane is calculated using equations from the theory of plates and shells for deflection of radially constrained thin circular disk under uniform load. The microwave frequency shift is calculated using the equations from the microwave cavity perturbation theory. The transducer model was validated with proof-of-concept experimental data of water flow sensing with a Brass 360 right circular cylinder with 22.2mm diameter and 203µm thick wall. The cavity was excited through a subwavelength hole in TE 011 mode with resonant frequency f ≈ 17.8GHz. Subsequently, we have performed a preliminary proof-of-principle test of flow sensing in high temperature liquid sodium in environment. For this test, we have developed a cylindrical resonator with the same dimensions as for the Brass 360 cavity. The resonator for high temperature liquid sodium test was machined from stainless steel 316 and electroplated with silver on the interior surfaces. We have also developed and insertion probe consisting of a 50cm WR-42 brass waveguide enclosed in a protective SS316 tube. The cavity was excited through a subwavelength hole on the side of the wall of the cylinder in the TE 011 mode with resonant frequency f ≈ 17.8GHz. The liquid sodium setup consists of a cylindrical vessel with a center feed line, where a transducer inserted through the top cap of the vessel measures velocity of the impinging liquid jet. The sodium flow sensing study was performed in a liquid sodium vessel at 340°C temperature and ambient pressure. The flow rate was changed by varying sodium pump power, and a frequency shift of several hundred kHz was observed. The transducer was removed after 70 days of testing. No visual damage was observed, and the electromagnetic response remained the same. As a calibration experiment, we have assembled a water vessel with a center feed, and with dimensions similar to those of the liquid sodium setup. Water flow sensing was performed before and after liquid sodium test. The responses in both cases are similar and agree with COMSOL computer simulations. Because long-term (multi-year) experimental tests of transducer resilience to harsh environment are not practical, we have developed a probabilistic model of creep to estimate transducer resilience to the harsh environment. The probabilistic model considers diffusion creep under the condition of high temperature and low stress, where the stress and temperature are allowed to be random variables with Gaussian distributions. Using the probabilistic model, we estimate inelastic membrane deflections due to creep for several temperature ranges. We conclude that for temperatures less than 650°C, creep has negligible long-term effect on the transducer performance.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2015

Revisited reaction-diffusion model of thermal desorption spectroscopy experiments on hydrogen retention in material

Desorption phase of thermal desorption spectroscopy (TDS) experiments performed on tungsten samples exposed to flux of hydrogen isotopes in fusion relevant conditions is analyzed using a reaction-diffusion model describing hydrogen retention in material bulk. Two regimes of hydrogen desorption are identified depending on whether hydrogen trapping rate is faster than hydrogen diffusion rate in material during TDS experiments. In both regimes, a majority of hydrogen released from material defects is immediately outgassed instead of diffusing deeply in material bulk when the evolution of hydrogen concentration in material is quasi-static, which is the case during TDS experiments performed with tungsten samples exposed to flux of hydrogen isotopes in fusion related conditions. In this context, analytical expressions of the hydrogen outgassing flux as a function of the material temperature are obtained with sufficient accuracy to describe main features of thermal desorption spectra (TDSP). These expressions are then used to highlight how characteristic temperatures of TDSP depend on hydrogen retention parameters, such as trap concentration or activation energy of detrapping processes. The use of Arrhenius plots to characterize retention processes is then revisited when hydrogen trapping takes place during TDS experiments. Retention processes are also characterized using the shape of desorption peaks in TDSP, and it is shown that diffusion of hydrogen in material during TDS experiment can induce long desorption tails visible aside desorption peaks at high temperature in TDSP. These desorption tails can be used to estimate activation energy of diffusion of hydrogen in material.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Impact of cement composition, brine concentration, diffusion rate, reaction rate and boundary condition on self-sealing predictions for cement-CO 2 systems

Geological CO 2 storage (GCS) plays an important role in curbing CO 2 emissions by reducing the carbon footprint of difficult to decarbonize operations and achieve negative CO 2 emissions through activities like Bioenergy with Carbon Capture and Storage (BECCS) and Direct Air Carbon Capture and Storage (DACCS). Leakage of CO 2 through wells is an important concern when it comes to deployment of large-scale GCS. Existing wells in sites that are otherwise suitable for GCS can act as conduits for stored CO 2 to escape the reservoir. There is broad consensus that the main risk of leakage through wellbores is via fractures/damaged pathways. The results from several studies evaluating the permeability evolution of cement fractures in wells upon leakage of CO 2 agree that smaller fracture apertures, slower brine velocities and higher brine residence times promote self-sealing of fractures by mineral precipitation. Quantitatively, however, the differences in sealing conditions are significant and are typically attributed to differences in experimental conditions or model assumptions. Here we examine the sensitivity of our model, describing CO 2 leakage through wellbores, to cement composition, brine concentration, diffusion rates, and reaction rates. We also evaluate the impact of the boundary condition to allow comparisons between observations from experiments performed at constant flow rate and model predictions made at constant pressure conditions. Our results show that diffusion and reactions rates have the most impact on the self-sealing criteria for cement-CO 2 systems. In addition, conditions associated with self-sealing of fractures at constant flow rate require longer fractures, smaller fracture apertures, and slower velocities than under constant pressure.

58 GEOSCIENCES↗

DOE OSTI2024

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↗

DOE OSTI2022

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↗

DOE OSTI2020

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↗

DOE OSTI2022

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↗

DOE OSTI2021

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↗

DOE OSTI2024

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↗

DOE OSTI2013

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↗

DOE OSTI2013

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↗

DOE OSTI2017

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↗

DOE OSTI2021

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↗

DOE OSTI2016

Diffusion-driven growth of nanowires by low-temperature molecular beam epitaxy

With ZnTe as an example, we use two different methods to unravel the characteristics of the growth of nanowires (NWs) by gold-catalyzed molecular beam epitaxy at low temperature. In the first approach, CdTe insertions have been used as markers, and the nanowires have been characterized by scanning transmission electron microscopy, including geometrical phase analysis and energy dispersive electron spectrometry; the second approach uses scanning electron microscopy and the statistics of the relationship between the length of the tapered nanowires and their base diameter. Axial and radial growth are quantified using a diffusion-limited model adapted to the growth conditions; analytical expressions describe well the relationship between the NW length and the total molecular flux (taking into account the orientation of the effusion cells), and the catalyst-nanowire contact area. A long incubation time is observed. This analysis allows us to assess the evolution of the diffusion lengths on the substrate and along the nanowire sidewalls, as a function of temperature and deviation from stoichiometric flux.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Current-induced spin torque resonance of magnetic insulators affected by field-like spin-orbit torques and out-of-plane magnetizations

The spin-torque ferromagnetic resonance (ST-FMR) in a bilayer system consisting of a magnetic insulator such as Y{sub 3}Fe{sub 5}O{sub 12} and a normal metal with spin-orbit interaction such as Pt is addressed theoretically. We model the ST-FMR for all magnetization directions and in the presence of field-like spin-orbit torques based on the drift-diffusion spin model and quantum mechanical boundary conditions. ST-FMR experiments may expose crucial information about the spin-orbit coupling between currents and magnetization in the bilayers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Diffusion-limited hydrolysis in polymeric materials

Polymer degradation under aggressive environmental stressors often develops heterogeneities due to diffusion-limited reaction phenomena. This is well established for diffusion-limited oxidation (DLO), which is known to occur for most polymeric materials at elevated temperatures but has been less summarized for the conditions of diffusion-limited hydrolysis (DLH). Here, an overview of hydrolysis for several materials and a computational model, analogous to the underlying equations for DLO, able to define this diffusion-reaction system is presented. A systematic study of the influence of various parameters, such as water diffusivity, reaction rate and order, and a more in-depth focus on residual isocyanate hydrolysis in a polymeric methylene diphenyl diisocyanate (pMDI) based polyurethane (PU) foam is given. For this system, we present experimental data for model ‘input’ parameters and discuss predictions for different conditions. We conceptually compare the behavior of diffusion-limited oxidation to that of diffusion-limited hydrolysis (DLH). With the mathematical framework and key material properties presented herein, any DLH phenomena following Fickian diffusion behavior can be understood, modeled, and predicted.

36 MATERIALS SCIENCE↗