High Performance Computing Architectures Research Articles

Characterizing structural features of two-dimensional particle systems through Voronoi topology

Abstract This paper introduces a new approach toward characterizing local structural features of two-dimensional particle systems. The approach can accurately identify and characterize defects in high-temperature crystals, distinguish a wide range of nominally disordered systems, and robustly describe complex structures such as grain boundaries. This paper also introduces two-dimensional functionality into the open-source software program VoroTop which automates this analysis. This software package is built on a recently-introduced multithreaded version of Voro++, enabling the analysis of systems with billions of particles on high-performance computer architectures.

TOLEDO: enhancing Maestro GUI for non-expert users to perform massive MD simulations

Classical Molecular Dynamics (MD) simulates the dynamical evolution of biological systems at the atomic level. Using MD in conjunction with high-performance computing (HPC) architectures, we can evaluate the possible interactions between a ligand library against one protein target to find a drug that can influence a protein target to cure a disease. Simultaneously, we can also obtain information about their dynamic evolution. One of the primary software packages for MD simulations is Desmond, which employs Maestro for the setup, execution, and analysis of MD through a graphical user interface (GUI), which is suitable even for non-expert users. However, using the GUI, users can typically run only one short (less than 1000 ns) MD each time. Our work aims to create a method/protocol to run several MD simulations simultaneously on a remote HPC cluster within Maestro-Desmond. In this work, we provide TOLEDO (Throughput Optimization of Ligand-Protein Systems Exploration through Dynamics simulation in Optimized HPC systems) to overcome such limitations and run several MD simulations simultaneously. The best feature of TOLEDO is its independence from the usual time constraints of many clusters, with storage space being the only limitation. To run TOLEDO, we prepare/set up the protein-ligand complex before running MD via Maestro GUI. Next, we run the main TOLEDO script for several MD simulations on a supercomputer. When TOLEDO finishes, users obtain reports and graphics. The obtained results are easily interpretable. In essence, TOLEDO significantly enhances MD throughput beyond the capabilities of the Maestro GUI.

Evaluating ARM and RISC-V Architectures for High-Performance Computing with Docker and Kubernetes

This paper thoroughly assesses the ARM and RISC-V architectures in the context of high-performance computing (HPC). It includes an analysis of Docker and Kubernetes integration. Our study aims to evaluate and compare these systems’ performance, scalability, and practicality in a general context and then assess the impact they might have on special use cases, like HPC. ARM-based systems exhibited better performance and seamless integration with Docker and Kubernetes, underscoring their advanced development and effectiveness in managing high-performance computing workloads. On the other hand, despite their open-source architecture, RISC-V platforms presented considerable intricacy and difficulties in working with Kubernetes, which hurt their overall effectiveness and ease of management. The results of our study offer valuable insights into the practical consequences of implementing these architectures for HPC, highlighting ARM’s preparedness and the potential of RISC-V while acknowledging the increased complexity and significant trade-offs involved at this point.

Analysis model of energy consumption variables for data processing in high-performance computing systems

One of the main challenges in the efficient operation of a high-performance computing (HPC) center is the energy consumption generated by the operation of the data center where the HPC equipment is housed, mainly because this consumption is reflected in very high accounts payable, and this may affect the level of service offered to users. The study of the different factors and elements that can make energy consumption more efficient in these data centers provides an opportunity to focus these resources on elements that favor the use of HPC. The design variables provided by manufacturers to manage HPC systems and monitoring systems provide an accurate view of the behavior of these variables according to how they are used. HPC architectures are configured in a very particular way for each HPC data center, creating particular scenarios of operation and performance in each implementation. Various proposals and technologies have been developed for the analysis of the energy consumption of a data center, and the processing elements include a series of indicators and technologies that manufacturers have developed to determine the energy efficiency. This article seeks to identify this series of processing and performance variables, which affect the energy consumption of HPC equipment, for the implemented computing architectures based on the analysis of performance models to obtain a general over-view of their effect on energy consumption in a case study to identify the behaviors of particular job assignment factors and provide an analysis of the energy consumption under particular conditions.

