Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fuel lattices”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Quantization effects in semiconductor nanostructures and singlet fission in molecular chromophores for photovoltaics and solar fuels

Hot carriers in semiconductors are electrons and/or holes that have energies greater than carriers that reside at the top and bottom of the conduction and valence bands, respectively; the latter carriers are in equilibrium with the lattice and have a temperature equal to the lattice (ambient) temperature. Hot carriers are created in semiconductors upon the absorption of photons with energies greater than the bandgap. The excess energy above the bandgap energy is in the form of kinetic energy. The hot carriers can cool to the lattice temperature via electron–phonon scattering and establish separate Boltzmann distributions for electrons and holes at the band edges. Thus, upon cooling, the excess kinetic energy of the hot carriers is transformed into heat and is unavailable to be converted into electrical or chemical free energy in a solar photon-converting cell. This hot-carrier cooling process significantly limits the maximum possible power conversion efficiency (PCE) of the free energy of solar photons into electrical free energy or chemical free energy in chemical bonds of fuels. However, if hot-carrier cooling can be slowed such that the hot carriers can be extracted before cooling to create higher photovoltages, or utilized to create additional electron-hole pairs through carrier multiplication, then beneficial use can be made of the excess energy of hot carriers and the maximum thermodynamic PCE can be dramatically increased. Quantization effects in semiconductor nanostructures, unique properties of some bulk semiconductors, and exciton multiplication in both semiconductors and molecular chromophores that undergo singlet fission have shown a path forward for potential dramatic increases in the PCEs of solar photons into electricity and solar fuels. The status, history, and future promise of the science and technology of these future approaches for different types of photovoltaic cells and cells for solar fuels will be discussed.

14 SOLAR ENERGY↗
DOE OSTI2025

A machine learning approach to quantify degradation of nuclear fuels and the effects of fission products

Nuclear fuel performance is critically dependent on understanding the evolution of fuel properties under operational conditions, a complex challenge driven by chemical changes and substantial radiation damage during fission. Traditionally, property evolution has been determined via empirical data collected following irradiation. However, these empirical correlations are limited in their applicability beyond the specific conditions in which they were obtained. This study explores a novel approach to address this challenge by applying materials informatics to develop a machine learning random forest (ML-RF) model that captures the effects of fission products on fuel compounds. The model predicts formation enthalpy (ΔH f ) by leveraging extensive quantum materials property data and correlating it with material descriptors such as composition, atomic and site features, and crystal lattice properties. This ML-RF model enables rapid interpolation across the compositional and structural spaces covered by the training data, thus supporting high-throughput screening and energetic ranking of candidate phases. The model demonstrates the ability to predict ΔH f with a mean absolute error (MAE) of approximately 0.1 to 0.2 eV/atom across a wide range of compounds, including key nuclear fuel systems (U-O, U-N, U-C, U-Si, and U-Mo). For example, it was used to assess shifts in stoichiometry for UO 2 (O/M) and UN (N/M) fuels, revealing their distinct tendencies in chemical potential variation and enabling preliminary convex hull analyses. Furthermore, the model provides insights into how individual fission products affect fuel properties. Results indicate that larger fission products (e.g., Nd, Pu, Ce) have a more pronounced impact on UO 2 , while lighter ones (e.g., Zr) strongly influence UN. Here, the model developed in this work can be used to support the Accelerated Fuel Qualification approach by facilitating preliminary evaluations prior to extensive materials modeling and experimentation. To this end, the trained model has been made available to the fuel community to support ongoing fuel development efforts.

Accelerated fuel qualification↗
DOE OSTI2025

A strategy for automated core design to increase economic viability and minimize fuel fragmentation, relocation, and dispersal susceptibility in high-burnup cores

