SEARCH · Engineering Papers

Results for “average neutron multiplicity”

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

All sources NASA NTRS DOE OSTI

NASA NTRS1972

Effect of multiple spin species on spherical shell neutron transmission analysis

A series of Monte Carlo calculations were performed in order to evaluate the effect of separated against merged spin statistics on the analysis of spherical shell neutron transmission experiments for gold. It is shown that the use of separated spin statistics results in larger average capture cross sections of gold at 24 KeV. This effect is explained by stronger windows in the total cross section caused by the interference between potential and J(+) resonances and by J(+) and J(-) resonance overlap allowed by the use of separated spin statistics.

Semler, T. T.↗

DOE OSTI2021

Which nuclear data can be validated with LLNL pulsed-sphere experiments?

It is shown here that 14-MeV D+T LLNL pulsed-sphere experiments bring complementary information into the process of validating nuclear data compared to experiments that are traditionally used for this purpose—such as critical assemblies. To be more specific, the 14-MeV D+T LLNL pulsed-sphere neutron-leakage spectra enable to validate scattering and fission nuclear data up to 15 MeV (compared to approximately up to 5 MeV when using criticality experiments) and employ to this end simple compound targets containing only few isotopes. In this work, sensitivity profiles of the spectra to nuclear data are calculated in order to understand in detail which isotopes, observables, and energy ranges of nuclear data contribute significantly to their simulation. These profiles are presented for a few selected spheres containing 16 O, 12 C, 56 Fe, and 239 Pu. It is shown that the neutron-leakage spectra of spheres containing light elements are mostly sensitive to elastic- and inelastic-scattering cross sections on discrete levels and corresponding angular distributions. Spheres of structural materials are sensitive to elastic- and inelastic-scattering cross sections, including scattering on discrete levels and the continuum, and double-differential cross sections. Actinide spheres are also strongly sensitive to the fission observables, in particular to the total-fission neutron spectrum. Thin spheres (in which neutrons experience on average less than one scatter) are mostly sensitive to data near the elastic peak, in the energy range from 12–15 MeV, while thicker ones can be sensitive to data at lower incident-neutron energies due to multiple-scattering effects. This information is brought together with simulations of 71 pulsed-sphere neutron-leakage spectra using the ENDF/B-VII.1 and ENDF/B-VIII.0 nuclear-data libraries. This analysis highlights ENDF/B-VIII.0 data that could be further investigated for potential shortcomings ( 6 Li, 12 C, 16 O, 24-26 Mg, 27 Al, 48 Ti, 56 Fe, and 208 Pb) or are likely reliable ( 1,2 H, 7 Li, 9 Be, 14 N, 235,238 U, and 239 Pu) as indicated by validating with LLNL pulsed-sphere experiments.

14-MeV D+T LLNL pulsed-sphere neutron-leakage spec↗

DOE OSTI2020

An Evaluation of a Neutron Time Correlated Interrogation Method for Measurement of Fissile Content in Research Reactor Spent Fuel Assemblies

There are many research reactors worldwide that have been, or are being converted from the use of high enrichment uranium (HEU) to low enrichment uranium (LEU, < 20% enriched). The verification of fissile content and initial enrichment of the spent fuel is needed for the effective safeguards of the fuel. The advanced experimental fuel counter (AEFC) was developed for the measurement of spent fuel rods and assemblies from research reactors for safeguards verification. This measurement system contains components for active neutron interrogation, passive neutron totals counting, neutron coincidence counting, and gross gamma-ray counting. This report presents the first application of the time correlated interrogation technique for the measurement of the 235 U content in research reactor spent fuel assemblies. The technique, called time correlated induced fission (TCIF), uses a 252 Cf neutron source to irradiate the fuel assembly, and the subsequent induced fission events in the fissile material are measured by coincidence counting. The doubles rates are enhanced by having the neutron trigger events from both the 252 Cf source and the induced fission neutrons in the same time gate in the coincidence analysis. The average neutrons per fission of the 252 Cf source is 3.76 and the induced fission neutrons for 235 U is 2.44, so the number of neutrons that are produced is higher than for random neutron interrogation. This high effective neutron number increases the multiplicity counting rates and reduces the statistical error. The background coincidence counts from the 252 Cf are reduced by the water in the sample cavity and the polyethylene surrounding the 3 He detector tubes. This method of active neutron interrogation has been applied to the measurement of spent research reactor (IRT) fuel assemblies. The advanced experimental fuel counter (AEFC) was used to compare the TCIF method with the typically used AmLi neutron interrogation source that emits neutrons that are random in time. The statistical uncertainty for the use of the random neutron source (AmLi) and the time correlated source ( 252 Cf) for spent fuel interrogations was evaluated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

