Search NASA⌕ Search

SEARCH · Search NASA

Results for “deterministic 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.

532 records · Page 30
All sourcesNASA NTRSDOE OSTI
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 OSTI2018

Upgrading Limiting Peak-Power Analysis Techniques with Modern Validation and Uncertainty Quantification for the Advanced Test Reactor

Here, this work demonstrates the acceptability of the 2D deterministic transport code, HELIOS, to replace the legacy diffusion code, PDQ, for computing the peak-power performance parameters of the Advanced Test Reactor (ATR). The 95% Confidence Rule, commonly used in the commercial reactor sector, is explored to develop the so-called “reliability factors” which provide statistical confidence that the peak-power limits within the hottest location along a fuel plate, referred to as the hot-stripe, will not be exceeded. Additionally, an alternative “legacy” methodology was explored that attempts to mimic the exact PDQ analysis process used for defining the peak-power limits. The legacy methodology, involves interpolating power between regions at azimuthal boundaries subtending the regions of interest. In order to apply the 95% Confidence Rule, a statistically significant calculation-to-measurement bias must first be established. Unlike the commercial world where thousands of power observations can be collected every cycle using on-line flux mapping instrumentation, the ATR power distribution must be measured during “depressurized” zero-power measurements using fission wires in polyethylene wands. In 2012, fission wire activation data was collected during a flux run in the Advanced Test Reactor – Critical facility. Also to improve statistical validity, archival data from ATR zero power flux runs from 1977, 1986, and 1994 were digitized from scanned reports and used to create new benchmark models. Borrowing from least-squares adjustment methods commonly used for neutron activation spectroscopy, adjusted fission wire powers were calculated for all four datasets. The mean and standard deviation of the bias between a priori calculated and adjusted wire-powers was then taken as the bias and uncertainty used in the 95% Confidence Rule.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Griffin Software Development Plan

Griffin is a MOOSE-based reactor physics application for advanced reactor multiphysics modeling and simulation. The application is developed in a consistent multiphysics environment with strong software quality assurance. Griffin inherited most of the capabilities of MAMMOTH/Rattlesnake and is adopting the capabilities from PROTEUS that are needed in the code. The toolset includes a variety of deterministic radiation transport solvers for fixed source, k-eigenvalue, ad-joint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. The code contains the cross-section preparation capabilities applicable to fast and thermal reactors, including TRISO-fueled reactors. Core management capabilities include core performance, fuel depletion and shuffling, equilibrium core calculation, pebble-bed reactor run-in and equilibrium core, molten-salt reactor delayed neutron precursor drift, and control rod and drum movement with cusping correction. This software development plan presents the current and future capabilities and features in Griffin for the design and analysis of non-light-water reactor systems in steady-state and transient conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Development of Critical Experiments to Provide Validation Data for Multiphysics Coupling Methods (Final Report)

