/ / Article

Research Article | Open Access

Volume 2015 |Article ID 256726 | https://doi.org/10.1155/2015/256726

Muhammad Arghand, Majid Amirfakhrian, "A Meshless Method Based on the Fundamental Solution and Radial Basis Function for Solving an Inverse Heat Conduction Problem", Advances in Mathematical Physics, vol. 2015, Article ID 256726, 8 pages, 2015. https://doi.org/10.1155/2015/256726

# A Meshless Method Based on the Fundamental Solution and Radial Basis Function for Solving an Inverse Heat Conduction Problem

Accepted11 Mar 2015
Published14 Apr 2015

#### Abstract

We propose a new meshless method to solve a backward inverse heat conduction problem. The numerical scheme, based on the fundamental solution of the heat equation and radial basis functions (RBFs), is used to obtain a numerical solution. Since the coefficients matrix is ill-conditioned, the Tikhonov regularization (TR) method is employed to solve the resulted system of linear equations. Also, the generalized cross-validation (GCV) criterion is applied to choose a regularization parameter. A test problem demonstrates the stability, accuracy, and efficiency of the proposed method.

#### 2. Mathematical Formulation of the Problem

In this section, we consider the following one-dimensional inverse heat conduction problem:with the following initial and boundary conditions:where and are considered as known functions and is a given positive constant, while and are regarded as unknown functions. So, in order to estimate and , we consider additional temperature measurements and heat flux given at a point , , as overspecified conditions:To solve the above problem, at first, we divide the problem (1)–(3) into two separate problems. The problem is as follows:and the problem is considered as follows: Obviously, the problems A and B are considered as IHCP, where , and , are unknown functions in the problems A and B, respectively.

#### 3. Method of Fundamental Solutions and Method of Radial Basis Functions

In this section, we introduce the numerical scheme for solving the problem (1)–(3) using the fundamental solutions and radial basis functions.

##### 3.1. Method of Fundamental Solutions

The fundamental solution of (1) is presented as below:where is Heaviside unit function. Assuming that is a constant, it can be demonstrated that the time shift functionis also a nonsingular solution of (1) in the domain .

In order to solve an IHCP by MFS, as the problem , since the basis function satisfies the heat equation (1) automatically, we assume that is a given set of scattered points on the boundary . An approximate solution is defined as a linear combination of as follows:where is given by (7) and ’s are unknown coefficients which can be determined by solving the following matrix equation:where and is a known vector. Also, as the fundamental functions are the solution of the heat equation, only the initial and boundary conditions are practiced to make the system of linear equations; that is, is a square matrix which is defined using the initial and the boundary conditions as follows:For more details, see [13, 18].

##### 3.2. Method of Radial Basis Functions

In this section, we consider RBF method for interpolation of scattered data. Suppose that and are a fixed point and an arbitrary point in , respectively. A radial function is defined via , where . That is, the radial function depends only on the distance between and . This property implies that the RBFs are radially symmetric about . Some well-known infinitely smooth RBFs are given in Table 1. As it is observed, these functions depend on a free parameter , known as the shape parameter, which has an important role in approximation theory using RBFs.

Let be a given set of distinct points in domain in . The main idea of using the RBFs is interpolation with a linear combination of RBFs of the same types as follows:where and ’s are unknown coefficients for . Assume that we want to interpolate the given values , . The unknown coefficients ’s are obtained, so that , , which results in the following system of linear equations:where , , and for with entries . Generally, the matrix has been shown to be positive definite (and therefore nonsingular) for distinct interpolation points for infinitely smooth RBFs by Schoenberg’s theorem . Also, using Micchelli’s theorem , it is shown that is an invertible matrix for a distinct set of scattered nodes in the case of MQ-RBF.

#### 4. Method of Fundamental Solutions-Radial Basis Functions

This section is specified to introduce the numerical scheme for solving the problem (1)–(3) using the fundamental solutions and radial basis functions. At first, we are going to figure out how radial basis functions and the fundamental solutions are applied to approximate the solution of the problem , (4). Similarly, this method will be used for solving the problem on domain . Because the one-dimensional heat conduction equation depends on both parameters of and , scattered nodes are considered in the domain . We assume that are two sets of scattered points in the domain , where and and are the interior and the boundary points in domains and , respectively. Also, we assume thatwhereso that .

