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 1,063 records · Page 59

Advanced modeling and simulation of research reactors using dynamic mode decomposition

Full text of publication follows. Due to the ever-increasing safety requirements, the current trend of nuclear reactor analysis is shifting towards high-fidelity multi-physics models, which have a very high computational cost and modelling complexity. As the cost of even a single model run makes it impossible to analyse the behaviour and performance of these models on large-scale commercial plants, it has become even more significant to provide suitable benchmarks to validate and test them extensively. In this sense, research reactors offer a promising solution for the initial validation of high-fidelity models, as they are significantly smaller than commercial reactors and their characteristics are well known. In particular, the reactors of the TRIGA family have been used to assess and validate models and methods for Generation-IV designs, as they have some similar features (such as the dominance of natural convection as cooling mechanism and the difficulties in performing sub-channel analysis using standard codes). Still, the computational requirements of high-fidelity models make them unsuitable for real-time analysis, even following their assessment on research reactors. In this sense, Model Order Reduction (MOR) techniques give an additional strategy to reduce the computational cost of high-fidelity models (whilst preserving sufficient accuracy). In particular, this work focuses on Dynamic Mode Decomposition (DMD), a non-intrusive MOR technique that aims at representing models with explicit temporal dynamics by extracting the time-varying characteristics and the governing structures based only on a set of available data, thus without needing any underlying knowledge of the governing equations. In addition, DMD also computes a low-dimensional surrogate of the dynamic matrix of the system, making it suited for stability analysis and real-time evaluations. This work focuses on the application and validation of the DMD method on the Computational Fluid-Dynamics (CFD) model TRIGA Mark II reactor, also discussing in detail the potentiality of this algorithm as an advanced modelling tool for nuclear reactor analysis. (author)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Uncertainty Propagation in Coupled Atmosphere–Wave–Ocean Prediction System: A Study of Hurricane Earl (2010)

This study focuses on understanding the evolution of Hurricane Earl (2010) with respect to random perturbations in the storm’s initial strength, size, and asymmetry in wind distribution. We rely on the Unified Wave Interface-Coupled Model (UWIN-CM), a fully coupled atmosphere–wave–ocean system to generate a storm realization ensemble, and use polynomial chaos (PC) expansions to build surrogate models for time evolution of both the maximum wind speed and minimum sea level pressure in Earl. The resulting PC surrogate models provide statistical insights on probability distributions of model responses throughout the simulation time span. Statistical analysis of rapid intensification (RI) suggests that initial perturbations having intensified and counterclockwise-rotated winds are more likely to undergo RI. In addition, for the range of initial conditions considered RI seems mostly sensitive to azimuthally averaged maximum wind speed and asymmetry orientation, rather than storm size and asymmetry magnitude; this is consistent with global sensitivity analysis of PC surrogate models. Finally, we combine initial condition perturbations with a stochastic kinetic energy backscatter scheme (SKEBS) forcing in the UWIN-CM simulations and conclude that the storm tracks are substantially influenced by the SKEBS forcing perturbations, whereas the perturbations in initial conditions alone had only limited impact on the storm-track forecast.

54 ENVIRONMENTAL SCIENCES↗

Automated model-predictive design of synthetic promoters to control transcriptional profiles in bacteria

Abstract Transcription rates are regulated by the interactions between RNA polymerase, sigma factor, and promoter DNA sequences in bacteria. However, it remains unclear how non-canonical sequence motifs collectively control transcription rates. Here, we combine massively parallel assays, biophysics, and machine learning to develop a 346-parameter model that predicts site-specific transcription initiation rates for any σ 70 promoter sequence, validated across 22132 bacterial promoters with diverse sequences. We apply the model to predict genetic context effects, design σ 70 promoters with desired transcription rates, and identify undesired promoters inside engineered genetic systems. The model provides a biophysical basis for understanding gene regulation in natural genetic systems and precise transcriptional control for engineering synthetic genetic systems.

42 ENGINEERING↗

Approximate Green’s Function Coupled Cluster A Method Employing Effective Dimension Reduction

