Determining An Unknown Boundary Condition by An Iteration Method

This paper investigates the boundary value in the heat conduction problem by a variational iteration method. Applying the iteration method, a sequence of convergent functions is constructed, the limit approximates the exact solution of the heat conduction equation in a few iterations using only the initial condition. This method does not require discretization of the variables. Numerical results show that this method is quite simple and straightforward for models that are currently under research.


Introduction
The problem of heat conduction with unknown boundary condition is a kind of inverse problem.The inverse problem of the heat conduction equation (IPHCE) arises in some important fields of engineering and science.For instance, because of the space program, starting about 1956, some researchers in References [1][2][3] researched the applications of IPHCE which were related to rocket nozzles, to nose cones of missiles and probes, and other devices; in References [4,5], the researchers researched the solutions of IPHCE for testing of nuclear reactor components; Howse in Reference [6] researched the solidification of glass by the IPHCE; the problem of periodic heating in combustion chambers of internal combustion engines was solved by the IPHCE in Reference [7].For other applications of the IPHCE the reader is referred to References [8,9].
In order to briefly introduce the IPHCE, consider the following one-dimensional heat conduction equation: with the initial condition: and the Dirichlet or Neumann boundary conditions: or Here, t denotes time variable, x denotes space variable, v = v(x, t) represents the temperature distribution function, ∂x 2 , F(x, t) represents heat source (sink), a 1 = a 1 (x, t) represents thermal conductivity, a 2 = a 2 (x, t) represents fluid velocity, and a 3 = a 3 (x, t) represents absorbtion (perfusion) coefficient.
There are four main types of IPHCEs, as follows: • The reversed-time problem (RTP): v(x, 0) is unknown; • the inverse boundary problem (IBP): the Dirichlet or Neumann condition is unknown; • the inverse coefficient problem (ICP): a 1 , a 2 , or a 3 is unknown; and • the inverse heat source (sink) problem (IHSP): F(x, t) is unknown.
Many authors have studied the IPHCE, and some types of IPHCE have been solved by numerical methods.For instance, the authors in Reference [10] applied a numerical method to solve an IBP; in Reference [11], the boundary integral method was applied to solve an IBP; and in Reference [12,13], the authors applied the method of finite volume or boundary element to an IBP.For other methods, the reader is referred to References [14,15].For this article, we use the variational iteration method (VIM) to investigate an IBP.
The structure of this article is as follows.We describe this inverse problem mathematically in Section 2. In Section 3, we briefly introduce the VIM and apply the method to IBP.Illustrative examples are presented in Section 4, and the conclusions are given in Section 5.

Problem Statement
We consider the following two IBPs in this Section.Firstly, we consider the Equations ( 1)-( 3) with a 1 = 1, a 2 = a 3 = F(x, t) = 0 as follows: with the overdetermined condition: Here, x * ∈ (0, 1) is an observation point; f (x) and H(t) are given, continuous functions; and g 2 (t) is a known infinitely differentiable function.The unknown functions v(x, t) and g 1 (t) will be determined.When v(x, 0), v(0, t), and v(1, t) are all given functions which satisfy certain conditions, then the problem in Equations ( 5)-( 8) is a forward problem.According to the existence and uniqueness theorem (see Reference [16]), only one bounded solution of the problem in Equations ( 5)-( 8) exists.
Secondly, we take the Equation (1) with a 1 = 1, a 2 = a 3 = 0, Equations ( 2) and (4) as follows: with the overdetermined condition: Here, x * ∈ (0, 1) is a fixed observation point; F(x, t), v(x, 0), v x (1, t), and H(t) are all given functions; and we will determine the unknown functions v(x, t) and h 1 (t).When F(x, t), v(x, 0), v x (0, t) and v x (1, t) are all given functions which satisfy certain conditions, then the problem in Equations ( 10)-( 13) is a straightforward problem.According to the existence and uniqueness theorem (See Reference [17]), under those conditions only one bounded solution of the problem in Equations ( 10)-( 13) exists.

Analysis of VIM
In [18], VIM was proposed.Using the iteration method, a sequence of convergent functions is constructed, the limit approximates the exact solution of the heat conduction equation in a few iterations by using only the initial condition.This method has been successfully applied to integro-differential equations, Bratu-like equation, and in other fields [19][20][21][22].In Reference [23], researchers studied the convergence of the VIM via Banach's fixed point theorem.
For a brief introduction to the basic theory and the main steps of VIM, consider the following differential equation: Here, L and N respectively represent linear and nonlinear operators, and f (x) represents an inhomogeneous term.
If v n−1 (n ≥ 1) is an approximate solution of Equation (15), it follows that: In order to improve its accuracy, a correction term can be added, which is written as: Here, the subscript n represents the n-th approximation of v, and λ(ξ) represents a generalized Lagrange multiplier [18], which can be optimally identified in as simple form as possible using variational theory, the non-linear part must be considered as a restricted variation term v n−1 [18], then the variation of the nonlinear part is zero (i.e., δ v n−1 = 0).The stationary condition of Equation ( 16) requires that: for arbitrary δv n−1 .
Considering VIM, we can construct the correct functional for the solution in the −t direction for for Equations ( 5) and (10) as follows: and: respectively.Determining λ via variational theory and taking the stationarity condition for Equation ( 19), we have: With integration by parts, we have: For arbitrary δv n−1 , this yields the following stationary conditions: Hence, we find: Then the iterative formulas for Equations ( 5) and ( 10) reduce to: and: respectively.By using Equation ( 23), a sequence of convergent of solutions v n for Equation ( 5) is found by selecting the appropriate v 0 (x, t), which limit the exact solution v.Then, we obtain the boundary condition in Equation (7).Using a similar procedure, a sequence of convergent solutions v n for Equation (10) are found from Equation (24) and then, we obtain the boundary condition in Equation (12).

Illustrative Examples
In this section, we give two examples of IBP with known exact solution to illustrate the efficiency and validity of VIM.
Thus, the problem in Equations ( 5)-( 9) has the solution v = exp(−t) sin(x) when n → +∞, which satisfies the given differential equation.Then, the boundary condition g 1 (t) = 0 is obtained.Hence, our iterative approach recovers the exact solution.

Example 2
Consider the IBP in Equations ( 10)-( 14) with the following known functions:
In References [10][11][12][13], similar IBPs have been solved by other numerical methods.Here, the advantage of the VIM is that it does not need to discretize the variables, so it is not subject to computational error and has no need to face computing time and large computer memory.Additionally, the VIM solves this problem using only the initial condition, and therefore is very practical and effective for the currently researched models.

Conclusions
In this work, we successfully utilized VIM in a heat conduction problem with unknown boundary conditions.All the calculations could be made with simple manipulations.The method provides the solution with high precision and high efficiency in the form of a rapidly converging series by using only the initial condition.The solutions obtained through VIM recover the exact solutions as the number of iterations becomes infinite.VIM solves the IBP easily, while it does not need to discretize the variables.For this reason, it is not subject to computational error and requires neither substantial computing time nor large amounts of computer memory.The results of the preceding examples illustrate that this method is easy to operate and that it is a practical, reliable, and effective iterative technique for the IBP.Thus, we can say that the VIM is very practical and effective for models currently under study.