Computational simulation has become one of the most important tools in the modern toolbox for nuclear reactor design. Reactor simulations, like the reactors themselves, involve multiple physical processes including neutronics, fluid flow and heat transfer (thermal-hydraulics), thermo-mechanics, fuel behavior, chemistry, and balance of plant. Each process can affect the behavior of the others, and advanced computational tools have been developed that couple separate physics models together to simulate the feedback effects. The objective of this work is to generate experimental data to support the validation of multi-physics reactor codes. It specifically targets the SHARP toolset from the NEAMS Reactor Product Line, but it is intended to be useful in any set of coded that couples thermal-hydraulics with neutronics simulations. The project takes advantage of the inherent flexibility of the RPI Reactor Critical Facility (RCF), a low power open-pool research reactor. The standard configuration of the RCF is used to explore the effects of moderator temperature on neutronic behavior and reactivity. Because the low power of the reactor generates no appreciable temperature change in the moderator, the moderator temperature is artificially changed with two 18 kW electric immersion heaters. The temperature feedback effects can be observed by recording the change in the reactor power as the reactor goes from a supercritical state, to critical, to a subcritical state. To expand the temperature range and temperature rate of change that can be investigated, an experimental apparatus consisting of a heated water loop was designed and installed in the RCF. In this setup, water from a heated reservoir tank is pumped through a pipe placed through the center of the reactor core. The water in the reservoir can be heated to around 70 °C and pumped through the system. The pump has a variable flow rate and electronic control, so a number of scenarios can be developed. A slug of water can be pushed through the system to generate a reasonably fast change in reactivity. The pump speed can be increased or decreased over the course of a measurement, or even oscillated at different periods, to observe the response of the reactor. To complement the experimental data generated by the RCF, a set of computational models have been developed. They are intended to illustrate the usefulness of the experimental data generated in validating coupled multi-physics simulations. A model previously constructed in MCNP has been used in conjunction with safety analysis at the RCF for many years and was modified to support the design of the heated water loop experiment. This model was translated into Serpent 2, a similar Monte Carlo code, which provided the added advantage of being able to produce collapsed cross-section data files for use in the deterministic neutron transport codes that are part of the SHARP toolset. A complete reconstruction of the RCF model, both with and without the heated fluid loop experiment, were necessary for use with the SHARP codes, PROTEUS (neutronics) and Nek5000 (thermal hydraulics), as these require an unstructured mesh format to operate. The computational requirements for fully coupled three dimensional simulations are significant, and challenges to full implementation and comparison with RCF experiments remain. A wide variety of experiments have been conducted for eventual benchmarking and validation of multi-physics reactor codes. A complete database of detailed time- and temperature-dependent measurements have been generated and stored. Several examples of the experimental results are presented here, with refinement and analysis of experiments ongoing. While these experiments still only address a finite range of operating conditions, in particular being limited to water in the sub-boiling regime, they represent a unique (to date) look at the feedback between temperature and neutronic behavior suited for validation of the coupling routines in multi-physics codes. It is hoped that these will be the foundation of future benchmarks that will provide confidence in the tools that will be used to develop the next generation of nuclear reactor systems.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Uncertainty quantification of a physics-informed model based on sparse identification of a Thermal Energy Distribution System

Integrated energy systems (IES)s are crucial for enhancing the economy and efficiency of power generation sources (e.g., nuclear energy) necessary to unleash American energy dominance. These systems can be integrated with thermal energy storage (TES) and intermittent renewable energies to optimize overall energy use, peak-load regulation, and demand-side responses. However, the stabilization of energy generation, transport, and utilization introduces operational complexities that exceed the challenges of managing each sub-component individually. Currently, though IESs rely on human operators for efficiency and stability, reducing human error risk and enhancing performance through automation is highly desirable. Recent advances at Idaho National Laboratory have demonstrated successful control of the Thermal Energy Distributed System (TEDS). However, the automatic control system depends on a deterministic Sparse Identification of Nonlinear Dynamics with Control (SINDyC) model, which are trained based on simulation data from physics-based simulations. Because of uncertainties in physics-based simulation, SINDyC model results in large discrepancies against experimental data and cannot be reliably used in automatic control. In this paper, we present an innovative approach to address these discrepancies by quantifying uncertainties and developing a more robust model. We first generated trajectories by using first-principles physics codes to encapsulate the experiment. Next, we trained thousands of models by randomly sampling these trajectories. We then collapsed all those models into one probabilistic SINDyC by fitting a multivariate Gaussian distribution onto the resulting coefficient’s distribution. Despite its simplicity, our approach successfully produced 95% confidence intervals that captured the experimental trajectories. It even did so with a higher probability and better U-pooling score across six of the seven relevant quantities of interest (QoIs), as compared to other classical approaches. In conclusion, ongoing research is focusing on generating new experimental trajectories to validate this approach, and on employing Bayesian calibration to refine parametric uncertainties and guide future model development efforts.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Distributions of the Mean for Lognormally Distributed Data - 19427