We suppose that the solution of the problem (4) in can be expressed as follows:where and is the fundamental solution of the heat equation which is described in Section 3.2. Also, and is a radial function which is defined in Section 3.1. To determine the coefficients in (16), we impose the approximate solution to satisfy the given partial differential equation with the other conditions at any point , so we achieve the following system of linear equations:where , , and 0 is zero matrix. Also, the matrix can be subdivided into two submatrices as follows:where is the obtained submatrix of applying the interior points whose entries are defined as follows:because as already mentioned above, the fundamental functions satisfy the heat equation and alsoUsing the boundary points of domain , namely, , the submatrix results, whereIt is essential to keep in mind that inherent errors in measurement data are inevitable. On the other hand, as it was mentioned, a radial basis function depends on the shape parameter , which has an effect on the condition number of . Therefore, the obtained system of linear equations (17) is ill-conditioned. So, it cannot be solved, directly. Since, small perturbation in initial data may produce a large amount of perturbation in the solution, we use the rand function in matlab in the numerical example presented in Section 6 and we produce noisy data as the follows:where is the exact data, is a random number uniformly distributed in , and the magnitude displays the noise level of the measurement data.

#### 5. Regularization Method

Solving the system of linear equations (17) usually does not lead to accurate results by most numerical methods, because condition number of matrix is large. It means that the ill-conditioning of matrix makes the numerical solution unstable. Now, in methods as the proposed method in this study, which is based on RBFs, the condition number of depends on some factors such as the shape parameter , as well. On the other hand, for fixed values of the shape parameter , the condition number increases with the number of scattered nodes . In practice, the shape parameter must be adjusted with the number of interpolating points. Also, the accuracy of radial basis functions relies on the shape parameter. So, in case a suitable amount of it is chosen, the accuracy of the approximate solution will be increased. Despite various research works which are done, finding the optimal choice of the shape parameter is still an open problem . Accordingly, some regularization methods are presented to solve such ill-conditioned systems. Tikhonov regularization (TR) method is mostly used by researchers . In this method, the regularized solution for the system of linear equations (17) is explained as the solution of the following minimization problem:where denotes the Euclidean norm and is called the regularization parameter. Some methods such as -curve , cross-validation (CV), and generalized cross-validation (GCV)  are carried out to determine the regularization parameter for the TR method. In this work, we apply (GCV) to obtain regularization parameter. In this method, regularization parameter minimizes the following (GCV) function:where .

The regularized solution (17) is shown by , in which is a minimizer of . Then the approximate solution for the problem (4) is written as follows:

#### 6. Numerical Experiment

In this section, we investigate the performance and the ability of the present method by giving a test problem. Therefore, in order to illustrate the efficiency and the accuracy of the proposed method along with the TR method, initially, we define the root mean square (RMS) error, the relative root mean square (RES) error, and the maximum absolute error as follows:where is the total number of testing points in the domain , and are the exact and the approximated values at these points, respectively. RMS, RES, and errors of the functions and are similarly defined, as well. In our computation, the Matlab code developed by Hansen  is used for solving the discrete ill-conditioned system of linear equations (17).

Example. For simplicity, we assume that and . By these assumptions, the exact solution of the problem (1)–(3) is given as follows :

