Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear capture”

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 991 records · Page 55

Depth-resolved thermal conductivity and damage in swift heavy ion irradiated metal oxides

We investigated thermal transport in swift heavy ion (SHI) irradiated insulating single crystalline oxide materials: yttrium aluminum garnet- Y 3 Al 5 O 12 (YAG), sapphire (Al 2 O 3 ), zinc oxide (ZnO) and magnesium oxide (MgO) irradiated by 167 MeV Xe ions at 10 12 – 10 14 ions/cm 2 fluences. Depth profiling of the thermal transport on nano- and micro- meter scales was assessed by time-domain thermoreflectance (TDTR) and modulated thermoreflectance (MTR) methods, respectively. This combination allowed us to isolate the conductivities of different sub-surface damage-regions characterized by their distinct microstructure evolution regimes. Thermal conductivity degradation in SHI irradiated YAG and Al 2 O 3 is attributed to formation of ion tracks and subsequent amorphization, while in ZnO and MgO it is mostly due to point defects. Additionally, notably lower conductivity when probed by very low penetrating thermal waves is consistent with surface hillock formation. An analytical model based on Klemens-Callaway method for thermal conductivity coupled with a simplified microstructure evolution capturing saturation in defect concentration was used to obtain depth dependent damage across the ion impacted region. The studies showed that YAG has the highest damage profile resulting in the less dependence of thermal conductivity with the depth, while MgO on the contrary has the strongest dependence. The presented work sheds new light on how SHI induced defects affect thermal transport degradation and recovery of oxide ceramics as promising candidates for next generation nuclear reactor applications.

36 MATERIALS SCIENCE↗

Nondestructive Evaluation (NDE) of Cable Anomalies using Frequency Domain Reflectometry (FDR) and Spread Spectrum Time Domain Reflectometry (SSTDR)

This report presents a comparative assessment of the performance of frequency domain reflectometry (FDR) and spread spectrum time domain reflectometry (SSTDR) in detecting a wide range of electrical cable anomalies. All tests and results reported herein were performed at the PNNL Accelerated and Real-Time Environmental Nodal Assessment (ARENA) cable and motor test bed. The primary objective of this work was to evaluate the effectiveness of SSTDR, a fledgling cable monitoring technique that shows promise for application in online monitoring of energized cable systems, against FDR, an offline technique widely employed in the nuclear power plant (NPP) industry. FDR tests are becoming more widely used in nuclear power plant cable aging management and test programs – particularly for low voltage cables. FDR capabilities for these kinds of tests have been reported by PNNL and others. The FDR test is performed on de-energized cables by connecting the FDR instrument to two of the cable conductors, or one conductor and the shield. A broad band low voltage (< 5 V) chirp is introduced in the cable, and any reflected response is captured in the frequency domain. The captured reflection is then processed by performing an inverse Fourier transform to a time domain response which can then be converted to a distance response based on the cable velocity of propagation (VoP). SSTDR measurements are functionally similar to FDR measurements in that a broad-band voltage signal composed of a square or sine wave modulated pseudo-random sequence of chips (< 5 volts), is injected onto one of the cable conductors. The injected signal will experience partial energy reflection and transmission at each impedance discontinuity along the transmission line. Any reflected response is detected by computing a cross-correlation between the reflected signals and a delayed copy of the incident SSTDR signal. the time delay for the reflected signal to experience the best matched correlation with the incident signal, indicates the travel time for the signal to reach a change in impedance. By knowing this time delay and velocity of propagation (VoP) of the signal, one can compute the physical distance. A big advantage that SSTDR measurements have over other methods is the ability to be connected to energized or live wires (currently up to 1kV) thereby enabling online monitoring of cables. SSTDR has been used successfully in several applications, e.g., aircraft, rail, and photovoltaic systems. In this work FDR and SSTDR cable assessment techniques were used to characterize a variety of cable anomalies and faults including: (1) Presence or absence of a motor; (2) Ground faults and short circuit faults; (3) Moist environments and water ingress faults; (4) Accelerated thermal aging. Both shielded and non-shielded cables were evaluated in this report. Offline measurements were made using FDR and online measurements were made by SSTDR for a range of test scenarios. Based on the results across all cable anomalies evaluated in this study, FDR displayed high sensitivity towards cable condition assessment, while SSTDR showed promise for future application in monitoring NPP cable systems. However, further developments are suggested to improve the resolution and sensitivity of SSTDR towards faults and anomalies in low voltage cables.rt presents a comparative

