SEARCH · Search NASA

Results for “Nuclides”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Coproduction of Light and Heavy r-process Elements via Fission Deposition

We apply for the first time fission yields determined across the chart of nuclides from the macroscopic-microscopic theory of the Finite Range Liquid Drop Model to simulations of rapid neutron capture (r-process) nucleosynthesis. With the fission rates and yields derived within the same theoretical framework utilized for other relevant nuclear data, our results represent an important step toward self-consistent applications of macroscopic-microscopic models in r-process calculations. The yields from this model are wide for nuclei with extreme neutron excess. We show that these wide distributions of neutron-rich nuclei, and particularly the asymmetric yields for key species that fission at late times in the r process, can contribute significantly to the abundances of the lighter heavy elements, specifically the light precious metals palladium and silver. Since these asymmetric yields correspondingly also deposit into the lanthanide region, we consider the possible evidence for coproduction by comparing our nucleosynthesis results directly with the trends in the elemental ratios of metal-poor stars rich in r-process material. We show that for r-process enhanced stars palladium over europium and silver over europium display mostly flat trends suggestive of coproduction and compare to the lanthanum over europium trend which is often used to justify robustness arguments in the lanthanide region. We find that such robustness arguments may be extendable down to palladium and heavier and demonstrate that fission deposition is a mechanism by which such a universality or robustness can be achieved.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

New Thermonuclear Rate of 7 Li(d,n)2 4 He Relevant to the Cosmological Lithium Problem

Accurate 7 Li(d,n)2 4 He thermonuclear reaction rates are crucial for precise prediction of the primordial abundances of lithium and beryllium and to probe the mysteries beyond fundamental physics and the standard cosmological model. However, uncertainties still exist in current reaction rates of 7 Li(d,n)2 4 He widely used in big bang nucleosynthesis (BBN) simulations. In this work, we reevaluate the 7 Li(d,n)2 4 He reaction rate using the latest data on the three near-threshold 9 Be excited states from experimental measurements. We present for the first time uncertainties that are directly constrained by experiments. Additionally, we take into account for the first time the contribution from the subthreshold resonance at 16.671 MeV of 9 Be. We obtain a 7 Li(d,n)2 4 He rate that is overall smaller than the previous estimation by about a factor of 60 at the typical temperature of the onset of primordial nucleosynthesis. We implemented our new rate in BBN calculations, and we show that the new rates have a very limited impact on the final light element abundances in uniform density models. Typical abundance variations are in the order of 0.002%. For nonuniform density BBN models, the predicted 7 Li production can be increased by 10% and the primordial production of light nuclides with mass number A > 7 can be increased by about 40%. Our results confirm that the cosmological lithium problem remains a long-standing unresolved puzzle from the standpoint of nuclear physics.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Bayesian Estimation of the D(p,γ) 3 He Thermonuclear Reaction Rate

Big bang nucleosynthesis (BBN) is the standard model theory for the production of the light nuclides during the early stages of the universe, taking place for a period of about 20 minutes after the big bang. Deuterium production, in particular, is highly sensitive to the primordial baryon density and the number of neutrino species, and its abundance serves as a sensitive test for the conditions in the early universe. The comparison of observed deuterium abundances with predicted ones requires reliable knowledge of the relevant thermonuclear reaction rates, and their corresponding uncertainties. Recent observations reported the primordial deuterium abundance with percent accuracy, but some theoretical predictions based on BBN are at tension with the measured values because of uncertainties in the cross section of the deuterium-burning reactions. In this work, we analyze the S-factor of the D(p,γ) 3 He reaction using a hierarchical Bayesian model. We take into account the results of eleven experiments, spanning the period of 1955--2021; more than any other study. We also present results for two different fitting functions, a two-parameter function based on microscopic nuclear theory and a four-parameter polynomial. Furthermore, our recommended reaction rates have a 2.2\% uncertainty at 0.8~GK, which is the temperature most important for deuterium BBN. Differences between our rates and previous results are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Estimation of the NiCu Cycle Strength and Its Impact on Type I X-Ray Bursts