Implementation and evaluation of diabatic advection in the Lagrangian transport model MPTRAC 2.6

Abstract. Diabatic transport schemes with hybrid zeta coordinates, which follow isentropes in the stratosphere, are known to greatly improve Lagrangian transport calculations compared to the kinematic approach. However, some Lagrangian transport calculations with a diabatic approach, such as the Chemical Lagrangian Transport Model of the Stratosphere (CLaMS), are not well prepared to run on modern high-performance computing (HPC) architectures. Here, we implemented and evaluated a new diabatic transport scheme in the Massive-Parallel Trajectory Calculations (MPTRAC) model. While MPTRAC can be used either with shared-memory multiprocessing on CPUs or with GPUs to offload computationally intensive calculations, making it flexible for many HPC applications, it has been limited to kinematic trajectories in pressure coordinates. The extended modelling approach now enables the use of either kinematic or diabatic vertical velocities and the coupling of different MPTRAC modules based on pressure or hybrid zeta coordinates. This study focus on the accuracy of the implementation in comparison to the CLaMS model. The evaluation of the new transport scheme in MPTRAC shows that, after 90 d of forward calculations, distributions of air parcels in the upper troposphere and lower stratosphere (UTLS) are almost identical for MPTRAC and CLaMS. No significant bias between the two Lagrangian models was found. Furthermore, after 1 d, internal uncertainties (e.g. due to interpolation or the numerical integration method) in the Lagrangian transport calculations are at least 1 order of magnitude smaller than external uncertainties (e.g. from reanalysis selection or downsampling of ERA5). Differences between trajectories using either CLaMS or MPTRAC are on the order of the combined internal uncertainties within MPTRAC. Since the largest systematic differences are caused by the reanalysis and the vertical velocity (diabatic vs. kinematic), the results support the development efforts for trajectory codes that can access the full resolution of ERA5 in combination with diabatic vertical velocities. This work is part of a larger effort to adapt Lagrangian transport in state-of-the-art models such as CLaMS and MPTRAC to current and future HPC architectures and exascale applications.

Edge HPC Architectures for AI-Based Video Surveillance Applications

The introduction of artificial intelligence (AI) in video surveillance systems has significantly transformed security practices, allowing for autonomous monitoring and real-time detection of threats. However, the effectiveness and efficiency of AI-powered surveillance rely heavily on the hardware infrastructure, specifically high-performance computing (HPC) architectures. This article examines the impact of different platforms for HPC edge servers, including x86 and ARM CPU-based systems and Graphics Processing Units (GPUs), on the speed and accuracy of video processing tasks. By using advanced deep learning frameworks, a video surveillance system based on YOLO object detection and DeepSort tracking algorithms is developed and evaluated. This study thoroughly assesses the strengths, limitations, and suitability of different hardware architectures for various AI-based surveillance scenarios.

TTCF4LAMMPS: A toolkit for simulation of the non-equilibrium behaviour of molecular fluids at experimentally accessible shear rates

We present TTCF4LAMMPS, a toolkit for performing non-equilibrium molecular dynamics (NEMD) simulations to study the fluid behaviour at low shear rates using the LAMMPS software. By combining direct NEMD simulations and the transient-time correlation function (TTCF) technique, we study the fluid response to shear rates spanning 15 orders of magnitude. We present two examples for simple monatomic systems: one consisting of a bulk liquid and another with a liquid layer confined between two solid walls. The small bulk system is suitable for testing on personal computers, while the larger confined system requires high-performance computing (HPC) resources. We demonstrate that the TTCF formalism can successfully detect the system response for arbitrarily weak external fields. We provide a brief mathematical explanation for this feature. Although we showcase the method for simple monatomic systems, TTCF can be readily extended to study more complex molecular fluids. Moreover, in addition to shear flows, the method can be extended to investigate elongational or mixed flows as well as thermal or electric fields. The high computational cost needed for the method is offset by the two following benefits: i) the cost is independent of the magnitude of the external field, and ii) the simulations can be made highly efficient on HPC architectures by exploiting the parallel design of the algorithm. We expect the toolkit to be useful for computational researchers striving to study the nonequilibrium behaviour of fluids under experimentally-accessible conditions. Program summaryProgram title: TTCF4LAMMPSCPC Library link to program files:https://doi.org/10.17632/hh2rkcxbrf.1Developer's repository link:https://github.com/edwardsmith999/TTCF4LAMMPSLicensing provisions: GNU General Public License 3Programming language: Python 3Nature of problem: Measuring the nonequilibrium behaviour of bulk and confined fluids under experimentally accessible strain rates in non-equilibrium molecular dynamics (NEMD) simulations.Solution method: Creating a Python-based code that utilises the transient-time correlation function method and the LAMMPS software to enable the bulk fluid properties (e.g. viscosity) and confined fluid interfacial properties (e.g. shear stress and slip velocity) to be computed at low shear rates with NEMD.

