Search NASA⌕ Search

SEARCH · Search NASA

Results for “Large ensemble”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Scalable Cloning on Large-Scale GPU Platforms with Application to Time-Stepped Simulations on Grids

Cloning is a technique to efficiently simulate a tree of multiple what-if scenarios that are unraveled during the course of a base simulation. However, cloned execution is highly challenging to realize on large, distributed memory computing platforms, due to the dynamic nature of the computational load across clones, and due to the complex dependencies spanning the clone tree. In this paper, we present the conceptual simulation framework, algorithmic foundations, and runtime interface of CloneX, a new system we designed for scalable simulation cloning. It efficiently and dynamically creates whole logical copies of a dynamic tree of simulations across a large parallel system without full physical duplication of computation and memory. The performance of a prototype implementation executed on up to 1,024 graphical processing units of a supercomputing system has been evaluated with three benchmarks—heat diffusion, forest fire, and disease propagation models—delivering a speed up of over two orders of magnitude compared to replicated runs. Finally, the results demonstrate a significantly faster and scalable way to execute many what-if scenario ensembles of large simulations via cloning using the CloneX interface.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

State-to-state, multi-collision, energy transfer in H–H{sub 2} gas ensembles

We use our recently developed computational model of energy flow in gas ensembles to study translation-to-internal energy conversion in an ensemble consisting of H{sub 2}(0; 0) in a bath of H atoms. This mixture is found in plasmas of industrial importance and also in interstellar clouds. The storage of energy of relative motion as rovibrational energy of H{sub 2} represents a potential mechanism for cooling translation. This may have relevance in astrophysical contexts such as the post-recombination epoch of the early universe when hydrogenic species dominated and cooling was a precondition for the formation of structured objects. We find that conversion of translational motion to H{sub 2} vibration and rotation is fast and, in our closed system, is complete within around 100 cycles of ensemble collisions. Large amounts of energy become stored as H{sub 2} vibration and a tentative mechanism for this unequal energy distribution is suggested. The “structured dis-equilibrium” we observe is found to persist through many collision cycles. In contrast to the rapidity of excitation, the relaxation of H{sub 2}(6; 10) in H is very slow and not complete after 10{sup 5} collision cycles. The quasi-equilibrium modal temperatures of translation, rotation, and vibration are found to scale linearly with collision energy but at different rates. This may be useful in estimating the partitioning of energy within a given H + H{sub 2} ensemble.

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

Bridging Hydrological Ensemble Simulation and Learning Using Deep Neural Operators

Ensemble-based simulation and learning (ESnL) has long been used in hydrology for parameter inference, but computational demands of process-based ESnL can be quite high. To address this issue, we propose a deep neural operator learning approach. Neural operators are generic machine learning algorithms that can learn functional mappings between infinite-dimensional spaces, providing a highly flexible tool for scientific machine learning. Our approach is built upon DeepONet, a specific deep neural operator, and is designed to address several common problems in hydrology, namely, model parameter estimation, prediction at ungaged locations, and uncertainty quantification. Here we demonstrate the effectiveness of our DeepONet-based workflow using an existing large model ensemble created for an eastern U.S. watershed that is instrumented with 10 streamflow gages. Results suggest DeepONet achieves high efficiency in learning an ML surrogate model from the model ensemble, with the modified Kling-Gupta Efficiency exceeding 0.9 on holdout test sets. Parameter inference, carried out using the trained DeepONet surrogate model and genetic algorithm, also yields robust results. Additionally, we formulate and train a separate DeepONet model for physics-informed, seq-to-seq streamflow forecasting, which further reduces biases in the pre-trained DeepONet surrogate model. While this study focuses primarily on a single watershed, our approach is general and may be extended to enable learning from model ensembles across multiple basins or models. Thus, this research represents a significant contribution to the application of hybrid machine learning in hydrology.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Calibration and Uncertainty Estimation Using the Ensemble Kalman Filter with a Large Subsurface Flow and Transport Model - 20321

