Numerical inversions of a source term in the FADE with a Dirichlet boundary condition using final observations

Abstract

This paper deals with an inverse problem of determining a source term in the one-dimensional fractional advection–dispersion equation (FADE) with a Dirichlet boundary condition on a finite domain, using final observations. On the basis of the shifted Grünwald formula, a finite difference scheme for the forward problem of the FADE is given, by means of which the source magnitude depending upon the space variable is reconstructed numerically by applying an optimal perturbation regularization algorithm. Numerical inversions with noisy data are carried out for the unknowns taking three functional forms: polynomials, trigonometric functions and index functions. The reconstruction results show that the inversion algorithm is efficient for the inverse problem of determining source terms in a FADE, and the algorithm is also stable for additional data having random noises.