Engineering Papers⌕ Search

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.

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

At least 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Grain boundary structure search by using an evolutionary algorithm with effective mutation methods

Grain boundaries (GBs) accommodate the misorientation between adjacent grains in a polycrystalline material. GBs are geometrically described by the macroscopic and microscopic degrees of freedom. Besides, at the atomistic level, GBs exhibit complicated behaviors under varying thermodynamic conditions. The complexity of atomistic GB structures demands stochastic searching for possible states. The effectiveness of stochastic search methods relies on techniques to recreate and select atomistic structures. In this work, we developed a new mutation operator that can induce direct and collective atomistic structure changes to boost the search efficiency of exploring GB structures with evolutionary algorithms (EA). We implemented the mutation methods along with innovative selection, crossover, boundary condition preprocessing methods to form an EA-based package to explore GB structures in grand canonical ensembles with atomistic simulations. We used this package to study the [001] symmetric tilt grain boundaries (STGBs) in FCC copper (Cu), the [110] STGBs in BCC tungsten (W), and the $[1\bar{2}10]$ STGBs in HCP magnesium (Mg). The results show that our design and implementation based on new mutation procedures, selection, and boundary conditions provide a high-quality search of atomistic GB structures in the grand canonical ensemble for different crystal lattices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Iterative sampling of expensive simulations for faster deep surrogate training

Deep neural network (DNN) surrogates of expensive physics simulations are enabling a rapid change in the way that common experimental design and analysis tasks are approached. Surrogate models allow simulations to be performed in parallel and separately from downstream tasks, thereby enabling analyses that would be impossible with the simulation in-the-loop; surrogates based on DNNs can effectively emulate diverse non-scalar data of the types collected in fusion and laboratory-astrophysics experiments. The challenge is in training the surrogate model, for which large ensembles of physics simulations must be run, preferably without wasting computational effort on uninteresting simulations. Here, in this paper, we present an iterative sampling scheme that can preferentially propose simulations in interesting regions of parameter space without neglecting unexplored regions, allowing high-quality and wide-ranging surrogate models to be trained using 2–3 times fewer simulations compare to space-filling designs. Our approach uses an explicit importance function defined on the simulation output space, balanced against a measure of simulation density which serves as a proxy for surrogate accuracy. It is easy to implement and can be tuned to find interesting simulations early in the study, allowing surrogates to be trained quickly and refined as new simulations become available; this represents an important step towards the routine generation of deep surrogate models quickly enough to be truly relevant to experimental work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Simulating the Transport and Rupture of Pollen in the Atmosphere

Pollen, one type of primary biological aerosol particle (PBAP), is emitted from the terrestrial biosphere and can undergo physical changes in the atmosphere via particle rupture. To examine the fate of pollen and its atmospheric processing, a pollen emission and transport scheme is coupled to the Weather Research and Forecasting Model with Chemistry (WRF-Chem). We simulate the emission of pollen and its impacts on the cloud properties and precipitation in the Southern Great Plains from 12 to 19 April 2013, a period with both high pollen emissions and convective activity. We conduct a suite of ensemble runs that simulate primary pollen and three different pollen rupture mechanisms that generate subpollen particles, including (a) high humidity-induced surface rupture, (b) high humidity-induced in-atmosphere plus surface rupture, and (c) lightning-induced rupture, where in-cloud and cloud-to-ground lightning strikes trigger pollen rupture events. When relative humidity is high (>80%), coarse primary pollen (~1 μg m –3 ) is converted into fine subpollen particles (~1.2e –4 μg m –3 ), which produces 80% more subpollen particles than lightning-induced rupture. The in-atmosphere humidity-driven rupture predominantly produces subpollen particles, which is further enhanced during a frontal thunderstorm. During strong convection, vertical updrafts lift primary pollen and subpollen particles (~0.5e –4 μg m –3 ) to the upper troposphere (~12 km) and laterally transports the ruptured pollen in the anvil top outflow. In regions of high pollen and strong convection, ruptured pollen can influence warm cloud formation by decreasing low cloud (<4 km) cloud water mixing ratios and increasing ice phase hydrometeors aloft (>10 km)

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Vegetation Representation Influences Projected Streamflow Changes in the Colorado River Basin

