Search NASA⌕ Search

SEARCH · Search NASA

Results for “ensemble simulation”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Modeling aerosol transmission spectra from n(λ) and k(λ) infrared optical constants measurements of organic liquids and solids

The effects of light scattering and refraction play significantly different roles for aerosols than for bulk materials, making it challenging to identify aerosolized chemicals using traditional spectral methods or spectral reference libraries. Due to a potentially infinite number of particle morphologies, sizes, and compositions, constructing a database of laboratory-measured aerosol spectra is not a practical solution. Here, as an alternative approach, the measured n / k optical vectors of two example organic materials (diethyl phthalate and D-mannitol) are used in combination with particle absorption / scattering theory (Mie theory and FDTD) and the Beer-Lambert law to generate a series of synthetic infrared transmission / scattered light spectra. The synthetic spectra show significant differences versus simple slab transmission spectra, even for small changes in particle size (e.g., 5 vs. 10 µm) for both single particles and ensembles, potentially serving as useful reference data for aerosol sensing. For spherical single particles with diameters of 1 to 10 µm, FDTD simulations predict changes in the magnitudes of spectral shifts and the shapes of the peaks vs. particle size with only small deviations from Mie theory predictions, yet reliably capture the direction of the shifts. Typical spectral peak shifts in the longwave infrared correspond to Δλ ∼0.20 µm (∼34 cm -1 ) when compared to corresponding slab transmission spectra. Additionally, synthetic spectra generated from the n / k values derived using two different methods (KBr pellet transmission and single-angle reflectance) are compared using the Mie theory model.

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

A Simulation Study of Quantum Clock Synchronization Using Teleportation

An important requirement in implementing distributed computing and sensing application is the synchronization of clocks at various locations. The Internet relies on the Network Time Protocol (NTP), which synchronizes clocks with accuracy in the order of milliseconds. More recently an ensemble of atomic clocks is used for navigation based on GPS. These clocks are highly accurate and provide time with very low uncertainty. Even so, many physics experiments such as distributed LIGO-based systems may require more accurate clock synchronization that is achievable using quantum entanglement. This requires the deployment of a network of quantum clocks synchronized by exploiting entangled atomic clock qubits. In this paper, we carry out a simulation study of synchronizing a network of quantum clocks interconnected by a fiber plant that supports the ESnet; the latter is used to support the classical communication needed for teleportation. We consider an existing protocol for synchro-nizing the atomic clock qubits that relies on the GHZ states. To assess the performance of the protocol we developed a discrete-event simulation of the network using IBM Qiskit framework for underlying quantum gate operations and measurements. The simulation results shed light on the resources required in terms of the number entangled qubits and the time needed to achieve the synchronization of different number of nodes in ESnet.

Kiran, Mariam↗
DOE OSTI2024

Recalibration of missing low-frequency variability and trends in the North Atlantic Oscillation

Abstract Multi-decadal trends in the wintertime North Atlantic Oscillation (NAO) are under-represented by coupled general circulation models (CGCMs), consistent with a lack of autocorrelation in their NAO index series. This study proposes and tests two simple “reddening” approaches for correcting this problem in simulated indices based on simple one parameter short-term (AR; Auto-Regressive order 1) and long-term (FD; Fractional-Difference) time series filters. Using CGCMs from the Coupled Model Intercomparison Project Phase 6 (CMIP6), the FD filter successfully improves the autocorrelation structure of the NAO, and in turn the simulation of extreme trends, while the AR filter is less successful. The 1963–1993 NAO trend is the maximum 31-year trend in the historical period. Raw CGCMs underestimate the likelihood of this trend by a factor of ten but this discrepancy is corrected after reddening. CMIP6 future projections show that long-term (2024–2094) NAO ensemble mean trends systematically increase with the magnitude of radiative forcing: -2.4 to 3.5 hPa/century for low-to-high forcing after reddening (more than double the range using raw output). The related likelihood of future maximum 31year trends comparable to 1963–1993 ranges from 3 to 7% whereas none of these CMIP6 projections simulate this without reddening. Near-term projections of the next 31 years (2024–2054) are less sensitive than long term trends to the future scenario, showing weak-to-no forced trend. However, reddening increases the ensemble range by 74% (to +/-1 standard deviation/decade), which could increase/decrease regional climate change signals in the Northern Hemisphere by magnitudes that are underestimated when using raw CGCM output.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

