SEARCH · Search NASA

Results for “ensemble dynamics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Molecular-Level Exploration of the Structure-Function Relations Underlying Interfacial Charge Transfer in the Subphthalocyanine/C 60 Organic Photovoltaic System

The arrangement of organic molecules at the donor-acceptor interface in an organic photovoltaic (OPV) cell can have a strong effect on the generation of charge carriers and thereby cell performance. In this paper, we report the molecular-level exploration of the ensemble of interfacial donor-acceptor pair geometries and the charge-transfer (CT) rates to which they give rise. Our approach combines molecular-dynamics simulations, electronic structure calculations, machine learning, and rate theory. This approach is applied to the boron subphthalocyanine chloride (donor) and C 60 (acceptor) OPV system. In this work we find that the interface is dominated by a previously unreported donor-acceptor pair edge geometry, which contributes significantly to device performance in a manner that depends on the initial conditions. Quantitative relations between the morphology and CT rates are established, which can be used to advance the design of more efficient OPV devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1984

Fracture, failure, and fragmentation

A general approach is required for describing matter of behavior when the failure is likely to involve growth and coalescence of a large number of fractures. Failures of this kind appear frequently in rapid dynamic processes, particularly in the formation of spall fragments. An approach to formulating constitutive relations that accounts for the opening, shear and growth of an ensemble of cracks is discussed. The approach accounts for plastic flow accompanying fragmentation. The resulting constitutive relations were incorporated into a Lagrangian computer program. A theoretical approach to coalescence is described. The simplest formulation uses a linear Liouville equation, with crack growth limited by the mean free path of cracks, assumed constant. This approach allows for an anisotropic distribution of cracks. An alternative approach in which the decrease of the mean free path with increasing crack size is accounted for, but the crack distribution is assumed isotropic is described. A reduction of the governing Liouville equation to an ordinary differential equation of third order is possible, and the result can be used to determine how mean free path decreases with increasing crack size.

Dienes, J. K.↗

NASA NTRS1985

A monthly forecast experiment: Preliminary report

An experiment on monthly forecasts with eight winter cases is being carried out, using a 1980 general circulation model (GCM), which incorporates a set of dubgrid-scale physics characterized by the Mellon-Yamada turbulence closure (hierarchy level 2.5), the Monin-Obukhov parameterization for the layer next to the ground surface, Manabe's cumulus parameterization, and the soil heat conduction. The sample cases adopted are for the month of January in the years from 1977 to 1983, which include the extraordinarily severe winter of 1977 and the most pronounced E1 Nino year of 1983. Each case is predicted by prescribing the climatologically normal sea surface temperature as the lower boundary conditions and by using an ensemble means of three individual integrations. These integrations start with three different initial conditions based on the Level data generated separately at the Geophysical Fluid Dynamics Laboratory, the National Meteorological Center, and the European Centre for Medium Range Weather Forecasts.

Miyakoda, K.↗

NASA NTRS1993

An X-ray method for detecting substructure in galaxy clusters - Application to Perseus, A2256, Centaurus, Coma, and Sersic 40/6

We use the moments of the X-ray surface brightness distribution to constrain the dynamical state of a galaxy cluster. Using X-ray observations from the Einstein Observatory IPC, we measure the first moment FM, the ellipsoidal orientation angle, and the axial ratio at a sequence of radii in the cluster. We argue that a significant variation in the image centroid FM as a function of radius is evidence for a nonequilibrium feature in the intracluster medium (ICM) density distribution. In simple terms, centroid shifts indicate that the center of mass of the ICM varies with radius. This variation is a tracer of continuing dynamical evolution. For each cluster, we evaluate the significance of variations in the centroid of the IPC image by computing the same statistics on an ensemble of simulated cluster images. In producing these simulated images we include X-ray point source emission, telescope vignetting, Poisson noise, and characteristics of the IPC. Application of this new method to five Abell clusters reveals that the core of each one has significant substructure. In addition, we find significant variations in the orientation angle and the axial ratio for several of the clusters.

Mohr, Joseph J.↗

NASA NTRS1999