Vegetation parameters for the Variable Infiltration Capacity (VIC) hydrologic model were recently updated using observations from the Moderate Resolution Imaging Spectroradiometer (MODIS). Previous work showed that these MODIS-based parameters improved VIC evapotranspiration simulations when compared to eddy covariance observations. Due to the importance of evapotranspiration within the Colorado River basin, this study provided a basin-by-basin calibration of VIC soil parameters with updated MODIS-based vegetation parameters to improve streamflow simulations. Interestingly, while both configurations had similar historic streamflow performance, end-of-century hydrologic projections, driven by 29 downscaled global climate models under the RCP8.5 emissions scenario, differed between the two configurations. The calibrated MODIS-based configuration had an ensemble mean that simulated little change in end-of-century annual streamflow volume (+0.4%) at Lees Ferry, Arizona, relative to the historical period (1960–2005). In contrast, the previous VIC configuration, which is used to inform decisions about future water resources in the Colorado River basin, projected an 11.7% decrease in annual streamflow. Both VIC configurations simulated similar amounts of evapotranspiration in the historical period. However, the MODIS-based VIC configuration did not show as much of an increase in evapotranspiration by the end of the century, primarily within the upper basin’s forested areas. Differences in evapotranspiration projections were the result of the MODIS-based vegetation parameters having lower leaf area index values and less forested area compared to previous vegetation estimates used in recent Colorado River basin hydrologic projections. These results highlight the need to accurately characterize vegetation and better constrain climate sensitivities in hydrologic models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Changes in Four Decades of Near‐CONUS Tropical Cyclones in an Ensemble of 12 km Thermodynamic Global Warming Simulations

We evaluate tropical cyclones (TCs) in a set of thermodynamic global warming (TGW) simulations over the continental United States (CONUS). A 12 km simulation forced by ERA5 provides a 40‐year historical (1980–2019) control. Four complimentary future scenarios are generated using thermodynamic deltas applied to lateral boundary, interior, and surface forcing. We curate a data set of 4,498 6‐hourly TC snapshots in the control and find a corresponding “twin” in each counterfactual, permitting a paired comparison. Warming results in an increase in mean dynamical TC intensity and moisture‐related quantities, with the latter being more pronounced. TC inner cores contract slightly but outer storm size remains unchanged. The frequency with which TCs become more intense is only moderately consistent, with snapshots having increased hazards ranging from 50% to 80% depending on warming level. The fractions of TCs undergoing rapid intensification and weakening both increase across all warming simulations, suggesting elevated short‐term intensity variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

High Pressure Infiltration/Expulsion of Aqueous NaCl in Planar Hydrophobic Nanopores

Pressure-driven permeation of water in a poorly wettable material results in a conversion of mechanical work into surface free energy representing a new form of energy storage, or absorption. When water is replaced by a concentrated electrolyte solution, the storage capacity of a nanoporous medium becomes comparable to high-end supercapacitors. The addition of salt can also reduce the hysteresis of the infiltration/expulsion cycle. Our molecular simulations provide a theoretical perspective into the mechanisms involved in the process, and underlying structures and interactions in compressed nanoconfined solutions. Specifically, we consider aqueous NaCl in planar confinements of widths of 1.0 nm and 1.64 nm and pressures of up to 3 kbar. Open ensemble Monte Carlo simulations utilizing fractional exchanges of molecules for efficient additions/removal of ions have been utilized in conjunction with pressure-dependent chemical potentials to model bulk phases under pressure. Confinements open to these pressurized bulk, aqueous electrolyte phases show reversibility at narrow pore sizes and irreversibility in wider ones, consistent with experiment, as well as strong hysteresis at both pore size. The addition of salt results in significant increases in the solid/liquid interfacial tension in narrower pores and associated infiltration and expulsion pressures. Furthermore, these changes are consistent with strong desalination effects at the lower pore size observed irrespective of external pressure and initial concentration.

25 ENERGY STORAGE↗
NASA NTRS2023

The Impact of ENSO and NAO Initial Conditions and Anomalies on the Modeled Response to Pinatubo-Sized Volcanic Forcing