The nuclear industry aims to increase the cycle length of pressurized water reactors from 18 to 24 months to increase power plant capacity factors and economic viability. These cycle length extensions will inherently require fuel rods to exceed the current peak rod average burnup limit of 62 GWd/MTU. A chief concern of operating beyond the current burnup limit is the fuel fragmentation, relocation, and dispersal (FFRD) phenomenon in which pulverized fuel fragments can axially relocate and escape through a burst in the cladding formed during a loss-of-coolant accident. In this work, we demonstrate an approach for automating core design employing an optimization tool based on a penalty-free, parallel simulated annealing algorithm to produce pressurized water reactor core designs with two different optimization objectives. The two objectives were to produce core designs with (1) mitigated FFRD susceptibility while achieving 24-month cycle lengths (2) maximum cycle length with no regard for the likelihood of FFRD. Batch size was considered in tandem with both cases to maximize economic viability. The PARCS nodal model was the primary reactor physics tool used in the optimizations and used nuclear cross sections calculated with 2D Polaris lattice physics models. Reactor performance and safety characteristics of the optimized cores were verified using high-fidelity Virtual Environment for Reactor Applications models. The core designs produced by the optimization tool are compared with each other and to a high-burnup core design produced and analyzed in previous works to highlight the fuel management strategies that may enhance high-burnup reactor safety and economic viability. The optimized cores satisfied their respective objective functions, producing a maximum cycle length of 720 effective full-power days in one core design and one that may reduce FFRD susceptibility by up to 50% based on the first-order approximation to FFRD risk formulated in this work. The optimized cores met most constraints but exceeded the hot channel factor limit, especially in FFRD cases where fresh fuel carried more power. Furthermore, this highlights the need for future lattice-level optimizations and broader assembly options.

Cycle length↗
DOE OSTI2022

Studsvik R2 materials test reactor Ad Hoc depletion strategy for the derivation of the fuel isotopic composition of the MPCMIV benchmark

The Ad Hoc Depletion Strategy elaborated by the NINE company, developed in support of the organization of the MPCMIV (Multi-physics Pellet Cladding Mechanical Interaction Validation) benchmark input and output specifications, is presented. This work aims at illustrating the strategy itself and then showing the results obtained with its application over the Studsvik R2 Testing Reactor, which is analyzed in the benchmark. The objective of the application of the strategy is to compute the fuel elements isotopic compositions at the beginning of some core loadings of interest for the benchmark. To this objective, it is necessary to implement first the simulation model of the three single assembly types and perform the infinite lattice depletions, then, to build the full core model and to perform the simulation of the core cycle. All the models and simulations were carried out with the use of the Monte Carlo particle transport code Serpent 2. Finally, the simulations results are assessed against Studsvik isotopic compositions of the fuel elements discharged from the R2 Reactor at the end of the core loading. Several assumptions were necessary during all the steps of the strategy, to overcome the lack of information regarding the core management. For this reason, the solution found at the end of the current analysis may not be completely optimized and further improvements regarding the model assumptions will be tested in a future work. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

MPACT 4.4 Theory Manual

MPACT is a three-dimensional (3D) full-core neutron transport code capable of calculating subpin power distributions. Calculations are based on the Boltzmann transport equation for neutron fluxes for problems in which the detailed geometrical configuration of fuel components such as the pellet and cladding are explicitly retained. The cross-section data needed for the neutron transport calculation are obtained directly from a multigroup cross section library, which has traditionally been used by lattice physics codes to generate few-group homogenized cross sections for nodal core simulators. Hence, MPACT assumes neither a priori homogenization nor group condensation for the full core spatial solution. The 3D MPACT transport solution can be obtained using the method of characteristics (MOC), which employs discrete ray tracing within each fuel pin. However, for practical reactor applications, the direct application of MOC to 3D core configurations requires an excessive amount of memory and computing time due to the very large number of rays. For practical 3D full-core calculations, MPACT commonly uses an approximate “2D/1D” method that treats the radial (x and y) variables differently from the axial (z) variable. In particular, the radial dependence of the solution is calculated using transport theory, and the axial dependence is calculated using diffusion or P 3 theory. The 2D/1D method requires the core to be divided into a vertical stack of axial slices with a thickness of Δ z ≈ 5–10 cm. Each axial slice is divided radially into coarse spatial cells with boundaries that usually constitute the pin cell boundaries, for which Δ x = Δ y ≈ 1.5 cm. Then, each coarse radial cell (pin cell) is divided into 50–100 fine radial cells, which resolve the angular flux in the fuel, cladding, and moderator regions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Deep Learning for Multigroup Cross-Section Representation in Two-Step Core Calculations

