Chinese Journal of Mathematics

Volume 2013 (2013), Article ID 761272, 9 pages

http://dx.doi.org/10.1155/2013/761272

Research Article

A Meshless Method to Determine a Source Term in Heat Equation with Radial Basis Functions

Baiyu Wang and Anping Liao

College of Mathematics and Econometrics, Hunan University, Changsha, Hunan 410082, China

Received 11 July 2013; Accepted 13 August 2013

Academic Editors: Y. Feng and G. Toth

Copyright © 2013 Baiyu Wang and Anping Liao. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This paper considers a numerical method based on the radial basis functions for the inverse problem of heat equation; the inverse problem is determining an unknown source term subject to the overdetermination along with the usual initial boundary conditions, and the unknown source term is only time-dependent. The radial basis functions method is a meshless method with high accuracy for the inverse problem. Some numerical experiments using this method are presented and discussed.

1. Introduction

The inverse problem of parabolic equations appears naturally in a wide variety of physical and engineering settings, such as elasticity, hydrology, material sciences, heat transfer, medical imaging, transport problems, and control theory.

In this paper, we consider an inverse problem for determining an unknown source term in a heat equation. There are many numerical methods which have been developed to solve the inverse problem of parabolic equations. One approach is the method of output least squares, which assumes that the unknown source function is a specific functional form depending on some parameters and then seeks to determine optimal parameter values which minimized an error function based on the overdetermined data. Other effective and broadly applicable techniques such as finite difference method, finite elements method, boundary elements method, and finite volume method are widely used in inverse problems; however, the finite elements method approximates the solution by using low-order piecewise polynomials, and in the finite difference method, the derivatives of the solution are approximated by difference quotients. These methods depend on a suitable generation of meshes, which is difficult for problems with very complicated and irregular geometries. The method we used in this paper is the meshless method based on the radial basis functions.

Meshless methods are a class of numerical methods for solving partial differential equations. In these methods, mesh generation on the spatial domain of the problem is not needed, and this property is the main advantage of these techniques over the mesh-dependent methods such as finite difference methods and finite element methods. In addition, the radial basis functions method is a technique for interpolation on scattered data; this method uses the distributed nodal points to approximate the unknown function; the distribution of nodes could be selected regularly or randomly in the analyzed domain, and there is no demand for the geometry of the domain. The meshless method based on the radial basis functions provides a global interpolation formula not only for the solution but also for the derivatives of the solution and possesses a high order of accuracy. A number of investigators have considered the meshless method based on the radial basis functions: Buhmann [1], Wendland [2, 3], Dehghan and Tatari [4, 5], Wu et al. [6, 7], and so forth.

In Section 2, we will give an outline of the radial basis functions. In Section 3, we solve the inverse problem using the meshless method based on the radial basis functions. Numerical experiments will be given in Section 4 to illustrate the application of the method.

2. An Outline of the Radial Basis Functions (RBF)

The radial basis functions method is a meshless method, it has been introduced for interpolation of scattered data, and it is very suitable for the numerical solution of partial differential equations.

Let be the distance between a fixed point and any ; that is, .

Define , and is a univariate function.

Let be a given set of distinct points. According to the definition of the function , we denote . For a given data , the interpolating function can be written as follows:

where is the unknown RBF coefficient and is the number of data points.

In order to compute , assume that we want to interpolate the values ; that is,

By substituting for in (1), we obtain the following linear system.

where

so we can get the unknown coefficient from the linear system . The matrix has been shown to be positive definite and is invertible for distinct set of scattered set of points. But in most cases, the matrix is very ill-conditioned, so we usually use the more stable method instead of the direct method to solve the coefficient .

The well-known radial basis functions are listed in Table 1.

Table 1: Some well-known functions that generate RBFs.

Definition

Gaussian (GA)

Hardy multiquadrics (MQ)

Inverse multiquadrics (IMQ)

Inverse quadric (IQ)

3. The Meshless Method for the Inverse Problem

In this section, we consider a class of inverse problems as follows.

Find and which satisfy

and the additional specification

where , , , and are known functions and is a fixed prescribed point in the admissible set.

If the functions , , and satisfy consistence, then the inverse problem (5)-(6) has a unique solution [8].

This kind of inverse problems of identifying unknown source control function has been studied by several authors [8, 9]. The approaches used in these papers are the finite difference method based on the trace type functional formulation and the boundary element method. In this paper, we use the meshless method based on radial basis functions combined with the trace type functional formulation to solve the inverse problem as follows.

If the functions pair solves the inverse problem (5)-(6), applying the trace type functional formulation, we have

From this, we get

Substituting (8) into (5), we have the following initial boundary value problem:

The approximate function of can be represented as

where is the radial basis functions described in Section 2.

Substituting (10) into (9), the system (9) can be rewritten as

For the solution of the system (11), we apply a simple one-step forward difference formula to the time and obtain

From the iterative method, we can obtain , , , and combining (10) and (8), we get the numerical solution as follows:

and the approximate solution as follows:

4. Numerical Experiments and Discussions

In this section, we give an example to illustrate the feasibility of the method. In the following experiments, we use the Gauss radial basis function and let .

