A parallel mixed–hybrid finite element method for two phase flow problems in porous media using MP

A parallel mixed–hybrid finite element method for two phase flow problems in porous media using MPI

Jakub Solovský , Radek Fučík

Faculty of Nuclear Sciences and Physical Engineering, Czech Technical University in Prague, Trojanova 13, 120 00 Prague, Czech Republic.




This paper deals with a numerical solution of two–phase flow problems in porous media. To solve this type of problems, we propose a numerical method based on mixed–hybrid finite element method. The method is briefly introduced for arbitrary dimension but only a 2D case is considered in this work. We implement several variations of this method using The accuracy and the computational efficiency of the method is verified using a problem with a known exact solution. Numerical experiments show that the errors are similar for all variations of the method. The method is convergent and the experimental order of convergence is slightly less than one. There are differences in the computational time in favor to the iterative solvers, especially when using finer meshes. For computations on fine meshes it is also advantageous to use the parallelism that significantly speeds-up the computation.different approaches to solve the resulting system of linear algebraic equations. Direct and iterative solvers are used and a parallel implementation of this method based on the domain decomposition method using MPI is described.

Cite as:

Solovský , J., Fučík, R. (2017). A parallel mixed–hybrid finite element method for two phase flow problems in porous media using MPI. Computer Methods in Materials Science, 17(2), 84 – 93. https://doi.org/10.7494/cmms.2017.2.0593

Article (PDF):


Mixed–Hybrid Finite Element Method, MPI, Parallel Implementation, Two Phase Flow in Porous Media


Brezzi, F. Fortin, M., 1991, Mixed and Hybrid Finite ElementMethods, Springer-Verlag.Brooks, R., Corey, A., 1964, Hydraulic Properties of PorousMedia, Colorado State University, Hydrology Paper, 3,27.

Burdine, N., 1953, Relative Permeability Calculations FromPore Size Distribution Data, Journal of PetroleumTechnology, 5, 71-78.

Davis, T. A., 2004, Algorithm 832: UMFPACK, anunsymmetric-pattern multifrontal method, ACMTransactions on Mathematical Software, 30, 196–199.

Fučík, R., Illangasekare, T. H., Beneš, M., 2016,Multidimensional self-similar analytical solutions oftwo-phase flow in porous media, Advances in WaterResources, 90, 51-56.

Fučík, R., Mikyška, J., Solovský, J., Klinkovský, J., Oberhuber,T., 2017, Multidimensional Mixed-Hybrid FiniteElement Method for Compositional Two-Phase Flow inHeterogeneous Porous Media and its Massively ParallelImplementation on GPU, in review in Computer PhysicsCommunications.Fučík, R., Mikyška, J., 2011, Mixed-hybrid finite elementmethod for modelling two-phase flow in porous media,Journal of Math for Industry, 3, 9-19.

Fučík, R., Mikyška, J., Beneš, M., Illangasekare, T. H., 2007,An Improved Semi-Analytical Solution for Verificationof Numerical Models of Two Phase Flow in PorousMedia, Vadose Zone Journal, 6(1), 93-104.

Geuzaine, C., Remacle, J. F., 2009, Gmsh: a three-dimensionalfinite element mesh generator with built-in pre- andpost-processing facilities, International Journal forNumerical Methods in Engineering, 79, 1309-1331.

LeVeque, R. J., 2002, Finite Volume Methods for HyperbolicProblems, Cambridge University Press.McWhorter, D. B., Sunada, D. K., 1990, Exact integral solutionsfor two-phase flow, Watter Resources Research, 26,399-413.

Mualem, Y., 1976, A new model for predicting the hydraulicconductivity of unsaturated porous media, WaterResources Research, 12, 513-522.

Saad, Y., 2003, Iterative Methods for Sparse Linear Systems,Society for Industrial and Applied Mathematics.The Open MPI Project, 2016, Open MPI: Open Source HighPerformance Computing, available online athttps://www.open-mpi.org/, Accessed: 10. 3. 2016.

Toselli, A., Windlund, O., 2005, Domain DecompositionMethods – Algorithms and Theory, Springer-Verlag.van Genuchten, M., 1980, A Closed-form Equation forPredicting the Hydraulic Conductivity of UnsaturatedSoils, Soil Science Society of America Journal, 40, 892-898.