Results for “parallel systems” — Page 79

DOE OSTI2024

Proposed crystal structure of carbadox, C 11 H 10 N 4 O 4

A model for the crystal structure of carbadox has been generated and refined using synchrotron X-ray powder diffraction data and optimized using density functional theory techniques. Carbadox crystallizes in space groupP2 1 (#4) witha= 13.8155(3),b= 21.4662(1),c= 16.3297(3) Å,β= 110.0931(7)°,V= 4548.10(3) Å 3 , andZ= 16. The crystal structure is characterized by approximately parallel stacking of the eight independent carbadox molecules parallel to thebc-plane. There are two different molecular configurations of the eight carbadox molecules; five are in the lower-energy configuration and three are in a ~10% higher-energy configuration. This arrangement likely achieves the lowest-energy crystalline packing via hydrogen bonding. Hydrogen bonds link the molecules both within and between the planes. Each of the amino groups forms a N–H⋯O hydrogen bond to an oxygen atom of the 1,4-dioxidoquinoxaline ring system of another molecule. The result is four pairs of hydrogen-bonded molecules, which form rings with graph setR2,2(14). Variation in specimen preparation can affect the preferred orientation of particles considerably. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®).

Materials Science↗

DOE OSTI2021

Toward active disruption avoidance via real-time estimation of the safe operating region and disruption proximity in tokamaks

This paper describes a real-time capable algorithm for identifying the safe operating region around a tokamak operating point. The region is defined by a convex set of linear constraints, from which the distance of a point from a disruptive boundary can be calculated. The disruptivity of points is calculated from an empirical machine learning predictor that generates the likelihood of disruption. While the likelihood generated by such empirical models can be compared to a threshold to trigger a disruption mitigation system, the safe operating region calculation enables active optimization of the operating point to maintain a safe margin from disruptive boundaries. The proposed algorithm is tested using a random forest disruption predictor fit on data from DIII-D. The safe operating region identification algorithm is applied to historical data from DIII-D showing the evolution of disruptive boundaries and the potential impact of optimization of the operating point. Real-time relevant execution times are made possible by parallelizing many of the calculation steps and implementing the algorithm on a graphics processing unit. Lastly, a real-time capable algorithm for optimizing the target operating point within the identified constraints is also proposed and simulated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Carbon Cycling, Environmental & Rural Economic Impacts of Collecting & Processing Specific Woody Feedstocks in Biofuels

Woody biomass will be an essential feedstock for a large-scale cellulosic biofuel industry. The life cycle carbon accounting for the production of biofuels from woody feedstocks is complex and has engendered significant controversy. Issues such as below ground carbon, carbon debt, varied regional forest practices, multiple parallel forest product lines, and development of realistic counterfactual scenarios all contribute to the complexity. DOE funded CORRIM to develop comprehensive and definitive lifecycle inventories and assessments on the production of fuels from woody feedstocks. CORRIM brought together expertise in forest practices, short rotation woody crops, bioconversion of woody biomass, process modeling, and life cycle assessment to successfully accomplish this goal. The CORRIM team has developed data on forest productivity, fuel usage, and fuel production for six regionally specific forest systems. These six forest systems include three current commercial systems; southern pine plantations, Douglas-fir plantations, naturally regenerated Northeastern (NE) spruce/fir, and three ‘short rotation woody crops’, poplar, eucalyptus and willow, which are at different stages of demonstration in the US. The fuel production systems include cellulosic ethanol and bio-oil based hydrocarbons. The project was successfully reviewed at the BETO Program review in March 2015, 2017, and 2019 and has resulted currently in 19 publications and reports and 35 presentations at both national and international conferences, with more in the pipeline. See publications and presentations for links to each document. Finally, and most importantly, CORRIM has gone beyond the scope of the original DOE proposal to work closely with GREET at Argonne National Laboratory (ANL) to incorporate all the life cycle data and scenario models into their modeling system. GREET is the most widely used and definitive information source for evaluating lifecycle carbon emissions for fuels. The incorporation of CORRIM data from this project guarantees the results of the research will be extensively used and widely disseminated. Technical process improvements and policy relevant accomplishments are detailed in the relevant programmatic sections in the full report including citations therein. Highlights are summarized here for easy reference.

09 BIOMASS FUELS↗

DOE OSTI2020

Optimizing the Accelerated Recursive Doubling Algorithm for Block Tridiagonal Systems of Equations

The need to solve block tridiagonal systems with hundreds or thousands of right-hand sides for the same block tridiagonal matrix is common in a variety of disciplines. To meet this need, the Accelerated Recursive Doubling Algorithm was developed. After a right-hand side independent phase, the algorithm allows for the quick, online calculation of solutions for different right-hand sides. In this work, we present methods to optimize the Accelerated Recursive Doubling Algorithm in memory usage and computation time in a hybrid parallelization model. The right-hand side independent phase of the naïve implementation takes ≥ 11/3 the amount of memory required to store the tridiagonal matrix, while our implementation reduces the fraction to ≈ 5/3 . The right-hand side dependent phase of the naïve implementation takes ≥ 6 times the amount of memory required to store the right-hand side, while our implementation reduces the fraction to ≈ 3. The computation time for the independent phase is reduced to ≈ 2/3 times that of the naïve implementation, while the computation time for the dependent phase is reduced to ≈ 5/9 . With increasing numbers of shared-memory threads q on every distributed processing element, we have O(q) theoretical speedup.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Assessment of Model Hub Height Wind Speed Performance Using DOE Lidar Buoy Data

As the offshore wind industry continues to gain momentum both in the U.S. and globally, reanalysis products remain essential tools for the lifecycle of a wind farm. From project siting and production estimates to determining construction and maintenance windows, reanalysis products play an important role. These regional and global long-term models also serve as inputs to higher resolution models. With such wide-ranging impacts on the offshore wind industry, the need for validation of reanalysis products is strong, particularly at typical turbine hub heights. Pacific Northwest National Laboratory (PNNL) operates two AXYS WindSentinelTM buoys for the U.S. Department of Energy (DOE) in order to collect meteorological and oceanographic data, including hub height wind speed, in areas of interest for offshore wind development. This work uses DOE buoy observations off the coasts of New Jersey and Virginia to evaluate the performance of Modern Era Retrospective Analysis for Research and Applications-2 (MERRA-2), the North American Regional Reanalysis (NARR), the European Center for Medium-Range Weather Forecasts (ECMWF) Reanalysis 5th generation (ERA5), and the National Oceanic and Atmospheric Administration (NOAA) assimilation system Rapid Refresh (RAP). Biases and degrees of correspondence are determined for each reanalysis product in order to provide insights on the performance and uncertainty for long-term wind resource characterization. With the baseline metrics determined, this work proceeds with an investigation into the sources of large deviations between the modeled and observed hub height wind speeds during the East Coast lidar buoy deployments. Consistent sources of reanalysis model error include stably stratified flow conditions, high wind shear, flow parallel to the coastline, and tropical and winter storms. For the near shore New Jersey location, offshore winds (i.e. blowing from land to water) tend to be associated with model overestimation of observed hub height wind speed, while winds parallel to the coastline are correlated with model underestimation of wind speed. Similarly, near the Virginia buoy, during southerly winds (parallel to the coastline) models underestimate the observed wind speed due to coastal upwelling. Reanalysis model wind speeds underestimate the hub height winds with increasing stable atmospheric stratification. Reanalysis model biases are seen when the observed wind speeds approach or exceed typical turbine cut-out speeds (around 20 m s-1). Large model errors during winter and tropical storm events, high wind speed events, and ramp events are often the result of wind speed timing and magnitude offsets. Overall, ERA5 provides the most successful representation of observed offshore hub height wind speeds at the U.S. East Coast buoy locations and is therefore best suited as an input boundary condition for model case studies. ERA5 performs well for a variety of atmospheric phenomena, including storms and sea breezes, however RAP, the next most successful reanalysis model, performs best in capturing the frequency of ramp events.

17 WIND ENERGY↗

DOE OSTI2023

Robust Decentralized Secondary Control Scheme for Inverter-based Power Networks

Inverter-dominated microgrids are quickly becoming a key building block of future power systems. They rely on centralized controllers that can provide reliability and resiliency in extreme events. Nonetheless, communication failures due to cyber-physical attacks or natural disasters can make autonomous operation of islanded microgrids challenging. This paper examines a unified decentralized secondary control scheme that is robust to inverter clock synchronization errors and can be seamlessly applied to grid-following or grid-forming control architectures. The proposed scheme overcomes the well-known stability problem that arises from parallel operation of local integral controllers. Theoretical guarantees for stability are provided along with criteria to appropriately tune the secondary control gains to achieve good frequency regulation performance while ensuring fair power sharing. The efficacy of our approach is demonstrated through simulations on a 5-bus microgrid with four grid-forming inverters.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2023

Editorial: Advanced water splitting technologies development: Best practices and protocols

As the level of deployment and utilization of renewable energy sources, including wind and solar, continues to rise, large-scale, long-term energy storage technologies that could accommodate weekly and seasonal energy fluctuations will play a significant role in the overall deployment of renewable energies in the future. Harnessing and storing renewable energy resources via electrochemical, photoelectrochemical, or thermochemical processes by converting renewable energy into sustainable (energy storage) fuels have the potential to meet the long-term, terawatt scale energy storage challenge. Renewable hydrogen production is the cornerstone for sustainable fuel production and deep decarbonization of multiple sectors in our society. Cost-competitive clean hydrogen provides value to applications, such as 1) in the transportation sector for fuel cell vehicles, 2) in the electric grid sector for system stability and load balancing, and 3) in the industrial sector with metal refineries, cement production, and biomass upgrading (carbon-free fertilizer production). In addition, coupling clean renewable hydrogen with the carbon and nitrogen cycles enables known and well-established thermal-chemical processes to generate renewable hydrocarbon fuels and ammonia. The Advanced Water Splitting Technologies (AWST): low temperature electrolysis (LTE), high temperature electrolysis (HTE), photoelectrochemical (PEC) and solar thermo-chemical hydrogen (STCH) provide four unique and parallel approaches to produce low cost, low greenhouse gas (GHG) emission hydrogen at scale (Figure 1). Cost competitive clean hydrogen production using these four technologies is a current high priority focus for governments and industry. In June of 2022, the U.S. Department of Energy (DOE) launched the first in a series of Earthshot Initiatives. The Hydrogen Shot, “1 1 1” aims to reduce the cost of clean hydrogen by more than 80% to one dollar per one kilogram in 1 decade ($\$$1/kg H 2 ). The European Green Deal and the International Energy Agency (IEA) have implemented a strong focus on green hydrogen production for a clean and secure energy future.