The Green’s function coupled cluster (GFCC) method, originally proposed in the early 1990s, is a powerful many-body tool for computing and analyzing the electronic structures of molecular structures and periodic systems, especially when electrons are strongly correlated. However, in order for the GFCC method to become an approach that can be routinely used in the electronic structure calculations, more robust linear algebraic techniques and approximations need to be employed to reduce its extremely high computational overhead. In our recent studies, we showed that the GFCC equation can be numerically solved by seeking solutions of multiple linear systems in the frequency regime of interest, which can be easily distributed in a massively parallel environment, and broaden the avenue for more state-of-the-art numerical linear algebraic tools associated with multiple linear systems to be kicked in to this field. In the present pursuit of this practice, we show a successful application of model order reduction (MOR) techniques in the GFCC framework. Briefly speaking, for a frequency regime of interest that requires a high resolution of spectral function, instead of solving GFCC linear equation in the original dimension for every single frequency point, an easily-solvable linear system model with reduced dimension can be built upon projecting the original GFCC linear system into a subspace, from which the interpolation and extrapolation can give a reasonable approximation to the results of the original GFCC linear equations. Here, we show that the subspace can be properly constructed in an iterative manner from the auxiliary vectors of the GFCC linear equations at some select frequencies inside the frequency regime. During the iterations, the quality of the subspace, as well as the linear system model, can be systematically improved. The approach is tested in this work in terms of the efficiency and accuracy of computing spectral functions for some typical molecular systems including carbon monoxide, 1,3-butadiene, benzene, and adenine molecules. As a byproduct, the high quality linear system model is also found to be able to provide high quality initial guess for the existing linear solver to improve the convergence performance.

Peng, Bo↗

Physics beyond colliders at CERN: beyond the Standard Model working group report

The Physics Beyond Colliders initiative is an exploratory study aimed at exploiting the full scientific potential of the CERN’s accelerator complex and scientific infrastructures through projects complementary to the LHC and other possible future colliders. These projects will target fundamental physics questions in modern particle physics. This document presents the status of the proposals presented in the framework of the Beyond the Standard Model physics working group, and explore their physics reach and the impact that CERN could have in the next 10-20 years on the international landscape.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

FAST-1.0.1 User Installation and Verification Guide

The purpose of this document is to provide the user information about the installation of FAST-1.0.1 on their computers or servers. General information about the code and supported operating systems is described in Section 1.0.1. Self-service oriented FAST-1.0.1 software licensing steps are described in Section 2.0. An installation verification test suite is provided with FAST-1.0.1 and described in Section 3.0. A convenience script for converting FRAPCON to FAST inputs is discussed in Section 4.0. FAST-1.0.1 was developed and released under a software quality assurance program based upon NQA-1-2017. FAST-1.0.1 is the latest baseline code and the result of a bug fixes to FAST-1.0 with other software and methodology developments. The installation verification test suite contains both steady state and transient Anticipated Operation Occurrences (AOOs). Capability to model accident conditions, such as Reactivity Initiated Accidents (RIAs) and Loss Of Coolant Accidents (LOCAs), are targeted for a later release of FAST. FRAPTRAN-2.0 will continue to be used for accident conditions until the release of FAST with accident condition modeling capabilities.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Convective-Permitting Hindcast Simulations during the North American Monsoon GPS Transect Experiment 2013: Establishing Baseline Model Performance without Data Assimilation

During the North American monsoon global positioning system (GPS) Transect Experiment 2013, daily convective-permitting WRF simulations are performed in northwestern Mexico and the southern Arizona border region using the operational Global Forecast System (GFS) and North American Mesoscale Forecast System (NAM) models as lateral boundary forcing and initial conditions. Compared to GPS precipitable water vapor (PWV), the WRF simulations display a consistent moist bias in the initial specification of PWV leading to convection beginning 3–6 h early. Given appreciable observed rainfall, days are classified as strongly and weakly forced based only on the presence of an inverted trough (IV); gulf surges did not noticeably impact the development of mesoscale convective systems (MCSs) and related convection in northwestern Mexico. Strongly forced days display higher modeled precipitation forecast skill than weakly forced days in the slopes of the northern Sierra Madre Occidental (SMO) away from the crest, especially toward the west where MCSs account for the greatest proportion of all monsoon-related precipitation. A case study spanning 8–10 July 2013 illustrates two consecutive days when nearly identical MCSs evolved over northern Sonora. Although a salient MCS is simulated on the strongly forced day (9–10 July 2013) when an IV is approaching the core monsoon region, a simulated MCS is basically nonexistent on the weakly forced day (8–9 July 2013) when the IV is farther away. Finally, the greater sensitivity to the initial specification of PWV in the weakly forced day suggests that assimilation of GPS-derived PWV for these types of days may be of greatest value in improving model precipitation forecasts.