Measuring the neutrino-oxygen neutral current quasielastic cross section using the accelerator neutrino neutron interaction experiment

The Accelerator Neutrino Neutron Interaction Experiment (ANNIE) is a 26-ton gadolinium-doped water Cherenkov detector located on-axis to Fermilab’s Booster Neutrino Beam (BNB). ANNIE is uniquely positioned to perform high-statistics measurements of neutrino-nucleus interactions in water, benefiting from a large neutrino flux due to a short (100-meter) baseline. A central focus of ANNIE’s physics program is the measurement of both charged current (CC) and neutral current (NC) cross sections on water, including neutral current quasielastic (NCQE) and CC-inclusive channels. The NCQE measurement is particularly critical for constraining uncertainties in rare-event searches such as the Diffuse Supernova Neutrino Background (DSNB), where atmospheric $\nu$NCQE interactions constitute a significant and poorly constrained background. This dissertation presents a measurement of the flux-averaged neutrino-oxygen neutral current quasielastic ($\nu$NCQE) cross section using $2.573 \times 10^{20}$~POT of BNB exposure from the 2022 and 2023 beam years. The $\nu$NCQE interaction is identified through the primary $\gamma$-rays produced by nuclear de-excitation of the residual $^{15}$N$^*$ or $^{15}$O$^*$ nucleus following nucleon knockout from $^{16}$O. A dedicated Monte Carlo (MC) re-tuning campaign was conducted using an americium-beryllium (AmBe) calibration source, Michel electrons from stopped muons, and throughgoing dirt muons originating upstream of the detector. This multi-sample approach provided a wide-ranging $\mathcal{O}(\text{MeV})$--$\mathcal{O}(\text{GeV})$ dataset for tuning the simulated detector response, which was subsequently validated against AmBe neutron and Michel electron data for use in the $\nu$NCQE analysis. A dedicated laser calibration campaign was carried out to reduce timing uncertainties across the PMT system, enabling reconstruction of the BNB bunch substructure with sufficient resolution to serve as a background rejection tool. By selecting events in-time with individual neutrino bunches, beam-correlated $\nu$NCQE events are separated from diffuse and accelerator-induced backgrounds, notably skyshine neutrons and externally-originating events, that would otherwise dominate traditional charge-based selections within a small-scale, surface-level, short-baseline detector. A data-driven estimation of the skyshine neutron and external background rates was performed and incorporated into the systematic uncertainty budget. The flux-averaged $\nu$NCQE cross section on oxygen is measured to be $1.57 \pm 0.06\,(\text{stat.})$ $^{+0.91}_{-0.67}\,(\text{syst.})$ $\times 10^{-38}\ \text{cm}^{2}$. A full systematic budget is constructed by propagating uncertainties in the secondary hadronic interaction modeling, background cross section normalizations, detector response, neutrino flux, and the primary $\gamma$-ray emission probabilities from oxygen nuclear de-excitation. An idealized de-excitation model, constructed from existing measurements in the literature is developed to benchmark the predictions of the \textsc{GENIE} event generator. A comparison reveals that \textsc{GENIE} systematically overpredicts the primary $\gamma$-ray emission probability from oxygen de-excitation by a factor of $1.49\times$ for $E_\gamma > 6$~MeV and $3.07\times$ in the $3$--$6$~MeV band. This comparison motivates the dominant systematic uncertainty in this analysis, where a conservative uncertainty of $^{+39.9\%}_{-0\%}$ on the primary $\gamma$-ray signal prediction is assigned. The ANNIE result is consistent with and complementary to existing flux-averaged $\nu$NCQE cross section measurements from T2K and Super-Kamiokande, providing an independent measurement with a different detector, neutrino beam, and analysis methodology. Looking ahead, an upgrade to the ANNIE DAQ infrastructure enabling continuous extended readout will allow a complementary $\nu$NCQE neutron multiplicity measurement, directly relevant to constraining the NCQE background in DSNB searches, competitive with the recent T2K measurement at SK-Gd. The planned Super-SANDI upgrade, deploying a large Water-based Liquid Scintillator (WbLS) volume, will further extend ANNIE's reach to hadronic final states and exclusive NC channels, and enable joint measurements with liquid argon detectors sharing the BNB beamline ahead of DUNE and Hyper-Kamiokande.

