Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model initialization”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A Dynamic Risk Framework for the Physical Security of Nuclear Power Plants

This paper describes ongoing work within the Light Water Reactor Sustainability pathway at Idaho National Laboratory (INL) to optimize the security and cost of nuclear power plants. It introduces the dynamic risk assessment tool developed at INL, Event Modeling Risk Assessment using Linked Diagrams (EMRALD). EMRALD is leveraged to optimize the security posture of a nuclear power plant by integrating force-on-force (FOF) simulations and operator mitigation actions, including dynamic and flexible coping strategies (FLEX). To illustrate the methodology, four attack scenarios are modeled in a commercially available FOF simulation tool using a hypothetical nuclear power plant facility. The simulation results provide valuable insights into possible attack outcomes, as well as the probabilistic risk of a core damage event given these outcomes. Safety mitigation procedures are modeled in EMRALD dependent on the attack outcomes by considering human operator uncertainties. The results demonstrate that the number of armed responders can be optimized, while still maintaining the same protection level as the initial security posture. The proposed modeling and simulation framework of integrating FLEX equipment with FOF models enables the nuclear power plants to credit FLEX portable equipment in the plant security posture, resulting in an efficient and optimized physical security system.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Collaborative Research: Improved Efficiency and Coupling of the Radiation Code in the ACME Earth System Model (Final Report)

The complexity of radiative transfer, its importance to the exchange of energy in the climate system, and its high computational cost establishes importance of an accurate and efficient radiative transfer parameterization for climate simulation. Previous work by the proposing team at AER led to the development of the radiation code, RRTMG, which has been widely accepted by the global modeling community as a fast and accurate advancement over the previous generation of radiation codes. It has been in use in the NCAR CESM for many years, and it has been implemented in the initial version the DOE E3SM model. However, its computational cost remains high relative to other components in part due to its complexity and to its inefficient use of modern optimization strategies, and this project helped address this limitation for the code’s application in E3SM. Under other funding, the Investigators of this project led an effort to develop a high-performance broadband radiation code, called RTE+RRTMGP, which is a completely restructured code that will take advantage of modern computational capabilities to enhance its performance while retaining the strengths and accuracy of the original code. Designed to perform over a range of computer architectures, RTE+RRTMGP makes extensive use of Fortran 2003 features to improve both its efficiency and its use of memory. RTE+RRTMGP is expected to be adopted in the next generation of RTE+RRTMGP. This project allowed for advancements to the code’s computational capabilities and adding new and enhanced features. This included revising RTE+RRTMGP to run on GPU processors, analysis of and improvements to the code’s timing, modifying the gas optics in RRTMGP to increase the code’s accuracy, extensive validation of computed fluxes, heating rates and forcings, generating a cloud optical property and vertical sampling capabilities for RTE+RRTMGP, implementing a fast and accurate longwave scattering capability, and groundwork for the inclusion of a capability to specify solar variability. Many of the accomplishment in this project necessitated significant collaboration with the E3SM development team. The result of this project was optimization of a key physical component (radiative transfer calculations) of E3SM, directly supporting E3SM’s overarching global modeling objectives. More broadly, this project provided overall advancements in the use of radiative transfer calculations in atmospheric modeling and simulation, particularly for climate.

58 GEOSCIENCES↗
DOE OSTI2022

A Dynamic Risk Framework for the Optimization of Physical Security Posture of Nuclear Power Plants

This paper describes an ongoing work within the Light Water Reactor Sustainability pathway at Idaho National Laboratory (INL) to optimize security and cost of nuclear power plants. It introduces the dynamic risk assessment tool developed at INL, Event Modeling Risk Assessment using Linked Diagrams (EMRALD). EMRALD was leveraged to optimize the security posture of a nuclear power plant by integrating force-on-force (FOF) simulations and operator mitigation actions including the dynamic and flexible coping strategies (FLEX). To illustrate the methodology, four attack scenarios were modeled in a commercially available FOF simulation tool using a hypothetical nuclear power plant facility. The simulation results provide valuable insights into possible attack outcomes, as well as the probabilistic risk of core damage event given these outcomes. Safety mitigation procedures were modeled in EMRALD dependent on the attack outcomes by considering human operator uncertainties. The results demonstrate that the number of armed responders can be optimized, while still maintaining the same protection level as the initial security posture. The proposed modeling and simulation framework of integrating FLEX equipment with FOF models enables the nuclear power plants to credit FLEX portable equipment in the plant security posture, resulting in an efficient and optimized physical security system.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Legacy of the Madden-Julian Oscillation Working Group

