Numerical modelling of porous flow in a low-permeability matrix with high-permeability inclusions is a challenging task because the large ratio of permeabilities ill-conditions the finite element system of equations. We propose a coupled model where Darcy flow is used for the porous matrix and potential flow is used for the inclusions. We discuss appropriate interface conditions in detail and show that the head drop in the inclusions can be prescribed in a very simple way. Algorithmic aspects are treated in full detail. Numerical examples show that this coupled approach precludes ill-conditioning and is more efficient than heterogeneous Darcy flow.
This is the pre-peer reviewed version of the following article: Rodríguez, A.; Sarrate, J.; Huerta, A. Numerical modelling of void inclusions in porous media. "International journal for numerical methods in engineering", Gener 2004, vol. 59, núm. 4, p. 577-596, which has been published in final form at http://www3.interscience.wiley.com/journal/106570919/abstract
CitationRodríguez, A.; Sarrate, J.; Huerta, A. Numerical modelling of void inclusions in porous media. "International journal for numerical methods in engineering", Gener 2004, vol. 59, núm. 4, p. 577-596.
All rights reserved. This work is protected by the corresponding intellectual and industrial property rights. Without prejudice to any existing legal exemptions, reproduction, distribution, public communication or transformation of this work are prohibited without permission of the copyright holder. If you wish to make any use of the work not provided for in the law, please contact: firstname.lastname@example.org