Search NASA⌕ Search

SEARCH · Search NASA

Results for “Large ensemble”

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

A Pacific Tropical Decadal Variability Challenge for Climate Models

Understanding and forecasting Tropical Pacific Decadal-scale Variability (TPDV) strongly rely on climate model simulations. Using a Linear Inverse Modeling (LIM) diagnostic approach, we reveal Coupled Model Intercomparison Project Phase 6 models have significant challenges in reproducing the spatial structure and dominant mechanisms of TPDV. Specifically, while the models' ensemble mean pattern of TPDV resembles that of observations, the spread across models is very large and most models show significant differences from observations. In observations, removing the coupling between extratropics and tropics reduces TPDV by ~60%–70%, and removing the tropical thermocline variability makes the central tropical Pacific a key center of action for TPDV and El Niño Southern Oscillation variability. These characteristics are only confirmed in a subset of models. Differences between observations and simulations are outside the range of natural internal TPDV noise and pose important questions regarding our ability to model the impacts of natural internal low-frequency variability superimposed on long-term climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Finite-size effects in lead scandium tantalate relaxor thin films

Large electromechanical effects in relaxor ferroelectrics are generally attributed to the collective response of an ensemble of correlated, nanometer-sized polar structures induced by chemical and charge disorder. Here, we study finite-size effects on such polar order (i.e., how it evolves when sample dimensions approach the polarization correlation length) in 7–70-nm-thick films of the relaxor ferroelectric PbSc 0.5 Ta 0.5 O 3 . Temperature-dependent polarization studies reveal a linear suppression of the polarization and nonlinearity associated with relaxor order as the film thickness decreases to ≈ 30 nm. Below this thickness, however, the suppression rapidly accelerates, and polarization is completely absent by film thicknesses of ≈ 7 nm, despite the continued observation of a broad peak in dielectric permittivity and frequency dispersion. Diffuse-scattering measurements reveal the diffuse-scattering symmetry, and analysis suggests the films have a polarization correlation length of ≈ 23 nm. Finally, taken together, it is apparent that reduction of sample size and the resulting distribution of polar structures drive suppression and eventual quenching of the electrical response of relaxors, which may be attributed to increasing dipole-dipole and dipole-interface interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Volcanic forcing degrades multiyear-to-decadal prediction skill in the tropical Pacific

Volcanic aerosol forcing can affect global climate, but its role in climate prediction remains poorly understood. We isolate the impact of volcanic eruptions on multiyear-to-decadal climate prediction skill by comparing two suites of initialized decadal hindcasts conducted with and without historical volcanic forcing. Unexpectedly, the inclusion of volcanic forcing in the prediction system significantly degrades the forecast skill of detrended multiyear-to-decadal sea surface temperature (SST) variability in the central-eastern tropical Pacific. The ensemble mean hindcasts produce multiyear-to-decadal tropical Pacific SST cooling in response to large tropical volcanic eruptions through thermodynamic and El Niño–Southern Oscillation (ENSO)–like dynamic processes. However, in observations, these eruptions coincided with tropical Pacific warming, which is well predicted by the no-volcano hindcasts and, hence, is likely related to internal climate variability. Improved model representation of volcanic response and its interaction with internal climate variability is required to advance prediction of tropical Pacific decadal variability and associated global impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Shallow Cumulus Entrainment Dynamics in a Sheared Environment

Abstract Vertical wind shear has long been known to tilt convective towers and reduce thermal ascent rates. The purpose of this study is to better understand the physical mechanisms responsible for reduced ascent rates in shallow convection. In particular, the study focuses on cloud-edge mass flux to assess how shear impacts mass-flux profiles of both the ensemble and individual clouds of various depths. A compositing algorithm is used to distill large-eddy simulation (LES) output to focus on up- and down-shear cloud edges that are not influenced by complex cloud geometry or nearby clouds. A direct entrainment algorithm is used to estimate the mass flux through the cloud surface. We find that the dynamics on the up- and down-shear sides are fundamentally different, with the entrainment of environmental momentum and dilution of buoyancy being primarily responsible for the reduced down-shear ascent rates. Direct estimates of fluid flow through the cloud interface indicate a counter-shear organized flow pattern that entrains on the down-shear side and detrains on the up-shear side, resulting from the subcloud shear being lifted into the cloud layer by the updraft. In spite of organized regions of entrainment and detrainment, the overall net lateral mass flux remains unchanged with respect to the no shear run, with weak detrainment present throughout cloud depth.