42 ENGINEERING↗

Updates and Validation for the n+ 63,65 Cu Cross Sections [Abstract]

The neutron induced total, elastic, and capture cross sections of 63,65 Cu isotopes were selected for evaluation in the resolved and unresolved resonance energy ranges by the National Criticality Safety Program to resolve discrepancies related to benchmark performance. This is especially evident for the series of ZEUS benchmarks in which copper is used as a reflector. Because copper is also used as structural material in both fission and fusion reactors, the need to address benchmark discrepancies linked to nuclear data deficiencies is a task of primary importance. The aim of this work is to describe the steps of evaluation work towards a consistent improvement of the benchmark performance. The R-matrix analysis with the SAMMY code focused on the 63 Cu(n,γ) reaction channel between 100-300 keV coupled to unresolved resonance region parameters up to 650 keV to fit average cross section data from a recent experiment. Due to the high sensitivity of many benchmarks to elastic scattering angular distribution data, especially for the 65 Cu isotope, the impact of these data was tested by generating Legendre coefficients from both resonance parameters and the Hauser-Feshbach model. Guided by the findings of Shaw et al., the performance of the current evaluation for 65 Cu was compared to that of ENDF/B-VII.1 and ENDF/B-VIII.0 by testing the reactivity coefficients corresponding to the validation suite of experimental criticality benchmarks for thermal, intermediate, and fast systems taken from the International Criticality Safety Benchmark Experiments Project Handbook. The benchmark performance is especially sensitive to 63 Cu(n,γ) and 65 Cu elastic scattering for neutron energies in the 100–500 keV region, whereas 100 keV is the upper limit of the resolved resonance region in the ENDF/B-VIII.0 evaluations for 63,65 Cu. The results highlight the need to handle the transition from the resolved resonance region to the high energy region carefully.

07 ISOTOPE AND RADIATION SOURCES↗

Implementation of a Drift-Flux Model in SAM for Modeling of Passively Transported GAS in Molten Salt Reactors

The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program is further developing the Multiphysics Object Oriented Simulation Environment (MOOSE) system thermal hydraulics (T/H) code, System Analysis Module (SAM), to include modeling of molten salt reactor (MSR) designs. MSR designs in which the fuel is dissolved in the salt coolant poses unique modeling and simulation (M&S) challenges because fission products—some of which will be noncondensable gasses—travel with the coolant, thus allowing for deposition and interaction throughout the entire system flow loop. This behavior must be modeled to successfully understand MSR system behavior during normal and off-normal operation events. In past work, SAM has been extended to include a passive species transport capability for modeling of MSR fission products that travel with the salt coolant. This work extends the model to include a drift-flux model to capture the tendency of the gas phase to move at a different velocity than the bulk liquid. Additional constitutive models have also been added for calculation of gas volume fraction, bubble radius, and interfacial area, which will be important for calculation of species phase migration in future coupling activities. Testing has been performed to ensure that the new model agrees with expected values for several simple models, but additional future work will be required to test the model for more complex geometries that will be needed for full-scale MSR simulations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Physical and Thermomechanical Properties of Yttrium Hydride from Large Scale Bulk Metal Hydriding Furnace