Structure-Based Identification of Novel Histone Deacetylase 4 (HDAC4) Inhibitors

Histone deacetylases (HDACs) are important cancer drug targets. Existing FDA-approved drugs target the catalytic pocket of HDACs, which is conserved across subfamilies (classes) of HDAC. However, engineering specificity is an important goal. Herein, we use molecular modeling approaches to identify and target potential novel pockets specific to Class IIA HDAC-HDAC4 at the interface between HDAC4 and the transcriptional corepressor component protein NCoR. These pockets were screened using an ensemble docking approach combined with consensus scoring to identify compounds with a different binding mechanism than the currently known HDAC modulators. Binding was compared in experimental assays between HDAC4 and HDAC3, which belong to a different family of HDACs. HDAC4 was significantly inhibited by compound 88402 but not HDAC3. Two other compounds (67436 and 134199) had IC50 values in the low micromolar range for both HDACs, which is comparable to the known inhibitor of HDAC4, SAHA (Vorinostat). However, both of these compounds were significantly weaker inhibitors of HDAC3 than SAHA and thus more selective, albeit to a limited extent. Five compounds exhibited activity on human breast carcinoma and/or urothelial carcinoma cell lines. The present result suggests potential mechanistic and chemical approaches for developing selective HDAC4 modulators.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Modeling and simulations of hydrodynamic shocks in a plasma flowing across randomized ICF scale laser beams

High-energy laser beams interacting with flowing plasmas can produce a plasma response that leads to deflection of the beam, beam bending. Such beams have usually a speckle structure generated by optical smoothing techniques that reduce the spatial and temporal coherence in the laser field pattern. The cumulative plasma response from laser speckles slows down the velocity of the incoming flow by momentum conservation. For slightly super-sonic flow the cumulative plasma response to the ponderomotive force exerted by the beam speckle ensemble is the strongest, such that slowing down the flow to subsonic velocities leads eventually to the generation of a shock around the cross section of the beam. This scenario has been predicted theoretically and is confirmed here by our hydrodynamic simulations in two dimensions with speckled beams and in one dimension with a reduced model. The conditions of shock generation are given in terms of the ponderomotive pressure, speckle size and the flow velocity. The nonlinear properties of the shocks are analyzed using Rankine–Hugoniot relations. According to linear theory, temporally smoothed laser beams exhibit a higher threshold for shock generation. Numerical simulations with beams that are smoothed by spectral dispersion compare well with the linear theory results, diverging from those produced by beams with only a random phase plates in the nonlinear regime. The conditions necessary for shock generation and their effects on the laser plasma coupling in inertial confinement fusion (ICF) experiments are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Impacts of stochastic coalescence variability on warm rain initiation using Lagrangian microphysics in box and large-eddy simulations

