Ph.D
Group : Parallel Systems
Solving Incompressible Navier-Stokes Equations on Heterogeneous Parallel Architectures
Starts on 03/10/2011
Advisor : BABOULIN, Marc
[LE MAITRE Olivier]
Funding : contrat doctoral DIGITEO
Affiliation : Université Paris-Saclay
Laboratory : LRI
Defended on 09/04/2015, committee :
Directeurs de thèse
M. Marc Baboulin, Professeur, Université Paris-Sud
M. Olivier Le Maître, Directeur de Recherche, LIMSI/CNRS
Rapporteurs
Mme. Fabienne Jézéquel, Maître de Conférences, LIP6
Mme. Masha Sosonkina, Professeur, Old Dominion University, USA
Examinateurs
M. Abdel Lisser, Professeur, Université Paris-Sud
M. Michel Kern, Chargé de Recherche, INRIA
Membre invité
M. Yann Fraigneau, Ingénieur de Recherche, LIMSI/CNRS
Research activities :
Abstract :
In this PhD thesis, we present our research in the domain of high performance software for computational fluid dynamics (CFD). With the increasing demand of high-resolution simulations, there is a need of numerical solvers that can fully take advantage of current manycore accelerated parallel architectures. In this thesis we focus more specifically on developing an efficient parallel solver for 3D incompressible Navier-Stokes (NS) equations on heterogeneous CPU/GPU architectures.
We first present an overview of the CFD domain along with the NS equations for incompressible fluid flows and existing numerical methods. We describe the mathematical model and the numerical method that we chose, based on an incremental prediction-projection method.
A balanced distribution of the computational workload is obtained by using a domain decomposition method. A two-level parallelization combined with SIMD vectorization is used in our implementation to take advantage of the current distributed multicore machines. Numerical experiments on various parallel architectures show that this solver provides satisfying performance and good scalability.
In order to further improve the performance of the NS solver, we integrate GPU computing to accelerate the most time-consuming tasks. The resulting solver can be configured for running on various heterogeneous architectures by specifying explicitly the numbers of MPI processes, threads and GPUs.
This thesis manuscript also includes simulation results for two benchmarks designed from real physical cases. The computed solutions are compared with existing reference results. The code developed in this work will be the base for a future CFD library for parallel CPU/GPU computations.
