A novel unstructured mesh finite element method for solving the time-space fractional wave equation on a two-dimensional irregular convex domain

摘要

Most existing research on applying the finite element method to discretize space fractional operators is studied on regular domains using either uniform structured triangular meshes, or quadrilateral meshes. Since many practical problems involve irregular convex domains, such as the human brain or heart, which are difficult to partition well with a structured mesh, the existing finite element method using the structured mesh is less efficient. Research on the finite element method using a completely unstructured mesh on an irregular domain is of great significance. In this paper, a novel unstructured mesh finite element method is developed for solving the time-space fractional wave equation on a two-dimensional irregular convex domain. The novel unstructured mesh Galerkin finite element method is used to discretize in space and the Crank-Nicolson scheme is used to discretize the Caputo time fractional derivative. The implementation of the unstructured mesh Crank-Nicolson Galerkin method (CNGM) is detailed and the stability and convergence of the numerical scheme are analysed. Numerical examples are presented to verify the theoretical analysis. To highlight the ability of the proposed unstructured mesh Galerkin finite element method, a comparison of the unstructured mesh with the structured mesh in the implementation of the numerical scheme is conducted. The proposed numerical method using an unstructured mesh is shown to be more effective and feasible for practical applications involving irregular convex domains.