Doran, Steven [Iowa State U.]↗

DOE OSTI2021

Extension of the Hauser-Feshbach fission fragment decay model to multichance fission

The Hauser-Feshbach fission fragment decay model, HF 3 D, which calculates the statistical decay of fission fragments, has been expanded to include multichance fission, up to neutron incident energies of 20 MeV. The deterministic decay takes as input prescission quantities—fission probabilities and the average energy causing fission—and postscission quantities—yields in mass, charge, total kinetic energy, spin, and parity. From these fission fragment initial conditions, the full decay is followed through both prompt and delayed particle emissions, allowing for the calculation of prompt neutron and γ properties, such as multiplicity and energy distributions, both independent and cumulative fission yields, and delayed neutron observables. Here, we describe the implementation of multichance fission in the HF 3 D model, and show an example of prompt and delayed quantities beyond first-chance fission, using the example of neutron-induced fission on 235 U. This expansion represents significant progress in consistently modeling the emission of prompt and delayed particles from fissile systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Revised 207 Pb Thermal Capture Gamma-ray Spectrum Evaluation

This report details the revision of the JENDL-4.0 thermal capture gamma-ray spectrum evaluation for 207 Pb and the comparison of the current evaluations and the revised spectrum. The history of the available reaction evaluations is described and the deficiencies are identified. The modified spectrum has an average energy of 7.329 MeV and an average multiplicity of 1.005, leading to a total energy released value of 7.366 MeV, just 0.03% off from the Q-value of the reaction, 7.368 MeV. Based on comparison to the measurements and structure libraries considered, the revised evaluation is closer to reality and should provide more accurate coupled neutron-gamma ray calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Measurement of the 239 Pu$(n, f$) prompt fission neutron spectrum from 10 keV to 10 MeV induced by neutrons of energy 1–20 MeV

Although the prompt fission neutron spectrum (PFNS) is an essential component of neutron-driven systems that has been measured for decades, there are still multiple glaring unknowns regarding the PFNS of major actinides in the fission neutron incident energy range, specifically with regard to multichance fission and pre-equilibrium neutron emission processes. The only impactful experimental 239 Pu PFNS measurements included in recent nuclear data evaluations were measured over a limited outgoing neutron energy range at thermal and 1.5-MeV average incident neutron energy, while other potentially impactful measurements have been shown to contain errors that resulted in either large uncertainty increases or in complete exclusion from nuclear data evaluation. Here, we report a measurement of the 239 Pu PFNS over a wide range of incident neutron energy (1–20 MeV) and three orders of magnitude in outgoing neutron energy (0.01–10 MeV) resulting from the Chi-Nu experiment at the Los Alamos Neutron Science Center. These results are the combination of separate PFNS measurements in the same experimental area, one using a Li-glass and the other a liquid scintillator detector array. Covariances between all PFNS data points from each detector and within each incident energy range were generated between all other data in both detector arrays and within all other incident neutron energy bins, yielding a single covariance matrix for all 1300 PFNS data points reported here. These covariances are based on a thorough assessment of systematic bias and uncertainties associated with the measurement, PFNS extraction technique, combination of data from each detector type, and other aspects of the analysis. The existence of covariances between PFNS data points in different incident neutron energy ranges yielded covariances between average PFNS energy values at each incident energy to be reported here as well, which allowed for firm statements to be made regarding a shape of a purely experimental mean PFNS energy trend for the first time. Although minor PFNS shape differences exist between the results reported here and recent nuclear data evaluations, the ENDF/B-VIII.0 and JEFF-3.3 PFNS evaluations agree reasonably well with the present results from 1-to 10-MeV incident neutron energy, which spans the well-measured 1.5-MeV incident neutron energy PFNS from Lestone and Shores as well as the onset of second-chance fission. However, while the pre-equilibrium component of the PFNS above 12-MeV incident neutron energy roughly agrees in position and magnitude with ENDF/B-VIII.0 and JEFF-3.3, clear differences relating to the relative magnitude of third-chance fission PFNS features are present in the PFNS shape and in the mean PFNS energy trends.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Evaluation of the spatial self-shielding impact for TRISO-based nuclear fuel depletion

