AAA Abstract and Applied Analysis 1687-0409 1085-3375 Hindawi Publishing Corporation 803615 10.1155/2014/803615 803615 Research Article A Two-Grid Finite Element Method for a Second-Order Nonlinear Hyperbolic Equation Chen Chuanjun 1 Liu Wei 2 Zhao Xin 1 Zhang Xinguang 1 School of Mathematics and Information Sciences Yantai University Yantai 264005 China ytu.edu.cn 2 School of Statistics and Mathematics Shandong University of Finance and Economics Jinan 250014 China sdu.edu.cn 2014 532014 2014 06 12 2013 28 01 2014 5 3 2014 2014 Copyright © 2014 Chuanjun Chen et al. 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.

We present a two-grid finite element scheme for the approximation of a second-order nonlinear hyperbolic equation in two space dimensions. In the two-grid scheme, the full nonlinear problem is solved only on a coarse grid of size H. The nonlinearities are expanded about the coarse grid solution on the fine gird of size h. The resulting linear system is solved on the fine grid. Some a priori error estimates are derived with the H1-norm O(h+H2) for the two-grid finite element method. Compared with the standard finite element method, the two-grid method achieves asymptotically same order as long as the mesh sizes satisfy h=O(H2).

1. Introduction

Let Ω2 be a bounded convex domain with smooth boundary Γ, and consider the initial-boundary value problem for the following second-order nonlinear hyperbolic equation (1)utt-·(A(u)u)=f(x,t),xΩ,0<tT,u(x,t)=0,xΓ,0<tT,u(x,0)=u0(x),ut(x,0)=u1(x),xΩ, where utt and ut denote 2u/t2 and u/t, respectively. x=(x1,x2). We assume that A(u) is a symmetric positive definite matrix. A(u) and Au(u) satisfy the Lipschitz continuous condition with respect to u, where Au=A/u and (2)|A(u1)-A(u2)|L|u1-u2|,(3)|Au(u1)-Au(u2)|L|u1-u2|,u1,u2, where L is a positive constant.

Two-grid method is a discretization technique for nonlinear equations based on two grids of different sizes. The main idea is to use a coarse-grid space to produce a rough approximation of the solution of nonlinear problems and then use it as the initial guess for the solution on the fine grid. This method involves a nonlinear solution on the coarse grid with grid size H and a linear solution on the fine grid with grid size h<H. Two-grid method was first introduced by Xu [1, 2] for linear (nonsymmetric or indefinite) and especially nonlinear elliptic partial differential equations. Later on, two-grid method was further investigated by many authors. Dawson and Wheeler [3, 4], Chen and Liu  have constructed the two-grid method by using finite difference method, mixed finite element method, and piecewise linear finite element method for nonlinear parabolic equations, respectively. Wu and Allen  have applied two-grid method combined with mixed finite element method to reaction-diffusion equations. Chen et al.  have constructed two-grid methods for expanded mixed finite-element solution of semilinear and nonlinear reaction-diffusion equations. Bi and Ginting  have studied two-grid finite volume element method for linear and nonlinear elliptic problems. Chen et al. , Chen and Liu [13, 14] have studied two-grid methods for semilinear parabolic and second-order hyperbolic equations using finite volume element method.

The finite element analysis for the second-order linear hyperbolic equations was discussed by Dupont  and Baker . They have obtained optimal L(L2) estimates for the error, O(hr), using subspaces of piecewise polynomial functions of degree r-1, for r1. Then Yuan and Wang [17, 18] have studied error estimates for the finite element method of the second-order nonlinear hyperbolic equations and proved the optimal error estimates in the L2 and H1 norm. Kumar et al.  presented and discussed semidiscrete piecewise linear finite volume approximations for a second-order wave equation and obtained optimal error estimates in L2, H1, and L norms. For second-order hyperbolic equations with a nonlinear reaction term, Chen and Liu  have presented a two-grid method using finite volume element method and obtained error estimate in the H1-norm.

However, as far as we know there is no two-grid finite element convergence analysis for the second-order nonlinear hyperbolic equations (1). In this paper, based on two conforming piecewise linear finite element spaces SH and Sh on one coarse grid with grid size H and one fine grid with grid size h, respectively, we consider the two-grid finite element discretization techniques for the second-order nonlinear hyperbolic problems. With the proposed techniques, solving the nonlinear problems on the fine-grid space is reduced to solving a linear system on the fine-grid space and a nonlinear system on a much smaller space. This means that solving a nonlinear problem is not much more difficult than solving one linear problem, since dimSHdimSh and the work for solving the nonlinear problem is relatively negligible. A remarkable fact about this simple approach is, as shown in , that the coarse mesh can be quite coarse and still maintain a good accuracy approximation.

