Search NASA⌕ Search

SEARCH · Search NASA

Results for “reference model”

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,081 records · Page 60

Discovering the nature of neutrinos: How nuclear theory will advance this grand endeavor. Final report

Why is there more matter than antimatter in the visible universe? Are neutrinos their own antiparticles? These are some of the most profound questions in science. The detection of a (0vββ) neutrinoless doublebeta decay and the measurement of its half-life is one of the most important experimental endeavors in science because it will answer these questions. Paramount to these measurements are accurate predictions with quantified uncertainties of the 0vββ-decay nuclear transitions, especially for the germanium-76 and xeon-136 nuclei, two candidate materials for building the next generation of DOE funded 0vββ-decay experiments. These predictions will guide the Office of Science’s down selection of experiments as well as aid in determining how much material is needed by these experiments. We studied a specific difficult to assess source of uncertainty that will impact coming microscopic predictions, uncertainty which arises from an incomplete yet systematically improvable model of interactions between nucleons and of nucleons with the other particles involved in decay process, referred to as the theory uncertainty. We employed a surrogate model of the germanium-76 nucleus, and its decay products, just detailed enough to adequately mimic the decaying nucleons need for a reasonable measure of this uncertainty and nothing more, avoiding the immense computational expense required in quantifying uncertainties of the full calculation. We found that microscopic calculations of this process can be expected to achieve relative theory uncertainties around 30%, including an unconstrained contribution that had the potential to spoil microscopic calculations. This results in a nearly order of magnitude on previously estimated uncertainties for this rare decay process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

HostSub_GP: Precise Galaxy Background Subtraction in Transient Long-slit Spectroscopy with Gaussian Processes

We present a novel host galaxy subtraction technique in long-slit spectroscopy for extragalactic transients. Unlike classic methods which generally estimate the background using simple interpolation of local galaxy flux in the 2D spectrum, our approach leverages multi-band archival images of the host galaxies to model the background emission from the galaxy in the 2D spectrum. Such imaging encodes the wavelength-dependent galaxy profile along the slit, and is readily accessible through wide-field imaging surveys. We construct a smooth prior for the 2D galaxy profile with a Gaussian process (GP) based on these reference images, and use another GP to model the correlated deviations from the prior in the observed spectrum. This enables accurate inference of the galaxy flux blended with the transient. On synthetic long-slit data of a spiral galaxy extracted from a Multi Unit Spectroscopic Explorer hyper-spectral cube, the GP method remains robust as long as the host galaxy is spatially resolved and consistently outperforms classic methods. We apply the method to archival Keck spectra of two real transients, SN 2019eix and AT 2019qiz, to further demonstrate how the method uniquely recovers weak spectral features amid strong galaxy contamination, enabling refined constraints on the properties of both transients. We have released the software implementation, HostSub_GP, a scalable toolkit that leverages JAX, with an MIT license.

79 ASTRONOMY AND ASTROPHYSICS↗

Benchmark of numerical modeling approaches on the systematic performance evaluation of wave energy converters

Different numerical modeling methods have been developed and applied to evaluate a variety of performance indicators of wave energy converters (WECs), including the power performance, structural loads, levelized cost of energy, etc. Based on the modeling fidelity, the commonly used numerical modeling approaches can be classified as linear modeling, weakly nonlinear modeling and fully nonlinear modeling approaches. Each method differs in accuracy and computational efficiency, making them suitable for different stages of WEC design. However, the selection of modeling approach could significantly impact evaluation outcomes. For instance, simplified linear models may underestimate structural loads or overestimate energy production in some operational conditions, potentially leading to less cost-effective designs. Given the widespread utilization of these models, it is essential to understand the uncertainties brought by them in performance evaluations. This work is dedicated to benchmarking different linear-potential-flow-based numerical models for evaluating the systematic performance of WECs. Three representative numerical modeling approaches are considered in this work, including linear frequency-domain modeling, statistically linearized spectral-domain modeling and Cummins equation-based nonlinear time-domain modeling. A generic point absorber WEC is considered as the research reference in this work, and different sea sites are taken into account. The numerical models are utilized to predict critical performance indicators, including power performance, the annual energy production, the capacity factor, the levelized cost of energy and the PTO fatigue loads. By comparing the results, this work identifies the uncertainties associated with different modeling approaches in evaluating WEC performance.

Fatigue↗

Model of Solar Wind in the Heliosphere at Low and High Latitudes