There is a long history of observational and model evaluation of the MJO. Even so, the underlying physics of the MJO remains poorly understood. Initial analyses of the MJO in models concentrated on large-scale diagnosis using such variables as velocity potential and/or zonal winds, and occasionally precipitation or outgoing longwave radiation. The first time a large number of models were systematically evaluated using a common set of diagnostics was in Slingo et al. (1996), who evaluated 15 general circulation models that were contributed to the Atmospheric Model Intercomparison Project (AMIP; Gates et al. 1992). However, the high- LLNL-TR-762421 frequency (daily) model data available to Slingo et al. (1996) were limited, and it was only with the advent of NWP-based reanalyses that a more comprehensive understanding of the 3-dimensional structure and the thermodynamics of the MJO became possible, thus helping to establish the variables relevant to the study of the MJO. International interest in the MJO began to blossom on the heels of the 2006 “Workshop on the Organization and Maintenance of Tropical Convection and the Madden Julian Oscillation,” held at the International Center for Theoretical Physics in Trieste, Italy. With its stated goals of understanding the influence of MJO in the context of the tropical circulation and its remote impacts, and discussion of future directions including the MJO impact on predictive capability, the time was right for a more coordinated attack on understanding and simulating the MJO.

58 GEOSCIENCES↗
DOE OSTI2024

Simplified inelastic constitutive models for ASME Section III, Division 5 design by inelastic analysis

This report describes the development of simplified, universal constitutive model that captures the high temperature monotonic and cyclic behavior of a range of commonly-used high temperature materials. The goal of the work is to provide a simple, universal constitutive model to replace the current bespoke models for Grade 91, 316H, and Alloy 617 included in Nonmandatory Appendix HBB-Z of the ASME Boiler & Pressure Vessel Code, and to extend this model to cover Alloy 800H. We initiated this work in response to feedback from reactor vendors and other Code users requesting simplified models, compared to the current models, that are easier to implement and use in commercial finite element analysis software. This report describes the completion of this effort by developing a model to correct the defects in standard model forms presently used for high temperature material modeling, described in past work, developing and implementing new numerical methods to train this model against test data, and then actually training the model for the four materials. The report provides a complete mathematical description of the model along with the tabulated material coefficients for the four materials. The final step will be to formulate an ASME Code change to introduce the new models into the Code.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Validation of a Hybrid Domain Overlapping Coupling Between SAM and CFD Against the TALL-3D Transients

The System Thermal Hydraulics (STH) code SAM has been coupled to the Computational Fluid Dynamics (CFD) code Simcenter STAR-CCM+ utilizing a hybrid domain overlapping method with an explicit coupling in time. The coupling aims to extend the STH code’s applicability to scenarios where local momentum and energy transfers are important yet difficult for STH codes to capture, such as three-dimensional (3D) mixing. The coupling method’s numerical stability was verified in the past against two closed-loop configurations, and it was validated against a double T-junction experiment with 3D scalar mixing. In the present work, the coupling method is validated against the TALL-3D STH/CFD coupling benchmark facility. TALL-3D is a three-legged, liquid-metal facility with a large, pool-type enclosure (test section) that exhibits 3D flow effects to be modeled by a CFD code. The rest of the system exhibits approximately 1D behavior well-predicted by an STH code. First, the present STAR-CCM+ CFD model of the 3D test section is validated against experimental data. Then, the SAM-STARCCM+ coupled model is validated against six different TALL-3D steady states, including SAM standalone model results for comparison. Lastly, the SAM-STARCCM+ coupled model is validated against two TALL-3D transients, one exhibiting flow reversal in the test section and one exhibiting nonlinear, Limit Cycle Oscillations (LCO). For the first transient, the SAM-STARCCM+ coupled model properly predicts an increase in the test section’s inlet temperature during flow reversal, and this is not predicted by the SAM standalone model. Following flow reversal, the SAM-STARCCM+ coupled model better-predicts the initial flow recovery and following oscillations as the system approaches a final natural circulation state. For the second transient, no true final steady state is observed due to LCO. Neither the SAM-STARCCM+ coupled model nor the SAM standalone model can perfectly capture the experiment’s changing oscillation frequency during the transient. However, the SAM-STARCCM+ coupled model does reproduce the oscillatory feedback observed in the system. This is a significant achievement as the SAM-STARCCM+ coupled model only uses an explicit coupling in time, as opposed to a semi-implicit coupling. In comparison, previous STH/CFD coupling efforts of the TALL-3D facility required semi-implicit coupling to obtain similar results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Bulb-style Kaplan Turbine - Run-of-the-river (ror) Hydro Governor And Turbine Model

