Search NASA⌕ Search

SEARCH · Search NASA

Results for “time-dependent data”

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.

1,706 records · Page 95

FLUX SENSOR EVALUATIONS AT THE ATR CRITICAL FACILITY

The Advanced Test Reactor (ATR) and the ATR Critical (ATRC) facilities lack real-time methods for detecting thermal neutron flux and fission reaction rates for irradiation capsules. Direct measurements of the actual power deposited into a test are now possible without resorting to complicated correction factors. In addition, it is possible to directly measure minor actinide fission reaction rates and to provide time-dependent monitoring of the fission reaction rate or fast/thermal flux during transient testing. A joint Idaho State University /Idaho National Laboratory ATR National Scientific User Facility (ATR NSUF) project was recently initiated to evaluate new real-time state-of-the-art in-pile flux detection sensors. Initially, the project is comparing the accuracy, response time, and long duration performance of French Atomic Energy Commission (CEA)-developed miniature fission chambers, specialized self-powered neutron detectors (SPNDs) by the Argentinean National Energy Commission (CNEA), specially developed commercial SPNDs, and back-to-back fission (BTB) chambers developed by Argonne National Laboratory (ANL). As discussed in this paper, specialized fixturing and software was developed by INL to facilitate these joint ISU/INL evaluations. Calculations were performed by ISU to assess the performance of and reduce uncertainties in flux detection sensors and compare data obtained from these sensors with existing integral methods employed at the ATRC. Ultimately, project results will be used to select the detector that can provide the best online regional ATRC power measurement. It is anticipated that project results may offer the potential to increase the ATRC’s current power limit and its ability to perform low-level irradiation experiments. In addition, results from this effort will provide insights about the viability of using these detectors in the ATR. Hence, this effort complements current activities to improve ATR software tools, computational protocols and in-core instrumentation under the ATR Modeling, Simulation and V&V Upgrade initiative, as well as the work to replace nuclear instrumentation under the ATR Life Extension Project (LEP) and provide support to the ATR NSUF.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Measurement of $D^0-overline{D}^0$ mixing parameters using $D^0 \to K_S^0 \pi^ \pi^-$ and $D^0 \to K_S^0 K^ K^-$ decays

We report a direct measurement of D{sup 0}-{bar D}{sup 0} mixing parameters through a time-dependent amplitude analysis of the Dalitz plots of D{sup 0} {yields} K{sub S}{sup 0}{pi}{sup +}{pi}{sup -} and, for the first time, D{sup 0} {yields} K{sub S}{sup 0}K{sup +}K{sup -} decays. The low-momentum pion {pi}{sub s}{sup +} in the decay D*{sup +} {yields} D{sup 0}{pi}{sub s}{sup +} identifies the flavor of the neutral D meson at its production. Using 468.5 fb{sup -1} of e{sup +}e{sup -} colliding-beam data recorded near {radical}s = 10.6 GeV by the BABAR detector at the PEP-II asymmetric-energy collider at SLAC, we measure the mixing parameters x = [1.6 {+-} 2.3 (stat.) {+-} 1.2 (syst.) {+-} 0.8 (model)] x 10{sup -3}, and y = [5.7 {+-} 2.0 (stat.) {+-} 1.3 (syst.) {+-} 0.7 (model)] x 10{sup -3}. These results provide the best measurement to date of x and y. The knowledge of the value of x, in particular, is crucial for understanding the origin of mixing.

43 PARTICLE ACCELERATORS↗

A COUPLED TH/NEUTRONICS/CRUD FRAMEWORK IN PREDICTION OF CIPS PHENOMENON

A coupled TH/Neutronics/CRUD framework, which is able to simulate the CRUD deposits impact on CIPS phenomenon, was described in this paper. This framework includes the coupling among three essential physics, thermal-hydraulics, CRUD and neutronics. The overall framework was implemented by using the CFD software STAR-CCM+, developing CRUD codes, and using the neutronics code DeCART. The coupling was implemented by exchanging data between softwares using intermediate exchange files. A typical 3 by 3 PWR fuel pin problem was solved under this framework. The problem was solved in a 12 months length period of time. Time-dependent solutions were provided, including CRUD deposits inventory and their distributions on fuels, boron hideout amount inside CRUD deposits, as well as power shape changing over time. The results clearly showed the power shape suppression in regions where CRUD deposits exist, which is a strong indication of CIPS phenomenon.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Dynamic modeling of a single-stage downward firing, entrained flow gasifier