benchmarking, low temperature electrolysis↗

DOE OSTI2022

Pele: An Exascale-Ready Suite of Combustion Codes

High fidelity simulations of realistic combustion devices are extremely demanding computationally because of the requirements to capture complex fuel chemical decomposition, its intricate interactions with turbulent, often multiphase, flows, and the wide separation of space and time scales between the thin flame and the device boundaries. Software required to carry out such computations tends to be extremely complex, particularly when designed to exploit hardware accelerators, and can be difficult to port and maintain. We present Pele, a performance portable suite of tools for the simulation of combustion systems, including codes to evolve reactive multiphase configurations in the low Mach number and compressible flow regimes, along with a set of inter-compatible post processing and in situ analysis tools. The Pele suite of tools is built on top of the AMReX framework for block-structured adaptive mesh refinement, which provides efficient data structures and algorithms that enable the development of a wide variety of efficient mesh and particle based PDE integration schemes. A hierarchical MPI+X parallelism scheme supports CPU-only and accelerated architectures, where X can be OpenMP, CUDA, and HIP based approaches for intra-node computational work distribution. The algorithms and data structures underlying the Pele simulation and analysis tools are highly scalable and performant across a wide variety of high-performance computing platforms, including DOEs newest exascale-class machines, Frontier and Aurora. The simulation and analysis tools are fully documented and freely distributed as open source via GitHub. We present key algorithmic and software challenges, solution strategies, performance and resulting set of capabilities.

