Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Thermal diffusion”

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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Summary Report on Ion Irradiation Study of Ceramic Coating on Suppressing FCCI

Advanced cladding is critical for advanced nuclear reactors with an enhanced performance in radiation tolerance and neutron transparency. Using advanced cladding will ensure the adequate thermal conductivity and mechanical stability of the cladding base material, corrosion resistance, high-temperature coolant compatibility of the cladding surface, and chemical stability in the cladding inner wall against fuel cladding chemical interaction (FCCI). An innovative cladding with a three-layer structure (i.e., a modified surface, a clad base material, and a modified inner wall) promises to meet all these requirements. Initial research and development (R&D) regarding this innovative cladding seek to demonstrate the effectiveness of a thin ceramic coating in suppressing FCCI under ion irradiation to high dose. This United States (U.S.) Department of Energy (DOE)–Office of Nuclear Energy (NE) project report summarizes the results of recent ion irradiation studies of diffusion-couple samples with details of an ion irradiation experiment, the characterization of the interface microstructure of cerium (Ce)/titanium nitride (TiN)/oxide-dispersion strengthened (ODS), Ce/TiN/iron (Fe), Ce/ODS, and Ce/Fe samples irradiated with 80 million electron-volts (MeV) of xenon (Xe) ions to 100 displacements per atom (dpa) at 500?C using the Argonne Tandem Linac Accelerator System (ATLAS) facility at Argonne National Laboratory (ANL). The results of the ODS-substrate sample are encouraging and demonstrate the effectiveness of a thin ceramic coating on FCCI mitigation under irradiation. The results from an Fe-substrate sample revealed a complex microstructure, and the root cause is discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A Code-Agnostic Driver Application for Coupled Neutronics and Thermal-Hydraulic Simulations

While the literature has numerous examples of Monte Carlo and computational fluid dynamics (CFD) coupling, most are hard-wired codes intended primarily for research rather than as standalone, general-purpose applications. In this work, we describe an open source application, ENRICO, that enables coupled neutronic and thermal-hydraulic simulations between multiple codes that can be chosen at runtime (as opposed to a coupling between two specific codes). The application has been designed such that the control flow logic, domain mapping, nonlinear fixed-point iteration, solution transfers, and convergence checks are all agnostic to the underlying physics solvers used. Special emphasis has also been placed on enabling efficient execution on distributed-memory computing environments. The transfer of solution fields between solvers is performed in memory rather than through filesystem I/O. Additionally, solvers can be configured to run on overlapping or disjoint sets of processes. To date, coupling with the OpenMC and Shift Monte Carlo codes, the Nek5000 CFD code, and a simplified heat diffusion and subchannel solver has been implemented in ENRICO. We present results for coupled simulations of a single light-water reactor fuel assembly based on the NuScale reactor using various combinations of the physics solvers. For this problem, the coupled simulations are shown to converge in about four Picard iterations. A comparison of the heat source and temperature distributions computed by ENRICO using OpenMC coupled with Nek5000 and Shift coupled with Nek5000 illustrates remarkable agreement between the codes.

42 ENGINEERING↗
DOE OSTI2024

Ab-initio molecular dynamics study of eutectic chloride salt: MgCl2–NaCl–KCl

Ionic liquid materials are viable candidates as a heat transfer fluid (HTF) in a wide range of applications, notably within concentrated solar power (CSP) technology and molten salt reactors (MSRs). For next-generation CSP and MSR technologies that strive for higher power generation efficiency, a HTF with wide liquid phase range and energy storage capabilities is crucial. Studies have shown that eutectic chloride salts exhibit thermal stability at high temperatures, high heat storage capacity, and are less expensive than nitrate and carbonate salts. However, the experimental data needed to fully evaluate the potential of eutectic chloride salts as a HTF contender are scarce and entail large uncertainties. Considering the high cost and potential hazards associated with the experimental methods used to determine the properties of ionic liquids, molecular modeling can be used as a viable alternative resource. In this study, the eutectic ternary chloride salt MgCl 2 –NaCl–KCl is modeled using ab-initio molecular dynamics simulations (AIMDs) in the liquid phase. Using the simulated data, the thermophysical and transport properties of eutectic chloride salt can be calculated: density, viscosity, heat capacity, diffusion coefficient, and ionic conductivity. For an initial model validation, experimental pair-distribution function data were obtained from X-ray total scattering techniques and compared to the theoretical pair-distribution function. Additionally, theoretical viscosity values are compared to experimental viscosity values for a similar system. The results provide a starting foundation for a MgCl 2 –NaCl–KCl model that can be extended to predict other fundamental properties.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Chemical structure and curing dynamics of bisphenol S, PEEK™-like, and resveratrol phthalonitrile thermoset resins

