Search NASA⌕ Search

SEARCH · Search NASA

Results for “choice modeling”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Late-time post-merger modeling of a compact binary: effects of relativity, r-process heating, and treatment of transport

Detectable electromagnetic counterparts to gravitational waves from compact binary mergers can be produced by outflows from the black hole-accretion disk remnant during the first 10 s after the merger. Two-dimensional axisymmetric simulations with effective viscosity remain an efficient and informative way to model this late-time post-merger evolution. In addition to the inherent approximations of axisymmetry and modeling turbulent angular momentum transport by a viscosity, previous simulations often make other simplifications related to the treatment of the equation of state and turbulent transport effects. In this paper, we test the effect of these modeling choices. By evolving with the same viscosity the exact post-merger initial configuration previously evolved in Newtonian viscous hydrodynamics, we find that the Newtonian treatment provides a good estimate of the disk ejecta mass but underestimates the outflow velocity. We find that the inclusion of heavy nuclei causes a notable increase in ejecta mass. An approximate inclusion of r-process effects has a comparatively smaller effect, except for its designed effect on the composition. Diffusion of composition and entropy, modeling turbulent transport effects, has the overall effect of reducing ejecta mass and giving it a speed with lower average and more tightly-peaked distribution. Also, we find significant acceleration of outflow even at distances beyond 10 000 km, so that thermal wind velocities only asymptote beyond this radius and at higher values than often reported.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Constraining computational modeling of indirect drive double shell capsule implosions using experiments

Double shell capsule implosions are an alternative approach to achieving alpha heating on the National Ignition Facility. Current machining techniques construct the outer shell as two hemispheres that are glued together, and the deuterium and tritium (DT) liquid inside the inner shell will be injected by a fill tube. These features introduce asymmetries and jetting that may disrupt the confinement of the DT fuel if not carefully controlled. Simulations indicate that in order to achieve high yields in the laboratory, these features as well as susceptibility to the Rayleigh–Taylor instability (RTI) must be mitigated. Furthermore, due to uncertainties in computational models and the expense of using the best physics models at adequate resolution in three dimensions, our computational modeling must be constrained by experiments. We report on the results of recent hydrogrowth radiography and dual-axis keyhole experiments with double shell targets that have been used to evaluate our modeling of the outer shell joint as well as the impacts of high-energy x-ray preheat that strongly impacts RTI growth. Our simulations show good agreement with the experimental data and inform several important modeling choices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Revisiting Power Systems Time-Domain Simulation Methods and Models

The changing nature of power systems dynamics is challenging present practices related to modeling and study of system-level dynamic behavior. While developing new techniques and models to handle the new modeling requirements, it is also critical to review some of the terminology used to describe existing simulation approaches and the embedded assumptions. This article provides a first-principles review of the simplifications and transformations commonly used in the formulation of time-domain simulation models. It introduces a taxonomy and classification of time-domain simulation models depending on their frequency bandwidth, network representation, and software availability. Furthermore, it focuses on the fundamental aspects of averaging techniques, and model reduction approaches that result in modeling choices, and discusses the associated challenges and opportunities of applying these methods in systems with large shares of Inverter Based Resources (IBRs). The article concludes with an illustrative simulation that compares the trajectories of an IBR-dominated system.

behavioral sciences↗
DOE OSTI2020

Nek5000 developments in support of industry and the NRC

