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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Mesoscale Modeling of the Effects of Accelerated Burnup on UO2 Microstructural Evolution

Accelerating the nuclear fuel qualification process will rely on some combination of advanced modeling and simulation techniques with accelerated irradiation testing and separate effects experiments to enable the development of new fuel concepts in a shorter time frame. One of the key challenges to successfully leveraging accelerated irradiation tests will be understanding the artifacts that may be introduced with accelerated accumulation of dose and/or burnup. This work presents phase field (MARMOT) simulations of the evolution of representative 2D UO2 microstructures up to 40 MWd/kgU. Simulations were performed under both commercial light water reactor fuel conditions as well as those that would be expected for highly accelerated (~10x) burnup conditions similar to those used in the MiniFuel irradiations in Oak Ridge National Laboratory’s High Flux Isotope Reactor. The phase field model was coupled with a discrete nucleation algorithm to model re- structuring at high burnup. The effect of the different fission rates in both microstructures was investigated at two temperatures: 650?C and 800?C. The lower temperature simulations both showed an onset of restructuring at nearly 60 MWd/kgU. More extensive restructuring was obtained in the MiniFuel microstructure compared with that of the PWR fuel. At 800?C, no restructuring was obtained as a result of the thermally activated diffusion of Xe atoms and U vacancies to fission gas bubbles, which reduces the nucleation driving force. These results highlight the importance of using modeling and simulation tools to inform the environmental conditions during targeted accelerated irradiation tests to extract the most useful fuel performance data.

accelerated fuel qualification, Phase Field, Restr↗
DOE OSTI2021

Radiation-Enhanced Anion Transport in Hematite

The influence of radiation-induced (1 MeV energy H+ to ~0.1 dpa at 450°C), non-equilibrium point defect populations on mass transport is studied with an integrated campaign of experimental and theoretical methods. Using epitaxial thin films of hematite with embedded 18O tracer layers and nanoscale atom probe tomography measurements, it is shown that anion self-diffusion is enhanced by at least 2 orders of magnitude under irradiation compared to thermal diffusion alone. Complementary scanning transmission electron microscopy of vacuum annealed specimens, reveals associated microstructural changes in the oxide films, including local phase transformation to Fe3O4 and the development of nanoscale voids from vacancy coalescence. Point defect formation and migration energies were computed from density functional theory and applied within the context of chemical rate theory to analyze contributions from both interstitial and vacancy mechanisms to self-diffusion in thermal and irradiation conditions. Comparisons are made between calculated, literature and newly measured self-diffusion values, revealing good agreement on the magnitude of radiation-enhanced anion diffusion. Further, the model suggests a transition from vacancy to interstitialcy mechanisms at low temperatures and high oxygen activity, consistent with the varied activation energies reported from prior studies.

iron oxides, isotopic tracers, atom probe tomograp↗
DOE OSTI2016

Inversion of geothermal heat flux in a thermomechanically coupled nonlinear Stokes ice sheet model

We address the inverse problem of inferring the basal geothermal heat flux from surface velocity observations using a steady-state thermomechanically coupled nonlinear Stokes ice flow model. This is a challenging inverse problem since the map from basal heat flux to surface velocity observables is indirect: the heat flux is a boundary condition for the thermal advection–diffusion equation, which couples to the nonlinear Stokes ice flow equations; together they determine the surface ice flow velocity. This multiphysics inverse problem is formulated as a nonlinear least-squares optimization problem with a cost functional that includes the data misfit between surface velocity observations and model predictions. A Tikhonov regularization term is added to render the problem well posed. We derive adjoint-based gradient and Hessian expressions for the resulting partial differential equation (PDE)-constrained optimization problem and propose an inexact Newton method for its solution. As a consequence of the Petrov–Galerkin discretization of the energy equation, we show that discretization and differentiation do not commute; that is, the order in which we discretize the cost functional and differentiate it affects the correctness of the gradient. Using two- and three-dimensional model problems, we study the prospects for and limitations of the inference of the geothermal heat flux field from surface velocity observations. The results show that the reconstruction improves as the noise level in the observations decreases and that short-wavelength variations in the geothermal heat flux are difficult to recover. We analyze the ill-posedness of the inverse problem as a function of the number of observations by examining the spectrum of the Hessian of the cost functional. Motivated by the popularity of operator-split or staggered solvers for forward multiphysics problems – i.e., those that drop two-way coupling terms to yield a one-way coupled forward Jacobian – we study the effect on the inversion of a one-way coupling of the adjoint energy and Stokes equations. Here, we show that taking such a one-way coupled approach for the adjoint equations can lead to an incorrect gradient and premature termination of optimization iterations. This is due to loss of a descent direction stemming from inconsistency of the gradient with the contours of the cost functional. Nevertheless, one may still obtain a reasonable approximate inverse solution particularly if important features of the reconstructed solution emerge early in optimization iterations, before the premature termination.