Since, for a large fixed number of scattered points , the matrix will be more ill-conditioned and also, smaller values of the shape parameter generate more accurate approximations, we suppose that and the number of interpolating points is . The obtained values of , , , , , , , , and , as well as condition number of for and and various values of , are given in Tables 2 and 4. Numerical results indicate that this method is not depended on parameter . By the assumptions and and with various levels of noise added into the data, Tables 3 and 5 illustrate the relative root mean square error of and at three points , 0.5, and of the interval , respectively. Figures 1 and 3 feature out a comparison between the exact solutions and the approximate solutions for and and various levels of noise added into the data. It is observed that as the noise level increases, the approximated functions will have acceptable accuracy. Figures 2 and 4 display the relationship between the accuracy and the parameter with noise of level added into the data. They elucidate not only that the numerical results are stable with respect to parameter , but also that they retain at the same level of accuracy for a wide range of values . Therefore, the accuracy of the numerical solutions is not relatively dependent on the parameter . In addition, the study of the presented numerical results, via MFS in , indicates that with noise of level , namely, with the noiseless data and for of the interval , with and for . Also, when the data are considered with noise, and , while by the proposed method in this work and with the same assumptions, we achieve and for . Also, with noise of level , we obtain and . Accordingly, the MFSRBF, which is based on the fundamental solutions of the heat equation and radial basis functions, is more accurate in comparison to MFS.

 1.1 0.017391 0.012603 0.044511 0.014968 0.012701 0.024588 0.005081 0.003682 0.019141 0.007528 0.006388 0.019141 0.005061 0.003668 0.014484 0.005682 0.004821 0.014485 0.002907 0.002107 0.008510 0.003587 0.003043 0.008510 0.003463 0.002510 0.013148 0.004864 0.004128 0.013147 0.004472 0.003241 0.021264 0.007954 0.006749 0.021264 0.003613 0.002618 0.012343 0.005781 0.004905 0.012343
 1.1 0.023701 0.011307 0.109059 0.053926 0.041309 0.109060 0.003687 0.001759 0.012035 0.005482 0.004199 0.008311 0.003769 0.001798 0.011163 0.006226 0.004769 0.011164 0.005477 0.002613 0.019829 0.001900 0.001455 0.003587 0.005138 0.002451 0.018525 0.003612 0.002767 0.006490 0.005022 0.002036 0.013977 0.005012 0.002917 0.008241 0.006734 0.003212 0.012163 0.004516 0.003459 0.006278
 0.1

Table 3 indicates the values of the obtained for various levels of noise, different values of and by MFSRBF.

#### 7. Conclusion

The present study successfully applies a new meshless scheme based on the fundamental solution of the heat equation and the radial basis function with the Tikhonov regularization method in order to solve a backward inverse heat conduction problem. At first, we have divided the inverse problem into two separate problems and then by applying the MFSRBF with the TR method on the obtained problems, two unknown functions in this IHCP are approximated simultaneously. Numerical results clarify that the presented method is an exact and reliable numerical technique to solve a BIHCP and also the accuracy of the numerical solution is not relatively depended on the parameter . Accordingly, not only this method can be applied to solve a BIHCP in higher dimensions, but also it is possible to practice it in other inverse problems. Hence, this will be contemplated more in future researches.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

