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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Examination of Legacy Metallic Fuel Pins (U-10Zr) Tested in FFTF

The MFF series of metallic fuel (U-10wt%Zr) tests performed in the Fast Fuel Test Facility (FFTF) were the beginning tests to qualify the fuel as a driver fuel for FFTF. They all performed very well, to relatively high burnup and with no pin breaches. Tests MFF-3 and MFF-5 were chosen to be destructively examined because they were run at high peak cladding temperatures, 643 and 649°C respectively, and to modest to high burnup (138 and 101 MWd/kgM respectively). They were the only sodium fast reactor (SFR) metallic tests, with long (91.4 cm) fuel columns and clad in HT9, to be operated at these high temperatures. Extensive operating condition reconstruction was performed to accompany the examination, providing detailed operating conditions axially along the pins. In addition, detailed ORIGEN calculations provided fuel burnup along the pin length. These calculations and detailed cross-section metallography allowed a PhD dissertation to be performed where the fuel/cladding chemical interaction (FCCI) was modeled using Fickian and Soret Effect driven diffusion of rare earth fission products to the cladding surface and interact with the cladding. The examination of the MFF-3 and MFF-5 pins included neutron radiography, axial profilometry including pin bow and length measurements, precision gamma scanning, plenum gas puncturing to measure fission gas release, detailed chemical/isotopic analyses of fuel samples to confirm burnup calculations, and metallography of pin cross-sections, including micro-hardness measurements of fuel and cladding.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Elucidating Hydrogen Reaction-Induced Water Desorption from Oxide-Passivated Metal Surfaces for Plasma Applications

Elucidating the mechanisms responsible for sub-microsecond desorption of water and other impurities from electrode surfaces at high heating rates is crucial for understanding pulsed power behavior. Ionization of desorbed impurities in the vacuum regions causes power or current loss; devising methods to limit such desorption during the short time scale of pulsed power is needed to improve corresponding applications. Previous molecular modeling studies have strongly suggested that, under high vacuum conditions, the amount of water impurity desorbing from oxide surfaces on metal electrodes is at a sub-monolayer level at room temperature, which appears insufficient to explain observed pulsed power energy losses at high current densities. In this work, we apply Density Functional Theory (DFT) techniques to show that hydrogen trapped inside iron metal can diffuse into hematite (α-Fe 2 O 3 ) on the metal surface, ultimately reacting with the oxide to form Fe(II) and H 2 O. The latter desorbs at elevated temperature and may explain the anomalous amount of desorbed impurity inferred from pulsed-power experiments. We also apply a suite of characterization techniques to demonstrate that when iron metal is heated to 650 °C, the dominant surface oxide component becomes α-Fe 2 O 3 . The oxide facets exposed are found to be a mixture of (0001), (10-10), and others, in agreement with the DFT models used.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Molecular Dynamics Study of the Effect of Grafting Density on Ion Diffusivity in a MARTINI Coarse-Grained Strong Polyelectrolyte Brush

Because surface-grafted polyelectrolyte brushes (PEBs) are responsive to external stimuli, such as electric fields and ionic strength, PEBs are attractive for applications ranging from drug delivery to separation technologies. Essential to PEB utilization is understanding how critical parameters like grafting density (σ) impact the PEB structure and the dynamics of the PEB and counterions. To study the effect of σ on PEB and the counterion structure and dynamics, we fine-tune a coarse-grained model that retains the chemical specificity of a strong polyelectrolyte, poly[(2-(methacryloyloxy)ethyl) trimethylammonium chloride] (PMETAC), using the MARTINI force field. Using “salt-free” conditions where the counterion concentration balances the charge on the brush, we build coarse-grained (CG) molecular dynamics simulations for MARTINI PMETAC brushes (N = 150 monomers; M W = 31.2 kg/mol) at experimentally relevant values of σ = 0.05, 0.10, 0.20, and 0.40 chains/nm 2 . Using 5 μs simulations, we investigate the effects of grafting density on the PEB structure, ion dissociation dynamics, polymer mobility, and counterion diffusivity. Results show that competition between electrostatic interactions, steric hindrance, and polymer mobility controls counterion diffusivity. Finally, the interplay of these factors leads to diffusivity that depends non-monotonically on σ, with counterion diffusivity peaking at an intermediate σ = 0.10 chains/nm 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2016