Here we investigate using deep learning, a type of machine-learning algorithm employing multiple layers of artificial neurons, for the mathematical representation of multigroup cross sections for use in the Griffin reactor multiphysics code for two-step deterministic neutronics calculations. A three-dimensional fuel element typical of a high-temperature gas reactor as well as a two-dimensional sodium-cooled fast reactor lattice are modeled using the Serpent Monte Carlo code, and multigroup macroscopic cross sections are generated for various state parameters to produce a training data set and a separate validation data set. A fully connected, feedforward neural network is trained using the open-source PyTorch machine-learning framework, and its accuracy is compared against the standard piecewise linear interpolation model. Additionally, we provide in this work a generic technique for propagating the cross-section model errors up to the k eff using sensitivity coefficients with the first-order uncertainty propagation rule. Quantifying the eigenvalue error due to the cross-section regression errors is especially practical for appropriately selecting the mathematical representation of the cross sections. We demonstrate that the artificial neural network model produces lower errors and therefore enables better accuracy relative to the piecewise linear model when the cross sections exhibit nonlinear dependencies; especially when a coarse grid is employed, where the errors can be halved by the artificial neural network. However, for linearly dependent multigroup cross sections as found for the sodium-cooled fast reactor case, a simpler linear regression outperforms deeper networks.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Preliminary modeling of triply periodic minimal surface (TPMS) structures using RELAP5-3D

With the United States Department of Energy (DOE)’s goal of quadrupling the nation’s nuclear energy supply by 2050, and with the Advanced Fuels Campaign pushing for new types of advanced reactor fuels and geometries, the need has arisen for new nuclear fuel designs. One such design is to swap out current nuclear fuel geometries in exchange for another type of geometry, called a Triply Periodic Minimal Surface (TPMS). TPMSs are self-supporting, infinitely repeating lattices—attributes that lend themselves well to additive manufacturing. These surfaces also possess enhanced heat transfer properties thanks to their internal area changes and large surface-area-to-volume ratios. Their drawback, however, is an increased pressure drop. Given the small amount of correlations and data (Reynolds numbers in the 2,000–8,000 range), and the minimal amount of experience so far obtained by modeling TPMS structures using 1D systems codes such as the Reactor Excursion and Leak Analysis Program (RELAP5-3D), further research into this topic was needed. Using data from the University of Wisconsin - Madison (UW), curve fits were created for both a Heat Transfer Coefficient (HTC) correlation and a Darcy friction factor empirical coefficient correlation. The curves’ coefficients and multipliers were then output and utilized in RELAP5-3D models of two upcoming experiments—Flow Loop for INFLUX Pressure drop (FLIP) and Microreactor Agile Non-nuclear Experimental Test (MAGNET)—aimed at increasing the available data for Reynolds numbers to the 16,000–36,000 range for TPMS structures. The models were run under the conditions utilized by a Computational Fluid Dynamics (CFD) analysis performed by another group at Idaho National Laboratory. Only CFD pressure drop values were obtained from the FLIP test, and those values showed that the RELAP5-3D models had a lower rate of pressure increase in comparison to the CFD values. In addition, there seemed to be a vertical shift upward in the pressure drop for both models whenever the TPMS porosity decreased, and the RELAP5-3D models showed a higher vertical shift in comparison to the CFD values. The MAGNET results did not correspond to any CFD or experimental results against which they could be compared, so they were instead compared against the proposed CFD input conditions. These values were then compared with each other to make sure the model seemed to be performing as expected, paving the way for future tests that can be run for the purpose of further analyses and comparisons. The pressure drop increased with temperature and mass flow rate independently. The temperature change would decrease with increasing mass flow rate and temperature, which was just as we expected based on the fact that the lower viscosity and decreased density would result in higher friction and churning losses. The last metric that was assessed was the enthalpy flow change, which increased with increasing mass flow rate and decreasing temperature.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Influence of Pt-Metal Alloy Catalysts with Various Ionomers on Oxygen Reduction Reaction in Fuel Cell Application

