Inverse Problems for Fractional Diffusion Equations

In recent decades, significant interest, based on physics and engineering applications, has developed on so-called anomalous diffusion processes that possess different spread functions with classical ones. The resulting differential equation whose fundamental solution matches this decay process is best modeled by an equation containing a fractional order derivative. This dissertation mainly focuses on some inverse problems for fractional diffusion equations.

After some background introductions and preliminaries in Section 1 and 2, in the third section we consider our first inverse boundary problem. This is where an unknown boundary condition is to be determined from overposed data in a time- fractional diffusion equation. Based upon the fundamental solution in free space, we derive a representation for the unknown parameters as the solution of a nonlinear Volterra integral equation of second kind with a weakly singular kernel. We are able to make physically reasonable assumptions on our constraining functions (initial and given boundary values) to be able to prove a uniqueness and reconstruction result. This is achieved by an iterative process and is an immediate result of applying a certain fixed point theorem. Numerical examples are presented to illustrate the validity and effectiveness of the proposed method.

In the fourth section a reaction-diffusion problem with an unknown nonlinear source function, which has to be determined from overposed data, is considered. A uniqueness result is proved and a numerical algorithm including convergence analysis under some physically reasonable assumptions is presented in the one-dimensional case. To show effectiveness of the proposed method, some results of numerical simulations are presented. In Section 5, we also attempted to reconstruct a nonlinear source in a heat equation from a number of known input sources. This represents a new research even for the case of classical diffusion and would be the first step in a solution method for the fractional diffusion case. While analytic work is still in progress on this problem, Newton and Quasi-Newton method are applied to show the feasibility of numerical reconstructions.

In conclusion, the fractional diffusion equations have some different properties with the classical ones but there are some similarities between them. The classical tools like integral equations and fixed point theory still hold under slightly different assumptions. Inverse problems for fractional diffusion equations have applications in many engineering and physics areas such as material design, porous media. They are trickier than classical ones but there are also some advantages due to the mildly ill-conditioned singularity caused by the new kernel functions.