Analysis of Two-Grid Methods for Nonlinear Parabolic Equations by Expanded Mixed Finite Element Methods

Transcription

1 Advances in Applied athematics and echanics Adv. Appl. ath. ech., Vol. 1, No. 6, pp DOI: /aamm.09-m09S09 December 009 Analysis of Two-Grid ethods for Nonlinear Parabolic Equations by Expanded ixed Finite Element ethods Yanping Chen 1,, Peng Luan and Zuliang Lu 3 1 School of athematical Science, South China Normal University, Guangzhou , Guangdong, China athematics Teaching and Research Group, Sanshui Experimental High School of Zhongshan University, Foshan 58145, Guangdong, China 3 Hunan Key Lab for Computation and Simulation in Science and Engineering, Department of athematics, Xiangtan University, Xiangtan , Hunan, China Received 15 April 009; Accepted in revised version) 08 October 009 Available online 18 November 009 Abstract. In this paper, we present an efficient method of two-grid scheme for the approximation of two-dimensional nonlinear parabolic equations using an expanded mixed finite element method. We use two Newton iterations on the fine grid in our methods. Firstly, we solve an original nonlinear problem on the coarse nonlinear grid, then we use Newton iterations on the fine grid twice. The two-grid idea is from Xu s work [SIA J. Numer. Anal., ), pp ] on standard finite method. We also obtain the error estimates for the algorithms of the two-grid method. It is shown that the algorithm achieve asymptotically optimal approximation rate with the two-grid methods as long as the mesh sizes satisfy h = OH 4k+1)/k+1) ). AS subject classifications: 65N30, 65N15, 651 Key words: Nonlinear parabolic equations, two-grid scheme, expanded mixed finite element methods, Gronwall s Lemma. 1 Introduction In this paper, we consider the following nonlinear parabolic equations p t Kp) p) = f p, p), x, t) Ω J, 1.1) Corresponding author. URL: info.asp?usernamesp=%b3%c%d1%de%c6%bc yanpingchen@scnu.edu.cn Y. Chen), 341lp@163.com P. Luan), zulianglux@16.com Z. Lu) c 009 Global Science Press

2 Y. Chen, P. Luan, Z. Lu / Adv. Appl. ath. ech., 6 009), pp with initial condition and boundary condition px, 0) = p 0 x), x Ω, 1.) Kp) p) ν = 0, x, t) Ω J, 1.3) where Ω R is a bounded and convex domain with C 1 boundary Ω, ν is the unit exterior normal to Ω, J = 0, T], K is a symmetric positive definite tensor and K : Ω R R. Eq. 1.1) can be rewritten as following p + Ψ = f p, p), t 1.4) Kp) 1 Ψ + p = ) Eqs. 1.1)-1.3) are the simplification of the modeling of groundwater through porous media [9]. In this particular case, p denotes the fluid pressure; and K is a symmetric, uniformly positive definite tensor with L Ω) components representing the permeability divided by the viscosity; Ψ represents the Darcy velocity of the flow; and f p, p) models the external flow rate. The two-grid methods was first introduced by Xu [16, 17] as a discretization technique for nonlinear and nonsymmetric indefinite partial differential equations. It based on the fact the nonlinearity, nonsymmetry and indefiniteness behaving like low frequencies are governed by coarse grid and the related high frequencies are governed by some linear or symmetric positive definite operators. The basic idea of the two-grid method is to solve a complicated problem nonlinear, nonsymmetric indefinite) on a coarse grid mesh size H) and then solve an easier problem linear, symmetric positive) on the fine grid mesh size h and h H) as correction. In many partial differential equations, the objective functional contains the gradient of the state variables. Thus, the accuracy of the gradient is important in numerical discretization of the coupled state equations. ixed finite element methods are appropriate for the equations in such cases since both the scalar variable and its flux variable can be approximated to the same accuracy by using such methods. Some specialists have made many important works on some topic of mixed finite element method for linear elliptic or reaction-diffusion equations. In [5, 6], Chen has studied the expanded mixed element methods for some quasilinear second order elliptic equations. However, there doesn t seem to exist much work on theoretical analysis for two-grid methods for mixed finite element approximation of quasilinear or nonlinear parabolic equations in the literature. any contributions have been done to the multi-grid schemes for finite element methods, see, for example [11,1]. In [8], the authors have studied a two-grid finite difference scheme for nonlinear parabolic equations. Xu and Zhou have considered some multi-scale schemes for finite element method of elliptic partial differential equations in [18]. Recently, we constructed a new two-grid method of expanded mixed finite

