Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reactor Physics Method”

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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Development of a continuous synthesis process for carbamazepine using validated in-line Raman spectroscopy and kinetic modelling for disturbance simulation

Mitigation of failure modes in the continuous synthesis (CS) of a drug substance (DS) has the potential to widen the adoption of continuous manufacturing (CM) technologies by the pharmaceutical industry. Here, this work demonstrates the development of a robust continuous process for the synthesis of carbamazepine (CBZ), an essential medicine as per the World Health Organization (WHO), facilitated by kinetic modelling and monitored by in-line Raman spectroscopy. Accurate kinetic modelling and the use of validated process analytical technology (PAT) models for quantitative measurement were found to play an important role in developing CS of drug substances. Kinetic data for the formation of CBZ from iminostilbene (ISB) were collected by batch reaction sampling and high-performance liquid chromatography (HPLC) analysis. A non-linear solver and iterative method was applied to determine two sets of Arrhenius parameters simultaneously for the reaction system by minimizing the standard error of the model fit. The start-up and dynamic equilibrium stages for the CS of CBZ using a continuous stirred tank reactor (CSTR) were modelled based on the batch kinetic data and employed to optimize conversion and simulate process disturbances. An in-line Raman spectroscopy method was successfully developed, validated, and integrated to determine the concentrations of CBZ and ISB within the operating range for the CS. The CS kinetic model was evaluated experimentally from startup to dynamic equilibrium over 10 residence times with monitoring by HPLC and in-line Raman spectroscopy. The developed kinetic model in tandem with in-line Raman spectroscopy successfully predicted disturbances due to changes in process variables and can serve as a useful tool in the future design of advanced process control strategies for the continuous synthesis of CBZ.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Investigations of Properties of {sup 89}Y Levels with Reactor Fast Neutrons

