Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 205169, 11 pages

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

Research Article

Application of the RBF Method to the Estimation of Temperature on the External Surface in Laminar Pipe Flow

Shiqin Lyu, Can Wu, and Sufang Zhang

College of Mathematics, Taiyuan University of Technology, Taiyuan 030024, China

Received 6 July 2013; Revised 26 August 2013; Accepted 2 September 2013

Academic Editor: Herb Kunze

Copyright © 2013 Shiqin Lyu et al. 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.

The inverse heat conduction problem on the heat transfer characteristics of cooled/heated laminar flows through finite length thick-walled circular tubes is studied, using temperature measurements taken at several different locations within the fluid in this paper. The method of radial basis functions is coupled with the boundary control technique to estimate the unknown temperature on the external surface of the circular pipe. The main idea of the proposed method is to solve the direct problem instead of solving the inverse problem directly. The temperature data obtained from the direct problem are used to simulate the temperature measurement for the inverse problem, and during the calculation, the Tikhonov regularization and L-curve methods are employed to solve the inverse problem. Therefore, this study also considers the influence of errors in these measurements upon the precision of the estimated results as well as the influence of the locations and number of sensors used upon the accuracy of the estimated results. The results indicate that the accuracy of the estimated results is improved by taking temperature measurements in locations close to the the unknown boundary. Finally, the results confirm that the proposed method is capable of yielding accurate results even when errors in the temperature measurements are present.

1. Introduction

The study of the inverse heat conduction problem on the heat transfer characteristics of cooled/heated laminar flows through a finite length thick-walled circular tubes has great practical significance in industrial production. It can provide a systemic and referential method to solve some common problems in industrial production such as the estimation of the fluid heat consumption and the choice of materials for tubes. So this kind of inverse heat conduction problem has been extensively investigated by many scholars. In 2009, Wei and Li [1] have presented the regularization method to identify the moving boundary in a multilayer domain based on Cauchy data under one-dimensional heat equations. Lu et al. [2] developed a conjugate gradient method to estimate the unknown temperature of fluid near the inner wall of pipe in the two-dimensional inverse heat conduction problem in 2010. In 2011, the conjugate gradient method combined with a finite element method was applied by Lu et al. [3] to study three-dimensional inverse heat conduction problem. A solution was proposed to investigate inverse problem by using a finite difference method and a least square method in literature [4]. Some other methods, such as the method of fundamental solutions [5], were also used to solve this kind of inverse heat conduction problem.

In this work, the inverse heat conduction problem we study is a typical Cauchy problem which is based on Laplace's equation. This class of problem is typically ill-posed. Some numerical methods such as the method of fundamental solutions (MFS) [6], the boundary element method (BEM) [7], and the finite difference method (FDM) [5, 8–10] were commonly used for solving ill-posed problems. All these methods rely heavily on the grid. Tiny changes to the grid structure could have significant impacts on calculations. Therefore, in recent years, many scholars have focussed their attention on meshless methods, especially meshless methods based on radial basis functions (RBF). It is a new and rapidly developing numerical method that follows traditional methods. The most prominent advantage of a meshless method is that it overcomes the reliance on grid, eliminating the division and reconstruction of a grid completely or partially, so the meshless method, like the RBF method, is widely used in the field of numerical calculation. It was proven that the RBF method was a great method for solving the Cauchy problem by Yao et al. [11]. Lai [12] introduced the application of the RBF method in computational electromagnetics, which reflected the wide applicability of the RBF method. Cheng and Cabral [13] applied the RBF method to get a direct solution to the Cauchy ill-conditioned boundary value problem, but the heat source has not been accounted for. Through these studies, it has been found that the RBF method is easier to implement than other methods and the initial conditions were not needed. In addition, compared with BEM, FDM, and finite element methods (FEM) which depended on the grid [8], the RBF method is faster in calculation and higher in accuracy.

This paper is organized as follows. In Section 2, the physical model is briefly introduced and developed. Based on this model, a numerical method is presented in Section 3 to solve the linear system. The temperature attained in the direct problem is used as a known condition for calculating the temperature on external surface of the tube, combining regularization method in inverse problem. In Section 4, two examples are presented to consider the influence of the measurement error, the locations of measurement points, and the number of measurement points upon the accuracy of the estimated results by analysing the temperature distributions and the relative root mean square error of two examples.

2. Physical Model and Governing Equations