3 83 Y. Chen, P. Luan, Z. Lu / Adv. Appl. ath. ech., 6 009), pp element method for semi-linear reaction-diffusion equation [4]. As a continued work of Chen [3], in this paper, we shall discuss the case that the coefficient matrix function K and the reaction term f are nonlinear, namely, K = Kp) and f = f p, p). In this paper, we adopt the standard notation W m,p Ω) for Sobolev spaces on Ω with a norm m,p given by and a semi-norm m,p given by v p m,p = v p m,p = D α v p L p Ω), α m D α v p L p Ω). α =m We set For p=, we denote W m,p 0 Ω) = { v W m,p Ω) : v Ω = 0 }. H m Ω) = W m, Ω), H m 0 Ω) = W m, 0 Ω), m = m,, = 0,. We denote by L s J; W m,p Ω)) the Banach space of all L s integrable functions from J into W m,p Ω) with norm T ) 1 v Ls J;W m,p Ω)) = v s W m,p Ω) dt s, for s [1, ), 0 and the standard modification for s =. The details can be found in [10]. In addition C or c denotes a general positive constant independent of h. At first, we make the following assumptions for the nonlinear parabolic equations. A1) There exist positive constants K and K, such that for z R, K z z t Kx, s)z K z, for x Ω, and that each element of K is twice continuously differential with derivatives up to second order bounded above by K. A) For some integer k 1, we assume that the solution function 1.4)-1.5) has the following regularity: ) p L J; W k+,4 Ω), Ψ L J; H Ω))) k+1. To analyze the discretization on a time interval 0, T), let N > 0, = T/N, and

4 Y. Chen, P. Luan, Z. Lu / Adv. Appl. ath. ech., 6 009), pp t n = n, and set φ n = φ, t n ), t φ n = φn φ n 1, N 1 φ l 0,T);X) = φ n X), φ l 0,T);X) = max 1 n N φn X, T ) 1 ϕ L 0,T);X) = ϕ n, t) X dt. 0 The plan of this paper is as follows. In next section, we introduce some notations and projections. The approximation properties of these projections will be recalled. In section 3, we construct a expanded mixed finite element discretization for the nonlinear parabolic equations 1.1)-1.3). A two-grid scheme and error estimates for twodimensional nonlinear parabolic equations using an expanded mixed finite element method are discussed in section 4. Finally, we give the conclusion and future work in section 5. Projection operators and approximation properties Let, ) denote the L Ω) inner product. Let V = Hdiv; Ω) = {v L Ω)), v L Ω)}, ṼV = V {v ν = 0}, W = L Ω). The Hilbert space V is equipped with the following norm: v Hdiv;Ω) = v + divv ) 1. Let Γ h denote a quasi-uniform partition of Ω into rectangles or triangles with the partition step h. We form V h and W h, discrete subspaces of V and W, using standard mixed finite element space such as the RT [13] spaces of order k, RT k, or Brezzi-Douglas- arini [1] spaces of order k, BD k. In addition, we assume that A3) The following inclusion hold for the RT spaces or BD spaces v h W h, v h V h. We shall employ some projection operators. Let Q h denote the L projection defined by φ, w h ) = Q h φ, w h ), w h W h,.1)