This effort provides a multifaceted analysis of the structural changes and material dynamics of thermally driven softening and curing of three distinct phthalonitrile (PN) resins that cross-link into thermally stable and oxidation-resistant thermosets. Although this material system has yielded a large subset of fire-retardant composites that require facile processing and low-temperature curing, to date, insufficient information had been available on the fundamental processes that drive their softening and curing stages. Our approach conducted a complementary analysis the chemistry, monomer mobility, and rheology of three PN polymers in order to correlate the curing processes with corresponding structural and behavioral transformations of thermosets. Here, we focused on PNs with a bisphenol S backbone, a bisphenol A (PEEK™-like) backbone, and a resveratrol backbone. We relied on quasi-elastic neutron scattering (QENS) in order to analyze the in situ dynamics and self-diffusion properties of PN monomers, and to track changes in their mobilities during cross-linking and staging. Our analysis facilitates proper control over the staging and final curing of these resins and enables more efficient processing of these thermosets.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Evaluation of high-temperature properties of binary, eutectic Ni-X alloys (X = B, Ca, Ce, La, Y, Zr)

We investigated the high-temperature behavior of binary eutectic nickel-base alloys with low-solubility elements, towards developing alloys with new combinations of structural and functional (e.g. thermal or corrosion) properties. Based on thermodynamic screening, we cast binary near-eutectic Ni-B, Ni-Ca, Ni-Ce, Ni-La, Ni-Y, and Ni-Zr alloys and characterized their microstructure and high-temperature properties. All alloys (except Ni-Ca) contain a fine eutectic microstructure, containing submicron alternating lamellae of γ-Ni and continuous intermetallic phases. The eutectic microstructure did not significantly coarsen at 700°C (except for Ni-B), but measurably coarsened in all alloys at 900°C, most quickly in Ni-B, followed by Ni-Ce, Ni-La, Ni-Y, and Ni-Zr. Coarsening in these alloys occurs by fault migration; diffusion within the intermetallic, and its interfacial energy, likely influence coarsening kinetics. At room-temperature, all alloys possess negligible ductility due to the continuous, brittle intermetallic. At 700 and 900°C, the Ni-Zr alloy showed the best combination of strength and ductility. Oxidation experiments at 900°C on Ni-Zr and Ni-Ce revealed heavy internal oxidation of the continuous intermetallic phase. Overall, the Ni-Zr alloy showed good high-temperature coarsening resistance and strength, though ductility and oxidation resistance require improvement for use in structural applications. Possible routes for improving alloy properties via additive manufacturing are discussed.

Ekaputra, Clement [ORNL] (ORCID:0000000307166479)↗
DOE OSTI2022

Dynamical chemistry: non-equilibrium effective actions for reactive fluids

Abstract We present two approaches for describing chemical reactions taking place in fluid phase. The first method mirrors the usual derivation of the hydrodynamic equations of motion by relating conserved—or to account for chemical reactions, non-conserved—currents to local-equilibrium parameters. The second method involves a higher-brow approach in which we attack the same problem from the perspective of non-equilibrium effective field theory (EFT). Non-equilibrium effective actions are defined using the in–in formalism on the Schwinger–Keldysh contour and are therefore capable of describing thermal fluctuations and dissipation as well as quantum effects. The non-equilibrium EFT approach is especially powerful as all terms in the action are fully specified by the symmetries of the system; in particular the second law of thermodynamics does not need to be included by hand, but is instead derived from the action itself. We find that the equations of motion generated by both methods agree, but the EFT approach yields certain advantages. To demonstrate some of these advantages we construct a quadratic action that is valid to very small distance scales—much smaller than the scales at which ordinary hydrodynamic theories break down. Such an action captures the full thermodynamic and quantum behavior of reactions and diffusion at quadratic order. Finally, taking the low-frequency and low-wavenumber limit, we reproduce the linearized version of the well-known reaction–diffusion equations as a final coherence check.

Mechanics↗
DOE OSTI2025

Code Coverage Status of the ARC Code PERSENT