At routinely monitored groundwater contamination sites, periodically measured environmental conditions such as groundwater levels and contaminant concentrations are used to inform and confirm a conceptual site model (CSM) and guide the development and calibration of a numerical groundwater flow and transport model. The calibration of groundwater flow and transport models after each measurement (sampling) event can illuminate deficiencies in a CSM, identify areas where additional monitoring is warranted, and predict the behavior of the system to guide decision making. However, manual and automated (e.g. PEST) model calibration tools can be time-consuming and computationally expensive to implement after each sampling event. Perhaps as a result, such calibration tools generally utilize all available monitoring data simultaneously rather than sequentially assimilating monitoring data one sampling event at a time as the results from sampling become available. A more real-time data assimilation approach may reduce parameter uncertainty, quantify the value of additional monitoring data, and produce a usable model more quickly and with less effort. To mitigate the potential time-consuming aspects of manual and widely applied automated calibration techniques, a data assimilation algorithm called the ensemble Kalman filter (EnKF) was evaluated as a relatively efficient method of model calibration and uncertainty assessment via the sequential integration of monitoring data into a model. The EnKF was able to successfully and efficiently assimilate monitoring and modeling data to calibrate a complex flow and transport model at a real-world site with significant subsurface heterogeneity, uncertainty, and 12 years of monitoring data (over 4,000 individual measurements of groundwater levels and over 2,500 measurements of contaminant concentrations). Starting with an uncalibrated model data from annual sampling events were sequentially assimilated, and the resultant predication errors and estimated parameter uncertainties were tracked. After all monitoring data were assimilated, both flow and transport residuals at the end of the EnKF process were comparable to those produced via a concurrent PEST calibration effort but required fewer model simulations. Both uncertainty and prediction errors decreased over time. In a real-time application, the adequacy of the model could be assessed after each sampling event. The benefits of such a real-time approach to utilizing monitoring data include reduced costs (in the form of model updates or site characterization efforts), early flagging of possible errors in the CSM, and a reduced risk of overfitting and corresponding increased confidence in model predictions. This tool may be particularly useful compared to other calibration techniques (e.g. manual, PEST) when model runtimes are long, calibration parameters are many, or parameter uncertainty is large. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2018

The multi-assumption architecture and testbed (MAAT v1.0): R code for generating ensembles with dynamic model structure and analysis of epistemic uncertainty from multiple sources

Abstract. Here, computer models are ubiquitous tools used to represent systems across many scientific and engineering domains. For any given system, many computer models exist, each built on different assumptions and demonstrating variability in the ways in which these systems can be represented. This variability is known as epistemic uncertainty, i.e. uncertainty in our knowledge of how these systems operate. Two primary sources of epistemic uncertainty are (1) uncertain parameter values and (2) uncertain mathematical representations of the processes that comprise the system. Many formal methods exist to analyse parameter-based epistemic uncertainty, while process-representation-based epistemic uncertainty is often analysed post hoc, incompletely, informally, or is ignored. In this model description paper we present the multi-assumption architecture and testbed (MAAT v1.0) designed to formally and completely analyse process-representation-based epistemic uncertainty. MAAT is a modular modelling code that can simply and efficiently vary model structure (process representation), allowing for the generation and running of large model ensembles that vary in process representation, parameters, parameter values, and environmental conditions during a single execution of the code. MAAT v1.0 approaches epistemic uncertainty through sensitivity analysis, assigning variability in model output to processes (process representation and parameters) or to individual parameters. In this model description paper we describe MAAT and, by using a simple groundwater model example, verify that the sensitivity analysis algorithms have been correctly implemented. The main system model currently coded in MAAT is a unified, leaf-scale enzyme kinetic model of C 3 photosynthesis. In the Appendix we describe the photosynthesis model and the unification of multiple representations of photosynthetic processes. The numerical solution to leaf-scale photosynthesis is verified and examples of process variability in temperature response functions are provided. For rapid application to new systems, the MAAT algorithms for efficient variation of model structure and sensitivity analysis are agnostic of the specific system model employed. Therefore MAAT provides a tool for the development of novel or toy models in many domains, i.e. not only photosynthesis, facilitating rapid informal and formal comparison of alternative modelling approaches.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Third-Party Supplier Risk Re-Classification Using Multi-Model Semantic Voting and External Web Augmentation

