Generating Multiphase Fluid Configurations in Fractures using Diffusion Models
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
This Career Award accomplished most of its goals, in demonstrating for the first time a low- threshold nuclear recoil detector based on the quantum evaporation of helium-4. Through this work, we trained several PhD students and evolved this niche R&D topic into a central pillar of the DOE-funded TESSERACT project, which now plans to take underground helium data in 2028.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Lightning talk presentation for the ignite off competition consisting of 20 picture-focused slides, as submitted to the national competition hosted by ORNL.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
To fulfill the design and licensing requirements for advanced reactors, the U.S. Department of Energy (DOE) is specifically interested in targeted research to quantify the mechanisms of fission product transport and retention in graphitic grades that will likely be used in next-generation very/high temperature reactors, fluoride salt-cooled high temperature reactors and molten salt reactors. We have examined an integrated experimental-computational approach to determine the diffusivities of selected fission products – ruthenium (Ru) and silver (Ag) in five nuclear graphite grades – POCO AXF-5Q, POCO ZXF-5Q, PCEA, IG110, and NBG-18. Experimental investigation for diffusion of other fission products (iodine, cesium, and krypton) in graphite could not be carried out due to extreme evaporation due to high vapor pressure, despite several attempts incorporating various experimental schemes. Experimental investigation with irradiated graphite could not be carried out because necessary approval to transport irradiated graphite could not be obtained due to procedural reasons at ORNL and national policy change in UK. An effort was made to “produce” by ion sputtering, but it was discontinued because the outcome did not adequately represent irradiated graphite. Thin film and “cup-and-cap” methods were employed to deposit fission products, diffusion annealed were carried out in an argon atmosphere, and concentration profiles were determined by dynamic secondary ion mass spectroscopy depth profiling technique. Regardless of graphite type, anomalous diffusion behavior of Ru (i.e., lower than expected based on Arrhenius temperature dependence) at an intermediate temperature, around 600 to 700 °C, was observed. Diffusion behavior of silver on the other hand agreed well to the Arrhenius temperature dependence although the excessive evaporation of silver hampered the accurate determination of diffusivities at temperature higher than 900 °C. In general, silver was observed to diffuse faster than Ru in the temperature range compared from 500 to 900 °C, had consistently had higher pre-exponential factor corresponding to higher jump frequency. To understand the abnormal diffusional behavior of Ru in graphite in intermediate temperature range, a detailed microstructural analysis was carried out. Clustering of Ru into particles on the surface of graphite was observed around 700 °C corresponding to de-wetting and spheroidization. This change would influence the boundary condition of diffusion, although features of clustering was much smaller (<1 mm) than the SIMS raster size (120 x 120 mm). Raman spectroscopy also demonstrated that Raman peak intensity for Ru decreased significantly for samples annealed higher than 700 °C and presence of Ru in graphite altered the crystallinity of graphite. To further elucidate the abnormal diffusion behavior observed in intermediate temperature investigated, the vibrational properties of graphene, bilayer graphene and Kr atoms in bilayer graphene were investigated with molecular dynamics. Krypton, which is a fission product, was introduced between a bilayer graphene system to evaluate (i) whether the phonon dispersion curves change with high temperature in bilayer graphene (ii) whether the addition of Kr atoms between the bilayer introduces any significant change in the vibrational properties, and (iii) whether any change in the vibrational properties can be correlated to the diffusive behavior of Kr atoms. The results show that the Kr atoms indeed exhibit an anomalous diffusive behavior at 550°C. Phonon analysis shows that a few phonon frequencies are enhanced as much as 20 to 30 % with the addition of just 0.2 at.% Kr. So, this may explain the unexpected self-diffusivity of fission products.
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.
The Nabarro-Herring (NH) diffusional creep theory postulates the vacancy-mediated transport of atoms under a stress gradient as the creep mechanism under low-stress and high-temperature conditions. In multicomponent alloys, we premise that this stress-assisted flow of vacancies to and from grain boundaries will produce elemental segregation. An observation of such segregation, validated with theoretical predictions, can provide the necessary experimental evidence for the occurrence of NH creep. Theoretical calculations of the segregation tendencies via analyzing the dominant solute diffusion mechanisms and the difference in diffusivities of the elements are therefore essential. To this end, this study applies density functional theory calculations of migration barriers and solute-vacancy binding energies as input to the self-consistent mean-field theory to assess the vacancy-mediated diffusion mechanisms, transport coefficients, and segregation tendencies of Co, Cr, Mo, Re, Ta, and W solutes in face-centered-cubic Ni. We find Co, Re, and W to be slow diffusers at high temperatures and Cr, Mo, and Ta to be fast diffusers. Further analysis shows that the slow diffusers tend to always enrich at vacancy sinks over a wide range of temperatures. In contrast, the fast diffusers show a transition from depletion to enrichment as the temperature lowers. Furthermore, our analysis of the segregation tendencies under tensile hydrostatic strains shows that slow diffusers are largely unaffected by the strain and favor enrichment. On the other hand, the fast diffusers exhibit high sensitivity to strain and their segregation tendency can transition from depletion to enrichment at a given temperature. Finally, the transport coefficients calculated in this work are expected to serve as input to mesoscale microstructure models to provide a more rigorous assessment of solute segregation under NH creep conditions.
In this work, a composite continuum-pore network formulation is presented to model single and two-phase transport thin porous layers, such as gas diffusion layers in polymer electrolyte fuel cells (PEFCs) and active electrodes in redox flow batteries (RFBs). The formulation can be integrated into CFD codes, thus combining the ease of implementation of continuum-based modeling and the computational power of pore-network modeling. The composite model includes a control volume (CV) mesh at the layer scale, which embeds a cubic pore network [1,2]. The pore-network model is used to determine analytically local anisotropic effective transport properties (local effective diffusivity and permeability), which are mapped onto the CV mesh to simulate transport in the porous transport layer [3,4]. Good agreement is found between the predicted global effective transport properties (global effective diffusivity and permeability) under dry and wet conditions and previous experimental data reported in the literature for Toray TGP-H series carbon paper. Water saturation distributions are also compared with results obtained using X-ray computed tomography [4]. [1] P.A. García-Salaberri, I.V. Zenyuk, J.T. Gostick, A.Z. Weber, Modeling Gas Diffusion Layer in Polymer Electrolyte Fuel Cells Using a Continuum-Based Pore-Network Formulation, ECS Trans. 97 (2020) 615. [2] P.A. García-Salaberri, Modeling diffusion and convection in thin porous transport layers using a composite continuum-network model: Application to gas diffusion layers in polymer electrolyte fuel cells, Int J. Heat Mass Transf. (2020), submitted. [3] P.A. García-Salaberri, J.T. Gostick, G. Hwang, A.Z. Weber, M. Vera, Effective diffusivity in partially-saturated carbon-fiber gas diffusion layers: Effect of local saturation and application to macroscopic continuum models, J. Power Sources 296 (2015) 440–453. [4] P.A. García-Salaberri, G. Hwang, M. Vera, A.Z. Weber, J.T. Gostick, Effective diffusivity in partially-saturated carbon-fiber gas diffusion layers: Effect of through-plane saturation distribution, Int. J. Heat Mass Transf. 86 (2015) 319–333.
Metal hydrides are a promising moderator material for high temperature fission reactors. Yttrium hydride can be loaded to a high hydrogen density with relatively high hydrogen stability at temperatures up to 800°C. This makes yttrium hydride a potential moderator material for microreactors as foreseen in the fission surface power program. However, during the operation of such advanced reactors temperature gradients are expected which can change the local hydrogen density in the moderator. Hydrogen diffusion in metals is driven by a concentration gradient (Fick’s law) and thermal diffusion (Soret diffusion). Thermal diffusion is the transport of hydrogen, or other species, due to a temperature gradient. For example, hydrogen might migrate from the hot side of a sample to the cold side of a sample. Measuring Fickian diffusion is achieved through various permeation or absorption experiments, however measuring thermal diffusion is challenging and has rarely been performed. The Hydrogen Experimental Apparatus for Thermal Diffusion (HEATD) experiment is designed to induce thermal diffusion in samples and quench those samples so that the hydrogen distribution can be analyzed using hot vacuum extraction (HVE). One side of the sample was heated to a high temperature e.g., 800°C, while the other side of the sample is at a lower temperature. The sample was held under the applied temperature gradient for a given time until the anticipated hydrogen diffusion has occurred. The actual time depends depend on the sample composition and hydrogen concentration. The results from HVE showed that thermal diffusion took place in the specimen and the Soret coefficient was calculated.
Abstract A model based on the code CrunchClay is presented for a fracture-clay matrix system that takes electrostatic effects on transport into account. The electrostatic effects on transport include those associated with the development of a diffusion potential as captured by the Nernst-Planck equation, and the formation of a diffuse layer bordering negatively charged clay particles within which partial anion exclusion occurs. The model is based on a dual continuum formulation that accounts for diffuse layer and bulk water pore space, providing a more flexible framework than is found in the classical mean electrostatic potential models. The diffuse layer model is obtained by volume averaging ion concentrations in the Poisson-Boltzmann equation, but also includes the treatment of longitudinal transport within this continuum. The calculation of transport within the bulk and diffuse layer porosity is based on a new formulation for the Nernst-Planck equation that considers averaging of diffusion coefficients and accumulation factors at grid cell interfaces. Equations for function residuals and the associated Jacobian matrix are presented such that the system of nonlinear differential-algebraic equations can be solved with Newton’s method. As an example, we consider a 2D system with a single discrete fracture within which flow and advective transport occurs that is coupled to diffusion in the clay-rich matrix. The simulation results demonstrate the lack of retardation for anions (e.g., 36 Cl − ) of the contaminant plume within the fracture flow system because they are largely excluded from the charged clay rock, while the migration of cations (e.g., 90 Sr ++ ) is more strongly attenuated. The diffusive loss of divalent cations in particular from the fracture is accentuated by their accumulation in the diffuse layer within the clay-rich matrix.