Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

LIM1TR (Lithium-Ion Modeling with 1-D Thermal Runaway) is an open-source code that uses the finite volume method to simulate heat transfer and chemical kinetics on a quasi 1-D domain. The target application of this software is to simulate thermal runaway in systems of lithium-ion batteries. The source code for LIM1TR can be found at https://github.com/ajkur/lim1tr. This user guide details the steps required to create and run simulations with LIM1TR starting with setting up the Python environment, generating an input file, and running a simulation. Additional details are provided on the output of LIM1TR as well as extending the code with custom reaction models. This user guide concludes with simple example analyses of common battery thermal runaway scenarios. The corresponding input files and processing scripts can be found in the “Examples” folder in the on-line repository, with select input files included in the appendix of this document.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

A cohesive phase-field model of ductile fracture in a finite-deformation setting is presented. The model is based on a free-energy function in which both elastic and plastic work contributions are coupled to damage. Using a strictly variational framework, the field evolution equations, damage kinetics, and flow rule are jointly derived from a scalar least-action principle. Particular emphasis is placed on the use of a rational function for the stress degradation that maintains a fixed effective strength with decreasing regularization length. The model is employed to examine crack growth in pure mode-I problems through the generation of crack growth resistance (J-R) curves. In contrast to alternative models, the current formulation gives rise to J-R curves that are insensitive to the regularization length. Numerical evidence suggests convergence of local fields with respect to diminishing regularization length as well.

Abstract not provided.

Performance assessment (PA) of geologic radioactive waste repositories requires three-dimensional simulation of highly nonlinear, thermo-hydro-mechanical-chemical (THMC), multiphase flow and transport processes across many kilometers and over tens to hundreds of thousands of years. Integrating the effects of a near-field geomechanical process (i.e. buffer swelling) into coupled THC simulations through reduced-order modeling, rather than through fully coupled geomechanics, can reduce the dimensionality of the problem and improve computational efficiency. In this study, PFLOTRAN simulations model a single waste package in a shale host rock repository, where re-saturation of a bentonite buffer causes the buffer to swell and exert stress on a highly fractured disturbed rock zone (DRZ). Three types of stress-dependent permeability functions (exponential, modified cubic, and Two-part Hooke's law models) are implemented to describe mechanical characteristics of the system. Our modeling study suggests that compressing fractures reduces DRZ permeability, which could influence the rate of radionuclide transport and exchange with corrosive species in host rock groundwater that could accelerate waste package degradation. Less permeable shale host rock delays buffer swelling, consequently retarding DRZ permeability reduction as well as chemical transport within the barrier system.

Ni-Cr alloys exhibit oscillatory segregation behaviors near low index surfaces, in which the preferred segregation species changes from Ni in the first layer to Cr in the second layer. In many dilute-alloy systems, this oscillatory pattern is attributed to the elastic release of stresses in the local lattice around the segregating solute or impurity atom. These stresses are mostly thought to originate from mismatches in the atomic size of the solute and host atoms. In Ni-Cr alloys, however, an appreciable mismatch in atomic size is not present, leading to questions about the origins of the oscillatory behavior in this alloy. Using density functional theory, we have modeled the segregation of a single Cr atom in the (100) and (111) surfaces of FCC Ni, an alloy which exhibits this oscillatory behavior. Using Bader charge analysis, we show that the negative energy correlates directly with the amount of charge on the Cr atom. As Ni atoms strip valence charge from the Cr, the Cr contracts slightly in size. The greatest contraction and highest positive charge for the Cr occurs when it is in the second layer of the surface where the system exhibits the oscillating negative segregation energy. We then find that this behavior persists in other alloy systems (Ag-Nb, Cu-Cr, Pt-Nb, and Pt-V), which exhibit similar atomic radii and electronegativity differences between host and solute to Ni-Cr. These represent alloys in which the host metal exhibits an FCC ground-state structure while the solute metal exhibits a BCC ground-state structure.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

U.S. nuclear power facilities face increasing challenges in meeting dynamic security requirements caused by evolving and expanding threats while keeping cost reasonable to make nuclear energy competitive. The past approach has often included implementing security features after a facility has been designed and without attention to optimization, which can lead to cost overruns. Incorporating security in the design process can provide robust, cost effective, and sufficient physical protection systems. The purpose of this work is both to develop a framework for the integration of security into the design phase of High Temperature Gas Reactors (HTGRs) that utilize pebble-based fuels and increase the use of modeling and simulation tools to optimize the design of physical protection systems. Specifically, this effort focuses on integrating security into the design phase of a model HTGR that meets current Nuclear Regulatory Commission (NRC) physical protection requirements and providing advanced solutions to improve physical protection and decrease costs. A suite of tools, including SCRIBE3D©, PATHTRACE© and Blender© were used to model a hypothetical, generic domestic HTGR facility. Physical protection elements such as sensors, cameras, barriers, and guard forces were added to the model based on best practices for physical protection systems. Multiple outsider sabotage scenarios were examined with four-to eight adversaries to determine security metrics. The results of this work will influence physical protection system designs and facility designs for U.S. domestic HTGRs. This work will also demonstrate how a series of experimental and modeling capabilities across the Department of Energy (DOE) Complex can impact the design of and complete Safeguards and Security by Design (SSBD) for SMRs. The conclusions and recommendations in this document may be applicable to all SMR designs.