Pt-M/C (M = Co, Ni, Mn, etc.) alloy catalysts exhibit superior oxygen reduction reaction (ORR) activity compared to pure Pt/C, leading to a high energy efficiency in hydrogen fuel cells. However, many Pt-M/C alloy catalysts were synthesized and evaluated at the lab scale in model test-bed systems like rotating disc electrodes, which don't always correlate to performance within a fuel cell system; there is a clear need to evaluate catalysts in electrodes that can be prepared at industrially relevant scales to evaluate how factors like ink formulation can greatly affect device-level of fuel cell performance. Herein, three commercial Pt-M/C alloy catalysts (two Pt-Co/C and one Pt-Ni/C) were comprehensively characterized by various techniques. The results show that the average particle sizes of the three catalysts are close to 5 nm; the atomic ratio of Pt/M is around 4; and the M was successfully embedded into Pt lattice, resulting in the positive shift of Pt 4f in XPS spectra and XRD patterns. These catalytic materials were incorporated into 9 different cathode catalyst layers (CCLs) with three kinds of ionomers (Nafion D2020, high oxygen permeability ionomer (HOPI), and Aquivion D79-25BS), and their performance in proton exchange membrane fuel cells (PEMFCs) were investigated. The results demonstrate that the Pt-Co/C catalysts possess a higher mass activity (MA) than Pt-Ni/C; the cathodes with Nafion ionomer provide the highest MA while electrodes with Aquivion ionomer showed the lowest activity, attributed to poor H+ conductivity resulting from suboptimal ionomer incorporation. Finally, these alloys were shown to exceed DOE targets for MA and H2/Air performance reported in the recent publications at beginning of life and after 90k cycle catalyst AST protocol. This study provides valuable performance benchmarks for these materials guiding future Pt-M/C catalyst design and material integration for heavy duty PEMFC applications.

08 HYDROGEN↗
DOE OSTI2024

Generation of Enrichment-Dependent Thermal Neutron Scattering Data

This work details the generation of enrichment-dependent thermal neutron scattering cross sections for several crucial uranium fuel compounds. The evaluations of the thermal scattering law (TSL) and associated cross sections for uranium dioxide (UO 2 ), uranium carbide (UC), and uranium nitride (UN) were performed using standard ab initio lattice dynamics (AILD) methods. The data for uranium metal was produced using a novel hybrid approach of molecular dynamics combined with lattice dynamics methods. 235 U enrichments of 5%, 10% (LEU+), 19.75% (HALEU), 93% (HEU), and 100% were considered, in addition to natural uranium. The enrichment-dependent masses and free atom cross sections were used in the generation of elastic and inelastic thermal neutron scattering cross sections, while the calculation of the phonon density of states (DOS) and resulting TSL considered only the natural isotopic composition of uranium. The use of an identical DOS for all enrichments is expected to have minimal impact on the final data, as the small change in uranium mass should not significantly affect lattice vibrations. The cross sections are shown to exhibit significant dependence on 235 U enrichment. The submission of this data to the National Nuclear Data Center (NNDC) for release in the ENDF/B-VIII.1 database should support the design of advanced reactor concepts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

White Paper: Research & Development for the Time at Temperature Approach