Given the superior thermal stability and highly attainable hydrogen density, yttrium hydride is an excellent high-temperature moderator material in advanced thermal neutron spectrum reactors that require small core volumes. Yttrium hydride has been selected as the moderator material for the Transformational Challenge Reactor, which was launched at Oak Ridge National Laboratory (ORNL) in 2019. However, fabrication of large-scale crack-free yttrium hydride is challenging and very limited efforts have been committed to the characterization of bulk yttrium hydride in response to the need to establish a complete database of the thermomechanical properties of YHx. In this report, the challenges associated with fabricating large-scale crack-free yttrium hydride are discussed herein. In response to those challenges, a hydriding system was designed and constructed at ORNL and was used to successfully fabricate crack-free yttrium hydride in complex geometries at large scales. This was accomplished by precisely controlling the hydrogen’s partial pressure and the retort temperature, which was informed by the well-established thermodynamic properties of the binary H-Y system. Hydrogen content in as-fabricated hydride was determined by the weight change method and vacuum hot extraction technique, complemented by the X-ray diffraction (XRD). In addition, significant efforts are being dedicated to establishing a complete database of the thermomechanical properties of as-fabricated yttrium hydride. In FY2020, we investigated the thermophysical properties of yttrium hydrides as a function of temperature (room temperature to 700°C) and hydrogen concentration (H/Y ratio ranges from 1.52 to 1.93). The results indicate that at the temperatures below 300 °C, the hydrogen content did not have a significant influence on the thermal expansion, while the specific heat capacity, the thermal diffusivity, and the calculated thermal conductivity were slightly higher for the higher H/Y ratio. Between 300°C and 700 °C, a reversible second-order endothermic transition in all measured thermal properties was observed. It was also found that the onset temperatures of the observed transition varied, with the composition having inverse dependence on the hydrogen content. An attempt was made to explain the behavior of the thermophysical properties at higher temperatures by considering the order– disorder transition as a result of hydrogen redistribution. In addition, nanoindentation was employed to determine the elastic modulus and hardness and to capture the crystal orientation dependence of these parameters. Vickers hardness was also reported. The final section of the report introduces ongoing neutron irradiation campaign of yttrium hydride.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

140,142 Ce Neutron Cross Section Resolved Resonance Region Evaluation [Slides]

ORNL has completed new RRR evaluations for 140,142 Ce. A lack of integral measurements made validation difficult, but the difference from Kadonis values suggests that more investigation of capture cross section is warranted. The complete file (with updated File2 and File32) has been submitted to NNDC repository for inclusion in next release of ENDF. A sponsor report publication is in progress.

07 ISOTOPE AND RADIATION SOURCES↗

Status of the RRR evaluations for 233,235 U and 239 Pu [Slides]

The presentation discusses the R -Matrix Analysis of 233 U, 235 U, and 239 Pu and how their interplay among PFNS, TNC, RRR, and neutron multiplicities has a significant impact on the benchmarks. The presentation also discusses the 233 U improved agreement on benchmarks, an additional investigation on the thickness of beryllium reflectors, work on URR (up to 40 keV) and/or extension of RRR (600 eV to 2 keV). Also discussed is how the 235 U progresses on understanding reactivity rates related to depletion calculations, improved agreement with measured (capture) data is quite insensitive to resolve the negative slope, and how the 239 Pu negative slope as a function of the temperature is not understood yet (bias over 2σ lower than measured).

07 ISOTOPE AND RADIATION SOURCES↗

Proposed Updates to Pu Evaluation for ENDF/B-VIII.0β 2 [Slides]

This presentation shows current status of in ENDF/B-VIII.0 and ENDF/B-VIII.0β 1 Pu evaluations. With a focus on the Pu thermal solutions (PST) and sanity check on fast metal assemblies (PMF). Pu ENDFβ 1 file in the low-energy region included updates in RRR (up to 2.5 keV), TNC and PFNS. Discussed are improvements over ENDF/B-VIII.0 in the C/E benchmarks (PST 81 cases) with preliminary ν p and cross section covariance generation for the resolved resonance energy range up to 2.25 keV. Highlighted are proposed updates for ENDFβ 2 with an extension of the RRR evaluation up to 5 keV, inclusion of Mosby’s data including proper resolution function (thanks to Marian J.), and a particular focus on η energy dependence from thermal up to 10 eV with the impact of the capture-to-fission ratio at 0.29 eV. This presentation concludes with an updated ν p and cross section covariance generation for the resolved resonance energy range up to 5 keV. And additional validation to test performance at about 2 eV and in the keV region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Automated Calculation of U-10Mo Fuel and Zr Cladding Thickness