Reactor physics analyses of nuclear cores with nuclear fuel concepts containing tristructural isotropic (TRISO) particles, such as pebbles or compact fuel elements, rely on various degrees of simplification to keep these highly heterogeneous problems computationally tractable. One such limitation regards the level of spatial discretization employed during burnup calculations, where traditionally only a limited number of spatial zones are modeled at the full core level and assume that the spectrum is constant within the fuel elements and TRISO particles in this depletion zone. This type of assumption neglects the impact of spatial self-shielding effect within the kernels (microscale level) as well as within the compact or pebbles (mesoscale level). Furthermore, the Monte Carlo code Serpent 2 contains many relevant features for efficiently modeling this type of geometry, including a collision-based domain decomposition intended for very large burnup calculations, which we leveraged for this work to quantify the impact of capturing neutron flux variations occurring at the micro- and mesoscale level on a series of high-temperature gas-cooled reactor fuel element depletion problems. While spatial self-shielding is observed at both scales, with differences from a volume-averaged burnup of ±7% within the kernels and ±2% between TRISO particles within the fuel element, the conjugated effect on nuclide inventories and multiplication factor are negligible, hence confirming that assuming a single average spectrum value may be sufficient for most applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1981

SAS 3 observations of GX 1 + 4

GX 1 + 4 is one of the brightest celestial sources of high-energy X-rays. It is a pulsar with a period of approximately 2 min (perhaps a multiple of 2 min), decreasing at a variable rate which, since 1971, has averaged approximately 2% per year, but which can be larger than 5% per year. This is the largest rate of decrease observed for any pulsar. The rate of decrease appears to be correlated with the luminosity, in support of the idea that the period decrease is produced by accretion torques acting upon a neutron star. No evidence is seen for a Doppler shift due to motion of the pulsar in a binary orbit; this is consistent with the results of optical observations which suggest that any orbital period is fairly long (months to years). The spectrum of GX 1 + 4 is measured as a function of pulse phase, as well as the phase-averaged total spectrum, and the average spectrum of the pulses alone. The shape of the average pulsed spectrum suggests that the pulsations may be produced by 'hot spots' which are a few hundred meters in extent, with temperatures of approximately 10 to the 8th K (kT being approximately equal to 8 keV).

Doty, J. P.↗

DOE OSTI2022

Model calibration of the liquid mercury spallation target using evolutionary neural networks and sparse polynomial expansions

The mercury constitutive model predicting the strain and stress in the target vessel plays a central role in improving the lifetime prediction and future target designs of the mercury targets at the Spallation Neutron Source. We leverage the experiment strain data collected over multiple years to improve the mercury constitutive model through a combination of large scale simulations of the target behavior and the use of machine learning tools for parameter estimation. We present two interdisciplinary approaches for surrogate-based model calibration of expensive simulations using evolutionary neural networks and sparse polynomial expansions. The newly calibrated simulations achieve 7% average improvement on the prediction accuracy and 8% reduction in mean absolute error compared to previously reported reference parameters, with some individual sensors experiencing up to 30% improvement. The calibrated simulations can aid in fatigue analysis to estimate the mercury target lifetime, which reduces abrupt failure and saves tremendous amount of costs.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2021

Progress on the reevaluation and validation of the n+233U neutron cross sections

