Abstract

A general multidimensional numerical scheme, primarily designed to simulate two-phase compositional flow in porous media, is presented with serial and parallel implementations suitable for solving problems with degenerate (capillary) diffusion or capillary barrier effect in heterogeneous porous materials. The numerical scheme is based on the mixed–hybrid finite element method with the semi-implicitapproach for the time discretization in order to obtain a system of linear equations in each time step. The scheme is implemented in serial for CPU and in parallel for CPU and GPU using TNL that provides an efficient abstraction layer for accessing various parallel hardware architectures. In order to demonstrate the applicability of the numerical scheme, a numerical analysis is presented for problems of two-phase flow in 1D, 2D, and 3D for which exact (semi-analytical) solutions are known and a series of benchmark problems for two-phase flow in heterogeneous porous media is discussed to show correct simulation of the capillary barrier effect. The efficiency and accuracy of the implementations on CPU and GPU are discussed. Moreover, we construct an analytical solution and use it to demonstrate convergence of the numerical scheme for two-phase compositional flow problems in porous media.

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.