Hand calculation and verification of cladding and fuel thickness in a hot-isostatically pressed uranium-molybdenum (U-10Mo) alloy can incur human errors and longer image processing time when analyzing a large set of cross-sectional images. To help alleviate both time cost and errors made in hand measurements, an automated image processing procedure was developed using Octave, an open-source MATLAB alternative, to repeatedly determine the thickness of Zr and U-10Mo layers. U-10Mo specimens were imaged using standard secondary electron and backscattered-electron imaging at 250× magnification to capture the various layers present in U-10Mo. Further image processing used pixel-by-pixel calculation of the Zr and U-10Mo layers to gather quantitative statistics on the thickness variations associated with each of the layers.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

New Design Yields Robust Large-Area Framing Camera

Manuscript to be submitted to Review of Scientific Instruments to determine interest in possible publication. Website: https://aip.scitation.org/rsi/authors/manuscript High-speed imaging is a core diagnostic capability for the Nevada National Security Site (NNSS). Data from experiments that measure the behavior of materials under extreme conditions are essential for the stockpile mission. For example, during subcritical experiments, materials that are shocked and compressed under high-explosive drive change behavior rapidly over timescales measured in millionths or even billionths of a second. Characterizing and understanding this behavior requires radiographic and/or optical imaging to track motion on these timescales. Framing cameras provide a discrete series of images over a short period of time in a manner that mimics a high-speed movie camera, but with individual shutter times that must be extremely short in order to capture freeze-frame images of rapidly evolving phenomena such as high-explosive shock-driven ejecta, dynamic compression of metals, and high-velocity fluid flow.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Kinetics of Radiation-induced Cr( II ) and Cr( III ) Redox Chemistry in Molten LiCl-KCl Eutectic

Chromium (Cr) is a frequent constituent of the metal alloys proposed for molten salt nuclear reactor (MSR) applications, and is typically the least noble metal ion present. Consequently, chromium is preferentially corroded into molten salt solutions. Here, the redox poise and redox cycling of chromium ions in the salt can greatly influence its corrosivity towards structural alloys, ultimately impacting the longevity of MSR systems. Radiation-induced chemistry is expected to play a significant role in determining the chromium oxidation state distribution during MSR operations. In the present research, electron pulse radiolysis techniques were employed to characterize the reactivity of Cr(II) and Cr(III) ions with primary radiolysis products in molten lithium chloride–potassium chloride (LiCl–KCl) eutectic over a temperature range of 400–600 °C. Both chromium oxidation states were found to rapidly react with the primary products of molten chloride salt radiolysis, i.e., the solvated electron (e S - ) and the dichlorine radical anion (Cl 2 ˙ - ). For reactions with the e S - , second-order rate coefficients (k) of k = (4.1 ± 0.2) and (6.1 ± 0.3) × 10 10 M -1 s -1 at 400 °C for Cr(II) and Cr(III), respectively, were determined. Temperature-dependent measurements allowed for the derivation of activation parameters for electron capture by Cr(II) and Cr(III). Both chromium ions also react with Cl 2 ˙ - , k = (7.2 ± 0.3) and (1.4 ± 0.1) × 10 9 M -1 s -1 at 400 °C for Cr(II) and Cr(III), respectively.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗

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↗

Ab-initio informed cluster dynamics simulation of self- and Xe diffusivity in uranium mononitride under irradiation

Uranium mononitride (UN) is one of the ceramic nuclear fuel alternatives to oxide fuels considered for light water reactor and advanced reactor designs, as it presents significant advantages such as high uranium density (better economics) and high thermal conductivity and melting point (increased safety). Self- and fission gas diffusivities need to be better understood, given that they influence key fuel performance phenomena like swelling and fission gas release. Recently, radiation enhanced diffusivity was investigated in UN by means of cluster dynamics simulations relying on empirical potential-based parameterizations, the reliability of which highly depends on the interatomic potential accuracy. Here, in this work, we refine this approach by determining, using ab-initio calculations, the properties of defect clusters containing vacancies, self-interstitials and Xe impurities. We also consider larger clusters than previous studies. The obtained dataset (formation enthalpies, entropies, and migration barriers) is used to parameterize a cluster dynamics model of mobile clusters, and to calculate the defect cluster concentrations under irradiation. This gives us access to the radiation enhanced self- and fission gas diffusivities. Although the resulting diffusivities are close to the values reported in the literature, we find important qualitative differences in the diffusion mechanisms. Capturing the correct mechanisms is crucial to properly describe the chemistry and fission rate dependence of the model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Step-loaded creep testing of Zircaloy-4 cladding at higher temperatures in the α-phase