The physical model under consideration consists of constant property flow with steady laminar motion in a circular tube. In this system, the pipe wall is assumed to be homogeneous with a constant thermal conductivity , and the fluid is assumed to be incompressible and homogeneous with a constant thermal conductivity . The fluid temperature at the entrance of the pipe is assumed to be uniform and equal to , and is the fluid velocity distribution. An adiabatic condition is applied in the upstream region and the downstream. As the fluid moves through the pipe, some of its heat escapes through the wall of the pipe via heat conduction. This leads to the uncertain and spatially nonuniform distribution of temperature along the surface of the outer wall. The pipe flow and conjugate heat transfer system are shown in Figure 1. Because of the symmetrical characteristics of our model, the domain can be considered one half of the pipe flow, and this problem can be assumed to be two-dimensional.

Figure 1: Physical model.

In the wall region, the governing equations for the temperature field can be described by the following partial differential equations:

According to our assumptions, the boundary conditions can be expressed as follows:

where and are the inner and outer radii of the pipe, respectively.

In the fluid region, the governing equations for the temperature field can be described by the following partial differential equations:

where is the fluid density and is the specific heat at constant pressure. In this model, is given by the following conventional Hagen-Poiseuille expression:

where is the average velocity, and it can be determined from the measured flow rate.

Similarly, according to our assumptions, we have

Notice that since the length of the downstream section is long enough to account for the wall axial conduction effects, at , the following boundary condition is appropriate:

In addition, at , the interface between the pipe wall and the fluid, the following interface conditions are given

In the heat conduction problem for this system, every dependent parameter would, more or less, affect the heat conduction. A parametric study of the individual parameters would require a vast set of results; however, they are not the principal objective of this work; hence, we assume that these parameters were kept constant in our work:

3. Numerical Method

3.1. Direct Problem

In order to reduce the amount of computation and simplify representation, the following scale parameters are introduced to reduce the governing equations and the boundary conditions into the dimensionless forms:

The RBF method is employed to analyze the direct problem. The main aim is to determine the temperature at some points when all the boundary conditions and thermal properties are known. This temperature data with added noise will then be used as measured data for the inverse problem.

For simplicity, the adopted notations are shown in Table 1.

Table 1: The adopted notations.

Firstly, we collocate some points to ensure that the governing equations and boundary conditions are satisfied; that is, take interior points on satisfying the governing equation (1), take interior points on satisfying the governing equation (3), and take collocation points on satisfying the boundary conditions (2) and (5)–(7), respectively.

Next, we assume that the approximate solution of is given by the series

where , , and are constant coefficients to be determined by collocation, and are basis functions. The multiquadric basis function (MQ)

is used throughout this paper, in which is called the shape parameter. It is well known that the value of leads to a method with exponential convergence properties for solving the PDE. Let

where . By fitting (1)–(7), it is easy to deduce a linear system of equations

where

is a coefficient matrix, in which

is a vector with as an vector of all the elements of 1. Consider the following:

Finally, we can determine the vector of unknown coefficients by the Gauss elimination method. Then the temperature at point can be approximated by (10).

3.2. Inverse Problem

The main purpose of the inverse problem is to determine the unknown outer wall temperature by using the temperature measurement value obtained within the fluid in the pipe.

Let be the measuring points in the fluid region, and let be the temperature observation with random noise added at these points. Although higher accuracy can be obtained by selecting a larger number of measurement points the computational and experimental costs increase as increases, so the number of measurement points is usually small. Hence, it is necessary to specify a number of measuring locations which yields an acceptable compromise between the precision and the associated process costs.

Suppose that the function can be approximated by

where are unknown coefficients are basis functions, usually taken to be triangular basis functions or power basis functions. At points , we can obtain

In the matrix notation, (18) can be rewritten as

where and is an matrix with as its th element.

From (13), we have , where represents the generalized inverse of , thus we can get

Let

then

that is, is the last rows of , and

where is by unit matrix; according to (17) and (19), which can be expressed in terms of as

To satisfy the conditions at the measurement points , we require

In matrix form, the values of the unknown coefficients can be obtained from solving the following matrix equation:

where is a by matrix with as its th element and

Let and . Then by substituting (24) into (26), and using some simple calculation, we can obtain that

