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 577 records · Page 32

Relation between pore size and the compressibility of a confined fluid

When a fluid is confined to a nanopore, its thermodynamic properties differ from the properties of a bulk fluid, so measuring such properties of the confined fluid can provide information about the pore sizes. Here, we report a simple relation between the pore size and isothermal compressibility of argon confined in such pores. Compressibility is calculated from the fluctuations of the number of particles in the grand canonical ensemble using two different simulation techniques: conventional grand-canonical Monte Carlo and grand-canonical ensemble transition-matrix Monte Carlo. Our results provide a theoretical framework for extracting the information on the pore sizes of fluid-saturated samples by measuring the compressibility from ultrasonic experiments.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Massively parallel Monte Carlo for many-particle simulations on GPUs

Current trends in parallel processors call for the design of efficient massively parallel algorithms for scientific computing. Parallel algorithms for Monte Carlo simulations of thermodynamic ensembles of particles have received little attention because of the inherent serial nature of the statistical sampling. In this paper, we present a massively parallel method that obeys detailed balance and implement it for a system of hard disks on the GPU. We reproduce results of serial high-precision Monte Carlo runs to verify the method. This is a good test case because the hard disk equation of state over the range where the liquid transforms into the solid is particularly sensitive to small deviations away from the balance conditions. On a Tesla K20, our GPU implementation executes over one billion trial moves per second, which is 148 times faster than on a single Intel Xeon E5540 CPU core, enables 27 times better performance per dollar, and cuts energy usage by a factor of 13. With this improved performance we are able to calculate the equation of state for systems of up to one million hard disks. These large system sizes are required in order to probe the nature of the melting transition, which has been debated for the last forty years. In this paper we present the details of our computational method, and discuss the thermodynamics of hard disks separately in a companion paper.

97 MATHEMATICS AND COMPUTING↗

Improving Simulated Soil Moisture Fields Through Assimilation of AMSR-E Soil Moisture Retrievals with an Ensemble Kalman Filter and a Mass Conservation Constraint

Model simulated soil moisture fields are often biased due to errors in input parameters and deficiencies in model physics. Satellite derived soil moisture estimates, if retrieved appropriately, represent the spatial mean of soil moisture in a footprint area, and can be used to reduce model bias (at locations near the surface) through data assimilation techniques. While assimilating the retrievals can reduce model bias, it can also destroy the mass balance enforced by the model governing equation because water is removed from or added to the soil by the assimilation algorithm. In addition, studies have shown that assimilation of surface observations can adversely impact soil moisture estimates in the lower soil layers due to imperfect model physics, even though the bias near the surface is decreased. In this study, an ensemble Kalman filter (EnKF) with a mass conservation updating scheme was developed to assimilate the actual value of Advanced Microwave Scanning Radiometer (AMSR-E) soil moisture retrievals to improve the mean of simulated soil moisture fields by the Noah land surface model. Assimilation results using the conventional and the mass conservation updating scheme in the Little Washita watershed of Oklahoma showed that, while both updating schemes reduced the bias in the shallow root zone, the mass conservation scheme provided better estimates in the deeper profile. The mass conservation scheme also yielded physically consistent estimates of fluxes and maintained the water budget. Impacts of model physics on the assimilation results are discussed.

Li, Bailing↗

Earth's energy imbalance since 1960 in observations and CMIP5 models: Earth's energy imbalance since 1960

Observational analyses of running 5 year ocean heat content trends (Ht) and net downward top of atmosphere radiation (N) are significantly correlated (r ~ 0.6) from 1960 to 1999, but a spike in Ht in the early 2000s is likely spurious since it is inconsistent with estimates of N from both satellite observations and climate model simulations. Variations in N between 1960 and 2000 were dominated by volcanic eruptions and are well simulated by the ensemble mean of coupled models from the Fifth Coupled Model Intercomparison Project (CMIP5). Here, we find an observation-based reduction in N of -0.31 ± 0.21 W m -2 between 1999 and 2005 that potentially contributed to the recent warming slowdown, but the relative roles of external forcing and internal variability remain unclear. Finally, while present-day anomalies of N in the CMIP5 ensemble mean and observations agree, this may be due to a cancelation of errors in outgoing longwave and absorbed solar radiation.

