SEARCH · Engineering Papers

Results for “ensemble simulations”

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 343 records · Page 19

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Molecular Simulations Probing the Adsorption and Diffusion of Ammonia, Nitrogen, Hydrogen, and Their Mixtures in Bulk MFI Zeolite and MFI Nanosheets at High Temperature and Pressure

Recent advances in the synthesis of MFI zeolite nanosheets have led to highly selective membranes that are promising candidates for small-scale ammonia separation from ammonia/nitrogen/hydrogen mixtures in distributed green ammonia production plants. Using force-field-based molecular simulations, we evaluate the performance of bulk all-silica MFI zeolite and a 3 nm thick all-silica MFI nanosheet for ammonia/nitrogen/hydrogen separations at moderate temperature and pressure conditions (T = 373 K, p = 5 bar) as well as conditions relevant for a membrane-based reactor–separator process (T = 523 or 623 K, p = 80 bar). Isobaric–isothermal Gibbs ensemble Monte Carlo (GEMC) simulations were carried out to understand the selective adsorption behavior for these mixtures. Molecular dynamics simulations in the canonical ensemble were performed to examine the transport behavior of ammonia in a mixture with nitrogen or hydrogen at loadings obtained from the GEMC simulations. Our results show that both bulk MFI and the nanosheet with explicit surface silanols are highly selective toward ammonia adsorption, but the adsorption selectivity decreases by factors of 4 and 10 for ammonia/nitrogen and ammonia/hydrogen mixtures as the temperature is increased from 373 to 623 K. Conversely, the diffusion selectivities toward ammonia are more favorable at process-relevant temperatures and for the MFI nanosheet. Here, at T = 523 K and p = 80 bar, our simulations predict overall separation factors of 3.8 ± 0.2 and 4.5 ± 0.3 for ammonia/nitrogen and ammonia/hydrogen mixtures, respectively, in the bulk MFI zeolite, and separation factors of 8.6 ± 0.4 and 12.5 ± 0.8, respectively, for the MFI nanosheet.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2003

Spectral Retrieval of Latent Heating Profiles from TRMM PR Data: Comparison of Look-Up Tables

The primary goal of the Tropical Rainfall Measuring Mission (TRMM) is to use the information about distributions of precipitation to determine the four dimensional (i.e., temporal and spatial) patterns of latent heating over the whole tropical region. The Spectral Latent Heating (SLH) algorithm has been developed to estimate latent heating profiles for the TRMM Precipitation Radar (PR) with a cloud- resolving model (CRM). The method uses CRM- generated heating profile look-up tables for the three rain types; convective, shallow stratiform, and anvil rain (deep stratiform with a melting level). For convective and shallow stratiform regions, the look-up table refers to the precipitation top height (PTH). For anvil region, on the other hand, the look- up table refers to the precipitation rate at the melting level instead of PTH. For global applications, it is necessary to examine the universality of the look-up table. In this paper, we compare the look-up tables produced from the numerical simulations of cloud ensembles forced with the Tropical Ocean Global Atmosphere (TOGA) Coupled Atmosphere-Ocean Response Experiment (COARE) data and the GARP Atlantic Tropical Experiment (GATE) data. There are some notable differences between the TOGA-COARE table and the GATE table, especially for the convective heating. First, there is larger number of deepest convective profiles in the TOGA-COARE table than in the GATE table, mainly due to the differences in SST. Second, shallow convective heating is stronger in the TOGA COARE table than in the GATE table. This might be attributable to the difference in the strength of the low-level inversions. Third, altitudes of convective heating maxima are larger in the TOGA COARE table than in the GATE table. Levels of convective heating maxima are located just below the melting level, because warm-rain processes are prevalent in tropical oceanic convective systems. Differences in levels of convective heating maxima probably reflect differences in melting layer heights. We are now extending our study to simulations of other field experiments (e.g. SCSMEX and ARM) in order to examine the universality of the look-up table. The impact of look-up tables on the retrieved latent heating profiles will also be assessed.

