Identification of a Time-Dependent Coefficient in Heat Conduction Problem by New Iteration Method

This paper investigates the problem of identifying unknown coefficient of time dependent in heat conduction equation by new iteration method. In order to use new iteration method, we should convert the parabolic heat conductive equation into an integral equation by integral calculus and initial condition. This method constructs a convergent sequence of function, which approximates the exact solution with a few iterations and does not need complex calculation. Illustrative examples are given to demonstrate the efficiency and validity.

Nonlocal boundary specifications like (6) arise from many important applications in heat transfer, thermoelasticity, control theory, life science, etc. For example, in a heat transfer process, if we let represent the temperature distribution, then (1)-(4) and (6) can be regarded as a control problem with source control. A source control parameter ( ) needs to be determined so that a desired thermal energy can be obtained for a portion of the spatial domain.
Determination of unknown coefficients in inverse heat conduction problems is well known as inverse coefficient problems (ICPs). Identification of physical properties such as conductivity using measured temperature or heat flux values at some sites is an important inverse problem, and these problems have been studied by many authors and some types of ICPs have been solved by some numerical or analysis methods; for instance, Cannon JR and Lin YP researched the parameter ( ) in some quasilinear parabolic differential equations in [1], Kerimov NB and Ismailov MI studied the existence, uniqueness and continuous dependence upon the data of the solution of the ICPs using the generalized Fourier 2 Advances in Mathematical Physics method in [2], Pourgholi R and Rostamian M solved an ICPs using Tikhonov regularization method in [3], Hussein MS and Lesnic D studied the problem of time and space-dependent coefficients in a parabolic equation using the finite difference method in [4], in [5] Kamynin VL researched the existence and uniqueness theorems of the solution of the inverse problem of simultaneous determination of the righthand side and the lowest coefficient in multidimensional parabolic equations with integral observation, Fatullayev AG and Cula S determined an unknown space-dependent coefficient in a parabolic equation using the finite difference method in [6], Ozbilge E and Demir A researched the inverse problem for a time-fractional parabolic equation using semigroup method in [7], variational iteration method was applied to solve an inverse parabolic equation in [8], and other methods can be referred to literature [9][10][11][12].
Daftardar-Gejji and Jafari [13] proposed the new iteration method (NIM) which is based upon the Adomian decomposition method. This method has been applied to various equations such as algebraic equations, integral equations, and ordinary and partial differential equations of integer and fractional order. This method was compared with other numerical methods by some authors; for instance, in [14] Bhalekar and Daftardar-Gejji applied the NIM to fractional-order logistic equation and compared with the Adomian decomposition and homotopy perturbation method, and the results showed that the NIM converges faster to the approximate solutions. Srivastava and Rai [15] used the NIM and modified Adomian decomposition method to solve the multiterm fractional diffusion equation for different conditions, the results also showed that the NIM is direct and straightforward, and it avoids the volume of calculations resulting from computing the Adomain polynomials.
In this paper, we will use the NIM to solve inverse heat conduction problems (1)- (6). This paper has been organized as follows. In Section 2 the basis idea of the NIM is simply described, and in Section 3 determination of unknown coefficient in inverse heat conduction problem has been verified using the NIM. Illustrative examples have been presented in Section 4 and the conclusion is given in Section 5.

The NIM
Consider the general functional equation where is a continuous nonlinear operator from → ( , ⊆ , is a Banach space), is a known function, and ‖ ( )‖ ≤ , where is a positive constant. Assume that (7) has the following series solution: The nonlinear operator N can be decomposed as [13] ( ) = ( In view of (8) and (9), (7) can be rewritten as Define the recursion relation as follows: From (11), we have and If is a contraction operator, that is, Since the numerical series ∑ ∞ =0 ‖ 0 ‖ is convergent, therefore, the series ∑ ∞ =0 uniformly converges to a solution of (7) [16], which is unique in view of the Banach fixed point theorem [17].

Apply NIM to the Inverse Problem
We begin our investigation with a pair of invertible transformations for (1)-(6): We can rewrite (1)-(4) and (5) or (6) as follows: such that or It is clear that the original inverse problem (1)- (6) In this paper, we will focus on the numerical approach. Now we proceed to approximate solution pair (V, ) by the NIM. By integrating both sides of (22) with respect to from 0 to and using (23), we obtain which of the same form as (17) can be solved by the NIM, where is a linear operator with respect to V( , ). We obtain the recursion relation from (15) and (18) Here, the notation V ( , ) = 2 V ( , )/ 2 . If (29) satisfies the conditions of (17), then the series V = ∑ ∞ =0 V is convergent and the pair (V, ) is found; we can obtain the solution of the inverse problem from (28).
We can summarize the procedure from the analysis above. Firstly, we change the inverse problem (1)

Illustrative Examples
In this section, several examples of inverse heat conduction problems are given to illustrate the efficiency and validity of the NIM.
with = 2, * = 1/2. Using the (31) and letting V 0 ( , ) = V( , 0) = exp( ), we obtain and so on. Generally we obtain the solution in a closed form is and using (26) Applying (28), the exact solution of this inverse problem is It can be seen that the same results are obtained using the variational iteration method [8]. At the same time, it is worth pointing out that the NIM does not need to approximately identify the general Lagrange multipliers via the variational theory. The overall results show the computation efficiency of the NIM for the studied model.
with = 2, * = 1/2. Using (31)and letting V 0 ( , ) = V( , 0) = sin(( /2) ) exp( ), we have and so on. Generally we obtain the solution in a closed form is and using (26) Applying the (28), the exact solution of this inverse problem is and using (26) and so on. Generally we obtain the solution in a closed form is and using (27) Applying (28), the exact solution of this inverse problem is ( , ) = (cos ( ) + ) exp ( ) , We obtain the following conclusion from the examples above.
Comparing the NIM with the variational iteration method, it does not need to approximately identify the general Lagrange multipliers via the variational theory and reduces the computational difficulties; comparing this method with other numerical method, for instance, the finite difference method [4], it does not require discretization of the variables, then it is not effected by computation round-off errors, and one is not faced with necessity of large computer memory and time. It provides the solution with high accuracy and minimal calculation in a rapidly convergent series which lead to the solution in a closed form by using the initial condition only. The solutions obtained are highly in agreement with the exact solutions.

Conclusion
In this work, we have successfully utilized the NIM to an inverse heat conduction problem. It is observed that the present method reduces the computational difficulties of variational iteration method, it does not need to approximately identify the general Lagrange multipliers via complex calculation, and all the calculation can be made in simple manipulations. It does not require discretization of the variables; it is not effected by computation round-off errors and not faced with necessity of large computer memory and time.
It provides the solution with high accuracy and minimal calculation in a rapidly convergent series where the series may lead to the solution in a closed form by using the initial condition only. The solutions obtained are highly in agreement with the exact solutions; thus we can say the NIM is very simple and straightforward for the studied model.

Data Availability
(1) The new iteration method in this paper was used to support this study and is available at doi: 10

Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.