Potential Predictability of the Monsoon Subclimate Systems

While El Nino/Southern Oscillation (ENSO) phenomenon can be predicted with some success using coupled oceanic-atmospheric models, the skill of predicting the tropical monsoons is low regardless of the methods applied. The low skill of monsoon prediction may be either because the monsoons are not defined appropriately or because they are not influenced significantly by boundary forcing. The latter characterizes the importance of internal dynamics in monsoon variability and leads to many eminent chaotic features of the monsoons. In this study, we analyze results from nine AMIP-type ensemble experiments with the NASA/GEOS-2 general circulation model to assess the potential predictability of the tropical climate system. We will focus on the variability and predictability of tropical monsoon rainfall on seasonal-to-interannual time scales. It is known that the tropical climate is more predictable than its extratropical counterpart. However, predictability is different from one climate subsystem to another within the tropics. It is important to understand the differences among these subsystems in order to increase our skill of seasonal-to-interannual prediction. We assess potential predictability by comparing the magnitude of internal and forced variances as defined by Harzallah and Sadourny (1995). The internal variance measures the spread among the various ensemble members. The forced part of rainfall variance is determined by the magnitude of the ensemble mean rainfall anomaly and by the degree of consistency of the results from the various experiments.

Yang, Song↗

NASA NTRS2023

Diverging Trends in Rain-on-Snow Over High Mountain Asia

Rain-on-snow (ROS) over snow-dominated regions such as High Mountain Asia (HMA) modulates snowmelt and runoff and is key contributor in influencing water availability and hazards (e.g., floods and landslides). We studied the trends in ROS in HMA over the past two decades from 2001 to 2018 using the land surface model Noah-MP driven by an ensemble precipitation data set. Our results show that changes in precipitation phase and rainfall are altering ROS. Because of the strong physical heterogeneity and atmospheric dynamics of HMA, ROS characteristics and trends are region-dependent and ROS occurs predominantly over the Indus, Ganges-Brahmaputra, and northwestern basins. In the Indus, ROS representing ∼5% of the annual precipitation and ∼20% of the annual snowmelt, has an increasing trend. This is contrary to the Ganges-Brahmaputra characterized by decreasing ROS trends, where it represents ∼11% of the annual precipitation and ∼60% of the annual snowmelt. In the northwestern basins, ROS has bidirectional trends due to elevation patterns and trends in rainfall, and it constitutes ∼5 to ∼10% of the annual precipitation. Increasing trends in ROS over Indus contribute to reducing the snowpack in late summer, with concerns of reduced water availability and increased groundwater exploitation. Similarly, because of its high amount and contribution to snowmelt, the decreasing ROS trends in the Ganges-Brahmaputra will have consequences of decreased recharge from the headwaters and exacerbated use of groundwater unless increasing trends in rainfall compensate for the decreasing snowmelt. These results provide new insights on ROS-driven changes in the hydrological cycle over HMA.

rain-on-snow↗

DOE OSTI2024

CMIP6-based Multi-model Streamflow Projections over the Conterminous US, Version 1.1

This dataset presents an ensemble of streamflow projections covering the conterminous United States (CONUS), developed to support the SECURE Water Act Section 9505 Assessment for the US Department of Energy (DOE) Water Power Technologies Office (WPTO). Multiple Coupled Models Intercomparison Project phase 6 (CMIP6) Global Climate Models (GCMs) were downscaled using either statistical (DBCCA) or dynamical (RegCM) downscaling methods, based on two meteorological reference datasets (Daymet and Livneh). Subsequently, the downscaled precipitation, temperature, and wind speed data were used to drive two calibrated hydrologic models (VIC and PRMS), with total runoff routed through the Routing Application for Parallel computatIon of Discharge (RAPID) routing model, producing an ensemble of streamflow projections across 2.7 million NHDPlusV2 stream reaches across the CONUS. Each ensemble member covers the 1980-2019 baseline and 2020-2059 near-future periods under the high-end (SSP585) emission scenario. Additionally, using only DBCCA and Daymet, the projections extend to the 2060-2099 far-future period and encompass three additional emission scenarios (SSP370, SSP245, and SSP126). This dataset is designed to support the SECURE Water Act Section 9505 Assessment for the US Department of Energy (DOE) Water Power Technologies Office (WPTO). For further details, refer to Kao et al. (2022), Rastogi et al. (2022), and Ghimire et al. (2023).