Risk decisions in many third-party risk management (TPRM) workflows rely on static inherent risk questionnaires (IRQ). These static forms provide a snapshot of the vendor from the business users’ perspective, as these requests are processed without cross-referencing for evidence. Consequently, responses can be misinformed or embellished with inaccuracies, thereby masking the vendor’s true risk to the enterprise. This paper presents a multi-stage verification framework to augment IRQs with web evidence and a deterministic ensemble of large language model assessors to reclassify risk. In a case study of 100 submissions previously misclassified as low risk, the proposed framework correctly identified 76% of the cases as high risk, while the existing workflow identified none. McNemar’s continuity corrected statistics of 74 were obtained with a two sided p-value of 2.65 × 10-23, indicating a significantly more effective workflow compared to the legacy model.

99 - GENERAL AND MISCELLANEOUS↗
NASA NTRS2017

Classifying Multi-Model Wheat Yield Impact Response Surfaces Showing Sensitivity to Temperature and Precipitation Change

Crop growth simulation models can differ greatly in their treatment of key processes and hence in their response to environmental conditions. Here, we used an ensemble of 26 process-based wheat models applied at sites across a European transect to compare their sensitivity to changes in temperature (minus 2 to plus 9 degrees Centigrade) and precipitation (minus 50 to plus 50 percent). Model results were analysed by plotting them as impact response surfaces (IRSs), classifying the IRS patterns of individual model simulations, describing these classes and analysing factors that may explain the major differences in model responses. The model ensemble was used to simulate yields of winter and spring wheat at four sites in Finland, Germany and Spain. Results were plotted as IRSs that show changes in yields relative to the baseline with respect to temperature and precipitation. IRSs of 30-year means and selected extreme years were classified using two approaches describing their pattern. The expert diagnostic approach (EDA) combines two aspects of IRS patterns: location of the maximum yield (nine classes) and strength of the yield response with respect to climate (four classes), resulting in a total of 36 combined classes defined using criteria pre-specified by experts. The statistical diagnostic approach (SDA) groups IRSs by comparing their pattern and magnitude, without attempting to interpret these features. It applies a hierarchical clustering method, grouping response patterns using a distance metric that combines the spatial correlation and Euclidian distance between IRS pairs. The two approaches were used to investigate whether different patterns of yield response could be related to different properties of the crop models, specifically their genealogy, calibration and process description. Although no single model property across a large model ensemble was found to explain the integrated yield response to temperature and precipitation perturbations, the application of the EDA and SDA approaches revealed their capability to distinguish: (i) stronger yield responses to precipitation for winter wheat than spring wheat; (ii) differing strengths of response to climate changes for years with anomalous weather conditions compared to period-average conditions; (iii) the influence of site conditions on yield patterns; (iv) similarities in IRS patterns among models with related genealogy; (v) similarities in IRS patterns for models with simpler process descriptions of root growth and water uptake compared to those with more complex descriptions; and (vi) a closer correspondence of IRS patterns in models using partitioning schemes to represent yield formation than in those using a harvest index. Such results can inform future crop modelling studies that seek to exploit the diversity of multi-model ensembles, by distinguishing ensemble members that span a wide range of responses as well as those that display implausible behaviour or strong mutual similarities.

Classification; Climate change; Crop model; Ensemb↗
DOE OSTI2017

Deterministic skill of ENSO predictions from the North American Multimodel Ensemble