A parallel and performance portable implementation of a full-field crystal plasticity model

We have developed a parallel implementation of an Elasto-Viscoplastic Fast Fourier Transform-based (EVPFFT) micromechanical solver to enable computationally efficient crystal plasticity modeling for polycrystalline materials. Our primary focus lies in achieving performance portability, allowing a single EVPFFT implementation to run optimally on various homogeneous architectures, including multi-core Central Processing Units (CPUs), as well as on heterogeneous computer architectures comprising multi-core CPUs and Graphics Processing Units (GPUs) from different vendors. To accomplish this goal, we have leveraged MATAR, a C++ software library that simplifies the creation and utilization of multidimensional dense or sparse matrix and array data structures. These data structures are designed to be portable across diverse architectures through the use of Kokkos, a performance-portable library. Additionally, we have employed the Message Passing Interface (MPI) to efficiently distribute the computational workload among processors. The heFFTe (Highly Efficient FFT for Exascale) library is used to facilitate the performance portability of the fast Fourier transforms (FFTs) computation. The computational performance of EVPFFT is evaluated and presented in terms of parallel scalability and simulation runtime on different high-performance computing (HPC) architectures. The utility of the developed framework to efficiently simulate the micro-mechanical fields in polycrystalline microstructures in engineering applications is discussed. Program summaryProgram Title: EVPFFTCPC Library link to program files:https://doi.org/10.17632/2k8579fyyv.1Developer's repository link:https://github.com/lanl/FierroLicensing provisions: BSD 3-Clause LicenseProgramming language: C++External routines/libraries: MPI, Kokkos, MATAR, HeFFTe, HDF5Nature of problem: EVPFFT is a crystal plasticity code designed to compute micro-mechanical fields within a polycrystalline representative volume element (RVE) and predict the macroscale response of the RVE.Solution method: EVPFFT uses the periodic Green's function method in Fourier space to solve the field equations of static stress equilibrium in a periodic spatial domain.

A universal parallel simulation framework for energy pipeline networks on high-performance computers

Energy distribution networks represent crucial infrastructures for modern society, and various simulation tools have been widely used by energy suppliers to manage these intricate networks. However, simulation calculations include a large number of fluid control equations, and computational overhead limits the performance of simulation software. This paper proposes a universal parallel simulation framework for energy pipeline networks that takes advantages of data parallelism and computational independence between network elements. A non-pipe model of an energy supply network is optimized, and the input and output of the network model in the proposed framework are modified, which can reduce the development burden during the numerical computations of the pipeline network and weaken the computational correlation between different simulated components. In addition, independent computations can be performed concurrently through periodic data exchange procedures between component instances, improving the parallelism and efficiency of simulation computations. Further, a parallel water pipelines network simulation computing paradigm based on a heterogeneous computer hardware architecture is used to evaluate the proposed framework’s performance. A series of tests are conducted to verify the accuracy of the proposed framework, and simulation errors of less than 5% are achieved. The results of multi-threaded simulation experiments have demonstrated the feasibility of the proposed framework in a parallel computing approach. Moreover, an Advanced Micro Devices (AMD) Deep Computing Unit (DCU)-parallel program is implemented into a water supply network simulation system; the computational efficiency of this system is compared with that of its serial counterpart. The experimental results show that the proposed framework is appropriate for high-performance computer architectures, and the 18x speed-up ratio demonstrates that the parallel program based on the proposed universal framework outperforms the serial program. That provides the basis for the application of pipe network simulation on high-performance computers.