54 ENVIRONMENTAL SCIENCES↗

Equilibration and Diffusion for a Dynamical Lorentz Gas

We consider a model of a dynamical Lorentz gaz: a single particle is moving in ℝ{sup d} through an array of fixed and soft scatterers each possessing an internal degree of freedom coupled to the particle. Assuming the initial velocity is sufficiently high and modelling the parameters of the scatterers as random variables, we describe the evolution of the kinetic energy of the particle by a Markov chain for which each step corresponds to a collision. We show that the momentum distribution of the particle approaches a Maxwell–Boltzmann distribution with effective temperature T such that k{sub B}T corresponds to an average of the scatterers’ kinetic energy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Using an Advanced Distribution Management System Test Bed to Evaluate the Impact of Model Quality on Volt/VAR Optimization: Preprint

In this paper, we present a test bed for evaluating existing and future advanced distribution management system (ADMS) applications in a realistic laboratory setting, including other utility management systems and field equipment. We present an example of using it to evaluate the impact of the ADMS network model quality on a Volt/VAR optimization (VVO) application. The test bed integrates a commercial ADMS with a real-time simulation model of a utility distribution feeder. Representative power and controller hardware are integrated through hardware-in-the-loop (HIL) techniques. The performance of the ADMS VVO application is also evaluated for different levels of measurement density. Initial results indicate that a higher model quality achieves the highest possible energy savings while avoiding voltage violations, whereas a lower model quality results in increased energy savings but at the expense of more voltage violations.

ADMS↗

A helical actuator driven by biased SMA: design, model, and experiment

A helical actuator driven by biased shape memory alloy (SMA) patterns embedded into a soft composite ribbon base is presented in this work. Instead of common U-shape SMA wires, a single SMA wire is woven into planar patterns, which enable helical deformation of the composite ribbon from an initially flat geometry. An analytical static model is established for accurate and rapid prediction of the helical reconfiguration arising from the shape memory effect of woven SMA patterns, followed by validation of the static model using the finite element method (FEM). The finite element results are compared with the analytical solutions given by this static model, which show a high agreement. Parametric study of the influences of eight independent design variables on the dependent helical parameters, such as combined curvatures, pitches, and helical angles, is completed. It is found that the helically deformed geometry is mainly dominated by diameters, biased positions, inclined angles, and numbers of skewed segments of the SMA wire. Fabrication and in-situ experimental test of a prototype of such helical actuators qualitatively demonstrate its dramatic three-dimensional (3D) spiral reconfiguration from a two-dimensional (2D) flat ribbon. In conclusion, such SMA patterns will allow more diverse designs of soft actuators for a wider range of robotic applications.

36 MATERIALS SCIENCE↗

Energy based fracture initiation criterion for strain-crystallizing rubber-like materials with pre-existing cracks

Fracture prediction is indispensable for polymers, like rubbers, which have a broad range of applications mainly due to their high extensibility. The phenomenon known as strain-induced crystallization further contributes to the fracture toughness of certain rubbers. Here in this study, a criterion based on internal bond energy, incorporating the effects of crystallization, is proposed to predict fracture initiation in rubber-like materials with pre-existing cracks. First, a multi-scale mechanical model is developed for characterizing the behavior of rubber when subjected to both uniaxial and biaxial deformation states. At the microscale, both the amorphous and crystalline chain segments are modeled as elastic in order to consider the energy contribution by the molecular bond distortions. This internal energy is considered along with the entropic and crystalline free energy for each chain. In the chain model, the effects of loading condition and the relative orientation of a chain on its crystallinity are taken into account. At the macroscopic scale, an existing crystallinity distribution function is adapted and a mixed finite element formulation with an augmented Lagrangian multiplier is utilized to impose the incompressibility constraint. A non-affine maximal advance path constraint based homogenization model is utilized for bridging the two scales. Its potential to account for anisotropy in the stretched network compels the model to be preferable due to its physical significance, for the purpose of fracture modeling. The rigidity of the crystallites is accounted for by proposing a crystallite distortion energy in addition to the critical bond dissociation energy, for fracture initiation to occur. The model is validated by comparison with existing experimental results for both crystallizing and non-crystallizing rubbers. In addition to its potential to predict the material behavior when subjected to uniaxial and biaxial loading, the capability of the model to quantitatively estimate the effect of crystallization on fracture initiation is also verified.