Meteorology & Atmospheric Sciences↗
DOE OSTI2021

Anthropogenic Influences on Tornadic Storms

The impact of climate change on severe storms and tornadoes remains uncertain, largely owing to inconsistencies in observational data and limitations of climate models. Here, we performed ensembles of convection-permitting climate model simulations to examine how three tornadic storms would change if similar events were to occur in pre-industrial and future climates. The choice of events includes winter, nocturnal, and spring tornadic storms to provide insight into how the timing and seasonality of storms may affect their response to climate change. Updraft helicity (UH), convective available potential energy (CAPE), storm relative helicity (SRH), and convective inhibition (CIN) were used to determine the favorability for the three tornadic storm events in the different climate states. We found that from the pre-industrial to present, the potential for tornadic storms decreased in the winter event and increased in the nocturnal and spring events. With future climate change, the potential for tornadic storms increased in the winter and nocturnal events in association with increased CAPE, and decreased in the spring event despite greater CAPE.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

An efficient hybrid downscaling framework to estimate high-resolution river hydrodynamics

Flow depth and velocity are the most important hydrodynamic variables that govern various river functions, including water resources, navigation, sediment transport, and biogeochemical cycling. Existing high-resolution flow depth simulations rely on either computationally expensive river hydrodynamic models (RHMs) or data-driven models with formidable training costs, whereas data-driven modeling of flow velocity has rarely been explored. Here, using the hybrid Low-fidelity, Spatial analysis, and Gaussian process learning (LSG) model, we developed a downscaling approach to construct high-resolution flow depth and velocity from a two-dimensional (2-D) RHM simulation at coarse resolution. The LSG models were trained and tested in an urban watershed in Houston using two different hurricane-driven flood events. The high-resolution (as fine as 30 m resolution) and low-resolution (mostly 1000 m resolution) meshes include 664 724 and 14 536 grid cells, respectively. The results showed that through downscaling, the simulation errors were reduced to less than one-fourth and one-third of the errors of the low-resolution 2-D RHM for flow depth and velocity, respectively. Our analysis further revealed that the dominant uncertainty sources of the downscaled hydrodynamics are different, with flow velocity dominated by the dimensionality reduction error, which we reduced by using a regionalized training procedure. The downscaling approach achieves an 84-fold acceleration in computational time compared to the high-resolution 2-D RHM, making high-fidelity ensemble flood modeling feasible. More importantly, the developed method provides an opportunity to couple large-scale hydrodynamical processes with local physical, chemical, and biological processes in river models.