Recent advancements in nuclear power research are greatly improving reactor safety and performance through the development of Accident Tolerant Fuel (ATF) and Low-Enriched Uranium Plus (LEU+). These innovations can address Departure from Nucleate Boiling (DNB) margins, which are vital for reactor safety. DNB happens when the coolant switches to film boiling, significantly decreasing heat transfer and posing a risk of fuel cladding failure. The U.S. Nuclear Regulatory Commission (NRC) employs conservative DNB criteria, which can potentially restrict the operational flexibility and efficiency of reactors. The Time at Temperature (TaT) approach could provide a more detailed and adaptable operational guideline by establishing acceptable time-temperature limits, accounting for the duration a material can withstand elevated temperatures without losing its integrity. This method allows reactors to operate more efficiently and safely, offering additional operational margins, faster power adjustments, and improved fuel cycle economics. TaT criteria allow for higher power levels and more flexible responses to operational transients, particularly applicable for anticipated operational occurrences (AOOs) that result in short durations of post-DNB conditions. It enhances plant operational flexibility, allows faster startup times, and enables quicker power level adjustments, optimizing fuel loading patterns and improving fuel cycle economics. Implementing TaT limits reduces core design constraints, lowers fuel usage, and reduces costs, essential for the long-term sustainability of Light Water Reactors (LWRs). TaT maximizes the use of advanced fuel technologies like ATF and LEU+, further enhancing their economic and environmental benefits. To apply the TaT approach in existing LWRs, collaborative research activities among various DOE-sponsored programs are essential. These efforts should incorporate fuel experiments, physics-based high-fidelity modeling, ML-based surrogate modeling, and optimization techniques. This whitepaper proposes four research and development areas: 1) Investigation of the feasibility of new operations of LWR with updated safety limits; 2) Assessment of reactor operation limits through uncertainty reduction; 3) Evaluation of power uprate in virtual environment; and 4) Lattice and reactor core design for power uprate. Each area includes why this research is in need and a suggested scope of work. These comprehensive research areas ensure practical and beneficial advancements for existing reactors, translating innovations in nuclear fuel and cladding technology into improved reactor performance and safety.

42 - ENGINEERING↗
DOE OSTI2023

Temperature-dependent elastic constants of thorium dioxide probed using time-domain Brillouin scattering

Herein we report the adiabatic elastic constants of single-crystal thorium dioxide over a temperature range of 77–350 K. Time-domain Brillouin scattering, an all-optical, non-contact picosecond ultrasonic technique, is used to generate and detect coherent acoustic phonons that propagate in the bulk perpendicular to the surface of the crystal. These coherent acoustic lattice vibrations have been monitored in two hydrothermally grown single-crystal thorium dioxide samples along the (100) and (311) crystallographic directions. The three independent elastic constants of the cubic crystal (C 11 , C 12 , and C 44 ) are determined from the measured bulk acoustic velocities. The longitudinal wave along the (100) orientation provided a direct measurement of C 11 . Measurement of C 44 and C 12 was achieved by enhancing the intensity of quasi-shear mode in a (311) oriented crystal by adjusting the polarization angle relative to the crystal axes. We find the magnitude of softening of the three elastic constants to be ~2.5% over the measured temperature range. Good agreement is found between the measured elastic constants with previously reported values at room temperature, and between the measured temperature-dependent bulk modulus with calculated values. We find that semi-empirical models capturing lattice anharmonicity adequately reproduce the observed trend. We also determine the acoustic Grüneisen anharmonicity parameter from the experimentally derived temperature-dependent bulk modulus and previously reported temperature-dependent values of volumetric thermal expansion coefficient and heat capacity. This work presents measurements of the temperature-dependent elasticity in single-crystal thorium dioxide at cryogenic temperature and provides a basis for testing ab initio theoretical models and evaluating the impact of anharmonicity on thermophysical properties.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2000

Fretting Stresses in Single Crystal Superalloy Turbine Blade Attachments