15 GEOTHERMAL ENERGY↗
DOE OSTI2018

An a priori Analysis of Mixture Fraction-Based Modeling of Coal Combustion

Due to the high computational cost associated with using detailed kinetic mechanisms, it is common practice to implement reaction models which rely on relatively small number of variables when simulating turbulent combustion systems. In this study, we perform an a priori assessment of the accuracy of three chemistry models: a Burke-Schumann chemistry model, an equilibrium chemistry model, and a steady laminar flamelet model. All chemistry modeling approaches considered include a mixture fraction for devolatilization products, and a normalized enthalpy deficit as parameters, with scalar dissipation rate as an additional input parameter for the flamelet model. Data used as a basis for the study is generated through One-Dimensional Turbulence (ODT) simulation of a coal combustion system in which high-fidelity models are used for gas phase chemistry and devolatilization. The flamelet model was found to be more accurate than both the equilibrium and Burke-Schumann models over a majority of parameter space. The performance of all models considered declined substantially for sufficiently fuel-rich conditions. Furthermore, the analysis in this work indicate that both differential diffusion and speciation of volatiles have a perceptible negative impact on flamelet reconstructions, though the effect of differential diffusion was more substantial.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Micro-continuum approach for mineral precipitation

Abstract Rates and extents of mineral precipitation in porous media are difficult to predict, in part because laboratory experiments are problematic. It is similarly challenging to implement numerical methods that model this process due to the need to dynamically evolve the interface of solid material. We developed a multiphase solver that implements a micro-continuum simulation approach based on the Darcy–Brinkman–Stokes equation to study mineral precipitation. We used the volume-of-fluid technique in sharp interface implementation to capture the propagation of the solid mineral surface. Additionally, we utilize an adaptive mesh refinement method to improve the resolution of near interface simulation domain dynamically. The developed solver was validated against both analytical solution and Arbitrary Lagrangian–Eulerian approach to ensure its accuracy on simulating the propagation of the solid interface. The precipitation of barite (BaSO 4 ) was chosen as a model system to test the solver using variety of simulation parameters: different geometrical constraints, flow conditions, reaction rate and ion diffusion. The growth of a single barite crystal was simulated to demonstrate the solver’s capability to capture the crystal face specific directional growth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Small to Full-Height Scale Comparisons of Cesium Ion Exchange Performance with Crystalline Silicotitanate

Crystalline silicotitanate (CST) ion exchanger is planned to be used to remove cesium (137Cs) from the aqueous phases of Hanford tank wastes in preparation for vitrification at the Waste Treatment and Immobilization Plant (WTP) Pretreatment Facility. Column scale up testing was conducted to evaluate performance of Cs exchange onto the CST. Batch contact testing was conducted to assess exchange kinetics of Cs exchange at four different CST PSDs. Column testing was conducted at three column sizes, small (2.5% full height), medium (12% full height), and full height, to assess Cs load performance behavior. Testing at the small scale was compared to actual Hanford tank waste testing to verify the validity of simulant tests to accurately represent full height column performance. A change in CST particle size was essential in scaling the small column dynamics up to full scale. Based on these results, a determination of intraparticle and film diffusion impacts on overall mass transfer coefficients will allow future modeling of the breakthrough performance at a range of process conditions.