In environmental sciences, a distribution of the mean is often needed for right-skewed data. For example, risk assessments are often based on distributions of the mean to obtain upper confidence limits (UCLs) for analyte concentrations, and, in probabilistic performance assessment (PPA) models, parameters are given distributions rather than deterministic values, and they are usually given distributions of their mean rather than the population distributions to account for the large spatial and temporal domains to which they are applied. How this distribution of the mean is developed depends on the characteristics of the parameter's full distribution and how the parameter relates to its use. In this presentation, I will discuss the implications of a right skew in distribution development. Parameters may vary over orders of magnitude and/or may be right-skewed. The material property of saturated hydraulic conductivity, Ks, is a prime example of such a parameter. Right-skewed parameters are often assumed to follow lognormal distributions. If X follows a lognormal distribution, then log(X) follows a normal distribution. For some right-skewed parameters, practitioners go back and forth between the two forms of the data depending on their needs. For example, it may be preferred to use the untransformed version of X in a model or equation, but it is easier to visualize the parameter values on a logarithmic scale. However, distribution development for the average value of the parameter will differ depending on whether X or log(X) is modeled. This difference is due to Jensen's Inequality, a theorem which states that the expected value of a concave function of a variable is less than or equal to the same function of the expected value of the variable; a logarithm is a concave function. In other words, the mean of log(X) is always less than or equal to the log of the mean of X and therefore the mean estimated from a distribution based on log(X) is less than or equal to the mean estimated from a distribution based on X. If a parameter is naturally considered on a log scale by subject matter experts, then distributions developed for the mean of log(X) may be more consistent with the use of the parameters in the literature than distributions developed for X. The statistical decisions of distribution development will be discussed within the context of right-skewed, lognormal data. These issues include method-of-moments versus maximum likelihood estimation of the distribution parameters and the use of bootstrapping to aid in the visualization of a distribution of the mean. The details surrounding Jensen's Inequality will be explained and visualized within the context of distribution development of a right-skewed parameter such as saturated hydraulic conductivity. (authors)

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

The Photometric LSST Astronomical Time-series Classification Challenge PLAsTiCC: Selection of a Performance Metric for Classification Probabilities Balancing Diverse Science Goals

Classification of transient and variable light curves is an essential step in using astronomical observations to develop an understanding of the underlying physical processes from which they arise. However, upcoming deep photometric surveys, including the Large Synoptic Survey Telescope (LSST), will produce a deluge of low signal-to-noise data for which traditional type estimation procedures are inappropriate. Probabilistic classification is more appropriate for such data but is incompatible with the traditional metrics used on deterministic classifications. Furthermore, large survey collaborations like LSST intend to use the resulting classification probabilities for diverse science objectives, indicating a need for a metric that balances a variety of goals. We describe the process used to develop an optimal performance metric for an open classification challenge that seeks to identify probabilistic classifiers that can serve many scientific interests. The Photometric LSST Astronomical Time-series Classification Challenge (PLASTICC) aims to identify promising techniques for obtaining classification probabilities of transient and variable objects by engaging a broader community beyond astronomy. Using mock classification probability submissions emulating realistically complex archetypes of those anticipated of PLASTICC, we compare the sensitivity of two metrics of classification probabilities under various weighting schemes, finding that both yield results that are qualitatively consistent with intuitive notions of classification performance. We thus choose as a metric for PLASTICC a weighted modification of the cross-entropy because it can be meaningfully interpreted in terms of information content. Finally, we propose extensions of our methodology to ever more complex challenge goals and suggest some guiding principles for approaching the choice of a metric of probabilistic data products.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Operator learning for predicting multiscale bubble growth dynamics