This year, the Nuclear Energy Advanced Modeling Simulation program (NEAMS) thermal-hydraulics verification and validation (V&V) work has focused in three areas of Nek5000 V&V-driven development. First, in a close collaborative effort with the U. S. Nuclear Regulatory Commission (NRC) staff, we have continued V&V efforts for the HYMERES-2 project using the OECD/NEA sponsored testing in the PSI PANDA facility. This year’s focus of ANL-NRC collaboration involves Nek5000 setups and validation for a range of problems relevant to and including the HYMERES-2 benchmark from PSI. The primary outcome of this year efforts is a more efficient geometry and inlet modeling simplification after a careful sensitivity study of the inlet profiles and pipe geometries. The resulting modeling choice of a short recycling/fully-developed turbulent inlet is within the experimental uncertainty estimate. This finding simplifies the next step of the cross-V&V HYMERES-2 project. In addition, the ANL team continue to provide assistance to the NRC staff in the form of Nek5000 application support in general and on the use of the HPC platforms of ALCF and INL in particular. This supports the NRC’s assessment of Nek5000 for use with the NRC Blue CRAB code suite. Second, we have implemented and tested more robust model of URANS, namely the k – τ model, a variant of the k-ω model, along with other improvements to RANS Nek5000 modeling in general. Because of its demonstrated robustness and stability, the k – τ model is the only RANS model that has been implemented in the new GPU version of the Nek5000 code, nekRS. Lastly, we report the initial implementation of Jacobian-free Newton Krylov approach to the direct Newton method for steady fluid solvers aimed at acceleration of RANS modeling and at IC improvement for LES campaigns. Also leveraging the Exascale Computing Project (ECP) ANL/CEED & SMR team’s software development effort to support NEAMS problems at large scale of the advanced computing architectures, NekRS, a GPU variant of Nek5000, built on top of kernels from libParanumal using OCCA for portability, has been successfully run on the full system of Summit (4608 nodes, 27648 GPUs).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Temperatures of Anvil Clouds and Radiative Tropopause in a Wide Array of Cloud-Resolving Simulations

Here, we present 123 cloud-resolving simulations to study how temperatures of anvil clouds and radiative tropopause (RT) change with surface warming. Our simulation results show that the RT warms at approximately the same rate as anvil clouds. This relationship persists across a variety of modeling choices, including surface temperature, greenhouse gas concentration, and the representation of radiative transfer. We further show that the shifting ozone profile associated with climate warming may give rise to a fixed RT temperature as well as a fixed anvil temperature. This result points to the importance of faithful treatment of ozone in simulating clouds and climate change; the robust anvil–RT relationship may also provide alternative ways to understand what controls anvil temperature.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Efficient Neural Network Approaches for Conditional Optimal Transport with Applications in Bayesian Inference

In this work, we present two neural network approaches that approximate the solutions of static and dynamic conditional optimal transport (COT) problems. Both approaches enable conditional sampling and conditional density estimation, which are core tasks in Bayesian inference—particularly in the simulation-based (“likelihood-free”) setting. Our methods represent the target conditional distribution as a transformation of a tractable reference distribution. Obtaining such a transformation, chosen here to be an approximation of the COT map, is computationally challenging even in moderate dimensions. To improve scalability, our numerical algorithms use neural networks to parameterize candidate maps and further exploit the structure of the COT problem. Our static approach approximates the map as the gradient of a partially input convex neural network. It uses a novel numerical implementation to increase computational efficiency compared to state-of-the-art alternatives. Our dynamic approach approximates the conditional optimal transport via the flow map of a regularized neural ODE; compared to the static approach, it is slower to train but offers more modeling choices and can lead to faster sampling. We demonstrate both algorithms numerically, comparing them with competing state-of-the-art approaches, using benchmark datasets and simulation-based Bayesian inverse problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

How to estimate soil organic carbon stocks of agricultural fields? perspectives using ex-ante evaluation