13 HYDRO ENERGY↗

NASA NTRS2013

Simultaneous Estimation of Model State Variables and Observation and Forecast Biases Using a Two-Stage Hybrid Kalman Filter

In this paper, we present a two-stage hybrid Kalman filter to estimate both observation and forecast bias in hydrologic models, in addition to state variables. The biases are estimated using the discrete Kalman filter, and the state variables using the ensemble Kalman filter. A key issue in this multi-component assimilation scheme is the exact partitioning of the difference between observation and forecasts into state, forecast bias and observation bias updates. Here, the error covariances of the forecast bias and the unbiased states are calculated as constant fractions of the biased state error covariance, and the observation bias error covariance is a function of the observation prediction error covariance. In a series of synthetic experiments, focusing on the assimilation of discharge into a rainfall-runoff model, it is shown that both static and dynamic observation and forecast biases can be successfully estimated. The results indicate a strong improvement in the estimation of the state variables and resulting discharge as opposed to the use of a bias-unaware ensemble Kalman filter. Furthermore, minimal code modification in existing data assimilation software is needed to implement the method. The results suggest that a better performance of data assimilation methods should be possible if both forecast and observation biases are taken into account.

Data Assimilation↗

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 OSTI2024

Insight into molecular basis and dynamics of full-length CRaf kinase in cellular signaling mechanisms

Raf kinases play key roles in signal transduction in cells for regulating proliferation, differentiation, and survival. Despite decades of research into functions and dynamics of Raf kinases with respect to other cytosolic proteins, understanding Raf kinases is limited by the lack of their full-length structures at the atomic resolution. Here, we present the first model of the full-length CRaf kinase obtained from artificial intelligence/machine learning algorithms with a converging ensemble of structures simulated by large-scale temperature replica exchange simulations. Our model is validated by comparing simulated structures with the latest cryo-EM structure detailing close contacts among three key domains and regions of the CRaf. Our simulations identify potentially new epitopes of intramolecule interactions within the CRaf and reveal a dynamical nature of CRaf kinases, in which the three domains can move back and forth relative to each other for regulatory dynamics. The dynamic conformations are then used in a docking algorithm to shed insight into the paradoxical effect caused by vemurafenib in comparison with a paradox breaker PLX7904. In this study, we propose a model of Raf-heterodimer/KRas-dimer as a signalosome based on the dynamics of the full-length CRaf.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Random Close Packing as a Dynamical Phase Transition

Sphere packing is an ancient problem. The densest packing is known to be a face-centered cubic (FCC) crystal, with space-filling fraction Φ FCC = π / √ 18 ≈ 0.74 . The densest “random packing,” random close packing (RCP), is yet ill defined, although many experiments and simulations agree on a value Φ RCP ≈ 0.64 . We introduce a simple absorbing-state model, biased random organization (BRO), which exhibits a Manna class dynamical phase transition between absorbing and active states that has as its densest critical point Φ cmax ≈ 0.64 ≈ Φ RCP and, like other Manna class models, is hyperuniform at criticality. The configurations we obtain from BRO appear to be structurally identical to RCP configurations from other protocols. This leads us to conjecture that the highest-density absorbing state for an isotropic biased random organization model produces an ensemble of configurations that characterizes the state conventionally known as RCP.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

High fidelity simulations of contaminant dispersion in an urban environment with comparison to magnetic resonance imaging measurements