Westesen, Amy M.↗
DOE OSTI2017

Antimony diffusion in CdTe

Group V dopants may be used for next-generation high-voltage cadmium telluride (CdTe) solar photovoltaics, but fundamental defect energetics and kinetics need to be understood. Here, antimony (Sb) diffusion is studied in single-crystal and polycrystalline CdTe under Cd-rich conditions. Diffusion profiles are determined by dynamic secondary ion mass spectroscopy and analyzed with analytical bulk and grain-boundary diffusion models. Slow bulk and fast grain-boundary diffusion are found. Density functional theory is used to understand formation energy and mechanisms. Lastly, the theory and experimental results create new understanding of group V defect kinetics in CdTe.

14 SOLAR ENERGY↗
DOE OSTI2015

Challenges of accelerated aging techniques for elastomer lifetime predictions

Elastomers are often degraded when exposed to air or high humidity for extended times (years to decades). Lifetime estimates normally involve extrapolating accelerated aging results made at higher than ambient environments. Several potential problems associated with such studies are reviewed, and experimental and theoretical methods to address them are provided. The importance of verifying time–temperature superposition of degradation data is emphasized as evidence that the overall nature of the degradation process remains unchanged versus acceleration temperature. The confounding effects that occur when diffusion-limited oxidation (DLO) contributes under accelerated conditions are described, and it is shown that the DLO magnitude can be modeled by measurements or estimates of the oxygen permeability coefficient (P Ox ) and oxygen consumption rate (Φ). P Ox and Φ measurements can be influenced by DLO, and it is demonstrated how confident values can be derived. In addition, several experimental profiling techniques that screen for DLO effects are discussed. Values of Φ taken from high temperature to temperatures approaching ambient can be used to more confidently extrapolate accelerated aging results for air-aged materials, and many studies now show that Arrhenius extrapolations bend to lower activation energies as aging temperatures are lowered. Furthermore, best approaches for accelerated aging extrapolations of humidity-exposed materials are also offered.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Simulations of self- and Xe diffusivity in uranium mononitride including chemistry and irradiation effects

A combination of density functional theory and empirical potential atomic scale simulations have been used to determine a model for defect stability and mobility in uranium mononitride (UN), as a function of temperature (T) and N 2 partial pressure (p N 2 ). Using the model, predictions of hypo-stoichiometry under U-rich conditions compare favorably to CALPHAD calculations using the TAF-ID database. Furthermore, our predictions of U and N self-diffusivity are in good agreement with experiments carried out as a function of T at specific partial pressures under thermal equilibrium. The validated atomic scale data have then been implemented within a cluster dynamics method to simulate irradiation-enhanced defect concentrations. All defects and clusters studied have significantly enhanced concentrations, with respect to thermal equilibrium, as T is lowered. The irradiation-enhanced Xe diffusivity is compared to post-irradiation annealing and in-pile experiments. In conclusion, the contributions of various defects and clusters to non-stoichiometry, self-diffusivity, and Xe diffusivity are discussed.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2016

Sorption Modeling and Verification for Off-Gas Treatment