The rest of this paper is organized as follows. In Section 2, we describe the finite element scheme for the nonlinear second-order hyperbolic problem (1). Section 3 contains the error estimates for the finite element method. Section 4 is devoted to the two-grid finite element and its error analysis. Throughout this paper, the letter C or with its subscript denotes a generic positive constant which does not depend on the mesh parameters and may be different at its different occurrences.

2. Standard Finite Element Method

We adopt the standard notation for Sobolev spaces Ws,p(Ω) with 1p consisting of functions that have generalized derivatives of order s in the space Lp(Ω). The norm of Ws,p(Ω) is defined by (4)us,p,Ω=us,p=(Ω|α|s|Dαu|pdx)1/p, with the standard modification for p=. In order to simplify the notation, we denote Ws,2(Ω) by Hs(Ω) and omit the index p=2 and Ω whenever possible; that is, us,2,Ω=us,2=us. Let H01(Ω) be the subspace of H1(Ω) of functions vanishing on the boundary Γ.

For the variational formulation we multiply (1) by a smooth function v, which vanishes on Γ and find, after integration over Ω and using Green’s formula, that u(·,t)H01(Ω), 0<tT such that (5)(utt,v)+a(u;u,v)=(f,v),vH01(Ω),u(x,0)=u0(x),ut(x,0)=u1(x),xΩ, where (·,·) denotes the L2(Ω)-inner product and the bilinear form a(·;·,·) is defined by (6)a(w;u,v)=ΩA(w)u·vdx.

Henceforth, it will be assumed that the problem (5) has a unique solution u, and in the appropriate places to follow, additional conditions on the regularity of u which guarantee the convergence results, will be imposed.

Let 𝒯h be a quasiuniform triangulation of Ω with h=maxhK, where hK is the diameter of the triangle K𝒯h. With the triangulation 𝒯h, we associate the function space Sh consisting of continuous, piecewise linear functions on 𝒯h, vanishing on Γ; that is, (7)Sh={vC(Ω¯):vlinearinKforeachK𝒯h,v=0onΓ}. Using the above assumptions on 𝒯h, it is easy to see that Sh is a finite-dimensional subspace of the Hilbert space H01(Ω) .

Thus, the continuous-time finite element approximation is defined as to find a solution uh(t)Sh, 0<tT, such that (8)(uh,tt,vh)+a(uh;uh,vh)=(f,vh),vhSh,uh(0)=u0,uh,t(0)=u1, where uh,tt=2uh/t2. Since we have discretized only in the space variables, this is referred to as a spatially semidiscrete problem. The existence and uniqueness of the solution of (8) have been proved by Yuan and Wang .

3. Error Analysis for the Finite Element Method

To describe the error estimates for the finite element scheme (8), we will give some useful lemmas. In [17, 21] it was shown that the bilinear form a(·;·,·) is symmetric and positive definite and the following lemma was proved, which indicates that the bilinear form a(·;·,·) is continuous and coercive on Sh.

Lemma 1.

For h sufficiently small, there exist two positive constants C1,C2>0 such that, for all uh,vh,whSh, the coercive property (9)a(wh;uh,uh)C1uh12 and the boundedness property (10)|a(wh;uh,vh)|C2uh1vh1 hold true.

Lemma 2.

Let u~Sh be the standard Ritz projection such that (11)a(u(x,t);(u~-u)(x,t),vh)=0,vhSh. Thus u~ is the finite element approximation of the solution of the elliptic problem whose exact solution is u. From , we have (12)u-u~+hu-u~1Ch2u2,(13)(u-u~)t+h(u-u~)t1Ch2ut2, for some positive constant C independent of h and u.

And there exists a positive constant C0 independent of h, such that  (14)u~+u~tC0,fortT.

We now turn to describe the estimates for the finite element method. We give the error estimates in the H1-norm and L2-norm between the exact solution and the semidiscrete finite element solution.

Theorem 3.

Let u and uh be the solutions of problem (1) and the semidiscrete finite element scheme (8), respectively. Under the assumptions given in Section 1, if uh(0)=u~0 and uh,t(0)=u~1, for 0<tT, one has (15)u(t)-uh(t)+hu(t)-uh(t)1𝒞h2,(u(t)-uh(t))t𝒞h2, where 𝒞=C(uL2(H2),uL(H2),utL2(H2),uttL2(H2)) independent of h.

