Regularized Reconstruction of Scalar Parameters in Subdiffusion with Memory via a Nonlocal Observation
Andrii Hulianytskyi, Sergei Pereverzyev, Sergii Siryk, Nataliya Vasylyeva
TL;DR
The paper addresses the inverse problem of identifying scalar parameters in two- and multi-term time-fractional operators governing subdiffusion with memory, using a nonlocal observation $\psi(t)=\int_{\Omega}u(x,t)\,dx$. It develops explicit reconstruction formulas for the leading fractional orders $\nu_1$ and $\nu_2$ (and extensions to $\nu_i$ in $M$-term operators), proves uniqueness and stability, and proposes a computational algorithm based on Tikhonov regularization with a quasi-optimality principle to handle noisy data. The approach is grounded in fractional Hölder spaces, Volterra-type reformulations, and Mittag-Leffler calculus, and it yields regularity results for the state $u$ and concrete procedures to recover additional coefficients such as $\rho_i$ when needed. The results are relevant for practical modeling in biological transport (e.g., oxygen distribution) and demonstrate effective parameter recovery from discrete, noisy measurements, highlighting the method’s applicability to nonautonomous subdiffusion with memory.
Abstract
In the paper, we propose an analytical and numerical approach to identify scalar parameters (coefficients, orders of fractional derivatives) in the multi-term fractional differential operator in time, $\mathbf{D}_t$. To this end, we analyze inverse problems with an additional nonlocal observation related to a linear subdiffusion equation $\mathbf{D}_{t}u-\mathcal{L}_{1}u-\mathcal{K}*\mathcal{L}_{2}u=g(x,t),$ where $\mathcal{L}_{i}$ are the second order elliptic operators with time-dependent coefficients, $\mathcal{K}$ is a summable memory kernel, and $g$ is an external force. Under certain assumptions on the given data in the model, we derive explicit formulas for unknown parameters. Moreover, we discuss the issues concerning to the uniqueness and the stability in these inverse problems. At last, by employing the Tikhonov regularization scheme with the quasi-optimality approach, we give a computational algorithm to recover the scalar parameters from a noisy discrete measurement and demonstrate the effectiveness (in practice) of the proposed technique via several numerical tests.