The Argonne Reactor Code (ARC) software system supports users in their fast reactor design goals by providing neutronic, thermal-hydraulic, and structural analysis capabilities. PERSENT fulfills the role of generating reactivity coefficients for a given time point of a REBUS calculation usable in a point kinetics based safety analysis capability. PERSENT also provides a sensitivity coefficient capability on eigenvalue, reactivity worth, and several other key coefficients that are used in the follow-on safety analysis. Given a co-variance matrix, PERSENT can carry out the uncertainty quantification to indicate the amount of error in the reactivity coefficients derived from the errors in the cross section measurements. With continued improvement of computational resources, many of the geometry modeling capabilities in DIF3D that were primarily used in low order schemes are not really needed anymore. Today, the diffusion and transport capabilities of DIF3D-VARIANT are primarily used in the reactor design process with some scattered usage of DIF3D-FD and DIF3D-Nodal. PERSENT is part of the ARC code system and is built around DIF3D-VARIANT and the flux solution it provides. The purpose of the present work is to identify a set of test problems for PERSENT and assess the code coverage of PERSENT for those test problems. PERSENT treats the DIF3D executable as an external executable and thus the code coverage considerations only need to focus on the PERSENT source code and only a fraction of the connected modules in the existing ARC software library. The goal is to document what parts of the existing PERSENT code are touched by the set of test problems and which are not. Because the verification work done on PERSENT was focused on the most common uses of PERSENT for fast reactor analysis, the code coverage assessment of those capabilities is the highest priority. This will ensure that nothing is being missed by the existing verification test problems that users of PERSENT rely upon. The code coverage analysis of PERSENT was performed with the Code Coverage Tool of the Intel Fortran compiler which requires modifications to the compilation of PERSENT. The detailed coverage tables are given for each submodule of PERSENT. Most of the uncovered parts/files could be easily ignored because they are either for error message and debugging output or not needed by PERSENT today. Only a few uncovered parts of PERSENT deserve extending the verification test suite.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Effect of ammonia addition on suppressing soot formation in methane co-flow diffusion flames

Due to issues surrounding carbon dioxide emissions from carbon-containing fuels, there is growing interest in ammonia (NH 3 ) as an alternative combustion fuel. One attractive method of burning NH 3 is to co-fire it with hydrocarbons, such as natural gas, and in this case soot formation is possible. To begin understanding the influence of NH 3 on soot formation when co-fired with hydrocarbons, soot volume fractions and mole fractions of gas-phase species were computationally and experimentally interrogated for CH 4 flames with up to 40% NH 3 by volumetric fuel fraction. Mole fractions of gas-phase species, including C 2 H 2 and C 6 H 6 , were measured with on-line electron impact mass spectrometry, and soot volume fractions were obtained via color-ratio pyrometry. The simulations employed a detailed chemical mechanism developed for capturing nitrogen interactions with hydrocarbons during combustion. The results are compared to findings in N 2 -CH 4 flames, in order to separate thermal and dilution effects from the chemical influence of NH 3 on soot formation. Experimentally, C 2 H 2 concentrations were found to decrease slightly for the NH 3 -CH 4 flames relative to N 2 -CH 4 flames, and a stronger suppression of C 6 H 6 was found for NH 3 relative to N 2 additions. The measured results show a strong suppression of soot with the addition of NH 3 , with soot concentrations reduced by over a factor of 10 with addition of up to 20% or more NH 3 by mole fraction. The model satisfactorily captured relative differences in maximum centerline C 2 H 2 , C 6 H 6 , and soot concentrations with addition of N 2 , but was unable to match measured differences in NH 3 -CH 4 flames. Furthermore, these results highlight the need for an improved understanding of fuel-nitrogen interactions with higher hydrocarbons to enable accurate models for predicting particulate emissions from NH 3 /hydrocarbon combustion.

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

Nanoscale imaging of phonon dynamics by electron microscopy

Spatially resolved vibrational mapping of nanostructures is indispensable to the development and understanding of thermal nanodevices, modulation of thermal transport and novel nanostructured thermoelectric materials. Through the engineering of complex structures, such as alloys, nanostructures and superlattice interfaces, one can significantly alter the propagation of phonons and suppress material thermal conductivity while maintaining electrical conductivity. There have been no correlative experiments that spatially track the modulation of phonon properties in and around nanostructures due to spatial resolution limitations of conventional optical phonon detection techniques. Here we demonstrate two-dimensional spatial mapping of phonons in a single silicon–germanium (SiGe) quantum dot (QD) using monochromated electron energy loss spectroscopy in the transmission electron microscope. Tracking the variation of the Si optical mode in and around the QD, we observe the nanoscale modification of the composition-induced red shift. We observe non-equilibrium phonons that only exist near the interface and, furthermore, develop a novel technique to differentially map phonon momenta, providing direct evidence that the interplay between diffuse and specular reflection largely depends on the detailed atomistic structure: a major advancement in the field. Our work unveils the non-equilibrium phonon dynamics at nanoscale interfaces and can be used to study actual nanodevices and aid in the understanding of heat dissipation near nanoscale hotspots, which is crucial for future high-performance nanoelectronics.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Dynamics of iodine geminate recombination in supercritical xenon solvent: Caging effect