36 MATERIALS SCIENCE↗

First high-resolution γ-ray spectroscopy of 41 Si

The first high-resolution in-beam γ -ray spectroscopy is reported for the neutron-rich nucleus 41 Si, a tenant of the N = 28 island of inversion. Excited states were populated in the direct one-proton removal reaction from 42 P projectiles and pn removal from 43 P. Seven γ-ray transitions were observed, only one of which had been reported previously in the literature. Furthermore, this makes 41 Si the most neutron-rich odd-even N = 27 isotone with high-resolution excited-state information. For the one-proton removal, the measured partial cross-section distribution to all observed bound final states is contrasted with results from direct one-proton removal calculations that combine eikonal reaction dynamics with SDPF-MU shell-model spectroscopic factors and assume various possible initial states for the poorly known 42 P projectile. Rather distinct calculated cross-section distributions emerge that, in comparison to the new data, imply that the initial state in 42 P is most likely 3 – or 2 – rather than 1 – or 0 – , the predicted shell-model ground state of 42 P. It is further shown that the level scheme from the novel VS-IMSRG calculation closely agrees with the one of SDPF-MU, the most successful phenomenological shell-model effective interaction in describing the much discussed neighboring isotope 42 Si, perhaps cross-validating these complementary approaches on the quest to model rapid shell evolution away from the valley of β stability.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Process Optimization and Modeling for Minerals Sustainability (PrOMMiS) v1.0.0

The U.S. Department of Energy's ("DOE") Process Optimization and Modeling for Minerals Sustainability project ("PrOMMiS Project") was initiated in 2023 to transform the national Critical Minerals and Rare Earth Elements (CM & REE) landscape, thereby supporting DOE's three enduring strategic objectives: security, economic competitiveness, and environmental responsibility. To address this initiative, the PrOMMiS Project will develop, demonstrate, and deploy an open-source, advanced system modeling, optimization, and analysis application ("PrOMMiS Software Application") that will support industry decision-makers by accelerating the scale-up of novel CM & REE technologies by de-risking the development and deployment of commercial scale processes and maximizing learning throughout the development cycle.