AMReX↗

DOE OSTI2023

Towards a Circular Economy for PET Bottles in the US

The United States generates the most plastic waste of any country. Along with that GHG emissions from the global plastic economy are expected to increase to 15% of the global carbon budget by 2050. It is imperative that plastic recycling is made a reality to reduce both plastic pollution in the environment and GHG emissions. A portfolio of end-of-life strategies must be implemented to minimize environmental impacts and retain valuable plastic material, but it is challenging to compare options that generate products with different utility and lifetime. Plastic use reduction, reuse and recycling are thus increasingly important, but making informed policy and research decisions within this space can be challenging given the diverse range of available solutions. The novel analysis framework, Plastic Parallel Pathways Platform (4P) has been equipped with consequential life cycle assessment, techno-economic analysis, and a plastic circularity indicator to estimate the greenhouse gas (GHG) emissions, circularity, and cost of polyethylene terephthalate (PET) down-cycling to lower-quality resin, closed-loop recycling to food-grade PET bottles, up-cycling to fiber-reinforced plastic (FRP), and conversion to non-plastic products (electricity, oil) on a United States economy-wide basis. Integrating system dynamics into this robust plastics model that already incorporates techno-economics, circularity, and environmental impacts will enable identification of key bottlenecks between manufacturers, waste sorters, and reclaimers that currently prevent rapid decarbonization of the plastics economy. System dynamics (SD) explore the evolution of activities and technologies based on changed macro parameters such as plastic demand and supply, market shifts, and cross-sectoral interactions. This project particularly aims to explore the interplay between, waste collection, plastic waste sorting, recycling, and manufacturing, as well as the effect of plastic bale quality and plastic reuse initiatives on the surrounding process stages. This functionality will facilitate combinatory analysis in which a portfolio of end-of-life pathways are assessed simultaneously, with the exact makeup of that portfolio affected by parameters such as technology scales, resource constraints, and waste mitigation efforts. Integrating SD with the 4P framework enables analyzing the effect of increased revenue and reinvestment into improving process efficiencies, sorting and collection quantities. Through that, market effects of increased recycled resin availability can be studied for the plastics systems model for the US. The results will help identify technical or economic bottlenecks that currently limit efforts to decarbonize the U.S. plastics economy.

