Massively parallel numerical simulation using up to 36,000 CPU cores of an industrial-scale polydispersed reactive pressurized fluidized bed with a mesh of one billion cells

For the last 30 years, experimental and modeling studies have been carried out on fluidized bed reactors from laboratory up to industrial scales. The application of developed models for predictive simulations has however been strongly limited by the available computational power and the capability of computational fluid dynamics software to handle large enough simulations. In recent years, both aspects have made significant advances and we thus now demonstrate the feasibility of a massively parallel simulation, on whole supercomputers using NEPTUNE_CFD, of an industrial-scale polydispersed fluidized-bed reactor. This simulation of an olefin polymerization reactor makes use of an unsteady Eulerianmulti-fluid approach and relies on a billion cellsmeshing. This is a worldwide premiere as the obtained accuracy is yet unmatched for such a large-scale system. The interest of this work is two-fold. In terms of High Performance Computation (HPC), all steps of setting-up the simulation, running it with NEPTUNE_CFD, and post-processing results induce multiple challenges due to the scale of the simulation. The simulation ran using 1260 up to 36,000 cores on supercomputers, used 15 millions of CPU hours and generated 200 TB of rawdata for a simulated physical time of 25s. This article details the methodology applied to handle this simulation, and also focuses on computation performances in terms of profiling, code efficiency and partitioning method suitability. Though being by itself interesting, the HPC challenge is not the only goal of this work as performing this highly-resolved simulation will benefit chemical engineering and CFD communities. Indeed, this computation enables the possibility to account, in a realistic way, for complex flows in an industrial-scale reactor. The predicted behavior is described, and results are post-processed to develop sub-grid models. These will allow for lower-cost simulations with coarser meshes while still encompassing local phenomenon