Beattie, KeithS [Lawrence Berkeley National Labora↗

Renewable Energy Potential Model: Geothermal Supply Curves

The Renewable Energy Potential (reV) model is a geospatial platform for estimating technical potential and developing renewable energy supply curves, initially developed for wind and solar technologies. The model evaluates deployment constraints, considering land use, environmental, and cultural factors, and estimates the distance to existing grid features to connect future plants (Maclaurin et al., 2021). A pressing deficiency in the reV model, however, is representation of geothermal electricity generation technologies. To address this gap, we developed a novel geothermal generation module for reV that allows for representation and analysis at the same level of detail as other renewable technologies. The included paper describes our process for evaluating data sources for the modeling, and presents five preliminary reV geothermal results. More specifically, we present two sets of resource data that represent upper and lower bounds for geothermal potential. We then present several sensitivity runs using the upper bound resource data; the results are encouraging that levelized cost of electricity (LCOE) can be reduced by optimizing the location and estimated capacity of the spatially diverse geothermal resource while considering the distance to existing grid infrastructure. Our preliminary supply curves and levelized cost of electricity (LCOE) results provided here should be considered with care due to the high uncertainty in geothermal resource potential data. We present median LCOE values for the conterminous U.S. for three scenarios: two hydrothermal (3.5km depth, USGS heat flow & SMU temperatures respectively) and one EGS (4.5km depth, SMU temperatures). The capital and operating costs for each respective technology are modeled. We also compare results using two different resource data sources.

15 GEOTHERMAL ENERGY↗

Spatio-temporal Estimates of Disease Transmission Parameters for COVID-19 with a Fully-Coupled, County-Level Model of the United States

Sandia National Laboratories has developed a capability to estimate parameters of epidemiological models from case reporting data to support responses to the COVID-19 pandemic. A differentiating feature of this work is the ability to simultaneously estimate county-specific disease transmission parameters in a nation-wide model that considers mobility between counties. The approach is focused on estimating parameters in a stochastic SEIR model that considers mobility between model patches (i.e., counties) as well as additional infectious compartments. The inference engine developed by Sandia includes (1) reconstruction and (2) transmission parameter inference. Reconstruction involves estimating current population counts within each of the compartments in a modified SEIR model from reported case data. Reconstruction produces input for the inference formulations, and it provides initial conditions that can be used in other modeling and planning efforts. Inference involves the solution of a large-scale optimization problem to estimate the time profiles for the transmission parameters in each county. These provide quantification of changes in the transmission parameter over time (e.g., due to impact of intervention strategies). This capability has been implemented in a Python-based software package, epi_inference, that makes extensive use of Pyomo [5] and IPOPT [10] to formulate and solve the inference formulations.

42 ENGINEERING↗

Development of a Geothermal Module in reV: Quantifying the Geothermal Potential while Accounting for the Geospatial Intersection of the Grid Infrastructure and Land Use Characteristics

The Renewable Energy Potential (reV) model is a geospatial platform for estimating technical potential and developing renewable energy supply curves, initially developed for wind and solar technologies. The model evaluates deployment constraints, considering land use, environmental, and cultural factors, and estimates the distance to existing grid features to connect future plants (Maclaurin et al., 2021). A pressing deficiency in the reV model, however, is representation of geothermal electricity generation technologies. To address this gap, we developed a novel geothermal generation module for reV that allows for representation and analysis at the same level of detail as other renewable technologies. This paper describes our process for evaluating data sources for the modeling, and presents five preliminary reV geothermal results. More specifically, we present two sets of resource data that represent upper and lower bounds for geothermal potential. We then present several sensitivity runs using the upper bound resource data; the results are encouraging that levelized cost of electricity (LCOE) can be reduced by optimizing the location and estimated capacity of the spatially diverse geothermal resource while considering the distance to existing grid infrastructure. Our preliminary supply curves and levelized cost of electricity (LCOE) results should be considered with care due to the highly uncertainty in geothermal resource potential data. We present median LCOE values for the conterminous U.S. for five scenarios: four hydrothermal (3.5km depth) and one EGS (4.5km depth). The capital and operating costs for each respective technology are modeled. We also compare results using two different resource data sources.

exclusions↗

Development of a Geothermal Module in reV: Quantifying the Geothermal Potential While Accounting for the Geospatial Intersection of the Grid Infrastructure and Land Use Characteristics: Preprint

The Renewable Energy Potential (reV) model is a geospatial platform for estimating technical potential and developing renewable energy supply curves, initially developed for wind and solar technologies. The model evaluates deployment constraints, considering land use, environmental, and cultural factors, and estimates the distance to existing grid features to connect future plants (Maclaurin et al., 2021). A pressing deficiency in the reV model, however, is representation of geothermal electricity generation technologies. To address this gap, we developed a novel geothermal generation module for reV that allows for representation and analysis at the same level of detail as other renewable technologies. This paper describes our process for evaluating data sources for the modeling, and presents five preliminary reV geothermal results. More specifically, we present two sets of resource data that represent upper and lower bounds for geothermal potential. We then present several sensitivity runs using the upper bound resource data; the results are encouraging that levelized cost of electricity (LCOE) can be reduced by optimizing the location and estimated capacity of the spatially diverse geothermal resource while considering the distance to existing grid infrastructure. Our preliminary supply curves and levelized cost of electricity (LCOE) results should be considered with care due to the highly uncertainty in geothermal resource potential data. We present median LCOE values for the conterminous U.S. for five scenarios: four hydrothermal (3.5km depth) and one EGS (4.5km depth). The capital and operating costs for each respective technology are modeled. We also compare results using two different resource data sources.

exclusions↗