Vast Parameter Space Exploration of the Virtual Brain: A Modular Framework for Accelerating the Multi-Scale Simulation of Human Brain Dynamics

Global neural dynamics emerge from multi-scale brain structures, with nodes dynamically communicating to form transient ensembles that may represent neural information. Neural activity can be measured empirically at scales spanning proteins and subcellular domains to neuronal assemblies or whole-brain networks connected through tracts, but it has remained challenging to bridge knowledge between empirically tractable scales. Multi-scale models of brain function have begun to directly link the emergence of global brain dynamics in conscious and unconscious brain states with microscopic changes at the level of cells. In particular, adaptive exponential integrate-and-fire (AdEx) mean-field models representing statistical properties of local populations of neurons have been connected following human tractography data to represent multi-scale neural phenomena in simulations using The Virtual Brain (TVB). While mean-field models can be run on personal computers for short simulations, or in parallel on high-performance computing (HPC) architectures for longer simulations and parameter scans, the computational burden remains red heavy and vast areas of the parameter space remain unexplored. In this work, we report that our HPC framework, a modular set of methods used here to implement the TVB-AdEx model for the graphics processing unit (GPU) and analyze emergent dynamics, notably accelerates simulations and substantially reduces computational resource requirements. The framework preserves the stability and robustness of the TVB-AdEx model, thus facilitating a finer-resolution exploration of vast parameter spaces as well as longer simulations that were previously near impossible to perform. Comparing our GPU implementations of the TVB-AdEx framework with previous implementations using central processing units (CPUs), we first show correspondence of the resulting simulated time-series data from GPU and CPU instantiations. Next, the similarity of parameter combinations, giving rise to patterns of functional connectivity, between brain regions is demonstrated. By varying global coupling together with spike-frequency adaptation, we next replicate previous results indicating inter-dependence of these parameters in inducing transitions between dynamics associated with conscious and unconscious brain states. Upon further exploring parameter space, we report a nonlinear interplay between the spike-frequency adaptation and subthreshold adaptation, as well as previously unappreciated interactions between the global coupling, adaptation, and propagation velocity of action potentials along the human connectome. Given that simulation and analysis toolkits are made public as open-source packages, this framework serves as a template onto which other models can be easily scripted. Further, personalized data-sets can be used for for the creation of red virtual brain twins toward facilitating more precise approaches to the study of epilepsy, sleep, anesthesia, and disorders of consciousness. These results thus represent potentially impactful, publicly available methods for simulating and analyzing human brain states.

Scalable Multi-FPGA HPC Architecture for Associative Memory System.

Associative memory is a cornerstone of cognitive intelligence within the human brain. The Bayesian confidence propagation neural network (BCPNN), a cortex-inspired model with high biological plausibility, has proven effective in emulating high-level cognitive functions like associative memory. However, the current approach using GPUs to simulate BCPNN-based associative memory tasks encounters challenges in latency and power efficiency as the model size scales. This work proposes a scalable multi-FPGA high performance computing (HPC) architecture designed for the associative memory system. The architecture integrates a set of hypercolumn unit (HCU) computing cores for intra-board online learning and inference, along with a spike-based synchronization scheme for inter-board communication among multiple FPGAs. Several design strategies, including population-based model mapping, packet-based spike synchronization, and cluster-based timing optimization, are presented to facilitate the multi-FPGA implementation. The architecture is implemented and validated on two Xilinx Alveo U50 FPGA cards, achieving a maximum model size of 200×10 and a peak working frequency of 220 MHz for the associative memory system. Both the memory-bounded spatial scalability and compute-bounded temporal scalability of the architecture are evaluated and optimized, achieving a maximum scale-latency ratio (SLR) of 268.82 for the two-FPGA implementation. Compared to a two-GPU counterpart, the two-FPGA approach demonstrates a maximum latency reduction of 51.72× and a power reduction exceeding 5.28× under the same network configuration. Compared with the state-of-the-art works, the two-FPGA implementation exhibits a high pattern storage capacity for the associative memory task.