circular economy↗

DOE OSTI2024

Towards a Circular Economy for PET Bottles in the U.S. - 4P Model

The United States generates the most plastic waste of any country. Along with that GHG emissions from the global plastic economy are expected to increase to 15% of the global carbon budget by 2050. It is imperative that plastic recycling is made a reality to reduce both plastic pollution in the environment and GHG emissions. A portfolio of end-of-life strategies must be implemented to minimize environmental impacts and retain valuable plastic material, but it is challenging to compare options that generate products with different utility and lifetime. Plastic use reduction, reuse and recycling are thus increasingly important, but making informed policy and research decisions within this space can be challenging given the diverse range of available solutions. The novel analysis framework, Plastic Parallel Pathways Platform (4P) has been equipped with consequential life cycle assessment, techno-economic analysis, and a plastic circularity indicator to estimate the greenhouse gas (GHG) emissions, circularity, and cost of polyethylene terephthalate (PET) down-cycling to lower-quality resin, closed-loop recycling to food-grade PET bottles, up-cycling to fiber-reinforced plastic (FRP), and conversion to non-plastic products (electricity, oil) on a United States economy-wide basis. Integrating system dynamics into this robust plastics model that already incorporates techno-economics, circularity, and environmental impacts will enable identification of key bottlenecks between manufacturers, waste sorters, and reclaimers that currently prevent rapid decarbonization of the plastics economy. System dynamics (SD) explore the evolution of activities and technologies based on changed macro parameters such as plastic demand and supply, market shifts, and cross-sectoral interactions. This project particularly aims to explore the interplay between, waste collection, plastic waste sorting, recycling, and manufacturing, as well as the effect of plastic bale quality and plastic reuse initiatives on the surrounding process stages. This functionality will facilitate combinatory analysis in which a portfolio of end-of-life pathways are assessed simultaneously, with the exact makeup of that portfolio affected by parameters such as technology scales, resource constraints, and waste mitigation efforts. Integrating SD with the 4P framework enables analyzing the effect of increased revenue and reinvestment into improving process efficiencies, sorting and collection quantities. Through that, market effects of increased recycled resin availability can be studied for the plastics systems model for the US. The results will help identify technical or economic bottlenecks that currently limit efforts to decarbonize the U.S. plastics economy.