A refined understanding of zirconium-based cladding thermomechanical performance during rapid transients is essential for enhancing the safety and operation of light-water reactors. Traditional models for zirconium alloys under accident conditions generally assume that creep dominates fuel cladding performance. Here, these historic models have largely remained unchanged and serve as the basis for safety criteria development. As the U.S. nuclear industry pursues higher burnup levels, the increased release of fission gases during transients raises the risk of cladding failure in the low-temperature hcp α-phase, making the fidelity of these models of greater importance. Creep testing was conducted from 550–700°C with 25–120 MPa applied hoop stresses to investigate Zircaloy-4 deformation at accident-relevant temperatures in the α-phase. Step-loading was employed to capture creep behavior across a wide stress range from a single sample. The stress-strain rate data at higher temperatures (650 and 700°C) were well-described by isotropic versions of the Erbacher and Kaddour models, while the lower temperature data (550 and 600°C) were underpredicted by both anisotropic and isotropic model variants. Greater strain rates during the initial loading step at 650 and 700°C were attributed to recrystallization and grain growth of sub-micron crystallites. Yet, texture analysis revealed the basal split texture remained after testing. These observations produced results suggesting Zircaloy-4 claddings experience higher creep rates across the α-phase than previously thought, possibly related to dynamic anisotropy due to temperature dependent activation of deformation mechanisms, effects of biaxial loading, and variation in material condition between the current testing used in previous model development.

36 MATERIALS SCIENCE↗

Estimating outflow masses and velocities in merger simulations: Impact of r -process heating and neutrino cooling

Here, the determination of the mass, composition, and geometry of matter outflows in black hole-neutron star and neutron star-neutron star binaries is crucial to current efforts to model kilonovae and to understand the role of neutron star merger in r-process nucleosynthesis. In this manuscript, we review the simple criteria currently used in merger simulations to determine whether matter is unbound and what the asymptotic velocity of ejected material will be. We then show that properly accounting for both heating and cooling during r-process nucleosynthesis is important to accurately predict the mass and kinetic energy of the outflows. These processes are also likely to be crucial to predict the fall-back timescale of any bound ejecta. We derive a model for the asymptotic velocity of unbound matter and binding energy of bound matter that accounts for both of these effects and that can easily be implemented in merger simulations. We show, however, that the detailed velocity distribution and geometry of the outflows can currently only be captured by full three-dimensional fluid simulations of the outflows, as nonlocal effect ignored by the simple criteria used in merger simulations cannot be safely neglected when modeling these effects. Finally, we propose the introduction of simple source terms in the fluid equations to approximately account for heating/cooling from r-process nucleosynthesis in future seconds-long three-dimensional simulations of merger remnants, without the explicit inclusion of out-of-nuclear statistical equilibrium reactions in the simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

The role of irradiation-enhanced interstitial diffusion in over-pressurizing fission gas bubbles in UO 2

Fission gas bubbles in UO 2 nuclear fuel have been observed to exhibit pressures in excess of the equilibrium bubble pressure; however, the cause of bubble over-pressurization has not yet been demonstrated. The mechanical interaction between a bubble and the surrounding matrix or grain boundary depends on the internal pressure of the bubble and local stress state, such that over-pressurized bubbles are thought to be responsible for fragmentation and pulverization, when exposed to a temperature ramp. Here, in this work, we investigate the role of U interstitials, produced through irradiation, in over-pressurizing bubbles by using a combined molecular dynamics (MD) and cluster dynamics approach. Firstly, the energies for the capture of interstitials and vacancies by bubbles have been determined from MD as a function of the ratio of gas atoms to vacancies that make up the bubble. Secondly, these reaction energies have been implemented in the cluster dynamics code Centipede to predict bubble over-pressurization as a function of temperature for typical fission rates. It was found that there is a transition from low pressure bubbles (at high temperatures) to high pressure bubbles (at lower temperatures). The cause of this behavior was shown to be the creation of irradiation-induced interstitials that are highly mobile relative to vacancies at low temperature; whereas, vacancies are sufficiently mobile at high temperatures to limit bubble pressures. This result supports the hypothesis that over-pressurized bubbles form during steady-state operation and that this behavior is highly sensitive to the local pellet temperature.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Evaluation of Hybrid FPOG Applications in Regulated and Deregulated Markets Using HERON

