Identifying source term in the subdiffusion equation with L^2-TV regularization

Fractional derivatives and integrals (26A33) Smoothness and regularity of solutions to PDEs (35B65) Fractional partial differential equations (35R11) Inverse problems for PDEs (35R30) Regularization by noise (60H50) Numerical solutions of ill-posed problems in abstract spaces; regularization (65J20) Finite difference methods for initial value and initial-boundary value problems involving PDEs (65M06) Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs (65M12) Error bounds for initial value and initial-boundary value problems involving PDEs (65M15) Numerical methods for ill-posed problems for initial value and initial-boundary value problems involving PDEs (65M30) Numerical methods for inverse problems for initial value and initial-boundary value problems involving PDEs (65M32) Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs (65M60) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30)

arXiv classification

cs.NA math.NA

Abstract: In this paper, we consider the inverse source problem for the time-fractional diffusion equation, which has been known to be an ill-posed problem. To deal with the ill-posedness of the problem, we propose to transform the problem into a regularized problem with L^2 and total variational (TV) regularization terms. Differing from the classical Tikhonov regularization with L^2 penalty terms, the TV regularization is beneficial for reconstructing discontinuous or piecewise constant solutions. The regularized problem is then approximated by a fully discrete scheme. Our theoretical results include: estimate of the error order between the discrete problem and the continuous direct problem; the convergence rate of the discrete regularized solution to the target source term; and the convergence of the regularized solution with respect to the noise level. Then we propose an accelerated primal-dual iterative algorithm based on an equivalent saddle-point reformulation of the discrete regularized model. Finally, a series of numerical tests are carried out to demonstrate the efficiency and accuracy of the algorithm.

This page was built for publication: Identifying source term in the subdiffusion equation with L^2-TV regularization

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6367104)