This project was successfully executed to provide valuable adsorption data and improve a comprehensive model developed in previous work by the authors. Data obtained were used in an integrated computer program to predict the behavior of adsorption columns. The model is supported by experimental data and has been shown to predict capture of off gas similar to that evolving during the reprocessing of nuclear waste. The computer program structure contains (a) equilibrium models of off-gases with the adsorbate; (b) mass-transfer models to describe off-gas mass transfer to a particle, diffusion through the pores of the particle, and adsorption on the active sites of the particle; and (c) incorporation of these models into fixed bed adsorption modeling, which includes advection through the bed. These models are being connected with the MOOSE (Multiphysics Object-Oriented Simulation Environment) software developed at the Idaho National Laboratory through DGOSPREY (Discontinuous Galerkin Off-gas SeParation and REcoverY) computer codes developed in this project. Experiments for iodine and water adsorption have been conducted on reduced silver mordenite (Ag0Z) for single layered particles. Adsorption apparatuses have been constructed to execute these experiments over a useful range of conditions for temperatures ranging from ambient to 250°C and water dew points ranging from -69 to 19°C. Experimental results were analyzed to determine mass transfer and diffusion of these gases into the particles and to determine which models best describe the single and binary component mass transfer and diffusion processes. The experimental results were also used to demonstrate the capabilities of the comprehensive models developed to predict single-particle adsorption and transients of the adsorption-desorption processes in fixed beds. Models for adsorption and mass transfer have been developed to mathematically describe adsorption kinetics and transport via diffusion and advection processes. These models were built on a numerical framework for solving conservation law problems in one-dimensional geometries such as spheres, cylinders, and lines. Coupled with the framework are specific models for adsorption in commercial adsorbents, such as zeolites and mordenites. Utilizing this modeling approach, the authors were able to accurately describe and predict adsorption kinetic data obtained from experiments at a variety of different temperatures and gas phase concentrations. A demonstration of how these models, and framework, can be used to simulate adsorption in fixed- bed columns is provided. The CO 2 absorption work involved modeling with supportive experimental information. A dynamic model was developed to simulate CO 2 absorption using high alkaline content water solutions. The model is based upon transient mass and energy balances for chemical species commonly present in CO 2 absorption. A computer code was developed to implement CO 2 absorption with a chemical reaction model. Experiments were conducted in a laboratory scale column to determine the model parameters. The influence of geometric parameters and operating variables on CO 2 absorption was studied over a wide range of conditions. Continuing work could employ the model to control column operation and predict the absorption behavior under various input conditions and other prescribed experimental perturbations. The value of the validated models and numerical frameworks developed in this project is that they can be used to predict the sorption behavior of off-gas evolved during the reprocessing of nuclear waste and thus reduce the cost of the experiments. They can also be used to design sorption processes based on concentration limits and flow-rates determined at the plant level.

42 ENGINEERING↗
DOE OSTI2021

Particle Size Distribution Dynamics Can Help Constrain the Phase State of Secondary Organic Aerosol

Particle phase state is a property of atmospheric aerosols that has important implications for the formation, evolution, and gas/particle partitioning of secondary organic aerosol (SOA). In this work, we use a size-resolved chemistry and microphysics model (SOM-TOMAS), updated to include an explicit treatment of particle phase state, to constrain the bulk diffusion coefficient (Db) of SOA produced from ??- pinene ozonolysis. By leveraging data from laboratory experiments performed in the absence of a seed and under dry conditions, we find that the Db for SOA can be constrained (1-7 ×10-15 cm2 s-1 in these experiments) by simultaneously reproducing the time-varying SOA mass concentrations and the evolution of the particle size distribution. Another version of our model that used the predicted SOA composition to calculate the glass transition temperature, viscosity, and, ultimately, Db (~10-15 cm2 s-1) of the SOA was able to reproduce the mass and size distribution measurements when we included oligomer formation (oligomers accounted for about a fifth of the SOA mass). Our work highlights the potential of a size resolved SOA model to constrain the particle phase state of SOA by utilizing historical measurements of the evolution of the particle size distribution.

He, Yicong↗
DOE OSTI2022

Modeling Ni redistribution in the hydrogen electrode of solid oxide cells through Ni(OH) 2 diffusion and Ni-YSZ wettability change

Ni redistribution in the hydrogen electrodes of solid oxide cells is an important degradation mechanism. Its driving force is still under debate. This work focuses on the interplay between Ni(OH) 2 diffusion and Ni-YSZ wettability change. With a Ni(OH) 2 diffusion model developed in our previous work, we further employ three models for the Ni-YSZ contact angle to quantify the Ni-YSZ wettability change. The microstructure evolutions in a reconstructed Ni-YSZ electrode are simulated under selected experimental operating conditions. It is shown that the phenomenological model can capture the driving force of Ni spreading/detachment on YSZ surface and the Ni migration reported in experiments can be qualitatively reproduced through the competition between Ni(OH) 2 diffusion and Ni-YSZ wettability change. It is also found that the initial microstructure of the active layer and the microstructure characteristics in the support layer can strongly affect the distribution of steam partial pressure and overpotential, hence the Ni redistribution, in the active layer, which may explain the inconsistencies in experiments. Further, our results show that there are two missing pieces in current theory of Ni redistribution: a mechanism of fast Ni diffusion under humid atmosphere and the physical reason behind the Ni spreading/detachment.

