160 likes | 259 Views
High performance Computing Applied to a Saltwater Intrusion Numerical Model. E. Canot IRISA/CNRS J. Erhel IRISA/INRIA Rennes C. de Dieuleveult IRISA/INRIA Rennes. Outline. Context Model and Test Cases Coupled Model Test Cases : Henry and Elder
E N D
High performance Computing Applied to a Saltwater Intrusion Numerical Model E. Canot IRISA/CNRS J. Erhel IRISA/INRIA Rennes C. de Dieuleveult IRISA/INRIA Rennes
Outline • Context • Model and Test Cases • Coupled Model • Test Cases : Henry and Elder • Parallel performances • Presentation • Results • Conclusion
Context Context Effect of pumping in coastal aquifers Necessity to predict the evolution of the water supply.
Coupled Model Presentation of the software • Based on TVDV-2D software developped at IMFS in Strasbourg • Simulation of density driven coupled flow and transport in a porous media • Originally sequential
Coupled Model Model FLOW TRANSPORT convection dispersion
Time n Coupled Model Strong Coupling Velocity FLOW TRANSPORT Density Concentration Time n+1
Test Cases Henry Stable test case
Unstable test case Test Cases Elder
Parallel performances First parallel version (v1) Parallel linear solvers
Parallel performances Parallel sparse solver Use of MUMPS (« MUltifrontal Massively Parallel Solver »), a free package for solving linear systems of equations Ax=b Adapted for sparse unsymmetric, and symmetric definite positive matrices.
Parallel performances Global parallelisation (v2) • Partitioning mesh thanks to METIS, a free package for partitioning graphs, meshes and for producing fill reducing orderings for sparse matrices. Partitioning example with 5 parts • Better data distribution.
Results Sensitivity to convergence criterion Elder (mesh 256x160) Due to the convergence criterion
Results Sensitivity to convergence criterion Elder (mesh 256x160), time step number 8
Results Results for the first version (v1) Henry (mesh 510x254) Wall clock (in seconds) Number of processors
Results Results for the second version (v2) Henry (mesh 510x254) Wall clock (in seconds) Number of processors
Conclusion Conclusion - Perspectives • 3D geometry • Adaptative mesh • Improved coupling