Single crystal nickel base superalloy turbine blades are being utilized in rocket engine turbopumps and turbine engines because of their superior creep, stress rupture, melt resistance and thermomechanical fatigue capabilities over polycrystalline alloys. Currently the most widely used single crystal nickel base turbine blade superalloys are PWA 1480/1493 and PWA 1484. These alloys play an important role in commercial, military and space propulsion systems. High Cycle Fatigue (HCF) induced failures in aircraft gas turbine and rocket engine turbopump blades is a pervasive problem. Blade attachment regions are prone to fretting fatigue failures. Single crystal nickel base superalloy turbine blades are especially prone to fretting damage because the subsurface shear stresses induced by fretting action at the attachment regions can result in crystallographic initiation and crack growth along octahedral planes. Furthermore, crystallographic crack growth on octahedral planes under fretting induced mixed mode loading can be an order of magnitude faster than under pure mode I loading. This paper presents contact stress evaluation in the attachment region for single crystal turbine blades used in the NASA alternate Advanced High Pressure Fuel Turbo Pump (HPFTP/AT) for the Space Shuttle Main Engine (SSME). Single crystal materials have highly orthotropic properties making the position of the crystal lattice relative to the part geometry a significant factor in the overall analysis. Blades and the attachment region are modeled using a large-scale 3D finite element (FE) model capable of accounting for contact friction, material orthotrophy, and variation in primary and secondary crystal orientation. Contact stress analysis in the blade attachment regions is presented as a function of coefficient of friction and primary and secondary crystal orientation, Stress results are used to discuss fretting fatigue failure analysis of SSME blades. Attachment stresses are seen to reach peak values at locations where fretting cracks have been observed. Fretting stresses at the attachment region are seen to vary significantly as a function of crystal orientation. Attempts to adapt techniques used for estimating fatigue life in the airfoil region, for life calculations in the attachment region, are presented. An effective model for predicting crystallographic crack initiation under mixed mode loading is required for life prediction under fretting action.

Arakere, Nagaraj K.↗
DOE OSTI2021

Laser-induced annealing of aged PuO 2

Plutonium dioxide (PuO 2 ) is an important compound used in nuclear fuel, irradiation targets, and heat sources. As such, improved understanding of its structural and spectroscopic properties has numerous applications. Alpha particle-induced damage of a PuO 2 crystal lattice modifies several properties of its Raman spectrum, including band intensities, positions, and widths. The decay also induces growth of new bands and creates electron-trapped defects with luminescent properties. Herein, we show for the first time that laser-induced heating can reverse damage to an aged and damaged PuO 2 lattice. Using automated instrumentation to heat a single 6–10 μm spot to temperatures above 1,300°C, we show that laser-induced annealing of aged PuO 2 results in restoration of T2g band intensity with a concomitant decrease in defect band intensity, a result that indicates laser annealing can be used to reverse age-damage in PuO 2 . This laser annealing approach permitted in situ observation of temperature-dependent Raman spectral changes, thereby providing insight into the thermodynamics of structural alterations in a radiolytically damaged PuO 2 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Polaris-PARCS Sensitivity Study on LWR Fuel Cycles: Polaris Input Options

This study is the first of a multi-phase effort to assess the sensitivity of light-water reactor (LWR) core-level prediction biases to changes in lattice-level calculation parameters. Prediction bias is the measured-to-predicted difference in a core-level quantity of interest (QOI) which can be estimated by comparing the simulation results with the plant-measured data for key nuclear parameters. The LWR two-step neutronics codes employed herein are the SCALE–Polaris lattice physics code (v6.3.1) and the Purdue Advanced Reactor Core Simulator (PARCS) nodal diffusion simulator (v3.4.2), both funded and used for confirmatory analysis to support licensing by the US Nuclear Regulatory Commission (NRC). Polaris–PARCS is used to model Watts Bar Unit 1 cycles 1–3 and Peach Bottom Unit 2 cycles 1–3. This study focuses on the impact of changes to Polaris input options such as scattering treatment or quadrature settings and how these input options induce changes in core-level quantities of interest (QOIs)bias. The report documents multiple bias assessments for different modeling choices and compares the bias magnitude to the QOI measurement uncertainties. Future companion reports will investigate the sensitivity of core-level LWR prediction bias to Polaris input options and Polaris-computed QOIs such as few-group assembly-homogenized cross sections to gain an understanding of the key drivers of prediction bias at lattice and core levels for application of a two-step LWR neutronics procedure in a licensing scenario.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Preliminary Study on TRISO Fuel Cross Section Generation

