Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Delayed Neutron Fraction”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Validation of KENO Delayed Neutron Fraction Capabilities

The delayed neutron fraction, which is necessary for the safe operation of reactors, varies by fissionable material and neutron spectrum. A prompt fission neutron spectrum has been included in the SCALE 6.2 release to account for the removal of delayed neutrons from precursor materials from the neutron balance. This paper describes how experimental measurements from criticality safety benchmarks were used to validate the SCALE calculational procedure for the delayed neutron fraction. Additionally, the effect of nuclear data uncertainty on the KENO calculation of β eff was examined with the SAMPLER sequence.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Data Assimilation using Non-invasive Monte Carlo Sensitivity Analysis of Reactor Kinetics Parameters [Slides]

Sensitivity coefficients of prompt neutron decay constant (α) and effective delayed neutron fraction (β eff ) to Pu-239 nuclear data were calculated using a newly available tool (ACEtk). Uncertainty calculations showed the following trends: (1) The prompt neutron decay constant can be used to reduce nuclear data-induced uncertainty in Pu-239(n,f) and (2) the effective delayed neutron fraction can be used to effectively reduce nuclear data-induced uncertainty in Pu-239 total fission $\tilde{v}$. Sensitivity/uncertainty analysis can be used to determine optimal experiment/detector setup to identify compensating errors in nuclear data. Upcoming work includes investigating other nuclide-reaction pairs and response sensitivities, such as Δ$\rho$, neutron leakage spectra, and reaction rate measurements to constrain nuclear data of interest.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Estimating β eff from static, late-time Feynman neutron noise measurements

The effective delayed neutron fraction, β eff , is a kinematic parameter describing the contribution of delayed neutrons to the effective multiplication factor, $\kappa$ eff , and is of great interest to the nuclear criticality safety community and beyond. In fact, $\kappa$ eff presently must be inferred (as opposed to directly measured) since variables like β eff must be obtained from simulation, reference tables, or a dynamic or pulsed measurement where the assembly of fissile material is perturbed and not in steady state. Several recent works have proposed neutron noise techniques to estimate β eff from static measurements (such as benchmark experiments). This work experimentally investigates an approach based on looking at late-time (up to thousands of seconds) correlated neutrons in Feynman histograms. This work utilizes measurements of weapons grade plutonium and highly enriched uranium to construct late-time Feynman histograms and investigate a proposed methodology of estimating β eff from a fit of the distribution. Data is shown from both He-3 detectors and an organic scintillator array. This work is preliminary in nature, demonstrating some experimental feasibility and proposing improvements for future investigations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Preservation of kinetics parameters generated by Monte Carlo calculations in two-step deterministic calculations

The generation of accurate kinetic parameters such as mean generation time Λ and effective delayed neutron fraction β eff via Monte Carlo codes is established. Employing these in downstream deterministic codes warrants another step to ensure no additional error is introduced by the low-order transport operator when computing forward and adjoint fluxes for bilinear weighting of these parameters. Another complexity stems from applying superhomogenization (SPH) equivalence in non-fundamental mode approximations, where reference and low-order calculations rely on a 3D full core model. In these cases, SPH factors can optionally be computed for only part of the geometry while preserving reaction rates and K-effective, but the impact of such approximations on kinetics parameters has not been thoroughly studied. This paper aims at studying the preservation of bilinearly-weighted quantities in the Serpent–Griffin calculation procedure. Diffusion and transport evaluations of IPEN/MB-01, Godiva, and Flattop were carried out with the Griffin reactor physics code, testing available modeling options using Serpent-generated multigroup cross sections and equivalence data. Verifying Griffin against Serpent indicates sensitivities to multigroup energy grid selection and regional application of SPH equivalence, introducing significant errors; these were demonstrated to be reduced through the use of a transport method together with a finer energy grid.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Utilization of ACE nuclear data file toolkit ACEtk to calculate relative sensitivity coefficients of point-kinetics parameters