The gasifier is the heart of the integrated gasification combined cycle (IGCC), a technology that has emerged as an attractive alternative to conventional coal-fired power plant technology due to its higher efficiency and cleaner environmental performance especially with the option of CO{sub 2} capture and sequestration. Understanding the optimal performance of the gasifier is therefore paramount for the efficient operation of IGCC power plants. Numerous gasifier models of varying complexity have been developed to study the various aspects of gasifier performance. These range from simple one-dimensional (1D) process-type models to rigorous higher order 2-3D models based on computational fluid dynamics (CFD). Whereas high-fidelity CFD models can accurately predict most key aspects of gasifier performance, they are computationally expensive and typically take hours to days to execute on high-performance computers. Therefore, faster 1D partial differential equation (PDE)-based models are required for use in dynamic simulation studies, control system analysis, and training applications. A number of 1D gasifier models can be found in the literature, but most are steady-state models that have limited application in the practical operation of the gasifier. As a result, 1D PDE-based dynamic models are needed to further study and predict gasifier performance under a wide variety of process conditions and disturbances. In the current study, a 1D transient model of a single-stage downward-fired GE/Texaco-type entrained-flow gasifier has been developed. The model comprises mass, momentum and energy balances for the gas and solid phases. The model considers the initial gasification processes of water evaporation and coal devolatilization. In addition, the key heterogeneous and homogeneous chemical reactions have been modeled. The resulting time-dependent PDE model is solved using the well-known method of lines approach in Aspen Custom Modeler®, whereby the PDEs in the spatial domain are discretized and the resulting differential algebraic equations (DAEs) are then integrated over time using a dynamic integrator. The dynamic response results of the gasifier performance parameters to certain disturbances commonly encountered during practical operation are presented. These disturbances include ramp and step changes to input variables such as coal flow rate, oxygen-to-coal ratio and water-to-coal ratio among others. Comparison of model predictions to available dynamic data will also be discussed.

01 COAL, LIGNITE, AND PEAT↗

Machine Learning and Economic Models to Enable Risk-Informed Condition Based Maintenance of a Nuclear Plant Asset

The primary objective of this research is to address challenges in the implementation of risk-informed, condition-based predictive maintenance (PdM), which reduces operating costs while still maintaining the safety and reliability of commercial nuclear power plants (NPPs). To achieve the objective, risk models are being developed by taking advantage of advancements in data analytics, deep learning, machine learning (ML), and artificial intelligence (AI). The notable outcomes presented in the report include ? Development of a ML models using heterogeneous plant process and vibration data collected at different spatial and temporal resolutions from the Salem?s CWS to diagnose a circulating water pump (CWP) failure based on salient fault signatures. The developed diagnostic models are extendable to other faults associated with CWPs and CWP motors given associated fault signatures. ? Development of a natural language processing (NLP) technique to automatically classify the WO data into different categories. The developed NLP technique was validated on independent WO data. This automates the tedious and time-consuming activity of mining and classifying WOs by subject matter experts. ? Estimation of mean time between downtime (i.e., time duration between time instances when 1 or more CWPs are not available) and developed an approach to establish reliability of CWS components using unstructured WO data along with CWS plant process data. ? Formulation of economic model based on Markov chain models. The parameters of associated with the transition rate between different states of Markov chain models were estimated using WO data. The economic model formulation and discussion captures both time-independent and time-dependent parameter variation, leading to risk-informed decision-making.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Examination of Digital-Delay Rossi-$\alpha$ for 252 Cf-Driven Highly Enriched Uranium Using Organic Scintillators

Neutron noise techniques constitute several analysis methods applicable to non-destructive assay. One technique is the Rossi-α method to calculate the prompt neutron decay constant (α) or its inverse, the prompt neutron period (1/α), for assemblies of fissionable material. This work evaluates a high-data-throughput measurement at the National Criticality Experiments Research Center (NCERC) with organic scintillators measuring a subcritical assembly of highly enriched uranium (HEU) metal (93% 235 U). The assembly comprises hemi shells stacked together to form fully closed shells and is driven by a Cf source at the center. The assembly is a total of 59.85 kg HEU and a k eff of 0.99, according to MCNP6.2 KCODE simulations. Measurements were acquired with a three-by-four array of 5.08-cm-diameter by 5.08-cm-length trans-stilbene crystals 166 cm from the assembly center. Two types of coincident binning methods are used to build the Rossi-α distribution of coincident neutron detections: A) type 1 binning, also known as any-and-all forward time differences, and B) type 1 binning with a digital-delay technique that is analogous to Orndoff’s use of delay cabling in early Rossi-alpha measurements. Method B) disregards same-detector coincidences and once all nearest time coincidences are collected between all detectors after a trigger, a time delay of 0.75 µs is implemented. The measured prompt neutron decay constants for both techniques are calculated from single exponential fits and the two methods are compared. Method A) shows an apparent timing discontinuity near 500 ns time differences. This timing discontinuity interferes with the fitting method used to calculate the prompt neutron decay constant. Method B) mitigates the artificial timing discontinuity and removes the disagreement of the fit. Future work will model the time-dependent detector response to discern why this timing discontinuity occurs, discern how Method B) reduces the amplitude of this timing discontinuity, and apply this method to critical assembly measurements to ultimately confirm the recommendation to use Method B) for high-data-throughput measurements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Level 1 transient model for a molybdenum-99 producing aqueous homogeneous reactor and its applicability to the tracy reactor