Cross section self-shielding methodologies for TRISO fuel were assessed to provide accurate multigroup cross sections for a high-fidelity reactor physics code so that the code is able to accurately model and simulate advanced reactors with TRISO fuel. Initially, the two existing methodologies (the SCALE method and the Sanchez-Pomraning method) were studied and implemented to MC2-3 for detailed performance tests. Additionally, a new spatial self-shielding method, named the iterative local spatial self-shielding (ILSS) method, for particulate fuels was developed based on the disadvantage factor and implemented to MC2-3 as well. The new method approximately accounts for the effect of randomly distributed particles on the particle shadowing effect using a homogenized compact region surrounding a particle of interest at the center. The self-shielded cross sections of the particle at the center are determined iteratively since the cross sections of the homogenized compact region are calculated using them. For the energy range above 100 keV where the fuel-to-moderator ratio is more important than the random distribution of particles, a single particle unit-cell model is used by preserving the average amount of moderator per fuel particle in the system. The three self-shielding methods implemented in MC2-3 were tested using numerical benchmark problems made based on fuel compact problems of a prismatic-type very high temperature reactor. Test results indicated that the ILSS method produced slightly better results than the SCALE and Sanchez-Pomraning methods, compared to the Serpent-2 Monte Carlo results obtained with 25 independent random particle configurations. The SCALE and Sanchez-Pomraning methods tend to underestimate the heterogeneity effect by 150 and 100 pcm, respectively, while the new ILSS method overestimates the heterogeneity effect by 70 pcm. In future, the new self-shielding method will be extended to perform pebble calculations and compare results with those from the SCALE and Sanchez-Pomraning methods. Furthermore, the new method will be optimized for practical applications to on-the-fly resonance treatment for lattice or whole-core calculations for advanced reactors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRSDate not supplied

Molten Salt Lattice Confinement Fusion (LCF) Fast Fission Reactor for Lunar and Planetary Surface Power

Molten salt fission reactors (MSR) have been suggested for lunar and planetary surface power systems. They have the advantage of operating at high temperature, for efficient thermal-electrical conversion, low pressure, long-lived with high nuclear fuel burnup. MSR are often designed to breed fissile 233U from natural 232Th by neutron capture and  decay via: 232Th(n,)233Th(,)233Pa(,)233U Unfortunately, this process requires 233Pa isotope separation and segregation to decay to 233U. This requirement prevents additional neutron capture that interferes with 233U breeding. Instead, Wooley’s sub-critical, fast fission, molten salt reactor would use externally generated tokamak fusion neutrons1,2 to fission all actinides.We propose a simpler fusion-fast-fission sub-critical reactor that generates fast neutrons in situ from lattice confinement fusion (LCF) to fission fertile and fissile actinides. This hybrid reactor doesn't require enriched 235U fissile pins to initiate fis-sion reactions, nor 233Pa separation and segregation during operation. Like Wooley’s, this hybrid reactor “burns” natural uranium (238U) or thorium (232Th) which avoids uranium enrichment and additional fissile material launch safety and security costs. The LCF neutron source is initiated by bremsstrahlung photoneutrons (Fig. 1)3 or isotopic neutron sources in electron-screened lattices (Fig. 2)4,5. Alternatively, the electrolytic Pd-deuterium co-deposition6protocol fast fissions7 both 232Th and 238U. However, an aqueous electrolyte-based system, without pressurization similar to conventional pressurized water fission reactors, is incapable of high temperatures due to the boiling point of the electrolyte slightly over 100 C. Molten salts can be used instead as was demonstrated at the University of Hawaii8 using a variety of Ni and Pd cathodes in lithiated, hydrided and deuterated salts. These salts have melting points often exceeding 500C making them suitable to efficiently produce electrical power9 through either Advanced Stirling Genera-tors (< 100 kWe) or closed-Brayton Cycle (> 100 kWe). This hybrid reactor could power a wide range of lunar or Martian applications from unmanned in-struments, to charging vehicles and entire facilities such as human habitats or in situ resource utilization. The power conversion cycles are Carnot Cycle limited, but generally 30% efficient at best. However, waste heat on the moon or Mars is important to surviving either two-week lunar nights or Martian nights as well as providing process heat for mineral extraction and “living off the land”

Lawrence Forsley↗