Numerical simulations for initial value inversion problem in a two-dimensional degenerate parabolic equation

: In this paper, we study the inverse problem of identifying the initial value of a two-dimensional degenerate parabolic equation, which often appears in the fields of engineering, physics


Introduction
In this work, we consider an inverse problem of recovering the initial value in a two-dimensional degenerate parabolic equation from terminal observation. The problem can be stated in the following form:  Problem P: Consider the following two-dimensional degenerate parabolic equation a X u X t Q T u X X X (1.1) where () aX is a given smooth function satisfying      ( ) | 0, ( ) 0, , a X a X X and   2 2 = (0, ) l R is a two-dimensional bounded plane region. () X is an unknown initial function in (1.1). Assume that an additional condition is given as: where  0 T is a fixed time, and  () X is the observation data. In this paper, we use X to denote the two-dimensional spatial variable ( , ).
xy The inverse problem is to determine the functions u and  satisfying (1.1) and (1.2).
In recent years, many inverse problems of degenerate parabolic equations have been proposed to meet actual needs (see [2,6,16,29,33]), which often appear in many engineering and financial fields. It is noted that the degenerate parabolic equation has attracted the attention of many scholars and has been successfully applied to porous media (see [1,[3][4][5]25,27]). These results have a profound impact on the development of related fields. Eq (1.1) also belongs to degenerate parabolic equations. Compared with the traditional parabolic equation, the degenerate one can get rid of the limitation of boundary conditions in some cases. According to the famous Fichera theory, whether or not the boundary conditions of the degenerate parabolic equation should be given is determined by the corresponding symbol of the Fichera function (see [23]). If the sign is nonnegative, the boundary condition is not necessary; otherwise, it is the opposite. Using the Fichera's theorem, we can easily calculate that the value of Fichera function of Eq (1.1) is zero at the corresponding degenerate boundary. Therefore, Eq (1.1) does not need to give boundary conditions.
For general heat conduction problems, the initial value function usually represents the initial temperature distribution field of the object. While for the inverse heat conduction problem, the main task is to reconstruct the initial temperature with some additional observation data. This is a kind of severly ill-posed problem, that is, the small change of observation data will lead to the big change of solution. This is easy to understand. In fact, most things in nature are irreversible. If we can easily reconstruct the past temperature by using the current measurement data, then we can easily obtain the temperature of the earth 65 million years ago. Maybe we can know how dinosaurs died out, but it's impossible. (The reason why dinosaurs died out has always been a mystery. One is that 65 million years ago, the earth's climate suddenly changed, the temperature dropped sharply, dinosaurs were frozen to death. Another view is that because of the eruption of volcanoes, a large amount of carbon dioxide spewed out, resulting in the earth's rapid greenhouse effect, making food death, and eventually the extinction of dinosaurs.) However, this is a problem that must be solved in some engineering fields. For example, in steel smelting, we need to obtain the ideal temperature distribution at a certain time, so how to accurately control the initial temperature is such a problem.
For the inverse initial value problems of nondegenerate parabolic equations, there have been a lot of research work, and abundant theoretical and numerical results have been obtained (see [9-15, 17-20,31]). In [17], the data assimilation problem for the soil water movement model is analysed carefully, where the numerical computation is realized through searching the minimum of the cost function iteratively by a descent gradient method.
For the case of degeneration, the literatures are quite few. These works mostly discuss the coefficient inversion of degenerate parabolic equation from the angle of theoretical analysis, such as source term coefficient or zero-order term coefficient, etc. Documents [7,8,22] are earlier ones involving degenerate parabolic equations. In these papers, the authors establish the Carleman estimate in the case of boundary degradation, and discuss the null controllability of related equations. Since then, these works have been widely extended and successfully applied to the inverse coefficient(s) problem of degenerate parabolic equations (see [2,6,29]). In [26], the uniqueness of the inverse source problem of one-dimensional degenerate heat equation is proved by the extremum principle, and the numerical simulation is also given. In [29], the authors consider the inverse problem of determining the insolation function, where the underlying equation is the nonlinear Sellers climate model. The uniqueness and Lipschitz stability of the solution are also obtained. In [2], the inverse source problem for multidimensional parabolic operators of Grushin type 2 12 | | ( ) ( , , ), ( , , ) (0, ), is a bounded open region that contains the origin and thus the principle coefficient is degenerate.
In [32], numerical inversion for an inverse initial value problem is studied carefully, where the mathematical model is governed by a one-dimensional degenerate parabolic equation. In [33], the authors consider an inverse problem of simultaneously reconstructing the initial value and source term in a degenerate parabolic equation. Uniqueness and stability of the solution are obtained by the Carleman estimate. It should be pointed out that the vast majority of literature in the previous work is about one-dimensional degenerate parabolic equation, and there are few relevant literature about higher-dimensional case. The main reason is that in the case of higher dimension, the degree and way of coefficient degeneration is much more complex than that in the case of one dimension. However, the application in high dimension is more extensive than that in one dimension. For example, in the field of computer signal processing, one-dimensional functions can only represent one-dimensional signals, but two dimensional functions can represent more complex information, such as digital images. If it is a color image or a video, the dimension of the function will be higher.
This work studies the inverse initial value problem for a two-dimensional degenerate diffusion equation. In a sense, this is an extension of the one-dimensional form, and its thermal conductivity is also extended from one-dimension to two-dimension. The change of thermal conductivity enables us to take into account both anisotropic diffusion and slow diffusion, but it also brought essential difficulty for theoretical analysis and numerical simulations. In particular, these results in the case of two-dimension can be easily applied to the field of digital image processing. The outline of the manuscript is as follows: In Section 2, based on the finite volume method, a difference scheme of the forward problem is proposed. The stability and convergence of the scheme are proved in Section 3. In Section 4, the Landweber iteration and conjugate gradient method (CGM) are used to solve the numerical solution for the inverse problem. Some typical numerical examples are presented in the last section to show the validity of the inversion method.