carbon↗

DOE OSTI2024

Towards a Circular Economy for PET Bottles in the U.S. - 4P Framework

The United States generates the most plastic waste of any country. Along with that GHG emissions from the global plastic economy are expected to increase to 15% of the global carbon budget by 2050. It is imperative that plastic recycling is made a reality to reduce both plastic pollution in the environment and GHG emissions. A portfolio of end-of-life strategies must be implemented to minimize environmental impacts and retain valuable plastic material, but it is challenging to compare options that generate products with different utility and lifetime. Plastic use reduction, reuse and recycling are thus increasingly important, but making informed policy and research decisions within this space can be challenging given the diverse range of available solutions. The novel analysis framework, Plastic Parallel Pathways Platform (4P) has been equipped with consequential life cycle assessment, techno-economic analysis, and a plastic circularity indicator to estimate the greenhouse gas (GHG) emissions, circularity, and cost of polyethylene terephthalate (PET) down-cycling to lower-quality resin, closed-loop recycling to food-grade PET bottles, up-cycling to fiber-reinforced plastic (FRP), and conversion to non-plastic products (electricity, oil) on a United States economy-wide basis. Integrating system dynamics into this robust plastics model that already incorporates techno-economics, circularity, and environmental impacts will enable identification of key bottlenecks between manufacturers, waste sorters, and reclaimers that currently prevent rapid decarbonization of the plastics economy. System dynamics (SD) explore the evolution of activities and technologies based on changed macro parameters such as plastic demand and supply, market shifts, and cross-sectoral interactions. This project particularly aims to explore the interplay between, waste collection, plastic waste sorting, recycling, and manufacturing, as well as the effect of plastic bale quality and plastic reuse initiatives on the surrounding process stages. This functionality will facilitate combinatory analysis in which a portfolio of end-of-life pathways are assessed simultaneously, with the exact makeup of that portfolio affected by parameters such as technology scales, resource constraints, and waste mitigation efforts. Integrating SD with the 4P framework enables analyzing the effect of increased revenue and reinvestment into improving process efficiencies, sorting and collection quantities. Through that, market effects of increased recycled resin availability can be studied for the plastics systems model for the US. The results will help identify technical or economic bottlenecks that currently limit efforts to decarbonize the U.S. plastics economy.

circular economy↗

DOE OSTI2025

Extending TOUGH + HYDRATE with a parallel particle transport simulator: numerical investigation of sand production during gas production from hydrate deposits