Since the original inverse problem belongs to this class of ill-posed problems, any small error in measured data can lead to a dramatic change to the solution. Although it is transformed into the form of direct problem, the ill-conditioning of the matrix in (28) still persists. Due to the large condition number of the matrix , most standard numerical methods are unfit to solve the matrix equation (28). Many regularization methods have been developed for solving these kinds of ill-conditioned problems. In our case, the Tikhonov regularization [14] is employed to solve the matrix equation (28). The Tikhonov regularized solution is defined to be the solution to the following penalized least square problem:

where denotes the usual Euclidean norm and is called a regularization parameter. The determination of a suitable value for the regularization parameter is crucial to the accuracy of the method; the L-curve [15] criterion is used in our computation to determine a suitable value of .

Once the coefficient vector is determined, the outer wall temperature of pipe can be evaluated by using (17).

4. Result and Discussion

To demonstrate the effectiveness and stability of our method for this kind of inverse problem, we present some examples. Without loss of generality, in our illustration, the value of the dimensionless length of the pipe, , is specified as equal to 20, which is sufficiently long to observe the thermal development of the flow. The dimensionless radius of the interface between the pipe wall and the fluid inside the pipe is set equal to 0.6, and the dimensionless temperature on the external surface of pipe is given as follows

Example 1. Consider the following:

Example 2. Consider the following:

In the present RBF collocation solution, the multiquadric given in (11) is used as the basis functions. The collocation points form a 11 11 grid as shown in Figure 2. Different types of collocation points are marked in different symbols. On the interior points, shown in symbol, the governing equation is collocated. On the four sides of the boundary and the interface, shown in symbol, the Dirichlet and the Neumann boundary conditions are collocated.

Figure 2: RBF collocation nodes.

The dimensionless temperature distribution and error distribution for the outer wall are shown in Figures 3 and 4, the left for Example 1 with the shape parameters and , and the right for Example 2 with the shape parameters and . The shape parameters are chosen by minimizing the sum of residual on the given Dirichlet boundary and . In Figure 3, The solid line in black represents the exact solution and the dotted line with a plus at each data point represents the approximate solution by the method outlined in Section 3.1. It can be seen that there is almost no difference between the two lines. In Figure 4, the error of   at each point is defined as the difference between the approximate solution and the exact solution. From the numerical results demonstrated in Figures 3 and 4, we can observe that the proposed scheme in Section 3.1 gives a reasonable approximation to the exact solution . The maximum values for the errors in Examples 1 and 2 are less than and , respectively.

Figure 3: The exact and approximate dimensionless temperature distribution for outer wall.

Figure 4: The error distribution for the outer wall.

For the inverse problem, generally speaking, the results are mainly affected by the following three factors: measurement error, the location of measured points, and the number of measured points. In order to validate the effectiveness and stability of the method discussed in this paper, the influence of these three factors for calculation results is analysed as follows.

For numerical error estimation, we define the relative root mean square error of the boundary on which the temperature data are unknown, by the following formula

There are no real temperature measurements taken in advance, so the temperature which is used for calculation in the inverse problem is obtained from direct problem. In fact, due to the different degree of accuracy of the measuring instruments and the use of different measured methods, the errors are produced more or less; therefore, the simulated temperature tend to consider the influence of the measurement error in the inverse problem. The temperature obtained in the direct problem is taken to be the exact temperature. We add a random error disturbance to these exact values and use this noisy data in the inverse problem. The temperature observation values with random noises at measurement points can be expressed as

where is a random number between and 1 with zero mean and indicates the error level. In our computations, the random variable has been calculated by employing the unifrnd function in MATLAB.

4.1. The Influence of the Measured Error for the Results

To illustrate the impact on the estimated results of different measured error on the estimated results, the value of measuring points is fixed at , where is the dimensionless radius of the interface between the pipe wall and the fluid inside the pipe. In our calculation, we take 9 measuring points on the line of  ; that is,

The basis functions are taken as power basis function with .

Figures 5 and 6 demonstrate the influence of different measured error level for Examples 1 and 2, respectively. The solid lines with no symbol represent the exact distribution of  , the lines with a plus symbol represent the approximate distribution of   without noise, the lines with a star symbol represent the approximate distribution on the level of , and the lines with a circle symbol represent the approximate distribution on the level of .

Figure 5: Estimated dimensionless outer-wall temperature distributions for Example 1.

Figure 6: Estimated dimensionless outer-wall temperature distributions for Example 2.

From Figures 5 and 6, we observe that the estimated lines are consistent with the actual shape. These show that this method can obtain accurate results even with an error rate of up to . Furthermore, we can easily see that the deviation from the exact solution increased with the increasing of measurement error.