Difference scheme of forward problem
The main purpose of this paper is to study the problem P from the perspective of numerical analysis. The conclusions on the theoretical aspects are only briefly introduced, and those proofs can be found in the references (see [20,23]). Let where  is an arbitrarily small positive constant. If a and  satisfy  [23]). In this paper, we always assume that the solution of the direct problem has sufficient regularity (at least belongs to . Considering the embedding property of space, it seems that  5 k is enough. Moreover, the uniqueness and stability of the solution for problem P can be proved by the logarithmic convexity method (see [20]).
In this article, we use the finite difference method to obtain the numerical solution of Eq (1.1). For the sake of simplicity, we assume that the domain  is a square. For the case of rectangle, we can discuss it similarly, which is not the key of the problem.
Assumed that the domain  is divided into a mesh JJ  with the spatial step size    Using the finite volume method (see [21,28]), one can derive the following implicit difference scheme of (1.1).
In this section, we always assume that () aX at least belongs to , the boundary conditions can be discretized as , is not zero, we can discuss it similarly.

Case 3. If at least one of the four values
The initial condition x y x l y l can be discretized by Without loss of generality, in the following discussions, the boundary conditions are assumed to be Case 2.
In order to write the difference equation conveniently, we introduce the following symbol /, rh the difference scheme of (2.1) can be rewritten as r a a a a ra ra r a a a a ra C ra r a a a a ra ra r a a a  .
Solving Eq (2.6), we can obtain the numerical solution of the forward problem.

Landweber iterative method
It can be easily seen that the Problem P is a linear problem. Let K be the parameter-to-data mapping, that is , K u T (4.1) where u is the weak solution of (1.1). In this article, we would like to use the Landweber iteration method to solve the operator Eq (4.1). Volume   The proof of stability and convergence for the CGM can be found in [24].

Numerical results
In this section, we try to reconstruct the initial temperature field through the terminal observation data of  tT , which is a serious ill-posed problem. For a given initial condition, the change of temperature will decrease exponentially with the increase of time. So the terminal moment T cannot be too large. We take  0.5 T in the process of numerical simulation.   As can be seen from the above Figures 1, 2, when the number of iterations reaches 1000, the reconstruction effect is remarkably good. The initial guess is taken as zero. Obviously, this initial guess is not good at all, but the convergence of the iterative algorithm is very stable and the error order is about 2 10 . After adding noise, the reconstruction effect is also very good (see the Figures 3, 4).    When the number of CGM iterations reaches 300, the error reaches 3 10 , which shows that the reconstruction effect is satisfactory (see the Figures 5, 6). Furthermore, it shows that the convergence speed of CGM is faster than that of the Landweber iterative method. For the case of adding noise, it can be seen from the figure that the reconstruction effect is also satisfactory (see the Figures 7, 8). . The initial value () X is the original clear image, while the blur image  () X is obtained from the clear image through the above diffusion process. Our task is to restore the original clear image from the blur image. We have tested the two algorithms, and the results are similar. The convergence speed of CGM algorithm is faster than that of Landweber iterative algorithm. The Landweber iterative algorithm needs 100 iterations, while CGM only needs about 50 iterations. The reconstruction effect of the two images is shown in Figures 9-12 (since the reconstructed results are similar, we only use the results obtained by the Landweber iterative algorithm).    From these four groups of pictures, we can see that the image restoration effect is very good. Because the fuzzy part is mainly concentrated in the middle of the image (this is true of most images in reality (see [30])), our diffusion function can describe the evolution process quite well. After a period of iterations, the middle blurred part of the image is almost completely restored, and the edge part of the image is also well preserved.

Conclusions
In this paper, the inverse initial value problem for a two-dimensional degenerate parabolic equation is discussed. The finite volume method is used to construct the difference scheme for the forward problem, and the stability and convergence of the difference scheme is also proved. Then, we use the Landweber iteration and CGM to construct the calculation process of the inverse problem and do numerical experiments. Numerical results show that our algorithm is stable and converges quickly. These methods and results can be widely applied in engineering heat conduction problems and computer graphics and image processing.