Hindcasts and real-time predictions of the east-central tropical Pacific sea surface temperature (SST) from the North American Multimodel Ensemble (NMME) system are verified for 1982–2015. Skill is examined using two deterministic verification measures: mean squared error skill score (MSESS) and anomaly correlation. Verification of eight individual models shows somewhat differing skills among them, with some models consistently producing more successful predictions than others. The skill levels of MME predictions are approximately the same as the two best performing individual models, and sometimes exceed both of them. A decomposition of the MSESS indicates the presence of calibration errors in some of the models. In particular, the amplitudes of some model predictions are too high when predictability is limited by the northern spring ENSO predictability barrier and/or when the interannual variability of the SST is near its seasonal minimum. The skill of the NMME system is compared to that of the MME from the IRI/CPC ENSO prediction plume, both for a comparable hindcast period and also for a set of real-time predictions spanning 2002–2011. Comparisons are made both between the MME predictions of each model group, and between the average of the skills of the respective individual models in each group. Acknowledging a hindcast versus real-time inconsistency in the 2002–2012 skill comparison, the skill of the NMME is slightly higher than that of the prediction plume models in all cases. This result reflects well on the NMME system, with its large total ensemble size and opportunity for possible complementary contributions to skill.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

High-Resolution Momentum Imaging—From Stern’s Molecular Beam Method to the COLTRIMS Reaction Microscope [Book Chapter]

Multi-particle momentum imaging experiments are now capable of providing detailed information on the properties and the dynamics of quantum systems in Atomic, Molecular and Photon (AMO) physics. Historically, Otto Stern can be considered the pioneer of high-resolution momentum measurements of particles moving in a vacuum and he was the first to obtain sub-atomic unit (a.u.) momentum resolution (Schmidt-Böcking et al. in The precision limits in a single-event quantum measurement of electron momentum and position, these proceedings). A major contribution to modern experimental atomic and molecular physics was his so-called molecular beam method, which Stern developed and employed in his experiments. With this method he discovered several fundamental properties of atoms, molecules and nuclei. As corresponding particle detection techniques were lacking during his time, he was only able to observe the averaged footprints of large particle ensembles. Today it is routinely possible to measure the momenta of single particles, because of the tremendous progress in single particle detection and data acquisition electronics. A “state-of-the-art” COLTRIMS reaction microscope can measure, for example, the momenta of several particles ejected in the same quantum process in coincidence with sub-a.u. momentum resolution. Such setups can be used to visualize the dynamics of quantum reactions and image the entangled motion of electrons inside atoms and molecules. This review will briefly summarize Stern’s work and then present in longer detail the historic steps of the development of the COLTRIMS reaction microscope. Furthermore, some benchmark results are shown which initially paved the way for a broad acceptance of the COLTRIMS approach. Finally, a small selection of milestone work is presented which has been performed during the last two decades.

Jahnke, T.↗
DOE OSTI2016

Benefits of mitigation for future heat extremes under RCP4.5 compared to RCP8.5

We investigate changes in statistics of extreme daily temperature using ensembles from the Community Earth System Model (CESM) under a high and a lower emission scenarios, . The ensembles provide large samples for a robust application of extreme value theory. We estimate return values and return periods for annual maxima of the daily high and low temperatures as well as the 3-day averages of the same variables in current and future climate. Results indicate statistically significant increases (compared to the reference period of 1996–2005) in extreme temperatures over all land areas as early as 2025 under both scenarios, with statistically significant differences between them becoming pervasive over the globe by 2050. The substantially smaller changes, for all indices, produced under the lower emission case translate into sizeable benefits from emission mitigation: By 2075, in terms of reduced changes in 1-day heat extremes, about 95 % of land regions would see benefits of 1 °C or more under the lower emissions scenario, and 50 % or more of the land areas would benefit by at least 2 °C. 6 % of the land area would benefit by 3 °C or more in projected extreme minimum temperatures and 13 % would benefit by this amount for extreme maximum temperature. Benefits for 3-day metrics are similar. The future frequency of current extremes is also greatly reduced by mitigation: by the end of the century, under RCP8.5 more than half the land area experiences the current 20-year events every year while only between about 10 and 25 % of the area is affected by such severe changes under RCP4.5.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

High precision statistical Landau gauge lattice gluon propagator computation vs. the Gribov–Zwanziger approach