Various coalescence methods for Lagrangian microphysics schemes are tested in box and large-eddy simulation (LES) models, including the stochastic all-or-nothing super-droplet method (SDM) and a version of SDM (dSDM) that applies a fractional approach similar to the average impact method. In LES, variabilities driven by microphysics and by flow realizations are separated using the “piggybacking” technique. Rain initiation averaged over many realizations of the box model is delayed and rain variability increases as the number of super-drops per collision volume ( N SD ) is decreased using SDM. In contrast, rain initiation time using SDM in LES is insensitive to N SD for 32 ≤ N SD ≤ 512. This is explained through the interaction between LES grid boxes, each acting as a separate collision volume. Variability across the ensemble of LES collision volumes using SDM results in rain quickly initiating in some of the LES grid cells at low N SD and leading to a similar overall timing of rain initiation from the cloud compared to simulations with high N SD . There is a ~20% decrease in the total rain mass and mean rain flux as N SD is increased from 32 to 256, with little additional change as N SD is increased from 256 to 512. The fractional coalescence approach in dSDM leads to reduced microphysical variability and a 15-18 min delay in rain initiation compared to SDM. Finally, an additional LES ensemble with microphysical variability feeding back to the dynamics shows that flow variability dominates the impact of microphysical variability on rain properties. Thus, flow variability must be constrained to isolate impacts of microphysical variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Determining the Ensemble N -Representability of Reduced Density Matrices

The N-representability problem for reduced density matrices remains a fundamental challenge in electronic structure theory. Following our previous work that employs a unitary-evolution algorithm based on an adaptive derivative-assembled pseudo-Trotter variational quantum algorithm to probe pure-state N-representability of reduced density matrices [J. Chem. Theory Comput. 2024, 20, 9968], in this work we propose a practical framework for determining the ensemble N-representability of a p-body matrix. This is accomplished using a purification strategy that embeds an ensemble state into a pure state defined on an extended Hilbert space, such that the reduced density matrices of the purified state reproduce those of the original ensemble. By iteratively applying variational unitaries to an initial purified state, the proposed algorithm minimizes the Hilbert-Schmidt distance between its p-body reduced density matrix and a specified target p-body matrix, which serves as a measure of the N-representability of the target. This methodology facilitates both error correction of defective ensemble reduced density matrices and quantum-state reconstruction on a quantum computer, offering a route for density-matrix refinement. We validate the algorithm with numerical simulations on systems of two, three, and four electrons in both simple models as well as molecular systems at finite temperature, demonstrating its robustness.

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

Process‐Oriented Calibration of a Turbulence Scheme in the DOE's Global Storm‐Resolving Model Using Machine Learning

A process‐oriented calibration framework is developed for the Simplified Higher‐Order Closure (SHOC) turbulence scheme in DOE's Simple Cloud Resolving E3SM Atmospheric Model (SCREAM). This framework leverages machine learning surrogates and observational constraints to efficiently calibrate SHOC adjustable parameters across two convective regimes: clear‐sky dry convective boundary layer and fair‐weather shallow cumulus clouds from ARM observations. We use perturbed‐parameter ensembles of a doubly periodic version of SCREAM to train surrogates and apply Markov Chain Monte Carlo sampling guided by cost functions based on benchmarking large‐eddy simulations and observations to identify optimized parameter sets that perform well in both regimes. The calibrated SHOC parameters substantially improve boundary‐layer turbulence and cloud boundaries, and modeled cloud fraction and radiative effects align better with observations than the default. These results demonstrate that combining multiple process‐specific convective regimes with machine‐learning surrogates can reduce parametric uncertainties and yield a model more faithful to cloud–turbulence interactions.

58 GEOSCIENCES↗
DOE OSTI2025

A Polynomial-Time Classical Algorithm for Noisy Quantum Circuits

We provide a polynomial-time classical algorithm for noisy quantum circuits. The algorithm computes the expectation value of any observable for any circuit, with a small average error over input states drawn from an ensemble (e.g., the computational basis). Our approach is based upon the intuition that noise exponentially damps nonlocal correlations relative to local correlations. This enables one to classically simulate a noisy quantum circuit by keeping track of only the dynamics of local quantum information. Our algorithm also enables sampling from the output distribution of a circuit in quasipolynomial time, so long as the distribution anticoncentrates. A number of implications are discussed, including a fundamental limit on the efficacy of noise mitigation strategies: For constant noise rates, any quantum circuit for which error mitigation succeeds in polynomial-time on most input states can also be classically simulated in polynomial-time on most input states. Our algorithms scale exponentially in the inverse noise rate, which is fundamental and makes them impractical for current quantum devices.