58 GEOSCIENCES↗

Enabling immersive engagement in energy system models with deep learning

Complex ensembles of energy simulation models have become crucial components of renewable energy research in recent years. Often the increasing computational cost, high-dimensional structure, and other complexities hinder researchers from fully utilizing these data sources for knowledge building. Researchers at National Renewable Energy Laboratory have determined an immersive visualization workflow to dramatically improve user engagement and analysis capability through a combination of low-dimensional structure analysis, deep learning, and custom visualization methods. We present case studies for two energy simulation platforms.

97 MATHEMATICS AND COMPUTING↗

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↗

CoRE MOF DB: A curated experimental metal-organic framework database with machine-learned properties for integrated material-process screening

Here, we present an updated version of the Computation-Ready, Experimental (CoRE) Metal-Organic Framework (MOF) database, which includes a curated set of computation-ready MOF crystal structures designed for high-throughput computational materials discovery. Data collection and curation procedures were improved from the previous version to enable more frequent updates in the future. Machine-learning-predicted properties, such as stability metrics and heat capacities, are included in the dataset to streamline screening activities. An updated version of MOFid was developed to provide detailed information on metal nodes, organic linkers, and topologies of an MOF structure. DDEC6 partial atomic charges of MOFs were assigned based on a machine-learning model. Gibbs ensemble Monte Carlo simulations were used to classify the hydrophobicity of MOFs. The finalized dataset was subsequently used to perform integrated material-process screening for various carbon-capture conditions using high-fidelity temperature-swing adsorption (TSA) simulations. Our workflow identified multiple MOF candidates that are predicted to outperform CALF-20 for these applications.

CoRE MOF database↗

High-Temperature Water Adsorption Isotherms and Ambient Temperature Water Diffusion Rates on Water Harvesting Metal–Organic Frameworks

Water adsorption isotherms from 25 to 125 °C were measured for three metal−organic frameworks (MOFs), MOF-303, MOF-LA2-1, and MIL-100(Fe), which are frequently studied for water harvesting applications. The results show how the step in the water adsorption isotherm varies as a function of temperature and detail the combination of pressure and temperature necessary to remove adsorbed water. Furthermore, isobaric−isothermal Gibbs ensemble Monte Carlo simulations performed for MOF-303 shed light on the change in occupation numbers of the different known water adsorption sites with increasing temperature. Additionally, the diffusion rates of water through these materials were measured using concentration swing frequency response, and micropore diffusion was identified as the controlling mechanism. The Darken relation was used to show the dependence of the diffusion rate on the concentration and the impact of the adsorption isotherm slope. The adsorption of water on MOF-LA2-1 is faster than that on MIL-100(Fe). These data show that MOF-LA2-1 with its high-water adsorption capacity, quick adsorption rate, and favorable desorption energetics is a leading candidate for atmospheric water harvesting.

Adsorption↗

Nanoscale Probing of Dynamics in Local Molecular Environments

Vibrational spectroscopy can provide information about structure, coupling, and dynamics underlying the properties of complex molecular systems. While measurements of spectral line broadening can probe local chemical environments, the spatial averaging in conventional spectroscopies limits insight into underlying heterogeneity, in particular in disordered molecular solids. Here, using femtosecond infrared scattering scanning near-field optical microscopy (IR s-SNOM), we resolve in vibrational free-induction decay (FID) measurements a high degree of spatial heterogeneity in polytetrafluoroethylene (PTFE) as a dense molecular model system. In nanoscopic probe volumes as small as 103 vibrational oscillators, we approach the homogeneous response limit, with extended vibrational dephasing times of several picoseconds, that is, up to 10 times the inhomogeneous lifetime, and spatial average converging to the bulk ensemble response. We simulate the dynamics of relaxation with a finite set of local vibrational transitions subject to random modulations in frequency. The combined results suggest that the observed heterogeneity arises due to static and dynamic variations in the local molecular environment. This approach thus provides real-space and real-time visualization of the subensemble dynamics that define the properties of many functional materials.