25 ENERGY STORAGE↗
DOE OSTI2016

Activation volumes of oxygen self-diffusion in fluorite structured oxides

In this study, fluorite structured oxides are used in numerous applications and as such it is necessary to determine their materials properties over a range of conditions. In the present study we employ molecular dynamics calculations to calculate the elastic and expansivity data, which are then used in a thermodynamic model (the cBΩ model) to calculate the activation volumes of oxygen self-diffusion coefficient in ThO 2 , UO 2 and PuO 2 fluorite structured oxides over a wide temperature range. We present relations to calculate the activation volumes of oxygen self-diffusion coefficient in ThO 2 , UO 2 and PuO 2 for a wide range of temperature (300–1700 K) and pressure (–7.5 to 7.5 GPa).

36 MATERIALS SCIENCE↗
DOE OSTI2014

Exponential time differencing methods with Chebyshev collocation for polymers confined by interacting surfaces

We present a fast and accurate numerical method for the self-consistent field theory calculations of confined polymer systems. It introduces an exponential time differencing method (ETDRK4) based on Chebyshev collocation, which exhibits fourth-order accuracy in temporal domain and spectral accuracy in spatial domain, to solve the modified diffusion equations. Similar to the approach proposed by Hur et al. [Macromolecules 45, 2905 (2012)], non-periodic boundary conditions are adopted to model the confining walls with or without preferential interactions with polymer species, avoiding the use of surface field terms and the mask technique in a conventional approach. The performance of ETDRK4 is examined in comparison with the operator splitting methods with either Fourier collocation or Chebyshev collocation. Numerical experiments show that our exponential time differencing method is more efficient than the operator splitting methods in high accuracy calculations. This method has been applied to diblock copolymers confined by two parallel flat surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Controlling kinetic and diffusive length-scales during absorptive hydrogen removal in methane dehydroaromatization on MoC x /H-ZSM-5 catalysts

Addition of Zr metal absorbent to MoC x /H-ZSM-5 in the form of staged-bed, stratified-bed, and interpellet physical mixtures effectively scavenges H 2 from catalyst proximity, enhancing maximum single-pass benzene + naphthalene yield during methane dehydroaromatization (DHA) reactions to 14-16% compared to 8% in formulations without zirconium. The coupling of spatially-distinct catalytic and absorptive functions is achieved by dispersive/diffusive transport which conveys H 2 to staged Zr both co- and counter-current to bulk advection, thereby suppressing axial H 2 partial pressure profiles along the catalyst bed and enhancing net aromatization rates. We evince hitherto unreported significance of dispersive hydrogen transport during methane DHA by measurement of Péclet number, Pe = 1.32, in H 2 tracer studies with step-change or impulse input to inert catalyst proxies. Kinetic limits to methane pyrolysis are quantified by Damköhler number, Da, for synthesis of benzene, Da B = 0.15, and naphthalene, Da N = 0.03, determined from kinetic studies which rigorously account for reversibility of DHA reactions. Here, detailed reaction-transport models synthesize interplay of kinetic, diffusive, and convective length-scales captured by Péclet and Damköhler number to predict influence of catalyst-absorbent proximity and process flow-conditions on aromatization rates. Systematic control of catalyst bed-length, L, or linear flow velocity, u, predictably alters Pe and Da to effect improvements in methane conversion with and without Zr metal, corroborating results from simulation of the reaction-transport model.

03 NATURAL GAS↗
DOE OSTI2014

FY 2014 Sodium-Supercritical CO 2 Interactions in the SNAKE Experiment Facility