On the advantages of radial variation in assembly design for a nuclear reactor core

The geometric flexibility of additively manufactured metals and ceramics generates a vast, open design space requiring advanced modeling and simulation tools for physics simulations and the rigorous definition of design problems. This effort deploys artificial intelligence (AI) and machine learning (ML) algorithms to inform the design space, to enhance evaluation of potential designs, and to generate optimized results more efficiently.This report documents efforts of the Transformational Challenge Reactor (TCR) program to leverage advanced modeling and simulation techniques driven by AI/ML algorithms on high-performance computing (HPC) systems to yield more optimized TCR core designs. A multiphysics ML surrogate model was developed to run on the HPC architectures. The surrogate model is trained on high-fidelity simulation data of coupled neutronics and thermofluidics and is used to quickly evaluate thousands of candidate core designs in parallel, thus driving the evolution of the cooling channel shapes to minimize temperature peaking and material stress. The flexibility of additive manufacturing allows for the design of unique assemblies for each radial ring of the reactor core, which drives an increase in the selected performance metrics.One of the unique contributions made by this work is the capability to computationally demonstrate a 26.5% improvement in the selected reactor core performance metric by exploiting the freedom of axially variable assembly design for a nuclear reactor core. At the same pressure, the drop in the core decreased by a factor of 2.6, and the maximum temperature dropped by 10.7%. Additionally, the radially optimized design uses less materials, with a 15% decrease in total fuel and a 17.5% decrease in total moderator volume. These are the first quantitative results to allow for variation in assembly design in the core’s different radial rings.

An efficient heterogeneous parallel algorithm of the 3D MOC for multizone heterogeneous systems

Recent innovations in powerful high-performance computing (HPC) architecture have enabled high-fidelity whole-core neutron transport simulations in a reasonable amount of time. Among high-fidelity numerical methods, the one-step 3D Method of Characteristic (MOC) is becoming increasingly popular in the reactor design community due to its numerical accuracy, geometrical flexibility and high degrees of parallelism. In this work, we propose a sophisticated heterogeneous parallel algorithm of the 3D MOC based on a supercomputer prototype consisting of brand-new multizone heterogeneous processors called the MT-3000, which has a peak double precision performance of 11.6 TFLOPS when operating at 1.2 GHz. In the 3D MOC calculation procedure, the total run time is dominated by the so-called MOC transport sweep where large amounts of tracks traverse the mesh grids along different directions. We significantly improve the performance of the MOC transport sweep through a three-level parallel algorithm with multiple hiding latency techniques. We have also implemented heterogeneous level-1 basic linear algebra subprograms (BLAS) that's the second computational hot spot. The 3D C5G7 benchmark problem is employed to verify the correctness and efficiency of the heterogeneous parallel 3D MOC algorithm. Numerical results show that it achieves very high performance on the MT-3000-based supercomputer prototype. Specifically, the three-layer parallel algorithm has proven to be highly scalable for each layer, and the Quarter-MT-3000 performance of the MOC transport sweep has achieved 23.6% of the peak performance, which is close to the theoretical upper limit.

Does Academic Research Drive Industrial Innovation in Computer Architecture?—Analyzing Citations to Academic Papers in Patents

Examines the number of citations to papers published in International Conference on Architectural Support for Programming Languages and Operating Systems (ASPLOS), International Symposium on High- Performance Computer Architecture (HPCA), International Symposium on Computer Architecture (ISCA), and International Symposium on Microarchitecture (MICRO). Given that these four conferences are the most prestigious in computer architecture, they are an excellent proxy for academic research as they likely contain the most innovative ideas. As such, papers from these conferences could provide important building blocks for a company’s patented inventions; if so, these papers should be cited in the company’s patents.

Data Locality in High Performance Computing, Big Data, and Converged Systems: An Analysis of the Cutting Edge and a Future System Architecture