Abstract Type I X-ray bursts (XRBs) are powered by thermonuclear burning on proton-rich unstable nuclides. The construction of burst models with accurate knowledge of nuclear physics is required to properly interpret burst observations. Numerous studies that have investigated the sensitivities of burst models to nuclear inputs have commonly extracted the strength of the NiCu cycle in the rp process, determined by the 59 Cu( p , α ) 56 Ni and 59 Cu( p , γ ) 60 Zn thermonuclear reaction rates, as critical in the determination of reaction flow in the burst. In this study, the strength of the cycle at the XRB temperature range was estimated based on published experimental data. The nuclear properties of the compound nucleus 60 Zn were evaluated for the 59 Cu( p , α ) 56 Ni and 59 Cu( p , γ ) 60 Zn reaction rate calculations. Monte Carlo rate calculations were conducted to include the large uncertainties of nuclear properties in the calculations. In the current work, a weak NiCu cycle is expected, whereas the rates adopted by the previous studies suggest a strong NiCu cycle. Model simulations were performed with the new rates to assess the impact on Type I XRBs. The results show that the estimated cycle strength does not strongly influence the model predictions of the burst light curve or synthesized abundances.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

The 2025 Evaluation of Experimental Thermonuclear Reaction Rates (ETR25)

This work describes the formalism for estimating thermonuclear reaction rates for astrophysical applications, emphasizing modern statistical approaches such as Monte Carlo sampling and Bayesian models. We discuss related topics including the calculation of resonance energies from nuclear Q values, indirect estimates of particle partial widths, and matching of reaction rates at elevated temperatures to statistical model results. We have evaluated available experimental data on cross sections, resonance energies and strengths, partial widths, lifetimes, spin-parities, and spectroscopic factors. Based on these results, we have estimated numerical values of 78 experimental charged-particle thermonuclear reaction rates for target nuclei in the A = 2–40 mass region, for temperatures ranging from 1 MK to 10 GK. For each reaction, three rate values are provided: low, median, and high, corresponding to the 16th, 50th, and 84th percentiles, respectively, of the cumulative reaction rate probability density distribution. Additionally, we present the factor uncertainty of each rate at each temperature grid point. These results enable users to sample the reaction rate probability density in nucleosynthesis calculations, facilitating uncertainty estimates of nuclidic abundances. The rates presented here refer to their laboratory values. For use in stellar model simulations, these values need to be corrected for the effects of thermal excitations of the interacting nuclei. For each reaction, we include graphs that illustrate the fractional contributions to the overall reaction rate along with the associated uncertainty. These visuals are designed to assist both stellar modelers and nuclear experimentalists by identifying the primary sources of rate uncertain=^texttx);ty at specific stellar temperatures. A graphical comparison with earlier Monte Carlo rates is also provided.

Nuclear astrophysics↗

DOE OSTI2020

Isotopic Evolution of the Inner Solar System Inferred from Molybdenum Isotopes in Meteorites

The fundamentally different isotopic compositions of non-carbonaceous (NC) and carbonaceous (CC) meteorites reveal the presence of two distinct reservoirs in the solar protoplanetary disk that were likely separated by Jupiter. However, the extent of material exchange between these reservoirs, and how this affected the composition of the inner disk, are not known. Here we show that NC meteorites display broadly correlated isotopic variations for Mo, Ti, Cr, and Ni, indicating the addition of isotopically distinct material to the inner disk. The added material resembles bulk CC meteorites and Ca–Al-rich inclusions in terms of its enrichment in neutron-rich isotopes, but unlike the latter materials is also enriched in s-process nuclides. The comparison of the isotopic composition of NC meteorites with the accretion ages of their parent bodies reveals that the isotopic variations within the inner disk do not reflect a continuous compositional change through the addition of CC dust, indicating an efficient separation of the NC and CC reservoirs and limited exchange of material between the inner and outer disk. Instead, the isotopic variations among NC meteorites more likely record a rapidly changing composition of the disk during infall from the Sun's parental molecular cloud, where each planetesimal locks the instant composition of the disk when it forms. Here, a corollary of this model is that late-formed planetesimals in the inner disk predominantly accreted from secondary dust that was produced by collisions among pre-existing NC planetesimals.