Proof.

For convenience, let u-uh=(u-u~)+(u~-uh)=:η+ξ. Then from (1), (8), and (11), we get the following error equation: (16)(ξtt,vh)+a(uh;ξ,vh)=-(ηtt,vh)+a(uh;u,vh)-a(u;u,vh),vhSh. Choosing vh=ξt in (16) and by (11), we get (17)(ξtt,ξt)+a(uh;ξ,ξt)=-(ηtt,ξt)+a(uh;u~,ξt)-a(u;u~,ξt). For the terms of (17), we have (18)(ξtt,ξt)=12ddt(ξt,ξt)=12ddtξt2.(19)a(uh;ξ,ξt)=ΩA(uh)ξ·ξtdx=12ddta(uh;ξ,ξ)-12ΩA(uh)uuhtξ·ξdx.(20)a(uh;u~,ξt)-a(u;u~,ξt)=Ω[A(uh)-A(u)]u~·ξtdx=ddt(Ω[A(uh)-A(u)]u~·ξdx)-Ω[A(uh)-A(u)]u~t·ξdx-Ω(A(uh)uuht-A(u)uut)u~·ξdx. Integrating (17) from 0 to t, combining with (18)–(20), and noting that ξ(0)=0 and ξt(0)=0, we have (21)12ξt2+12a(uh;ξ,ξ)=-0t(ηtt,ξt)dt+120t(ΩA(uh)uuhtξ·ξdx)dt+Ω[A(uh)-A(u)]u~·ξdx-0t(Ω[A(uh)-A(u)]u~t·ξdx)dt-0t(Ω(A(uh)uuht-A(u)uut)u~·ξdx)dti=15Ti. Now let us estimate the right-hand side terms of (21); for T1, there is (22)|T1|C0tηttξtdtC0t(ηtt2+ξt2)dt. For T2, by (2), we obtain (23)|T2|C0t|A(uh)uuht|ξ2dtCL0t|uht|ξ2dtC0tξ12dt, where we used the fact that |uh/t| is bounded by a positive constant .

For T3 by (14), Schwarz inequality, and (3), we get (24)|T3|Cu~A(uh)-A(u)ξCu~Lξ+ηξ1C(η2+ξ2)+ϵξ12C(η2+0tξt2dt)+ϵξ12, with ϵ being a small positive constant. For T4, similarly we have (25)|T4|C0tu~tA(uh)-A(u)ξdtC0tu~tLξ+ηξ1dtC0t(ξ2+η2+ξ12)dtC0t(η2+ξ12)dt.

For T5, by Lemma 2, we obtain (26)|T5|C0tu~uht-utξdtC0t(ξt+ηt)ξ1dtC0t(ξt2+ηt2+ξ12)dt. By Lemma 1, from (21)–(26), we get (27)ξt2+C0ξ12C1[0t(ηtt2+ηt2+η2)dt+η2]+C20t(ξt2+ξ12)dt+ϵξ12. Choosing proper ϵ and kicking the last term into the left-hand side of (27), and applying Gronwall’s lemma, for tT, we have (28)ξt2+ξ12C1[0T(ηtt2+ηt2+η2)dt+η2]Ch4[0T(utt22+ut22+u22)dt+u22]. Together with (12) and (13), this yields (15).

4. Two-Grid Finite Element Method

In this section, we will present a two-grid finite element algorithm for problem (1) based on two different finite element spaces. The idea of the two-grid method is to reduce the nonlinear problem on a fine grid into a linear system on the fine grid by solving a nonlinear problem on a coarse grid. The basic mechanisms are two quasiuniform triangulations of Ω, 𝒯H and 𝒯h, with two different mesh sizes H and h (H>h), and the corresponding piecewise linear finite element spaces SH and Sh which will be called the coarse-grid and the fine-grid spaces, respectively.

To solve problem (1), we introduce two-grid algorithms into finite element method. This method involves a nonlinear solution on the coarse grid space and a linear solution on the fine grid space. We present the two-grid finite element method with two steps.

Algorithm 4.

Consider the following.

Step  1. On the coarse grid 𝒯H, find uHSH, such that (29)(uH,tt,vH)+a(uH;uH,vH)=(f,vH),vHSH,uH(0)=u~0,uH,t(0)=u~1.Step  2. On the fine grid 𝒯h, find uhSh, such that (30)(uh,tt,vh)+a(uH;uh,vh)=(f,vh),  vhSh,uh(0)=u~0,uh,t(0)=u~1.