Big data has revolutionized science and technology leading to the transformation of our societies. High-performance computing (HPC) provides the necessary computational power for big data analysis using artificial intelligence and methods. Traditionally, HPC and big data had focused on different problem domains and had grown into two different ecosystems. Efforts have been underway for the last few years on bringing the best of both paradigms into HPC and big converged architectures. Designing HPC and big data converged systems is a hard task requiring careful placement of data, analytics, and other computational tasks such that the desired performance is achieved with the least amount of resources. Energy efficiency has become the biggest hurdle in the realization of HPC, big data, and converged systems capable of delivering exascale and beyond performance. Data locality is a key parameter of HPDA system design as moving even a byte costs heavily both in time and energy with an increase in the size of the system. Performance in terms of time and energy are the most important factors for users, particularly energy, due to it being the major hurdle in high-performance system design and the increasing focus on green energy systems due to environmental sustainability. Data locality is a broad term that encapsulates different aspects including bringing computations to data, minimizing data movement by efficient exploitation of cache hierarchies, reducing intra- and inter-node communications, locality-aware process and thread mapping, and in situ and transit data analysis. This paper provides an extensive review of cutting-edge research on data locality in HPC, big data, and converged systems. We review the literature on data locality in HPC, big data, and converged environments and discuss challenges, opportunities, and future directions. Subsequently, using the knowledge gained from this extensive review, we propose a system architecture for future HPC and big data converged systems. To the best of our knowledge, there is no such review on data locality in converged HPC and big data systems.

P3 problem and Magnolia language: Specializing array computations for emerging architectures

The problem of producing portable high-performance computing (HPC) software that is cheap to develop and maintain is called the P3 (performance, portability, productivity) problem. Good solutions to the P3 problem have been achieved when the performance profiles of the target machines have been similar. The variety of HPC architectures is, however, large and can be expected to grow larger. Software for HPC therefore needs to be highly adaptable, and there is a pressing need to provide developers with tools to produce software that can target machines with vastly different profiles. Multi-dimensional array manipulation constitutes a core component of numerous numerical methods, such as finite difference solvers of Partial Differential Equations (PDEs). The efficiency of these computations is tightly connected to traversing and distributing array data in a hardware-friendly way. The Mathematics of Arrays (MoA) allows for formally reasoning about array computations and enables systematic transformations of array-based programs, e.g., to use data layouts that fit to a specific architecture. This paper presents a programming methodology aimed for tackling the P3 problem in domains that are well-explored using Magnolia, a language designed to embody generic programming. The Magnolia programmer can restrict the semantic properties of abstract generic types and operations by defining so-called axioms. Axioms can be used to produce tests for concrete implementations of specifications, for formal verification, or to perform semantics-preserving program transformations. We leverage Magnolia's semantic specification facilities to extend the Magnolia compiler with a term rewriting system. We implement MoA's transformation rules in Magnolia, and demonstrate through a case study on a finite difference solver of PDEs how our rewriting system allows exploring the space of possible optimizations.

HPC-enabling technologies for high-fidelity combustion simulations

With the increase in computational power in the last decade and the forthcoming Exascale supercomputers, a new horizon in computational modelling and simulation is envisioned in combustion science. Considering the multiscale and multiphysics characteristics of turbulent reacting flows, combustion simulations are considered as one of the most computationally demanding applications running on cutting-edge supercomputers. Exascale computing opens new frontiers for the simulation of combustion systems as more realistic conditions can be achieved with high-fidelity methods. However, an efficient use of these computing architectures requires methodologies that can exploit all levels of parallelism. The efficient utilization of the next generation of supercomputers needs to be considered from a global perspective, that is, involving physical modelling and numerical methods with methodologies based on High-Performance Computing (HPC) and hardware architectures. This review introduces recent developments in numerical methods for large-eddy simulations (LES) and direct-numerical simulations (DNS) to simulate combustion systems, with focus on the computational performance and algorithmic capabilities. Due to the broad scope, a first section is devoted to describe the fundamentals of turbulent combustion, which is followed by a general description of state-of-the-art computational strategies for solving these problems. These applications require advanced HPC approaches to exploit modern supercomputers, which is addressed in the third section. The increasing complexity of new computing architectures, with tightly coupled CPUs and GPUs, as well as high levels of parallelism, requires new parallel models and algorithms exposing the required level of concurrency. Advances in terms of dynamic load balancing, vectorization, GPU acceleration and mesh adaptation have permitted to achieve highly-efficient combustion simulations with data-driven methods in HPC environments. Therefore, dedicated sections covering the use of high-order methods for reacting flows, integration of detailed chemistry and two-phase flows are addressed. Final remarks and directions of future work are given at the end.