58 GEOSCIENCES↗

DOE OSTI2022

The Timing of Potential Last Nucleosynthetic Injections into the Protosolar Molecular Cloud Inferred from 41 Ca– 26 Al Systematics of Bulk CAIs

Short-lived radionuclides (SLRs) provide important information about the chronology of the early solar system. Among them, 41 Ca, due to its decay to 41 K with a half-life of only 0.1 Ma, is particularly valuable in constraining the timescales and origins of both SLRs and the formation of the oldest solar system materials, the Ca–Al-rich inclusions (CAIs). The initial abundance of 41Ca in the solar system, expressed as the ( 41 Ca/ 40 Ca)I ratio, is the key to unveiling the origin of this nuclide. Here, we report a new solar system ( 41 Ca/ 40 Ca)I ratio of 2.0 × 10 −8 derived from the K isotope compositions of two CAIs. This new ratio is about four times higher than the previous value inferred from a mineral isochron. Such a high ( 41 Ca/ 40 Ca)I ratio in the CAIs exceeds that expected for the protosolar molecular cloud by ∼1000×, implying very late injection of the 41 Ca (and possibly other SLRs) into the protosolar molecular cloud. The correlated enrichments of 41 Ca and 26 Al in the bulk CAI samples hint at a common stellar origin of both SLRs. The injection time estimated from our new data depends on the stellar source—it ranges from 0.6 Ma for a Wolf–Rayet wind to 1.0 Ma for a TP-AGB star ejecta.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Strontium-84 Enrichments in Presolar Grains Provide First Evidence of p -process Nucleosynthesis in Core-collapse Supernovae

This study reports detection of rare p-process isotopes within presolar grains. Presolar grains are relic dust grains from dying stars. These microscopic dust particles are found in primitive solar system materials. Their distinct isotopic compositions record the nucleosynthetic processes in their parent stars and the Galactic chemical environment in which these stars formed. We studied presolar graphite grains of high-density type from the Murchison meteorite and found five grains with subgrains that show enrichments in 84 Sr compared to the solar abundance. 84 Sr is the neutron-deficient isotope of strontium that can be produced in the deep oxygen-rich interior of high-mass stars that end their lives as core-collapse supernovae. The observed 84 Sr excesses cannot be produced in low-mass asymptotic giant branch stars, the source of most high-density presolar graphites found in meteorites. High-density graphites with embedded 84 Sr excesses are, instead, compatible with a core-collapse supernovae origin. The graphite subgrains condensed from carbon-rich materials in the outer layers of core-collapse supernovae, where 84 Sr was destroyed by neutron captures during hydrostatic evolution of the stars and their final explosion. Based on current theoretical stellar models, a few percent of contribution from the inner regions of core-collapse supernovae, which are enriched in p-process nuclides, to the outer carbon-rich regions is the most likely explanation for the observed enrichment of 84 Sr in the subgrains of the high-density graphites. In this study, we present the first observational evidence that core-collapse supernovae produce and eject isotopes made by the p-process.

Astronomy and AstroPhysics↗

DOE OSTI2025

East Antarctic Ice Sheet variability in the central Transantarctic Mountains since the mid Miocene