Sensitivity and uncertainty methods are quintessential for nuclear criticality safety and experiment design. This type of analysis relies on calculations of sensitivity coefficients; sensitivity coefficients of the effective neutron multiplication factor with respect to some nuclear data are predominantly calculated and used. As a part of the Laboratory Directed Research & Development project EUCLID (Experiments Underpinned by Computational Learning for Improvements in nuclear Data) at Los Alamos National Laboratory, sensitivity coefficients of many radiation detector measurement responses with respect to nuclear data were investigated. Specifically, this paper outlines a method to calculate point-kinetics parameters relative sensitivity coefficients with respect to nuclear data. Point-kinetics parameters such as the prompt neutron decay constant, effective delayed neutron fraction, and neutron generation time are especially important to experimenters and reactor operators designing systems with dynamic neutron populations. This method couples capabilities of the ACE (A Compact ENDF) nuclear data file toolkit, ACEtk, with the ability to load cross sections into the radiation transport code Monte Carlo N-Particle (MCNP). In conclusion, key aspects of optimizing this method for a particular application and sensitivity profiles of the Jezebel criticality experiment are examined and discussed.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2022

Nuclear data uncertainty propagation applied to the versatile test reactor conceptual design

We report the Versatile Test Reactor (VTR) currently under development is a 300 MWth sodium-cooled fast reactor (SFR) fueled with ternary metal alloy fuel, which aims to accelerate the testing of advanced nuclear fuels, materials, instrumentation, and sensors in high flux environments that are necessary to license the next generation of advanced reactor concepts. To support the VTR design process, uncertainties associated with the nuclear data has been propagated through the reactor core neutronics calculation to global parameters of interest, such as the core multiplication factor, kinetic parameters, and various reactivity feedback coefficients, following the sensitivity based uncertainty propagation approach. By folding the sensitivity coefficients, separately computed by the generalized perturbation theory code PERSENT and Monte Carlo code Serpent 2, with the variance-covariance matrices from COMMARA-2.0, we obtain the reaction-wise, isotope-wise, and overall uncertainties for each response of interest due to nuclear data uncertainty. With Serpent 2, the statistical error of the uncertainty is obtained by propagating the statistical error of the sensitivity coefficients through the same process using a newly developed uncertainty propagation method. From both codes, the overall top uncertainty contributors are found to be the cross section of Fe-56 elastic scattering, Na-23 elastic scattering, and U 238 inelastic scattering. The large contributions of the Fe-56 elastic scattering cross sections to global parameters are due to its relatively large relative uncertainty of 5–10% in nuclear data and the large volume of Fe-containing reflector assemblies in the fairly compact VTR core design. Both codes agreed well for the overall uncertainty estimates of all responses of interest, except the delayed neutron fraction, prompt neutron generation time, and the coolant density feedback coefficient, where Serpent 2 yielded a much larger value than PERSENT due to the large statistical error of sensitivity coefficients. The calculated uncertainties are also compared to those associated with other SFR cores. Another outcome of this study is a variance-covariance matrix of reactivity coefficients, which can be used in the subsequent uncertainty propagation to the system level to investigate the impact of identified uncertainties on system responses in the safety analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Using Prompt Neutron Decay Constant Measurements to Obtain Additional Kinetics Information

Prompt neutron decay constant measurements are typically used to infer the prompt neutron lifetime in a chain-reacting system, but these measurements include much more information about the kinetics of the system. This work focuses specifically on the Rossi-α method of measuring the prompt neutron decay constant. The Rossi-α system of equations consists of five unknown parameters. With five unknowns and three equations, the system of equations solved by the Rossi-α measurement is under-determined. It has been postulated that a well posed problem can be formulated using independent measurement configurations on a single system. The external source strength, the effective delayed neutron fractions, and Diven’s parameter do not typically change for different configurations of a single experiment. If good detector geometry is chosen, the detector efficiency does not change for multiple measurements of a single experiment using Rossi-α measurements. By making these assumptions, there are only two unique variables. This knowledge can be used to solve for most of the unique parameters using two independent, unique Rossi-α measurements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Data Assimilation using Non-invasive Monte Carlo Sensitivity Analysis of Reactor Kinetics Parameters