A new parallel code for simulating particle transport in porous media is integrated with the TOUGH + HYDRATE simulator to investigate sand production associated with gas production from unconsolidated gas hydrate-bearing sediments (HBS). Here, the parallel coupled simulator is named THMPT and uses the integral finite difference method to describe the Darcian and non-Darcian flow of fluids and heat transport, the finite element method to describe the associated geomechanical changes, and the discrete element method to track the trajectory of individual sand particles within the HBS. The THMPT simulator is written in Fortran, incorporates multiple optimized algorithms, and can comprehensively address the coupled flow, thermal, chemical, geomechanical, and particle transport processes that characterize the system behaviors during gas production from HBS. The simulator can capture all processes involved in sand particle transport in porous media, including sand detachment, collision, clogging (i.e., bridging), and migration. A benchmark case study of sand production in the course of depressurization-induced gas production from a representative HBS reveals various distinct microscopic particle migration mechanisms and the adverse impact of sand particle detachment, transport, and clogging. The numerical investigation also examines the effect of bottomhole pressure on mitigating sand production. The simulation results indicate that sand clogging near the wellbore significantly reduces permeability, decreasing gas production by at least 50%. Lastly, the efficiency of gravel packing in mitigating sand production is numerically evaluated, revealing that the structure of the porous media appears to profoundly influence the macroscopic motion behavior of sand particles and sand clogging characteristics.

discrete element method↗

DOE OSTI2021

Encoder–decoder neural network for solving the nonlinear Fokker–Planck–Landau collision operator in XGC

An encoder–decoder neural network has been used to examine the possibility for acceleration of a partial integro-differential equation, the Fokker–Planck–Landau collision operator. This is part of the governing equation in the massively parallel particle-in-cell code XGC, which is used to study turbulence in fusion energy devices. The neural network emphasizes physics-inspired learning, where it is taught to respect physical conservation constraints of the collision operator by including them in the training loss, along with the ℓ 2 loss. In particular, network architectures used for the computer vision task of semantic segmentation have been used for training. A penalization method is used to enforce the ‘soft’ constraints of the system and integrate error in the conservation properties into the loss function. During training, quantities representing the particle density, momentum and energy for all species of the system are calculated at each configuration vertex, mirroring the procedure in XGC. This simple training has produced a median relative loss, across configuration space, of the order of 10 –4 , which is low enough if the error is of random nature, but not if it is of drift nature in time steps. The run time for the current Picard iterative solver of the operator is O(n 2 ), where n is the number of plasma species. As the XGC1 code begins to attack problems including a larger number of species, the collision operator will become expensive computationally, making the neural network solver even more important, especially since its training only scales as O(n). Here, a wide enough range of collisionality has been considered in the training data to ensure the full domain of collision physics is captured. An advanced technique to decrease the losses further will be subject of a subsequent report. Eventual work will include expansion of the network to include multiple plasma species.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

A supernodal all-pairs shortest path algorithm

We show how to exploit graph sparsity in the Floyd-Warshall algorithm for the all-pairs shortest path (Apsp) problem. Floyd-Warshall is an attractive choice for Apsp on high-performing systems due to its structural similarity to solving dense linear systems and matrix multiplication. However, if sparsity of the input graph is not properly exploited, Floyd-Warshall will perform unnecessary asymptotic work and thus may not be a suitable choice for many input graphs. To overcome this limitation, the key idea in our approach is to use the known algebraic relationship between Floyd-Warshall and Gaussian elimination, and import several algorithmic techniques from sparse Cholesky factorization, namely, fill-in reducing ordering, symbolic analysis, supernodal traversal, and elimination tree parallelism. When combined, these techniques reduce computation, improve locality and enhance parallelism. We implement these ideas in an efficient shared memory parallel prototype that is orders of magnitude faster than an efficient multi-threaded baseline Floyd-Warshall that does not exploit sparsity. Our experiments suggest that the Floyd-Warshall algorithm can compete with Dijkstra's algorithm (the algorithmic core of Johnson's algorithm) for several classes sparse graphs.

Sao, Piyush↗

DOE OSTI2021

Dakota A Multilevel Parallel Object-Oriented Framework for Design Optimization Parameter Estimation Uncertainty Quantification and Sensitivity Analysis (V.6.14) (Theory Manual)