Recent changes in the U.S. energy market, such as low natural gas prices and increased electricity production for variable renewable energy (VRE) sources, have led to an economic crisis for existing light-water reactor (LWR) nuclear power plants (NPP). Many owners and operators of LWRs have elected to decommission these plants rather than continue using them as consistent sources of clean baseload power. This has led to exploration of various possibilities to increase the economic viability of these units, including market restructuring to monetize benefits LWRs already provide to the grid through ancillary markets, load following and economic dispatch, and possible integration of secondary systems directly to the NPP for production of additional products through technologies such as hydrogen electrolysis or water desalination. Previous studies have considered the technologies associated with these Integrated Energy Systems (IES) activities, and the analysis of markets for these secondary products. To analyze the economic viability of various system configurations including IES, especially given the uncertainty surrounding load demand, electricity prices, and the availability of VRE resources, the stochastic technoeconomic analysis package HERON (Heuristic Energy Resource Optimization Network) was released earlier this year as an extension of the risk analysis framework RAVEN (Risk Analysis Virtual Environment). HERON focuses foremost on making the complex uncertainty quantification analysis tools approachable for energy systems analysts, also providing general dispatch optimization algorithms for those workflows. HERON continues to be improved and tested as a significant part of the IES viability analyses performed in this work. HERON is not a capacity expansion model. To consider market and grid energy system development in a variety of scenarios, HERON is best used in coupling with modelling tools such as US-REGEN, which sacrifice some of the uncertainty analysis and resolution of HERON's modelling for the ability to efficiently predict the change in the grid energy system's profile due to economic drivers over decades. HERON can then use this information to explore the economic viability of introducing changes to the predicted outcomes, such as the introduction of an IES. In this work, experts at EPRI using US-REGEN provide six projection scenarios for use in HERON stochastic technoeconomic analysis (STEA) in considering the options available for increasing LWR economic viability through introduction of a hydrogen-centric IES using a high-temperature steam electrolysis plant (HTSE), hydrogen storage, and a constant-rate contracted hydrogen consumer. The results obtained are differential in nature; they do not report expected profits for any configuration, but rather report on the possible increase in the NPV of a configuration with respect to a baseline no-IES configuration. Due to the uncertainty captured in the variable net load of the systems, there is likewise uncertainty in the mean values reported. We consider this viability both in terms of a regulated market, where the energy producers and IES are owned and operated by single entity, as well as a deregulated market, where the IES chooses its bid for electricity generation and is then dispatched by the grid system operator. Results indicate that for deregulated markets, the inclusion of the IES is often statistically beneficial. This is especially true in policies that are not favorable towards nuclear, as nuclear is less often dispatched and is forced to deal with frequent idle capacity. In the nominal case as well as the case of carbon tax policies, inclusion of the IES clearly benefited the economic performance of the NPP. In the regulated case, however, there was a trend towards minimizing the IES, likely due to the optimal sizing performed by US-REGEN of the NPP within the system as well as the lack of penalty for idle capacity at the NPP in the regulated market analyses.

99 GENERAL AND MISCELLANEOUS↗

An ab initio exciton model for singlet fission

Here we present an ab initio exciton model that extends the Frenkel exciton model and includes valence, charge-transfer, and multiexcitonic excited states. It serves as a general, parameter-free, yet computationally efficient and scalable approach for simulation of singlet fission processes in multichromophoric systems. A comparison with multiconfigurational methods confirms that our exciton model predicts consistent energies and couplings for the pentacene dimer and captures the correct physics. Calculations of larger pentacene clusters demonstrate the computational scalability of the exciton model and suggest that the mixing between local and charge-transfer excitations narrows the gap between singlet and multiexcitonic states. Local vibrations of pentacene molecules are found to facilitate singlet–multiexcitonic state-crossing and hence are important for understanding singlet fission. The exciton model developed in this work also sets the stage for further implementation of the nuclear gradients and nonadiabatic couplings needed for first principles nonadiabatic quantum molecular dynamics simulations of singlet fission.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