cDPD: A new dissipative particle dynamics method for modeling electrokinetic phenomena at the mesoscale

In this work, we develop a “charged” dissipative particle dynamics (cDPD) model for simulating mesoscopic electrokinetic phenomena governed by the stochastic Poisson-Nernst-Planck and the Navier-Stokes equations. In particular, the transport equations of ionic species are incorporated into the DPD framework by introducing extra degrees of freedom and corresponding evolution equations associated with each DPD particle. Diffusion of ionic species driven by the ionic concentration gradient, electrostatic potential gradient, and thermal fluctuations is captured accurately via pairwise fluxes between DPD particles. The electrostatic potential is obtained by solving the Poisson equation on the moving DPD particles iteratively at each time step. For charged surfaces in bounded systems, an effective boundary treatment methodology is developed for imposing both the correct hydrodynamic and electrokinetics boundary conditions in cDPD simulations. To validate the proposed cDPD model and the corresponding boundary conditions, we first study the electrostatic structure in the vicinity of a charged solid surface, i.e., we perform cDPD simulations of the electrostatic double layer and show that our results are in good agreement with the well-known mean-field theoretical solutions. Furthermore, we simulate the electrostatic structure and capacity densities between charged parallel plates in salt solutions with different salt concentrations. Moreover, we employ the proposed methodology to study the electro-osmotic and electro-osmotic/pressure-driven flows in a micro-channel. In the latter case, we simulate the dilute poly-electrolyte solution drifting by electro-osmotic flow in a micro-channel, hence demonstrating the flexibility and capability of this method in studying complex fluids with electrostatic interactions at the micro- and nano-scales.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Uncovering the predictive pathways of lithium and sodium interchange in layered oxides

Ion exchange is a powerful method to access metastable materials with advanced functionalities for energy storage applications. However, high concentrations and unfavourably large excesses of lithium are always used for synthesizing lithium cathodes from parent sodium material, and the reaction pathways remain elusive. Here, using layered oxides as model materials, we demonstrate that vacancy level and its corresponding lithium preference are critical in determining the accessible and inaccessible ion exchange pathways. Taking advantage of the strong lithium preference at the right vacancy level, we establish predictive compositional and structural evolution at extremely dilute and low excess lithium based on the phase equilibrium between Li 0.94 CoO 2 and Na 0.48 CoO 2 . Such phase separation behaviour is general in both surface reaction-limited and diffusion-limited exchange conditions and is accomplished with the charge redistribution on transition metals. Guided by this understanding, we demonstrate the synthesis of Na y CoO 2 from the parent Li x CoO 2 and the synthesis of Li 0.94 CoO 2 from Na y CoO 2 at 1-1,000 Li/Na (molar ratio) with an electrochemical assisted ion exchange method by mitigating the kinetic barriers. Our study opens new opportunities for ion exchange in predictive synthesis and separation applications. Ion exchange is a powerful method to access metastable materials for energy storage, but identifying lithium and sodium interchange in layered oxides remains challenging. Using such model materials, vacancy level and corresponding lithium preference are shown to be crucial for ion exchange pathway accessibility.

25 ENERGY STORAGE↗
DOE OSTI2019

Realizability-preserving DG-IMEX method for the two-moment model of fermion transport