Estimating soil organic carbon (SOC) stocks of agricultural fields has a range of important applications from development of sustainable management practices to monitoring carbon stocks. There are many estimation strategies with the potential for more reliable estimates of SOC stock and more efficient use of soil sampling and analysis resources, especially by leveraging readily available auxiliary information such as remote sensing. However, concrete guidance for strategy selection is lacking. This study narrows this gap with a comparison of strategies for estimating deep SOC stock (0–60 cm) in a prototypical field. Using high density SOC stock measurements and simulation, we built on past studies by 1) ex-ante evaluating a large number of strategy options, 2) using a Bayesian approach to quantify the uncertainty of the comparison, and 3) considering multiple Bayesian models to assess sensitivity to this modeling choice. We found that, using readily available auxiliary information, both balanced and stratified sampling offer substantial improvements over simple random sampling. The auxiliary information most important for this improvement is a Sentinel-2 SOC index = blue / (green × red), followed by the topographic wetness index. We found that these results are robust to the choice of mapping method, but that there is uncertainty in the magnitude of improvement. Here, we recommend future studies implement this Bayesian approach for simulated ex-ante evaluation of SOC stock estimation strategies across more fields to investigate the generalizability of these findings.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Learning interpretable surface elasticity properties from bulk properties via neural network equation learners

Surface elasticity is central to understanding the mechanics and stability of surfaces and interfaces. It is characterized by quantities such as surface tension, residual surface stress, and surface stiffness. However their analytical expressions are typically difficult to derive from atomistic data, and depend strongly on modeling choices. This work presents a neural network-based equation learner which combines customized activation functions and connection-based pruning to discover parsimonious, closed-form equations for surface elasticity from atomistic simulations. Applying the method to seven face-centered cubic (FCC) metals, our equation learner uncovers interpretable equations that describe both low-Miller index and high-Miller index surface properties, capturing long-tail property distributions accurately. The discovered expressions are decoupled into two components: a universal, geometry-driven orientation function, and material-specific baseline coefficients. We find that lower-order properties such as surface tension are fundamentally geometry dependent, while higher-order properties such as surface stress and elasticity show more complex geometry and material dependence. We also relate material dependent coefficients to bulk properties, forming a clear map from bulk material properties to surface elasticity. Overall, this approach demonstrates that interpretable neurosymbolic machine learning can bridge the gap between atomistic simulations and physical laws, enabling the discovery of generalizable structure–property relationships for materials science phenomena such as surface elasticity.

Equation learning↗
DOE OSTI2025

Poor reliability of public charging stations can impede the growth of the electric vehicle market

How does the reliability of public charging infrastructure affect electric vehicle (EV) adoption? Substantial public and private investments are expanding EV charging networks, but concerns are growing about the poor reliability of existing chargers and its potential impacts on EV adoption. Using data from a nationwide survey, we employ a choice model to quantify the effects of perceived charging reliability on Americans’ intentions to purchase new or used EVs. By randomly assigning participants to receive information characterizing public charging as either very reliable or very unreliable, we show a causal effect of reliability perceptions on EV purchase intentions. In conclusion, we find that differences in perceived reliability are equivalent to changing price by 32 % of purchasing budget or changing range by 366 miles, underscoring the importance of reliable public charging.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Value-added catalog of physical properties for more than 1.3 million galaxies from the DESI survey

We present an extensive catalog of the physical properties of more than a million galaxies investigated with the Dark Energy Spectroscopic Instrument (DESI), one of the largest spectroscopic surveys to date. Spanning a full range of target types, including emission-line galaxies, luminous red galaxies, and quasars, our survey encompasses an unprecedented range of spectroscopic redshifts, all the way from 0 to 6. The physical properties, such as stellar masses and star formation rates, were derived via the CIGALE spectral energy distribution (SED) fitting code accounting for the contribution coming from active galactic nuclei (AGNs). Based on the modeling of the optical-mid-infrared (grz supplemented with WISE photometry) SEDs, we studied the galaxy properties with respect to their location on the main sequence. We have revised the dependence of stellar mass estimates on model choices and on the availability of WISE photometry. Indeed, the WISE data are required to minimize the misclassification of star-forming galaxies as AGNs. The lack of WISE bands in SED fits leads to elevated AGN fractions for 68% of star-forming galaxies identified using emission line diagnostic diagrams, but this does not significantly affect their stellar mass or star formation estimates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Considerations for the modeling of the laser ablation region of ICF targets with Lagrangian simulations

