Abstract
In this paper, we design and analyze the iterative two-grid methods for the discontinuous Galerkin (DG) discretization of semilinear elliptic partial differential equations (PDEs). We first present an iterative two-grid method that is just like the classical iterative two-grid methods for nonsymmetric or indefinite linear elliptic PDEs, namely, to solve a semilinear problem on the coarse space and then to solve a symmetric positive definite problem on the fine space. Secondly, we designed another iterative two-grid method, which replace the semilinear term by using the corresponding first-order Taylor expansion. Specifically, we need to construct a suitable initial value, which can be sorted out from an auxiliary variational problem, for the second iterative method. We also provide the error estimates for the second iterative algorithm and present numerical experiments to illustrate the theoretical result.
Similar content being viewed by others
References
Arnold, D.N.: An interior penalty finite element method with discontinuous elements. SIAM J. Numer. Anal. 19(4), 742–760 (1982)
Arnold, D., Brezzi, F., Cockburn, B., Marini, L.: Unified analysis of discontinuous Galerkin methods for elliptic problems. SIAM J. Numer. Anal. 39(5), 1749–1779 (2002)
Bi, C.J., Ginting, V.: Two-grid discontinuous Galerkin method for quasi-linear elliptic problems. J. Sci. Comput. 49(3), 311–331 (2011)
Bi, C.J., Wang, C., Lin, Y.P.: Pointwise error estimates and two-grid algorithms of discontinuous Galerkin method for strongly nonlinear elliptic problems. J. Sci. Comput. 67(1), 153–175 (2016)
Bi, C.J., Wang, C., Lin, Y.P.: A posteriori error estimates of two-grid finite element methods for nonlinear elliptic problems. J. Sci. Comput. 74(1), 23–48 (2018)
Brenner, S.C., Cui, J., Gudi, T., Sung, L.-Y.: Multigrid algorithms for symmetric discontinuous Galerkin methods on graded meshes. Numer. Math. 119(1), 21–47 (2011)
Congreve, S., Houston, P., Wihler, T.P.: Two-grid hp-version discontinuous Galerkin finite element methods for second-order quasilinear elliptic PDEs. J. Sci. Comput. 55(2), 471–497 (2013)
Logg, A., Mardal, K.-A., Wells, G.: Automated solution of differential equations by the finite element method. Springer (2012)
Scott, L.R., Zhang, S.: Finite element interpolation of nonsmooth functions satisfying boundary conditions. Math. Comp. 54(190), 483–493 (1990)
Xu, J.C.: A new class of iterative methods for nonselfadjoint or indefinite problems. SIAM J. Numer. Anal. 29(2), 303–319 (1992)
Xu, J.C.: A novel two-grid method for semilinear elliptic equations. SIAM J. Sci. Comput. 15(1), 231–237 (1994)
Xu, J.C.: Two-grid discretization techniques for linear and nonlinear PDEs. SIAM J. Numer. Anal. 33(5), 1759–1777 (1996)
You, X.: On existence and uniqueness of solutions of a class of second order semilinear elliptic boundary value problems. J. Sichuan Normal Univ. (Nat. Sci.) 26(5), 479–481 (2003)
Zhan, J.J., Zhong, L.Q., Peng, J.: Discontinuous Galerkin methods for semilinear elliptic boundary value problem. Adv. Appl. Math. Mech. 15(2), 450–467 (2023)
Zhang, W.F., Fan, R.H., Zhong, L.Q.: Iterative two-grid methods for semilinear elliptic equations. Comput. Math. Appl. 80(3), 522–530 (2020)
Zhong, L.Q., Zhou, L.L., Liu, C.M., Peng, J.: Two-grid IPDG discretization scheme for nonlinear elliptic PDEs. Commun. Nonlinear Sci. Numer. Simul. 95, 105587 (2021)
Funding
The second and third authors are supported by the National Natural Science Foundation of China (No. 12071160). The third author is also supported by National Natural Science Foundation of China (No. 12101250) and the Science and Technology Projects in Guangzhou (No. 202201010644).
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare no competing interests.
Additional information
Communicated by: Long Chen
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendix
Appendix
In this appendix, we present the proof of 4.24 in Lemma 4.5.
First, we define the projection operator \(P_H: H^1(\Omega )\rightarrow V_H\): for \(w \in H^1(\Omega )\),
Next, we demonstrate the error estimates of \(P_H\) as follows.
Lemma 2.1
For the projection operator \(P_H\) and the solution \(u \in H^{r+1} \cap H_0^1(\Omega )\) of 2.4, we have
Proof
For any \(v_H \in V_H\), using the coerciveness and continuity of \(a_h(\cdot , \cdot )\) and (6.1), we could obtain that
which implies that
Using the triangle inequality and (6.4), we have
Further using (6.5), the arbitrary of \(v_H\) and 4.6, we have
which derives (6.2).
For any fixed \(g \in L^2(\Omega )\), we introduce the following auxiliary problem: find \(w \in H^2(\Omega ) \cap H_0^1(\Omega )\) such that
and assume that the following regularity result holds
Using (6.6), similar technique in 2.8 with \(w \in H^2(\Omega ) \cap H_0^1(\Omega )\), (6.1), the continuity of \(a_h(\cdot , \cdot )\), 4.6 and (6.2), it is obtained that
At last, using (6.8), (6.7) and the arbitrary of g, we could obtain (6.3). \(\square \)
Assume \(u \in H^{r+1}(\Omega )\) is the solution of the problem 2.4, for any \(w_H \in V_H\), we define operator \(\Phi _H: V_H \rightarrow V_H\) by
where
By the coerciveness and continuity of \(a_h(\cdot , \cdot )\), and Assumption 2.1, we could prove that \(b_h(\cdot , \cdot )\) is continuous and coercive for \(V_H\), which implies the operator \(\Phi _H\) is well defined.
For the given solution u of the problem 2.4, we define a space
where the projection operator \(P_H\) is defined in (6.1), \(\delta _H = C_1 \Vert P_H u - u\Vert _0\) and \(C_1 \) is a constant which can be sufficiently large and does not depend on the mesh size.
Using highly similar technique in Section 2.3.4 of [14], we could prove that there is a fixed point \(\bar{w}_H \in B_H\) in (6.9) by Brouwer fixed point theorem, and here is omitted. Further, we can see that the fixed point \(\bar{w}_H\) is exactly the solution \(\tilde{u}_H\) of 3.3, which implies \(\tilde{u}_H \in B_H\). Then using the triangle inequality, (6.10) and (6.3), we have
which completes the proof of 4.24 in Lemma 4.5.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Zhan, J., Zhong, L. & Peng, J. Iterative two-grid methods for discontinuous Galerkin finite element approximations of semilinear elliptic problem. Adv Comput Math 49, 84 (2023). https://doi.org/10.1007/s10444-023-10082-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s10444-023-10082-8