Accurately predicting the criticality of an experiment before interacting with the experimental components is very important for criticality safety. Radiation transport software can be utilized to calculate the effective neutron multiplication factor of a nuclear system. Because of the integral nature of the effective neutron multiplication factor, the value calculated contains various sources of nuclear-data induced uncertainty. The sensitivity analysis and data assimilation technique presented in this paper exhibit one possible method of identifying and reducing the effective neutron multiplication factor nuclear-data induced uncertainty. The results presented in this work show that it is possible to use relative sensitivity coefficients of the prompt neutron decay constant and the effective delayed neutron fraction to 239 Pu nuclear data to reduce nuclear-data induced uncertainties in the effective neutron multiplication factor. This work has been utilized by members of the Los Alamos National Laboratory project EUCLID (Experiments Underpinned by Computational Learning for Improvements in Nuclear Data) for optimally designing a new experiment, which will be used to reduce compensating errors in 239 Pu nuclear data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Effects of Uranium Impurities in Downblended HEU on HTGR Performance

Many advanced reactor designs require fuel enriched between 5% and 20% 235 U. To assist in producing fuel at these enrichment levels, government-owned inventories of highly enriched uranium can be downblended. However, fuel produced from these inventories contain uranium impurities that are not often found when enriching natural uranium or accounted for when modeling reactor cores. To address this concern, this work models reactor designs like the X-energy Xe-100 and the Ultra Safe Nuclear Company’s Micro Modular Reactor, and compares their performance with fuel from enriching natural uranium to fuel from downblended highly enriched uranium. This paper evaluates the models based on the effective neutron multiplication factor, k eff , effective delayed neutron fraction, β eff , and energy- and spatially dependent neutron flux, ϕ, as well as the fuel, coolant, moderator, and total reactivity temperature feedback coefficients, α F , α C , α M , and α T . The results show that the fuel from downblended highly enriched uranium inventories leads to differences in each of the metrics, especially in the keff values. In the Xe-100–like and Micro Modular Reactor–like models, k eff changes by about 1400 pcm and up to 1200 pcm, respectively. Total reactivity feedback coefficients α T are negative with the impure fuels and the keff values remain above 1 for each core configuration and fuel composition. These results show that the impure fuel compositions do not necessarily prevent achieving key design parameters, such as cycle length, or from operating in a safe condition.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Demonstrate Capability of NEAMS Tools to Generate Reactor Kinetics Parameters for Pebble-Bed HTGRs Transient Modeling