5 834 Y. Chen, P. Luan, Z. Lu / Adv. Appl. ath. ech., 6 009), pp for any φ L Ω), and φ, v h ) = Q h φ, v h ), v h V h,.) for any vector valued function φ L Ω)). We also assume that the following approximation properties hold. For 1 < q, there exists a positive constant C independent of h such that [1]: for all φ W k+1,q Ω) or φ W k+1,q Ω)) ) Q h φ 0,q C φ 0,q, q <,.3) φ Q h φ 0,q C φ r,q h r, 0 r k ) We use the well known Π h projection for mixed finite element approximation spaces. We shall assume that there exists a projection operator Π h : H 1 Ω)) V h such that for q Hdiv; Ω), Π h q, w h ) = q, w h ), w h W h..5) The following approximation properties [1] hold for the projection Π h : We also assume that q Π h q 0,q C q r,q h r, 1 q < r k + 1,.6) q Π h q) 0,q C q r,q h r, 0 r k ) φ Q h φ C φ h r, 0 r k + 1,.8) v h C v h h 1, v h V h..9) Now, we recall the discrete Gronwall s Lemma see, e.g., [14]): Lemma.1. Assume that k n is a non-negative sequence, and that the sequence ϕ n satisfies ϕ 0 g 0, n 1 n 1.10) ϕ n g 0 + p s + k s ϕ s, n 1. s=0 s=0 Then ϕ n satisfies ϕ 1 g k 0 ) + p 0, n 1 ϕ n g 0 s=0 n 1 + k s ) + s=0 n 1 p s r=s k r ) + p n 1, n..11) oreover, if g 0 0 and p n 0 for n 0, it follows ) ) n 1 n 1 ϕ n g 0 + p s exp k s, n 1..1) s=0 s=0

6 Y. Chen, P. Luan, Z. Lu / Adv. Appl. ath. ech., 6 009), pp Expanded mixed finite element discretization First, we set three variables: the pressure p, the gradient Υ= p, and the flux Ψ=Kp)Υ. We suppose that the following assumption is satisfied. A4) The reaction term f is a sufficiently smooth function with bounded derivatives through the second order on Ω. For all ˆp, ˆΥ) W V there exists F > 0 such that f p ˆp, ˆΥ) F, f Υ ˆp, ˆΥ) F. Now, we define the weak form of the nonlinear parabolic equations 1.1)-1.3) as follows: Find p, Υ, Ψ) W V ṼV such that p ) t, w + Ψ, w) = f p, Υ), w ), w W, 3.1) Υ, v) = p, v), v ṼV, 3.) Ψ, v) = Kp)Υ, v ), v V. 3.3) The discrete time mixed finite element approximation to 3.1)-3.3) can be defined as follows: Find p n h, Υn h, Ψn h ) W h V h ṼV h such that p n h p n 1 ) h, w h + Ψ n h, w h) = f p n h, Υn h ), w h), wh W h, 3.4) Υ n h, v h ) = p n h, v h ), v h ṼV h, 3.5) Ψ n h, v h ) = Kp n h )Υn h, v h ), vh V h. 3.6) At the time of t = t n, we rewrite 3.1)-3.3) in the following form p n ) t, w h + Ψ n, w h ) = f p n, Υ n ) ), w h, wh W h, 3.7) Υ n, v h ) = p n, v h ), v h ṼV h, 3.8) Ψ n, v h ) = Kp n )Υ n, v h ), vh V h. 3.9) Using the definition of Π h and Q h, and the assumption A3) that V h W h, we can obtain Qh p n Q h p n 1 ), w h + Π h Ψ n, w h ) = f p n, Υ n ) ), w h + t p n p n t, w h ), 3.10) Q h Υ n, v h ) = Q h p n, v h ), 3.11) Π h Ψ n, v h ) = Kp n )Υ n, v h ) + Πh Ψ n Ψ n, v h ). 3.1) Set µ n = Q h p n p n h, ξn = Q h Υ n Υ n h, λn = Π h Ψ n Ψ n h.