Atkin, Joanna M.↗

Enhanced Luminescent Stability through Particle Interactions in Silicon Nanocrystal Aggregates

Close-packed assemblies of ligand-passivated colloidal nanocrystals can exhibit enhanced photoluminescent stability, but the origin of this effect is unclear. Here, we use experiment, simulation, and ab initio computation to examine the influence of interparticle interactions on the photoluminescent stability of silicon nanocrystal aggregates. The time-dependent photoluminescence emitted by structures ranging in size from a single quantum dot to agglomerates of more than a thousand is compared with Monte Carlo simulations of noninteracting ensembles using measured single-particle blinking data as input. In contrast to the behavior typically exhibited by the metal chalcogenides, the measured photoluminescent stability shows an enhancement with respect to the noninteracting scenario with increasing aggregate size. We model this behavior using time-dependent density functional theory calculations of energy transfer between neighboring nanocrystals as a function of nanocrystal size, separation, and the presence of charge and/or surface-passivation defects. In conclusion, our results suggest that rapid exciton transfer from “bright” nanocrystals to surface trap states in nearest-neighbors can efficiently fill such traps and enhance the stability of emission by promoting the radiative recombination of slowly diffusing excited electrons.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Role of Arctic Sea Ice in the 2014–2015 Eurasian Warm Winter

Abstract Because of large internal variabilities, the role of Arctic sea ice concentration (SIC) on regional climate cannot be easily separated. This study uses two groups of large ensemble atmospheric model simulations under climatological and observed SIC boundary conditions to investigate the role of SIC in the 2014–2015 December to February Eurasian warm winter (DJF15). It is shown that the SIC has large impact on the probability distribution function of the DJF15 temperature and pressure fields. The anomalous high Barents Sea ice during the 2014–2015 autumn and winter leads to significant shift in the probability distribution function skewness of the DJF15 surface temperature (from −0.13 to −0.48) and the related sea level pressure (from −0.18 to 0.32) that favor more occurrences of warm temperature anomaly and positive North Atlantic Oscillation‐like pattern. This asymmetry is consistent with anomalous forcing in phase with the anomalies of the sea level pressure field.

Xie, Jinbo↗

The Influence of Volcanic Aerosol Meridional Structure on Monsoon Responses over the Last Millennium

Monsoon responses to eruptions over the last millennium (LM) are examined in an ensemble of climate simulations as a function of eruption hemisphere. In this work, a composite analysis reveals a particularly strong sensitivity of monsoon rainfall in the year following Northern Hemisphere (NH) extratropical eruptions. Additional analysis focusing on the 18th century eruption of Mt. Laki and idealized simulations representing an analogue Southern Hemisphere eruption (SH-Laki) reveal monsoon responses that are approximately symmetric across hemispheres, despite exhibiting asymmetries in other aspects of the climate response. We conclude that 1) latitudinally mirrored eruptions result in approximately symmetric monsoon responses, 2) disproportionate weakening (strengthening) of NH (SH) monsoons by NH eruptions over the LM resulted in part from their relatively high latitudes, and 3) uncertainty in eruption latitude fundamentally limits our ability to accurately simulate associated monsoon and tropical precipitation responses in nature.

54 ENVIRONMENTAL SCIENCES↗

Uncertainty in Reconstructing Paleo‐Elevation of the Antarctic Ice Sheet From Temperature‐Sensitive Ice Core Records