This software package includes hydro governor and turbine models developed in Simulink (Mathworks Inc.) and RSCAD (RTDS Technologies) for low-head bulb-style Kaplan turbine hydropower plants. The model developed in Simulink has been designed and tuned to match the governor-turbine response for a real-world hydropower unit owned by Idaho Falls Power (IFP). The Simulink and RSCAD models enable real-time testing in a hardware-in-the-loop (HIL) implementation using OPAL-RT and RTDS digital real-time simulators, respectively. An automatic initialization for dynamic simulation has also been integrated to each model in both Simulink and RSCAD. The RSCAD initialization involves a novel python-based interfacing for automatic modification of the initial conditions in the model. The parameters of these models can be tuned to match other hydropower plants in this class.

Alam, SMShafiul↗
DOE OSTI2023

Robust scalable initialization for Bayesian variational inference with multi-modal Laplace approximations

Predictive modeling typically relies on Bayesian model calibration to provide uncertainty quantification. Variational inference utilizing fully independent (“mean-field”) Gaussian distributions are often used as approximate probability density functions. This simplification is attractive since the number of variational parameters grows only linearly with the number of unknown model parameters. However, the resulting diagonal covariance structure and unimodal behavior can be too restrictive to provide useful approximations of intractable Bayesian posteriors that exhibit highly non-Gaussian behavior, including multimodality. High-fidelity surrogate posteriors for these problems can be obtained by considering the family of Gaussian mixtures. Gaussian mixtures are capable of capturing multiple modes and approximating any distribution to an arbitrary degree of accuracy, while maintaining some analytical tractability. Unfortunately, variational inference using Gaussian mixtures with full-covariance structures suffers from a quadratic growth in variational parameters with the number of model parameters. The existence of multiple local minima due to strong nonconvex trends in the loss functions often associated with variational inference present additional complications, These challenges motivate the need for robust initialization procedures to improve the performance and computational scalability of variational inference with mixture models. In this work, we propose a method for constructing an initial Gaussian mixture model approximation that can be used to warm-start the iterative solvers for variational inference. The procedure begins with a global optimization stage in model parameter space. In this step, local gradient-based optimization, globalized through multistart, is used to determine a set of local maxima, which we take to approximate the mixture component centers. Around each mode, a local Gaussian approximation is constructed via the Laplace approximation. Finally, the mixture weights are determined through constrained least squares regression. The robustness and scalability of the proposed methodology is demonstrated through application to an ensemble of synthetic tests using high-dimensional, multimodal probability density functions. Here, the practical aspects of the approach are demonstrated with inversion problems in structural dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Progress Report on Model Development for the Transport of Aerosol through Microchannels

This report summarizes the current progress in the development of a phenomenological model of aerosol transport, deposition, and plugging through microchannels. The purpose is to introduce a generic, reliable numerical model for the prediction of aerosol transport, deposition, and plugging in leak paths while accounting for potential plugging formation, to a user community involving researchers, regulators, and industry. In that regard, a Graphical User Interface (GUI) was generated by integrating the individual MATLAB scripts that make up the model and adding additional features to aid the general understanding of the problem to a user. This report focuses on the model development and features included. Further, predictions from the GUI are compared with experimental data for validation. The strength of the GUI is the ability of a user to plug in basic parameters such as the initial pressure conditions and canister model specifications in order to obtain a first principal approximation of vital information such as the blowdown pressure differential, aerosol penetration, and deposition as a function of time. The user can obtain this without the know-how of the underlying MATLAB scripts and thus enables the model to be readily applied by regulators, industry, and shareholders to reduce the uncertainty in the off-site radiological consequences. The report also lists future tasks under this work scope. This includes ongoing efforts to include additional crack geometries (divergent and divergent-convergent slots) to the model; GUI development and improvement based on feedback from users, and integration of aerosol source term data from Sibling Pin tests into the model as we approach realistic canister stress corrosion cracking–induced aerosol release scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Development of Saturated Zone Three-Dimensional Initial Condition Plumes for the Composite Analysis and Cumulative Impacts Evaluation Modeling