In this article we report on lattice results for the Landau gauge gluon propagator using large statistical ensembles. In particular, the compatibility of the lattice data with the tree level predictions of the Refined Gribov–Zwanziger and the Very Refined Gribov–Zwanziger actions is analysed, complementing earlier work using small-scale statistics. Our results show that the data is well described by the tree level estimate only up to momenta p≲1 GeV and clearly favours the so-called Refined Gribov–Zwanziger scenario propagator. For the Very Refined Gribov–Zwanziger action, these results imply a particular relation between d=2 condensates. Furthermore, we also provide a global fit of the lattice data that interpolates between the above scenario at low momenta and the usual continuum one-loop renormalization improved perturbation theory after introducing an infrared log-regularizing term.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Structure and reactivity of Pt–In intermetallic alloy nanoparticles: Highly selective catalysts for ethane dehydrogenation

The structure of silica supported Pt and Pt-In bimetallic catalysts with nominal In: Pt atomic ratios of 0.7 and 1.4 were here determined by in situ synchrotron XAS and XRD. It was seen that the addition of In led to the formation of two different intermetallic alloy phases. At an In: Pt ratio of 0.7 the Pt 3 In phase with a Cu 3 Au structure was formed. When the ratio was increased to 1.4 a shell of PtIn 2 having a CaF 2 structure formed around a core of Pt 3 In. The catalysts were tested for ethane dehydrogenation at 600° to determine the effect of alloying on ethylene selectivity and turnover rate (TOR). The monometallic Pt catalysts was 73% selective for ethylene and had an initial TOR of 0.7 s -1 . Both alloy catalysts were approximate to 100% selective for dehydrogenation and had higher initial TOR, 2.8 s -1 and 1.6 s -1 for In: Pt ratio of 0.7 and 1.4, respectively. The increase in selectivity is attributed to the elimination of large Pt ensembles resulting from geometric changes to the catalyst surface upon alloying. Electronic changes due to the formation of Pt-In bonds are thought to be responsible for the increases in TOR in the alloy catalysts.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A theoretical calculation of stacking fault energy of Ni alloys: The effects of temperature and composition

Combining cluster expansion (CE) method with one dimensional axial Ising model, this work investigated the effects of alloying elements and configurational variations due to temperature on the stacking fault energy (SFE) of FCC Ni binary alloys. Ensembles of large numbers of atomistic structures, each with more than ~400 atoms, were generated to consider sufficient long-range chemical disorder and temperature effects due to configurational entropy. Here, a Monte Carlo Metropolis algorithm was used to generate these structures, whose energies were then evaluated based on the effective cluster interactions obtained from CE. As a baseline, this work had shown the SFE of pure Ni and Al to be 127 mJ/m 2 and 137 mJ/m 2 , respectively, which agreed with the experimental values of 125 mJ/m 2 and 150 mJ/m 2 reported in the literature. Additions of Al, Ti, Cr and Co to pure Ni were found to decrease the SFE to different extents. Although temperature does not strongly influence the SFE of the FCC Ni-Al and Ni-Cr binary alloys, it can lead to significant changes to the SFE of the FCC Ni-Ti and Ni-Co alloys. While effects of temperature and composition on SFE observed in this work were calculated from binary Ni alloys, the general trends are nonetheless expected to be valid in the γ phase of multicomponent Ni superalloys.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Scalable tuning of building models to hourly data

Energy models of existing buildings are unreliable unless calibrated so they correlate well with actual energy usage. Manual tuning requires a skilled professional, is prohibitively expensive for small projects, imperfect, non-repeatable, non-transferable, and not scalable to the dozens of sensor channels that smart meters, smart appliances, and cheap/ubiquitous sensors are beginning to make available today. A scalable, automated methodology is needed to quickly and intelligently calibrate building energy models to all available data, increase the usefulness of those models, and facilitate speed-and-scale penetration of simulation-based capabilities into the marketplace for actualized energy savings. The "Autotune'' project is a novel, model-agnostic methodology which leverages supercomputing, large simulation ensembles, and big data mining with multiple machine learning algorithms to allow automatic calibration of simulations that match measured experimental data in a way that is deployable on commodity hardware. This paper shares several methodologies employed to reduce the combinatorial complexity to a computationally tractable search problem for hundreds of input parameters. Furthermore, accuracy metrics are provided which quantify model error to measured data for either monthly or hourly electrical usage from a highly-instrumented, emulated-occupancy research home.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2017

