Application of physics-based pulverization models in FFRD simulations [Slides]
Summary presentation on application of physics-based models in FFRD analysis for EPRI/INL/DOE workshop. Contains content already cleared through LRS.
SEARCH · Search NASA
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.
Summary presentation on application of physics-based models in FFRD analysis for EPRI/INL/DOE workshop. Contains content already cleared through LRS.
Presentation summarizing work done by INL and ORNL investigating loss-of-coolant accidents involving high burnup fuel. This includes developing LOCA test conditions for TREAT and SATS experiments and highlighting potential differences between nuclear and electrical heating LOCA experiments.
Six reactivity-initiated-accident (RIA) commissioning tests have been performed at Idaho National Laboratory’s (INL) Transient Reactor Test Facility. Five of these tests were performed using fresh fuel rodlets, and the sixth test was performed using a previously irradiated rodlet from the ATF-2 irradiation experiment. These experiments demonstrate the ability to perform RIA testing including the ability to perform experiments using previously irradiated materials. Post-transient examinations revealed pellet cladding interactions and cladding ballooning depending on the initial boundary conditions of the test. Increased thickness of both zirconium oxide and alpha-zirconium were found in the cladding with tests performed at higher levels of total energy deposition. No oxide or alpha-zirconium was found in the test with the greatest amount of ballooning, which was also the test with the highest initial rod internal pressure. The highest energy deposition test (1110 J/g) resulted in failure of the rodlet. Fragmentation of the fuel pellets occurred in this test with the size of the fragments decreasing inversely to the pellet radius.
Title (20 words): Investigation of the Collision-Induced Dissociation Mechanism of Protonated TODGA with IRIS Introduction (120 words): One of the challenges facing wide-spread adoption of nuclear power is the development of efficient separation processes for used nuclear fuel. The molecules in separation processes are subjected to an extreme environment due to the high radiation fields from the used fuel and highly acidic media used for fuel dissolution, which results in significant molecular degradation, leading to reduced process efficiency. These degradation products must be identified and studied so mitigation strategies can be developed to maintain process efficiency. However, complex systems can have many degradation products, complicating identification. Untargeted analysis tools could be used to understand radiation chemistry in complex systems. However, this would necessitate improved understanding of the gas-phase fragmentation mechanisms of fuel cycle molecules like tetraoctyldiglycolamide (TODGA). Methods (120 words): The gas-phase fragmentation of protonated TODGA was investigated using collision-induced dissociation (CID), resonance ejection, and infrared ion spectroscopy (IRIS). CID and resonance ejection experiments were conducted using a Bruker Daltonics (Bremen, Gemany) SolariX XR fourier transform ion cyclotron resonance (FT-ICR) mass spectrometer. IRIS spectra of protonated TODGA and its two CID fragmentation products were measured using a modified Bruker amaZon Speed ETD 3D quadrupole ion trap mass spectrometer coupled to the Free Electron Lasers for Infrared eXperiments (FELIX) free electron laser. Measured spectra were compared with density functional theory (DFT) calculations using the Gaussian 16, Revision C.02 software package with the ?B97X-D functional and def2-TZVPP basis sets. Candidate structures were generated using the CREST 3.0 conformational sampling software tool. Preliminary Data (300 words): Collision-induced dissociation of protonated TODGA ([C36H73N2O3]+, m/z=581.562) results two fragment ions, one at m/z=340.285 assigned as [C20H38NO3]+ and the other at m/z=312.290, assigned as [C19H38NO2]+. Based on the assigned formula and the structure of protonated TODGA, the fragment at m/z=340.285 is likely formed from elimination of neutral dioctylamine. Comparison of the IRIS spectrum of m/z=340.285 with DFT predictions suggests it contains a ring structure, and is assigned as N-octyl-N-(6-oxo-1,4-dioxan-2-ylidene)octan-1-aminium. Based on this structure and the structure of protonated TODGA, we hypothesize this fragment formed from elimination of neutral dioctylamine followed by a ring closure mechanism. Comparison of the IRIS spectrum of the fragment at m/z=312.290 with DFT predictions also indicated the presence of a ring structure, assigned as N-(1,3-dioxolan-4-ylidene)-N-octyloctan-1-aminium. This product could be formed from elimination of carbon monoxide from the ring of m/z=340.285 as a sequential fragmentation or formed directly from protonated TODGA via elimination of neutral N,N-dioctylformamide followed by a ring closure. Resonance ejection experiments where m/z=340 was continuously ejected from the IRC cell showed no decrease in intensity of m/z=312.290 across several collision energies, suggesting that the later, direct formation mechanism, dominates. The location of the ionizing proton in protonated TODGA is important for modeling the fragmentation mechanisms. DFT calculations suggested that the position of bands involving the coupled vibrations of the amide C—N and C=O bonds in TODGA are the most sensitive to proton location. Evaluation of the IRIS spectrum of protonated TODGA suggests that the ionizing proton is located between the two amid oxygens. This protonation location was calculated to lie approximately 30 kJ/mol lower in energy than the next lowest energy location, with the proton located solely on one of the amide oxygens. Novel aspect (20 words): Infrared ion spectroscopy combined with resonance ejection experiments and density functional theory to probe the collision-induced dissociation mechanism of tetraoctyldiglycolamide.
Studies on high burnup UO 2 subjected to loss-of-coolant accident conditions have shown that restructured regions of the fuel are susceptible to pulverization and eventual dispersal. Due to a lack of pre-test characterization, the distinct microstructural features rendering the fuel prone to fragmentation remain ambiguous. Four samples of commercially irradiated light-water reactor UO 2 have been characterized utilizing electron backscatter diffraction to assess the susceptible microstructure. The microscopy focused on determining the burnup and temperature conditions responsible for the formation of the different microstructural regions where the regions were denoted as the high-burnup structure (HBS), HBS transition, mid-radial, restructured central, and central region. Previous works have outlined the specific conditions required for the restructuring of the microstructure into the HBS, but the conditions responsible for the restructuring in the central region of the fuel are not well understood. The four analyzed samples confirm a burnup threshold of 61 GWd/tU, and an unknown temperature range is needed to facilitate the formation of the restructured central region. In conclusion, additional fuel performance evaluations are needed to quantify the temperature range promoting restructuring in the central region.
The work described in this report has developed a preliminary conceptual flowsheet based on limited published literature for the head-end of a TRISO oxycarbide UNF reprocessing plant. The essential function of the head-end is to prepare the UNF for subsequent separation of actinides. The main steps of the conceptual TRISO UNF head-end reprocessing flowsheet are: 1. Fragmentation of the graphite moderator to expose the fuel particles. 2. Separation of fuel particles from fragmented graphite moderator by sieving and fluidization. 3. Fragmentation of the fuel particles to expose the fuel kernel. 4. Dissolution of the fuel kernel in nitric acid. Noteworthy at this step is the formation of organic acids (mainly oxalic and mellitic acids) from the carbon associated with the fuel and TRISO coatings. These acids can interfere with actinide separations and solvent extraction hydraulic performance and carbon dioxide gas generation. 5. Clarification of the dissolved fuel solution to separate the coating fragments. 6. Dissolver off-gas treatment with specific emphasis on managing generation of carbon dioxide gas containing carbon-14. The quantity of carbon dioxide is projected to be 100 to 1000 times greater than that arising from dissolving LWR UNF. Pre or post treatment of the exposed fuel kernels could be undertaken to mitigate the formation of organic acids but there still exists a significant off-gas treatment challenge. The literature on reprocessing uranium nitride UNF was reviewed but a flowsheet was not specifically developed since its scope and uncertainties are expected to be encompassed by the oxycarbide flowsheet. Overall, the predominant feature of processing TRISO UNF is carbon management to reduce waste volume and mitigate the formation of carbon dioxide and organic acids. This report recommends broad areas of research that will help to further define the head-end flowsheet.
Safety/transient testing to evaluate performance under off-normal conditions is an essential pillar for both the development of Accident Tolerant Fuels (ATF) and the optimization of fuel operation economics beyond current discharge burnups. Among other factors, the successful interpretation of the transient testing results relies upon the knowledge of the initial conditions of the test, including the characteristics of the fuel system under scrutiny. When testing pre-irradiated material, the assumptions that the fuel and the cladding still have the same properties as in the pre-irradiation stage is obviously wrong and could affect the results of the test. This is particularly true the more burnup accumulates in the fuel rod and irradiation progresses. The knowledge of the initial microstructure of both fuel and cladding allows a clearer interpretation of the subsequent transient testing results, provides validation of the physical phenomena underlying the model predictions and eliminates the uncertainties related to the limited knowledge of the sample status before the test. One example is the phenomenon of fine fragmentation that occurs in Light Water reactor (LWR) fuel. During a Loss of Coolant Accident (LOCA) or Reactivity Initiated Transient (RIA) the fuel can severely fragment. During LOCA, high burnup fuel tend to finely fragment, which has raised safety concerns due to the increased likelihood of dispersal of such small particles once the cladding has burst and due to the increased fission gas release. Therefore, efforts have been devoted to the assessment of a pulverization threshold that could determine the conditions under which fine fragmentation is predominant. However, the lack of information regarding the initial conditions of the fuel, and the connections between those conditions and the pre-transient irradiation history, have hindered the development of a fully mechanistic fragmentation and pulverization criterion. The empirical relationships rely on conservative estimations, due to the lack of information on critical material properties and characteristics. More generally, experimental evidence of the irradiation-induced modifications at microstructural scale are necessary to determine the behavior of the material at the macroscopic scale, with the latter being the one of technological interest. Significant progress has been made in the last two decades in the developments of analytical materials science techniques that can be applied to highly radioactive materials, such as high burnup fuels. The availability of new techniques and the improvement of existing ones has enabled investigations previously not possible that can deepen the understanding of the fuel characteristics and properties at high burnup. The better knowledge of material behavior and irradiation-induced phenomena could help the prediction of its performance. In this context, the scope of the present work is to apply a wide portfolio of advanced characterization techniques to determine properties that are relevant for safety and performance. The results are interpreted in the context of engineering scale post-irradiation examinations and available information on the irradiation conditions.
Biorefinery processes for converting lignocellulosic biomass to fuels and chemicals proceed via an integrated series of steps. Biomass is first pretreated and deconstructed using chemical catalysts and/or enzymes to liberate sugar monomers and lignin fragments. Deconstruction is followed by a conversion step in which engineered host organisms assimilate the released sugar monomers and lignin fragments, and produce value-added fuels and chemicals. Over the past couple of decades, a significant amount of work has been done to develop innovative biomass deconstruction and conversion processes that efficiently solubilize biomass, separate lignin from the biomass, maximize yields of bioavailable sugars and lignin fragments and convert the majority of these carbon sources into fuels, commodity chemicals, and materials. Herein, we advocate that advanced in silico approaches provide a theoretical framework for developing efficient processes for lignocellulosic biomass valorization and maximizing yields of sugars and lignin fragments during deconstruction and fuel and chemical titers during conversion. This manuscript surveys the latest developments in lignocellulosic biomass valorization with special attention given to highlighting computational approaches used in process optimization for lignocellulose pretreatment; enzyme engineering for enhanced saccharification and delignification; and prediction of the genome modification necessary for desired pathway fine-tuning to upgrade products from biomass deconstruction into value-added products. Physics-based modeling approaches such as density functional theory calculations and molecular dynamics simulations have been most impactful in studies aimed at exploring the molecular level details of solvent-biomass interactions, reaction mechanisms occurring in biomass-solvent systems, and the catalytic mechanisms and engineering of enzymes involved in biomass degradation. More recently, with ever increasing amounts of data from, for example, advanced mutli-omics experiments, machine learning approaches have begun to make important contributions in synthetic biology and optimization of metabolic pathways for production of biofuels and chemicals.
To improve the economics of commercial nuclear energy generation, U.S. utilities are currently seeking licensing approval to operate UO2 fuel to higher burnups. One significant safety issue that must be addressed to obtain approval is the potential for fine fragmentation/pulverization of the fuel during a loss-of-coolant accident (LOCA). The cause of pulverization has been hypothesized to be the rapid increase of pressure in fission gas bubbles in the high burnup region of the fuel during a LOCA. To better understand this phenomenon, a novel phase-field model of the fission gas bubble microstructure in UO2 has been developed and implemented in Idaho National Laboratory (INL)'s Marmot application for phase-field simulation of nuclear materials. Simulations of the bubble response to steady-state and transient conditions were conducted. Simulation results were used to inform a mechanistic model of pulverization in BISON, INL’s fuel performance simulation code.
This report summarized work carried out to investigate mechanisms relevant to fragmentation of high burnup UO 2 fuel. High pressure bubbles that form in high burnup structure are thought to be responsible for fragmentation and pulverization, when exposed to a temperature ramp during a transient. The over-pressurization of the bubbles during normal operation is a pre-requisite for this mechanism to occur. Therefore, in this work we investigate the role of interstitials, produced through irradiation, in over-pressurizing the bubbles by using a combined molecular dynamics (MD) and cluster dynamics approach. Firstly, the energies for the annihilation of interstitials and vacancies at bubbles have been determined from MD as a function of bubble pressure. Secondly, these reaction energies have been implemented in the cluster dynamics code Centipede to determine the steady-state defect concentrations and bubble pressure under irradiation. It was found that there is a transition from low pressure bubbles, at high temperature, to high pressure bubbles, at low temperatures. This indicates the formation of over-pressurized bubbles in the low temperature periphery of the pellet, where high burnup structure forms. This result supports the hypothesis that over-pressurized bubbles form during steady-state operation and can then contribute to fragmentation during a transient. Future work will examine the influence of bubble pressure on Xe diffusion and study the impact of other sinks (e.g. grain boundaries and dislocations) that are relevant to high burnup structure. Ultimately, the goal is to provide predictions of the bubble pressure as a function of burnup that can be used to derive a more mechanistic fragmentation threshold for use in BISON.
This report summarizes lower length scale computational research conducted to improve the pulverization criterion for high-burnup UO 2 fuel in the BISON fuel performance code. This research was sponsored by the NEAMS program during FY22. Efforts to improve the model primarily focus on calculating the current pressure of bubbles in the high-burnup structure (HBS) region, as well as the critical bubble pressure at which pulverization occurs in the HBS. The phase-field model for predicting initial bubble pressure in the HBS was improved by implementing a more realistic model for defect production, and by coupling the phase-field model for inter-granular bubble evolution with the spatially-resolved cluster dynamics code Xolotl to simulate intra-granular fission gas evolution. An evolution equation for current bubble pressure was added to BISON, an improvement over the previous model that assumed that bubble pressure was static following HBS formation. To improve calculations of critical bubble pressure for pulverization, 3-D phase-field fracture simulations were performed, and a function for critical bubble pressure was fit—based on the simulation results—to replace the previous function, which had been determined using 2-D simulations. The impacts of these modifications to the existing pulverization assessment cases will be reported in the forthcoming engineering scale milestone on high-burnup UO 2 pulverization. A poromechanics-based approach was used to include the effect of bubble over pressurization on the stress state of the pellet at the engineering scale, and the initial strategy for integration of the pulverization criterion with pellet-scale mechanical degradation using a smeared cracking model was developed. This report also describes the initial research that was conducted to inform a forthcoming transient fission gas release (tFGR) model in BISON. Focusing on the fission gas release caused by pulverization of the outer pellet rim, a model for the amount of fission gas release was developed as a function of fuel porosity, bubble size, and fragment size. Implementation of this model in BISON and its impact of engineering-scale predictions will also be described in the aforementioned engineering-scale milestone on high-burnup UO 2 pulverization.
Reactivity Initiated Accidents (RIA) are design basis accidents that could have adverse impacts on the core coolability. In the unlikely event that sufficient reactivity is inserted into the reactor core by the ejected/dropped control rod, prompt energy deposition into the fuel can occur, which when sufficiently high can lead to fuel rod failure or, at large energy deposition levels, expulsion of UO2 fragments or molten UO2 material from the fuel rod. This results in a release of fission product and fuel into the coolant and has the potential to compromise core coolability and threatening the pressure boundary of the primary coolant system. The design basis RIA is one of the industry challenging problems that the CASL aimed to address; the CASL RIA Challenge Problem Charter [1] states, ?The Pressurized Water Reactor (PWR) Rod Ejection Accident (REA) and Boiling Water Reactor (BWR) Control Rod Drop Accident (CRDA) are postulated accidents with consequences that are important to nuclear safety (fuel rod integrity and core coolability). Currently each reload core design must be analyzed to meet regulatory acceptance criteria. The goal of CASLs ModSim capability for RIA is to model the event at a higher fidelity, with validation to existing tests, to better model the transient neutronics and the progression of the fuel and cladding thermal-mechanical behavior. These improved analytical capabilities can be used to better inform reload core design, limits on fuel assembly discharge burnup, restrictions on placement of fuel in the reactor, control rod insertion limits, operating margin, and performance sensitivities.? In support of that charter, BISON, the fuel performance code has been used to demonstrate the simulation of thermal-mechanical behavior of LWR fuels during a reactivity-initiated accident. The combination of mechanical, thermal, and thermal-hydraulic phenomena all present during an RIA makes a multi-physics code such as BISON a valuable tool for modeling these scenarios. This paper will highlight many of the activities associated with BISON relevant to reactivity-initiated accident capability development and validation efforts. These efforts have been performed under both CASL and the Nuclear Energy Advanced Modeling and Simulation (NEAMS) programs.
Zirconium carbide (ZrC) is commonly used for energy sector research, as well as a surrogate for the proposed advanced nuclear fuel candidate uranium carbide. Here this study investigates structural modifications to nanocrystalline and microcrystalline ZrC resulting from dense electronic excitations induced by swift heavy ion exposure. Samples were irradiated with 946 MeV Au ions to various fluences up to 6 × 10 13 ions cm –2 and characterized using synchrotron-based x-ray diffraction. The evolution of the unit-cell parameter and heterogeneous microstrain were evaluated as a function of fluence and compared with those of nanocrystalline and microcrystalline CeO 2 (a surrogate for UO 2 fuel) irradiated under identical conditions. Distinct differences were observed in the radiation responses of the carbide and oxide across both grain sizes. Most notably, microcrystalline ZrC exhibits swelling characterized by two distinct regimes, which does not result in saturation at the ion fluences achieved. This contrasts with CeO 2 , which exhibits the well-documented direct-impact defect accumulation mechanism, reaching a steady-state saturation of swelling at higher fluences. Nanocrystalline CeO 2 undergoes more pronounced swelling compared with microcrystalline CeO 2 , in contrast to nanocrystalline ZrC, which exhibits only minimal unit-cell changes. These results demonstrate that swift heavy ion-induced structural changes can be quite different in carbides and oxides, which must be considered when extrapolating fission-fragment type damage in current fuels to advanced fuels.
Radiation of fuel rods with neutrons during a fission reaction in a nuclear reactor is known to affect the material microstructure of the fuel as a result of fission fragment damage, and lead to the formation of atomic-level defects such as voids, dislocation loops, and lattice swelling from fission-gas release. A consequence of the radiation-induced damage in nuclear fuels is the drastic alteration of material properties, in particular, the reduction of the thermal conductivity – a key parameter that governs the transport of thermal energy released from fissile fuel atoms to the surrounding coolant. Swelling can induce high stresses and ultimate failure of the cladding. A fundamental understanding of the role of radiation-induced damage on the thermal transport and mechanical properties of nuclear fuels is therefore critical for efficient, reliable and safe operation of a nuclear power plant. As opposed to conventional post-irradiation examination techniques that are often time-consuming and require significant sample preparation, laser-based characterization methods have emerged as promising non-contact, non-destructive tools to measure the evolution of microstructure-induced material property changes. Moreover, laser methods offer the promise of in-situ characterization while the material is being irradiated. This project aims to develop an improved understanding of the impact of radiation-induced changes to material microstructure on the thermal and elastic properties of nuclear materials. A laser ultrasonic technique, known as the transient grating (TG) spectroscopy method, is used to simultaneously measure thermal diffusivity and elastic properties in unirradiated and ion-irradiated oxide nuclear fuel samples. The results will complement on-going investigations on electron- and phonon-mediated thermal transport in nuclear materials and will provide foundational work for incorporating the influence of defects in fuel performance codes.
The separation and characterization of two irradiated uranium targets, a depleted uranium (DU) and a highly enriched uranium (HEU) target as well as a plutonium (Pu) target, was conducted in April of 2022. The three targets were assembled at Los Alamos National Laboratory (LANL) and irradiated using the GODIVA critical assembly at the National Criticality Experiments Research Center (NCERC). Splits of the dissolved targets were received by Pacific Northwest National Laboratory (PNNL) after which the PNNL and LANL teams chemically separated the solutions using independent separation schemes and analyzed the separated fractions for short lived actinides and fission products. Chemical separations were traced with stable or radioactive tracers to allow for the determination of chemical yields, analyzing using either inductively coupled plasma optical emission spectroscopy (ICP-OES), inductively coupled plasma mass spectrometry (ICP-MS) or gamma emission analysis (GEA) depending on the nature of the tracer. The Pu target solution was traced with stable elements at LANL to follow elemental fractionation during a Pu removal step. Many analytical techniques were used by PNNL including kinetic phosphorescence analysis (KPA), ICP-OES, ICP-MS, GEA, and thermal ionization mass spectrometry (TIMS) depending on the analyte’s need.
This paper presents a TRACE model of a postulated large-break loss-of-coolant accident (LBLOCA) event in a four-loop pressurized water reactor (PWR). The input files have been made publicly available and serve as a starting point for researchers to improve upon or apply to their purposes. The model, based on a legacy PWR model, contains numerous improvements and modifications consistent with US Nuclear Regulatory Commission LBLOCA analysis guidelines. The model is applied to analyze high-burnup (>62 GWD/MTU) fuel behavior during LBLOCA, improving on previous work to provide realistic thermal hydraulic predictions for future fuel performance analyses and experiments related to fission fragment, relocation, and dispersal (FFRD). Sensitivity studies are presented to quantify the impact of the vessel-modeling approach, core radial discretization, and other phenomena. In conclusion, the model is intended for research purposes and contains sufficient plant geometric, operational, and safety system details to represent the main physical phenomena pertaining to LBLOCA.
A 3D numerical model for the ablation of pellets and shattered pellet injection (SPI) fragments in tokamaks in the plasma disruption mitigation and fueling parameter space has been developed based on the Lagrangian particle code [R. Samulyak, X. Wang, H.-S. Chen, Lagrangian Particle Method for Compressible Fluid Dynamics, J. Comput. Phys., 362 (2018), 1-19]. The pellet code implements the low magnetic Reynolds number MHD equations, kinetic models for the electronic heating, a pellet surface ablation model, an equation of state that supports multiple ionization states, radiation, and a model for grad-B drift of the ablated material across the magnetic field. The Lagrangian particle algorithm is highly adaptive, capable of simulating a large number of fragments in 3D while eliminating numerical difficulties of dealing with the tokamak background plasma. The code has achieved good agreement with theory for spherically symmetric ablation flows. Axisymmetric simulations of neon and deuterium pellets in magnetic fields ranging from 1 to 6 Tesla have been compared with previous simulations using the FronTier code, and very good agreement has also been obtained. Furthermore, the main physics contribution of the paper is a detailed study of the influence of 3D effects, in particular grad-B drift, on pellet ablation rates and properties of ablation clouds. Smaller reductions of ablation rates in magnetic fields compared to axially symmetric simulations have been demonstrated because the ablated material is not confined to narrowing channels in the presence of grad-B drift. Contribution of various factors in the grad-B drift model has also been quantified.
Microstructural defects formed as a result of fission fragment damage are known to drastically alter thermal conductivity, a fuel property that governs the efficiency of a nuclear reactor. A fundamental understanding of radiation-induced effects on thermal transport is critical for the development of advanced fuels. Ion-irradiation has widely been used to simulate the effects of neutron-irradiation by seeding atomic-to-nanoscale defects of controllable size and density. In this work, we investigate the influence of microstructural defects, induced by ion-irradiation, on the thermal properties of single crystal thorium dioxide samples grown using the hydrothermal technique. The samples were irradiated using 2 MeV protons at room temperature. A laser-based modulated thermoreflectance technique was used to measure the thermal diffusivity and conductivity within the damaged ThO2 region. The experimentally measured thermal properties were compared with Boltzmann Transport Equation predictions, to determine the impact of various defect types on the thermal properties of ThO2.