Strong, strato-volcanic eruptions are a substantial, intermittent source of natural climate variability. Initial atmospheric and oceanic conditions, such as El Niño Southern Oscillation (ENSO) and the North Atlantic Oscillation (NAO), also naturally impact climate on interannual timescales. We examine how initial conditions of ENSO and NAO contribute to the evolution of climate in the period following a Pinatubo-type eruption using a large (81-member) ensemble of model simulations in GISS model E2.1-G. Simulations are initialized from sampled conditions of ENSO and NAO using the protocol of the coordinated CMIP6 Volcanic Model Intercomparison Project (VolMIP) – where aerosols are forced with respect to time, latitude, and height. We analyze paired anomalous variations (perturbed – control) to understand changes in global and regional climate responses under positive, negative, and neutral ENSO and NAO conditions. In particular, we find that for paired anomalies there is a high probability of strong (∼1.5 °C) warming of northern Eurasia surface air temperature in the first winter after the volcanic eruption for negative NAO ensembles coincident with decreased lower stratospheric temperature at the poles, decreased geopotential height, and strengthening of the stratospheric polar vortex. Climate anomalies (relative to average conditions across the control period), however, show no mean warming and suggest that the strength of this response is impacted by conditions present in the selected period of the control run. Again using paired anomalies, we also observe that under both +ENSO and −ENSO ensembles sea surface temperature decreases in the first post-eruptive boreal winter coinciding with surface cooling from volcanic aerosols. Neutral ENSO ensembles, on the other hand, show variability in their response with no clear trend in post-eruptive warming or cooling. In general, paired anomalies from unperturbed simulations give insight into the evolution of the climate response to volcanic forcing; however, when compared with anomalies from climatological conditions, it is clear that paired anomalies are significantly affected by sampled initial conditions occurring at the time of the volcanic eruption.

volcanic eruptions↗
NASA NTRS2019

Evaluating the Uncertainty of Terrestrial Water Budget Components over High Mountain Asia

This study explores the uncertainties in terrestrial water budget estimation over High Mountain Asia (HMA) using a suite of uncoupled land surface model (LSM) simulations. The uncertainty in the water balance components of precipitation (P), evapotranspiration (ET), runoff (R), and terrestrial water storage (TWS) is significantly impacted by the uncertainty in the driving meteorology, with precipitation being the most important boundary condition. Ten gridded precipitation datasets along with a mix of model-, satellite-, and gauge-based products, are evaluated first to assess their suitability for LSM simulations over HMA. The datasets are evaluated by quantifying the systematic and random errors of these products as well as the temporal consistency of their trends. Though the broader spatial patterns of precipitation are generally well captured by the datasets, they differ significantly in their means and trends. In general, precipitation datasets that incorporate information from gauges are found to have higher accuracy with low Root Mean Square Errors and high correlation coefficient values. An ensemble of LSM simulations with selected subset of precipitation products is then used to produce the mean annual fluxes and their uncertainty over HMA in P, ET, and R to be 2.11 ± 0.45, 1.26 ± 0.11, and 0.85 ± 0.36 mm per day, respectively. The mean annual estimates of the surface mass (water) balance components from this model ensemble are comparable to global estimates from prior studies. However, the uncertainty/spread of P, ET, and R is significantly larger than the corresponding estimates from global studies. A comparison of ET, snow cover fraction, and changes in TWS estimates against remote sensing-based references confirms the significant role of the input meteorology in influencing the water budget characterization over HMA and points to the need for improving meteorological inputs.

Terrestrial water budget↗
DOE OSTI2020

Data-guided Multi-Map variables for ensemble refinement of molecular movies

Driving molecular dynamics simulations with data-guided collective variables offer a promising strategy to recover thermodynamic information from structure-centric experiments. In this study, the three-dimensional electron density of a protein, as it would be determined by cryo-EM or x-ray crystallography, is used to achieve simultaneously free-energy costs of conformational transitions and refined atomic structures. Unlike previous density-driven molecular dynamics methodologies that determine only the best map-model fits, our work employs the recently developed Multi-Map methodology to monitor concerted movements within equilibrium, non-equilibrium, and enhanced sampling simulations. Construction of all-atom ensembles along the chosen values of the Multi-Map variable enables simultaneous estimation of average properties, as well as real-space refinement of the structures contributing to such averages. Using three proteins of increasing size, we demonstrate that biased simulation along the reaction coordinates derived from electron densities can capture conformational transitions between known intermediates. The simulated pathways appear reversible with minimal hysteresis and require only low-resolution density information to guide the transition. The induced transitions also produce estimates for free energy differences that can be directly compared to experimental observables and population distributions. The refined model quality is superior compared to those found in the Protein Data Bank. We find that the best quantitative agreement with experimental free-energy differences is obtained using medium resolution density information coupled to comparatively large structural transitions. Practical considerations for probing the transitions between multiple intermediate density states are also discussed.

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

