ISRN Applied Mathematics

Volume 2013 (2013), Article ID 498383, 7 pages

http://dx.doi.org/10.1155/2013/498383

## A Broken -Nonconforming Finite Element Method for Incompressible Miscible Displacement Problem in Porous Media

^{1}Department of Mathematics, Shandong Jiaotong University, Jinan 250023, China^{2}College of Mathematical Science, Shandong Normal University, Jinan 250014, China

Received 10 June 2013; Accepted 31 July 2013

Academic Editors: G. Kyriacou and M. Ohlberger

Copyright © 2013 Fengxin Chen and Huanzhen Chen. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

An approximate scheme is defined for incompressible miscible displacement in porous media. This scheme is constructed by using immersed interface finite element method for the pressure equation which is based on the broken -nonconforming piecewise linear polynomials on interface triangular elements and utilizing finite element method for the concentration equation. Error estimates for pressure in broken norm and for concentration in norm are presented.

#### 1. Introduction

Miscible displacement of one incompressible fluid in a porous medium over time interval is modeled by the system

In the formulation, we have utilized the fact that in applications the vertical dimension of a subsurface geological formation is often much smaller than the horizontal dimension. We thus simplify the problem via a vertical average to rewrite the problem as a two-dimensional problem in the horizontal dimension, so the physical domain . , is the permeability tensor of the medium and is the viscosity of the fluid mixture that may be discontinuous across some interfaces; is the porosity of the medium; is the Darcy velocity of the mixture; represents flow rates at wells, commonly a linear combination of Dirac measures; is the diffusion-dispersion tensor, with , , and being the molecular diffusion, the transverse, and longitudinal dispersivities, respectively, is the identity tensor; is specified at sources and at sinks. is the initial concentration. The dependent variable is is the pressure in the fluid mixture, and is the concentration of a solvent injected into resident reservoir.

For (1) we assume that

Let ; for some restrict the variable to lie between and , so we can suppose that the following regularity for , , , and holds:

There exist much literature concerning numerical method and numerical analysis of miscible placement problem in porous media, for example, [1–4], and so on. As we know the porous media equations used to model the interface between oil and an injected fluid in simulations of secondary recovery in oil reservoirs. For the pressure equation (1)(a), the coefficient often changes rapidly across fluid interfaces, and this sharp change is accompanied by large changes in the pressure gradient, as a compensatory, yielding a fairly smooth Darcy velocity , so we can trade (1)(a) as an interface problem. When the interface is smooth enough, the solution of the interface problem is also very smooth in individual regions where the coefficient is smooth, but due to the jump of the coefficient across the interface, the global regularity is usually very low and has order of , . Due to the low global regularity and the irregular geometry of the interface, it seems to be difficult for the standard finite element method to achieve high accuracy.

For better approximation, the fitted finite element methods whose mesh depends on the smooth interface are developed [5]. However, this method using fitted grids is costly for more complicated time dependent problems in which the interface moves with time and repeated grid generation is called for. Compared with the fitted finite element methods, the immersed interface method proposed by LeVeque and Li [6] allows the mesh to be independent of the interface, such as a Cartesian mesh. In recent years, Li et al. [7] studied an immersed finite element method using uniform grid, and they proved the approximation property of the finite element space of this scheme. On the other hand, Kwak et al. [8] introduced an immersed finite element method based on the broken -nonconforming piecewise linear polynomials on interface triangular elements; this method uses edge averages as degrees of freedom, and the basis functions are C-R type [9]. Theory analysis and numerical experiments also show the optimal-order convergence of the method.

In this paper, we apply -nonconforming finite element method (FEM) to pressure equation, while a FEM was used to approximate the concentration equation. Other methods could also be employed for the discretization of the concentration equation, for example, characteristic Galerkin method and so on. However, since the main point of this paper is to show the feasibility of the use of FEM for pressure, we will discuss the concentration equation in the single case.

The rest of the paper is organized as follows: in Section 2, we first introduce some preliminaries. In Section 3, we briefly describe the FEM-FEM scheme. In Section 4, we present some projections and lemmas used in the following error analysis. In Section 5, we prove the main error estimate.

#### 2. Preliminaries

In this paper, without loss of generality, we consider the case in which is a rectangular domain and the interface is a smooth curve separating into two subdomains ** and such that ; see Figure 1 for an illustration.

The proper jump conditions for elliptic (1)(a) are given by where . We assume is a positive function bounded below and above by two positive constants.

Let , , be the Sobolev spaces consisting of functions whose derivatives up to order- are th integrable on , and and .

For the analysis, we introduce the space equipped with the norm

We integrate (1)(a) multiplied by any test function over the domain and apply the divergence theorem; then the variational formulation of the pressure equation is as follows Find such that where the bilinear formation .

