Recovering an unknown source in a fractional diffusion problem

A standard inverse problem is to determine a source which is supported in an unknown domain D from external boundary measurements. Here we consider the case of a time-independent situation where the source is equal to unity in an unknown subdomain D of a larger given domain Omega and the boundary of D has the star-like shape, i.e. partial derivative D = {q(theta)(cos theta, sin theta)(T) : theta is an element of [0, 2 pi]}. Overposed measurements consist of time traces of the solution or its flux values on a set of discrete points on the boundary partial derivative Omega. The case of a parabolic equation was considered in [6]. In our situation we extend this to cover the subdiffusion case based on an anomalous diffusion model and leading to a fractional order differential operator. We will show a uniqueness result and examine a reconstruction algorithm. One of the main motives for this work is to examine the dependence of the reconstructions on the parameter alpha, the exponent of the fractional operator which controls the degree of anomalous behaviour of the process. Some previous inverse problems based on fractional diffusion models have shown considerable differences between classical Brownian diffusion and the anomalous case. (C) 2018 Elsevier Inc. All rights reserved.