By employing a beam of reactor fast neutrons, the spectrum of gamma rays up to an energy of 4.6MeV and their angular distributions with respect to the neutron-beamaxis aremeasured in the reaction {sup 89}Y(n, n'γ). The multipolarities and multipole-mixture parameters for 34 gamma transitions and the spin–parities J{sup π} of states excited in this reaction are determined. The lifetimes of the lowest 32 levels of 89Y were measured by the Doppler shift attenuation method, and the reduced probabilities for the respective gamma transitions were calculated. Levels of the K{sup π} = +5/2{sup +} and K{sup π} = −7/2{sup +} bands associated with, respectively, prolate and oblate deformation shapes are found in {sup 89}Y at low excitation energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Determination of plasma density from data on the ion current to cylindrical and planar probes

To improve probe methods of plasma diagnostics, special probe measurements were performed and numerical models describing ion transport to a probe with allowance for collisions were developed. The current–voltage characteristics of cylindrical and planar probes were measured in an RF capacitive discharge in argon at a frequency of 81 MHz and plasma densities of 10{sup 10}–10{sup 11} cm{sup –3}, typical of modern RF reactors. 1D and 2D numerical models based on the particle-in-cell method with Monte Carlo collisions for simulating ion motion and the Boltzmann equilibrium for electrons are developed to describe current collection by a probe. The models were used to find the plasma density from the ion part of the current–voltage characteristic, study the effect of ion collisions, and verify simplified approaches to determining the plasma density. A 1D hydrodynamic model of the ion current to a cylindrical probe with allowance for ion collisions is proposed. For a planar probe, a method to determine the plasma density from the averaged numerical results is developed. A comparative analysis of different approaches to calculating the plasma density from the ion current to a probe is performed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Rancor-HUNTER: A Virtual Plant and Operator Environment for Predicting Human Performance

Advances in simulation capabilities to model physical systems have outpaced the development of simulations for humans using those physical systems. There is an argument that the infinite span of potential human behaviors inherently render human modeling more challenging than physical systems. Despite this challenge, the need for modeling humans interacting with these complex systems is paramount. As technologies have improved, many of the failure modes originating from the physical systems have been solved. This means the overall proportion of human errors has increased, such that it is not uncommon to be the primary driver of system failure in modern complex systems. Moreover, technologies such as automated systems may introduce emerging contexts that can cause new, unanticipated modes of human error. Therefore, it is now more important than ever to develop models of human behavior to realize overall system error reductions and achieve established safety margins. To support new and novel concepts of operations for the anticipated wave of advanced nuclear reactor deployments, human factors and human reliability analysis researchers need to develop advanced simulation-based approaches. This talk presents a simulation environment suitable to both collect data and then perform Monte Carlo simulations to evaluate human performance and develop better models of human behavior. Specifically, the Rancor Microworld Simulator models a complex energy production system in a simplified manner. Rancor includes computer-based procedures, which serve as a framework to automatically classify human behaviors without manual, subjective experimenter coding during scenarios. This method supports a detailed level of analysis at the task level. It is feasible for collecting large sample sizes required to develop quantitative modelling elements that have historically challenged traditional full-scope simulator study approaches. Additionally, the other portion of this experimental platform, the Human Unimodel for Nuclear Technology to Enhance Reliability (HUNTER), is presented to show how the collected data can be used to evaluate novel scenarios based on the contextual factors, or performance shaping factors, derived from Rancor simulations. Rancor-HUNTER is being used to predict operator performance with new procedures, such as results from control room modernization or new-build situations. Rancor-HUNTER is also proving a useful surrogate platform to model human performance for other complex systems.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

(Doublon) Benchmarking of Different Inverse Point Kinetics Implementations for an Autocorrected Reactimeter Algorithm

In November 2017, the Transient Reactor Test Facility returned to operation. Since that time, many transient test series have been completed, such as the Transient Heatsink Overpower Response capsule (THOR), the Transient Water Irradiation System for TREAT (TWIST), and Sirius. Each has provided valuable data for materials performance and reactor safety that can be applied in future designs. During each experimental series, detector count rates provided important information on the core behavior during transients. However, a limitation of these data is that variations in the neutron distribution during experiments can cause errors when attempting to infer reactivity evolution from detector signals. Neutron physics codes can be used to compute the flux shape variations. However, this is a poor solution when the experimental data is used for code verification, validation and uncertainty quantification. Indeed, if the output of the code is used both as a reference and to correct what the reference is compared to, the circular dependency limits the quality of the verification, validation and uncertainty quantification approach. To overcome this problem, the autocorrected reactimeter algorithm (ACRA) has been developed. This approach infers a time-dependent reactivity evolution by testing different spatial corrections and selecting the one that minimizes reactivity variations when the core is in a frozen configuration (i.e., when there is no variation in parameters affecting reactivity). However, the scope of this method was limited to transients where there were negligible thermal feedback. Indeed, the core is never in a frozen configuration when the fuel temperature varies during the whole transient. This is our motivation for developing an improved version of the ACRA that does not require frozen configurations. To develop this new algorithm, we need a precise and unbiased implementation of the inverse point kinetic equations (IPKEs) as any error in the reactivity evaluation will be propagated into the choice of the optimal spatial correction. Indeed, the previous reactimeter algorithm would use approximations, such as a negligible flux amplitude derivative, to focus on rapidity. For the numerical validation of ACRA, we aim at absolute error under for reactivity derived from signals similar to the one of this study. In this summary, we test eight different IPKE implementations. Each will process a mockup signal built for this study, similar to those that the future ACRA will process. Each reactivity output will be compared to the reference reactivity that has been used to generate the mockup signal. The implementation minimizing the difference with the reference reactivity will be used in the development of a new ACRA formulation.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

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

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

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Analysis of the Daya Bay Reactor Antineutrino Flux Changes with Fuel Burnup

We investigate the recent Daya Bay results on the changes in the antineutrino flux and spectrum with the burnup of the reactor fuel. We find that the discrepancy between current model predictions and the Daya Bay results can be traced to the original measured 235 U/ 239 Pu ratio of the fission beta spectra that were used as a base for the expected antineutrino fluxes. An analysis of the antineutrino spectra that is based on a summation over all fission fragment beta-decays, using nuclear database input, explains all of the features seen in the Daya Bay evolution data. However, this summation method still predicts an anomaly. Thus, we conclude that there is currently not enough information to use the antineutrino flux changes to rule out the possible existence of sterile neutrinos.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Radiological Safety Analysis Computer (RSAC) Program Version 7.4 Users’ Manual

The Radiological Safety Analysis Computer (RSAC) Program Version 7.4 (RSAC-7) is the newest version of the RSAC legacy code. RSAC-7 calculates the consequences of a release of radionuclides to the atmosphere. Users generates a fission product inventory from either reactor operating history or a nuclear criticality event. RSAC-7 models the effects of high-efficiency particulate air filters or other cleanup systems and calculates the decay and ingrowth during transport through processes, facilities, and the environment. Doses are calculated for inhalation, air immersion, ground surface, ingestion, and cloud gamma pathways. RSAC-7 is used as a tool to evaluate accident conditions in emergency response scenarios, radiological sabotage events, and safety basis accident consequences. This users’ manual contains the mathematical models and operating instructions for RSAC-7. Instructions, screens, and examples are provided to guide the user through the functions provided by RSAC-7. This program is designed for users who are familiar with radiological dose assessment methods.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Robust Online Monitoring for Calibration Assessment of Transmitters and Instrumentation

Robust online monitoring (OLM) technologies are expected to enable the extension or elimination of periodic sensor calibration intervals in operating and new reactors. These advances in OLM technologies will improve the safety and reliability of current and planned nuclear power systems through improved accuracy and increased reliability of sensors used to monitor key parameters. In this article, we discuss an overview of research being performed within the Nuclear Energy Enabling Technologies (NEET)/Advanced Sensors and Instrumentation (ASI) program, for the development of OLM algorithms to use sensor outputs and, in combination with other available information, 1) determine whether one or more sensors are out of calibration or failing and 2) replace a failing sensor with reliable, accurate sensor outputs. Algorithm development is focused on the following OLM functions: • Signal validation • Virtual sensing • Sensor response-time assessment These algorithms incorporate, at their base, a Gaussian Process-based uncertainty quantification (UQ) method. Various plant models (using kernel regression, GP, or hierarchical models) may be used to predict sensor responses under various plant conditions. These predicted responses can then be applied in fault detection (sensor output and response time) and in computing the correct value (virtual sensing) of a failing physical sensor. The methods being evaluated in this work can compute confidence levels along with the predicted sensor responses, and as a result, may have the potential for compensating for sensor drift in real-time (online recalibration). Evaluation was conducted using data from multiple sources (laboratory flow loops and plant data). Ongoing research in this project is focused on further evaluation of the algorithms, optimization for accuracy and computational efficiency, and integration into a suite of tools for robust OLM that are applicable to monitoring sensor calibration state in nuclear power plants.

