Computing Sciences

Over the academic year 2022–23, we discussed the teaching of software performance engineering with more than a dozen faculty across North America and beyond. Our outreach was centered on research-focused faculty with an existing interest in this course material. These discussions revealed an enthusiasm for making software performance engineering a more prominent part of a curriculum for computer scientists and engineers. Here, we discuss how MIT’s longstanding efforts in this area may serve as a launching point for community development of a software performance engineering curriculum, challenges in and solutions for providing the necessary infrastructure to universities, and future directions.

The processing and analysis of computed tomography (CT) imaging is important for both basic scientific development and clinical applications. In AutoCT, we provide a comprehensive pipeline that integrates an end-to-end automatic preprocessing, registration, segmentation, and quantitative analysis of 3D CT scans. The engineered pipeline enables atlas-based CT segmentation and quantification leveraging diffeomorphic transformations through efficient forward and inverse mappings. The extracted localized features from the deformation field allow for downstream statistical learning that may facilitate medical diagnostics. On a lightweight and portable software platform, AutoCT provides a new toolkit for the CT imaging community to underpin the deployment of artificial intelligence-driven applications.

The Gaussian process (GP) is a popular statistical technique for stochastic function approximation and uncertainty quantification from data. GPs have been adopted into the realm of machine learning (ML) in the last two decades because of their superior prediction abilities, especially in data-sparse scenarios, and their inherent ability to provide robust uncertainty estimates. Even so, their performance highly depends on intricate customizations of the core methodology, which often leads to dissatisfaction among practitioners when standard setups and off-the-shelf software tools are being deployed. Arguably, the most important building block of a GP is the kernel function, which assumes the role of a covariance operator. Stationary kernels of the Matérn class are used in the vast majority of applied studies; poor prediction performance and unrealistic uncertainty quantification are often the consequences. Non-stationary kernels show improved performance but are rarely used due to their more complicated functional form and the associated effort and expertise needed to define and tune them optimally. In this perspective, we want to help ML practitioners make sense of some of the most common forms of non-stationarity for Gaussian processes. We show a variety of kernels in action using representative datasets, carefully study their properties, and compare their performances. Based on our findings, we propose a new kernel that combines some of the identified advantages of existing kernels.

Predictive high-fidelity modeling of wind turbines with computational fluid dynamics, wherein turbine geometry is resolved in an atmospheric boundary layer, is important to understanding complex flow accounting for design strategies and operational phenomena such as blade erosion, pitch-control, stall/vortex-induced vibrations, and aftermarket add-ons. The biggest challenge with high-fidelity modeling is the realization of numerical algorithms that can capture the relevant physics in detail through effective use of high-performance computing. For modern supercomputers, that means relying on GPUs for acceleration. In this paper, we present ExaWind, a GPU-enabled open-source incompressible-flow hybrid-computational fluid dynamics framework, comprising the near-body unstructured grid solver Nalu-Wind, and the off-body block-structured-grid solver AMR-Wind, which are coupled using the Topology Independent Overset Grid Assembler. Turbine simulations employ either a pure Reynolds-averaged Navier–Stokes turbulence model or hybrid turbulence modeling wherein Reynolds-averaged Navier–Stokes is used for near-body flow and large eddy simulation is used for off-body flow. Being two-way coupled through overset grids, the two solvers enable simulation of flows across a huge range of length scales, for example, 10 orders of magnitude going from O(μm) boundary layers along the blades to O(10 km) across a wind farm. In this paper, we describe the numerical algorithms for geometry-resolved turbine simulations in atmospheric boundary layers using ExaWind. We present verification studies using canonical flow problems. Validation studies are presented using megawatt-scale turbines established in literature. Additionally presented are demonstration simulations of a small wind farm under atmospheric inflow with different stability states.

A comprehensive understanding of human-induced changes to rainfall is essential for water resource management and infrastructure design. However, at regional scales, existing detection and attribution studies are rarely able to conclusively identify human influence on precipitation. Here we show that anthropogenic aerosol and greenhouse gas (GHG) emissions are the primary drivers of precipitation change over the United States. GHG emissions increase mean and extreme precipitation from rain gauge measurements across all seasons, while the decadal-scale effect of global aerosol emissions decreases precipitation. Local aerosol emissions further offset GHG increases in the winter and spring but enhance rainfall during the summer and fall. Our results show that the conflicting literature on historical precipitation trends can be explained by offsetting aerosol and greenhouse gas signals. At the scale of the United States, individual climate models reproduce observed changes but cannot confidently determine whether a given anthropogenic agent has increased or decreased rainfall.

Global warming increases available sensible and latent heat energy, increasing the thermodynamic potential wind intensity of tropical cyclones (TCs). Supported by theory, observations, and modeling, this causes a shift in mean TC intensity, which tends to manifest most clearly at the greatest intensities. The Saffir-Simpson scale for categorizing damage based on the wind intensity of TCs was introduced in the early 1970s and remains the most commonly used metric for public communication of the level of wind hazard that a TC poses. Because the scale is open-ended and does not extend beyond category 5 (70 m/s windspeed or greater), the level of wind hazard conveyed by the scale remains constant regardless of how far the intensity extends beyond 70 m/s. This may be considered a weakness of the scale, particularly considering that the destructive potential of the wind increases exponentially. Here, we consider how this weakness becomes amplified in a warming world by elucidating the past and future increases of peak wind speeds in the most intense TCs. A simple extrapolation of the Saffir-Simpson scale is used to define a hypothetical category 6, and we describe the frequency of TCs, both past and projected under global warming, that would fall under this category. We find that a number of recent storms have already achieved this hypothetical category 6 intensity and based on multiple independent lines of evidence examining the highest simulated and potential peak wind speeds, more such storms are projected as the climate continues to warm.

The trade-off between robustness and tunability is a central challenge in the pursuit of quantum simulation and fault-tolerant quantum computation. In particular, quantum architectures are often designed to achieve high coherence at the expense of tunability. Many current qubit designs have fixed energy levels and consequently limited types of controllable interactions. Here by adiabatically transforming fixed-frequency superconducting circuits into modifiable Floquet qubits, we demonstrate an XXZ Heisenberg interaction with fully adjustable anisotropy. This interaction model can act as the primitive for an expressive set of quantum operations, but is also the basis for quantum simulations of spin systems. To illustrate the robustness and versatility of our Floquet protocol, we tailor the Heisenberg Hamiltonian and implement two-qubit iSWAP, CZ and SWAP gates with good estimated fidelities. In addition, we implement a Heisenberg interaction between higher energy levels and employ it to construct a three-qubit CCZ gate, also with a competitive fidelity. Our protocol applies to multiple fixed-frequency high-coherence platforms, providing a collection of interactions for high-performance quantum information processing. It also establishes the potential of the Floquet framework as a tool for exploring quantum electrodynamics and optimal control.

We demonstrate an indirect, rather than direct, role of quasi-resonant amplification of planetary waves in a summer weather extreme. We find that there was an interplay between a persistent, amplified large-scale atmospheric circulation state and soil moisture feedbacks as a precursor for the June 2021 Pacific Northwest Heat Dome event. An extended resonant planetary wave configuration prior to the event created an antecedent soil moisture deficit that amplified lower atmospheric warming through strong nonlinear soil moisture feedbacks, favoring this unprecedented heat event.