The response of the East Antarctic Ice Sheet to warmer-than-present climate conditions has direct implications for projections of future sea level, ocean circulation, and global radiative forcing. Nonetheless, it remains uncertain whether the ice sheet is likely to undergo net loss due to amplified melting coupled with dynamic instabilities or whether such losses will be balanced, or even offset, by enhanced accumulation under a higher-precipitation regime. The glacial depositional record from the central Transantarctic Mountains (TAM) provides a robust geologic means to reconstruct the past behaviour of the East Antarctic Ice Sheet, including during periods thought to have been warmer than today, such as the mid-Pliocene Warm Period (~3.3–3.0 Ma). This study describes a new surface-exposure-dated moraine record from Otway Massif in the central TAM spanning the last ~9 Myr and synthesises these data in the context of previously published moraine chronologies constrained with cosmogenic nuclides. The resulting record, although fragmentary, represents the majority of direct and unambiguous terrestrial evidence for the existence and size of the East Antarctic Ice Sheet during the last 14 Myr, and it thus provides new insight into the long-term relationship between the ice sheet and global climate. At face value, the existing TAM moraine record does not exhibit a clear signature of the mid-Pliocene Warm Period, thus precluding a definitive verdict on the East Antarctic Ice Sheet's response to this event. In contrast, an apparent hiatus in moraine deposition both at Otway Massif and the neighbouring Roberts Massif suggests that the ice sheet surface in the central TAM was potentially lower than present during the late Miocene and earliest Pliocene.

58 GEOSCIENCES↗

DOE OSTI2022

The collapse of the Cordilleran–Laurentide ice saddle and early opening of the Mackenzie Valley, Northwest Territories, Canada, constrained by 10 Be exposure dating

Deglaciation of the northwestern Laurentide Ice Sheet in the central Mackenzie Valley opened the northern portion of the deglacial Ice-Free Corridor between the Laurentide and Cordilleran ice sheets and a drainage route to the Arctic Ocean. In addition, ice sheet saddle collapse in this section of the Laurentide Ice Sheet has been implicated as a mechanism for delivering substantial freshwater influx into the Arctic Ocean on centennial timescales. However, there is little empirical data to constrain the deglaciation chronology in the central Mackenzie Valley where the northern slopes of the ice saddle were located. Here, we present 30 new 10Be cosmogenic nuclide exposure dates across six sites, including two elevation transects, which constrain the timing and rate of thinning and retreat of the Laurentide Ice Sheet in the area. Our new 10Be dates indicate that the initial deglaciation of the eastern summits of the central Mackenzie Mountains began at ~15.8 ka (17.1–14.6 ka), ~1000 years earlier than in previous reconstructions. The main phase of ice saddle collapse occurred between ~14.9 and 13.6 ka, consistent with numerical modelling simulations, placing this event within the Bølling–Allerød interval (14.6–12.9 ka). Our new dates require a revision of ice margin retreat dynamics, with ice retreating more easterly rather than southward along the Mackenzie Valley. In addition, we quantify a total sea level rise contribution from the Cordilleran–Laurentide ice saddle region of ~11.2 m between 16 and 13 ka.

58 GEOSCIENCES↗

DOE OSTI2023

Reversible ice sheet thinning in the Amundsen Sea Embayment during the Late Holocene

Cosmogenic-nuclide concentrations in subglacial bedrock cores show that the West Antarctic Ice Sheet (WAIS) at a site between Thwaites and Pope glaciers was at least 35 m thinner than present in the past several thousand years and then subsequently thickened. This is important because of concern that present thinning and grounding line retreat at these and nearby glaciers in the Amundsen Sea Embayment may irreversibly lead to deglaciation of significant portions of the WAIS, with decimeter- to meter-scale sea level rise within decades to centuries. A past episode of ice sheet thinning that took place in a similar, although not identical, climate was not irreversible. We propose that the past thinning–thickening cycle was due to a glacioisostatic rebound feedback, similar to that invoked as a possible stabilizing mechanism for current grounding line retreat, in which isostatic uplift caused by Early Holocene thinning led to relative sea level fall favoring grounding line advance.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Assessing the suitability of sites near Pine Island Glacier for subglacial bedrock drilling aimed at detecting Holocene retreat–readvance