The set of 233 U resonance parameters of the ENDF/B-VIII.0 nuclear data library was adopted from the previous ENDF/B-VII.1 evaluation using the external levels to update the thermal values. Adoption of IAEA 2017 thermal standards ( σ f = 533.0 ± 2.2 b, σ c = 44.9 ± 0.9 b, and ν ‾ tot = 2.487 ± 0.011 ) and of the IAEA-recommended thermal-neutron induced prompt fission neutron spectrum (PFNS) with average PFNS energy of 2.030 ± 0.013 MeV requires a re-evaluation of 233 U neutron cross sections in the resolved resonance region. A newly produced evaluation is being tested on benchmarks carefully selected from the Handbook of International Criticality Safety Benchmark Experiments (ICSBEP) which are highly sensitive to 233 U data. An important goal of this work was to eliminate the strong negative gradient of the calculated effective multiplication factors with respect to the epithermal fission fraction observed in the validation of the ENDF/B-VIII.0 library for those assemblies. A significant improvement in integral performance of critical 233 U solutions is observed for the newly proposed evaluation. Further work addressing the fast neutron region is needed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Correlations between energy and $\gamma$-ray emission in 239 Pu($n, f$)

Here, we study $\gamma$-ray emission following $^{239}\mathrm{Pu}(n,\mathrm{f})$ over an incident neutron energy range of $2 < E_i < 40$ MeV. We present the first experimental evidence for positive correlations between the total angular momentum generated in fission and the excitation energy of the compound nucleus prior to fission. The $\gamma$-ray multiplicity increases linearly with incident energy below the 2\textsuperscript{nd}-chance fission threshold with a slope of $0.085 \pm 0.010$~MeV$^{-1}$. This linear trend appears to hold for the average excitation energy of the compound nucleus between $9 < \langle E_x \rangle < 19$ MeV. Most of the multiplicity increase comes from an enhancement around a $\gamma$-ray energy of 0.7 MeV, which we interpret as stretched quadrupole $\gamma$ rays that indicate an increase in total fission-fragment angular momentum with excitation energy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Multiple Critical Unresolved Region Integral Experiment (MCURIE): A Proposed Integral Critical Experimental Framework for Unresolved Region and Intermediate Energies

The Critical Unresolved Region Integral Experiment (CURIE) critical experiment was performed at the National Criticality Experiments Research Center (NCERC) at the Device Assembly Facility (DAF) at the Nevada Nuclear Security Sites (NNSS) in 2020. The objective of CURIE was to improve the quality of integral nuclear data in the uranium-235 ( 235 U) unresolved resonance region (URR) by performing benchmark integral experiments that were sensitive to the URR energy ranges. The CURIE experiment was evaluated for the International Criticality Safety Benchmark Evaluation Project (ICSBEP) handbook and the benchmark evaluation was accepted in 2022. The URR is a region within the intermediate neutron energy range (the intermediate energy ranges from 0.7 eV to 100 keV). The observed resonance structure in neutron cross sections is due to discrete energy levels in the nucleus and are characterized by resonance parameters. In the URR region, the resonance parameters are only partially resolved as the resolution of the experiment techniques becomes comparable to the average width of the resonances themselves, and the resonances are so close to one another that the structure cannot be determined empirically. The ENDF/B-VIII.0 and JEFF-3.3 nuclear data libraries define the URR as beginning at 2.25 keV and continuing until 25 keV. There are minimal intermediate neutron energy benchmarks available in the ICSBEP benchmark handbook, and aside from CURIE there are none that are highly sensitive in the URR energy region. There is a current need for additional experiments sensitive to the URR energy region. There is a new proposed subgroup for the Organization for Economic Co-operation and Development and Nuclear Energy Agency (OECD NEA) Working Party on International Nuclear Data Evaluation Cooperation (WPEC), so any new models or information will need experimental validation and testing. The NEA Working Party on Nuclear Criticality Safety (WPNCS) recent experimental needs and priority list includes intermediate energy 235 U and 238 U experiments, as does the recent Integral Experiments to Address Nuclear Criticality Safety Needs Meeting (May 2023). A variant on the original CURIE experiment called the Multiple Critical Unresolved Region Integral Experiment (MCURIE) is proposed to provide further investigation of the intermediate and URR energy region for uranium. MCURIE will utilize existing fuels at NCERC but use alternative moderators and reflectors to modify the neutron absorption, scattering, and fission spectra of the experiments, allowing for precise targeting of nuclear data sensitivities. The overall goal of MCURIE is to develop a framework for designing and performing integral benchmark experiments with high sensitivities in the intermediate and URR energy regions for nuclear data validation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