fault detection↗
DOE OSTI2021

Prediction of the Thermal Runaway Limit and Optimal Operation of Heat Transfer-Limited, Fixed-Bed Reactor Systems

Here, we derive a new prediction for thermal runaway starting from the alpha model for fixed-bed reactor systems. This method accounts for thermal resistance internal to the reactor tube and the radial temperature gradients that result. To showcase our method, we compare its predictions to other common criteria for thermal runaway using o-xylene oxidation as the example chemistry. Even in systems where internal heat transfer is negligible, the empirical practical design criterion for thermal runaway is inaccurate. For cases where internal heat transfer is relevant, our runaway limit is more stringent than limits derived from simpler 1-D models. To augment our work, we optimize the product yield with the thermal runaway constraint using orthogonal collocation. Using the alpha model, the results illustrate that the thermal runaway limit can be accurately determined using either numerical or analytical methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Transport and kinetics properties of LaF 3 in FLiNaK molten salt determined by electrochemical methods

It is found that lanthanum species can exist in the form of LaF 6 3- in FLiNaK (46.5 LiF-11.5 NaF-42 KF, mol%) molten salt in previous studies. The understanding of transport and reaction kinetics parameters of LaF 6 3- ions in FLiNaK molten salt is essential for the electrochemical separation of lanthanum from FLiNaK melt to achieve the coolant clean-up in fluoride salt cooled high-temperature reactor (FHR). Through conducting chronopotentiometry tests, the diffusion coefficients of LaF 6 3- ions in FLiNaK molten salt at the temperature range of 923 K–1023 K were determined. By conducting potentiodynamic polarization measurements at different concentrations and temperatures, the exchange current density, reaction rate constant, and charge transfer coefficient were obtained. The experiment data of potentiodynamic polarization were analyzed using a non-simplified electrode kinetics equation which incorporates both mass transfer and reaction kinetics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