Shige, Shoichi↗

DOE OSTI2020

An analytic asymmetric-piston model for the impact of mode-1 shell asymmetry on ICF implosions

For many years, low mode asymmetry in inertially confined fusion (ICF) implosions has been recognized as a potential performance limiting factor, but analysis has been limited to using simulations and searching for data correlations. In this paper, an analytically solvable model based upon the simple picture of an asymmetric piston is presented. Asymmetry of the shell driving the implosion, as opposed to asymmetry in the hot-spot, is key to the model. The model provides a unifying framework for the action of mode-1 shell asymmetry and the resulting connections between various diagnostic signatures. A key variable in the model is the shell asymmetry fraction, f, which is related to the areal density variation of the shell surrounding the hot-spot. It is shown that f is simply related to the observed hot-spot mode-1 velocity and to the concept of residual energy in an implosion. The model presented in this paper yields explicit expressions for the hot-spot diameter, stagnation pressure, hot-spot energy, inertial confinement-time, Lawson parameter, hot-spot temperature, and fusion yield under the action of mode-1 asymmetry. Agreement is found between the theory scalings when compared to ICF implosion data from the National Ignition Facility and to large ensembles of detailed simulations, making the theory a useful tool for interpreting data. The theory provides a basis for setting tolerable limits on asymmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

The role of SAXS and molecular simulations in 3D structure elucidation of a DNA aptamer against lung cancer

Aptamers are short, single-stranded DNA or RNA oligonucleotide molecules that function as synthetic analogs of antibodies and bind to a target molecule with high specificity. Aptamer affinity entirely depends on its tertiary structure and charge distribution. Therefore, length and structure optimization are essential for increasing aptamer specificity and affinity. Here, we present a general optimization procedure for finding the most populated atomistic structures of DNA aptamers. Based on the existed aptamer LC-18 for lung adenocarcinoma, a new truncated LC-18 (LC-18t) aptamer LC-18t was developed. A three-dimensional (3D) shape of LC-18t was reported based on small-angle X-ray scattering (SAXS) experiments and molecular modeling by fragment molecular orbital or molecular dynamic methods. Molecular simulations revealed an ensemble of possible aptamer conformations in solution that were in close agreement with measured SAXS data. The aptamer LC-18t had stronger binding to cancerous cells in lung tumor tissues and shared the binding site with the original larger aptamer. The suggested approach reveals 3D shapes of aptamers and helps in designing better affinity probes.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Hydrodynamic computations of high-power laser drives generating metal ejecta jets from surface grooves

Understanding dynamic fragmentation in shock-loaded metals and predicting properties of the resulting ejecta are of considerable importance for both basic and applied science. The nature of material ejection has been shown to change drastically when the free surface melts on compression or release. In this work, we present hydrodynamic simulations of laser-driven microjetting from micron-scale grooves on a tin surface. We study microjet formation across a range of shock strengths from drives that leave the target solid after release to drives that induce shock melting in the target. The shock-state particle velocity (Up) varies from 0.3 to 3 km/s and the shock breakout pressure is 3–120 GPa. The microjet tip velocity is 1–8 km/s and the free-surface velocity varies from 0.1 to 5 km/s. Two tin equations of state are examined: a “soft” model (LEOS 501) where the target melts for U p > 1 km/s and a more detailed multiphase model (SESAME 2161) that melts for U p > 1.4 km/s. We use these two models to examine the influence of phase change and the choice of the material model on microjet formation and evolution. We observe in our computational results that jet formation can be classified into three regimes: a low-energy regime where material strength affects jet formation, a moderate-energy regime dominated by the changing phase of tin material, and a high-energy regime where results are insensitive to the material model and jet formation is described by an idealized steady-jet theory. Using an ensemble of 2D simulations, we show that these trends hold across a wide range of drive energies and groove angles.