Now we consider the error estimates in the H1-norm for the two-grid finite element method Algorithm 4.

Theorem 5.

Let u and uh be the solutions of problem (1) and the two-grid scheme Algorithm 4, respectively. Under the assumptions given in Section 1, if uh(0)=u~0 and uh,t(0)=u~1, for 0<tT, we have (31)u(t)-uh(t)1𝒞(h+H2), where 𝒞=C(uL2(H2),uL(H2),uL(W1,),utL2(H2),utL(W1,),uttL2(H2)) independent of h.

Proof.

Once again, we set u-uh=(u-u~)+(u~-uh)=:η+ξ and choose vh=ξt. Then for Algorithm 4, we get the error equation (32)(ξtt,ξt)+a(uH;ξ,ξt)=-(ηtt,ξt)+a(uH;u~,ξt)-a(u;u~,ξt). Similarly as the proof of Theorem 3, we get (33)12ξt2+12a(uH;ξ,ξ)=-0t(ηtt,ξt)dt+120t(ΩA(uH)uuHtξ·ξdx)dt+Ω[A(uH)-A(u)]u~·ξdx-0t(Ω[A(uH)-A(u)]u~t·ξdx)dt-0t(Ω(A(uH)uuHt-A(u)uut)u~·ξdx)dti=15Ti. For T1 and T2, we can estimate them similarly as in Theorem 3. So our main task is to deal with T3T5. By (3), we have (34)|T3|Cu~A(uH)-A(u)ξCu~LuH-uξ1CuH-u2+ϵξ12, with ϵ being a small positive constant (35)|T4|C0tu~tA(uH)-A(u)ξdtC0tu~tLuH-uξ1dtC0t(uH-u2+ξ12)dt,(36)|T5|C0tu~uHt-utξdtC0t(uH-u)tξ1dtC0t((uH-u)t2+ξ12)dt. Substituting the estimates of Ti in (33) and by Lemma 1, we obtain (37)ξt2+C0ξ12C10t(ηtt2+uH-u2+(uH-u)t2)dt+C2uH-u2+C30t(ξt2+ξ12)dt+ϵξ12. Choosing proper ϵ and kicking the last term into the left-hand side of (33), and applying Gronwall’s lemma, for tT, we have (38)ξt2+ξ12C10t(ηtt2+uH-u2+(uH-u)t2)dt+C2uH-u2. By Theorem 3, we obtain (39)ξt2+ξ12𝒞(h4+H4), where 𝒞=C(uL2(H2),uL(H2),uL(W1,),utL2(H2),utL(W1,),uttL2(H2)) independent of h. Thus, (40)ξt+ξ1𝒞(h2+H2). By (12) and the triangular inequality, the proof is complete.

Remark 6.

In order to give the fully discrete scheme, we further discretize time t of the semidiscrete two-grid finite element method in this section. We consider a time step Δt and approximate the solutions at tn=nΔt, Δt=T/N, n=0,1,,N. Denote uhn=uh(tn), uh,ttn=(uhn+1-2uhn+uhn-1)/(Δt)2, uh,tn=(uhn+1-uhn)/Δt, we can get the fully discrete two-grid finite element scheme for (1). For simplicity and convenience, we only give the fully discrete scheme for Algorithm 4.

Algorithm  4 . Consider the following.

Step  1. On the coarse grid 𝒯H, find uHnSH(n=1,2,), such that (41)(uH,ttn,vH)+a(uHn+1;uHn+1,vH)=(fn+1,vH),vHSH,uH0=u~0,uH,t0=u~1.

Step  2. On the fine grid 𝒯h, find uhnSh(n=1,2,), such that (42)(uh,ttn,vh)+a(uHn+1;uhn+1,vh)=(fn+1,vh),vhSh,uh0=u~0,uh,t0=u~1.

We can get the same kind of estimate as Theorem 5 with the result un-uhn1𝒞((Δt)2+h+H2).

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (no. 11301456) and Shandong Province Natural Science Foundation (nos. ZR2010AQ010 and ZR2011AQ021).