The Flash ADC system and PMT waveform reconstruction for the Daya Bay experiment

To better understand the energy response of the Antineutrino Detector (AD), the Daya Bay Reactor Neutrino Experiment installed a full Flash ADC readout system on one AD that allowed for simultaneous data taking with the current readout system. This paper presents the design, data acquisition, and simulation of the Flash ADC system, and focuses on the PMT waveform reconstruction algorithms. For liquid scintillator calorimetry, the most critical requirement to waveform reconstruction is linearity. Several common reconstruction methods were tested but the linearity performance was not satisfactory. A new method based on the deconvolution technique was developed with 1% residual non-linearity, which fulfills the requirement. Here, the performance was validated with both data and Monte Carlo (MC) simulations, and 1% consistency between them has been achieved.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Advanced Modeling and Simulation Methods for Evaluation of Thermal Neutron Scattering Materials

With the rise of interest in thermal neutron scattering data for advanced reactor, criticality safety, and shielding applications, new experimental data are required for evaluation of new materials or for re-evaluation (or validations) of previously evaluated materials. New experimental data are evaluated in a three-step process: (1) computing the phonon characteristics, (2) computing the dynamic structure factor (DSF) from the data, and (3) using the experimental setup to simulate the experimental data. All three steps have challenges, ranging from the need for a sufficiently general material simulation code—a processing code that can compute the corresponding DSF—to having a detailed layout of the instrument/beamline/facility where the data were measured. Whereas phonon characteristics of materials can be computed using various methods (molecular dynamics, density functional theory, etc.), a high-fidelity computation of the DSF and the simulation of the experiment based on the DSF is vital to the accuracy of the evaluation. The latter two steps can be achieved by using the two corresponding code systems developed by instrument scientists at the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory: (1) OCLIMAX, a program that calculates the dynamic structure factor from DFT and MD simulation results, and (2) MCViNE, a Monte Carlo neutron ray-tracing program designed to simulate neutron scattering experiments. Recently, polyethylene and yttrium hydride were measured at the Wide Angular-Range Chopper (ARCS) and SEQUOIA instrument stations of the SNS. These experiments are simulated using the density functional theory code, the Cambridge Serial Total Energy Package (CASTEP), to compute its phonon characteristics (eigenvalues/vectors and PDOS), which is then processed using OCLIMAX to yield the DSF, and finally the data at each instrument station are simulated by the MCViNE for comparison to the measured data for evaluation. For comparison to conventional evaluation methods, the scattering data processed from OCLIMAX are compared against those processed from the LEAPR module of NJOY, and the results from MCViNE simulations are compared against previously used simplified beamline models implemented in the Monte Carlo N-Particle (MCNP) code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Combining CO 2 Reduction with Ethane Oxidative Dehydrogenation by Oxygen-Modification of Molybdenum Carbide

The surface properties that determine the selectivity of Mo 2 C catalysts in ethane oxidative dehydrogenation with CO 2 as a soft oxidant were investigated using a combination of pulse experiments and in-situ spectroscopic methods. Oxygen modification was discovered to be crucial for inhibiting the cleavage of the C–C bond in ethane and enhancing the production of ethylene. Here, the addition of the Fe promoter accelerated the formation of surface oxygen species and stabilized them from reduction by ethane, leading to a shorter induction period, higher ethylene yield, and improved stability.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