Tan, Zeli [Pacific Northwest National Laboratory (↗
DOE OSTI2023

Large-scale molecular dynamics simulations of bubble collapse in water: Effects of system size, water model, and nitrogen

Molecular dynamics simulations in the microcanonical ensemble are performed to study the collapse of a bubble in liquid water using the single-site mW and the four-site TIP4P/2005 water models. To study system size effects, simulations for pure water systems are performed using periodically replicated simulation boxes with linear dimensions, L, ranging from 32 to 512 nm with the largest systems containing 8.7 × 10 6 and 4.5 × 10 9 molecules for the TIP4P/2005 and mW water models, respectively. The computationally more efficient mW water model allows us to reach converging behavior when the bubble dynamics results are plotted in reduced units, and the limiting behavior can be obtained through linear extrapolation in L –1 . Qualitative differences are observed between simulations with the mW and TIP4P/2005 water models, but they can be explained by the models’ differences in predicted viscosity and surface tension. Although bubble collapse occurs on time scales of only hundreds of picoseconds, the system sizes used here are sufficiently large to obtain bubble dynamics consistent with the Rayleigh–Plesset equation when using the models’ thermophysical properties as input. For the conditions explored here, extreme heating of the interfacial water molecules near the time of collapse is observed for the larger mW water systems (but the model underpredicts the viscosity), whereas heating is less pronounced for the TIP4P/2005 water systems because its larger viscosity contribution slows the collapse dynamics. The presence of nitrogen within the bubble only starts to affect bubble dynamics near the very end of the initial collapse, leading to an incomplete collapse and strong rebound for the mW water model. Although nitrogen is non-condensable at 300 K, it becomes highly compressed and reaches a liquid-like density near the collapse point. We find that the dissolution of nitrogen is much slower than the movement of the collapsing water front, and the re-expansion of the dense nitrogen droplet gives rise to bubble rebound. The incompatibility of the collapse and dissolution time scales should be considered for continuum-scale modeling of bubble dynamics. Finally, we also confirm that the diffusion coefficient for dissolved nitrogen is insensitive to pressure as the liquid transitions from a compressed to a stretched state.

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

Mesospheric temperature and circulation response to the Hunga Tonga-Hunga-Ha'apai volcanic eruption

The Hunga Tonga Hunga-Ha'apai (HTHH) volcanic eruption on 15 January 2022 injected water vapor and SO 2 into the stratosphere. Several months after the eruption, significantly stronger westerlies, and a weaker Brewer-Dobson circulation developed in the stratosphere of the Southern Hemisphere and were accompanied by unprecedented temperature anomalies in the stratosphere and mesosphere. In August 2022, the Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) satellite instrument observed record-breaking temperature anomalies in the stratosphere and mesosphere that alternate signs with altitude. Ensemble simulations carried out with the Whole Atmosphere Community Climate Model (WACCM6) indicate that the strengthening of the stratospheric westerlies explains the mesospheric temperature changes. The stronger westerlies cause stronger westward gravity wave drag in the mesosphere. Although the enhanced gravity wave drag is partly balanced by a weakening of planetary wave forcing, the net result is an acceleration of the mesospheric mean meridional circulation. The stronger mesospheric circulation, in turn, plays a dominant role in driving the changes in mesospheric temperatures. This study highlights the impact of large volcanic eruptions on middle atmospheric dynamics and provides insight into their long-term effects in the mesosphere. On the other hand, we could not discern a clear mechanism for the observed changes in stratospheric circulation. In fact, an examination of the WACCM ensemble reveals that not every member reproduces the large changes observed by SABER. We conclude that there is a stochastic component to the stratospheric response to the HTHH eruption.

58 GEOSCIENCES↗
DOE OSTI2019

CRAGE enables rapid activation of biosynthetic gene clusters in undomesticated bacteria

Advances in the prediction of biosynthetic gene clusters (BGCs) from genomic and metagenomic data have made the large-scale identification of candidate secondary metabolite biosynthesis pathways routine. However, most predicted BGCs cannot be successfully expressed in the limited range of expression hosts currently available, which prevents the identification and characterization of their respective products. Here we show that parallelized use of a phylogenetically and physiologically diverse panel of purpose-engineered chassis strains offers a versatile strategy for the rapid activation of BGCs. We domesticated an ensemble of 25 diverse ?-proteobacteria from 11 different genera by integrating a unified landing pad that enables efficient introduction of large BGCs into their genome. We deployed this chassis-independent recombinase-assisted genome engineering (CRAGE) strategy to express several previously posited but experimentally elusive non-ribosomal peptide synthase (NRPS) and NRPSpolyketide synthase (PKS) hybrid BGCs from Photorhabdus luminescens. Successful activation of 6 BGCs identified 21 products, including several with likely roles in host-microbe interactions. Taken together, we illustrate how the regulatory and physiological diversity naturally present across bacteria can be harnessed for the exploration of novel bioproducts and provide a sizable panel of domesticated chassis strains to enable this strategy.

applied microbiology, genetic engineering, industr↗
DOE OSTI2022

Shock ignition and deflagration growth in plastic-bonded TATB (1, 3, 5-trinitro-2, 4, 6-triaminobenzene) microstructures

TATB (1,3,5-triamino-2,4,6-trinitrobenzene) plastic-bonded explosives (PBX) were shocked with laser-launched flyer plates. The spectral radiance of the emitted light from a small portion of the microstructure (a “microenvironment”) containing a small number of TATB particles with an estimated mass of 150 ng was measured every 0.8 ns from 1 ns to 200 μs and was analyzed to give radiance and time-dependent graybody temperatures. By fabricating an array with 186 PBX charges, we could obtain ≥15 shots at each of 12 velocities between 1.8 and 4.7 km/s. We found that every microenvironment generated a unique radiance fingerprint. Some of these microenvironments were much more reactive than average. The radiance has two peaks around 20 ns and 5 μs, associated with shock ignition and deflagration growth. In our interpretation, the shock creates an ensemble of hot spots of various sizes and temperatures. Of those hot spots that ignite, only a small portion, at about 2200 K, was large enough and hot enough to survive long enough (>100 ns) to ignite individual TATB particles, leading to deflagration. Integrating various time intervals of the radiance can quantify the strength of the shock–PBX interaction, and the decay and growth of the hot spot ensemble and the deflagration.

Akhtar, Meysam (ORCID:0000000243897837)↗
DOE OSTI2024

Approaching periodic systems in ensemble density functional theory via finite one-dimensional models

Ensemble density functional theory (EDFT) is a generalization of ground-state DFT, which is based on an exact formal theory of finite collections of a system's ground and excited states. EDFT in various forms has been shown to improve the accuracy of calculated energy level differences in isolated model systems, atoms, and molecules, but it is not yet clear how EDFT could be used to calculate band gaps for periodic systems. We extend the application of EDFT toward periodic systems by estimating the thermodynamic limit with increasingly large finite one-dimensional 'particle in a box' systems, which approach the uniform electron gas (UEG). Using ensemble-generalized Hartree and local spin density approximation exchange-correlation functionals, we find that corrections go to zero in the infinite limit, as expected for a metallic system. However, there is a correction to the effective mass, with results comparable to other calculations on 1D, 2D, and 3D UEGs, which indicates promise for non-trivial results from EDFT on periodic systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Multifidelity validation of digital surrogates using variable-density turbulent mixing models

In this study, ensembles of experimental data are presented and utilized to compare and validate two models used in the simulation of variable-density (Atwood = 0.22), compressible turbulent mixing. Though models of this kind (Reynolds averaged NavierStokes and large-eddy simulations) have been validated extensively with more canonical flows in previous studies, here the present approach offers novelty in the complexity of the geometry, the ensemble-based validation, and the uniformity of the computational framework on which the models are tested. Moreover, all experimental and computational tasks were completed by the authors which has led to a tightly coupled experimental configuration with its “digital twin.” The experimental divergent-shock-tube facility and its data acquisition methods are described and replicated in simulation space. A 2D Euler model which neglects the turbulent mixing at the interface is optimized to experimental data using a Gaussian process. This model then serves as the basis for both the 2D RANS and 3D LES studies that make comparisons to the mixing-layer data from the experiment. A relatively simple RANS model is shown to produce good agreement with experimental data only at late flow development times. The LES ensembles generally show good agreement with experimental data but display sensitivity to the characterization of initial conditions. Resolution-dependent behavior is also observed for certain higher-order statistics of interest. Overall, the LES model successfully captures the effects of divergent geometry, compressibility, and combined nonlinear instabilities inherent to the problem. The successful prediction of mixing width and its growth rate highlight the existence of three distinct regimes in the development of the instability, each with similarities to previously studied instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Finite-Temperature Structures of Supported Subnanometer Catalysts Inferred via Statistical Learning and Genetic Algorithm-Based Optimization

Single-atom catalysts (SACs) minimize noble metal utilization and can alter the activity and selectivity of supported metal nanoparticles. However, the morphology of active centers, including single atoms and subnanometer clusters of a few atoms, remains elusive due to experimental challenges. The computational cost to describe numerous cluster shapes and sizes makes direct first-principles calculations impractical. We present a computational framework to enable structure determination for single-atom and subnanometer cluster catalysts. As a case study, we obtained the low energy structures of Pd n (n = 1-21) clusters supported on CeO 2 (111), which are critical components of automobile three-way catalysts. Trained on density functional theory data, a three-dimensional cluster expansion is established using statistical learning to describe the Hamiltonian and predict energies of supported Pdn clusters of any structure. Low energy stable and metastable structures are identified using a Metropolis Monte Carlo-based genetic algorithm in the canonical ensemble at 300 K. We observe that supported single atoms sinter to form bilayer clusters and large cluster isomers share similarities in both shape and energy, and elucidate the significance of the support and microstructure on cluster stability. We discovered a simple surrogate structure-energy model, where the energy per atom scales with the square root of the average first coordination number, which can be used to estimate energies and compare the stability of clusters. Our framework, applicable to any metal/support system, fills an important methodological gap to predict the stability of supported metal catalysts in the subnanometer regime.

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

Vegetation biogeography is a main source of uncertainty in modelling the land carbon cycle

The terrestrial biosphere exchanges a large amount of CO 2 with the atmosphere through photosynthesis and respiration, determining the magnitude of land carbon sink and consequently influencing the rate of global warming. The magnitudes of global photosynthesis and respiration, however, vary widely across models (100-200 PgC/year), constituting a key and persistent source of uncertainty in carbon cycle and climate modelling. Here, we argue that the uncertainty in the land carbon cycle modelling is largely attributable to the uncertainty in biogeography – the distribution of plant functional types (PFTs). Using an ensemble of dynamic global vegetation models (DGVMs), we find a strong dependence of total photosynthesis on total area for each PFT. The dependence allows us to reduce the spread of land carbon cycle estimates by ~75% using remote sensing-based PFT maps. We further find that 56 ± 21% of climate-driven changes in global photosynthesis modelled by DGVMs are caused by changes in PFT distribution in the last two decades. Our study identifies vegetation biogeography as a main controlling factor of uncertainty in land carbon cycle modelling and highlights the importance of biogeography-climate interactions in carbon cycle and climate studies.

Zhao, Ruiying [National Univ. of Singapore (Singap↗
DOE OSTI2021

Cosmic-ray transport and gamma-ray emission in M31

Here we study the possibility that an extended cosmic-ray leptonic and/or hadronic halo is at the origin of the large-scale gamma-ray emission detected from the Andromeda galaxy (M31). We consider a broad ensemble of nonhomogeneous diffusion scenarios and of cosmic-ray injection sources. We find that cosmic-ray electrons and protons could be, and very likely are, responsible for part, or all, of the gamma-ray emission from M31, including out to more than 100 kpc from the center of the galaxy. We also simulate possible emission from pulsars in M31 and consider the effect of regions of highly inefficient diffusion around cosmic-ray acceleration sites, as suggested by recent TeV halo observations with Cherenkov telescopes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Tree-based algorithms for weakly supervised anomaly detection

Weakly supervised methods have emerged as a powerful tool for model-agnostic anomaly detection at the Large Hadron Collider (LHC). While these methods have shown remarkable performance on specific signatures such as dijet resonances, their application in a more model-agnostic manner requires dealing with a larger number of potentially noisy input features. In this paper, we show that using boosted decision trees as classifiers in weakly supervised anomaly detection gives superior performance compared to deep neural networks. Boosted decision trees are well known for their effectiveness in tabular data analysis. Our results show that they not only offer significantly faster training and evaluation times, but they are also robust to a large number of noisy input features. By using advanced gradient boosted decision trees in combination with ensembling techniques and an extended set of features, we significantly improve the performance of weakly supervised methods for anomaly detection at the LHC. This advance is a crucial step toward a more model-agnostic search for new physics. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2020

The vertical structure of convective mass-flux derived from modern radar systems: Data analysis in support of cumulus parametrization (Final Report)

The project delivered against all its major aims. We first developed a new calibration technique that enabled the construction of a long-term (17 years) calibrated radar dataset for the Darwin region. We documented both the technique and the data set in publications. We then developed a corresponding long-term dataset that characterises the large-scale state of the atmosphere. Equipped with both, we examined the relationship of key convective ensemble characteristics, such as cloud number and size, with the convective environment. In addition to the stated goals, we also applied the data set to i) estimate rainfall efficiency; ii) derive a new metric for convective organisation and iii) derive a longterm data set of convective mass-flux in the Darwin region and apply it to derive estimates of entrainment and detrainment rates in a convective cloud ensemble. Most of the research has been published or submitted for publication and will therefore only briefly summarised here, with links to or copies of the respective papers supplied.

58 GEOSCIENCES↗
DOE OSTI2022

Hybrid quantum-classical algorithms for approximate graph coloring

We show how to apply the recursive quantum approximate optimization algorithm (RQAOA) to MAX- k -CUT, the problem of finding an approximate k -vertex coloring of a graph. We compare this proposal to the best known classical and hybrid classical-quantum algorithms. First, we show that the standard (non-recursive) QAOA fails to solve this optimization problem for most regular bipartite graphs at any constant level p : the approximation ratio achieved by QAOA is hardly better than assigning colors to vertices at random. Second, we construct an efficient classical simulation algorithm which simulates level- 1 QAOA and level- 1 RQAOA for arbitrary graphs. In particular, these hybrid algorithms give rise to efficient classical algorithms, and no benefit arising from the use of quantum mechanics is to be expected. Nevertheless, they provide a suitable testbed for assessing the potential benefit of hybrid algorithm: We use the simulation algorithm to perform large-scale simulation of level- 1 QAOA and RQAOA with up to 300 qutrits applied to ensembles of randomly generated 3 -colorable constant-degree graphs. We find that level- 1 RQAOA is surprisingly competitive: for the ensembles considered, its approximation ratios are often higher than those achieved by the best known generic classical algorithm based on rounding an SDP relaxation. This suggests the intriguing possibility that higher-level RQAOA may be a potentially useful algorithm for NISQ devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