Abstract. Unambiguous identification of past episodes of ice sheet thinning below the modern surface and grounding line retreat inboard of present requires recovery and exposure dating of subglacial bedrock. Such efforts are needed to understand the significance and potential future reversibility of ongoing and projected change in Antarctica. Here we evaluate the suitability for subglacial bedrock drilling of sites in the Hudson Mountains, which are located in the Amundsen Sea sector of West Antarctica. We use an ice sheet model and field data – geological observations, glaciological observations and bedrock samples from nunataks, and ground-penetrating radar from subglacial ridges – to rate each site against four key criteria: (i) presence of ridges extending below the ice sheet, (ii) likelihood of increased exposure of those ridges if the grounding line was inboard of present, (iii) suitability of bedrock for drilling and geochemical analysis, and (iv) accessibility for aircraft and drilling operations. Our results demonstrate that although no site in the Hudson Mountains is perfect for this study when assessed against all criteria, the accessibility, N–S orientation and basaltic bedrock lithology of Winkie Nunatak's southernmost ridge (74.86° S, 99.77° W) make it a feasible site both for drilling and subsequent cosmogenic nuclide analysis. Furthermore, the ridge is strewn with glacial erratics at all elevations, providing valuable constraints on its early Holocene deglacial history. Based on our experiences during this study, we conclude with a series of recommendations for assessing site suitability for future bedrock drilling campaigns. We emphasise the importance of consulting a range of expertise prior to drilling and ensuring that sufficient field reconnaissance is undertaken (including obtaining detailed grids of radar survey data and bedrock samples).

58 GEOSCIENCES↗

DOE OSTI2026

Antarctic ice sheet model comparison with uncurated geological constraints shows that higher spatial resolution improves deglacial reconstructions

Accurately reconstructing past changes to the shape and volume of the Antarctic ice sheet relies on the use of physically based and thus internally consistent ice sheet modeling, benchmarked against spatially limited geologic data. The challenge in model benchmarking against geologic data is diagnosing whether model-data misfits are the result of an inadequate model, inherently noisy or biased geologic data, and/or incorrect association between modeled quantities and geologic observations. In this work we address this challenge by (i) the development and use of a new model-data evaluation framework applied to an uncurated data set of geologic constraints, and (ii) nested high-spatial-resolution modeling designed to test the hypothesis that model resolution is an important limitation in matching geologic data. While previous approaches to model benchmarking employed highly curated datasets, our approach applies an automated screening and quality control algorithm to an uncurated public dataset of geochronological observations (specifically, cosmogenic-nuclide exposure-age measurements from glacial deposits in ice-free areas). This optimizes data utilization by including more geological constraints, reduces potential interpretive bias, and allows unsupervised assimilation of new data as they are collected. We also incorporate a nested model framework in which high-resolution domains are downscaled from a continent-wide ice sheet model. We highlight the application of this framework by applying these methods to a small ensemble of deglacial ice-sheet model simulations, and demonstrate that the nested approach improves the ability of model simulations to match exposure age data collected from areas of complex topography and ice flow. We develop a range of diagnostic model-data comparison metrics to provide more insight into model performance than possible from a single-valued misfit statistic, showing that different metrics capture different aspects of ice sheet deflation.

Geosciences↗

DOE OSTI2026

On the Statistical Uncertainty of Monte Carlo-Calculated Scattering Sensitivities