First observation of unlocking the locked mode by electrode biasing on J-TEXT tokamak

Major disruptions, often caused by locked m/n = 2/1 modes (m/n is the poloidal/toroidal mode number), are great threat to a tokamak fusion reactor and should be mitigated or avoided. The locked 2/1 modes have been unlocked by electrode biasing (EB) for the first time in J-TEXT tokamak experiments. Furthermore, the application of a sufficiently negative EB voltage in plasma edge region quickly changes the phase and amplitude of the locked 2/1 mode and drives the mode to rotate in about 10–30 ms. For a larger EB current or a smaller locked mode amplitude, the mode is more easily to be unlocked, revealing a new method for mode unlocking and avoiding plasma disruptions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Automated X-Ray Diffraction of Irradiated Materials

Synchrotron-based X-ray diffraction (XRD) and small-angle Xray scattering (SAXS) characterization techniques used on unirradiated and irradiated reactor pressure vessel steels yield large amounts of data. Machine learning techniques, including PCA, offer a novel method of analyzing and visualizing these large data sets in order to determine the effects of chemistry and irradiation conditions on the formation of radiation induced precipitates. In order to run analysis on these data sets, preprocessing must be carried out to convert the data to a usable format and mask the 2-D detector images to account for experimental variations. Once the data has been preprocessed, it can be organized and visualized using principal component analysis (PCA), multi-dimensional scaling, and k-means clustering. In conclusion, from these techniques, it is shown that sample chemistry has a notable effect on the formation of the radiation induced precipitates in reactor pressure vessel steels.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

A consistent, Bayesian, approach to the cross section probability distribution in the unresolved resonance region

The cross sections of neutron-induced reactions can be divided into three energy ranges: the resolved resonance region (RRR), the unresolved resonance region (URR), and the fast region. In general, the cross sections in the URR show significant fluctuations that cannot be predicted and cannot be experimentally resolved, thus, it is commonly assumed that the cross section at a specific energy is given by a probability distribution function (PDF) over a range of values that can span several orders of magnitude. The current methodology used to describe such behavior is to construct the PDF by stochastically generating resonance ladders and numerically measuring the PDF. The resonance ladders are sampled using known resonance statistical properties and average resonance widths and spacings extrapolated from the RRR. Although this is a standard and widely used technique, it is computationally very expensive, therefore, an alternative, analytical, approach would be preferable due to the considerable speed up of the computational time in real life applications. Moreover, the current methodology does not take into account existing experimental data, such for total and capture cross sections, that are available for many nuclei. Finally, this approach was developed to be used in reactor-scale applications and it is not suited for use in single-event applications. In this work we will rethink the entire approach to the PDF construction using a Bayesian mindset. This will allow us to provide a different definition of the PDF that allows a much faster calculation of the higher-temperature PDFs and a proper combination of theoretical and experimental PDFs following the probability theory. We will also show that our definition is well suited for single-event applications and we will make an explicit connection between our method and the standard approach. We do this by showing that the central limit theorem applies and our method leads to the same PDF obtained with the standard methodology, for a large number of events per history.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Wet to Dry Controls Lanthanide Scandate Synthesis

The choice of temperature and gas conditions used in a water pressure-controlled reactor is guided by density functional theory (DFT) to synthesize nearly phase pure lanthanide scandate nanoparticles (LnScO 3 , Ln = La, Nd, Sm, Gd). In this synthetic method, low water vapor partial pressures, well below water’s gas liquidus, inhibit particle growth, while an excess of water vapor results in undesired rare earth hydroxide and oxyhydroxide secondary phases. The optimal humidity for high-purity LnScO 3 particle synthesis is shown to vary with the lanthanide; DFT is used to calculate the thermodynamics of secondary phase formation for each lanthanide tested, such that the role of water vapor may be quantified and used to maintain phase purity (greater than 96 mol%) across the series. Furthermore, the combination of thermodynamic calculation and experimental confirmation with this pressure-controlled reactor provides an opportunity to explore analogous syntheses of other inorganic perovskite nanoparticles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