Recently, much effort has been dedicated to the improvement of models and modeling choices utilized in radiation hydrodynamic simulations of direct drive inertial confinement fusion experiments in an effort to improve their predictive capability. In this paper, we consider the choice in mesh for the simulation of the laser ablation of a direct-drive-like target and compare Lagrangian simulations with various mesh zoning choices with Eulerian simulations with fixed resolution in the laser energy deposition region. Using these simulations, we demonstrate how errors in ablation pressure, laser deposition rate, shock speed, and density profile arise from insufficient zoning following from the conservation of mass of Lagrangian zones. These considerations place stringent requirements on the initial t = 0 zoning in the solid density shell for simulations aiming at resolving the ablation and laser absorption region. However, with sufficiently fine zoning in the t = 0 shell, agreement with Eulerian simulations and analytic scaling laws can be recovered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Winter roosting ecology of tricolored bats ( Perimyotis subflavus ) in trees and bridges

Abstract Tricolored bats (Perimyotis subflavus) that roost in subterranean hibernacula have experienced precipitous declines from white-nose syndrome (WNS); however, understudied populations also use during winter non-subterranean roosts such as tree cavities, bridges, and foliage. Our objectives were to determine winter roost use by tricolored bats in an area devoid of subterranean roosts, determine roost microclimates to relate them to growth requirements of the fungal causal agent of WNS, and determine habitat factors influencing winter tree selection. From November to March 2017–2019, we used radiotelemetry to track 15 bats to their day roosts in the upper Coastal Plain of South Carolina and recorded microclimates in accessible tree cavities and bridges. We also characterized habitat and tree characteristics of 24 used trees and 153 random, available trees and used discrete choice models to determine selection. Roost structures included I-beam bridges, cavities in live trees, and foliage. Bridges were warmer and less humid than cavities. Roost temperatures often were amenable to fungal growth (< 19.5°C) but fluctuated widely depending on ambient temperatures. Bats used bridges on colder days (8.7°C ± 5.0 SD) and trees on warmer days (11.3°C ± 5.4). Bats selected low-decay trees closer to streams in areas with high canopy closure and cavity abundance. Bats also appeared to favor hardwood forests and avoid pine forests. Our results suggest that access to multiple roost microclimates might be important for tricolored bats during winter, and forest management practices that retain live trees near streams and foster cavity formation in hardwood forests likely will benefit this species. Our results also suggest tricolored bats using bridge and tree roosts might be less susceptible to WNS than bats using subterranean hibernaculum roosts. Thus, forests in areas without subterranean hibernacula in the southeastern United States that support bats during winter might represent important refugia from WNS for multiple species.

Zoology↗
DOE OSTI2021

The Dark Energy Survey supernova programme: modelling selection efficiency and observed core-collapse supernova contamination

ABSTRACT The analysis of current and future cosmological surveys of Type Ia supernovae (SNe Ia) at high redshift depends on the accurate photometric classification of the SN events detected. Generating realistic simulations of photometric SN surveys constitutes an essential step for training and testing photometric classification algorithms, and for correcting biases introduced by selection effects and contamination arising from core-collapse SNe in the photometric SN Ia samples. We use published SN time-series spectrophotometric templates, rates, luminosity functions, and empirical relationships between SNe and their host galaxies to construct a framework for simulating photometric SN surveys. We present this framework in the context of the Dark Energy Survey (DES) 5-yr photometric SN sample, comparing our simulations of DES with the observed DES transient populations. We demonstrate excellent agreement in many distributions, including Hubble residuals, between our simulations and data. We estimate the core collapse fraction expected in the DES SN sample after selection requirements are applied and before photometric classification. After testing different modelling choices and astrophysical assumptions underlying our simulation, we find that the predicted contamination varies from 7.2 to 11.7 per cent, with an average of 8.8 per cent and an r.m.s. of 1.1 per cent. Our simulations are the first to reproduce the observed photometric SN and host galaxy properties in high-redshift surveys without fine-tuning the input parameters. The simulation methods presented here will be a critical component of the cosmology analysis of the DES photometric SN Ia sample: correcting for biases arising from contamination, and evaluating the associated systematic uncertainty.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Stringent σ8 constraints from small-scale galaxy clustering using a hybrid MCMC + emulator framework