A suite of neutron time-of-flight detectors to measure hot-spot motion in direct-drive inertial confinement fusion experiments on OMEGA

A suite of six neutron time of flight detectors capable of measuring the primary DT neutron energy spectrum has recently been completed on the OMEGA Laser System. The detectors are positioned along multiple quasi-orthogonal lines of sight (LOS’s) in either a stand-alone, collinear, or antipodal configuration. The collinear detector configuration enables one to measure the neutron velocity along the detector LOS independent of the neutron-production history, while the antipodal detector configuration enables one to directly measure of the hot-spot velocity along the detector axis and Gamow velocity shift due to the plasma ion temperature. Using these six measurements of the primary DT neutron energy spectrum, the neutron-averaged hot-spot velocity and Gamow velocity shift of DT fusion neutrons have been measured in cryogenic experiments on OMEGA for the first time. Finally, the detector suite is operable between DT neutron yields of 1 × 10 13 to 2 × 10 14 with an uncertainty of ~14 km/s in the neutron-averaged hot-spot velocity measurement.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2026

Taming nuclear mass models with Gaussian processes

We propose a new set of nuclear mass predictions based on multiple theoretical mass models. By employing Gaussian process regression with the Matérn kernel, we achieved root-mean-square (rms) deviations below 100 keV for the training dataset. The best-performing mass models achieved rms deviations below 150 keV for the new precise mass data from AME2020, whereas the ensemble average showed robust performance across the nuclear chart. Our approach uniquely combines: (1) systematic refinement of eight mass models through their residuals, (2) physics-informed features, including magic numbers, nucleon parity numbers, neutron excess, and nuclear collectivity, and (3) theory-to-theory validation demonstrating robust extrapolation capability. We find that the Matérn kernel provides superior uncertainty quantification compared to the RBF kernel, with a length-scale analysis revealing enhanced inter-nuclei correlations. We provide complete mass predictions for all unknown nuclides in AME2020, offering valuable constraints for nuclear structure studies and astrophysical modeling when used with proper uncertainty propagation.

Gaussian processes↗

NASA NTRS2001

The X-Ray Globular Cluster Population in NGC 1399

We report on X-ray sources detected in the Chandra images of the elliptical galaxy NGC 1399 and identified with globular clusters (GCs). The 8'x 8' Chandra image shows that a large fraction of the 2-10 keV X-ray emission is resolved into point sources, with a luminosity threshold of 5 x 10 (exp 37) ergs s-1. These sources are most likely Low Mass X-ray Binaries (LMXBs). More than 70% of the X-ray sources, in a region imaged by Hubble Space Telescope (HST), are located within GCs. Many of these sources have super-Eddington luminosity (for an accreting neutron star) and their average luminosity is higher than the remaining sources. This association suggests that, in giant elliptical galaxies, luminous X-ray binaries preferentially form in GCs. The spectral properties of the GC and non-GC sources are in most cases similar to those of LMXBs in our galaxy. Two of the brightest sources, one of which is in GC, have a much softer spectra as seen in the high state black hole. The "apparent" super-Eddington luminosity in many cases may be due to multiple LMXB systems within individual GC, but with some of the most extreme luminous systems containing massive black holes.

Angelini, Lorella↗

DOE OSTI2023

The Extension of the Hauser-Feshbach Fission Fragment Decay Model to Multi-chance Fission and its Application to 239 Pu