Building on the framework of Zhang & Shu [1], [2], we develop a realizability-preserving method to simulate the transport of particles (fermions) through a background material using a two-moment model that evolves the angular moments of a phase space distribution function f . The two-moment model is closed using algebraic moment closures; e.g., as proposed by Cernohorsky & Bludman [3] and Banach & Larecki [4]. Variations of this model have recently been used to simulate neutrino transport in nuclear astrophysics applications, including core-collapse supernovae and compact binary mergers. We employ the discontinuous Galerkin (DG) method for spatial discretization (in part to capture the asymptotic diffusion limit of the model) combined with implicit-explicit (IMEX) time integration to stably bypass short timescales induced by frequent interactions between particles and the background. Appropriate care is taken to ensure the method preserves strict algebraic bounds on the evolved moments (particle density and flux) as dictated by Pauli's exclusion principle, which demands a bounded distribution function (i.e., f ϵ [0,1]). This realizability-preserving scheme combines a suitable CFL condition, a realizability-enforcing limiter, a closure procedure based on Fermi-Dirac statistics, and an IMEX scheme whose stages can be written as a convex combination of forward Euler steps combined with a backward Euler step. The IMEX scheme is formally only first-order accurate, but works well in the diffusion limit, and — without interactions with the background — reduces to the optimal second-order strong stability-preserving explicit Runge-Kutta scheme of Shu & Osher [5]. Numerical results demonstrate the realizability-preserving properties of the scheme. As a result, we also demonstrate that the use of algebraic moment closures not based on Fermi-Dirac statistics can lead to unphysical moments in the context of fermion transport.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

A mechanism of defect-enhanced phase transformation kinetics in lithium iron phosphate olivine

Abstract Antisite defects are a type of point defect ubiquitously present in intercalation compounds for energy storage applications. While they are often considered a deleterious feature, here we elucidate a mechanism of antisite defects enhancing lithium intercalation kinetics in LiFePO 4 by accelerating the FePO 4 → LiFePO 4 phase transformation. Although Fe Li antisites block Li movement along the [010] migration channels in LiFePO 4 , phase-field modeling reveals that their ability to enhance Li diffusion in other directions significantly increases the active surface area for Li intercalation in the surface-reaction-limited kinetic regime, which results in order-of-magnitude improvement in the phase transformation rate compared to defect-free particles. Antisite defects also promote a more uniform reaction flux on (010) surface and prevent the formation of current hotspots under galvanostatic (dis)charging conditions. We analyze the scaling relation between the phase boundary speed, Li diffusivity and particle dimensions and derive the criteria for the co-optimization of defect content and particle geometry. A surprising prediction is that (100)-oriented LiFePO 4 plates could potentially deliver better performance than (010)-oriented plates when the Li intercalation process is surface-reaction-limited. Our work suggests tailoring antisite defects as a general strategy to improve the rate performance of phase-changing battery compounds with strong diffusion anisotropy.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Adaptive Interface-PINNs (AdaI-PINNs) for transient diffusion: Applications to forward and inverse problems in heterogeneous media

We model transient diffusion in heterogeneous materials using a novel physics-informed neural networks framework (PINNs) termed Adaptive interface physics-informed neural networks or AdaI-PINNs (Roy et al. arXiv preprint arXiv:2406.04626, 2024). AdaI-PINNs utilize different activation functions with trainable slopes tailored to each material region within the computational domain, allowing for a fully automated and adaptive PINNs approach to model interface problems with strongly and weakly discontinuous solutions. To enhance its performance in highly heterogeneous transient diffusion systems, we prescribe a suite of robust practices, including appropriate non-dimensionalization of equations, a biased sampling method, Glorot initialization, and the hard enforcement of boundary and initial conditions. Here we evaluate the efficacy of the proposed method on several benchmark forward and inverse problems. Comparative studies on one-dimensional and two-dimensional benchmark problems reveal that the modified AdaI-PINNs outperform its unmodified counterpart, achieving root-mean-square errors that are at least two orders of magnitude better in forward problems. For inverse problems, the maximum errors in the approximated diffusion coefficients by modified AdaI-PINNs are four orders of magnitude better than those of the unmodified version. Additionally, modified AdaI-PINNs demonstrate improved stability in problems with large material mismatches.