Understanding the dynamics of chemical reactions in solutions is vital, as their rates and kinetics are significantly affected by the solvent environment. Supercritical solvents offer extensive applications in chemical reactions by enabling the manipulation of the solution environment. Here, in this study, we investigate the geminate recombination of iodine in a supercritical xenon solvent by using ReaxFF-based molecular dynamics simulations. Our findings reveal that the highest iodine recombination rate occurs near supercritical conditions, while lower-pressure conditions lead to reduced collision rates and unstable recombination, and higher-pressure conditions hinder iodine diffusion, resulting in a lower recombination rate. Our analysis shows that the xenon local density at the time of recombination is at least 2.5 times higher than the global density, confirming the presence of xenon clusters surrounding the Iodine atoms. This observation is further supported by coordination number analysis, which confirms an elevated xenon local density during recombination. In addition, the correlation between the total energy of xenon atoms within a cluster and recombined iodine atoms underscores the kinetic energy transfer process, validating the occurrence of geminate recombination. The excess kinetic energy from the recombining iodine atoms is transferred to the surrounding xenon atoms. Our examination of geminate recombination demonstrates that iodine atoms confined within xenon clusters—whether through manual insertion of atoms or the fast dissociation of an iodine molecule within xenon clusters—are more likely to recombine as primary geminate recombination. However, extending the iodine molecule dissociation time allows iodine atoms to diffuse out of the cluster, and the recombination to shift toward secondary geminate recombination.

Cage effect↗
DOE OSTI2023

Nanoscale mapping of point defect concentrations with $\mathrm{4D-STEM}$

Vacancies are missing atoms in a crystalline material, and occur both at equilibrium (varying with temperature) and out of equilibrium such as when crystalline materials are damaged with radiation or corrosion. While we know of their importance, particularly regarding diffusive mechanisms, it is not straightforward to experimentally measure their concentration or to visualize them directly. Traditionally, measurements such as positron annihilation spectroscopy and to some degree X-ray diffraction can measure average concentrations, but generally lack the ability to visualize or quantify a heterogenous concentration of vacancies that can occur at the level of individual defects and microstructural features. Here, we present a method to map vacancy concentrations and their distribution using local lattice parameter measurements with a high-resolution electron microscope. Our method utilizes a Au thin film as a model to demonstrate the method via four-dimensional scanning transmission electron microscopy (4D-STEM) by correlating the differences between changes in lattice parameter and the volumetric thermal expansion during in situ heating experiments. The vacancy mapping methodology is also applied to non-equilibrium defects accumulated in pure Al via knock-on electron beam irradiation. Our method demonstrates the ability to map point defect concentrations in heterogeneous systems in situ with nanometer spatial resolution. The result is a technique that can provide direct measurements of vacancy concentrations at the level of individual defects in studies of materials in and out of equilibrium.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Brazing characteristics, microstructure, and wettability of laser powder bed fusion additive manufactured GRCop-84 compared to CuCrZr and OFC, and brazing to titanium-zirconium-molybdenum alloy limiters

Laser Powder Bed Fusion (L-PBF) of Glenn Research Copper 84 (GRCop-84), a Cr 2 Nb (8 at. % Cr, 4 at. % Nb) precipitation hardened alloy, produces a fully dense, high conductivity alloy with a yield strength of 500 MPa and ultimate tensile strength (UTS) of 740 MPa with 20% elongation; superior to other competing copper alloys. Braze wetting characteristics of GRCop-84 with Ag-Cu-X, and Au-Cu brazes were similar to CuCrZr, but less than oxygen free copper. No difference in wetting was observed between infill and surface contour areas in L-PBF GRCop-84. Wet sanding to 240 grit (R a =0.24 µm) was considered the optimal surface condition. Silver diffusing through GRCop-84 depleted Cr 2 Nb precipitates from the copper grain and deposited agglomerations of coarsened precipitates within silver-rich regions of intergranular diffusion once a density threshold was reached. Microstructure modification was minimized with 50Au-50Cu braze implying that silver caused precipitate coarsening and agglomeration, and not high temperature exposure. Coarsened precipitates were observed on the surface within braze pools implying a contribution to braze wetting. Palcusil-25, Ticusil, CuSil-ABA, and 50Au-50Cu brazes were suitable for brazing to unplated Titanium-Zirconium-Molybdenum (TZM), while sulfamate nickel plating to allows wetting with CuSil or other non-active brazes. Additionally, vacuum brazing techniques were developed to join a 1 mm thick layer of TZM to the front of additive manufactured GRCop-84 waveguides considering the brazing characteristics of both GRCop-84, TZM, and internal stress from the difference in coefficient in thermal expansion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Development of MOSCATO: A CFD-Level Electrochemistry and Corrosion Simulator for Molten Salt Systems