Equipartition and the Temperature of Maximum Density of TIP4P/2005 Water

Here, we simulate TIP4P/2005 water in the temperature range of 257 to 318 K with time-steps δ = 0.25, 0.50, 1.00, 2.00, and 4.00 fs. The density–temperature behavior obtained using 0.25 or 0.50 fs is in excellent agreement with each other but differs from those obtained using time steps that have been shown earlier to lead to a breakdown of equipartition. For δt = 0.25 or 0.50 fs, the temperature of maximum density (TMD) is 277.15 K and the density value is in close agreement with experiments. For δt = 1.00 fs, the TMD is 277.15 K, but the density value is shifted higher. For the other time steps considered here, the TMD is shifted to progressively lower values for longer time steps, a trend that holds for different thermostat/barostat combinations. Enhancing the water–water dispersion interaction, as has been recommended for simulating disordered proteins in TIP4P/2005, degrades the description of the liquid–vapor phase envelope. We present a simple physically transparent explanation that highlights the separation of the time scales between translational and rotational motion. We also develop a metric, χ, that we term the equipartition anomaly, to detect equipartition violations in simulations that include molecules that are treated as rigid objects. Calculating χ is shown to be straightforward and sensitive to equipartition violations. A key takeaway from this study is that using sufficiently short time steps (≤0.5 fs) to preserve equipartition is essential for obtaining meaningful liquid water properties and for producing reliable simulation data, as correct ensemble sampling is fundamental to ensure reproducibility across codes and simulation algorithms.

Asthagiri, Dilipkumar N. [Oak Ridge National Labor↗
DOE OSTI2026

High-Performance Semiempirical Excited-State Molecular Dynamics Powered by Graphics Processing Units

Here, this Letter introduces excited-state molecular dynamics in PYSEQM, a GPU-accelerated semiempirical quantum chemistry engine implemented in PyTorch. The new module enables Born–Oppenheimer molecular dynamics (BOMD) using configuration-interaction singles and random phase approximation for excited states, allowing long trajectories and large statistical ensembles to be simulated efficiently on a single GPU. We also implement an extended Lagrangian excited-state BOMD (XL-ESMD) scheme that propagates auxiliary electronic variables, enabling relaxed ground and excited-state convergence thresholds without compromising energy conservation. The excited-state BOMD implementation scales smoothly from small chromophores to a nearly 900-atom dendrimer (taking 6.5 s per MD step). PYSEQM also supports batched execution, allowing many geometries or trajectories to be evaluated in a single GPU launch, substantially increasing throughput and making ensemble-based protocols routine. As a demonstration, we compute absorption, emission, and infrared spectra from trajectories propagated on the ground and first excited states. The XL-ESMD scheme yields identical spectra at significantly lower computational cost, establishing the role of extended Lagrangian based dynamics for efficient excited-state BOMD simulations. Beyond raw performance, PYSEQM’s PyTorch foundation provides automatic differentiation for forces, efficient GPU batching, and seamless interfacing with machine learning models. These capabilities position PYSEQM as a practical platform for machine learning-augmented excited-state dynamics and lay the foundation for future data-driven nonadiabatic excited-state dynamics modeling of ultrafast spectroscopic probes.

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

Scaling Analysis of Two–Phase Flow in Fractal Permeability Fields

Fluid mixing in permeable media is essential in many practical applications. The mixing process is a consequence of velocity fluctuations owing to geological heterogeneities and mobility contrast of fluids. Heterogeneities in natural rocks are often spatially correlated, and their properties, such as permeability, may be described using fractal distributions. This work models the fractal characteristics of such permeability fields in which the covariance function is expressed as a power-law function. A generalized scaling relation is derived relating various fractal permeability fields using the magnitude of their fluctuations. Here, this relation reveals the self-similar behavior of two-phase flow in such permeable media. To that end, a recently developed, high-resolution numerical simulator is employed to validate the analytically derived scaling relations. Two flow problems are considered in which flow is governed by 1) a linear, and 2) a nonlinear transport equation. Due to the probabilistic representation of the fractal permeability fields, a sensitivity study is conducted for each flow scenario to determine the number of realizations required for statistical convergence. Scaling analysis is performed using ensemble averages of simulated saturation profiles and their mixing lengths. Results support the validity of the developed scaling relation across the range of investigated flow conditions at intermediate times. The dynamics of linear flow in the asymptotic regime is affected by the correlation structure of heterogeneity. In nonlinear flow, scaling behavior appears to be dominated by the degree of nonlinearity.

