Abstrac tFlow modelling in heterogeneous media is a great challenge, especially in the presence of welk or hydraulic fracturen . The most used numerical scheme, finite volume method, allows taking heterogeneities into account . But finite volume meshes adapted to the geometry of welk or fracturen cannot be built efficiently . Conversely, Boundary element method (BEM) can be used to model complex geometries in homogeneous media . The use of BEM requires working on homogeneized area and accepting to loose the information about heterogeneities . In this paper, we propose an innovative method to solve this conflict and to use the combined advantages of both methods : we develop an approach using the coupling of finite volume scheme and boundary element method . The BEM is only used in reduced areas of complex geometries . The existence and uniqueness of the solution of our coupling scheme is mathematically proved and its performance is numerically studied . Finally, we show how our coupling scheme can bring a solution to the problem of modelling objects of complex geometry in heterogeneous media .
Introductio nBEM han been used in well tests in reservoir engineering [5,9] . As example, Pecher and Stanislas [9] proposed a resolution of a multi-domain integral methods to simulate well tests . The main interest of the boundary integral method is its ability to describe complex geometry . So it seems to be tempted to use boundary element methods around objects of complex geometry such as advanced welk or fractures and to use more classical discrete numerical scheme far from the well . In this paper, we consider a steady-state problem . One application of such procedure is the calculus of well performance in an heterogeneous reservoir . Another application could be the homogeneization of geological objects . Moumas et al .[7] propose a method to discretize welk using integral methods . In this paper, we deal with the problem of coupling a finite volume scheme and a boundary integral representation of an intertor domaio . The first part describes the simplified problem considered and presents the discretization in space . In the second part, we generalize this approach to several bounded and simply connected domains, where the solution is given by an integral representation . Part 3 is devoted to numerical tests .