For both coolant and fueled variants of molten salt reactors (MSRs), the corrosion of structural materials is a significant challenge. The corrosion stems from chemical and electrochemical reactions initiated by fissile material, fission products, and impurities in the salt. Lower-fidelity models rely on empirical correlations for mass transfer, simplified lumped temperature profiles, and similar assumptions. They do not capture detailed spatial variations in complex geometries, creating the need for high-fidelity modeling to bridge this gap.As we approach the demonstration and possible deployment of MSRs in this decade, the development of a high-fidelity, high-performance simulator becomes imperative. To simulate the complex electrochemical environment and corrosion within molten salt systems, we have developed the Molten Salt Chemistry And TranspOrt (MOSCATO) code. This endeavor is comprised of three essential components. First, mass transfer equations are coupled with the Navier-Stokes equations in order to account for the transport of species in the salt. Second, the diffusion of alloy constituents, such as Cr, Fe, Ni, etc. is simulated within the structural metals. Third, the alloy and salt domains are coupled to account for the heterogeneous chemical and electrochemical reactions that occur at the salt-alloy interface.MOSCATO manages all three components within the framework of the highly scalable, open-source spectral element method computational fluid dynamics code Nek5000/NekRS. This integration enables MOSCATO to harness the immense computational power of modern high-performance computing resources, ensuring both high fidelity and computational speed.In addition to code development, we have initiated a comprehensive verification and validation campaign, utilizing data from diverse sources. First, MOSCATO's electrochemical solver was verified with reference numerical data. Then validation occurred against experiments: one of a thermal galvanic cell and the other for corrosion in flowing molten salt of FLiNaK (LiF-NaF-KF). This campaign verified and validated MOSCATO as a reliable tool for simulating electrochemical environments and corrosion in molten salt systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

DECOVALEX-2023: Task F2 Salt Final Report

The subject of Task F of DECOVALEX-2023 concerns performance assessment modelling of radioactive waste disposal in deep mined repositories. The primary objectives of Task F are to build confidence in the models, methods, and software used for performance assessment (PA) of deep geologic nuclear waste repositories, and/or to bring to the fore additional research and development needed to improve PA methodologies. In Task F2- (salt), these objectives have been accomplished through staged development and comparison of the models and methods used by participating teams in their PA frameworks. Coupled-process submodels and deterministic simulations of the entire PA model for a reference scenario for waste disposal in domal salt have been conducted. The task specification has been updated continuously since the initiation of the project to reflect the staged development of the conceptual repository model and performance metrics. Thermal, hydrological, mechanical, and chemical properties of individual components of the engineered and natural system were chosen for relevance by participating teams. The salt reference case system was characterized using data and measurements collected at relevant underground research laboratories (URLs), field sites, and simulation results from teams with specialized modelling capability. Participating teams made a wide range of model assumptions from compartmentalized networks to full 3D models of the salt formation. No single contributed model includes full-fidelity representation of all the features, events, and processes (FEPs) detailed in the task specification, but almost all features and processes are represented in at least one model. Despite differences in the modelling strategies developed by participating teams, all models indicate that salt compaction and radionuclide diffusion are key processes in the repository, and for the FEPs and model scenario considered, little of the disposed radionuclides will migrate beyond the repository seal over the 100,000 year simulations. In general, the model output quantities have the largest differences over the short term and near the waste. The models tend to be more similar further from waste and at later time. Disparities between the models are believed to be due to differing simplifications from the task specification, some of which are chosen simplifications to reduce complexity, and some are restrictions imposed by the modelling tools. A second round of this task has been accepted for DECOVALEX-2027 in conjunction with Task F1 on crystalline PA modelling. The future round includes waste package heating, improved modelling of salt creep closure, additional comparisons of coupled-process sub-models, and the impact of repository engineering design on radionuclide migration in the repository. Participants will also propose and finalize a set of uncertain inputs for the reference case simulations, propagate these uncertainties in a set of realizations, and conduct sensitivity analyses on the simulation results.

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

