In this paper,we develop a two-grid method(TGM)based on the FEM for 2D nonlinear time fractional two-term mixed sub-diffusion and diffusion wave equations.A two-grid algorithm is proposed for solving the nonlinear sys...In this paper,we develop a two-grid method(TGM)based on the FEM for 2D nonlinear time fractional two-term mixed sub-diffusion and diffusion wave equations.A two-grid algorithm is proposed for solving the nonlinear system,which consists of two steps:a nonlinear FE system is solved on a coarse grid,then the linearized FE system is solved on the fine grid by Newton iteration based on the coarse solution.The fully discrete numerical approximation is analyzed,where the Galerkin finite element method for the space derivatives and the finite difference scheme for the time Caputo derivative with orderα∈(1,2)andα1∈(0,1).Numerical stability and optimal error estimate O(h^(r+1)+H^(2r+2)+τ^(min{3−α,2−α1}))in L^(2)-norm are presented for two-grid scheme,where t,H and h are the time step size,coarse grid mesh size and fine grid mesh size,respectively.Finally,numerical experiments are provided to confirm our theoretical results and effectiveness of the proposed algorithm.展开更多
基金This work is supported by the State Key Program of National Natural Science Foundation of China(11931003)National Natural Science Foundation of China(41974133,11971410)+2 种基金Project for Hunan National Applied Mathematics Center of Hunan Provincial Science and Technology Department(2020ZYT003)Hunan Provincial Innovation Foundation for Postgraduate,China(XDCX2020B082,XDCX2021B098)Postgraduate Scientific Research Innovation Project of Hunan Province(CX20210597).
文摘In this paper,we develop a two-grid method(TGM)based on the FEM for 2D nonlinear time fractional two-term mixed sub-diffusion and diffusion wave equations.A two-grid algorithm is proposed for solving the nonlinear system,which consists of two steps:a nonlinear FE system is solved on a coarse grid,then the linearized FE system is solved on the fine grid by Newton iteration based on the coarse solution.The fully discrete numerical approximation is analyzed,where the Galerkin finite element method for the space derivatives and the finite difference scheme for the time Caputo derivative with orderα∈(1,2)andα1∈(0,1).Numerical stability and optimal error estimate O(h^(r+1)+H^(2r+2)+τ^(min{3−α,2−α1}))in L^(2)-norm are presented for two-grid scheme,where t,H and h are the time step size,coarse grid mesh size and fine grid mesh size,respectively.Finally,numerical experiments are provided to confirm our theoretical results and effectiveness of the proposed algorithm.