Babcock and Wilcox Technical Services Group (B and W) has identified aqueous homogeneous reactors (AHRs) as a technology well suited to produce the medical isotope molybdenum 99 (Mo-99). AHRs have never been specifically designed or built for this specialized purpose. However, AHRs have a proven history of being safe research reactors. In fact, in 1958, AHRs had 'a longer history of operation than any other type of research reactor using enriched fuel' and had 'experimentally demonstrated to be among the safest of all various type of research reactor now in use [1].' A 'Level 1' model representing B and W's proposed Medical Isotope Production System (MIPS) reactor has been developed. The Level 1 model couples a series of differential equations representing neutronics, temperature, and voiding. Neutronics are represented by point reactor kinetics while temperature and voiding terms are axially varying (one-dimensional). While this model was developed specifically for the MIPS reactor, its applicability to the Japanese TRACY reactor was assessed. The results from the Level 1 model were in good agreement with TRACY experimental data and found to be conservative over most of the time domains considered. The Level 1 model was used to study the MIPS reactor. An analysis showed the Level 1 model agreed well with a more complex computational model of the MIPS reactor (a FETCH model). Finally, a significant reactivity insertion was simulated with the Level 1 model to study the MIPS reactor's time-dependent response. (authors)

07 ISOTOPE AND RADIATION SOURCES↗

Comparison of MCNP6 and experimental results for neutron counts, Rossi-{alpha}, and Feynman-{alpha} distributions

MCNP6, the general-purpose Monte Carlo N-Particle code, has the capability to perform time-dependent calculations by tracking the time interval between successive events of the neutron random walk. In fixed-source calculations for a subcritical assembly, the zero time value is assigned at the moment the neutron is emitted by the external neutron source. The PTRAC and F8 cards of MCNP allow to tally the time when a neutron is captured by {sup 3}He(n, p) reactions in the neutron detector. From this information, it is possible to build three different time distributions: neutron counts, Rossi-{alpha}, and Feynman-{alpha}. The neutron counts time distribution represents the number of neutrons captured as a function of time. The Rossi-a distribution represents the number of neutron pairs captured as a function of the time interval between two capture events. The Feynman-a distribution represents the variance-to-mean ratio, minus one, of the neutron counts array as a function of a fixed time interval. The MCNP6 results for these three time distributions have been compared with the experimental data of the YALINA Thermal facility and have been found to be in quite good agreement. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A variable-order time-dependent neutron transport method for nuclear reactor kinetics using analytically-integrated space-time characteristics

A new time-dependent neutron transport method based on the method of characteristics (MOC) has been developed. Whereas most spatial kinetics methods treat time dependence through temporal discretization, this new method treats time dependence by defining the characteristics to span space and time. In this implementation regions are defined in space-time where the thickness of the region in time fulfills an analogous role to the time step in discretized methods. The time dependence of the local source is approximated using a truncated Taylor series expansion with high order derivatives approximated using backward differences, permitting the solution of the resulting space-time characteristic equation. To avoid a drastic increase in computational expense and memory requirements due to solving many discrete characteristics in the space-time planes, the temporal variation of the boundary source is similarly approximated. This allows the characteristics in the space-time plane to be represented analytically rather than discretely, resulting in an algorithm comparable in implementation and expense to one that arises from conventional time integration techniques. Furthermore, by defining the boundary flux time derivative in terms of the preceding local source time derivative and boundary flux time derivative, the need to store angularly-dependent data is avoided without approximating the angular dependence of the angular flux time derivative. The accuracy of this method is assessed through implementation in the neutron transport code DeCART. The method is employed with variable-order local source representation to model a TWIGL transient. The results demonstrate that this method is accurate and more efficient than the discretized method. (authors)

97 MATHEMATICS AND COMPUTING↗

Calculating infinite-medium {alpha}-eigenvalue spectra with a transition rate matrix method

The time-dependent behavior of the energy spectrum in neutron transport was investigated with a formulation, based on continuous-time Markov processes, for computing {alpha}-eigenvalues and eigenvectors in an infinite medium. For this, a research Monte Carlo code called TORTE was created and used to estimate elements of a transition rate matrix. TORTE is capable of using both multigroup and continuous-energy nuclear data, and verification was performed. Eigenvalue spectra for infinite homogeneous mixtures were obtained and an eigenfunction expansion was used to investigate transient behavior of the neutron energy spectrum. (authors)