The purpose of this environmental calculation file (ECF) is to document and present the methodology, input data, and results of extending two-dimensional (2-D) plume depictions documented in DOE/RL- 2017-66, Hanford Site Groundwater Monitoring for 2017, into three dimensions (3-D) for select contaminants of interest (COIs) in the 200-BP-5, 200-UP-1, 200-ZP-1, and 200-PO-1 Groundwater Operable Units at the Hanford Site. Plumes with sufficient available vertical profile data were interpolated in 3-D using the Leapfrog Geo®1 geologic modeling platform, which includes an interpolation engine for analyzing and visualizing 3-D data. Results of these interpolations are intended for use in developing inputs for transport modeling to support forecasts of Central Plateau contaminant plume fate and transport.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Doubling protected land area may be inefficient at preserving the extent of undeveloped land and could cause substantial regional shifts in land use

Abstract Projection of land use and land‐cover change is highly uncertain yet drives critical estimates of carbon emissions, climate change, and food and bioenergy production. We use new, spatially explicit land availability data in conjunction with a model sensitivity analysis to estimate the effects of additional land protection on land use and land cover. The land availability data include protected land and agricultural suitability and is incorporated into the Moirai land data system for initializing the Global Change Analysis Model. Overall, decreasing land availability is relatively inefficient at preserving undeveloped land while having considerable regional land‐use impacts. Current amounts of protected area have little effect on land and crop production estimates, but including the spatial distribution of unsuitable (i.e., unavailable) land dramatically shifts bioenergy production from high northern latitudes to the rest of the world, compared with uniform availability. This highlights the importance of spatial heterogeneity in understanding and managing land change. Approximately doubling the current protected area to emulate a 30% protected area target may avoid land conversion by 2050 of less than half the newly protected extent while reducing bioenergy feedstock land by 10.4% and cropland and grazed pasture by over 3%. Regional bioenergy land may be reduced (increased) by up to 46% (36%), cropland reduced by up to 61%, pasture reduced by up to 100%, and harvested forest reduced by up to 35%. Only a few regions show notable gains in some undeveloped land types of up to 36%. Half of the regions can reach the target using only unsuitable land, which would minimize impacts on agriculture but may not meet conservation goals. Rather than focusing on an area target, a more robust approach may be to carefully select newly protected land to meet well‐defined conservation goals while minimizing impacts to agriculture.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Complex Fracture Nucleation and Evolution with Nonlocal Elastodynamics

A mechanical model is introduced for predicting the initiation and evolution of complex fracture patterns without the need for a damage variable or law. The model, a continuum variant of Newton’s second law, uses integral rather than partial differential operators where the region of integration is over finite domain. The force interaction is derived from a novel nonconvex strain energy density function, resulting in a nonmonotonic material model. The resulting equation of motion is proved to be mathematically well-posed. The model has the capacity to simulate nucleation and growth of multiple, mutually interacting dynamic fractures. In the limit of zero region of integration, the model reproduces the classic Griffith model of brittle fracture. As a result, the simplicity of the formulation avoids the need for supplemental kinetic relations that dictate crack growth or the need for an explicit damage evolution law.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Select Transient Metallic Fuel Bison Benchmarks

To explore the current capabilities of Bison with respect to transient metallic fuel modeling, four transient experiments are compared with their respective Bison model. The four experiments are the OPT-1 experiment of the Experimental Breeder Reactor II~(EBR-II) and the M5, M6, and M7 experiments of the Transient Reactor Test facility (TREAT). Each benchmark is a 2D-RZ geometry of a single pin employing coupled multiphysics. Bison models of the benchmark agree well and illustrate areas needing development. The areas of suggested development deal with better material models, better coolant channel module for transients, proper treatment of material melting, improved feedback with failure models, and an ability to set initial conditions from post irradiation examination (PIE) data. Two Fast Flux Test Facility (FFTF) pins are also modeled with M7 transient conditions to illustrate Bison being applied to investigate potential experiments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

A Bayesian approach to evaluation of soil biogeochemical models

Abstract. To make predictions about the carbon cycling consequences of rising global surface temperatures, Earth system scientists rely on mathematical soil biogeochemical models (SBMs). However, it is not clear which models have better predictive accuracy, and a rigorous quantitative approach for comparing and validating the predictions has yet to be established. In this study, we present a Bayesian approach to SBM comparison that can be incorporated into a statistical model selection framework. We compared the fits of linear and nonlinear SBMs to soil respiration data compiled in a recent meta-analysis of soil warming field experiments. Fit quality was quantified using Bayesian goodness-of-fit metrics, including the widely applicable information criterion (WAIC) and leave-one-out cross validation (LOO). We found that the linear model generally outperformed the nonlinear model at fitting the meta-analysis data set. Both WAIC and LOO computed higher overfitting risk and effective numbers of parameters for the nonlinear model compared to the linear model, conditional on the data set. Goodness of fit for both models generally improved when they were initialized with lower and more realistic steady-state soil organic carbon densities. Still, testing whether linear models offer definitively superior predictive performance over nonlinear models on a global scale will require comparisons with additional site-specific data sets of suitable size and dimensionality. Such comparisons can build upon the approach defined in this study to make more rigorous statistical determinations about model accuracy while leveraging emerging data sets, such as those from long-term ecological research experiments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Comparative Analysis of HEATNETS for Geothermal Network Performance: Preprint