42 ENGINEERING↗
DOE OSTI2016

Determination of krypton diffusion coefficients in uranium dioxide using atomic scale calculations

We present a study of the diffusion of krypton in UO 2 using atomic scale calculations combined with diffusion models adapted to the system studied. The migration barriers of the elementary mechanisms for interstitial or vacancy assisted migration are calculated in the DFT + U framework using the nudged elastic band method. The attempt frequencies are obtained from the phonon modes of the defect at the initial and saddle points using empirical potential methods. The diffusion coefficients of Kr in UO 2 are then calculated by combining this data with diffusion models accounting for the concentration of vacancies and the interaction of vacancies with Kr atoms. We determined the preferred mechanism for Kr migration and the corresponding diffusion coefficient as a function of the oxygen chemical potential μ O or nonstoichiometry. For very hypostoichiometric (or U-rich) conditions, the most favorable mechanism is interstitial migration. For hypostoichiometric UO 2 , migration is assisted by the bound Schottky defect and the charged uranium vacancy, V U 4– . Around stoichiometry, migration assisted by the charged uranium–oxygen divacancy (V UO 2– ) and V U 4– is the favored mechanism. Finally, for hyperstoichiometric or O-rich conditions, the migration assisted by two V U 4– dominates. Kr migration is enhanced at higher μ O , and in this regime, the activation energy will be between 4.09 and 0.73 eV depending on nonstoichiometry. The experimental values available are in the latter interval. Since it is very probable that these values were obtained for at least slightly hyperstoichiometric samples, our activation energies are consistent with the experimental data, even if further experiments with precisely controlled stoichiometry are needed to confirm these results. Finally, the mechanisms and trends with nonstoichiometry established for Kr are similar to those found in previous studies of Xe.

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

Quantifying Prediction Fidelity in Multiscale Multiphysics Simulations (Final Technical Report)

This is a collaborative proposal that aims to establish theoretical foundations and computational tools that enable uncertainty quantification (UQ) in tightly coupled atomistic-to-continuum multiscale simulations. The program emphasizes the following three research thrusts: 1. UQ and its propagation in atomistic simulations, whether through intrusive or nonintrusive approaches; 2. Extraction of macroscale observables from atomistic simulations and propagation across scales; and 3. Uncertainty quantification and propagation in continuum simulations for macroscale properties tightly coupled with instantaneous states of the atomistic systems. Thus, the project offers to enable the use of multiscale multiphysics simulations as predictive design tools for complex systems.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Evaluation of Engineered Barrier Systems FY20 Report

This report describes research and development (R&D) activities conducted during fiscal year 2020 (FY20) specifically related to the Engineered Barrier System (EBS) R&D Work Package in the Spent Fuel and Waste Science and Technology (SFWST) Campaign supported by the United States (U.S.) Department of Energy (DOE). The R&D activities focus on understanding EBS component evolution and interactions within the EBS, as well as interactions between the host media and the EBS. A primary goal is to advance the development of process models that can be implemented directly within the Generic Disposal System Analysis (GDSA) platform or that can contribute to the safety case in some manner such as building confidence, providing further insight into the processes being modeled, establishing better constraints on barrier performance, etc. The FY20 EBS activities involved not only modeling and analysis work, but experimental work as well. Despite delays to some planned activities due to COVID-19 precautions, progress was made during FY20 in multiple research areas and documented in this report as follows: (1) EBS Task Force: Task 9/FEBEX Modeling Final Report: Thermo-Hydrological Modeling with PFLOTRAN, (2) preliminary sensitivity analysis for the FEBEX in-situ heater test, (3) cement-carbonate rock interaction under saturated conditions: from laboratory to modeling, (4) hydrothermal experiments, (5) progress on investigating the high temperature behavior of the uranyl-carbonate complexes, (6) in-situ and electrochemical work for model validation, (7) investigation of the impact of high temperature on EBS bentonite with THMC modeling, (8) sorption and diffusion experiments on bentonite, (9) chemical controls on montmorillonite structure and swelling pressure, (10) microscopic origins of coupled transport processes in bentonite, (11) understanding the THMC evolution of bentonite in FEBEX-DP—coupled THMC modeling, (12) modeling in support of HotBENT, an experiment studying the effects of high temperatures on clay buffers/near-field, and (13) high temperature heating and hydration column test on bentonite.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Investigation of Scale-Dependent Groundwater/Surface-water Exchange in Rivers by Gradient Self-Potential Logging: Numerical Modeling and Field Experiments