58 GEOSCIENCES↗
DOE OSTI2024

Avoiding pitfalls in molecular simulation of vapor sorption: Example of propane and isobutane in metal–organic frameworks for adsorption cooling applications

This study introduces recommendations for conducting molecular simulations of vapor adsorption, with an emphasis on enhancing the accuracy, reproducibility, and comparability of results. The first aspect we address is consistency in the implementation of some details of typical molecular models, including tail corrections and cutoff distances, due to their significant influence on generated data. We highlight the importance of explicitly calculating the saturation pressures at relevant temperatures using methods such as Gibbs ensemble Monte Carlo simulations and illustrate some pitfalls in extrapolating saturation pressures using this method. For grand canonical Monte Carlo (GCMC) simulations, the input fugacity is usually calculated using an equation of state, which often requires the critical parameters of the fluid. We show the importance of using critical parameters derived from the simulation with the same model to ensure internal consistency between the simulated explicit adsorbate phase and the implicit bulk phase in GCMC. We show the advantages of presenting isotherms on a relative pressure scale to facilitate easier comparison among models and with experiment. Extending these guidelines to a practical case study, we evaluate the performance of various isoreticular metal–organic frameworks (MOFs) in adsorption cooling applications. This includes examining the advantages of using propane and isobutane as working fluids and identifying MOFs with a superior performance.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Spectral Gap Filling Method Using the Spectral Information from CLARREO-Pathfinder Measured Earth-Reflected Solar Spectra to Support CPF-CERES Intercalibration

The high-accuracy Earth-reflected solar spectral radiances from CLARREO Pathfinder (CPF) will serve as in-orbit benchmark measurements for intercalibrating other satellite imagers. The CPF team is developing a state-of-the-art intercalibration method to calibrate the shortwave channel (300-5000 nm) of the CERES instrument against the CPF reference with an uncertainty of 0.3%. Because the CPF instrument only covers spectral range from 350 nm to 2300 nm, the team has developed a spectral extension method to fill in the gaps between the spectral domains of CPF and the CERES shortwave channel. The method solely relies on using the hyperspectral information from the CPF wavelengths to predict the radiances below 350 nm and above 2300 nm. To quantify the relationship between the CPF spectrum and the non-overlapping portions of the spectrum between CPF and CERES, the team has simulated a large ensemble of spectra that cover the spectral range for both CPF and CERES. A hyperspectral bidirectional reflectance model for land surface based on thousands of measured land surface reflectance has been developed. This land surface model as well as the Cox-Munk model for ocean surface are used in the spectral gap filling. Various water/ice clouds, aerosols, atmosphere profiles and surfaces were used to simulate Earth-reflected radiances from 200 nm to 15 um under different solar and sensor geometries. The targeted uncertainty for the CPF-extended broadband CERES radiances is less than 0.1%. The previously derived training spectra were updated to replace the Kurucz solar irradiance model with the TSIS-1 Hybrid Solar Reference Spectrum. The influence of the scalar assumption in radiative transfer model equations on the uncertainty in estimating the broadband radiances was also investigated.

Hyperspectral↗
DOE OSTI2021

DART-PFLOTRAN: An ensemble-based data assimilation system for estimating subsurface flow and transport model parameters