By the Sobolev embedding theorem, for any , we can get for any . Then we have the following regularity lemma for the weak solution of the variational formulation (8).

Lemma 1 (see [10]). *Assume that . Then there exists a unique solution to the variational formulation (8) such that
*

#### 3. Formation of the Method

In this section, firstly, we define a broken -nonconforming finite element method for pressure equation, and then, we use a finite element method to approximate the concentration equation.

In order to construct a broken -nonconforming finite element procedure for pressure equation (1)(a), we assume the following situation.

Let be the usual regular triangulation of the domain such that the elements have diameters bounded by . Without loss of generality, we assume that the triangles in the partition used have the following features.()If meets one edge at more than two points, then this edge is one part of . ()If meets a triangle at two points, then the two points must be on the different edges of this triangle. ()The interface curve is defined by a piecewise function, and the mesh is formed such that the subset of in any interface element is .

Then we can separate the triangles of partition into two classes. For an element , if the interface passes through the interior of , we call it an interface element and denote it by ; otherwise, we call it a noninterface element and denote it by , respectively.

As usual, we will define linear finite element spaces on each element of the partition , so we have to construct the local basis functions. For a noninterface element , we can simply use the standard linear shape functions on with degrees of freedom at average values along edges of and construct the linear finite element space as follows:

For this space, we define the interpolation operator , where the following well-known approximation property is satisfied: and we use to denote the space of piecewise -nonconforming finite element space on the domain .

For interface element , since the finite element space on a general element can be obtained from the counterpart of a reference element through an affine mapping, we consider a typical interface element whose geometric configuration is given in Figure 2 in which the three vertices are given by , , and , and the curve between points and is part of the interface with for and for . Let denote the line segment connecting points and , which divides into two parts and with .

Let , be the edges of , and let denote the average of along , that is, , , and then we would like to construct a new function which is linear on and , respectively, and satisfies the jump condition (5) on . For this purpose, we write the modified basis function on an interface element as follows: with the following constraints: where , are averages along and , , , are given values, and is the unit normal vector on the line , that is, .

By similar calculation to that given in Theorem 2.2 of [8], we can rewrite (12) and (13) in the matrix form: where the coefficient matrix is defined by and the determinant of the matrix is

Therefore, the coefficients of (12) are uniquely determined by conditions , , respectively, and when , , have the same value , we can get the piecewise linear function by uniqueness.

Therefore we can define the following finite element space on an interface element :

Also we use to denote the finite element space defined on the domain .

*Remark 2. *Although for functions in , the flux jump condition is enforced on line segments, they actually satisfy a weak flux jump condition along the interface when is a piecewise constant such that
which is proved by an application of divergence theorem as in [11].

Combining the definitions of and ,we can describe the immersed finite element space on the whole domain , endowed with the broken norm,

Now, assume that the approximation of concentration is known. Then, the pressure can be determined by the following system: where is the bilinear formulation defined on with .

The question at hand is to discretize the concentration equation.

Multiplying the concentration equation (1)(b) by a test function and integrating over which leads to the following weak formulation for concentration :

Let denote the standard piecewise linear finite element space associated with , and denote the bound of elements diameters. Then the discrete procedure of (22) is to find such that where is the elliptic projection of .

#### 4. Projection and Some Lemmas

For convergence analysis conducted in the following section, we need an interpolation operator using the average of on each edge by and when , we can define by . Then the following interpolation estimate hold for the FEM approximation [8]

Next, let be the elliptic projection of defined by

The function will be chosen to assure coercivity of the form. Then, the following estimates hold [12]:

For convergence analysis conducted in Section 5, we need to apply quote some lemmas from [9, 13] and references therein.

Lemma 3 (see [9] (the second Strang lemma)). *Let and be the solution of (8) and (21). Then there exists a constant such that
*

Lemma 4 (see [13]). *Let be an edge of . Then there exists a constant such that for all :
**
where .*

#### 5. Convergence Analysis

In this section, we will present convergence analysis for the pressure and concentration.

Theorem 5. *Let , be the solution of (8) and (21); then there exists a constant independent of and the location of the interface, such that the following result holds:
*

*Proof. *Since the immersed finite element formulation (21) is nonconforming, we can use Lemma 3 to prove the error bound. The first part in (28) is an approximation error, which can be estimated simply:

For the second part, it is a consistency error estimate. By the definition of and Green’s formula, we can get
where is the unit outward normal vector on each , , and by the construction of the space we have a well-defined property on the interior edges, that is, for is the common edge of adjacent element and . By boundary condition, we note that vanishing average on the boundary. So we rewrite and use Lemma 4 to derive that

Combining (28) with (31)–(33) yields the theorem result.

We now are in the position to prove the error estimate for concentration.

Theorem 6. *Assume that the true solution of (1) satisfies , . Let be the projection of satisfying (26) and let be the solution of (23). Suppose that the mesh parameters satisfy . Then the following result holds:
*