Paleotemperature reconstructions from ice cores are mixed signals of changes in climate and ice‐surface elevation. A common, temperature‐based paleoaltimetry method suggests these signals can be disentangled by comparing two proxy locations with similar climates. The difference between the records is assumed to be due to elevation, which is estimated by scaling the temperature difference by a lapse rate. We investigate the uncertainty associated with this approach using a case study of the Antarctic Ice Sheet during the Last Glacial Maximum. From an ensemble of climate simulations, we extract modeled temperatures at locations of real ice cores. We find uncertainty on the order of hundreds of meters that results from spatial heterogeneity in non‐adiabatic temperature change, which itself stems in part from elevation‐induced atmospheric circulation change. Our findings suggest that caution is needed when interpreting temperature‐based paleoaltimetry results for ice sheets.

58 GEOSCIENCES↗

Convective Updraft and Downdraft Characteristics of Continental Mesoscale Convective Systems in the Model Gray Zone

Here, the "gray zone" of convective modeling is defined as the range of horizontal grid spacings (Δx) at which turbulent transport processes are only partially resolved by the dynamics of the numerical model. This zone typically covers Δx from a few kilometers to several hundred meters, wherein the realistic representation of convective cloud processes can be challenging. This study characterizes the convective draft behaviors at multiple Δx across the gray zone and determines the appropriate Δx that can reliably capture these salient convective properties. We perform an ensemble of idealized simulations of mesoscale convective systems (MCS) using the Weather Research and Forecasting model at various Δx from 4 km to 250 m over the central U.S. An evaluation of key MCS kinematic properties is constrained using unique, long-term vertical velocity estimates obtained by radar wind profilers deployed by the Department of Energy Atmospheric Radiation Measurement user facility. MCS simulations for all Δx tested overestimate (underestimate) the probabilities of convective updrafts (downdrafts) compared to the observations. In terms of the convective draft intensity, finer-Δx models overestimate the updraft intensity, while the opposite is found for downdrafts. Moving from Δx = 4 km to 250 m, downdrafts become stronger and more frequent especially at middle and upper levels, attributed to additional drag from an increasing graupel frequency and compensation of enhanced updraft velocity. Simulated draft characteristics, including core size, intensity, and probability of occurrence, exhibit pronounced changes at Δx greater than 500 m, but suggest less sensitivity when Δx is reduced below 500 m.

54 ENVIRONMENTAL SCIENCES↗

Future Projections of EURO‐CORDEX Raw and Bias‐Corrected Daily Maximum Wind Speeds Over Scandinavia

Twenty historical and future Coordinated Regional Climate Downscaling Experiment ensemble for Europe simulations are bias corrected to investigate the future changes in the daily maximum wind speed over Scandinavia. We use quantile mapping to adjust the wind for the historical period (1985–2014) and quantile delta mapping for two future periods (2041–2070, 2071–2100, RCP8.5) with the 3-km spatial resolution NORA3 hindcast as the reference data set. Decomposing the variance, we find that most of the inter-model spread in the bias and the response to climate change is due to the regional climate model over land and mainly to the general climate model over sea. On average, the mean daily maximum wind speed is projected to increase everywhere except along the western coast of Norway and Denmark. In summer, we see an opposite sign over Sweden and Finland. The Norwegian Sea experiences weaker mean and high wind whereas the Baltic Sea experiences a strengthening. Bias correction influences the amplitude of the response, not the response pattern. Wind speed distributions can have different shapes in the future, with, for example, a flattening of the distribution off the coast of Norway with more frequent weak and very strong winds. Apart from a few locations in Norway, we find an increase in the number of days in the local highest historical wind category. Overall, summer exhibits opposite signals compared to the three other seasons. At country scale, the sign of the change in the mean and 98th percentile varies greatly depending on the region and among the simulations, especially for Norway.

54 ENVIRONMENTAL SCIENCES↗

Scenario Storyline Discovery for Planning in Multi‐Actor Human‐Natural Systems Confronting Change