1. J. V. Beck, B. Black-Well, and C. R. St-Clair, Inverse Heat Conduction Ill-Posed Problems, John Wiley & Sons, 1985.
2. J. R. Cannon, The One-Dimensional Heat Equation, Addison-Wesley, 1984. View at: Publisher Site | MathSciNet
3. O. M. Alifanov, Inverse Heat Conduction Problems, Springer, 1994.
4. M. N. Ozisik, Boundary Value Problems of Heat Conduction, Dover Publications, New York, NY, USA, 1989.
5. J. R. Cannon and J. A. Vander Hoke, “Implicit finite difference scheme for the diffusion of mass in porous media,” in Numerical Solution of Partial Differential Equations, vol. 527, pp. 527–539, North-Holland, 1982. View at: Google Scholar
6. B. T. Johansson and D. Lesnic, “Determination of a spacewise dependent heat source,” Journal of Computational and Applied Mathematics, vol. 209, no. 1, pp. 66–80, 2007. View at: Publisher Site | Google Scholar | MathSciNet
7. 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. View at: Publisher Site | Google Scholar | MathSciNet
8. M. D. Buhmann, Radial Basis Functions, Cambridge University Press, Cambridge, UK, 2003.
9. M. D. J. Powell, “Radial basis functions for multivariable interpolation: a review,” in Numerical Analysis, D. F. Griffths and G. A. Watson, Eds., pp. 223–241, Longman Scientific & Technical, Harlow, UK, 1987. View at: Google Scholar
10. Z. M. Wu, “Hermite-Birkhoff interpolation of scattered data by radial basis functions,” Approximation Theory and Its Applications, vol. 8, no. 2, pp. 1–10, 1992. View at: Google Scholar | MathSciNet
11. M. A. Golberg and C. S. Chen, “The method of fundamental solutions for potential hemholtz and diffusion problems,” in Boundary Integral Methods: Numerical and Mathematical Aspects, pp. 103–176, Computational Mechanics Publication, Southampton, UK, 1999. View at: Google Scholar
12. B. Jin and L. Marin, “The method of fundamental solutions for inverse source problems associated with the steady-state heat conduction,” International Journal for Numerical Methods in Engineering, vol. 69, no. 8, pp. 1570–1589, 2007. View at: Publisher Site | Google Scholar | MathSciNet
13. L. Yan, C.-L. Fu, and F.-L. Yang, “The method of fundamental solutions for the inverse heat source problem,” Engineering Analysis with Boundary Elements, vol. 32, no. 3, pp. 216–222, 2008. View at: Publisher Site | Google Scholar
14. S. Chantasiriwan, “Methods of fundamental solutions for time-dependent heat conduction problems,” International Journal for Numerical Methods in Engineering, vol. 66, no. 1, pp. 147–165, 2006. View at: Publisher Site | Google Scholar
15. W. Chen and M. Tanaka, “A meshless, integration-free, and boundary-only RBF technique,” Computers & Mathematics with Applications, vol. 43, no. 3–5, pp. 379–391, 2002. View at: Publisher Site | Google Scholar | MathSciNet
16. V. D. Kupradze and M. A. Aleksidze, “The method of functional equations for the approximate solution of certain boundary-value problems,” USSR Computational Mathematics and Mathematical Physics, vol. 4, pp. 82–126, 1964. View at: Google Scholar | MathSciNet
17. Y. C. Hon and T. Wei, “A fundamental solution method for inverse heat conduction problem,” Engineering Analysis with Boundary Elements, vol. 28, no. 5, pp. 489–495, 2004. View at: Publisher Site | Google Scholar
18. N. S. Mera, “The method of fundamental solutions for the backward heat conduction problem,” Inverse Problems in Science and Engineering, vol. 13, no. 1, pp. 65–78, 2005. View at: Publisher Site | Google Scholar | MathSciNet
19. L. Marin and D. Lesnic, “The method of fundamental solutions for the Cauchy problem in two-dimensional linear elasticity,” International Journal of Solids and Structures, vol. 41, no. 13, pp. 3425–3438, 2004. View at: Publisher Site | Google Scholar
20. L. Marin and D. Lesnic, “The method of fundamental solutions for the Cauchy problem associated with two-dimensional Helmholtz-type equations,” Computers & Structures, vol. 83, no. 4-5, pp. 267–278, 2005. View at: Publisher Site | Google Scholar | MathSciNet
21. E. J. Kansa, “Multiquadrics—a scattered data approximation scheme with applications to computational fluid-dynamics. I. Surface approximations and partial derivative estimates,” Computers & Mathematics with Applications, vol. 19, no. 8-9, pp. 127–145, 1990. View at: Publisher Site | Google Scholar | MathSciNet
22. E. J. Kansa, “Multiquadrics—a scattered data approximation scheme with applications to computational fluid-dynamics—II. Solutions to parabolic, hyperbolic and elliptic partial differential equations,” Computers & Mathematics with Applications, vol. 19, no. 8-9, pp. 147–161, 1990. View at: Publisher Site | Google Scholar | MathSciNet
23. E. J. Kansa, “Exact explicit time integration of hyperbolic partial differential equations with mesh free radial basis functions,” Engineering Analysis with Boundary Elements, vol. 31, no. 7, pp. 577–585, 2007. View at: Publisher Site | Google Scholar
24. G. E. Fasshauer, “Solving partial differential equations by collocation with radial basis functions,” in Surface Fitting and Multires-Olution Methods, A. Le Mhaur, C. Rabut, and L. L. Schumaker, Eds., Vanderbilt University Press, Nashville, Tenn, USA, 1997. View at: Google Scholar
25. M. Tatari and M. Dehghan, “On the solution of the non-local parabolic partial differential equations via radial basis functions,” Applied Mathematical Modelling, vol. 33, no. 3, pp. 1729–1738, 2009. View at: Publisher Site | Google Scholar | MathSciNet
26. 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. View at: Publisher Site | Google Scholar | MathSciNet
27. M. Dehghan and V. Mohammadi, “The numerical solution of Fokker-Planck equation with radial basis functions (RBFs) based on the meshless technique of Kansa's approach and Galerkin method,” Engineering Analysis with Boundary Elements, vol. 47, pp. 38–63, 2014. View at: Publisher Site | Google Scholar | MathSciNet
28. J. A. Kolodziei and A. UsciLowska, “Application of MFS for determination of effective thermal conductivity of unidirectional composites with linearly temperature dependent conductivity of constituents,” Engineering Analysis with Boundary Elements, vol. 36, no. 3, pp. 293–302, 2012. View at: Publisher Site | Google Scholar | MathSciNet
29. L. Hui and L. Jijun, “Solution of Backward Heat problem by Morozov principle and conditional stability,” Numerical Mathematics, vol. 14, no. 2, pp. 180–192, 2005. View at: Google Scholar | MathSciNet
30. I. J. Schoenberg, “Metric spaces and completly monotonic functions,” Annals of Mathematics, vol. 39, pp. 811–841, 1938. View at: Google Scholar
31. H. Wendland, “Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree,” Advances in Computational Mathematics, vol. 4, no. 4, pp. 389–396, 1995. View at: Publisher Site | Google Scholar | MathSciNet
32. G. E. Fasshauer and J. G. Zhang, “On choosing ‘optimal’ shape parameters for RBF approximation,” Numerical Algorithms, vol. 45, no. 1–4, pp. 345–368, 2007. View at: Publisher Site | Google Scholar | MathSciNet
33. R. E. Carlson and T. A. Foley, “The parameter R2 in multiquadric interpolation,” Computers & Mathematics with Applications, vol. 21, no. 9, pp. 29–42, 1991. View at: Publisher Site | Google Scholar | MathSciNet
34. A. H. Cheng, “Multiquadric and its shape parameter—a numerical investigation of error estimate, condition number, and round-off error by arbitrary precision computation,” Engineering Analysis with Boundary Elements, vol. 36, no. 2, pp. 220–239, 2012. View at: Publisher Site | Google Scholar | MathSciNet
35. L.-T. Luh, “The shape parameter in the Gaussian function II,” Engineering Analysis with Boundary Elements, vol. 37, no. 6, pp. 988–993, 2013. View at: Publisher Site | Google Scholar | MathSciNet
36. H. W. Engl, M. Hanke, and A. Neubauer, Regularization of Inverse Problems, vol. 375 of Mathematics and its Applications, Kluwer Academic, New York, NY, USA, 1996. View at: MathSciNet
37. P. C. Hansen, Rank-Defficient and Discrete Ill-Posed Problems, SIAM, Philadelphia, Pa, USA, 1998.
38. G. H. Golub, M. Heath, and G. Wahba, “Generalized cross-validation as a method for choosing a good ridge parameter,” Technometrics, vol. 21, no. 2, pp. 215–223, 1979. View at: Publisher Site | Google Scholar | MathSciNet
39. P. C. Hansen, “Regularization tools: a Matlab package for analysis and solution of discrete ill-posed problems,” Numerical Algorithms, vol. 6, no. 1-2, pp. 1–35, 1994. View at: Publisher Site | Google Scholar | MathSciNet
40. A. Shidfar, Z. Darooghehgimofrad, and M. Garshasbi, “Note on using radial basis functions and Tikhonov regularization method to solve an inverse heat conduction problem,” Engineering Analysis with Boundary Elements, vol. 33, no. 10, pp. 1236–1238, 2009. View at: Publisher Site | Google Scholar | MathSciNet

#### More related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.