The relative root mean square error on the boundary , obtained for Examples 1 and 2 using various levels of noisy data, is presented in Figure 7, on the left for Example 1 and on the right for Example 2. It can be seen that the values for are very small, even for a relatively high amount of noise () added into the measurement data. In other words, the numerical result for the temperature on the outer wall represents a good approximation when compared to its analytical value. Hence our method provides stable numerical solutions to this kind of inverse boundary value problem.

Figure 7: The relative root mean square error of the outer wall for different noise level.

4.2. The Influence of the Location of Measured Points on the Results

In this subsection, to investigate the relationship between the locations of the measurement point and the accuracy of the corresponding estimated results, in our simulation, the error level is fixed at , and the number of measurement points is fixed at .

The following three different types of locations of sensors are considered:

(i) is uniformly distributed in at intervals of 0.1, and ;(ii) is uniformly distributed in at intervals of 0.1, and ;(iii) is uniformly distributed in at intervals of 0.1, and .

Similar to the previous subsection, the basis functions are taken also as power basis function with .

Figures 8 and 9 indicate that under Examples 1 and 2, respectively, the different values of measuring points have an influence on the deviation of the estimated results and the exact values; that is, the position of the measured points has an influence on the estimated results. In the left figure, the solid lines with no symbol represent the exact distribution of  , the lines with a plus symbol represent the approximate distribution of  , which was obtained by the first type of locations of the measurement points, the lines with a star symbol represent the approximate distribution on the level of , and the lines with a circle symbol represent the approximate distribution on the level of  .

Figure 8: The relationship between the radial distance from the measuring points to the centerline of the pipe and the accuracy of the estimated results for Example 1.

Figure 9: The relationship between the radial distance from the measuring points to the centerline of the pipe and the accuracy of the estimated results for Example 2.

The relative error distribution of the three types of locations of the measurement points is given on the right of Figures 8 and 9. In addition, the relative root mean square error on the boundary , which was obtained for Examples 1 and 2 using each of the three types of locations of measuring points is presented in Table 2. From this table, we can also find that the variations of the deterministic bias are influenced by the radial distance from the measuring points to the center line of the pipe, and we can conclude that the accuracy of the inverse method is improved when the sensors are located closer to the interface between solid and fluid. As the locations of the sensors approach the outer wall whose unknown boundary conditions are to be predicted, the accuracy of the estimated results increases correspondingly.

Table 2: The relationship between the relative root mean square error and the value of .

Example 1

Example 2

4.3. The Influence of the Number of Measured Points on Results

During the process of estimating the temperature on the external surface of the laminar pipe flow, the number of selected points also affected our results. To analyse the influence of the number of measurement points on the estimated results, the error level is fixed at and the basis functions are still taken as power basis function with . Let denotes the number of measuring points. The following situations regarding the number of sensors are considered:

(i), is uniformly distributed in at intervals of 0.1, and or 0.3;(ii), is uniformly distributed in at intervals of 0.1, and ;(iii), is uniformly distributed in at intervals of 0.2, and .

Figures 10 and 11 present a comparison between the estimated results and the exact solutions of the temperature distributions obtained using the three situations above for a measurement error of . In the left figure, the solid lines with no symbol represent the exact distribution of , the lines with a plus symbol represent the approximate distribution of , which was obtained by 18 measurement points, the lines with a star symbol represent the approximate distribution on the level of , and the lines with a circle symbol represent the approximate distribution on the level of . The corresponding relative root mean square errors for different number of measuring points are given on the right of Figures 10 and 11.

Figure 10: The relationship between the number of the measuring points and the accuracy of the estimated results for Example 1.

Figure 11: The relationship between the number of the measuring points and the accuracy of the estimated results for Example 2.

We observe that estimations are accurate and robust for the first type (i.e., 18 measuring points) and that the accuracy of the results using the second type temperature measurements (i.e., 9 measuring points) is superior to that obtained using the third type (i.e., 5 measuring points). This indicates that the accuracy of the estimated results increases as the amount of input information (number of sensors) is increased. It also can be concluded that the proposed inverse method is capable of yielding accurate results despite the fewer measurement points.

From Table 3, it is again confirmed that the accuracy of the inverse method is improved with an increase in the number of measuring points. Even if the number of sensors is small, provided that (where denotes the number of unknown coefficients), we also can get very good accuracy. Thus, in practical applications, we can install fewer sensors to save costs without sacrificing accuracy. When the number of sensors is less than the number of unknown coefficients (i.e. ), the estimation accuracy drops dramatically; hence, the number of sensors installed must be at least .