The dispersion of a contaminant in an urban environment has the potential to impact a large population of people. In this work, a complex urban canopy flow based on the Oklahoma City downtown business district circa 2003 is studied using Magnetic Resonance Imaging (MRI) and high-fidelity Large Eddy Simulations (LES). MRI is a novel experimental technique that can provide high-resolution measurements in four dimensions (three spatial and temporal) for lab scale models. The experiments and simulations use the same geometry and boundary conditions providing a one-to-one comparison of the two methods. Results are presented on the time-averaged velocity and concentration fields, the temporal dynamics of the concentration plumes for a transient release, and a novel Cloud Identification Algorithm that can separate plumes produced by periodic contaminant releases used for ensemble averaging over many releases. The MRI and LES datasets both include millions of measurement voxels and the comparisons highlight the complex 3D nature of the flow including strong vertical velocities in spanwise street canyons and flow acceleration in streamwise street canyons. The concentration fields are qualitatively similar albeit the LES shows larger dispersion. A quantitative analysis with performance measures compares the datasets pointwise and demonstrates that the two 3D datasets are similar with respect to many measures including a fractional bias of 0.02 (ideal=0.0), correlation coefficient of 0.87 (ideal = 1.0), and the fraction points within a factor of 2 is 0.98 (ideal = 1.0). Plume analysis compares the arrival and residence time of contaminant and is found to vary significantly with location within the urban environment with arrival times between 0 and 1.25 and differences within the contaminant cloud less than 10% at most locations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Operator learning for energy-efficient building ventilation control with computational fluid dynamics simulation of a real-world classroom

Energy-efficient ventilation control plays an important role in reducing building energy consumption while ensuring occupant health and comfort. While Computational Fluid Dynamics (CFD) simulations provide detailed and physically accurate representations of indoor airflow, their high computational cost limits their use in real-time building control. In this work, we present a neural operator learning framework that combines the physical accuracy of CFD with the computational efficiency of machine learning to enable building ventilation control with the high-fidelity fluid dynamics models. Our method jointly optimizes the airflow supply rates and vent angles to reduce energy use and adhere to air quality constraints. We train an ensemble of neural operator transformer models to learn the mapping from building control actions to airflow fields using high-resolution CFD data. This learned neural operator is then embedded in an optimization-based control framework for building ventilation control. Experimental results show that our approach achieves significant energy savings compared to maximum airflow rate control, rule-based control, as well as data-driven control methods using spatially averaged CO 2 prediction and deep learning–based reduced-order models, while consistently maintaining safe indoor air quality. These results highlight the practicality and scalability of our method in maintaining energy efficiency and indoor air quality in real-world buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2023

Robust scalable initialization for Bayesian variational inference with multi-modal Laplace approximations

Predictive modeling typically relies on Bayesian model calibration to provide uncertainty quantification. Variational inference utilizing fully independent (“mean-field”) Gaussian distributions are often used as approximate probability density functions. This simplification is attractive since the number of variational parameters grows only linearly with the number of unknown model parameters. However, the resulting diagonal covariance structure and unimodal behavior can be too restrictive to provide useful approximations of intractable Bayesian posteriors that exhibit highly non-Gaussian behavior, including multimodality. High-fidelity surrogate posteriors for these problems can be obtained by considering the family of Gaussian mixtures. Gaussian mixtures are capable of capturing multiple modes and approximating any distribution to an arbitrary degree of accuracy, while maintaining some analytical tractability. Unfortunately, variational inference using Gaussian mixtures with full-covariance structures suffers from a quadratic growth in variational parameters with the number of model parameters. The existence of multiple local minima due to strong nonconvex trends in the loss functions often associated with variational inference present additional complications, These challenges motivate the need for robust initialization procedures to improve the performance and computational scalability of variational inference with mixture models. In this work, we propose a method for constructing an initial Gaussian mixture model approximation that can be used to warm-start the iterative solvers for variational inference. The procedure begins with a global optimization stage in model parameter space. In this step, local gradient-based optimization, globalized through multistart, is used to determine a set of local maxima, which we take to approximate the mixture component centers. Around each mode, a local Gaussian approximation is constructed via the Laplace approximation. Finally, the mixture weights are determined through constrained least squares regression. The robustness and scalability of the proposed methodology is demonstrated through application to an ensemble of synthetic tests using high-dimensional, multimodal probability density functions. Here, the practical aspects of the approach are demonstrated with inversion problems in structural dynamics.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Peatland fires in Alaska will double by the end of the century