The Dakota toolkit provides a flexible and extensible interface between simulation codes and iterative analysis methods. Dakota contains algorithms for optimization with gradient and nongradient-based methods; uncertainty quantification with sampling, reliability, and stochastic expansion methods; parameter estimation with nonlinear least squares methods; and sensitivity/variance analysis with design of experiments and parameter study methods. These capabilities may be used on their own or as components within advanced strategies such as surrogate-based optimization, mixed integer nonlinear programming, or optimization under uncertainty. By employing object-oriented design to implement abstractions of the key components required for iterative systems analyses, the Dakota toolkit provides a flexible and extensible problem-solving environment for design and performance analysis of computational models on high performance computers. This report serves as a theoretical manual for selected algorithms implemented within the Dakota software. It is not intended as a comprehensive theoretical treatment, since a number of existing texts cover general optimization theory, statistical analysis, and other introductory topics. Rather, this manual is intended to summarize a set of Dakota-related research publications in the areas of surrogate-based optimization, uncertainty quantification, and optimization under uncertainty that provide the foundation for many of Dakota's iterative analysis capabilities.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Dakota, A Multilevel Parallel Object-Oriented Framework for Design Optimization, Parameter Estimation, Uncertainty Quantification, and Sensitivity Analysis: Theory Manual (V.6.15)

The Dakota toolkit provides a flexible and extensible interface between simulation codes and iterative analysis methods. Dakota contains algorithms for optimization with gradient and nongradient-based methods; uncertainty quantification with sampling, reliability, and stochastic expansion methods; parameter estimation with nonlinear least squares methods; and sensitivity/variance analysis with design of experiments and parameter study methods. These capabilities may be used on their own or as components within advanced strategies such as surrogate-based optimization, mixed integer nonlinear programming, or optimization under uncertainty. By employing object-oriented design to implement abstractions of the key components required for iterative systems analyses, the Dakota toolkit provides a flexible and extensible problem-solving environment for design and performance analysis of computational models on high performance computers. This report serves as a theoretical manual for selected algorithms implemented within the Dakota software. It is not intended as a comprehensive theoretical treatment, since a number of existing texts cover general optimization theory, statistical analysis, and other introductory topics. Rather, this manual is intended to summarize a set of Dakota-related research publications in the areas of surrogate-based optimization, uncertainty quantification, and optimization under uncertainty that provide the foundation for many of Dakota's iterative analysis capabilities.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Simultaneous imaging and diffraction in the dynamic diamond anvil cell

The ability to visualize a sample undergoing a pressure-induced phase transition allows for the determination of kinetic parameters, such as the nucleation and growth rates of the high-pressure phase. For samples that are opaque to visible light (such as metallic systems), it is necessary to rely on x-ray imaging methods for sample visualization. Here, we present an experimental platform developed at beamline P02.2 at the PETRA III synchrotron radiation source, which is capable of performing simultaneous x-ray imaging and diffraction of samples that are dynamically compressed in piezo-driven diamond anvil cells. This setup utilizes a partially coherent monochromatic x-ray beam to perform lensless phase contrast imaging, which can be carried out using either a parallel- or focused-beam configuration. The capabilities of this platform are illustrated by experiments on dynamically compressed Ga and Ar. Melting and solidification were identified based on the observation of solid/liquid phase boundaries in the x-ray images and corresponding changes in the x-ray diffraction patterns collected during the transition, with significant edge enhancement observed in the x-ray images collected using the focused-beam. These results highlight the suitability of this technique for a variety of purposes, including melt curve determination.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Massively parallel GPU enabled third-order cluster perturbation excitation energies for cost-effective large scale excitation energy calculations

We present here a massively parallel implementation of the recently developed CPS(D-3) excitation energy model that is based on cluster perturbation theory. The new algorithm extends the one developed in Baudin et al. [J. Chem. Phys., 150, 134110 (2019)] to leverage multiple nodes and utilize graphical processing units for the acceleration of heavy tensor contractions. Furthermore, we show that the extended algorithm scales efficiently with increasing amounts of computational resources and that the developed code enables CPS(D-3) excitation energy calculations on large molecular systems with a low time-to-solution. More specifically, calculations on systems with over 100 atoms and 1000 basis functions are possible in a few hours of wall clock time. This establishes CPS(D-3) excitation energies as a computationally efficient alternative to those obtained from the coupled-cluster singles and doubles model.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