We report simulating and predicting multiscale problems that couple multiple physics and dynamics across many orders of spatiotemporal scales is a great challenge that has not been investigated systematically by deep neural networks (DNNs). Herein, we develop a framework based on operator regression, the so-called deep operator network (DeepONet), with the long-term objective to simplify multiscale modeling by avoiding the fragile and time-consuming “hand-shaking” interface algorithms for stitching together heterogeneous descriptions of multiscale phenomena. To this end, as a first step, we investigate if a DeepONet can learn the dynamics of different scale regimes, one at the deterministic macroscale and the other at the stochastic microscale regime with inherent thermal fluctuations. Specifically, we test the effectiveness and accuracy of the DeepONet in predicting multirate bubble growth dynamics, which is described by a Rayleigh–Plesset (R–P) equation at the macroscale and modeled as a stochastic nucleation and cavitation process at the microscale by dissipative particle dynamics (DPD). First, we generate data using the R–P equation for multirate bubble growth dynamics caused by randomly time-varying liquid pressures drawn from Gaussian random fields (GRFs). Our results show that properly trained DeepONets can accurately predict the macroscale bubble growth dynamics and can outperform long short-term memory networks. We also demonstrate that the DeepONet can extrapolate accurately outside the input distribution using only very few new measurements. Subsequently, we train the DeepONet with DPD data corresponding to stochastic bubble growth dynamics. Although the DPD data are noisy and we only collect sparse data points on the trajectories, the trained DeepONet model is able to predict accurately the mean bubble dynamics for time-varying GRF pressures. Taken together, our findings demonstrate that DeepONets can be employed to unify the macroscale and microscale models of the multirate bubble growth problem, hence providing new insight into the role of operator regression via DNNs in tackling realistic multiscale problems and in simplifying modeling with heterogeneous descriptions.

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

Nanoindentation of Micrograins in Polycrystals under Multiaxial Stress: Control of Abrupt & Stochastic Plastic Events