During recent summers, warm and dry conditions have increased the occurrence of wildfires and potentially peat-fires across Alaska. Limitations in resolving the fine-scale distribution of peatlands and climate observations have constrained our ability to accurately predict peat-fire dynamics. Using a new high-resolution peatland map of Alaska, we evaluated the climate and environmental controls of past and future peat-fire activity. Ensemble machine learning models identified reduced soil moisture, higher temperatures, and evapotranspiration as key predictors of annual total burned peatland area (tenfold CV R 2 = 0.62, RMSE = 221.1 km 2 ). By the end of the twenty-first century, models forced with climate datasets from representative concentration pathways (RCPs) 4.5, 6.0, and 8.5 emission scenarios project a statewide doubling of burned peatlands (increasing 61–121%), with regional increases ranging from 25–165% in polar, 61–95% in boreal, and 102–106% in maritime ecoregions. These projections indicate that wildfires will progressively encroach further into organic-rich moist and wet peaty soils, potentially amplifying soil carbon release across Alaska.

climate-change ecology↗

DOE OSTI2022

Artificial spin ice phase-change memory resistors

Abstract We present a proposal for realization of an electrical memory reminiscent of a memristor in connected Kagome artificial spin ice. We show that current flowing through the system alters the magnetic ensemble, which in turns controls the overall resistance thus leaving memory of current passage in the system. This introduces a current-dependent effect for a dynamic resistive state. We simulate a spin-induced thermal phase-change mechanism, and an athermal domain-wall spin inversion. In both cases we observe electrical memory behavior with an I – V hysteretic pinched loop, typical of memristors. These results can be extended to the more complex geometries in which artificial spin ice can be designed to engineer the hysteresis curve.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Extending the Dynamic Range of Galaxy Outflow Scaling Relations: Massive Compact Galaxies with Extreme Outflows

We investigate galactic winds in the HizEA galaxies, a collection of 46 late-stage galaxy mergers at z = 0.4–0.8, with stellar masses of $\mathrm{log}({M}_{* }/{M}_{\odot })=10.4\mbox{--}11.5$ , star formation rates (SFRs) of 20–500 M ⊙ yr -1 , and ultra-compact (a few 100 pc) central star-forming regions. We measure their gas kinematics using the Mg II λλ 2796,2803 absorption lines in optical spectra from MMT, Magellan, and Keck. We find evidence of outflows in 90% of targets, with maximum outflow velocities of 550–3200 km s -1 . We combine these data with ten samples from the literature to construct scaling relations for outflow velocity versus SFR, star formation surface density (Σ SFR ), M * , and SFR/M * . The HizEA galaxies extend the dynamic range of the scaling relations by a factor of ~2–4 in outflow velocity and an order of magnitude in SFR and Σ SFR . The ensemble scaling relations exhibit strong correlations between outflow velocity, SFR, SFR/R, and Σ SFR , and weaker correlations with M * and SFR/M * . The HizEA galaxies are mild outliers on the SFR and M* scaling relations, but they connect smoothly with more typical star-forming galaxies on plots of outflow velocity versus SFR/R and Σ SFR . These results provide further evidence that the HizEA galaxies' exceptional outflow velocities are a consequence of their extreme star formation conditions rather than hidden black hole activity, and they strengthen previous claims that Σ SFR is one of the most important properties governing the velocities of galactic winds.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1981

A numerical model study of long-term planetary wave predictability

A numerical model was developed to study planetary wave behavior over periods up to 4 months under the influences of stochastic and diffusive forcing by synoptic scale motions and of sea surface temperature anomalies. The model is suitable for Monte Carlo experiments which provide ensemble mean and variance predictions and permit the response to heating anomalies to be partly separated from random variability. The model incorporates the linear-balanced dynamical equations for wind at two tropospheric levels and temperature at one level. The model atmosphere is forced by semiempirical January heating linearly dependent on local air temperature. Topography and the lower boundary layer are represented by vertical motion and drag functions, respectively, of the lower level wind. The linear phase speeds of Rossby waves are modified to reflect the effect of the stratosphere.

Wobus, R. L.↗