One appealing feature of the supercritical carbon dioxide Brayton cycle energy conversion system is the small footprint that the hardware requires, which is in part due to the use of compact diffusion-bonded heat exchangers such as Printed Circuit Heat Exchanger (PCHEs) or Hybrid Heat Exchangers (H 2 X) as the heat source heat exchanger (sodium-to-CO 2 ) as well as the recuperator and cooler modules. Although PCHEs and H 2 Xs have a high degree of structural integrity, the potential for leaks to develop between the sodium and CO 2 coolant channels in the heat exchanger cannot be ruled out, and this would lead to discharge of high pressure CO 2 into the secondary coolant circuit. Due to the robustness of the PCHE design, catastrophic failure leading to CO 2 jet blowdown into the secondary sodium loop is not deemed likely. Rather, small cracks (or micro-leaks) may develop in which CO 2 will bleed into the secondary system at a relatively low rate and chemically react with the sodium. The goal of the sodium-CO 2 interaction tests is to gain a fundamental understanding of sodium-CO 2 interactions under prototypical conditions of compact diffusion-bonded heat exchanger failure, a fundamental understanding of self-plugging if it occurs, and the development of one-dimensional phenomenological models for the interactions between high-pressure CO 2 issuing into liquid sodium from a micro-leak across a stainless steel pressure boundary. These models will be validated using experiment data. Therefore, an experiment program at Argonne was initiated in Fiscal Year 2010 to investigate the reaction behavior between sodium and CO 2 under micro-leak conditions. Several reports have described the facility scaling rationale and design. Assembly of the SNAKE (S-CO 2 , Na Kinetics Experiment) began in Fiscal Year 2011 and was completed in July 2012. The sodium inventory in this facility consists of approximately 44 lbs (20 kg; ~21 L/5.5 gal).

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2013

Description of the First Observed Sodium-C0 2 Reactions in the Sodium CO 2 Interaction Experiment (SNAKE)

One appealing feature of the supercritical carbon dioxide Brayton cycle energy conversion system is the small footprint that the hardware requires, which is in part due to the use of Printed Circuit Heat Exchangers (PCHEs) as the heat source heat exchanger (sodium-to-CO 2 ) as well as the recuperator and cooler modules. Although PCHEs have a high degree of structural integrity, the potential for leaks to develop between the sodium and CO 2 coolant channels in the secondary heat exchanger cannot be ruled out, and this would lead to discharge of high pressure CO 2 into the secondary coolant circuit. Due to the robustness of the PCHE design, catastrophic failure leading to CO 2 jet blowdown into the secondary sodium loop is not deemed likely. Rather, small cracks (or micro-leaks) may develop in which CO 2 will bleed into the secondary system at a relatively low rate and chemically react with the sodium. The goal of the sodium-CO 2 interaction tests is to gain a fundamental understanding of sodium-CO 2 interactions under prototypical conditions of compact diffusion-bonded heat exchanger failure, a fundamental understanding of self-plugging if it occurs, and the development of one-dimensional phenomenological models for the interactions between high-pressure CO 2 issuing into liquid sodium from a micro-leak across a stainless steel pressure boundary. These models will be validated using experiment data. Therefore, an experiment program at Argonne was initiated in Fiscal Year 2010 to investigate the reaction characteristics between sodium and CO 2 under micro-leak conditions. Several reports have described the facility scaling rationale and design. Assembly of the SNAKE (S-CO 2 , Na Kinetics Experiment) began in Fiscal Year 2011 and was completed in July 2012. Approximately 44 lbs (20 kg; ~21 L/5.5 gal) of sodium was transferred into the SNAKE dump tank from a drum of clean sodium in July 2012.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

Laboratory Experiments on Bentonite Samples: FY16 Progress

The primary goal of this study is to improve the understanding of U(VI) sorption and diffusion behavior in sodium-montmorillonite in order to support the development of realistic conceptual models describing these processes in performance assessment models while (1) accounting for potential changes in system conditions over time and space, (2) avoiding overly conservative transport predictions, and (3) using a minimum number of fitting parameters.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