7 836 Y. Chen, P. Luan, Z. Lu / Adv. Appl. ath. ech., 6 009), pp Subtracting 3.10)-3.1) from 3.4)-3.6), we derive the error equations: µ n µ n 1 ), w h + λ n, w h ) = f p n, Υ n ) f p n h, Υn h ), w h) + t p n p n t, w h ), 3.13) ξ n, v h ) = µ n, v h ), 3.14) λ n, v h ) = Kp n )Υ n Kp n h )Υn h, v h ) + Πh Ψ n Ψ n, v h ). 3.15) Letting w h = µ n, v h = λ n, v h = ξ n in Eqs. 3.13)-3.15) respectively, we can get µ n µ n 1 ), µ n + λ n, µ n ) = f p n, Υ n ) f p n h, Υn h ), µn) + t p n p n t, µ n ), 3.16) ξ n, λ n ) = µ n, λ n ), 3.17) λ n, ξ n ) = Kp n h )Υn Q h Υ n ), ξ n) Kp n h ) Kpn ))Υ n, ξ n) + Kp n h )ξn, ξ n) + Π h Ψ n Ψ n, ξ n ). 3.18) Noticing that 1 µ n 1 + µ n ) µ n 1, µ n ), 3.19) Kp n h )ξ n, ξ n) = Kp n h ) 1 ξ n. 3.0) We have t µ n, µ n ) + Kp n h )ξn, ξ n) 1 Applying the Taylor expansions, we have µ n µ n 1 ) + Kp n h ) 1 ξ n. 3.1) f p n, Υ n ) f p n h, Υn h ) = f p p, Ῡ)p n p n h ) f Υ p, Ῡ)Υ n Υ n h ). Then, f p n, Υ n ) f p n h, Υn h ), µn) Fp n Q h p n ), µ n) + Fµ n, µ n) FΥ n Q h Υ n ), µ n) Fξ n, µ n ). 3.) oreover, t p n p n t, µ n ) ) = t t n 1 )p tt, s)ds, µ n t n t n p tt, s) ds µ n ) p, s) t n 1 t n 1 t ) + µ n, 3.3)

8 Y. Chen, P. Luan, Z. Lu / Adv. Appl. ath. ech., 6 009), pp and Kp n h ) Kp n ))Υ n, ξ n) K Υ n p n h pn ξ n CK p n h Q hp n + Q h p n p n ) ξ n = CK µ n + Q h p n p n ) ξ n C µ n + Q h p n p n ) + δ ξ n, 3.4) where we have used the assumptions A1)-A). Hölder s inequality gives From 3.16)-3.4), applying the 1 µn µ n 1 ) + Kp n h ) 1 ξ n µ n µ n 1 ), µ n + Kp n h ) 1 ξ n t p n p n t, µ n ) + Π h Ψ n Ψ n, ξ n ) + Fp n Q h p n ), µ n ) + Fµ n, µ n ) FΥ n Q h Υ n ), µ n ) Fξ n, µ n ) + Kp n h )Υn Q h Υ n ), ξ n) + Kp n h ) Kpn ))Υ n, ξ n) C µ n + C + δ ξ n + C[ Π h Ψ n Ψ n + Υ n Q h Υ n + p n Q h p n ]. 3.5) We next multiply both sides of 3.5) by and sum over n, n = 1,, N, applying the discrete Gronwall s Lemma, to obtain C µ N + N N Kp n h ) 1 ξ n [ Π h Ψ n Ψ n + Υ n Q h Υ n + p n Q h p n ] + µ 0 + C). With a proper choice of the initial function p 0 h = Q hp 0, we have µ 0 = 0, combining.4) and.6), we get Then, we have µ N + µ N + N Kp n h ) 1 ξ n C ) + h k+). N Kp nh ) 1 ξ n ) 1 C + h k+1). 3.6) Applying the triangle inequality, we have the following result.