Exchanges of groundwater and surface-water are fundamental to a wide range of water-supply and water-quality management issues but challenging to map beyond the reach scale. Waterborne gradient self-potential (SP) measurements are directly sensitive to water flow through riverbed sediments and can be used to infer exchange locations, direction (gain versus loss), scale, and relative changes, but to date applications to river corridor hydrology are limited. Numerical modeling and field experiments were therefore performed herein, each emphasizing waterborne gradient SP logging for identifying and locating focused vertical groundwater discharge (surface-water gain) and recharge (surface-water loss) in a river. Two and three-dimensional numerical models were constructed to simulate the polarities, appearances, and peak amplitudes of streaming-potential and electric-field anomalies on a riverbed and in the surface-water that were attributable to steady-state vertical fluxes of groundwater through high-permeability conduits in the riverbed. Effects of varied hydraulic length-scale of exchange and surface-water depth were tested through numerical modeling. Modeling results aided in data acquisition and interpretation for three repeated field experiments performed along a 1.5–2.0 km reach of the Quashnet River in Cape Cod, Massachusetts, where focused, meter-scale groundwater discharges occur at discrete locations within otherwise ubiquitous and more diffuse groundwater upwelling conditions. Strong gradient SP anomalies were repeatedly measured in the Quashnet River at previously confirmed locations of focused groundwater discharge, showing the efficacy of waterborne gradient SP logging in identifying and characterizing groundwater/surface water exchange dynamics at multiple river network scales.

Engineering↗
DOE OSTI2018

Experimental and modeling study of calcium carbonate precipitation and its effects on the degradation of oil well cement during carbonated brine exposure

Decalcification of cement in solutions of carbonated brine is important to a host of engineering applications, especially in subsurface service environments where cementitious materials are frequently utilized as engineered barriers for wellbore seals, as well as shaft and drift seals and waste forms for nuclear waste disposal. Analysis of leaching simulations and experiments shows that, depending on solution compositions (dissolved CO{sub 2} concentration, pH, Ca ion concentration), calcite precipitation occurring during leaching of cement in contact with carbonated brine can have a significant impact on cement reactivity, in some instances leading to complete arrest of reactivity via calcium carbonate ''pore-clogging''. We present modeling and experimental results that examine the range of solution conditions that can lead to pore-clogging. Analysis of the results shows that distinct regimes of leaching behavior, based on pH and pCO2, can be used to form a framework to better understand the occurrence of pore-clogging.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2013

Is the acceleration of anomalous cosmic rays affected by the geometry of the termination shock?

Historically, anomalous cosmic rays (ACRs) were thought to be accelerated at the solar-wind termination shock (TS) by the diffusive shock acceleration process. When Voyager 1 crossed the TS in 2004, the measured ACR spectra did not match the theoretical prediction of a continuous power law, and the source of the high-energy ACRs was not observed. When the Voyager 2 crossed the TS in 2007, it produced similar results. Several possible explanations have since appeared in the literature, but we follow the suggestion that ACRs are still accelerated at the shock, only away from the Voyager crossing points. To investigate this hypothesis closer, we study ACR acceleration using a three-dimensional, non-spherical model of the heliosphere that is axisymmetric with respect to the interstellar flow direction. We then compare the results with those obtained for a spherical TS. A semi-analytic model of the plasma and magnetic field backgrounds is developed to permit an investigation over a wide range of parameters under controlled conditions. The model is applied to helium ACRs, whose phase-space trajectories are stochastically integrated backward in time until a pre-specified, low-energy boundary, taken to be 0.5 MeV n{sup –1} (the so-called injection energy), is reached. Our results show that ACR acceleration is quite efficient on the heliotail-facing part of the TS. For small values of the perpendicular diffusion coefficient, our model yields a positive intensity gradient between the TS and about midway through the heliosheath, in agreement with the Voyager observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

