Isoparametric multigrid method for reaction-diffusion equations on two-dimensional domains (Q850188)

The authors consider the finite element discretization of reaction-diffusion problems on 2D curved domains with a Lipschitz-continuous boundary. The curved boundary is approximated by using isoparametric finite elements. For applying multigrid methods a sequence of finite element triangulations is constructed, which gives a sequence of domains \(\Omega _k\), \(k=0,1,\dots,l,\) approximating the domain \(\Omega \). The discretization is constructed by means of higher order finite elements on a sequence of corresponding nonnested triangulations. The authors define an intergrid transfer operator, which is necessary for a multigrid nonnested algorithm and analyze its stability. An optimal convergence rate for the two-level method with the Jacobi-like smoother is obtained with respect to the energy norm. An upper bound for the number of iterations of the coarse-level residual problem is found such that the cost of the computational work of one iteration step of the multigrid algorithm to be proportional to the number of unknowns.