Adjoint equation method for solving the inverse diffusion source problem

V. A. Litvinov ; V. V. Uchaikin

Vestnik Ûžno-Uralʹskogo gosudarstvennogo universiteta. Seriâ, Matematika, mehanika, fizika, Tome 11 (2019) no. 3, pp. 20-27

Résumé

The objects of the research are differential equations of diffusion (or thermal conductivity) kind. The subject of research is the algorithm for determining the function of the source or the initial conditions of the problem as per the experimentally measured values. The approach is based on a dual representation of functionals corresponding to experimentally observed quantities in the processes of mass and heat transfer. The inverse problem is formulated in the form of integral equations of the first kind, the core of which is the adjoint function (importance function) obtained as a solution of the adjoint (in the Lagrange sense) diffusion (thermal conductivity) equation with the detector sensitivity function in the right-hand side. Meanwhile, solving the adjoint equations by changing the variables is reduced to solving direct equations. To regularize the solution of the Volterra equation of the first kind corresponding to the problem of recovering the dependence of the boundary condition on time, the residual minimization for an overdetermined system of linear equations has been proposed to be used. The problem of reconstructing the dependence of the initial condition on the coordinate is formulated as a Fredholm equation of the first kind, the solution to which has been obtained using the Tikhonov regularization method. The results of model calculations are presented for the restoration of the time dependence of the sources set by a smooth function, a step function and a function with a harmonic component in the problem of one-dimensional diffusion in a homogeneous medium. These results prove that with the selected calculation parameters, the solutions obtained by the proposed method behave regularly and have quite an acceptable accuracy, even though the values of the sought-for function within the set search interval change by six orders of magnitude. This is seen by the authors as the main difference of their method from other approaches to solving this problem.

Détail

Keywords: inverse problem, thermal conductivity, importance-function, sensitivity.
Mots-clés : diffusion

Bibliographie

[1] Sobolev S., “Functionally-invariant solutions of wave equation”, Travaux Inst. Physico-Math. Stekloff, 5, 1934, 259–264 (in Russ.) | Zbl

[2] Kolmogorov A.N., “On the analytic methods of probability theory”, Uspekhi Mat. Nauk, 1938, no. 5, 5–41 (in Russ.)

[3] Lavrent'ev M.M., Reznitskaya K.G., Yakhno V.G., One-dimensional inverse problems of mathematical physics, Nauka Publ., Novosibirsk, 1982, 88 pp. (in Russ.) | MR

[4] Samarskiy A.A., Vabishchevich P.N., Numerical methods for solving inverse problems of mathematical physics, Izd-vo LKI publ., M., 2009, 478 pp. (in Russ.)

[5] Kabanikhin S.I., Inverse and ill-posed problems, Sibirskoe nauchnoe izd-vo Publ., Novosibirsk, 2009, 457 pp. (in Russ.)

[6] Goncharskiy A.V., Cherepashchuk A.M., Yagola A.G., Numerical methods for solving inverse problems of astrophysics, Nauka Publ., M., 1978, 335 pp. (in Russ.) | MR

[7] Kalinina A.E., “Numerical study of inverse problem of identification of source density for the twodimensional non-stationary convection-diffusion equation”, Dal'nevost. Mat. Zh., 5:1 (2004), 89–99 (in Russ.)

[8] Yaparova N.M., “Numerical Method for Solving some Inverse Heat Conduction Problems with Unknown Initial Conditions”, Bulletin of the South Ural State University. Ser. Computer Technologies, Automatic Control, Radio Electronics, 15:2 (2015), 55–65 (in Russ.) | DOI | MR

[9] Yaparova N.M., “Method for Solving an Inverse Term Source Problem Based on the Laplace Transform”, Bulletin of the South Ural State University. Series: Computational Mathematics and Software Engineering, 5:3 (2016), 20–35 (in Russ.) | DOI

[10] Kozhanov A.I., “Inverse problems of recovering the right-hand side of a special type of parabolic equations”, Mathematical notes of NEFU, 23:4 (2016), 31–44 (in Russ.)

[11] A. Hasanov, M. Slodicka, “An analysis of inverse source problems with final-time measured output data for the heat conduction equation: A semigroup approach”, Applied Mathematics Letters, 26:2 (2013), 207–214 | DOI | MR | Zbl

[12] A. Hazanee, M.I. Ismailov, D. Lesnic, N.B. Kerimov, “An inverse time-dependent source problem for the heat equation”, Applied Numerical Mathematics, 69 (2013), 13–33 | DOI | MR | Zbl

[13] L. Yang, M. Dehghan, J.-N. Yu, G.-W. Luo, “Inverse problem of time-dependent heat sources numerical reconstruction”, Mathematics and Computers in Simulation, 81:8 (2011), 1656–1672 | DOI | MR | Zbl

[14] Marchuk G.I., Computational Mathematics Methods, Mir Publ., M., 1980, 430 pp. (in Russ.)

[15] G.I. Marchuk, Adjoint equations and analysis of complex systems, Mathematics and Its Applications (MAIA), 295, Springer, Dordrecht, 1995, 468 pp. | MR

[16] Litvinov B.A., Uchaykin B.B., Izvestiya vuzov. Fizika, 29:2 (1986), 128 (in Russ.)

[17] Litvinov B.A., Uchaykin B.B., Izvestiya vuzov. Fizika, 29:12 (1986), 96 (in Russ.)

[18] Lavrent'ev M.M., Doklady AN SSSR, 127:1 (1959), 31–33 (in Russ.) | Zbl

[19] Magnitskiy N.A., Vestnik Moskovskogo universiteta. Seriya 15: “Vychislitel'naya matematika i kibernetika”, 1978, no. 1, 91–98 (in Russ.)