The spatial distributions of the magnetic field, plasma density, and current at distances of (20–400)R{sub S} from the Sun (where R{sub S} is the solar radius) are investigated within a stationary axisymmetric MHD model of the solar wind (SW) at all latitudes in the inertial frame of reference with the origin at the center of the Sun. The model takes into account differential (with respect to the heliolatitude) rotation of the Sun and full corotation of plasma inside a boundary sphere of radius 20RS, which breaks down beyond this sphere. Self-consistent distributions of the plasma density, current, and magnetic field in the SW are obtained by numerically solving a set of time-independent MHD equations in spherical coordinates. It is demonstrated that the calculated results do not contradict observational data and describe a gradual transition from the fast SW at high heliolatitudes to the slow SW at low heliolatitudes, as well as the steepening of the profiles of the main SW characteristics with increasing distance from the Sun. The obtained dependences extend understanding of the SW structure at low and high latitudes and agree with the well-known Parker model in the limit of a small Ampère force.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

THE BENEFITS OF VLBI ASTROMETRY TO PULSAR TIMING ARRAY SEARCHES FOR GRAVITATIONAL RADIATION

Precision astrometry is an integral component of successful pulsar timing campaigns. Astrometric parameters are commonly derived by fitting them as parameters of a timing model to a series of pulse times of arrival (TOAs). TOAs measured to microsecond precision over spans of several years can yield position measurements with sub-milliarcsecond precision. However, timing-based astrometry can become biased if a pulsar displays any red spin noise or a red signal produced by the stochastic gravitational wave background. We investigate how noise of different spectral types is absorbed by timing models, leading to significant estimation biases in the astrometric parameters. We find that commonly used techniques for fitting timing models in the presence of red noise (Cholesky whitening) prevent the absorption of noise into the timing model remarkably well if the time baseline of observations exceeds several years, but are inadequate for dealing with shorter pulsar data sets. Independent of timing, pulsar-optimized very long baseline interferometry (VLBI) is capable of providing position estimates precise to the sub-milliarcsecond levels needed for high-precision timing. In order to make VLBI astrometric parameters useful in pulsar timing models, the transformation between the International Celestial Reference Frame (ICRF) and the dynamical solar system ephemeris used for pulsar timing must be constrained to within a few microarcseconds. We compute a transformation between the ICRF and pulsar timing frames and quantitatively discuss how the transformation will improve in coming years. We find that incorporating VLBI astrometry into the timing models of pulsars for which only a couple of years of timing data exist will lead to more realistic assessments of red spin noise and could enhance the amplitude of gravitational wave signatures in post-fit timing residuals by factors of 20 or more.

79 ASTRONOMY AND ASTROPHYSICS↗

Discovering a Transferable Charge Assignment Model Using Machine Learning

Partial atomic charge assignment is of immense practical value to force field parametrization, molecular docking, and cheminformatics. Machine learning has emerged as a powerful tool for modeling chemistry at unprecedented computational speeds given accurate reference data. However, certain tasks, such as charge assignment, do not have a unique solution. We use a machine learning algorithm to discover a new charge assignment model by learning to replicate molecular dipole moments across a large, diverse set of nonequilibrium conformations of molecules containing C, H, N, and O atoms. The new model, called Affordable Charge Assignment (ACA), is computationally inexpensive and predicts dipoles of out-of-sample molecules accurately. Furthermore, dipole-inferred ACA charges are transferable to dipole and even quadrupole moments of much larger molecules than those used for training. We apply ACA to dynamical trajectories of biomolecules and produce their infrared spectra. Additionally, we find that ACA assigns similar charges to Charge Model 5 but with greatly reduced computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

ISMabs: A COMPREHENSIVE X-RAY ABSORPTION MODEL FOR THE INTERSTELLAR MEDIUM

We present an X-ray absorption model for the interstellar medium, to be referred to as ISMabs, that takes into account both neutral and ionized species of cosmically abundant elements, and includes the most accurate atomic data available. Using high-resolution spectra from eight X-ray binaries obtained with the Chandra High Energy Transmission Grating Spectrometer, we proceed to benchmark the atomic data in the model particularly in the neon K-edge region. Compared with previous photoabsorption models, which solely rely on neutral species, the inclusion of ions leads to improved spectral fits. Fit parameters comprise the column densities of abundant contributors that allow direct estimates of the ionization states. ISMabs is provided in the appropriate format to be implemented in widely used X-ray spectral fitting packages such as XSPEC, ISIS, and SHERPA.

79 ASTRONOMY AND ASTROPHYSICS↗

Adaptive Learning for Concept Drift in Application Performance Modeling