Xu J. Two grid finite element discretizations for linear and nonlinear elliptic equations 1992 AM105 University Park, Pa, USA Department of Mathematics, Pennsylvania State University Xu J. Two-grid discretization techniques for linear and nonlinear PDEs SIAM Journal on Numerical Analysis 1996 33 5 1759 1777 10.1137/S0036142992232949 MR1411848 Dawson C. N. Wheeler M. F. Two-grid methods for mixed finite element approximations of nonlinear parabolic equations Contemporary Mathematics 1994 180 191 203 Dawson C. N. Wheeler M. F. Woodward C. S. A two-grid finite difference scheme for nonlinear parabolic equations SIAM Journal on Numerical Analysis 1998 35 2 435 452 10.1137/S0036142995293493 MR1618822 Chen C. Liu W. A two-grid method for finite element solutions of nonlinear parabolic equations Abstract and Applied Analysis 2012 2012 11 10.1155/2012/391918 391918 MR2975310 Wu L. Allen M. B. Two-grid methods for mixed finite-element solution of coupled reaction-diffusion systems Numerical Methods for Partial Differential Equations 1999 15 5 589 604 10.1002/(SICI)1098-2426(199909)15:5&lt;589::AID-NUM6&gt;3.0.CO;2-W MR1708087 Chen Y. Huang Y. Yu D. A two-grid method for expanded mixed finite-element solution of semilinear reaction-diffusion equations International Journal for Numerical Methods in Engineering 2003 57 2 193 209 10.1002/nme.668 MR1970318 Chen Y. Liu H.-W. Liu S. Analysis of two-grid methods for reaction-diffusion equations by expanded mixed finite element methods International Journal for Numerical Methods in Engineering 2007 69 2 408 422 10.1002/nme.1775 MR2284576 Chen Y. Luan P. Lu Z. Analysis of two-grid methods for nonlinear parabolic equations by expanded mixed finite element methods Advances in Applied Mathematics and Mechanics 2009 1 6 830 844 MR2570473 Chen L. Chen Y. Two-grid method for nonlinear reaction-diffusion equations by mixed finite element methods Journal of Scientific Computing 2011 49 3 383 401 10.1007/s10915-011-9469-3 MR2853157 Bi C. Ginting V. Two-grid finite volume element method for linear and nonlinear elliptic problems Numerische Mathematik 2007 108 2 177 198 10.1007/s00211-007-0115-9 MR2358002 Chen C. Yang M. Bi C. Two-grid methods for finite volume element approximations of nonlinear parabolic equations Journal of Computational and Applied Mathematics 2009 228 1 123 132 10.1016/j.cam.2008.09.001 MR2514269 Chen C. Liu W. Two-grid finite volume element methods for semilinear parabolic problems Applied Numerical Mathematics 2010 60 1-2 10 18 10.1016/j.apnum.2009.08.004 MR2566074 Chen C. Liu W. A two-grid method for finite volume element approximations of second-order nonlinear hyperbolic equations Journal of Computational and Applied Mathematics 2010 233 11 2975 2984 10.1016/j.cam.2009.11.043 MR2592271 Dupont T. L 2 -estimates for Galerkin methods for second order hyperbolic equations SIAM Journal on Numerical Analysis 1973 10 880 889 MR0349045 10.1137/0710073 Baker G. A. Error estimates for finite element methods for second order hyperbolic equations SIAM Journal on Numerical Analysis 1976 13 4 564 576 MR0423836 10.1137/0713048 Yuan Y. R. Wang H. Error estimates for finite element solutions of nonlinear hyperbolic equations Journal of Systems Science and Mathematical Sciences 1985 5 3 161 171 MR847810 Wang H. Stability and convergence of the finite element method for a class of second order nonlinear hyperbolic equations Mathematica Numerica Sinica 1986 8 150 158 Kumar S. Nataraj N. Pani A. K. Finite volume element method for second order hyperbolic equations International Journal of Numerical Analysis and Modeling 2008 5 1 132 151 MR2479124 Larsson S. Thomée V. Partial Differential Equations with Numerical Methods 2003 45 Berlin, Germany Springer MR1995838 Thomée V. Galerkin Finite Element Methods for Parabolic Problems 2006 25 2nd Berlin, Germany Springer MR2249024 Wheeler M. F. A priori L2 error estimates for Galerkin approximations to parabolic partial differential equations SIAM Journal on Numerical Analysis 1973 10 723 759 MR0351124 10.1137/0710062 Hayes L. J. A modified backward time discretization for nonlinear parabolic equations using patch approximations SIAM Journal on Numerical Analysis 1981 18 5 781 793 10.1137/0718052 MR629663