Sensitivity coefficients calculated with Monte Carlo codes are widely used for nuclear data uncertainty quantification in the modeling and simulation of complex 3D reactor systems. This study systematically compares sensitivity coefficients and associated statistical uncertainties for the multiplication factor and fuel temperature reactivity across multiple Monte Carlo codes (SCALE/KENO, SCALE/Shift, MCNP, and Serpent) using simple models representing light-water reactors and advanced reactor concepts. For multiplication factor sensitivities, statistical uncertainties are generally acceptable, although scattering sensitivities show significantly larger statistical uncertainties than, for example, fission and capture reactions. Fuel temperature reactivity sensitivities show significantly larger statistical uncertainties across all reactions. Elastic scattering sensitivities are the most problematic: all Monte Carlo codes fail to resolve energy-dependent coefficients, and they produce dramatically different energy-collapsed values. Critically, the use of these sensitivity coefficients in nuclear data uncertainty propagation leads to reduced statistical uncertainties in individual uncertainty contributions. This can lead to the masking of unusable sensitivity coefficients and producing misleading uncertainty results. The findings of this study show that new or enhanced methods are needed to improve Monte Carlo elastic scattering sensitivity calculations. Additionally, this study shows the relevance of verifying sensitivity coefficients through direct perturbation calculations for individual nuclide reactions, instead of only for total cross sections as commonly done.

Bostelmann, Rike [ORNL] (ORCID:0000000165968088)↗

DOE OSTI2026

Automatic Volume Balancing for Online Refueling in Molten Salt Reactor Simulations with SCALE

This work introduces recently implemented capabilities in the SCALE code system’s TRITON reactor physics sequence that improve molten salt reactor (MSR) modeling: (1) a volume balancing option, which automatically balances the volumes of the fed material with a corresponding material removal to maintain fixed mixture volumes in the neutron transport model and accurate densities, and (2) continuous feed from mixtures, which enables users to define material feed streams directly from salt mixture definitions rather than individual nuclides. The new capabilities were demonstrated via SCALE/TRITON simulations of two representative MSR concepts. Simulations of the 180 MWth molten chloride fast reactor—a system with a fixed fuel salt volume in the reactor core—applied continuous refueling with fresh fuel salt. The results confirm approximately constant reactivity when the new volume balancing capability is used. Additionally, excellent agreement with a manual feed-and-drain method for volume balancing further verified the new implementation. Simulations of the 400 MWth EIRENE reactor, an integral MSR in which the salt volume may grow over time within the reactor core, demonstrated the capability to represent growing salt volume within the TRITON depletion calculation. Compared with reference solutions, the results show consistent trends in reactivity and isotopic evolution, confirming that these new user-friendly capabilities provide accurate, physically consistent approaches for modeling MSRs in SCALE.

Elzohery, Rabab [ORNL] (ORCID:0000000160043633)↗

DOE OSTI2020

Heterogeneity effects on nondestructive assay measurements of enrichment in UF 6 cylinders

A set of measurements were performed using multiple mechanically cooled HPGe detectors at six positions around a cylinder uranium hexafluoride to assess if matrix inhomogeneity is detectible and effects the determined apparent uranium enrichment. These measurements were performed on nine 30B cylinders of uranium hexafluoride. Six of the cylinders appeared to be homogeneous and had similar apparent measured uranium enrichment, as calculated by the FRAM software, at all positions. However, three cylinders appeared to have local inhomogeneity based on the observed results. This was manifested as very low apparent enrichments, often <10% of the declared enrichment, on select positions. The side with a low apparent enrichment shows elevated daughter isotope activity (e.g., elevated 1001 keV and 766 keV peak count rate) and slightly reduced 235 U activity (lower 185.7 keV count rate). We hypothesize that these heterogenous cylinders may have experienced asymmetric solar heating, which caused volatile UF 6 to sublime preferentially away from the warmed side. A similar phenomenon was observed at a second facility where dose rate measurements confirm that the “low-enriched side” of cylinders show elevated gamma ray dose rates, likely from the removal of UF 6 attenuation and concentration of daughter-products. Care should be taken during uranium enrichment verification when applying methods that include gamma-rays associated with daughter nuclides to the determination of uranium enrichment on cylinders that are stored in sunlight and have asymmetric dose readings around the cylinder.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Sensitivity Coefficients Calculated for the Prompt Neutron Decay Constant at or Near Delayed Critical [Abstract]