Example. Consider the equation , with the conditions

and we let , , and .

The exact solutions are

At first, in order to get the approximation effect of according to different nodes of , we let and , respectively. The exact solution , the numerical solution , and the error are plotted in Figures 1, 2, and 3.

(a)

(b)

Figure 1: The exact solution and the numerical solution of (a) and the error (b) with , .

(a)

(b)

Figure 2: The exact solution and the numerical solution of (a) and the error (b) with , .

(a)

(b)

Figure 3: The exact solution and the numerical solution of (a) and the error (b) with , .

As seen from the Figures 1–3, the error decreases and the approximation effect becomes better as decreases.

In Figure 4, we plot the exact solution and the numerical solutions at , with and .

(a)

(b)

(c)

(d)

Figure 4: The exact solution and numerical solution at , with , .

Figure 4 shows clearly that the approximate degree of numerical solution is very good for different value of .

Figures 5 and 6 give the numerical solution and the error with , , respectively.

(a)

(b)

Figure 5: The numerical solution of (a) and the error (b) with , .

(a)

(b)

Figure 6: The numerical solution of (a) and the error (b) with , .

As seen from the Figures 5 and 6, the effect of different to the approximate solution is the same as ; that is, the error decreases according to decreasing .

Next, in order to illustrate the effect of the noisy data on numerical solution, we introduce artificial error function into the additional specification data as follows:

where represents the noisy parameter, and we consider two cases in this paper:

(i), (ii) is a random function of uniformly distributed on .

We first consider case (i), Figures 7 and 8 plot the exact solution , the numerical solution , and the error with , , and , respectively.

(a)

(b)

Figure 7: The exact solution and the numerical solution of (a) and the error (b) with , , and .

(a)

(b)

Figure 8: The exact solution and the numerical solution of (a) and the error (b) with , , and .

In case (ii), let the parameter and , respectively; we plot the exact solution , the numerical solution , and the error with , , in Figures 9 and 10.

(a)

(b)

Figure 9: The exact solution and the numerical solution of (a) and the error (b) with , , and .

(a)

(b)

Figure 10: The exact solution and the numerical solution of (a) and the error (b) with , , and .

The results shown in Figures 7–10 illustrate sensitivity of the inverse problem solution with respect to the noisy data, but when the noisy parameter is const, the result is relatively stable, with the parameter decreasing and the error decreasing too. When the noisy parameter is a random function, the result we obtain is relatively worse, but from the figures we can see that the error can be controlled in a certain range.

In the last experiment, in order to illustrate the approximate effect, we define

where and are the exact solution and and are the approximate numerical solution. The error results with no noisy data are listed in Table 2. In Tables 3, 4, and 5, we list the error results with different noisy data in both cases.

Table 2: The errors for different .

0.05

0.01

0.005

0.001

0.0005

0.0001

Table 3: The errors for with different .

0.01

0.001

Table 4: The errors for with different .

0.01

0.001

Table 5: The errors for with different .

0.002

0.001

5. Conclusion

In this paper, the meshless method based on the radial basis functions is used for solving the inverse problem of heat equation with the time-dependent source term. From the numerical experiment, we can see that this method is high-efficiency and stable with no noisy data. When the noisy parameter is constant, the result is worsening, but not much. The worst result is the case that the noisy parameter is a random function; however, the error can be controlled in a certain range.

Acknowledgments

The project is supported by the Hunan Provincial Innovation Foundation for Postgraduate CX2011B156 and by NNSF of China Grants nos. 11271117 and 11201136.

References

M. D. Buhmann, Radial Basis Functions Theory and Implementations, Cambridge University Press, Cambridge, UK, 2003.

H. Wendland, “Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree,” Advances in Computational Mathematics, vol. 4, no. 1, pp. 389–396, 1995.

H. Wendland, Scattered Data Approximation, Cambridge University Press, Cambridge, UK, 2005.

M. Dehghan and M. Tatari, “Determination of a control parameter in a one-dimensional parabolic equation using the method of radial basis functions,” Mathematical and Computer Modelling, vol. 44, no. 11-12, pp. 1160–1168, 2006.

M. Tatari and M. Dehghan, “A method for solving partial differential equations via radial basis functions: application to the heat equation,” Engineering Analysis with Boundary Elements, vol. 34, no. 3, pp. 206–212, 2010.

Z.-M. Wu and R. Schaback, “Local error estimates for radial basis function interpolation of scattered data,” IMA Journal of Numerical Analysis, vol. 13, no. 1, pp. 13–27, 1993.

L. Ma and Z. Wu, “Radial basis functions method for parabolic inverse problem,” International Journal of Computer Mathematics, vol. 88, no. 2, pp. 384–395, 2011.

A. Farcas and D. Lesnic, “The boundary-element method for the determination of a heat source dependent on one variable,” Journal of Engineering Mathematics, vol. 54, no. 4, pp. 375–388, 2006.

A. Fatullayev and E. Can, “Numerical procedures for determining unknown source parameter in parabolic equations,” Mathematics and Computers in Simulation, vol. 54, no. 1–3, pp. 159–167, 2000.