Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Code of Accounts”

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 271 records · Page 15

Center for Tokamak Transient Simulations (RPI Unstructured Mesh Developments for FES SciDAC4 Partnerships) (Final Report)

The goal of this project was the development of unstructured mesh technologies for fusion simulation codes” for FES SciDAC partnerships and to integrate those technologies into the simulation workflows of those partnerships. Specific developments were executed in support of the following FES SciDAC4 partnerships: Partnership Center for High‐fidelity Boundary Plasma Simulation (HBPS), Center for Integrated Simulation of Fusion Relevant RF (RF‐SciDAC), Center for Plasma Surface Interactions: Predicting the Performance and Impact of Dynamic PFC Surfaces (PSI2), and Center for Tokamak Transient Simulations (CTTS). The key unstructured mesh development areas addressed in this project include (i) methods to most effectively perform PIC calculations on unstructured meshes; (ii) creating meshes for fusion systems accounting for any desired level of geometric complexity and providing physics aware mesh configurations, (iii) adapting unstructured meshes to control the discretization errors, (iv) executing unstructured mesh calculations on GPU accelerated systems, (v) supporting physics/application‐specific PIC operations including surface/wall interactions of particles, (vi) providing infrastructure tools to support the interactions of solvers with unstructured meshes, and (vii) providing advanced methods for coupling plasma physics codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A GPU ‐Accelerated 3D Unstructured Mesh Based Particle Tracking Code for Multi‐Species Impurity Transport Simulation in Fusion Tokamaks

ABSTRACT This paper presents the multi‐species global impurity transport capability developed in a GPU‐accelerated fully 3D unstructured mesh‐based code, GITRm, to simultaneously track multiple impurity species and handle interactions of these impurities with mixed‐material surfaces. Different computational approaches to model particle‐surface interaction or surface response have been developed and compared. Sheath electric field is taken into account by employing a fast distance‐to‐boundary calculation, which is carried out in parallel on distributed or partitioned meshes on multiple GPUs without the need for any inter‐process communication during the simulation. Several example cases, including two for the DIII‐D tokamak, that is, one with the SAS‐V divertor and the other with the collector probes, are used to demonstrate the utility of the current multi‐species capability. For the DIII‐D probe case, the capability of GITRm to resolve the spatial distribution of particles in localized regions, such as diagnostic probes, within non‐axisymmetric tokamak geometries is demonstrated. These simulations involve up to 320 million particles and utilize up to 48 GPUs.