Experimenters at Los Alamos National Laboratory (LANL) measure the prompt neutron decay constant for many experiments at the National Criticality Experiments Research Center (NCERC) to infer reactivity and the effective neutron multiplication factor. These quantities are very important for nuclear criticality safety and validating nuclear data. Uncertainty in measures of criticality of an experimental configuration can be determined prior to physically performing the experiment by applying first order perturbation theory to Monte Carlo codes, such as MCNP®. The first order perturbation theory produces first derivatives of some nuclear parameter to nuclear data (e.g., cross section data). This first derivative is commonly referred to as a sensitivity coefficient. Currently, the MCNP® software has the capability of computing effective neutron multiplication factor sensitivity coefficients to cross section data. This work builds off of this MCNP® capability and the first order perturbation theory to provide a method of calculating sensitivity coefficients for the prompt neutron decay constant at or near delayed critical to cross section data. The prompt neutron decay constant sensitivity coefficient calculated in this work does not depend on any modification of the MCNP® source code. Prompt neutron decay constant sensitivity coefficient calculations can be used to infer reactivity and effective neutron multiplication factor sensitivity coefficient values as well. By investigating the trends of prompt neutron decay constant sensitivity coefficients for nuclide-reaction pairs across energy spectra, experiments can be designed to maximize or minimize the uncertainty in the prompt neutron decay constant in a particular energy region, which can lead to further optimization studies. Prompt neutron decay constant sensitivity coefficients will be calculated for the Jezebel benchmark. Subsequently, these sensitivity coefficients will be used in a data assimilation process to determine if there is or are optimal experiments that can be performed to provide insight into adjustments of uncertain/inaccurate cross section data. Specifically, the effect of the prompt neutron decay constant sensitivity coefficients on the nuclear data-induced uncertainty in the effective neutron multiplication factor will be examined

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Modeling and Simulation of Fuel Burnup in Pebble Bed Reactors

Modeling and simulation of fuel burnup plays important roles in reactor design, operation, safety, and security as well as nuclear material control and accounting (MC&A) [1]. This task is uniquely challenging for pebble bed reactors (PBR) because the pebbles are continuously added and recycled into the reactor core, and their paths through the core are random. To address this problem, we present two simulation models in this paper. Brookhaven National Laboratory (BNL) developed a simple lattice model of a PBR in Serpent software to generate used pebble isotopic concentrations. The benefit of using Serpent software in this specific application is that it helps streamline the data generation process without having to use too many independent software codes in combination to achieve a simple task. For example, transport, burnup and zero power decay can be implemented in a single pass. Three-dimensional core models were developed using Serpent to simulate the burnup process of 5 subject pebbles starting from fresh till they reach nearly target burnup, with each pebble placed in one of the five artificially designated radial channels to capture the changing neutron spectra along the core radius. Equilibrium isotopic concentrations were assumed in all other pebbles in the core. To provide a verification for the Serpent isotope transmutation and decay results, Oak Ridge National Laboratory (ORNL) performed simple SCALE/ORIGEN calculations using the average neutron spectra calculated by Serpent for each of the 5 pebbles. The 252-group neutron spectra from Serpent were then used by ORIGEN to produce the one-group library for depletion and decay calculations. The isotopic concentrations of a few nuclides of interest and neutron and gamma source terms produced from the ORIGEN calculations were compared with the ones from the Serpent calculations. The model simulated in this work was based on the Pebble Bed Modular Reactor (PBMR)-400 design because many data needed for the simulation such as core power profiles, fuel and reflector temperatures, and equilibrium core composition are publicly available. In this paper, we will compare the results between these two approaches and benchmark the results against a set of well-established simulation results for PBMR-400.

Dim, Odera↗