9 838 Y. Chen, P. Luan, Z. Lu / Adv. Appl. ath. ech., 6 009), pp Theorem 3.1. Suppose that assumptions A1)-A4) are valid. Let p n h, Υn h, Ψn h ) W h V h ṼV h be the solution of the mixed finite element equations 3.4)-3.6) for n 1. If we choose the initial function p 0 h = Q hp 0, for 1 N and k 1, then there exists a positive constant C, independent of h such that ) 1 p ph + Kp n h ) 1 ) 1 Υ n Υ n h ) + Ψ n Ψ nh C + h k+1). 3.7) 4 A two-grid method and error estimates Note that the two-grid algorithms used in Wu and Allen [15] involve only one Newton iteration on the fine grid. Nevertheless, some more interesting results can be derived if one more Newton iteration is carried out on the fine grid. The idea is from Xu s work [16]. Thus, based on the works of Dawson [7], we can construct an algorithm of the two-grid method for Eqs. 1.1)-1.3) discretized by expanded mixed finite element methods. To iteratively solve the nonlinear system 3.4)-3.6), we find p n h, Υ n h, Ψ n h ) W h V h ṼV h in three steps as follows: Step 1: On the coarse grid Γ H, compute p n H, Υn H, Ψn H ) W H V H ṼV H to satisfy the following original nonlinear system: ) p n H p n 1 H, w H + Ψ n H, w H) = f n, w H ), w H W H, 4.1) Υ n H, v H ) p n H, v H ) = 0, v H ṼV H, 4.) Ψ n H, v H ) = Kp n H)Υ n H, v H ), v H V H. 4.3) Step : On the fine grid Γ h, compute p n h, Υ n h, Ψ n h ) W h V h ṼV h to satisfy the following linear system: ) p n h p n 1 h, w h + Ψ n h, w h) = f n, w h ), w h W h, 4.4) Υ n h, v h ) p n h, v h ) = 0, v h ṼV h, 4.5) Ψ n h, v h ) = Kp n H) Υ n h, v h ) + K p n H)Υ n H p n h pn H), v h ), v h V h. 4.6)

10 Y. Chen, P. Luan, Z. Lu / Adv. Appl. ath. ech., 6 009), pp Step 3: On the fine grid Γ h, compute p n h, Υ n h, Ψ n h ) W h V h ṼV h to satisfy the following linear system: p n h ) p n 1 h, w h + Ψ n h, w h) = f n, w h ), w h W h, 4.7) Υ n h, v h ) p n h, v h ) = 0, v h ṼV h, 4.8) Ψ n h, v h ) = K p n h ) Υ n h, v h ) + K p n h ) Υ n h p n h pn h ), v h ), v h V h. 4.9) Remark 4.1. Eq. 4.9) is motivated by the Taylor expansion K p n h ) Υ n h = K pn h ) Υ n h + K p n h ) Υ n h p n h pn h ) + K x n ) for some x n between p n h and pn h. Υ n h p n h pn h ), 4.10) Lemma 4.1. Suppose that assumptions A1)-A4) are fulfilled. Let p n h, Υ n h, Ψ n h ) W h V h ṼV h be the solution of the mixed finite element equations 4.1)-4.6) for n 1. If we choose the initial function p 0 h = Q hp 0, for 1 N and k 1, then there exists a positive constant C, independent of h such that ) 1 p p h + Kp n H) 1 Υ n Υ n h ) C + h k+1 + H k+1 ). The proof can be found in [7]. oreover, by an argument similar to that in Lemma 4.1, we can also prove the following result: Remark 4.. If we choose the initial function p 0 h = Q hp 0, then for 1 n N, q <, and k 1, we have p n p n h 0,q C + h k+1 + H k+1 ). 4.11) Now, we can prove the following main theorem. Theorem 4.1. Suppose that assumptions A1)-A4) are fulfilled. Let p n h, Υ n h, Ψ n h ) W h V h ṼV h be the solution of the mixed finite element equations 4.1)-4.9) for n 1. If we choose the initial function p 0 h = Q hp 0, 4.1) then for 1 N and k 1, we have ) 1 p p h + K p n h ) 1 Υ n Υ n h ) C + h k+1 + H 4k+1 ). 4.13)