Table 3: The relationship between relative root mean square error and the number of the measuring points.

Example 1

Example 2

5. Conclusions

This paper has successfully studied the use of RBF collocation method and boundary control technique in estimating the unknown outer surface temperature of thick-walled circular tubes with cooled/heated laminar flows.

Experiments indicate that when the measured points are closer to the unknown boundary, the accuracy of our estimation is greater. For a larger measurement error, the increasing number of measured points can improve the accuracy of the estimated value. Even if the measurement error is as large as , we can still get very accurate results. One advantage of this approach is that it does not require the initial temperature. Computationally, this method is more simple and efficient than traditional methods. The unknown boundary can be directly obtained by only one operation, and the repeated iteration process is also avoided. The experiment suggests that the RBF method is an accurate, stable, and effective method for estimating the unknown outer surface temperature of circular tubes.

Acknowledgment

The authors are very indebted to the editors for their valuable comments. The work described in this paper was partially supported by Shanxi Science and Technology Infrastructure Project (2012091003-0101).

References

T. Wei and Y. S. Li, “An inverse boundary problem for one-dimensional heat equation with a multilayer domain,” Engineering Analysis with Boundary Elements, vol. 33, no. 2, pp. 225–232, 2009.

T. Lu, B. Liu, P. X. Jiang, Y. W. Zhang, and H. Li, “A two-dimensional inverse heat conduction problem in estimating the fluid temperature in a pipeline,” Applied Thermal Engineering, vol. 30, no. 13, pp. 1574–1579, 2010.

T. Lu, B. Liu, and P. X. Jiang, “Inverse estimation of the inner wall temperature fluctuations in a pipe elbow,” Applied Thermal Engineering, vol. 31, no. 11-12, pp. 1976–1982, 2011.

R. Pourgholi and M. Rostamian, “A stable numerical algorithm for solving an inverse parabolic problem,” Journal of Information and Computing Science, vol. 4, pp. 290–298, 2009.

Y. S. Li and T. Wei, “Identification of a moving boundary for a heat conduction problem in a multilayer medium,” Heat and Mass Transfer, vol. 46, no. 7, pp. 779–789, 2010.

F. L. Yang and L. Ling, “On numerical experiments for Cauchy problems of elliptic operators,” Engineering Analysis with Boundary Elements, vol. 35, no. 7, pp. 879–882, 2011.

N. S. Mera, L. Elliott, D. B. Ingham, and D. Lesnic, “A comparison of different regularization methods for a Cauchy problem in anisotropic heat conduction,” International Journal of Numerical Methods for Heat and Fluid Flow, vol. 13, no. 5, pp. 528–546, 2003.

T. P. Fredman, “A boundary identification method for an inverse heat conduction problem with an application in ironmaking,” Heat and Mass Transfer/Waerme- und Stoffuebertragung, vol. 41, no. 2, pp. 95–103, 2004.

K. Iijima, “Numerical solution of backward heat conduction problems by a high order lattice-free finite difference method,” Journal of the Chinese Institute of Engineers, vol. 27, no. 4, pp. 611–620, 2004.

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.

G. M. Yao, Siraj-ul-Islam, and B. Sarler, “Assessment of global and local meshless methods based on collocation with radial basis functions for parabolic partial differential equations in three dimensions,” Engineering Analysis with Boundary Elements, vol. 36, no. 11, pp. 1640–1648, 2012.

S. J. Lai, The Radial Basis Meshless Method in Computational Electromagnetics, University of Electronic Science and Technology of China, Chengdu, China, 2010.

A. H.-D. Cheng and J. J. S. P. Cabral, “Direct solution of ill-posed boundary value problems by radial basis function collocation method,” International Journal for Numerical Methods in Engineering, vol. 64, no. 1, pp. 45–64, 2005.

L. Y. Zhang, Y. H. Zhang, and K. M. Qian, “On the iterated Tikhonov regularization for ill-posed problems,” Journal of Shandong University, vol. 4, no. 46, pp. 29–33, 2011.

P. C. Hansen and D. P. OLeary, “The use of the L-curve in the regularization of discrete ill-posed problems,” SIAM Journal on Scientific Computing, vol. 14, pp. 1487–1503, 1993.