🤖 AI Summary
To address the computational challenges of multiphysics flow simulations in engineering, medicine, and fundamental sciences, this work introduces MFC 5.0—a high-performance CFD solver. Methodologically, it unifies support for immersed boundary methods, multiphase flow with phase change, fluid–structure interaction, and reactive flows; features a novel Pyrometheus-driven thermochemical automatic code generation framework; and integrates advanced numerical techniques including relaxed characteristic boundary conditions, Strang splitting, and low-Mach-number formulations. Architecturally, it achieves end-to-end scalability—from single GPU/APU nodes to exascale systems (e.g., Frontier and El Capitan)—leveraging WENO schemes, Euler–Euler/Lagrangian subgrid modeling, and heterogeneous acceleration. Strong and weak scaling efficiencies exceed 90%, with measured performance reaching exaFLOPS-level throughput. These advances significantly enable large-scale simulations critical to aerospace, energy, and biomedical applications.
📝 Abstract
Engineering, medicine, and the fundamental sciences broadly rely on flow simulations, making performant computational fluid dynamics solvers an open source software mainstay. A previous work made MFC 3.0 a published open source source solver with many features. MFC 5.0 is a marked update to MFC 3.0, including a broad set of well-established and novel physical models and numerical methods and the introduction of GPU and APU (or superchip) acceleration. We exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC's single-GPU/APU performance, MFC achieves exascale computation in practice. With these capabilities, MFC has evolved into a tool for conducting simulations that many engineering challenge problems hinge upon. New physical features include the immersed boundary method, $N$-fluid phase change, Euler--Euler and Euler--Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Frontier and LLNL El Capitan see efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-time increase in device count show parallel efficiencies over 90% on OLCF Frontier. Other MFC improvements include ensuring code resilience and correctness with a continuous integration suite, the use of metaprogramming to reduce code length and maintain performance portability, and efficient computational representations for chemical reactions and thermodynamics via code generation with Pyrometheus.