The system analysis of anticipated operating occurrences and design basis accident for pebblebed reactor systems requires knowledge of the neutron kinetic parameters. These parameters take into account various types of feedback from the pebble-bed core and are typically integrated into system analysis tools such as the System Analysis Module (SAM) in the Nuclear Energy Advanced Modeling and Simulation (NEAMS) tool package. These parameters are typically generated by a higher fidelity full-core coupled neutronics/thermal fluids analysis, such as using Griffin and Pronghorn, in 2D or 3D. The delayed neutron fractions (beta) and the neutron lifetimes (lambda) are typically generated by applying the adjoint solution of the neutron flux. Temperature reactivity coefficients are the most important feedback for all anticipated operating occurrences and design basis accident, then depending on the transient other effects might be important as rector control system worth, s-curves, and Xenon generation. This work package will establish a methodology and process for calculating the reactor kinetics parameters for high-temperature gas-cooled reactor and provide them in a suitable form for system level analyses of typical anticipated operating occurrences and design basis accident.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Implementation of Perturbation Theory and Sensitivity Capabilities in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor Multiphysics analysis application, jointly developed by Argonne and Idaho National Laboratories under the DOE-NE NEAMS program. This fiscal year, capabilities for reactivity and sensitivity evaluation using perturbation methods were implemented and verified. The First Order Perturbation Method (FOPT) was employed to compute reactivity worth resulting from small perturbations in input parameters, while the Generalized Perturbation Theory (GPT) was used to evaluate sensitivities of a range of response types, including reaction rate ratio, k-eigenvalue, neutron generation time, and effective delayed neutron fraction. These perturbation methods enable users to quantify how response quantities change due to a perturbation in a input parameter without explicitly performing an additional transport simulation for each perturbed state. In particular, the GPT formulation accounts for indirect effects arising from flux changes by solving generalized inhomogeneous equations, for which a Neumann series-based iterative solution method was developed and implemented in Griffin. The implemented reactivity and sensitivity evaluation capabilities were verified using two test problems: an infinite homogeneous system and a two-dimensional hexagonal core. The results showed excellent agreement with reference solutions obtained by a direct method based on finite difference approximation as well as GPT-based results from the PERSENT code, confirming the accuracy of both reactivity and sensitivity evaluations. Additionally, preliminary uncertainty quantification (UQ) results were obtained by combining the sensitivity values computed using GPT and external covariance data, demonstrating that the implemented sensitivity results can be reliably used for uncertainty calculations. To further demonstrate the generality and practical strength of the implementation, the sensitivity evaluation capability was successfully applied to the Empire microreactor with a geometrically complex design that poses significant modeling challenges. The results confirm that Griffin enables sensitivity evaluations even for irregular and highly heterogeneous reactor configurations, thereby establishing a foundation for UQ applications in advanced reactor designs and analyses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Pile Oscillator Evaluation for Flattop [Poster]

Useful method for high-precision measurements of integral cross section in various samples. The availability of benchmark-quality effective delayed neutron fraction (ßeff) values for nuclear data and code validation remains limited.

MCNP↗
DOE OSTI2025

Uncertainty quantification and sensitivity analysis of a nuclear thermal propulsion reactor startup sequence

The research presented in this article describes progress in applying stochastic methods, uncertainty quantification, parametric studies, and variance-based sensitivity analysis (also known as Sobol sensitivity analysis) to a full-core model of a nuclear thermal propulsion (NTP) system simulated via the radiation transport code Griffin to simulate neutronics. Our goal is to develop a reduced-order (surrogate) model that can be rapidly sampled with perturbations to multiple input parameters. In this NTP system, reactivity and power feedback affect the rotation of control drums (CDs), which is itself controlled by a hybrid proportional-integral-derivative (PID) controller actuated by the power demand and reactivity feedback from the numerical model. This model uses reactor kinetic feedback (mean generation time [Λ] and effective delayed neutron fraction [ β eff ] from a transient Griffin simulation executed via Griffin’s improved quasi-static solver to provide the kinetic parameters) as inputs to functions that control the CD rotation angle. By investigating numerous stochastic approaches, we developed a dual-purpose surrogate model of the NTP system, using polynomial regression in the Multiphysics Object-Oriented Simulation Environment (MOOSE) Stochastic Tools Module (STM). The trained model can be rapidly sampled while simultaneously perturbing various input parameters, such as coefficients on the PID control or temperature (directly affecting the neutron cross section). The surrogate model delivers accurate (within 5%) results at speeds orders of magnitude faster (minutes, not days of computational time) than the base model. Once the surrogate model has been trained, distributions of the uncertain parameters can be changed at will to investigate the effects of perturbing multiple inputs as well as the effects of these inputs on the model output. For example, coefficients used in the PID control system may vary due to some type of physical interference, or uncertainty may exist in the temperature of the neutron cross sections in various regions of the reactor. A distribution can be placed on these parameters, and operational boundaries can be determined. The goal of this work is to support development of an advanced control system for operating CDs in a functioning NTP system. This work is a scoping study of the MOOSE STM.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2022