*Proof. *We decompose the error to , and . Then, combining (22), (23), and (26) results in the following error equation for the concentration:

We choose the test function as in (35) and integrate from 0 to on both sides to derive

The terms on the left side can be easily bounded by

Next, we should estimate terms on the right side one by one.

By (27), it is trivial that

Similarly, we can bound the second and third terms as follows:
provided that is Lipschitz continuous.

For the fifth term, we can derive the estimate by using the boundness of and (27), as well as (30),

Similarly, we can derive

In order to prove the forth term, we need the following induction hypothesis:

As the statement in [1], we assume that for some ,
and with this hypothesis and (4), we can get

Therefore, collecting all the bounds derived above and using Gronwall inequality leads to

To complete the argument, we have to check hypothesis (42). We use (25), (27), (30), (45), inverse inequality, and boundness of to get that

Combining (27) with (30) and (45), we deduce the error estimates for the pressure and the concentration.

Theorem 7. *Assume that the true solution of (1) satisfies , . Let , be the solution of (23) and (21). Suppose that the mesh parameters satisfy . Then the following result holds:
*

*Remark 8. *In this paper, we just get a priori error estimates for the coefficient . However, some mathematical models for miscible displacements, which are currently being used by oil companies, make the assumption that the coefficient and it does not depend on ; with this assumption, we can prove the optimal order convergence similarly.

#### Acknowledgments

This work is supported in part by the National Natural Science Foundation of China (no. 10971254), National Natural Science Foundation of Shandong Province (no. Y2007A14), and Young Scientists Fund of Shandong Province (no. 2008BS01008).

#### References

- R. E. Ewing and M. F. Wheeler, “Galerkin methods for miscible displacement problems in porous media,”
*SIAM Journal on Numerical Analysis*, vol. 17, no. 3, pp. 351–365, 1980. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Q. Li and H. Z. Chen, “Finite element method for incompressible miscible displacement in porous media,”
*Acta Mathematicae Applicatae Sinica*, vol. 9, no. 3, pp. 204–212, 1993. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Douglas Jr., R. E. Ewing, and M. F. Wheeler, “A time-discretization procedure for a mixed finite element approximation of miscible displacement in porous media,”
*RAIRO Analyse Numérique*, vol. 17, no. 3, pp. 249–265, 1983. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. E. Ewing, T. F. Russell, and M. F. Wheeler, “Convergence analysis of an approximation of miscible displacement in porous media by mixed finite elements and a modified method of characteristics,”
*Computer Methods in Applied Mechanics and Engineering*, vol. 47, no. 1-2, pp. 73–92, 1984. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - I. Babuška, “The finite element method for elliptic equations with discontinuous coefficients,”
*Computing*, vol. 5, pp. 207–213, 1970. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Z. Li, T. Lin, Y. Lin, and R. C. Rogers, “An immersed finite element space and its approximation capability,”
*Numerical Methods for Partial Differential Equations*, vol. 20, no. 3, pp. 338–367, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Z. Li, T. Lin, and X. Wu, “New Cartesian grid methods for interface problems using the finite element formulation,”
*Numerische Mathematik*, vol. 96, no. 1, pp. 61–98, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. Y. Kwak, K. T. Wee, and K. S. Chang, “An analysis of a broken ${P}_{1}$-nonconforming finite element method for interface problems,”
*SIAM Journal on Numerical Analysis*, vol. 48, no. 6, pp. 2117–2134, 2010. View at Publisher · View at Google Scholar · View at MathSciNet - P. G. Ciarlet,
*The Finite Element Method for Elliptic Problems*, vol. 4, North-Holland Publishing, Amsterdam, The Netherlands, 1978. View at MathSciNet - J. H. Bramble and J. T. King, “A finite element method for interface problems in domains with smooth boundaries and interfaces,”
*Advances in Computational Mathematics*, vol. 6, no. 2, pp. 109–138, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - T. F. Russell and M. F. Wheeler, “Finite element and finite difference methods for continuous flows in porous media,” in
*The Mathematics of Reservoir Simulation*, R. E. Ewing, Ed., pp. 35–106, SIAM, Philadelphia, Pa, USA, 1984. View at Google Scholar - J. Douglas Jr., R. E. Ewing, and M. F. Wheeler, “The approximation of the pressure by a mixed method in the simulation of miscible displacement,”
*RAIRO Analyse Numérique*, vol. 17, no. 1, pp. 17–33, 1983. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Crouzeix and P.-A. Raviart, “Conforming and nonconforming finite element methods for solving the stationary Stokes equations. I,”
*RAIRO-Analyse Numérique*, vol. 7, no. 3, pp. 33–75, 1973. View at Google Scholar · View at MathSciNet