ABSTRACT We present a novel simulation-based hybrid emulator approach that maximally derives cosmological and Halo Occupation Distribution (HOD) information from non-linear galaxy clustering, with sufficient precision for DESI Year 1 (Y1) analysis. Our hybrid approach first samples the HOD space on a fixed cosmological simulation grid to constrain the high-likelihood region of cosmology + HOD parameter space, and then constructs the emulator within this constrained region. This approach significantly reduces the parameter volume emulated over, thus achieving much smaller emulator errors with fixed number of training points. We demonstrate that this combined with state-of-the-art simulations result in tight emulator errors comparable to expected DESI Y1 LRG sample variance. We leverage the new abacussummit simulations and apply our hybrid approach to CMASS non-linear galaxy clustering data. We infer constraints on σ8 = 0.762 ± 0.024 and fσ8(zeff = 0.52) = 0.444 ± 0.016, the tightest among contemporary galaxy clustering studies. We also demonstrate that our fσ8 constraint is robust against secondary biases and other HOD model choices, a critical first step towards showcasing the robust cosmology information accessible in non-linear scales. We speculate that the additional statistical power of DESI Y1 should tighten the growth rate constraints by at least another 50–60 ${{\ \rm per\ cent}}$, significantly elucidating any potential tension with Planck. We also address the ‘lensing is low’ tension, which we find to be in the same direction as a potential tension in fσ8. We show that the combined effect of a lower fσ8 and environment-based bias accounts for approximately $50{{\ \rm per\ cent}}$ of the discrepancy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Hyper Suprime-Cam Year 3 results: Cosmology from cosmic shear power spectra

Here, we measure weak lensing cosmic shear power spectra from the 3-year galaxy shear catalog of the Hyper Suprime-Cam (HSC) Subaru Strategic Program imaging survey. The shear catalog covers 416 deg 2 of the northern sky, with a mean i-band seeing of 0.59 arcsec and an effective galaxy number density of 15 arcmin —2 within our adopted redshift range. With an i-band magnitude limit of 24.5 mag, and four tomographic redshift bins spanning 0.3 ≤ z ph ≤ 1.5 based on photometric redshifts, we obtain a high-significance measurement of the cosmic shear power spectra, with a signal-to-noise ratio of approximately 26.4 in the multipole range 300 < ℓ < 1800. The accuracy of our power spectrum measurement is tested against realistic mock shear catalogs, and we use these catalogs to get a reliable measurement of the covariance of the power spectrum measurements. We use a robust blinding procedure to avoid confirmation bias, and model various uncertainties and sources of bias in our analysis, including point spread function systematics, redshift distribution uncertainties, the intrinsic alignment of galaxies and the modeling of the matter power spectrum. For a flat Λ CDM model, we find S 8 ≡ σ 8 (Ω m /0.3) 0.5 = ${0.776}_{—0.033}^{+0.032}$, which is in excellent agreement with the constraints from the other HSC Year 3 cosmology analyses, as well as those from a number of other cosmic shear experiments. This result implies a ~2σ-level tension with the Planck 2018 cosmology. We study the effect that various systematic errors and modeling choices could have on this value, and find that they can shift the best-fit value of S 8 by no more than ~ 0.5σ, indicating that our result is robust to such systematics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

SPT clusters with DES and HST weak lensing. I. Cluster lensing and Bayesian population modeling of multiwavelength cluster datasets