Improved lithium storage capacity and high rate capability of nitrogen-doped graphite-like electrode materials prepared from thermal pyrolysis of graphene quantum dots

Adopting a solid-phase microwave-assisted technique followed by thermal pyrolysis of N-functionalized graphene quantum dots, novel nitrogen-doped graphite-like (NGL) electrode materials were synthesized in this work and served as the anode for Li-ion batteries. The NGL anode demonstrated reversible capacity of 530 mAh g -1 at 0.1C, superior rate capability at high C rate operation (420 mAh g -1 at 5C), remarkable initial coulombic efficiency (>95.7%), and excellent cyclic stability along with high efficiency (>99.1%) during entire cycling. The NGL anode nanostructure enables improved lithium ion mobility and reversible Li + storage during cycling. The analysis of the Ragone plots revealed that the specific energy of NGL anode reaches to ca. 840 Wh kg -1 at the power density of 4200 W kg -1 . The diffusion coefficient of Li ions was measured as 1.69 × 10 -9 cm 2 s -1 for the NGL anode material, substantially improving over commonly used graphite electrodes (15–26 times higher Li + diffusivity). The high-rate cyclability as well as the cyclic stability of the NGL anodes were also confirmed via long-term cycling of full pouch cells assembled with ternary cathode and NGL anode. The robust design of the NGL anode materials introduced in this work, paves the way for designing next-generation lithium-ion batteries operating at ultra-high C rates.

25 ENERGY STORAGE↗
DOE OSTI2021

ORNL Input to GDSA Repository Systems Analysis FY21

This document satisfies the M3 milestone M3SF-21OR010304082 titled “ORNL Input to GDSA Repository Systems.” This document describes the current status of Oak Ridge National Laboratory’s (ORNL’s) efforts related to analyses of dual-purpose canister (DPC) disposal in unsaturated alluvium with a focus on thermal hydrological constraints on criticality timing and power output. This analysis updates previous analyses of thermal hydrological constraints on timing and power output of a potential criticality event in Dual-Purpose Canisters (DPCs) in a hypothetical repository in unsaturated alluvium using a more realistic representation of heat transport inside dry DPCs. PFLOTRAN was used for the simulations of multiphase thermal hydrology near a single DPC. The scenario considers a DPC failure at 9000 years, allowing water to enter and eventually fill the DPC. Once the DPC is filled to a level that could support a criticality event, different values for criticality power output were added to the decay heat. The main objective is to bound the power output that could be produced by a criticality event without driving water out of the package. For the conditions analyzed here, following a package breach, the alluvial formation could supply enough water to allow enough accumulation in the DPC to support a criticality event. However, the power output that would be generated is limited to modest values by loss of water moderator due to evaporation and vapor diffusion. In the reference case scenario, the DPC would not start to fill with water until about 16,000 years post closure and would not fill to a level that allows a criticality event until approximately 25,000 years post closure. The long-term average power output that could be sustained without driving off the water and terminating the criticality event is limited to approximately 100 W. Sensitivity to assumed conditions and parameters in the reference case, especially the deep percolation rate, need to be addressed but could not be undertaken because of numerical failures of the PFLOTRAN code running in general mode in the dry conditions of an unsaturated alluvium repository

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

Phase Field Modeling of Chemical Reaction Related Damage Evolution in Environmental Barrier Coatings

The advent of next-generation engines necessitates materials capable of withstanding temperatures beyond the reach of current superalloys. SiC-based ceramic matrix composites, augmented with environmental barrier coatings (EBCs), present a promising materials solution. Given the active search for effective and durable EBCs, there is a pressing need for modeling tools to understand and predict damage evolution in these materials to help accelerate their development. This study introduces a phase-field model (PFM) designed to simulate the thermally grown oxides (TGO) and phase transformation in the degradation and failure of EBCs. The model accounts for the severe volume expansion due to oxidation, alongside phase transformations and porosity evolution during thermal cycling, offering a comprehensive view of the damage processes. Simulation results are validated against experimental findings reported in the literature, establishing the model's potential as a significant tool for understanding and improving the resilience of EBCs in cyclic oxidative environments.

fast-diffusion path↗