decoherence↗
DOE OSTI2026

Decadal Spiciness Variability in the Subtropical-Tropical Pacific in the CESM2 Large Ensemble

Tropical Pacific decadal variations impact weather and climate around the world and are also connected to variations in the global warming trend. The mechanisms driving these long-term modulations, particularly the role of subsurface ocean dynamics, are still debated. Here, we investigate the dynamics of spiciness (density-compensated temperature and salinity) anomalies in the tropical and subtropical Pacific, which are hypothesized as a possible driving mechanism of decadal climate variability. Based on the analysis of 100 realizations from the Community Earth System Model Version 2 Large Ensemble (CESM2-LE), we demonstrate a coupling between the subtropics and the equatorial Pacific by propagating spiciness anomalies at decadal time scales. The CESM2-LE simulates spiciness variability along a subduction path from the subtropics to the equator with frequency spectra that show the highest power at low frequencies and a power decay proportional to a −4 slope for frequencies greater than 0.01 cycles per months, corresponding to periods smaller than ∼8.5 years. Signals that originate in the Southern Hemisphere (SH) dominate and arrive with a larger magnitude at the equator compared to spiciness anomalies from the Northern Hemisphere (NH). Spiciness anomalies from the SH have shorter propagation times and are strengthened along their pathway as stochastic wind stress curl forcing generates anomalous baroclinic ocean pressure gradients. These pressure gradients generate spiciness anomalies via anomalous advection across climatological spiciness gradients in the SH. We conclude that the observed spiciness variance at decadal time scales is consistent with a forcing by stochastic wind variations that are low-pass filtered by ocean dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Probabilistic Forward Modeling of Galaxy Catalogs with Normalizing Flows

Abstract Evaluating the accuracy and calibration of the redshift posteriors produced by photometric redshift (photo- z ) estimators is vital for enabling precision cosmology and extragalactic astrophysics with modern wide-field photometric surveys. Evaluating photo- z posteriors on a per-galaxy basis is difficult, however, as real galaxies have a true redshift but not a true redshift posterior. We introduce PZFlow, a Python package for the probabilistic forward modeling of galaxy catalogs with normalizing flows. For catalogs simulated with PZFlow, there is a natural notion of “true” redshift posteriors that can be used for photo- z validation. We use PZFlow to simulate a photometric galaxy catalog where each galaxy has a redshift, noisy photometry, shape information, and a true redshift posterior. We also demonstrate the use of an ensemble of normalizing flows for photo- z estimation. We discuss how PZFlow will be used to validate the photo- z estimation pipeline of the Dark Energy Science Collaboration, and the wider applicability of PZFlow for statistical modeling of any tabular data.

Astronomy & Astrophysics↗
DOE OSTI2025

Deforestation Increases Vegetation Vulnerability to Drought Across Biomes

Land use and land cover changes have altered terrestrial ecosystem carbon storage, but their impacts on ecosystem sensitivity to drought and temperature fluctuations have not been evaluated spatially over the globe. We estimate drought and temperature sensitivities of ecosystems using vegetation greenness from satellite observations and vegetation biomass from dynamic global vegetation model (DGVM) simulations. Using a space-for-time substitution with satellite data, we first illustrate the effects of vegetation cover changes on drought and temperature sensitivity and compare them with the effects estimated from DGVMs. We also compare simulations forced by scenarios with and without land cover changes to estimate the historical land cover change effects. Satellite data and vegetation models both show that converting forests to grasslands results in a more negative or decreased positive sensitivity of vegetation greenness or biomass to drought. Significant variability exists among models for other types of land cover transitions. We identify substantial effects of historical land cover changes on drought sensitivity from model simulations with a generally positive direction globally. Deforestation can lead to either an increased negative sensitivity, as drought-tolerant forests are replaced by grasslands based on model ensemble mean, or a decreased negative sensitivity, since forests under current land cover are predicted to exhibit greater drought resistance compared to those under pre-industrial land cover. Overall, our findings emphasize the critical role of forests in maintaining ecosystem stability and resistance to drought and temperature fluctuations, thereby implying their importance in stabilizing the carbon stock under increasingly extreme climate conditions.