Deforming objects in the world of the small has proven a much more complex task than just bending a common-sized object like a spoon. Over the last decade, it has become clear that the application of stress on micron-sized objects, such as micropillars, leads to abrupt and stochastic plastic strain jumps that do not satisfy Weibull statistics (as expected for uncorrelated defect yielding) but instead seem to follow highly wide, power-law distributions, pointing towards cooperative dislocation plasticity. Are power-law distributed plastic events the rule or the exception in crystals? If the rule, can we design protocols to control and suppress them in order to achieve deterministic and smooth microforming processes? We propose to combine experiments and theory to unravel the mysteries of abrupt microplastic deformation in a natural dislocation environment: micron-sized grains (\micrograins") of polycrystals. We propose a novel experimental design that utilizes nanoindentation and 4-point bending techniques on the substrate of polycrystalline thin films of two typical FCC metals (Al and Cu), and make concrete theoretical conjectures, based on intuition built from previous theoretical and experimental efforts in micropillar compression and nanoindentation. Based on a wide set of observations over the last decade, we put forward two distinct hypotheses that we propose to investigate in both theory and experiments: Our first conjecture is that crystal plasticity in the central region of micron-sized grains has micropillar-like stochastic features: power-law plastic strain jumps at stresses above the macroscale yield strength, displaying strong size and rate effects. We will demonstrate this similarity in spherical nanoindentation experiments in the central-region of multi-orientation grains of Al and Cu thin films. Al and Cu are selected for their disparity in stacking fault energy and consequently, their different propensity for cross-slip relaxation; however, our protocol may be applied to any crystalline material of technological interest. We will test our hypothesis by performing two-dimensional extended discrete dislocation plasticity (2.5D-DDP) simulations, inspired by comparable successful studies of micropillar plasticity, and extended to model nanoindentation with two slip systems and on a surface of a bicrystal with microhard grain-boundaries (GB), ie. not allowing dislocation absorption (116). Our model will include effective three dimensional dislocation mechanisms such as dislocation multiplication and double-cross-slip assisted glide. Second, we hypothesize that the application of uniform in-plane stress on our thin films below the yield point can smoothen the micrograin plastic response during nanoindentation. On a technological level, we propose that indentation-induced localization (sink-in) may be diminished by the application of concurrent in-plane stress that activates slow, viscoplastic relaxations near GBs. We aim to demonstrate our hypothesis through a 4-point bending apparatus that introduces a homogeneous in-plane stress state before indentation. We will perform simulations of a coarse-grained mesoscale theory of plasticity in a system with microhard GBs, benchmarked through 2.5D-DDP simulations, which will capture the effects of slow viscoplastic relaxations. Our nanoindentation abrupt event statistics studies, either in the micrograin interior or near GBs, will be correlated to direct observations of induced surface dislocation microstructures. For we will use TEM-based methods (PACOM/TKD) in combination with EBSD. Indents placed in the center of grains will be analyzed based on the orientation of the grain, while those near to GBs will be analyzed based on the orientation of each grain, the distance from the GB and the 5 degrees of freedom of the GB.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The Influence of Soil Properties on Sea-Breeze Circulations in the Southeast U.S.

Sea-breeze circulations (SBCs) are common weather phenomena at and near coastal regions. They form because of a thermal gradient between the land surface at the coast and the sea surface. In a mid-day regime, a “thermal low” generated at the warm coast will lead to rising air motion, creating a wind shift coming from the sea near the surface displacing the coastal air. A “return flow” moving back towards the sea is generated by upper-level divergence because of the rising motion from the thermal low. SBCs propagate and serve as a method of urban pollutant dispersion in the Los Angeles region of California and are constrained to the coast due to the topography of surrounding mountains serving as a boundary for further inland propagation. Within the northeast U.S. SBCs are seen in the warm season but tend to remain coastally bound due to Coriolis distortion over long distances. Within the southeast U.S. (SEUS), the paradigmatic example of SBCs occurs over the Florida peninsula, where thunderstorms form on a nearly daily occurrence due to the convergence of SBCs from the east and west sides of the peninsula. However, there are further examples of sea-breezes in the SEUS that warrant study. Within the region bordering the SEUS and the Mid-Atlantic, just east of the southern Appalachian Mountains, warm-season SBCs form at the coast of Georgia and the Carolinas. Relatively flat topography ~150-200km inland allows for mostly unimpeded inland SBC propagation. Through visual analysis, Viner et al. catalogued several SBCs that propagated as far inland as the Central Savannah River Area surrounding Augusta, Georgia. Wermter et al. found that while the land-sea thermal gradient at the coast can influence coastal SBC genesis, the inland-coastal thermal gradient over the land is the primary influencer on the speed and depth of inland propagation of SBCs in this region. Additionally, soil moisture itself is a known correlative factor to sea-breeze formation, as it influences the soil temperature and the thermal gradient needed for SBC genesis and inland penetration. Physick determined that higher latent heat fluxes associated with wetter soil dampen the land-sea thermal gradient and suppress the formation of a SBC. Physick determined that higher latent heat fluxes associated with wetter soil dampen the land-sea thermal gradient and suppress the formation of a SBC. Conversely, drier soil enhances the thermal gradient and promotes SBC formation. However, while there is an inverse relationship between soil moisture and SBCs, higher soil moisture can actually promote more convective rainfall following a SBC if it does not significantly impact the thermal gradient. While the relationship between soil moisture and SBC formation has been conceptually explored and modeled numerically, there is a research gap in observed connections. The Soil Moisture Active Passive (SMAP) satellite mission has been operational since 2015 and has been used to create high-resolution re-analytical Level 4 (L4) datasets of soil moisture and soil temperature at different soil depths: the surface (0-5cm) and rootzone (0-1m). The surface soil temperature effectively acts as the “skin temperature” of the surface at these levels, and a spatial map of the land-sea as well as the coastal-inland thermal gradients can be represented. SMAP data are also assimilated in some atmospheric models such at the High Resolution Rapid Refresh (HRRR) mesoscale model. We propose leveraging the use of SMAP products to fill in spatial gaps left by weather and mesonet stations within the SEUS region, as well as assessing the effectiveness of utilizing SMAP products towards SBC forecasting in both deterministic and machine learning (ML) models.

58 GEOSCIENCES↗