11 840 Y. Chen, P. Luan, Z. Lu / Adv. Appl. ath. ech., 6 009), pp Proof. We now derive an estimate for this two-grid scheme. We set α n = Q h p n p n h, βn = Q h Υ n Υ n h, γn = Π h Ψ n Ψ n h. Using 3.10)-3.1), 4.7)-4.9), and the definition of projection Q h, Π h, we get the error equations: and αn α n 1, w h ) + Π h Ψ n Ψ n h ), w h) = E p, w h ), w h W h, 4.14) Q h Υ n Υ n h, v h ) Q h p n p n h, v h ) = 0, v h ṼV h, 4.15) Π h Ψ n Ψ n h, v h ) ) =Π h Ψ n Ψ n, v h ) Kp n )Q h Υ n Υ n ), v h ) + K p n h )Q hυ n Υ n h ), v h ) ) + α n + p n Q h p n ))K p n h ) Υ n h, v h + K p n h )pn p n h )Q hυ n Υ n h ), v h + 1 ) p n p n h ) K x n )Q h Υ n, v h, v h V h, 4.16) for some x n between p n h and pn h and E p = t p n p n t. Choosing w h = α n in 4.14), v h = γ n in 4.15), v h = β n in 4.16), adding 4.14) and 4.15), and subtracting 4.16), yield α n α n 1 ), α n + K p n h )βn, β n) = E p, α n ) Π h Ψ n Ψ n, β n ) + Kp n )Q h Υ n Υ n ), β n) α n + p n Q h p n ))K p n h ) Υ n h, βn) K p n h )pn p n h )Q hυ n Υ n h ), βn) 1 p n p n h ) K x n )Q h Υ n, β n), 4.17) for some x n between p n h and pn h. Note that p n p n h ) K x n )Q h Υ n, β n) [ = p n h p n H) p n h pn )p n p n H) p n p n H) ] K x n )Q h Υ n Υ n + Υ n Υ n H), β n) + p n p n h ) K x n )Υ n H Υ n + Υ n ), β n), 4.18) and using the Hölder s inequality, we get p n p n h ) K x n )Υ n, β n) C p n p n h 0,4 βn. 4.19)

12 Y. Chen, P. Luan, Z. Lu / Adv. Appl. ath. ech., 6 009), pp Using the Hölder s inequality and 3.3), multiplying by and summing n = 1 to, N, using the approximation properties, we see that where α α 0 + C) + δ +C C) + δ K p n h ) 1 β n β n + C [ α n + h k+ ] ] [ p n Q h p n ) Υ n h + p n p n h )Q hυ n Υ n h ) + T 1 β n + C [ α n + h k+ + T ], 4.0) T 1 = p n H p n h p n H p n h Q h Υ n Υ n + p n p n h p n p n H Q h Υ n Υ n + p n p n H p n p n H Q h Υ n Υ n + p n p n h p n p n h Υ n H Υ n + p n p n h 4 0,4, 4.1) ] T =C [ p n Q h p n ) Υ n h + p n p n h )Q hυ n Υ n h ) + T 1. 4.) Using 4.11) and the Hölder s inequality, we have T 1 = p n H p n h p n H p n h Q h Υ n Υ n + p n p n h p n p n H Q h Υ n Υ n + p n p n H p n p n H Q h Υ n Υ n + p n p n h p n p n h Υ n H Υ n + p n p n h 4 0,4 h h k+1 + H k+1 + ) 4 h k+ + h h k+1 + H k+1 + ) H k+1 + ) h k+ + h h k+1 + H k+1 + ) 4 Υ n Υ n H + h k+1 + H k+1 + ) ) Combining 4.)-4.3), Lemma 4.1 and Theorem 3.1, we can obtain T =C C ] [ p n Q h p n ) Υ n h + p n p n h )Q hυ n Υ n h ) + T 1 [ h k+ + h h k+1 + H k+1 + ) 4 h k+ + h h k+1 + H k+1 + ) H k+1 + ) h k+ + h k+1 + H k+1 + ) 4]

13 84 Y. Chen, P. Luan, Z. Lu / Adv. Appl. ath. ech., 6 009), pp h h k+1 + H k+1 + ) 4 Υ n Υ n H. From.4),.8)-.9), and Theorem 3.1, we have and C Υ n Υ n H C Υ n Q H Υ n + Q H Υ n Υ n H ) H k+ + C H 1 Q H Υ n Υ n H ) CH H k+1 + ), 4.4) Υ n Υ n H Υ n Q H Υ n + Q H Υ n Υ n H H 1 H k+1 + ). 4.5) Then, we obatin T C h k+ + H 8k+ + H 10k+4 h + ) ), 4.6) where h H k+1, we have T C h k+ + H 8k+ + ) ). 4.7) Noticing that α 0 = 0, and applying the discrete Gronwall s Lemma, we deduce that α + K p nh ) 1 β n ) 1 Ch k+1 + H 4k+1 + ). 4.8) Applying the triangle inequality and the approximation properties, we can prove the results of the Theorem. 5 Conclusion and future works In this paper, we have presented a new two-grid scheme for expanded mixed finite element solution of nonlinear parabolic equations. The theory demonstrates a remarkable fact about the two-grid scheme: we can iterate on a very coarse grid Γ H and still get good approximations by taking one iteration on the fine grid Γ h. We have derived