We present a Bayesian population modeling method to analyze the abundance of galaxy clusters identified by the South Pole Telescope (SPT) with a simultaneous mass calibration using weak gravitational lensing data from the Dark Energy Survey (DES) and the Hubble Space Telescope (HST). We discuss and validate the modeling choices with a particular focus on a robust, weak-lensing-based mass calibration using DES data. For the DES Year 3 data, we report a systematic uncertainty in weak-lensing mass calibration that increases from 1% at z = 0.25 to 10% at z = 0.95 , to which we add 2% in quadrature to account for uncertainties in the impact of baryonic effects. We implement an analysis pipeline that joins the cluster abundance likelihood with a multiobservable likelihood for the Sunyaev-Zel’dovich effect, optical richness, and weak-lensing measurements for each individual cluster. We validate that our analysis pipeline can recover unbiased cosmological constraints by analyzing mocks that closely resemble the cluster sample extracted from the SPT-SZ, SPTpol ECS, and SPTpol 500d surveys and the DES Year 3 and HST-39 weak-lensing datasets. This work represents a crucial prerequisite for the subsequent cosmological analysis of the real dataset.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Masses of Sunyaev-Zel’dovich galaxy clusters detected by the Atacama Cosmology Telescope: Stacked lensing measurements with Subaru HSC year 3 data

We present a stacked lensing analysis of 96 galaxy clusters selected by the thermal Sunyaev-Zel’dovich (SZ) effect in maps of the cosmic microwave background (CMB). We select foreground galaxy clusters with a 5 σ -level SZ threshold in CMB observations from the Atacama Cosmology Telescope, while we define background source galaxies for the lensing analysis with secure photometric redshift cuts in Year 3 data of the Subaru Hyper Suprime Cam survey. We detect the stacked lensing signal in the range of 0.1 < R [ h - 1 Mpc ] < 100 in each of three cluster redshift bins, 0.092 < z ≤ 0.445 , 0.445 < z ≤ 0.695 , and 0.695 < z ≤ 1.180 , with 32 galaxy clusters in each bin. The cumulative signal-to-noise ratios of the lensing signal are 14.6, 12.0, and 6.6, respectively. Using a halo-based forward model, we then constrain statistical relationships between the mass inferred from the SZ observation (i.e. SZ mass) and the total mass derived from our stacked lensing measurements. At the average SZ mass in the cluster sample ( 2.1 - 2.4 × 10 14 h - 1 M ⊙ ), our likelihood analysis shows that the average total mass differs from the SZ counterpart by a factor of 1.3 ± 0.2 , 1.6 ± 0.2 , and 1.6 ± 0.3 (68%) in the aforementioned redshift ranges, respectively. Our limits are consistent with previous lensing measurements, and we find that the cluster modeling choices can introduce a 1 σ -level difference in our parameter inferences.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Probing the Pulsar Explanation of the Galactic-Center GeV Excess Using Continuous Gravitational-Wave Searches

Over 10 years ago, Fermi observed an excess of GeV gamma rays from the Galactic Center whose origin is still under debate. One explanation for this excess involves annihilating dark matter, another requires an unresolved population of millisecond pulsars concentrated at the Galactic Center. In this work, we use the results from LIGO and Virgo’s most recent all-sky search for quasimonochromatic, persistent gravitational-wave signals from isolated neutron stars, which is estimated to be about 20%–50% of the population, to determine whether unresolved millisecond pulsars could actually explain this excess. First, we choose a luminosity function that determines the number of millisecond pulsars required to explain the observed excess. Then, we consider two models for deformations on millisecond pulsars to determine their ellipticity distributions, which are directly related to their gravitational-wave radiation. Lastly, based on null results from the O3 frequency-Hough all-sky search for continuous gravitational waves, we find that a large set of the parameter space in the pulsar luminosity function can be excluded. We also evaluate how these exclusion regions may change with respect to various model choices. Further, our results are the first of their kind and represent a bridge between gamma-ray astrophysics, gravitational-wave astronomy, and dark-matter physics.

79 ASTRONOMY AND ASTROPHYSICS↗