97 MATHEMATICS AND COMPUTING↗

Pseudo-dynamic Aspen Plus{sup R} Model of the Defense Waste Processing Facility - 18475

The Defense Waste Processing Facility (DWPF) at the Savannah River Site converts legacy high-level waste from the Cold War into borosilicate glass for permanent disposal in a geological repository. The waste sludge is chemically adjusted by the addition of acid, mixed with frit, and concentrated in the DWPF Chemical Process Cell (CPC) to ensure that the final glass has the desired properties. The resulting mixture is fed into the melter where the waste is melted and poured into stainless steel storage canisters. We present an Aspen Plus{sup R} model of the DWPF CPC, which is being developed to aid chemical processing decisions pertaining to the current flowsheet and to aid decisions regarding optimizing future processing. Given a waste feed composition with acid and frit additions, the model will predict the sludge and off-gas compositions throughout the processing, as well as the final melter feed composition. This allows different processing scenarios to be investigated easily, reducing the reliance on experimental testing, and hence saving time and money. The model has multiple potential applications, including calculating the acid and frit additions required to give acceptable glass properties, off-gas permitting, tracking mercury steam-stripped from the waste in the off-gas system, tracking flammable gases in the offgas system, investigating the effect of reducing air purges in the condensers, and general process optimization. The waste is processed in the CPC in batches, requiring dynamic simulation to track the time-dependent operations. Accurate modeling of the processing also requires calculations of phase and electrolyte chemical equilibria and parallel rate-dependent reactions. Due to existing technical limitations with dynamic simulation of the complex chemical system, we have adopted a pseudo-dynamic approach in which time intervals are treated as consecutive steady-state operations in Aspen Plus{sup R}. Thus, the model consists of ∼100 replicas of the processing vessel and its off-gas system, with the sludge output by one replica (representing one time interval) becoming the feed to the next (representing the next time interval). By examining the results from these replicas in sequence, one can see the time evolution of the model. This general approach can be used wherever simultaneous calculation of phase and electrolyte chemical equilibria in a dynamically evolving system is needed. In addition to describing the general structure and chemistry of the model, we present results from tests of the model against experimental data. These data are from laboratory-scale process simulations of the CPC with a nonradioactive waste sludge simulant. These tests include the model's predictions of the offgas composition, the sludge pH, and the final sludge composition. We will also discuss planned future developments of the model. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

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; each has provided valuable data for materials performance, 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 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 to do verification, validation and uncertainty quantification (VVUQ) on codes. 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 VVUQ approach. To overcome this problem, an Autocorrected Reactimeter Algorithm (ACRA) has been developed. This approach infers 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 the development of an improved version of the ACRA which does not require frozen configurations

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

(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↗

Constraints on Dynamical Dark Energy from Multiple Probes in the Full Dark Energy Survey

We present results on dark energy evolution, assuming a time-dependent equation of state $w(a)=w_0+w_a(1-a)$, from growth and geometric probes using the full six-year Dark Energy Survey dataset: type Ia supernovae, baryon acoustic oscillations, and weak gravitational lensing and galaxy clustering (3$\times$2pt). The combination yields $w_0=-0.84^{+0.10}_{-0.10}$ and $w_a=-0.44^{+0.60}_{-0.55}$, the tightest constraints ever obtained from a single survey, with $2.2σ$ deviation from a cosmological constant. Adding the DESI DR2 BAO data yields $w_0=-0.84^{+0.06}_{-0.07}$ and $w_a=-0.53^{+0.33}_{-0.28}$, representing the most stringent low-redshift-only test of dynamical dark energy to date, with a $2.3σ$ deviation. In this combination, adding 3$\times$2pt doubles the constraining power. Finally, when combined with primary CMB information, we obtain $w_0=-0.82^{+0.05}_{-0.05}$, $w_a=-0.63^{+0.21}_{-0.18}$, with a $3.0σ$ deviation. We find that including 3$\times$2pt in the previously studied SN + DESI BAO + CMB combination leaves the significance essentially unchanged ($3.2 σ$ to $3.0σ$) while improving the figure of merit by $\sim$10%. We systematically investigate the impact of leaving out each one of the probes and find that the significance of the deviation from a cosmological constant ranges from 2.3 to 3.2$σ$, with best-fit parameters consistently in the region $w_0 >-1$ and $w_a <0$. Excluding SN from the all data combination yields a $2.6σ$ departure from $Λ$CDM, providing a cross-check independent of supernova photometric calibration. These results support the weak preference for evolving dark energy reported by several recent cosmological analyses. By combining growth and geometric probes from a single survey, this work realizes the multi-probe dark energy program envisioned at the inception of DES.

Abbott, T. M.C. [Cerro-Tololo InterAmerican Obs.] ↗