14 Y. Chen, P. Luan, Z. Lu / Adv. Appl. ath. ech., 6 009), pp an error estimate for the algorithms of the two-grid method. It is shown that the algorithm achieve asymptotically optimal approximation as long as the mesh sizes satisfy ) h = O H 4k+1 k+1. The next step of our work is to carry out some the numerical experiment with the proposed two-grid scheme. Furthermore, we shall consider the two-grid schemes for expanded mixed finite element solution of nonlinear hyperbolic equations. Acknowledgements This work is supported by Guangdong Province Universities and Colleges Pearl River Scholar Funded Scheme 008), National Science Foundation of China , the National Basic Research Program under the Grant 005CB31703 and Hunan Provincial Innovation Foundation For Postgraduate CX009B119 References [1] F. BREZZI, J. DOUGLAS JR. AND L. D. ARINI, Two families of mixed finite elemetns for second order elliptic problems, Numer. ath., ), pp [] Y. CHEN, Y. HUANG AND D. YU, A two-grid method for expanded mixed finite element solution of semilinear reaction-diffusion equations, Internat. J. Numer. ethods in Engineering, ), pp [3] Y. CHEN AND L. LI, L p error estimates of two-grid schemes of expanded mixed finite element methods, Appl. ath. Comp., ), pp [4] Y. CHEN, D. YU AND S. LIU, Analysis of two-grid method for reaction-diffusion equations by expanded mixed finite element methods, Internat. J. Numer. ethods in Engineering, ), pp [5] Z. CHEN, Expanded mixed element methods for quasilinear second-order elliptic problemsi), RAIRO: odèl. ath. Anal. Numer., ), pp [6] Z. CHEN, Expanded mixed element methods for quasilinear second-order elliptic problemsii), RAIRO: odèl. ath. Anal. Numer., ), pp [7] C. N. DAWSON AND. F. WHEELER, Two-grid methods for mixed finite element approximations of non-linear parabolic equations, Contemporary ath., ), pp [8] C. N. DAWSON,. F. WHEELER AND C. S. WOODWARD, A two-grid finite difference scheme for nonlinear parabolic equations, SIA J. Numer. Anal., ), pp [9] J. DOUGLAS, JR., R. E. EWING AND. F. WHEELER, The approximation of the pressure by a mixed method in the simulation of miscible displacement, RAIRO Anal. Numér., ), pp [10] L. EVANS, Partial differential equations, Amer. ath. Soc., Providence, RI, [11] Z. LU AND H. ZHANG, V cycle multi-grid-cycle method of viscoelastic fluid flow obeying an Oldroyd B type constitutive law, Numer. Anal. Appl., ), pp [1] F. PRETORIUS AND. W. CHOPTUIK, Adaptive mesh refinement for coupled elliptichyperbolic systems, J. Comp. Phys., ), pp

15 844 Y. Chen, P. Luan, Z. Lu / Adv. Appl. ath. ech., 6 009), pp [13] P. A. RAVIART AND J.. THOAS, A mixed finite element method for nd order elliptic problems, ath. Aspects of the Finite Element ethod, Lecture Notes in ath., 606, Springer, Berlin, 1977), pp [14] V. THOÉE, Galerkin Finite Element ethods for Parabolic Problems, Springer-Verlag, [15] L. WU AND. B. ALLEN, A two-grid method for mixed finite-element solution of reactiondiffusion equations, Numer. ethods for Partial Diffenrential Equations, ), pp [16] J. XU, Two-grid discretization techniques for linear and non-linear PDEs, SIA J. Numer. Anal., ), pp [17] J. XU, A novel two-grid method for semilinear equations, SIA J. Sci. Comp., ), pp [18] J. XU AND A. H. ZHOU, Some multiscale methods for partial differential equations, Contemporary ath., ), pp. 1 7.