Xiao, Chenwei [Max Planck Institute for Biogeochem↗
DOE OSTI2025

$B\rightarrow D^{(*)}$ decays from $N_f=2+1+1$ highly improved staggered quarks and clover $b$-quark in the Fermilab interpretation

We present an update on the analysis of semileptonic $B\rightarrow D^{(*)}$ decays at non-zero recoil. Our computation employs $2 + 1 + 1$ FNAL-MILC ensembles with highly improved staggered quark (HISQ) action for sea and light valence quarks, while the bottom quark is treated using the clover action in the Fermilab interpretation. Simulations are performed across several lattice spacings, ranging approximately from $\sim 0.15$ fm to $\sim 0.06$ fm, and for various quark masses. We will present an overview of the analysis and show some preliminary results for the form factors.

DeTar, Carleton [Utah U.] (ORCID:0000000202166771)↗
DOE OSTI2025

Evaluating Mean State Cloud Properties in the Simple Cloud‐Resolving E3SM Atmosphere Model (SCREAM)

Accurately simulating clouds remains a key challenge in global climate models, primarily because cloud formation involves sub-grid processes that are parameterized and crudely represented in models. This study examines the performance of DOE's Simple Cloud-Resolving Energy Exascale Earth System (E3SM) Atmosphere Model (SCREAM) in simulating cloud properties and their spatio-temporal distribution by comparing against satellite observations. Two horizontal resolutions of SCREAM (3 and 12 km) are examined, and both depict a realistic spatial structure of mean-state cloud cover but underestimate its global mean magnitude. SCREAM 3 km reasonably reproduces the distribution of mean-state cloud properties across various cloud optical thickness and cloud-top pressure regimes, with performance comparable to CMIP5 and CMIP6 ensemble and marginally outperforming SCREAM 12 km. Still, SCREAM 3 km tends to underpredict low clouds and optically thin clouds, highlighting the need for continued improvement in representing unresolved processes. This study provides a basis for confidence in the representation of clouds in SCREAM, as simulating mean-state clouds is a necessary prerequisite for trusting its cloud responses to changes in aerosols and greenhouse gases.

Geosciences↗
DOE OSTI2025

A Parametric, Data-Driven, Non-Intrusive Reduced-Order Model Framework for Crystal Plasticity Simulations of Voids

The influence of the internal structure at micrometer length scales on the deformation of polycrystalline materials can be effectively captured using crystal plasticity finite element methods (CPFEM). However, the complexity and nonlinearity of the deformation equations CPFEM solves demand significant computational power and resources to achieve accurate predictions, limiting its broader application. To address this challenge, we have identified a reduced-order representation of the complex data in order to establish a computationally efficient reduced-order models (ROM) and drastically reduce the computational expense of CPFEM. Specifically, in this work, we developed a parametric, data-driven, and non-intrusive ROM framework for CPFEM using proper orthogonal decomposition (POD) and sparse variational Gaussian process (SVGP) regression for single-crystal microstructures under tensile loading conditions. The developed protocol enables one to compress field into a latent/low-dimensional space described by principal component analysis (PCA) via the singular value decomposition (SVD) algorithm. As a result, the high-dimensional data are reduced to a significantly smaller amount of dimensions with POD bases and POD coefficients. Furthermore, we deployed an ensemble of SVGPs—extended from the classical Gaussian process (GP) regression for scalability and handling big data—in a massively parallel manner to train and predict latent POD coefficients using known POD bases from a set of previously obtained simulations results. Lastly, using the predicted POD coefficients, we reconstructed the full-field results and showed reasonable agreement compared with the true values obtained from running CPFEM. The developed framework is validated with a set of CPFEM simulations of a single embedded void in single-crystal aluminum alloy. While the framework is broadly applicable, this work specifically focuses on single-crystal microstructures, a single load case (e.g., tensile), and a specific void geometry (spherical).

Anisotropy↗
DOE OSTI2025

How Representative Are Uncrewed Aircraft System Measurements of the Convective Boundary Layer?

Abstract Uncrewed aircraft systems (UAS) demonstrate significant potential for filling data gaps in the atmospheric boundary layer. However, the extent to which UAS observations—typically vertical profiles taken over 15 min—are representative of the boundary layer as a whole remains poorly characterized. Using large eddy simulations (LES) of the daytime convective boundary layer (CBL), we quantify random errors in UAS measurements that occur due to insufficient statistical convergence of the time average to the true ensemble mean. Random errors in first‐order moments increase as the CBL becomes increasingly unstable, and are largest near the surface for most quantities. Errors are on the order of 2–6 m for wind speed, 15–60 for wind direction, 0.2–3 K for potential temperature, and 0.1–1 g for specific humidity, with errors in turbulent fluxes on the order of 50%–100%. Sampling strategies that mitigate random errors are discussed in light of our results.

Greene, Brian R. [Now at Verisk Extreme Event Solu↗
DOE OSTI2024

HSW-V v1.0: localized injections of interactive volcanic aerosols and their climate impacts in a simple general circulation model

Abstract. A new set of standalone parameterizations is presented for simulating the injection, evolution, and radiative forcing by stratospheric volcanic aerosols against an idealized Held–Suarez–Williamson (HSW) atmospheric background in the Energy Exascale Earth System Model version 2 (E3SMv2). In this model configuration (HSW with enabled volcanism, HSW-V), sulfur dioxide (SO2) and ash are injected into the atmosphere with a specified profile in the vertical, and they proceed to follow a simple exponential decay. The SO2 decay is modeled as a perfect conversion to a long-living sulfate aerosol which persists in the stratosphere. All three species are implemented as tracers in the model framework and are transported by the dynamical core's advection algorithm. The aerosols contribute simultaneously to local heating of the stratosphere and cooling of the surface by a simple plane-parallel Beer–Lambert law applied on two zonally symmetric radiation broadbands in the longwave and shortwave ranges. It is shown that the implementation parameters can be tuned to produce realistic temperature anomaly signatures of large volcanic events. In particular, results are shown for an ensemble of runs that mimic the volcanic eruption of Mt. Pinatubo in 1991. The design requires no coupling to microphysical subgrid-scale parameterizations and thus approaches the computational affordability of prescribed aerosol forcing strategies. The idealized simulations contain a single isolated volcanic event against a statistically uniform climate, where no background aerosols or other sources of externally forced variability are present. HSW-V represents a simpler-to-understand tool for the development of climate source-to-impact attribution methods.

Hollowed, Joseph P. (ORCID:0000000286581672)↗
DOE OSTI2025

Stochastic Ensemble Generation for Improved Characterization of Representing Geologic Variability in a Reservoir: IBDP Case Study for SMART Initiative

This document is a poster covering the findings from activities on training data generation, specifically geologic ensemble generation. The generated geologic realizations captured the range of possible permeability distributions of the subsurface at the Illinois Basin - Decatur Project (IBDP) site, based on available well log variabilities. The percentages of reservoirs and baffles in the injection zone and a truncation of baffle permeability led to more variance in the simulations. This will be used to build forward modeling, history matching, and optimization workflows. The geologic realizations were also ranked according to dynamic measures of hydraulic diffusivity, and simulations confirm a greater contrast between the reservoir and the baffles during injection.

stochastic ensemble generation↗