Ensemble-based Data Assimilation (EDA), based on the Monte Carlo approach, has been effectively applied to estimate model parameters through inverse modeling in subsurface flow and transport problems. However, implementation of EDA approach involves a complicated workflow that include setting up and executing ensemble forward model simulations, processing observations and model simulation results for parameter updates, and repeat for sequential or iterative EDA. To facilitate the management of such workflow and lower the barriers for adopting EDA-based parameter estimation in subsurface science, we develop a generic software frame-work linking the Data Assimilation Research Testbed (DART) with a massively parallel subsurface FLOw and TRANsport code PFLOTRAN. The new DART-PFLOTRAN leverages both the core data assimilation engines in DART and the computational power afforded by PFLOTRAN. In addition to the standard smoother and filtering options, DART-PFLOTRAN enables an iterative EDA workflow based on the Ensemble Smoother for Multiple Data Assimilation method (ES-MDA) to improve estimation accuracy for nonlinear forward problems. Here, we verify the implementation of ES-MDA in DART-PFLOTRAN using two synthetic cases designed to estimate static permeability and dynamic exchange fluxes across the riverbed, respectively, from continuous temperature measurements made across a depth profile. One-dimensional hydro-thermal simulations are performed in both cases to relate temperature responses with the parameters of interest. In the case of estimating dynamic parameters, we demonstrate the flexibility of DART-PFLOTRAN in automating sequential ES-MDA workflow, which will significantly reduce the time researchers spend on managing complex workflows in similar applications. Both studies yield accurate estimations of the parameters compared to their synthetic truth, while ES-MDA leads to more accurate estimation when a high level of nonlinearity exist between observed responses and unknown parameters. With a code base in Python and Fortran, DART-PFLOTRAN paves the way for applications in large-scale subsurface inverse modeling by automating the complex workflow of sequential ES-MDA that can be executed on various computing platforms.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1997

Evaluation of Tropospheric Water Vapor Simulations from the Atmospheric Model Intercomparison Project

Simulations of humidity from 28 general circulation models for the period 1979-88 from the Atmospheric Model Intercomparison Project are compared with observations from radiosondes over North America and the globe and with satellite microwave observations over the Pacific basin. The simulations of decadal mean values of precipitable water (W) integrated over each of these regions tend to be less moist than the real atmosphere in all three cases; the median model values are approximately 5% less than the observed values. The spread among the simulations is larger over regions of high terrain, which suggests that differences in methods of resolving topographic features are important. The mean elevation of the North American continent is substantially higher in the models than is observed, which may contribute to the overall dry bias of the models over that area. The authors do not find a clear association between the mean topography of a model and its mean W simulation, however, which suggests that the bias over land is not purely a matter of orography. The seasonal cycle of W is reasonably well simulated by the models, although over North America they have a tendency to become moister more quickly in the spring than is observed. The interannual component of the variability of W is not well captured by the models over North America. Globally, the simulated W values show a signal correlated with the Southern Oscillation index but the observations do not. This discrepancy may be related to deficiencies in the radiosonde network, which does not sample the tropical ocean regions well. Overall, the interannual variability of W, as well as its climatology and mean seasonal cycle, are better described by the median of the 28 simulations than by individual members of the ensemble. Tests to learn whether simulated precipitable water, evaporation, and precipitation values may be related to aspects of model formulation yield few clear signals, although the authors find, for example, a tendency for the few models that predict boundary layer depth to have large values of evaporation and precipitation. Controlled experiments, in which aspects of model architecture are systematically varied within individual models, may be necessary to elucidate whether and how model characteristics influence simulations.

Gaffen, Dian J.↗
DOE OSTI2022

Impacts of Lake Surface Temperature on the Summer Climate Over the Great Lakes Region

The surface of the Great Lakes interacts with the atmosphere, influencing the weather and climate over the Great Lakes Region (GLR). However, most climate models were not designed with sufficient emphasis on lake–atmosphere interactions, which could potentially cause model biases over the GLR. To understand how lake surface temperature (LST) affects the regional summer climate over the GLR, we conducted twin experiments using the Weather Research and Forecasting model at a spatial resolution of 4 km using two different LST data sets as the bottom boundary condition over the Great Lakes. In this work, our simulations include 10 ensemble members for the summer of 2018 and a single multiyear run for the summers of 2014–2020. Results show that variations in LST influence atmospheric temperature and moisture at a local scale—while affecting the convective environment and precipitation processes over a much larger spatial scale. In particular, an LST that is only 1°C–3°C warmer (depending on the lake) increases near-surface air temperature by 1.93°C and 0.97°C over Lake Superior and Lake Erie, respectively, and increases evaporation over the lakes by 0.23 and 1.1 mm day –1 . The warmer LST reduces mesoscale convective precipitation upstream of the GLR; however, it increases isolated deep convective precipitation and nonconvective precipitation downstream of GLR due to increased local instability and enhancement of moisture transport. Our analyses confirm the robustness of these impacts, which is at least 2 times larger than the model internal variability and is seen across all simulated summer seasons.

54 ENVIRONMENTAL SCIENCES↗
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↗