NonInvasive Multigrid For SemiStructured Grids

Multigrid solvers for hierarchical hybrid grids (HHG) have been proposed to promote the efficient utilization of high performance computer architectures. These HHG meshes are constructed by uniformly refining a relatively coarse fully unstructured mesh. While HHG meshes provide some flexibility for unstructured applications, most multigrid calculations can be accomplished using efficient structured grid ideas and kernels. This paper focuses on generalizing the HHG idea so that it is applicable to a broader community of computational scientists, and so that it is easier for existing applications to leverage structured multigrid components. Specifically, we adapt the structured multigrid methodology to significantly more complex semi-structured meshes. Further, we illustrate how mature applications might adopt a semi-structured solver in a relatively non-invasive fashion. To do this, we propose a formal mathematical framework for describing the semi-structured solver. This formalism allows us to precisely define the associated multigrid method and to show its relationship to a more traditional multigrid solver. Additionally, the mathematical framework clarifies the associated software design and implementation. Numerical experiments highlight the relationship of the new solver with classical multigrid. We also demonstrate the generality and potential performance gains associated with this type of semi-structured multigrid.

Turbulence Modulation by Slender Fibers

In this paper, we numerically investigate the turbulence modulation produced by long flexible fibres in channel flow. The simulations are based on an Euler–Lagrangian approach, where fibres are modelled as chains of constrained, sub-Kolmogorov rods. A novel algorithm is deployed to make the resolution of dispersed systems of constraint equations, which represent the fibres, compatible with a state-of-the-art, Graphics Processing Units-accelerated flow-solver for direct numerical simulations in the two-way coupling regime on High Performance Computing architectures. Two-way coupling is accounted for using the Exact Regularized Point Particle method, which allows to calculate the disturbance generated by the fibers on the flow considering progressively refined grids, down to a quasi-viscous length-scale. The bending stiffness of the fibers is also modelled, while collisions are neglected. Results of fluid velocity statistics for friction Reynolds number of the flow Reτ=150 and fibers with Stokes number St = 0.01 (nearly tracers) and 10 (inertial) are presented, with special regard to turbulence modulation and its dependence on fiber inertia and volume fraction (equal to ϕ=2.12·10−5 and 2.12·10−4). The non-Newtonian stresses determined by the carried phase are also displayed, determined by long and slender fibers with fixed aspect ratio λtot=200, which extend up to the inertial range of the turbulent flow.

Responsibly Reckless Matrix Algorithms for HPC Scientific Applications

High-performance computing (HPC) achieved an astonishing three orders of magnitude performance improvement per decade for three decades, thanks to hardware technology scaling resulting in exponential improvement in the rate of floating point executions, though slowing in the most recent. Captured in the Top500 list, this hardware evolution cascaded through the software stack, triggering changes at all levels, including the redesign of numerical linear algebra libraries. HPC simulations on massively parallel systems are often driven by matrix computations, whose rate of execution depends on their floating point precision. Referred to by Jack Dongarra, the 2021 ACM A.M. Turing Award Laureate, as “responsibly reckless” matrix algorithms, we highlight the implications of mixed-precision (MP) computations for HPC applications. Introduced 75 years ago, long before the advent of HPC architectures, MP numerical methods turn out to be paramount for increasing the throughput of traditional and artificial intelligence (AI) workloads beyond riding the wave of the hardware alone. Reducing precision comes at the price of trading away some accuracy for performance (reckless behavior) but in noncritical segments of the workflow (responsible behavior) so that accuracy requirements of the application can still be satisfied. They offer a valuable performance/accuracy knob and, just as they are in AI, they are now indispensable in the pursuit of knowledge and discovery in simulations. In particular, we illustrate the MP impact on three representative HPC applications related to seismic imaging, climate/environment geospatial predictions, and computational astronomy.