210 likes | 437 Views
Parallel Simulations of Underground Flow in Porous and Fractured Media. H. Mustapha 1,2 , A. Beaudoin 1 , J. Erhel 1 and J.R. De Dreuzy 2. 1 - IRISA – INRIA de Rennes 2 - Géosciences de Rennes. Parallel Simulations of Underground Flow in Porous and Fractured Media. Outline.
E N D
Parallel Simulations of Underground Flow in Porous and Fractured Media H. Mustapha1,2, A. Beaudoin1, J. Erhel1 and J.R. De Dreuzy2 1 - IRISA – INRIA de Rennes 2 - Géosciences de Rennes
Parallel Simulations of Underground Flow in Porous and Fractured Media Outline • Problem presentation • Geometrical and flow fluid model • Numerical method and parallel matrix generation • Parallel linear solver • Direct solver - PSPASES • Iterative solver – HYPRE • Results and performances • Tests and architecture • Complexity and scalability analysis with PSPASES • Comparison between PSPASES and HYPRE • Conclusions and future work H. Mustapha INRIA, Parco sept 2005
Parallel Simulations of Underground Flow in Porous and Fractured Media Outline • Problem presentation • Geometrical and flow fluid model • Numerical method and parallel matrix generation • Parallel linear solver • Direct solver - PSPASES • Iterative solver – HYPRE • Results and performances • Tests and architecture • Complexity and scalability analysis with PSPASES • Comparison between PSPASES and HYPRE • Conclusions and future work H. Mustapha INRIA, Parco sept 2005
Parallel Simulations of Underground Flow in Porous and Fractured Media Geometrical model 2D porous media 3D fracture network Heterogeneous hydraulic conductivity field Discrete fracture network H. Mustapha INRIA, Parco sept 2005
Equations Q = - K*grad (h) div (Q) = 0 Fixed head Nul flux Parallel Simulations of Underground Flow in Porous and Fractured Media Flow fluid model • Boundary conditions 2D porous media 3D fracture network Nul flux Fixed head Fixed head Nul flux H. Mustapha INRIA, Parco sept 2005
Parallel Simulations of Underground Flow in Porous and Fractured Media Outline • Problem presentation • Geometrical and flow fluid model • Numerical method and parallel matrix generation • Parallel linear solver • direct solver – PSPASES • iterative solver – HYPRE • Result and performances • tests and architecture • complexity and scalability analysis with PSPASES • comparison between PSPASES and HYPRE • Conclusions and future work H. Mustapha INRIA, Parco sept 2005
Parallel Simulations of Underground Flow in Porous and Fractured Media Numerical method and parallel matrix generation • Mixed Hybrid Finite Element Method • Type of mesh • Global linear system Conforming triangular mesh Regular grid Assembled by all corresponding sub-domains linear systems H. Mustapha INRIA, Parco sept 2005
Parallel Simulations of Underground Flow in Porous and Fractured Media Outline • Problem presentation • Geometrical and flow fluid model • Numerical method and parallel matrix generation • Parallel linear solver • Direct solver – PSPASES • Iterative solver – HYPRE • Results and performances • Tests and architecture • Complexity and scalability analysis with PSPASES • Comparison between PSPASES and HYPRE • Conclusions and future work H. Mustapha INRIA, Parco sept 2005
Parallel Simulations of Underground Flow in Porous and Fractured Media Parallel Linear solver • Direct solver – PSPASES • sparse matrices • symmetric positive definite matrices • distributed-memory paradigm and MPI library • P (number of processors) = 2^N • Cholesky factorization • Iterative solver – HYPRE • sparse matrices • V-cycle multi-grid algorithm (SMG) • distributed-memory paradigm and MPI library H. Mustapha INRIA, Parco sept 2005
Parallel Simulations of Underground Flow in Porous and Fractured Media Outline • Problem presentation • Geometrical and flow fluid model • Numerical method and parallel matrix generation • Parallel linear solver • Direct solver – PSPASES • Iterative solver - HYPRE • Results and performances • Parallel architecture • Complexity and scalability analysis with PSPASES • Comparison between PSPASES and HYPRE • Conclusions and future work H. Mustapha INRIA, Parco sept 2005
Parallel Simulations of Underground Flow in Porous and Fractured Media Parallel architecture Parallel architecture distributed memory 2 nodes of 32 bi – processors (Proc AMD Opteron 2Ghz with 2Go of RAM) H. Mustapha INRIA, Parco sept 2005
Parallel Simulations of Underground Flow in Porous and Fractured Media Outline • Problem presentation • Geometrical and flow fluid model • Numerical method and parallel matrix generation • Parallel linear solver • Direct solver – PSPASES • Iterative solver - HYPRE • Results and performances • Parallel architecture • Complexity and scalability analysis with PSPASES • Comparison between PSPASES and HYPRE • Conclusions and future work H. Mustapha INRIA, Parco sept 2005
Parallel Simulations of Underground Flow in Porous and Fractured Media Complexity analysis with PSPASES • Test 1 : CPU time of matrix generation, linear solving and flow computation obtained with two processors H. Mustapha INRIA, Parco sept 2005
Parallel Simulations of Underground Flow in Porous and Fractured Media Complexity analysis with PSPASES • Test 2 : Memory requirements for both applications H. Mustapha INRIA, Parco sept 2005
Parallel Simulations of Underground Flow in Porous and Fractured Media Scalability analysis with PSPASES • Test 1 : Efficiency 2D medium 3D fracture network H. Mustapha INRIA, Parco sept 2005
Parallel Simulations of Underground Flow in Porous and Fractured Media Scalability analysis with PSPASES • Test 2 : Speed-up H. Mustapha INRIA, Parco sept 2005
Parallel Simulations of Underground Flow in Porous and Fractured Media Outline • Problem presentation • Geometrical and flow fluid model • Numerical method and parallel matrix generation • Parallel linear solver • Direct solver - PSPASES • Iterative solver - HYPRE • Results and performances • Parallel architecture • Complexity and scalability analysis with PSPASES • Comparison between PSPASES and HYPRE • Conclusions and future work H. Mustapha INRIA, Parco sept 2005
Parallel Simulations of Underground Flow in Porous and Fractured Media Comparison between PSPASES and HYPRE • Test 1 : Linear solving time HYPRE PSPASES H. Mustapha INRIA, Parco sept 2005
Parallel Simulations of Underground Flow in Porous and Fractured Media Comparison between PSPASES and HYPRE • Test 2 : Speed-up PSPASES HYPRE H. Mustapha INRIA, Parco sept 2005
Parallel Simulations of Underground Flow in Porous and Fractured Media Outline • Problem presentation • Geometrical and flow fluid model • Numerical method and parallel matrix generation • Parallel linear solver • Direct solver – PSPASES • Iterative solver – HYPRE • Results and performances • Parallel architecture • Complexity and scalability analysis with PSPASES • Comparison between PSPASES and HYPRE • Conclusions and future work H. Mustapha INRIA, Parco sept 2005
Parallel Simulations of Underground Flow in Porous and Fractured Media • Conclusions Scalable in the case of 2D medium Direct solver PSPASES Not scalable in the case of 3D fracture network Iterative solver HYPRE Efficient for large linear systems • Future work • A 3D extension of parallel software in the case of 2D porous media • Parallel software for simulating solute migration H. Mustapha INRIA, Parco sept 2005