Numerical inversions for diffusion coefficients in two-dimensional space fractional diffusion equation

This article deals with an inverse problem of determining the diffusion coefficients in 2D fractional diffusion equation with a Dirichlet boundary condition by the final observations at the final time. The forward problem is solved by the alternating direction implicit finite-difference scheme with the discrete of fractional derivative by shift Grunwald formula and a numerical text which is to prove its numerically stability and convergence is given. Furthermore, the homotopy regularization algorithm with the regularization parameter chosen by a Sigmoid-type function is introduced to solve the inversion problem numerically. Numerical inversions both with accurate data and noisy data are carried out for the unknown diffusion coefficients of constant and variable with polynomials, trigonometric and index functions. The reconstruction results show that the inversion algorithm is efficient for the inverse problem of determining diffusion coefficients in 2D space fractional diffusion equation, and the algorithm is also numerically stable for additional date having random noises.