Simultaneous inversion for the fractional exponents in the space-time fractional diffusion equation ∂t β u = −(− Δ) α/2 u − (− Δ) γ/2 u

Abstract

Abstract In this article, we consider the space-time fractional (nonlocal) equation characterizing the so-called “double-scale” anomalous diffusion ∂ t β u ( t , x ) = − ( − Δ ) α / 2 u ( t , x ) − ( − Δ ) γ / 2 u ( t , x ) , t > 0 , − 1 < x < 1 , $$\begin{array}{} \displaystyle \partial_t^\beta u(t, x) = -(-\Delta)^{\alpha/2}u(t,x) - (-\Delta)^{\gamma/2}u(t,x), \, \, t \gt 0, \, -1 \lt x \lt 1, \end{array}$$ where ∂ t β $\begin{array}{} \displaystyle \partial_t^\beta \end{array}$ is the Caputo fractional derivative of order β ∈ (0, 1) and 0 < α ≤ γ < 2. We consider a nonlocal inverse problem and show that the fractional exponents β, α and γ are determined uniquely by the data u(t, 0) = g(t), 0 < t ≤ T. The existence of the solution for the inverse problem is proved using the quasi-solution method which is based on minimizing an error functional between the output data and the additional data. In this context, an input-output mapping is defined and its continuity is established. The uniqueness of the solution for the inverse problem is proved by means of eigenfunction expansion of the solution of the forward problem and some basic properties of fractional Laplacian. A numerical method based on discretization of the minimization problem, namely the steepest descent method and a least squares approach, is proposed for the solution of the inverse problem. The numerical method determines the fractional exponents simultaneously. Finally, numerical examples with noise-free and noisy data illustrate applicability and high accuracy of the proposed method.