Scenarios have emerged as valuable tools in managing complex human-natural systems, but the traditional approach of limiting focus on a small number of predetermined scenarios can inadvertently miss consequential dynamics, extremes, and diverse stakeholder impacts. Exploratory modeling approaches have been developed to address these issues by exploring a wide range of possible futures and identifying those that yield consequential vulnerabilities. However, vulnerabilities are typically identified based on aggregate robustness measures that do not take full advantage of the richness of the underlying dynamics in the large ensembles of model simulations and can make it hard to identify key dynamics and/or storylines that can guide planning or further analyses. This study introduces the FRamework for Narrative Storylines and Impact Classification (FRNSIC; pronounced “forensic”): a scenario discovery framework that addresses these challenges by organizing and investigating consequential scenarios using hierarchical classification of diverse outcomes across actors, sectors, and scales, while also aiding in the selection of scenario storylines, based on system dynamics that drive consequential outcomes. We present an application of this framework to the Upper Colorado River Basin, focusing on decadal droughts and their water scarcity implications for the basin's diverse users and its obligations to downstream states through Lake Powell. We show how FRNSIC can explore alternative sets of impact metrics and drought dynamics and use them to identify drought scenario storylines, that can be used to inform future adaptation planning.

54 ENVIRONMENTAL SCIENCES↗

Urban Versus Lake Impacts on Heat Stress and Its Disparities in a Shoreline City

Abstract Shoreline cities are influenced by both urban‐scale processes and land‐water interactions, with consequences on heat exposure and its disparities. Heat exposure studies over these cities have focused on air and skin temperature, even though moisture advection from water bodies can also modulate heat stress. Here, using an ensemble of model simulations covering Chicago, we find that Lake Michigan strongly reduces heat exposure (2.75°C reduction in maximum average air temperature in Chicago) and heat stress (maximum average wet bulb globe temperature reduced by 0.86°C) during the day, while urbanization enhances them at night (2.75 and 1.57°C increases in minimum average air and wet bulb globe temperature, respectively). We also demonstrate that urban and lake impacts on temperature (particularly skin temperature), including their extremes, and lake‐to‐land gradients, are stronger than the corresponding impacts on heat stress, partly due to humidity‐related feedback. Likewise, environmental disparities across community areas in Chicago seen for skin temperature are much higher (1.29°C increase for maximum average values per $10,000 higher median income per capita) than disparities in air temperature (0.50°C increase) and wet bulb globe temperature (0.23°C increase). The results call for consistent use of physiologically relevant heat exposure metrics to accurately capture the public health implications of urbanization.

54 ENVIRONMENTAL SCIENCES↗

Estimating the CO 2 Fertilization Effect on Extratropical Forest Productivity From Flux‐Tower Observations

Abstract The land sink of anthropogenic carbon emissions, a crucial component of mitigating climate change, is primarily attributed to the CO 2 fertilization effect on global gross primary productivity (GPP). However, direct observational evidence of this effect remains scarce, hampered by challenges in disentangling the CO 2 fertilization effect from other long‐term confounding drivers, particularly climatic changes. Here, we introduce a novel statistical approach to separate the CO 2 fertilization effect on photosynthetic carbon uptake using eddy covariance (EC) records across 38 extratropical forest sites. We find the median stimulation rate of GPP to be 3.2 ± 0.9 gC m −2 yr −1 ppm −1 (or 16.4 ± 4.2% per 100 ppm) under increasing atmospheric CO 2 across these sites, respectively. To validate the robustness of our findings, we test our statistical method using factorial simulations of an ensemble of process‐based land surface models. We address additional factors, including nitrogen deposition and land management, that may impact plant productivity, potentially confounding the attribution to the CO 2 fertilization effect. Assuming these site‐specific effects offset to some extent across sites as random factors, the estimated median value still reflects the strength of the CO 2 fertilization effect. However, disentanglement of these long‐term effects, often inseparable by timescale, requires further causal research. Our study provides direct evidence that the photosynthetic stimulation is maintained under long‐term CO 2 fertilization across multiple EC sites. Such observation‐based quantification is key to constraining the long‐standing uncertainties in the land carbon cycle under rising CO 2 concentrations.

Environmental Sciences & Ecology↗