The Hauser-Feshbach fission fragment decay model, HF3D, calculates the statistical decay of fission fragments through both prompt and delayed neutron and γ-ray emissions in a deterministic manner. While previously limited to the calculation of only first-chance fission, the model has recently been extended to include multi-chance fission, up to neutron incident energies of 20 MeV. The deterministic decay takes as input prescission quantities–fission probabilities, pre-fission neutron energies, and the average energy causing fission– and post-scission quantities–yields in mass, charge, total kinetic energy, spin, and parity. From those fission fragment initial conditions, the full decay is followed through both prompt and delayed particle emissions. The evaporation of the prompt neutrons and γ rays is calculated through the Hauser-Feshbach statistical theory, taking into account the competition between neutron and γ-ray emission, conserving energy, spin, and parity. The delayed emission is taken into account using time-independent calculation using decay data. This whole formulation allows for the calculation of prompt neutron and γ-ray properties, such as multiplicities and energy distributions, both independent and cumulative fission yields, and delayed neutron observables, in a consistent framework. Here, we describe the implementation of multi-chance fission into the HF 3 D model, and show an example of prompt and delayed quantities beyond first-chance fission, using the example of neutron-induced fission on 239 Pu. This expansion represents significant progress in consistently modeling the emission of prompt and delayed particles from fissile systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Automated Control for Nuclear Thermal Propulsion Start-Up using MOOSE-based Applications

This report presents a Griffin/Bison/RELAP-7 numerical model of a prototypical NTP system that features fuel assemblies arranged in rings, and which was designed to simulate rapid startup transients. The physics modeled include full-core neutronics, assembly-wise heat conduction, and conjugate heat transfer, with the balance of plant mainly imposed through boundary conditions. In addition, various forms of automated reactivity control were deployed by using the MOOSE to autonomously drive the model and simulate the reactor transitioning from assumed initial conditions to nominal power in a fraction of a minute. To generate the cross-sections of the neutronics model, and in an effort to simultaneously account for the tremendous axial temperature gradients in the reactor and to limit the number of state points required for cross-section generation, the average component temperatures and hydrogen densities in the cooling channels were correlated to the average fuel and moderator temperatures, and fixed axial profiles were derived for nominal conditions and then used during the transient. With this approximation, a tractable cross-section library tabulated with fuel/moderator temperatures and CD angles was generated using Serpent. The full-core SPH correction procedure and the CD decusping technology in Griffin, respectively, ensure preservation of the multiplication factor and reaction rates at state points, along with a reasonably accurate reactivity worth between tabulated CD angles, despite using a coarse mesh. Feedback from other physics was calculated by modeling one representative fuel assembly per ring, along with the corresponding fuel and moderator cooling channels. To limit power overshoots during startup, another layer of multiphysics coupling was added to the model in order to automatically control the drums. Two different technologies presented herein showed outstanding performance in this regard: (1) a novel hybrid PID controller based on both power and reactivity signals, and (2) a PGC that relies on kinetics parameters and reactivity coefficients to predict future behavior and adjust the desired signal accordingly. A challenging benchmark was devised, featuring a power demand curve that exponentially increases by a factor of 500 within 30 seconds, then levels out after that. Both control approaches create a simulated power curve that closely follows the power demand curve and limits power overshoots to 1% or less. While the former approach requires more tuning of the internal parameters, the latter requires additional knowledge of the reactivity feedback coefficients and rates of change of the corresponding variables, including fuel and moderator temperature, which could be difficult to dynamically measure for a real NTP system. Fortunately, some inaccuracy in these quantities will not drastically degrade the PGC performance. Subsequently, a more realistic startup sequence was considered, in which the mass flow rate and outlet pressures are ramped up to model bootstrap and thrust build-up phases prior to reaching steady-state conditions, demonstrating the ability of the hybrid PID and PGCs to handle such transients, with both types of controllers exhibiting very similar behavior. Nevertheless, a significant chamber temperature overshoot was observed, caused by the demanded power signal and assumed mass flow rate. This issue could be mitigated by deploying a reactor controller that follows the chamber temperature signal and actuates both the control valves and drums (rather than using a power signal based solely on the drums to control reactivity). Enhancement of the hydrogen fluid properties available in MOOSE, as well as a better understanding of prototypical initial conditions, are also needed to further enhance this startup model. Finally, a study was performed to model decay heat post-shutdown, and to prepare for extending this model to predict shutdown behavior and post-shutdown pulsed cooling requirements.

33 ADVANCED PROPULSION SYSTEMS↗