A comprehensive modeling procedure for estimating statistical properties of forced ignition

In this work, a comprehensive modeling procedure for estimating the probability of ignition with application to high altitude relights of aircraft combustors is developed. In these configurations, an ignitor is used to introduce high-enthalpy discharge into a fuel-laden but stratified flow. Because of the inherent variabilities inflow conditions, kernel discharge process, and the chaotic turbulent flow, ignition cannot be described deterministically, but only as a probabilistic measure. The proposed modeling framework consists of three components. The turbulent flow is represented using the large eddy simulation (LES) framework. The ignition process is modeled using a manifold approach, where the initial kernel evolution is represented using a homogeneous reactor while the latter part of the evolution is represented as a competition between diffusion and chemical reactions using a flamelet-type mapping. A combined lookup table that can track the evolution of the kernel through these two distinct reaction stages is developed. The table lookup variables are solved within the LES framework. The resulting simulation tool is then embedded within an uncertainty quantification approach, where variations in the turbulent flow, as well as operating and kernel properties, are simulated using a Monte-Carlo-based sampling approach. Techniques to reduce computational cost are used to obtain a robust, numerically accurate, and physically representative model for engine relight. The method is validated using experimental data for ignition of methane/air mixtures. Due to the comprehensive nature of the modeling procedure, it is found that the simulation tool reproduces experimentally observed ignition probabilities over a wide range of operating conditions.

42 ENGINEERING↗
DOE OSTI2021

Drying model of a high salt content cementitious waste form: Effect of capillary forces and salt solution

Highlights: • Drying model for a high salt content cementitious waste form is developed. • Water vapor diffusion and capillary liquid flow are distinguished. • Capillary and salt solution effects are considered in vapor-liquid equilibrium. A water transport model coupling capillary liquid flow with vapor diffusion is developed to describe the drying process for a cementitious waste form with high salinity porewater. Vapor-liquid equilibrium is formulated as the driving force for vapor diffusion and the model accounts for pore capillary and high salinity effects on water thermodynamic activity. Pore filling and porewater surface tension as a function of pore size distribution and water saturation have been quantified for the material. Geochemical speciation modeling is used to simulate porewater activity as a function of composition over the range of saturation. The theoretical relationship between relative humidity and water saturation generally agrees with experimental measurement, and the developed model is capable of predicting drying rates under various external relative humidity conditions. The model was developed to be incorporated into reactive transport models considering the effects of drying such as salt redistribution and efflorescence.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modelling of edge plasma dynamics with active wall boundary conditions

A self-consistent 2D model is presented for transport in boundary plasma and plasma-facing material walls. Plasma dynamics in the domain is represented by a 2D collisional plasma fluid model in the edge-plasma code UEDGE, and transport of hydrogen and heat in the wall is represented by a system of reaction–diffusion equations in the 1D wall code FACE. To account for variation of parameters along the wall, in the coupled model multiple instances of the FACE code run in parallel. Here, the coupled model provides a tool for investigating a range of dynamic plasma–material interactions phenomena in 2D. For demonstration of its capability, one application of particular interest is the role of active wall in tokamak strike point sweeping proposed for mitigation of divertor heat loads. In the present study, the coupled calculations are applied to investigation of the impact of heat and hydrogen transport in the material wall on the divertor plasma and target heat load during sweeping of the target strike point for parameters of a high-power tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