36 MATERIALS SCIENCE↗

NASA NTRS2006

Radiative Forcing by Well-Mixed Greenhouse Gases: Estimates from Climate Models in the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report (AR4)

The radiative effects from increased concentrations of well-mixed greenhouse gases (WMGHGs) represent the most significant and best understood anthropogenic forcing of the climate system. The most comprehensive tools for simulating past and future climates influenced by WMGHGs are fully coupled atmosphere-ocean general circulation models (AOGCMs). Because of the importance of WMGHGs as forcing agents it is essential that AOGCMs compute the radiative forcing by these gases as accurately as possible. We present the results of a radiative transfer model intercomparison between the forcings computed by the radiative parameterizations of AOGCMs and by benchmark line-by-line (LBL) codes. The comparison is focused on forcing by CO2, CH4, N2O, CFC-11, CFC-12, and the increased H2O expected in warmer climates. The models included in the intercomparison include several LBL codes and most of the global models submitted to the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report (AR4). In general, the LBL models are in excellent agreement with each other. However, in many cases, there are substantial discrepancies among the AOGCMs and between the AOGCMs and LBL codes. In some cases this is because the AOGCMs neglect particular absorbers, in particular the near-infrared effects of CH4 and N2O, while in others it is due to the methods for modeling the radiative processes. The biases in the AOGCM forcings are generally largest at the surface level. We quantify these differences and discuss the implications for interpreting variations in forcing and response across the multimodel ensemble of AOGCM simulations assembled for the IPCC AR4.

Collins, W. D.↗

NASA NTRS2011

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↗

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 OSTI2025

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↗

DOE OSTI2024

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↗

DOE OSTI2022

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↗

DOE OSTI2022

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↗

DOE OSTI2023

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↗

DOE OSTI2024

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↗

DOE OSTI2023

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↗

DOE OSTI2024

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↗

DOE OSTI2021

Increasing large wildfires over the western United States linked to diminishing sea ice in the Arctic

Abstract The compound nature of large wildfires in combination with complex physical and biophysical processes affecting variations in hydroclimate and fuel conditions makes it difficult to directly connect wildfire changes over fire-prone regions like the western United States (U.S.) with anthropogenic climate change. Here we show that increasing large wildfires during autumn over the western U.S. are fueled by more fire-favorable weather associated with declines in Arctic sea ice during preceding months on both interannual and interdecadal time scales. Our analysis (based on observations, climate model sensitivity experiments, and a multi-model ensemble of climate simulations) demonstrates and explains the Arctic-driven teleconnection through regional circulation changes with the poleward-shifted polar jet stream and enhanced fire-favorable surface weather conditions. The fire weather changes driven by declining Arctic sea ice during the past four decades are of similar magnitude to other leading modes of climate variability such as the El Niño-Southern Oscillation that also influence fire weather in the western U.S.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Increased Asian aerosols drive a slowdown of Atlantic Meridional Overturning Circulation

Observational evidence and climate model experiments suggest a slowdown of the Atlantic Meridional Overturning Circulation (AMOC) since the mid-1990s. Increased greenhouse gases and the declined anthropogenic aerosols (AAs) over North America and Europe are believed to contribute to the AMOC slowdown. Asian AAs continue to increase but the associated impact has been unclear. Using ensembles of climate simulations, here we show that the radiative cooling resulting from increased Asian AAs drives an AMOC reduction. The increased AAs over Asia generate circumglobal stationary Rossby waves in the northern midlatitudes, which shift the westerly jet stream southward and weaken the subpolar North Atlantic westerlies. Consequently, reduced transport of cold air from North America hinders water mass transformation in the Labrador Sea and thus contributes to the AMOC slowdown. The link between increased Asian AAs and an AMOC slowdown is supported by different models with different configurations. Thus, reducing emissions of Asian AAs will not only lower local air pollution, but also help stabilize the AMOC.

54 ENVIRONMENTAL SCIENCES↗