Supervised learning is a promising approach for modeling the performance of applications running on large HPC systems. A key assumption in supervised learning is that the training and testing data are obtained under the same conditions. However, in production HPC systems these conditions might not hold because the conditions of the platform can change over time as a result of hardware degradation, hardware replacement, software upgrade, and configuration updates. These changes could alter the data distribution in a way that affects the accuracy of the predictive performance models and render them less useful; this phenomenon is referred to as concept drift. Ignoring concept drift can lead to suboptimal resource usage and decreased efficiency when those performance models are deployed for tuning and job scheduling in production systems. To address this issue, we propose a concept-drift-aware predictive modeling approach that comprises two components: (1) an online Bayesian changepoint detection method that can automatically identify the location of events that lead to concept drift in near-real time and (2) a moment-matching transformation inspired by transfer learning that converts the training data collected before the drift to be useful for retraining. We use application input/output performance data collected on Cori, a production supercomputing system at the National Energy Research Scientific Computing Center, to demonstrate the effectiveness of our approach. The results show that concept-drift-aware models obtain significant improvement in accuracy; the median absolute error of the best-performing Gaussian process regression improved by 58.8% when the proposed approaches were used.

adaptive learning↗

Finite deformation implementation of a mixed-mode single-integral type cohesive zone with reorienting surfaces of separation

To model material ductile failure and crack propagation, cohesive zone elements can be embedded along potential fracture paths in a finite element simulation. When damage criteria are met, elements in the mesh decohere, simulating the formation and propagation of a crack. In this paper, we present a novel computational algorithm based on finite deformation theory, essential to modeling crack initiation and growth in solids undergoing large deformations. This new algorithm was formulated within a Lagrangian frame of reference to extend previous cohesive zone algorithms to include modeling crack growth in finite deformation contexts. The local coordinate system, necessary for defining an embedded cohesive zone, is constructed based upon the current configuration and is updated within the nonlinear iteration process, thereby resulting in the convergence of the solution for a growing crack in a large deformation quasi-static setting. The model’s accuracy was demonstrated by comparing finite element model simulation results with the analytic case of a constant surface separation, as shown in the verification examples. The power and efficacy of the algorithm to capture large deformations during crack growth were then demonstrated with a double cantilever beam example case. It indicates that the model can be applied to a variety of physical circumstances for predicting crack initiation and growth with delamination and fracture.

42 ENGINEERING↗

Thermal Pressure in the Thermal Equation of State for Solid and a Proposed Substitute

Abstract The thermal equation of state (TEOS) for solids is a mathematic model among pressure, temperature and density, and is essential for geophysical, geochemical, and other high pressure–temperature (high P–T) researches. However, in the last few decades, there has been a growing concern about the accuracy of the pressure scales of the calibrants, and efforts have been made to improve it by either introducing a reference standard or building new thermal pressure models. The existing thermal equation of state, P ( V , T ) = P ( V , T 0 ) + P th ( V , T ), consists of an isothermal compression and an isochoric heating, while the thermal pressure is the pressure change in the isochoric heating. In this paper, we demonstrate that, for solids in a soft pressure medium in a diamond anvil cell, the thermal pressure can neither be determined from a single heating process, nor from the thermal pressure of its calibrant. To avoid the thermal pressure, we propose to replace the thermal pressure with a well-known thermal expansion model, and integrate it with the isothermal compression model to yields a Birch–Murnaghan-expansion TEOS model, called VPT TEOS. The predicted pressure of MgO and Au at ambient pressure from Birch–Murnaghan-expansion VPT TEOS model matches the experimental pressure of zero (0) GPa very well, while the pressure prediction from the approximated Anderson PVT TEOS exhibit a big deviation and a wrong trend.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

SU-E-J-148: Tools for Development of 4D Proton CT

Purpose: To develop tools for performing 4D proton computed tomography (CT). Methods: A suitable patient with a tumor in the right lower lobe was selected from a set of 4D CT scans. The volumetric CT images formed the basis for calculating the parameters of a breathing model that allows reconstruction of a static reference CT and CT images in each breathing phase. The images were imported into the TOPAS Monte Carlo simulation platform for simulating an experimental proton CT scan with 45 projections spaced by 4 degree intervals. Each projection acquired data for 2 seconds followed by a gantry rotation for 2 seconds without acquisition. The scan covered 180 degrees with individual protons passing through a 9-cm slab of the patient’s lung covering the moving tumor. An initial proton energy sufficient for penetrating the patient from all directions was determined. Performing the proton CT simulation, TOPAS provided output of the proton energy and coordinates registered in two planes before and after the patient, respectively. The set of projection data was then used with an iterative reconstruction algorithm to generate a volumetric proton CT image set of the static reference image and the image obtained under breathing motion, respectively. Results: An initial proton energy of 230 MeV was found to be sufficient, while for an initial energy of 200 MeV a substantial number of protons did not penetrate the patient. The reconstruction of the static reference image set provided sufficient detail for treatment planning. Conclusion: We have developed tools to perform studies of proton CT in the presence of lung motion based on the TOPAS simulation toolkit. This will allow to optimize 4D reconstruction algorithms by synchronizing the acquired proton CT data with a breathing signal and utilizing a breathing model obtained prior to the proton CT scan. This research has been supported by the National Institute Of Biomedical Imaging And Bioengineering of the National Institutes of Health under Award Number R01EB013118.