Evaluation of Oak Ridge National Laboratory Health Physics Research Reactor Operation Data for Critical Benchmark Creation [Abstract]

The Oak Ridge National Laboratory (ORNL) Health Physics Research Reactor (HPRR) was a research reactor designed and built at ORNL in 1961. The critical assembly was using a highly enriched uranium and molybdenum alloy as the fuel, and it could be operated in steady-state or burst modes. The reactor was used for about 25 years to produce a lot of publications related to dosimetry, radiobiology and radiation detectors testing before its decommissioning in 1987. In recent years, the idea of using legacy operation data from the to create a valuable critical accident alarm system shielding benchmark arose. Such a benchmark has been submitted to the International Criticality Safety Benchmark Experiment Project (ICSBEP) Technical Review Group for a potential inclusion in the 2022 version of the handbook. Another way to use the valuable data from the operation of the HPRR is to evaluate the feasibility of the creation of a subcritical or prompt supercritical benchmark for inclusion in the ICSBEP or the International Reactor Physics Experiments Evaluation Project (IRPhEP) handbooks. To initiate a burst, the HPRR had to be operated in a slightly subcritical state for a few minutes. Then, the insertion of the burst control rod would greatly increase the reactivity of the system and start the burst. No critical configuration of the HPRR critical assembly could be located. The only information available concerns stable subcritical and prompt supercritical states, found in a burst experiments’ logbook. In the recovered logbook pages, information about 8 different bursts is available. The information includes the rods positions before and during a burst, the recorded subcritical reactor period and reactivity, and the burst fission yield derived from the temperature elevation sulfur pellet irradiation analysis. By using the HPRR logbook information and the as-built drawings of the critical assembly, a highly detailed model of the HPRR was created with SCALE 6.2.4/KENO-VI. Eight KENO-VI models were created to replicate the sub-critical assembly configurations described in the eight bursts from the recovered logbook pages. KENO-VI calculates k eff and it can be linked to a reactivity value in cents by using the delayed neutron fraction B eff , also calculated by KENO-VI. KENO-VI can also be used to model the prompt super-critical configurations of the HPRR and to assess the similarity with the burst measurements by comparing the calculated k eff and the measured fission yields between each burst. Unfortunately, high uncertainty exist and the obtained discrepancies between experiments and calculation results are high, compromising the creation of a valuable critical benchmark from HPRR operation data. The reasons of the discrepancies and potential ways to solve them are explored.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Investigation of the OECD/NEA PWR MOX/UO{sub 2} core transient benchmark using a coupled whole core pin-by-pin route in WIMS

This paper investigates the new CAMELOT coupling route developed in WIMS, which uses prepared tabulated cross-sections combined with the flux solver MERLIN module and the integrated sub-channel thermal hydraulics solver ARTHUR module to solve for the coupled core state. This coupling route is applied to the OECD/NEA pressurized water reactor (PWR) core transient benchmark exercise, involving 3-D neutronics and thermal-hydraulics analysis of a full-sized PWR core featuring MOX fuel at various core states, as well as simulating a control rod ejection transient. Analysis of the results for each part of the benchmark demonstrate that the high-resolution pin-by-pin analysis used in WIMS closely matches the other participants in the benchmark, with differences of between 20-200 pcm for k-effective, 1-2.5% RMS difference for assembly-averaged powers, 2-8% difference in rod worth calculations, < 200 ppm for calculations of critical boron concentration, zero difference in delayed neutron fraction calculation and a very close match for power response to the rod ejection transient. This is the first application of the WIMS-CAMELOT approach for the transient analysis of a full-sized reactor and as such further improvements such as reducing computational cost through parallelization, optimisation and memory reduction, are currently in active development. This paper demonstrates the accuracy and flexibility of the CAMELOT coupling route, providing a straightforward and consistent user image that can be easily used for a variety of modelling problems involving coupled neutronics and thermal-hydraulics. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