| Peer-Reviewed

Two-scale Finite Element Discretizations for Semilinear Parabolic Equations

Received: 5 February 2020    Accepted: 25 September 2020    Published: 16 November 2020
Views:       Downloads:
Abstract

In this paper, to reduce the computational cost of solving semilinear parabolic equations on a tensor product domain Ω⊂ℝd with d = 2 or 3, some two-scale finite element discretizations are proposed and analyzed. The time derivative in semilinear parabolic equations is approximated by the backward Euler finite difference scheme. The two-scale finite element method is designed for the space discretization. The idea of the two-scale finite element method is based on an understanding of a finite element solution to an elliptic problem on a tensor product domain. The high frequency parts of the finite element solution can be well captured on some univariate fine grids and the low frequency parts can be approximated on a coarse grid. Thus the two-scale finite element approximation is defined as a linear combination of some standard finite element approximations on some univariate fine grids and a coarse grid satisfying H = O (h1/2), where h and H are the fine and coarse mesh widths, respectively. It is shown theoretically and numerically that the backward Euler two-scale finite element solution not only achieves the same order of accuracy in the H1 (Ω) norm as the backward Euler standard finite element solution, but also reduces the number of degrees of freedom from O(h-d×τ-1) to O(h-((d)+1)/2×τ-1) where τ is the time step. Consequently the backward Euler two-scale finite element method for semilinear parabolic equations is more efficient than the backward Euler standard finite element method.

Published in Applied and Computational Mathematics (Volume 9, Issue 6)
DOI 10.11648/j.acm.20200906.12
Page(s) 179-186
Creative Commons

This is an Open Access article, distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution and reproduction in any medium or format, provided the original work is properly cited.

Copyright

Copyright © The Author(s), 2024. Published by Science Publishing Group

Keywords

Two-scale, Finite Element, Combination, Semilinear Parabolic Equation