60 APPLIED LIFE SCIENCES↗

Non-Kyoto Radiative Forcing in Long-Run Greenhouse Gas Emissions and Climate Change Scenarios

Climate policies designed to achieve climate change objectives must consider radiative forcing from the Kyoto greenhouse gas, as well as other forcing constituents, such as aerosols and tropospheric ozone. Net positive forcing leads to global average temperature increases. Modeling of non-Kyoto forcing is a relatively new component of climate management scenarios. Five of the nineteen models in the EMF-27 Study model both Kyoto and non-Kyoto forcing. This paper describes and assesses current non-Kyoto radiative forcing modeling within these integrated assessment models. The study finds negative forcing from aerosols masking significant positive forcing in reference non-climate policy projections. There are however large differences across models in projected non-Kyoto emissions and forcing, with differences stemming from differences in relationships between Kyoto and non-Kyoto emissions and fundamental differences in modeling structure and assumptions. Air pollution and non-Kyoto forcing decline in the climate policy scenarios. However, non-Kyoto forcing appears to be influencing mitigation results, including allowable carbon dioxide emissions, and further evaluation is merited. Overall, there is substantial uncertainty related to non-Kyoto forcing that must be considered.

Rose, Steven K.↗

Logistic Regression in Clinical Studies

• A logistic regression model is used when the outcome of interest is binary. The term “logistic” refers to the underlying “logit” (log odds) function that is used to model the binary outcome. • Odds ratios are produced from a logistic regression model and have a useful interpretation. • Tips, tricks and concepts used to fit logistic regression models are similar to those used in linear regression models. • Modeling building that is knowledge-based rather than automatic is preferred in most applications of logistic regression. • A logistic regression model that is overparameterized (ie, too many variables for too few events) can result in odds ratios that are implausibly large and confidence intervals that are wide and uninterpretable. These types of “overfitted” models should be avoided. • Logistic regression models can be fit using most standard statistical software.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

MELCOR 2.2 Benchmarks of Peach Bottom NUREG/CR 7155 Uncertainty Analysis

The U.S. Nuclear Regulatory Commission (NRC) performed a first-of-a-kind uncertainty analysis (UA) of the accident progression, radiological releases, and offsite consequences for the State-of- the-Art Reactor Consequence Analyses (SOARCA) of an unmitigated long-term station blackout (LTSBO) severe accident scenario at the Peach Bottom Atomic Power Station. The objective of the UA was to evaluate the robustness of the SOARCA deterministic "best estimate results and conclusions documented in NUREG-1935, and to develop insight into the overall sensitivity of the SOARCA results to uncertainty in key modeling inputs. The study was completed in 2015 and documented in NUREG/CR-7155. Since 2015, two other SOARCA UAs were completed for two pressurized water reactor (PWR) plants. The PWR UAs incrementally updated the approach and methodology, including using the latest release of the MELCOR 2.2 computer code. There were also advances made in the state-of-the-art modeling related to NRC efforts using the Peach Bottom model in NUREG-2206, which provided the technical basis for the containment protection and release reduction rulemaking for boiling water reactors with Mark I and Mark II containments. This report documents the input model changes from the NUREG/CR-7155 study and performs a small number of reference calculations to assess the changes of the new computer code and the model input updates. The objective of the work is to verify whether the updated Peach Bottom MELCOR model and updated version of MELCOR support the conclusions formed in the Peach Bottom SOARCA UA by performing these representative calculations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A State-Space Model for Stability Boundary Analysis of Grid-Following Voltage Source Converters Considering Grid Conditions