Mechanisms of Lipid Scrambling by the G Protein-Coupled Receptor Opsin

Several class-A G protein-coupled receptor (GPCR) proteins act as constitutive phospholipid scramblases catalyzing the transbilayer translocation of >10,000 phospholipids per second when reconstituted into synthetic vesicles. To address the molecular mechanism by which these proteins facilitate rapid lipid scrambling, we carried out large-scale ensemble atomistic molecular dynamics simulations of the opsin GPCR. We report here that, in the process of scrambling, lipid head groups traverse a dynamically revealed hydrophilic pathway in the region between transmembrane helices 6 and 7 of the protein while their hydrophobic tails remain in the bilayer environment. We present quantitative kinetic models of the translocation process based on Markov State Model analysis. As key residues on the lipid translocation pathway are conserved within the class-A GPCR family, our results illuminate unique aspects of GPCR structure and dynamics while providing a rigorous basis for the design of variants of these proteins with defined scramblase activity.

59 BASIC BIOLOGICAL SCIENCES↗
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 OSTI2019

Identification of a Pt 3 Co Surface Intermetallic Alloy in Pt–Co Propane Dehydrogenation Catalysts

Bimetallic Pt-Co nanoparticles (NPs) were prepared and characterized by scanning transmission electron microscopy, in situ X-ray absorption spectroscopy, in situ synchrotron X-ray diffraction, and catalytic conversion for propane dehydrogenation with and without added H-2. In addition, the surface extended X-ray absorption fine structure (EXAFS) obtained by fitting the difference spectrum between the fully reduced and room-temperature-oxidized catalysts suggest that the surface structure remains Pt 3 Co, although the core changes from Pt to Pt 3 Co and to PtCo. At low Co loading, the bimetallic nanoparticles form a Pt3Co intermetallic surface alloy with Pt-rich core. With increasing Co loading, a full alloy forms where both the surface and NP compositions are Pt 3 Co. A further increase in Co loading leads to a Co-rich NP core, likely PtCo, with a surface of Pt 3 Co. Although Pt-Co intermetallic alloys form two different phases and several morphologies, the surface structures are similar in all catalysts. Although both monometallic Pt and Co are active for alkane dehydrogenation, all bimetallic Pt-Co catalysts are significantly more olefin selective than either single metal. The turnover rates of the bimetallic catalysts indicate that Pt is the active atom with little contribution from Co atoms. The high olefin selectivity is suggested to be due to Co acting as a less active structural promoter to break up large Pt ensembles in bimetallic NPs.

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

Accelerating and Parallelizing Lagrangian Simulations of Mixing-Limited Reactive Transport

Recent advances in random walk particle tracking have enabled direct simulation of mixing and reactions by allowing the particles to interact with each other using a multipoint mass transfer scheme. The mass transfer scheme allows separation of mixing and spreading processes, among other advantages, but it is computationally expensive because its speed depends on the number of interacting particle pairs. This note explores methods for relieving the computational bottleneck caused by the mass transfer step, and we use these algorithms to develop a new parallel, interacting particle model. The new model is a combination of a sparse search algorithm and a novel domain decomposition scheme, both of which offer significant speedup relative to the reference case—even when they are executed serially. We combine the strengths of these methods to create a parallel particle scheme that is highly accurate and efficient with run times that scale as 1/P for a fixed number of particles, where P is the number of computational cores (equivalently, subdomains, in this work) being used. Here, the new parallel model is a significant advance because it enables efficient simulation of large particle ensembles that are needed for environmental simulations and also because it can naturally pair with parallel geochemical solvers to create a practical Lagrangian tool for simulating mixing and reactions in complex chemical systems.

54 ENVIRONMENTAL SCIENCES↗