Nath, Dhyanjyoti D. [Scientific Computation Resear↗

The changing mouse embryo transcriptome at whole tissue and single-cell resolution

During mammalian embryogenesis, differential gene expression gradually builds the identity and complexity of each tissue and organ system 1 . Here we systematically quantified mouse polyA-RNA from day 10.5 of embryonic development to birth, sampling 17 tissues and organs. The resulting developmental transcriptome is globally structured by dynamic cytodifferentiation, body-axis and cell-proliferation gene sets that were further characterized by the transcription factor motif codes of their promoters. We decomposed the tissue-level transcriptome using single-cell RNA-seq (sequencing of RNA reverse transcribed into cDNA) and found that neurogenesis and haematopoiesis dominate at both the gene and cellular levels, jointly accounting for one-third of differential gene expression and more than 40% of identified cell types. By integrating promoter sequence motifs with companion ENCODE epigenomic profiles, we identified a prominent promoter de-repression mechanism in neuronal expression clusters that was attributable to known and novel repressors. Focusing on the developing limb, single-cell RNA data identified 25 candidate cell types that included progenitor and differentiating states with computationally inferred lineage relationships. We extracted cell-type transcription factor networks and complementary sets of candidate enhancer elements by using single-cell RNA-seq to decompose integrative cis -element (IDEAS) models that were derived from whole-tissue epigenome chromatin data. These ENCODE reference data, computed network components and IDEAS chromatin segmentations are companion resources to the matching epigenomic developmental matrix, and are available for researchers to further mine and integrate.

59 BASIC BIOLOGICAL SCIENCES↗

Atomic forces from Dirac-Kohn-Sham equations: Implementation in flexible (APW+lo/LAPW)+LO basis set

The atomic forces formulation based on the Dirac-Kohn-Sham equation and the flexible (APW+lo/LAPW)+LO basis set is presented here. The formulation was implemented in the code FlapwMBPT and allows a user to easily switch between different basis functions of the augmentation type (APW or LAPW) and between different kind of local orbitals. Similar to the work (Phys.Rev.B 91 (2015) 035105), the implementation takes into account small discontinuities of the wave functions, density, and potential at the muffin-tin sphere boundaries. Applications to the materials with strong relativistic effects, such as $\alpha$-Uranium, PuCoGa$_{5}$, and FePt, demonstrate robustness of the method. The comparison of the calculated forces with the ones obtained by numerical differentiation of free energy shows close agreement with deviations about 0.1% or less.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Transfer functions for Q A /Q B international regulatory limits for the safe transport of radioactive materials

This paper presents a proposed revision of the International Atomic Energy Agency transport regulations, related to the A 1 and A 2 limit values used to determine the radioactive transport classification. Based on the 'Q system', a novel methodology was introduced to derive Q A and Q B values related to scenarios involving external exposure from a distant source. These values are key parameters that respectively represent the total effective dose and total equivalent dose to the skin, from all primary and secondary particles contributing to radiation exposure. The International Working Group (WG A 1 /A 2 ) is established and associated with the TRANSSC Technical Expert Group on Radiation Protection. A review of the A 1 and A 2 values is performed in response to identified limitations within the existing Q system. The followed approach is based on Monte Carlo simulations that enabled the development of transfer functions aimed at reducing computational time and increasing the flexibility of dose evaluations for any radionuclide with known particle emission spectra. This method allows updating the Q A and Q B values to account for future data evolutions (decay data, fluence-to-dose conversion coefficients) and standardizing the calculation of regulation limits across all referenced radionuclides and scenarios related to external exposure. The transfer functions are established using three Monte Carlo simulation codes—FLUKA, Geant4, and MCNP—and address the previous limitations of the 'Q system', reflecting the latest International Commission for Radiation Protection recommendations and improvements in calculation techniques. The results of the WG show consistent agreement across the codes, with minor discrepancies observed at low primary energies due to statistical uncertainties and different handling of stopping power for electrons/positrons in the codes. This revised approach aligns with current standards and recommendations, ensuring that the radiological consequences of transport accidents are acceptable for the new A 1 and A 2 limits from a radiological protection perspective.

61 RADIATION PROTECTION AND DOSIMETRY↗

Monte Carlo Driven MDI Optimization at a Muon Collider

A Muon Collider represents a very interesting possibility for a future machine to explore the energy frontier in particle physics. However, to reach the needed luminosity, beam intensities of the order of 10⁹–10¹² muons per bunch are needed. In this context, the Beam-Induced Background must be taken into account for its effects on magnets and detector. Several mitigation strategies can however be conceived. In this view, it is of crucial importance to develop a flexible tool that allows to easily reconstruct the machine geometry in a Monte Carlo code, allowing to simulate in detail the interaction of muon decay products in the machine, while being able to change the machine optics itself to find the best configuration. In this contribution, a possible approach to such a purpose is presented, based on FLUKA for the Monte Carlo simulation and on LineBuilder for the geometry reconstruction. Results based on the 1.5 TeV machine optics developed by the MAP collaboration are discussed, as well as a first approach to possible mitigation strategies.

43 PARTICLE ACCELERATORS↗

Radiation Damage Calculation Methodology

This report outlines the consensus radiation damage calculation methodology used by the Neutronics Analysis Group at Idaho National Laboratory. This methodology includes calculation of radiation damage rates via the Monte Carlo N-Particle (MCNP) code, and calculation of radiation damage in mixed materials, such as alloys. A consistent method for calculating the average displacement threshold energy of a mixed material has been established by the Neutronics Analysis Group. No standard method currently exists to account for radiation damage caused by neutron-induced transmutations. Only one documented case addresses transmutation radiation damage in materials containing nickel, and that method was incorporated by the Neutronics Analysis Group into their radiation damage calculation methodology.

36 MATERIALS SCIENCE↗

HFIR Activity Workbook Generator (HAWK) User Guide

The HFIR Activity WorkbooK generator (HAWK) is a Python code that automates and streamlines the activity calculation of samples after irradiation in the High Flux Isotope Reactor (HFIR). HAWK’s results provide estimates of the activity and nuclide inventory of irradiated specimens before they are moved to hot cell facilities, where they undergo post-irradiation examination. The samples’ activity results guide the packing of shipping containers and inform the accountable inventories for the hot cell facilities. The toolkit was originally developed by Charles Daily, a former R&D staff member at Oak Ridge National Laboratory (ORNL). As of May 2025, HAWK is developed by the Radiation Transport & HPC Methods Group (Nuclear Energy and Fuel Cycle Division) at ORNL. Figure 1 presents HAWK’s workflow. To use HAWK, users need to: 1. Develop an Excel input workbook (i.e., XLSX extension) containing data from the experiment’s materials, irradiation history (cycles), and irradiation positions. 2. Make minor edits to an existing template JSON file (i.e., auxiliary_data.JSON) and to the Python driver. The driver sets the necessary environment variables, defines the material compositions, and ultimately calls HAWK. Once configured, HAWK runs the Oak Ridge Isotope Generation code (ORIGEN) to calculate the masses, activities, and heat load at the end of irradiation for each isotope in the specimen. ORIGEN is part of SCALE, ORNL’s in-house computational tool for performing nuclear safety and design calculations. Following this step, HAWK postprocesses the results and generates three output workbooks summarizing the activity calculations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Excitation of fishbone-like mode in tokamaks due to bounce resonances of trapped thermal ions

When the drift kinetic effect of thermal ions is taken into account, a high-frequency fishbone-like mode (FLM) is found to be driven unstable by trapped thermal ions (TTIs) in tokamak plasmas, according to self-consistent magneto-hydrodynamic (MHD)-kinetic hybrid modeling utilizing the MARS-K code [Liu et al., Phys. Plasmas 15, 112503 (2008)] as well as an analytic theory. It is found that, similar to energetic particles, TTIs can also stabilize the internal kink mode, whereas the FLM is excited when the effective beta of TTIs exceeds a threshold value. The real frequency of the FLM is comparable to the bounce frequency of TTIs. The mode structure of the FLM can be significantly different from the conventional step-like function for the associated plasma radial displacement. This drift kinetic induced modification of the mode structure near the q = 1 surface is captured by non-perturbative MHD-kinetic hybrid computations with MARS-K. Furthermore, the FLM can only be triggered by TTIs at sufficiently high thermal temperatures. Both the FLM and the internal kink can be stabilized by sufficiently fast plasma toroidal rotation and parallel sound wave damping. These two conditions of high thermal temperature and (fast) flow stabilization, though making it challenging to observe the TTI-driven FLM in present day experiments, are favorable for the mode excitation in future reactor scale devices.

36 MATERIALS SCIENCE↗

Multi-scale fission product release model with comparison to AGR data

TRistructural ISOtropic (TRISO) particle fuel is central to several advanced, high-temperature reactor designs. Each particle consists of a fuel kernel encapsulated by three layers of carbon and ceramics that prevent the release of fission products and ensure physical integrity. Despite outstanding retention properties, fission product release has been observed from intact particles. To better understand and quantify fission product release from TRISO particles, a multiscale, mechanistic model of fission product transport is being developed by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. Previous work focused on silver (Ag) transport and improved Ag release predictions. The work described in this report builds on this experience to better understand cesium (Cs) transport in silicon carbide (SiC), the main barrier to the release of fission products. Atomistic simulations provide bulk and grain boundary (GB) Cs diffusivities in SiC, which are used by phase field simulations in the mesoscale code Marmot to determine the temperature, microstructure, and irradiation-dependent Cs diffusivity at the mesoscale in SiC. This approach attributes the different temperature regimes experimentally observed for Cs diffusivities in SiC to a transition from bulk-dominated diffusivity at high temperatures to a GB-dominated regime at low temperatures, providing new insight. The multiscale, mechanistic effective diffusivity is then implemented in the fuel performance code BISON and further validated by comparing Cs release predictions from Advanced Gas Reactor (AGR)-1 and AGR-2 post-irradiation measurements. The new model improves BISON’s predictability. This document also reports improvements made on Ag transport modeling by accounting for different GB types having different diffusivities. Moreover, this report details preliminary efforts to model palladium (Pd) attack of the SiC at the mesoscale using a phase field approach. Pd attack and its impact on accelerated Ag transport remains a misunderstood phenomenon, and we use the model to demonstrate that the formation of lamellae that has been observed in experiments can be explained by the reaction of Pd with SiC to form alternating layers of graphite and Pd 2 Si. This effort aims to improve our understanding of the reaction and eventually provide a model for BISON to account for Pd penetration and its effects on fission product release.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling of High-Temperature Corrosion of Zirconium Alloys Using the eXtended Finite Element Method (X-FEM)

Oxidation modeling in modern nuclear fuel performance codes is currently limited by the lack of coupling with mechanics, thus preventing proper description of how high-temperature oxidation impacts mechanical properties. This is mostly due to the fact that the finite difference formalism adopted in corrosion models is incompatible with the direct coupling with mechanics in the finite element modeling employed in modern nuclear fuel performance codes. In this study, a physically based zirconium alloy corrosion model called the Coupled-Current Charge Compensation (C4) model, which was initially developed for operating temperature conditions, has been updated to include high-temperature corrosion in order to provide additional critical information (e.g., oxygen concentration profile) under loss-of-coolant accident (LOCA) conditions—information lacking in existing empirical models. The C4 model was implemented in the MOOSE finite-element framework developed at Idaho National Laboratory, enabling it to be used in the BISON nuclear fuel performance code based on the MOOSE framework. To precisely track the different interfaces at a relatively low computational cost, the eXtended Finite Element Method (X-FEM) was applied in MOOSE. The model’s results were compared to those of existing empirical models as well as metallographic analysis of high-temperature oxidized Zircaloy-4 coupons. Oxygen diffusivities in the a and ß phases resulting from this comparison closely agree with those found in the literature. The C4 model implemented with X-FEM in MOOSE now has the capability to accurately predict oxide, oxygen-stabilized a, and prior ß phase layer growth kinetics under isothermal exposure at high temperature (1000–1500°C). Furthermore, in contrast with the empirical models, the C4 model accounts for the finite thickness of the fuel cladding. It can predict the oxygen concentration profile evolution through the whole cladding, enabling evaluation of the remaining ductile thickness—a crucial variable for modeling the mechanical behavior of the fuel cladding under LOCA. Furthermore, this implementation allows direct coupling with mechanics, at a low computing cost, using finite-element-based nuclear fuel performance codes such as BISON.

36 MATERIALS SCIENCE↗

Comprehensive Margin Assessment of the ASME Section III, Division 5, Class A Primary Load Design Rules

This report provides a comprehensive margin assessment of the ASME Section III, Division 5, Subsection HB, Subpart B rules for the design of high temperature reactor components against load controlled stress limits. These rules, often called the Code primary load design rules, provide protection against creep rupture and plastic collapse under steady conditions. The method adopted here is to compare the deterministic ASME design life, according to the primary load design rules, against the expected, statistical service life of the component, considering creep rupture as the relevant failure mode. The report provides the design margin of the ASME rules in terms of the probability of premature failure -- the probability that the component fails in service before reaching the ASME design life. To complete this assessment, the report describes the development of a statistical creep life assessment procedure accounting for variations in the component loading, material creep rate, and the material rupture time. This includes a novel method for correlating creep rupture data using a Gaussian process accounting for heat-to-heat variation in the Larson-Miller time-temperature parameter. The report applies the complete margin assessment process to quantify the design margin inherent in the ASME primary load design rules as a function of temperature for relevant component geometries and materials. The overall conclusion is that the ASME rules are conservative, providing designs with very small probability of premature failure, but some optimization of the Code design rules is possible to provide a more uniform design margin.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Design of a prototypical natural circulation water-based reactor cavity cooling system (RCCS) for a pebble-bed generic FHR

A prototypical natural circulation water-based reactor cavity cooling system (RCCS) for the UC Berkeley Mark-1 pebble-bed gFHR is designed based on one-dimensional thermal hydraulics modeling and optimization implemented in an in-house MATLAB code. The model employs a lumped core consisting of fuel pebble and graphite pebble regions, while the graphite reflector, vessel, and RCCS are represented as separate but energy-coupled regions. The model is derived based on steady state energy balance equations accounting for conductive heat transfer from the lumped core to the vessel and subsequent radiative heat transfer to the RCCS and convective heat transfer in the water. Further, mass flow rate is calculated based on momentum and thermal energy balance in the RCCS. Starting with ANL’s water NSTF as baseline, effects of RCCS design parameters such as size of the plate between pipes, pipe diameter, source-to-sink distance, surface emissivity, and pitch from the core are examined in an integral effects framework which accounts for subsequent variation of the mass flow rate, number of pipes, convection coefficient, mean region and surface temperatures, radiative view factors, and water outlet temperature. The importance of including a prototypical reactor design in RCCS design calculations is emphasized as the calculations show that temperature drops from the core to the RCCS fluid to transfer the heat tend to be the limitation not the heat removal capacity of the fluid. A closed-loop design is obtained with capability to safely remove up to 0.72% of the nominal reactor power at maximum estimated peak conditions and 0.40% at shutdown. The design is based on physics calculations and does not account for economic optimization.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Cosmological parameter estimation with a joint-likelihood analysis of the cosmic microwave background and big bang nucleosynthesis

Here, we present a joint-likelihood analysis of big bang nucleosynthesis (BBN) and cosmic microwave background (CMB) data, consistently combining likelihoods and taking into account uncertainties in nuclear reaction rates for the first time. Bayesian inference is performed on the baryon abundance and the effective number of neutrino species, 𝑁 eff , using a CMB Boltzmann solver in combination with LINX , a new flexible and efficient BBN code. We marginalize over Planck nuisance parameters and nuclear rates to find 𝑁 eff =3.0⁢8$^{+0.15}_{−0.14}$, 2.9⁢4$^{+0.16}_{−0.15}$, or 2.96$^{+0.13}_{−0.14}$, for three separate reaction networks. This framework enables robust testing of the lambda cold dark matter paradigm and its variants with CMB and BBN data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A New CIERRA Gridded Product Climatology for LIS/OTD

The GLM-CIERRA processing generates two products: A level-2 cluster feature dataset that merges flashes that should have been a single flash following the clustering algorithm flash definition (for GLM, this is due to the LCFA thresholds); A set of level-3 gridded products including Flash Extent Density, Convective Probability, etc. The CIERRA Level-3 products are generated on a 0.1 degree grid. This is suboptimal given the spatial variations in GLM pixel size AND the size of gridpoints in the output grid. As a stepping-stone towards developing improved GLM-CIERRA grids, I am testing pixel matching techniques on the LIS/OTD data. The goal is to maintain accuracy while vastly improving computational efficiency. The CIERRA reclustering codes are also being applied to LIS / OTD to merge flashes split by the “first fit” clustering technique. The end result will be standardized grids at the nominal resolution of each instrument (5 km for LIS, 10 km for OTD) that take into account the pixel geometry for each event that comprises a given reclustered flash.

54 ENVIRONMENTAL SCIENCES↗

Developing Upscaling Approach for Swarming Hydraulic Fractures Observed at Hydraulic Fracturing Test Site through Multiscale Simulations

We show that this work aims to address a challenge posed by recent observations of tightly spaced hydraulic fractures in core samples from the hydraulic fracturing test site (HFTS) in the Middle Wolfcamp Formation. Many fractures in retrieved cores have subfoot spacing, which is at odds with conventional models in which usually one hydraulic fracture is initiated per cluster. Models assuming a single fracture at each cluster, although a common practice, often predict excessive fracture propagation that is inconsistent with microseismic observation. Here, we aim to develop a numerical approach to effectively account for densely spaced hydraulic fractures in field-scale simulations. Because it is impractical to explicitly model all aforementioned fractures, we develop a new upscaling law that enables existing simulation tools to predict reservoir response to fracture swarms. The upscaling law is derived based on an energy equivalence argument and validated through multiscale simulations using a high-fidelity code, GEOS. The swarming fractures are first modeled with a spacing that is much smaller than the cluster spacing; these fractures are then approximated by an upscaled, single fracture based on the proposed upscaling law. The upscaled fracture is shown to successfully match the energy input rate and produce the total fracture aperture and average propagation length of the explicitly simulated swarm. Afterward, the upscaling approach is further implemented in 3D field-scale simulations and validated against the HFTS microseismic data of a horizontal well. Our results show that hydraulic fracture swarming can significantly affect fracture propagation behaviors compared with the propagation of single fractures as assumed by conventional modeling approaches. Under the considered situations, the conventional treatment yields fast propagation speed that far exceeds that indicated by the microseismic data. We also illustrate that this discrepancy can be reduced readily through the implementation of the upscaling law. Our results demonstrate the importance of accounting for the fracture swarming effect in field-scale simulations and the efficacy of this approach to enable realistic predictions of reservoir responses to fracture swarms, without the need to model tightly spaced fractures individually.

02 PETROLEUM↗

Demonstrate new plasticity models for doped UO 2 that capture dislocation mechanisms

In light water reactors, fuel vendors are investigating the use of dopants to modify the properties of UO 2 pellets, with the goal of improving pellet-cladding mechanical interactions during operation. Dopants are expected to ‘soften’ the pellets; that is, the doped pellets have higher plastic deformation than conventional UO 2 . This leads to a reduction in the severity of mechanical pellet-cladding interactions, helping to reduce the hoop strain on the cladding. By minimizing the strain exerted by the pellet on the cladding, it is anticipated that cladding performance under accident conditions can be enhanced (i.e., lowering the risk of burst during a LOCA). Dopants such as chromium (Cr) promote grain growth during pellet fabrication, leading to larger grains; therefore, understanding the link between chemistry, microstructure and mechanical deformation (enhanced creep rates) behavior of UO 2 is critical to helping operators further substantiate the benefits of doping UO 2 . Historically, the nuclear energy industry has relied on empirical models to make assessments of performance. Compared to empirical models, mechanistic physics-based models provide benefits, such as, fewer data points for validation and better extrapolation where experimental data is scarce or non-existent. In this report, Bayesian inference techniques have been applied to a previously developed lower length-scale-informed diffusional creep model. The objective is to i) infer lower-length-scale parameter distributions from available experiment and then ii) determine the uncertainties in the measurable quantity (in this case creep rates) after propagating the inferred lower length scale parameter uncertainties. The approach requires many evaluations of the model, which becomes computationally insurmountable; therefore, a neural-network model is trained to data obtained by sampling the full model over the most important parameters. This neural-network is then used in the Bayesian inference approach to determine probability distributions in the parameter values that represent the uncertainty in the model given what is known from the experiments (posterior). A significant reduction compared to conservative initial (prior) uncertainties is achieved through inference against the experimental data, demonstrating the efficacy of this approach. Furthermore, by accounting for uncertainties in the experimental conditions and sample non-stoichiometry, it is possible to resolve apparent discrepancies in experimental measurements within a self-consistent grain boundary (Coble) creep model that is sensitive to chemistry. This work has been written up and submitted to Nuclear Technology for a special issue on accelerated fuel qualification (AFQ). This uncertainty quantification (UQ) work not only improves the diffusional model, while accounting for uncertainty, but also establishes a framework which can readily be applied to the mechanistic models of dislocation deformation developed in this study. The most likely values from the Bayesian analysis are incorporated into our UO 2 diffusional creep model and a lower length scale-informed irradiation UO 2 creep mechanistic model to generate a dataset. This dataset has been provided to our INL collaborators for training an artificial neural network surrogate model, which will be implemented in the BISON fuel performance code to assess how the results differ from those currently obtained using a fully empirical model and that of using the nominal (uncalibrated) atomic scale parameters in our mechanistic model. Plastic deformation (creep and glide) in UO 2 is a complex phenomenon, governed by multiple underlying processes such as local defect concentrations, applied stresses, and microstructural characteristics. Consequently, there is a need for a meso-scale model with polycrystalline resolution capable of extrapolating to large grain sizes applicable to doped UO 2 , where data is limited and the model can help bridge the knowledge gap. By integrating atomistic data into the polycrystal LApx code, it becomes possible to predict dislocation climb and glide plasticity that simple analytical models cannot accurately represent. The application of atomic-scale data within LApx demonstrated the importance of climb and glide mechanisms in reproducing high-stress UO 2 behavior. Behaviors such as this are crucial to capture and implement in BISON, as parts of the fuel pellet can reach temperatures where glide can occur before pellet cracking. This model which captures dislocation based mechanisms for UO 2 is then used to stand up the doped model accounting for larger grain sizes. It was found that larger grain sizes can lead to enhanced deformation rates in the glide regime, and therefore can help with the pellet cladding mechanical interaction. Therefore if the fuel pellet reaches conditions (stress/temperature) where glide is active, the enhanced creep rates for larger grains in the glide regime (doped UO 2 ) can help with pellet cladding mechanical interactions. Plastic deformation in UO 2 involves multiple mechanisms, including diffusional creep, dislocation climb, and glide. This milestone contains two parts: (1) UQ of a pre-existing lower length scale informed mechanistic diffusional creep model, and (2) development of a new LApx based model for dislocation-mediated creep mechanisms in UO 2 , with application to large-grain doped UO 2 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Studies of the influence of heat pipe model and coupling strategy in heat pipe microreactor simulations using Sockeye

This work demonstrates various modeling approaches for thermally coupling heat pipes to a graphite monolithic block in a microreactor context, accounting for the presence of a gap. We show that heat pipe model selection has a significant impact on the transient behavior of a three-dimensional coupled microreactor assembly problem by comparing results obtained using different heat pipe models from the heat pipe code Sockeye. Additionally, we show that the choice of coupling strategy has significant consequences for both accuracy and performance. Furthermore, we found that our heat flux transfer strategy exhibited greater accuracy than our temperature transfer strategy, even when comparing a loose coupling of the heat flux transfer strategy to a tight coupling of the temperature transfer strategy. Although this research uses only theoretical test problems, it provides important insights in modeling thermal fluids phenomena in heat pipe microreactors.

Direwolf↗