With the growing significance of renewable energy resources and energy storage systems, the number of grid-connected inverters has been rising at an increasingly rapid pace. Generally, these inverters are directly integrated with the distribution network by synchronizing with the grid voltage at the point of common coupling. However, the low grid strength and varying R/X ratios, as the common characteristics of most distribution networks or weak grids, can lead to dynamic interactions that comprise stability and limit the power transfer capacity of grid-connected inverters. To ensure stable operation of the inverters, researchers must determine the stability boundary, described as the maximum power transfer capacity of grid-connected inverters under the premise of maintaining system small-signal stability. For this purpose, we propose to formulate a state-space model of the system in the synchronously rotating dq-frame of reference and perform eigenvalue analysis to determine the stability boundary. With a detailed model of the control structure and parameters of the grid-connected inverters, the stability boundary is identified as a surface with respect to different grid strengths and R/X ratios. Case study results of proposed eigenvalue analysis are compared with those of admittance model-based stability analysis as well as time-domain simulation using a switching model in Matlab/Simulink, validating the effectiveness and accuracy of the proposed eigenvalue analysis for stability boundary identification.

grid-connected inverters↗

Computational fluid dynamics investigations of flow, heat transfer, and oxidation in heat recovery steam generator

Modern heat recovery steam generators (HRSGs) operate at elevated temperatures, leading to the formation of oxides inside the tubes of heat exchangers (HXs). This oxide growth reduces the heat recovery efficiency. Moreover, after reaching a certain critical thickness, some oxide scales detach from the tube surface (exfoliation), causing erosion damage to the components downstream. Predicting the metal temperature distribution and associated oxide thickness in the HXs of an HRSG can aid in mitigating these problems. Here, a computational fluid dynamics (CFD) model was developed within the commercial code STAR-CCM+ for the prediction of fluid flow, conjugate heat transfer, and associated oxidation in HRSGs. Moreover, a new Porous Media Model (PMM) method was developed to model the fin effect on the heat transfer in HX, which can substantially reduce the prohibitive computational costs of fin meshes. The developed CFD model was used to conduct a high-fidelity simulation of a real-scale HRSG to investigate flow, heat transfer, and oxide growth. The calculated oxide thickness on different tubes can be used to identify HX regions that require oxide-resistant coatings to prevent exfoliation and ensuing damages. Furthermore, this CFD framework can serve as a reference for future studies that intend to model and investigate high-temperature oxidation in HXs used for any applications.

42 ENGINEERING↗

Wind load analysis of a new linear Fresnel receiver assembly design

A linear Fresnel collector includes a low-profile reflector array and a receiver assembly with one or more absorber tubes and an optional secondary reflector. This combined optical system concentrates sunlight and converts it into thermal energy. The design of a receiver assembly is critical to the performance of a linear Fresnel collector. A position deviation of a few centimeters for the receiver assembly can result in notably reduced performance, thus leading to a direct loss in revenue associated with thermal power production. Wind load is one of the most significant environmental factors that can alter the optical - and therefore thermal - performance of a solar power system due to displacements after installation. At the same time, an over-designed receiver assembly may add unnecessary construction cost to a typically high-cost-constrained system. Thus, wind load analysis is particularly important when considering optimal engineering design of a receiver assembly and its supporting structure to cost-effectively mitigate the impacts of wind. In this study, a detailed computational fluid dynamics (CFD) model is adopted to derive the wind load of a commercial linear Fresnel receiver assembly. This wind load is then used as a reference to optimize the detailed engineering design. The CFD model is first carefully developed and benchmarked within a critical regime toward turbulence. The drag force, lift force, and vortex-shedding frequencies are derived at both the operating and survival wind-speed limits for target project deployment locations. The wind load analysis results provide a valuable reference for future engineering design and prototyping.

17 WIND ENERGY↗

System of Systems Model Development for Evaluating EMP Resilient Grid Mitigation Strategies

This Laboratory Directed Research and Development (LDRD) project focused on understanding the mathematical relationships that can be used in assessing the value of executing various EMP mitigation strategies on the grid. This is referred to as the EMP Resilient Grid Value Model. Because the range of mitigation strategies can contain widely differing characteristics (operational vs. technological), it is necessary to compute functions of many interrelated metrics at varying levels of fidelity that will be used to provide feedback as to the cost/benefit relationship of any proposed strategy. The value model is a hierarchical decomposition of a system of systems (SoS) model down to a grid circuit model. The model is intended to be suitable for use in subsequent decision support optimization for resilience to EMP events. The metric set goes beyond direct, technical impacts on the electrical grid to include ancillary impacts on dependent infrastructure and enterprise concerns (water, DoD, transportation, etc.).

24 POWER TRANSMISSION AND DISTRIBUTION↗