References
[1] Bank R E. Hierarchical bases and the finite element method. Acta Numer, Cambridge University Press, Cambridge, UK, 1996, 5: 1-43.
[2] Blum H, Lin Q, Rannacher R. Asymptotic error expansion and Richardson extrapolation for linear finite elements. Numer Math, 1986, 49: 11-37.
[3] Brenner S C, Scott L R. The Mathematical Theory of Finite Element Methods, 3rd ed.. Texts in Applied Mathematics, Springer, New York, 2008, 15.
[4] Bungartz H-J, Griebel M. Sparse grids. Acta Numer, 2004, 13: 1-123.
[5] Bungartz H-J, Griebel M, Rüde U. Extrapolation, combination, and sparse grid techniques for elliptic boundary value problems. Comput Methods Appl Mech Engrg, 1994, 116: 243-252.
[6] Chen C, Yang M, Bi C. Two-grid methods for finite volume element approximations of nonlinear parabolic equations. J Comput Appl Math, 2009, 228: 123-132.
[7] Chen C, Liu W. Two-grid finite volume element methods for semilinear parabolic problems. Appl Numer Math, 2010, 60 (1): 10-18.
[8] Chen Y, Chen L, Zhang X. Two-grid method for nonlinear parabolic equations by expanded mixed finite element methods. Numer Methods Partial Differential Equations, 2013, 29 (4): 1238-1256.
[9] Ciarlet P G. The Finite Element Method for Elliptic Problems. Studies in Mathematics and Its Applications, North-Holland Publishing Co., Amsterdam, 1978, 4.
[10] Dawson C, Wheeler M. Two-grid methods for mixed finite element approximations of nonlinear parabolic equations. Contemp Math, 1994, 180: 191-203.
[11] Dawson C, Wheeler M, Woodward C. A two-grid finite difference scheme for nonlinear parabolic equations. SIAM J Numer Anal, 1998, 35 (2): 435-452.
[12] Delvos F. d-Variate Boolean interpolation. J Approx Theory, 1982, 34: 99-114.
[13] Gao X, Liu F, Zhou A. Three-scale finite element eigenvalue discretizations. BIT, 2008, 48: 533-562.
[14] Garcke J, Griebel M. On the computation of the eigenproblems of hydrogen and helium in strong magnetic and electric fields with the sparse grid combination technique. J Comput Phys, 2000, 165: 694-716.
[15] Griebel M, Schneider M, Zenger C. A combination technique for the solution of sparse grid problem. Iterative Methods in Linear Algebra (P. de Groen and P. Beauwens, eds). IMACS, North Holland, Amsterdam, 1992, 263-281.
[16] Hegland M. Adaptive sparse grids. ANZIAM J, 2003, 44: C335-C353.
[17] Hegland M, Garcke J, Challis V. The combination technique and some generalisations. Linear Algebra Appl, 2007, 420: 249-275.
[18] Hennart J P, Mund E H. On the h- and p-versions of the extrapolated Gordon's projector with applications to elliptic equations. SIAM J Sci Statist Comput, 1988, 9: 773-791.
[19] Jin J, Wei N, Zhang H. A two-grid finite-element method for the nonlinear Schrödinger equation. J Comput Math, 2015, 33 (2): 146-157.
[20] Lasis A, Süli E. hp-version discontinuous Galerkin finite element method for semilinear parabolic problems. SIAM J Numer Anal, 2007, 45 (4): 1544-1569.
[21] Liao X, Zhou A. A multi-parameter splitting extrapolation and a parallel algorithm for elliptic eigenvalue problem. J Comput Math, 1998, 16: 213-220.
[22] Liem C B, Lü T, Shih T M. The Splitting Extrapolation Method. Series on Applied Mathematics, World Scientific Publishing Co., River Edge, NJ, 1995, 7.
[23] Lin Q, Lin J. Finite Element Methods: Accuracy and Improvement. Science Press, Beijing, 2006.
[24] Lin Q, Lü T. Asymptotic expansions for finite element eigenvalues and finite element solution. Extrapolation procedures in the finite element method, Bonner Math Schriften, Univ Bonn, Bonn, 1984, 158: 1-10.
[25] Lin Q, Yan N, Zhou A. A sparse finite element method with high accuracy. Part I. Numer Math, 2001, 88: 731-742.
[26] Lin Q, Zhu Q. Undirectional extrapolations of finite difference and finite elements. J Engrg Math, 1984, 1: 1-12 (in Chinese).
[27] Lin Q, Zhu Q. The Preprocessing and Postprocessing for the Finite Element Method. Shanghai Sci Tech Press, 1994 (in Chinese).
[28] Liu W, Rui H, Bao Y. Two kinds of two-grid algorithms for finite difference solutions of semilinear parabolic equations. J Sys Sci Math Sci, 2010, 30 (2): 181-190.
[29] Liu F, Stynes M, Zhou A. Postprocessed two-scale finite element discretizations, part I. SIAM J Numer Anal, 2011, 49: 1947-1971.
[30] Liu F, Zhou A. Two-scale finite element discretizations for partial differential equations. J Comput Math, 2006, 24: 373-392.
[31] Liu F, Zhou A. Localizations and parallelizations for two-scale finite element discretizations. Commun Pure Appl Anal, 2007, 6: 757-773.
[32] Liu F, Zhou A. Two-scale Boolean Galerkin discretizations for Fredholm integral equations of the second kind. SIAM J Numer Anal, 2007, 45: 296-312.
[33] Pflaum C. Convergence of the combination technique for second-order elliptic differential equations. SIAM J Numer Anal, 1997, 34: 2431-2455.
[34] Pflaum C, Zhou A. Error analysis of the combination technique. Numer Math, 1999, 84: 327-350.
[35] Shi D, Yang H. Unconditional optimal error estimates of a two-grid method for semilinear parabolic equations. Appl Math Comput, 2017, 310: 40-47.
[36] Shi D, Mu P, Yang H. Superconvergence analysis of a two-grid method for semilinear parabolic equations. Appl Math Lett, 2018, 84: 34-41.
[37] Thomée V. Galerkin Finite Element Methods for Parabolic Problems. Springer-Verlag, Berlin, 1984.
[38] Thomée V, Xu J, Zhang N. Superconvergence of the gradient in piecewise linear finite-element approximation to a parabolic problem. SIAM J Numer Anal, 1989, 26 (3): 553-573.
[39] Xu J. Two-grid discretization techniques for linear and nonlinear PDEs. SIAM J Numer Anal, 1996, 33: 1759-1777.
[40] Xu J, Zhou A. Local and parallel finite element algorithms based on two-grid discretizations. Math Comp, 2000, 69: 881-909.
[41] Xu, J, Zhou, A. A two-grid discretization scheme for eigenvalue problems. Math Comput, 2001, 70: 17-25.
[42] Yang J, Xing X. A two-grid discontinuous Galerkin method for a kind of nonlinear parabolic problems. Appl Math Comput, 2019, 346: 96-108.
[43] Yserentant H. On the multi-level splitting of finite element spaces. Numer Math, 1986, 49: 379-412.
[44] Yserentant, H. Old and new convergence proofs for multigrid methods. Acta Numerica, Cambridge University Press, Cambridge, UK, 1993, 285-326.
[45] Zenger, C. Sparse grids. Parallel Algorithms for Partial Differential Equations (Kiel, 1990), Notes Numer Fluid Mech, Vieweg, Braunschweig, 1991, 31: 241-251.
[46] Zhang H, Jin J, Wang J. Two-grid finite-element method for the two-dimensional time-dependent Schrödinger equation. Adv Appl Math Mech, 2013, 5 (2): 180-193.
[47] Zhou A, Li J. The full approximation accuracy for the stream function-vorticity-pressure method. Numer Math, 1994, 68: 427-435.
[48] Zhou A, Liem C, Shih T, Lü T. A multi-parameter splitting extrapolation and a parallel algorithm. Systems Sci Math Sci, 1997, 10: 253-260.
Cite This Article
  • APA Style

    Fang Liu. (2020). Two-scale Finite Element Discretizations for Semilinear Parabolic Equations. Applied and Computational Mathematics, 9(6), 179-186. https://doi.org/10.11648/j.acm.20200906.12

    Copy | Download

    ACS Style

    Fang Liu. Two-scale Finite Element Discretizations for Semilinear Parabolic Equations. Appl. Comput. Math. 2020, 9(6), 179-186. doi: 10.11648/j.acm.20200906.12

    Copy | Download

    AMA Style

    Fang Liu. Two-scale Finite Element Discretizations for Semilinear Parabolic Equations. Appl Comput Math. 2020;9(6):179-186. doi: 10.11648/j.acm.20200906.12

    Copy | Download

  • @article{10.11648/j.acm.20200906.12,
      author = {Fang Liu},
      title = {Two-scale Finite Element Discretizations for Semilinear Parabolic Equations},
      journal = {Applied and Computational Mathematics},
      volume = {9},
      number = {6},
      pages = {179-186},
      doi = {10.11648/j.acm.20200906.12},
      url = {https://doi.org/10.11648/j.acm.20200906.12},
      eprint = {https://article.sciencepublishinggroup.com/pdf/10.11648.j.acm.20200906.12},
      abstract = {In this paper, to reduce the computational cost of solving semilinear parabolic equations on a tensor product domain Ω⊂ℝd with d = 2 or 3, some two-scale finite element discretizations are proposed and analyzed. The time derivative in semilinear parabolic equations is approximated by the backward Euler finite difference scheme. The two-scale finite element method is designed for the space discretization. The idea of the two-scale finite element method is based on an understanding of a finite element solution to an elliptic problem on a tensor product domain. The high frequency parts of the finite element solution can be well captured on some univariate fine grids and the low frequency parts can be approximated on a coarse grid. Thus the two-scale finite element approximation is defined as a linear combination of some standard finite element approximations on some univariate fine grids and a coarse grid satisfying H = O (h1/2), where h and H are the fine and coarse mesh widths, respectively. It is shown theoretically and numerically that the backward Euler two-scale finite element solution not only achieves the same order of accuracy in the H1 (Ω) norm as the backward Euler standard finite element solution, but also reduces the number of degrees of freedom from O(h-d×τ-1) to O(h-((d)+1)/2×τ-1) where τ is the time step. Consequently the backward Euler two-scale finite element method for semilinear parabolic equations is more efficient than the backward Euler standard finite element method.},
     year = {2020}
    }
    

    Copy | Download

  • TY  - JOUR
    T1  - Two-scale Finite Element Discretizations for Semilinear Parabolic Equations
    AU  - Fang Liu
    Y1  - 2020/11/16
    PY  - 2020
    N1  - https://doi.org/10.11648/j.acm.20200906.12
    DO  - 10.11648/j.acm.20200906.12
    T2  - Applied and Computational Mathematics
    JF  - Applied and Computational Mathematics
    JO  - Applied and Computational Mathematics
    SP  - 179
    EP  - 186
    PB  - Science Publishing Group
    SN  - 2328-5613
    UR  - https://doi.org/10.11648/j.acm.20200906.12
    AB  - In this paper, to reduce the computational cost of solving semilinear parabolic equations on a tensor product domain Ω⊂ℝd with d = 2 or 3, some two-scale finite element discretizations are proposed and analyzed. The time derivative in semilinear parabolic equations is approximated by the backward Euler finite difference scheme. The two-scale finite element method is designed for the space discretization. The idea of the two-scale finite element method is based on an understanding of a finite element solution to an elliptic problem on a tensor product domain. The high frequency parts of the finite element solution can be well captured on some univariate fine grids and the low frequency parts can be approximated on a coarse grid. Thus the two-scale finite element approximation is defined as a linear combination of some standard finite element approximations on some univariate fine grids and a coarse grid satisfying H = O (h1/2), where h and H are the fine and coarse mesh widths, respectively. It is shown theoretically and numerically that the backward Euler two-scale finite element solution not only achieves the same order of accuracy in the H1 (Ω) norm as the backward Euler standard finite element solution, but also reduces the number of degrees of freedom from O(h-d×τ-1) to O(h-((d)+1)/2×τ-1) where τ is the time step. Consequently the backward Euler two-scale finite element method for semilinear parabolic equations is more efficient than the backward Euler standard finite element method.
    VL  - 9
    IS  - 6
    ER  - 

    Copy | Download

Author Information
  • School of Statistics and Mathematics, Central University of Finance and Economics, Beijing, China

  • Sections