Thermal energy networks (TENs), also known as 5th generation district energy systems, or more specifically geothermal networks when exchanging heat with geothermal boreholes, are an important technology for decarbonization. In these networks an ambient loop connects buildings and thermal sources, such as a borehole field, to exchange energy and maintain a desired loop temperature. Water-source heat pumps are used at the buildings to connect to the ambient or thermal loop to meet to the building heating and cooling loads and maintain comfort. A semi-transient, reduced-order technical model and techno-economic model, called HEATNETS, has been developed at NREL that captures the flow of energy around a TEN. In this work, a comparison of the HEATNETS technical model and a well-known coding platform used for modeling geothermal networks, TRNSYS, has been completed for a proposed geothermal network as a verification process. Hourly data provided from the TRNSYS simulation included building loads, pumping power, heat pump power, temperature entering and leaving the borehole field, and mass flow rates. The hourly borehole temperatures were used to create a linear regression model utilized in HEATNETS to estimate the borehole field heat exchange. The building loads and mass flow rates were direct inputs to HEATNETS while the pumping power, heat pump power, borehole temperatures, and coefficients of performance were all simulated and calculated by HEATNETS, allowing for direct comparison of the thermal energy transfer, rather than also comparing control systems responses. HEATNETS considers the full process from design inputs to economic outputs and can provide modeling options for high-level initial system design and operational optimization. This study shows that HEATNETS, while not intended to replace other modeling tools, can be a unique modeling tool for the performance of a full geothermal network system.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Modeling creep behavior in ceramic matrix composites

Here, in this work, a three-dimensional viscoplasticity formulation with progressive damage is developed and used to investigate the complex time-dependent constituent load transfer and progressive damage behavior in ceramic matrix composites (CMCs) subjected to creep. The viscoplasticity formulation is based on Hill's orthotropic plastic potential, an associative flow rule, and the Norton-Bailey creep power law with Arrhenius temperature dependence. A fracture mechanics-informed isotropic matrix damage model is used to account for CMC brittle matrix damage initiation and propagation, in which two scalar damage variables capture the effects of matrix porosity as well as matrix property degradation due to matrix crack initiation and propagation. The Curtin progressive fiber damage model is utilized to simulate progressive fiber failure. The creep-damage formulation is subsequently implemented as a constitutive model in the generalized method of cells (GMC) micromechanics formulation to simulate time-dependent deformation and material damage under creep loading conditions. The developed framework is used to simulate creep of single fiber SiC/SiC microcomposites. Simulation results are in excellent agreement with experimental and numerical data available in the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Impact of Lateral Flow on Surface Water and Energy Budgets over the Southern Great Plains – A Modeling Study

As we see the horizontal grid spacing decrease, treatment of hydrologic processes in land surface models, such as the lateral flow of surface and subsurface flow, need to be explicitly represented. Unlike previous studies that mainly focused on the mountainous regions, in this study the offline WRF-Hydro model is employed to study the impact of lateral flow on soil moisture and energy fluxes over the relatively flat southern Great Plains (SGP). The vast amount of measurements over the SGP provide an unique opportunity to assess the model behavior. In addition, newly developed land surface properties and input forcing are ingested into the model, in an attempt to reduce uncertainties associated with the initial and boundary forcing and help to identify model deficiencies. Our results show that the more realistic inputs (parameters, soil types, forcing) lead to larger underestimation of latent heat flux and dry bias, indicating the existence of model structural uncertainty (embedded errors) in WRF-Hydro that need to be characterized to inform future model development efforts. Including lateral flow processes partly mitigates the model deficiencies in representing hydrologic processes and alleviates the dry bias. In particular, both surface and subsurface lateral flow increase soil moisture mainly over the lower elevations, except that subsurface flow also affects soil moisture over steeper terrains. Additional simulations are performed to assess the effect of routing resolution on model results. When LSM resolution is high, noticeable differences in soil moisture are produced between different routing resolutions especially over steep terrain. Whereas when LSM resolution is coarse, differences between routing resolutions become negligible, especially over flat terrain.

54 ENVIRONMENTAL SCIENCES↗