Abstract

Based on the boundary element method and the decentralized fuzzy inference algorithm, the thermal conductivity in the two-dimensional unsteady-state heat transfer system changing with the temperature is deduced. The more accurate inversion results are obtained by introducing the variable universe method. The concrete method is as follows: using experimental means to obtain the instantaneous temperature in the material or on the boundary, to determine the thermal conductivity of the material by solving the inversion problem. The boundary element method is used to calculate the regional boundary and internal temperature in the direct problem. With the inversion problem, the decentralized fuzzy inference algorithm is used to compensate for the initial guess of the thermal conductivity by using the difference between the temperature measurement and the temperature calculation. In the inversion problem, the influence of the initial guess of different thermal conductivities, different numbers of measuring points, and the existence of measurement errors on the results is discussed. The example calculation and analysis prove that, with different initial guesses, existence of measurement errors, and the number of boundary measurements decrease, the methods adopted in this paper still maintain good validity and accuracy.

1. Introduction

Inversion heat transfer problems refer to the fact that some of the output information of the heat transfer system is obtained through experimental methods to invert some structural features or input information in the system. For example, the inversion of information such as the shape of the temperature boundary layer, the thermal conductivity of the material, and the heat flux density are all typical inversion heat conduction problems. Inversion heat transfer problems have been widely used in many fields such as nondestructive testing, geometrical shape optimization, aerospace engineering, power engineering, mechanical engineering, constructional engineering, bioengineering, metallurgical engineering, material processing, biological medicine, and food engineering, all of which achieve great success [113]. For inversion heat conduction problems, a lot of researches have been done by scholars at home and abroad.

The boundary element method and the complex variable derivation method are applied by Yu to invert the thermal conductivity of heterogeneous materials, which can effectively identify the thermal conductivity of single or multiple parameters [14]. When the heat conduction boundary value of the stability boundary is inverted by Yaparova, Laplace and Fourier transforms [15] are applied. The boundary element method is used to analyze two-dimensional transient conduction problems by Zhou et al. and the conjugate gradient method is introduced to solve the thermal conduction coefficient, which verify the effectiveness and stability of this method [16]. Mera et al. use iterative BEM to generate a stable numerical solution, which increases the number of boundary elements and reduces the amount of noise added in the input data [17]. Chen and Tanaka use a coupling application of the dual reciprocity boundary element method and dynamic programming filter to some inversion heat conduction problems [18]. A new and simple boundary element method is presented by Gao and Wang; this method is called interface integral boundary element method for solving heat conduction problems consisting of multiple media [19]. Wang et al. apply a nonsingular indirect boundary element method for the solution of three-dimensional inversion heat conduction problems. The exact geometrical representation of computational domain is adopted by parametric equations to eliminate the errors in traditional approaches of polynomial shape functions [20].

The differential transformation is studied and a stable differential calculus method is proposed to solve the inversion heat conduction problem [21] by Baranov et al. A new method to invert the thermal conductivity of material with temperature is proposed Miao et al., by which the temperature of measurement point is obtained by finite element method, and the residual between calculated value and measured value of temperature at the measurement point is minimized to get numerical solution, proving the effectiveness and accuracy of the algorithm [22].

The thermal conductivity of material with temperature changes is piecewise discrete by temperature range by Tang et al., and the genetic algorithms and the adjoint equation are used to carry out the inversion [23] of the thermal conductivity of the full temperature range. Based on the semi-infinite one-dimensional thermal model, the thermal conductivity inversion algorithm is studied by Lei; by changing the mathematical model, different intensity of noise is simulated, and the impact of noise on the accuracy of inversion is observed and studied; the method to improve the accuracy is proposed [24].

The decentralized fuzzy inference algorithm is successfully applied to the unsteady-state heat transfer system by Ran, which shows good anti-ill-posedness and obvious advantages and effectiveness [25]. In this paper, the boundary element method is used to solve the boundary and internal temperature in the two-dimensional unsteady-state heat transfer system, and the decentralized fuzzy inference algorithm is used to compensate for the initial guess of the thermal conductivity in order to minimize the residual between the calculated and the measured values of the temperature, and the true thermal conductivity is obtained.

2. Direct Heat Conduction Problem

2.1. The Boundary Integral Equation

The mathematical model of the two-dimensional unsteady-state heat transfer problem [26]:

In the mathematic model formula, there are , and the boundary of domain , which meets . And is the thermal conductivity coefficient , is the density of the object, and is the specific heat capacity of the object. And is the heat transfer coefficient of the object changing with temperature, is the temperature, and is the outer normal vector of the boundary. And is the ambient temperature, and is the heat flux density. The letter with “−” denotes the known quantity.

Weight function is introduced into the expression of weighted residual of governing equation [27].

The left side of the equation is decomposed to get:

In Green’s theorem for the Laplacian, of which is the boundary curve of plane closed Region and is the arc differential.

According to Laplace Green function, the following is obtained:

Equation (4) is taken into (3) and (2) becomes

It is further simplified, and is transformed to get

One Integration by parts () is carried out in the equation for to get

Equation (6) is taken into (5) to get

Because of and , , is discomposed and taken into (8) to get

The corresponding basic solution to this formula iswhere is the dimensionality of space, for two-dimensional problem, .

Differential derivation of (10) is done to get

In the formula, is the vertical distance from the source point “” to the boundary element line.

The basic solution has the following characteristics:

Equations (12) and (11) are taken into (10) and a good merger of similar items is done to get

2.2. The Boundary Element Equation

The change of the function , over time is small enough to be negligible compared to that of , , which can be reasonably approximated as a constant over small time intervals, and (13) can be segmented into time integration [28].

And the interval integral for is

In the formula, is the exponential integral function, which can be calculated by the series, which is

In the formula, is Euler function, , for ; generally the first five approximations are taken.

According to the above formula, (17) can be written as

For the spatial domain division, the boundary is divided into units and the domain is divided into units. Equation (17) becomes

Linear element interpolation is adopted and the interpolation function of linear element is . Therefore, the boundary curve is approximated by a straight line. The values of and in the unit are approximated by the linearity with two endpoint values.

Equation (18) is done as

In the formula, . The formula is written in matrix form:

and at the boundary node can be obtained by (21). Take and the inner point temperature is obtained by (13) and (18).

3. Decentralized Fuzzy Inversion

3.1. Inversion of Thermal Conductivity

The thermal conductivity inverted in this paper varies with the temperature of the material, the function of which is known. By the known measured temperature at a particular measurement point, the inversion algorithm is used to determine the constant coefficients of the function.

The difference between the temperature measurement and the temperature calculation is taken as the objective function, which is minimized as

is the inversion parameter in the objective function; and , respectively, represent the temperature measurement and the temperature calculation at the measuring point at time ; is the number of temperature measuring points; is the number of future time steps; and the minimum value of the objective function is calculated as the parameter vector in the inversion problem.

3.2. Decentralized Fuzzy Inference Method

The difference between the temperature measurement and the temperature calculation is used to correspondingly compensate for the initial guess of the thermal conductivity. Therefore, a multiple-input multiple-output fuzzy inference system is established. Each independent measurement point is a single fuzzy inference unit (FIU). Fuzzy inference unit is shown in Figure 1. The independent fuzzy inference controller is divided into four parts: fuzzy interface, knowledge base, inference engine, and fuzzy decision interface (defuzzy).

The input to fuzzy inference unit is the error of temperature calculation and measurement at the time of at the measuring point on the known boundary condition :

The independent fuzzy inference unit output is a fuzzy inference result corresponding to that of inputting and is a numerical value for compensating for the guess value of the inversion parameter with only one independent measurement point .

Linguistic values of each linguistic variable are defined: seven fuzzy sets are defined on the universe of inputting variable and output variable of and these are and . The linguistic values corresponding to the fuzzy sets are, respectively, , PM, PS, ZO, NS, NM, . There are a lot of files and the rules are formulated flexibly and detailed. However, rules are too many and complex, and the programming is difficult, accounting for more internal storage; rare files correspond to less rules, which can be easy to implement. The disadvantage is that control function becomes less detailed, whose effect cannot be satisfactory. So setting the fuzzy rule base is to take into account both the simplicity and accuracy.

The membership function of each language value and triangle membership function are defined; the shape and distribution of membership functions are shown in Figures 2 and 3.

The knowledge base is mainly composed of two parts, the database and language control rule base. The language control rules are based on the difference between the temperature calculation and the temperature measurement. If , it is proved that the temperature calculation is smaller than the temperature measurement; it is necessary to raise the guess value of the inversion parameter to eliminate the temperature error , and the larger the is, the larger the range ability of the guess value of the inversion parameter will be. When , it is proved that the temperature calculation is larger than the temperature measurement; it is necessary to reduce the guess value of the inversion parameter to eliminate the temperature error . The fuzzy control rules are shown in Table 1.

The fuzzy inference engine is based on the fuzzy input and language control rules and the fuzzy relational equation is solved to obtain the fuzzy output. Mamdani Maximum - Minimum Fuzzy Inference Algorithm is used to determine the fuzzy set of output variables. The set of output variables is from the following formula:

In the fuzzy decision interface, the fuzzy output is done with defuzzification to get a precise control. In the fuzzy set output by the fuzzy inference engine, the center of gravity is used to solve defuzzification:

3.3. Variable Universe

In variable universe, the appropriate universe extension factor is to be selected; the error is changed and some changes are also made to the universe, and through the universe changing with the error changes, the precise control effect can be achieved [29, 30].

In this paper, the thermal conductivity changing with temperature is inverted; the error of the input information has a more important influence on the inversion result. Therefore, effectively reducing the sensitivity of the inversion results to the error information is the prerequisite for obtaining a stable inversion result. Universe increases with decreasing, which makes the partition of universe rough, and the corresponding output becomes more detailed. Therefore, the following variable universe formula is applied:in which is the initial input universe of ; is the universe of after being changed. In this paper, two different sets , are taken based on different assumptions.

3.4. Inversion Process

The process of inversing thermal conductivity is as follows.

(1) Set the number of iterations and take the initial guess of thermal conductivity .

(2) Calculate the direct problem to get temperature calculation at the temperature at point .

(3) From (23), the deviation can be calculated to determine whether the convergence condition is satisfied. If it is satisfied, the iteration is stopped. The value is assumed to be the thermal conductivity; otherwise, the next calculation is performed.

(4) Calculate the real-time universe from (25).

(5) Determine by one-dimensional fuzzy inference unit .

(6) Calculate the new guess of thermal conductivity and return to step (2).

4. The Instance Calculation and Analysis

The schematic diagram is shown in Figure 4. A transient heat conduction problem in a 10 × 8 quadrilateral region is considered and the thermal conductivity meets

is the initial temperature, is the thermal conductivity at temperature , and is an experimentally determined constant. For simplicity of description, the physical properties of the material are set to . The temperature of the four boundaries is 1°C, the initial temperature in the domain is 0°C, the boundary is divided into 36 boundary elements, and the domain is divided into 16 quadrilateral elements. The material thermal conductivity inversion is done under the premise of the temperature at the measuring point in the domain is known. The function of thermal conductivity is known, , the correct thermal conductivity is introduced into the direct problem to get the temperature at the measuring point at different time, and the correct temperature is defined as the temperature measurement. In the inversion process , is the temperature at this position at the previous moment. The number of measuring points is selected as , at each measuring four point temperature calculations of different time are taken, and by inversion the three groups of data of the temperature and the corresponding thermal conductivity are obtained.

When there is a measurement error, temperature measurement at the measuring point is

In this formula, is the random number of the normal distribution and is the standard deviation of the measurement.

4.1. The Impact of the Number of Measuring Points

The initial guess value is taken as , the standard deviation of measurement is taken as , the number of temperature measurement points is taken as , , , and the temperature at the measuring point and the corresponding thermal conductivity are obtained by inversion. The inversion result is shown in Figure 5.

At the measuring point , , , by the least square method and are calculated, and the average relative errors are shown in Table 2. When , the average relative errors of and are 0.28% and 3.5%; when , the average relative errors of and are 0.14% and 3.0%, respectively; when , the average relative errors of and are 0.05% and 2.1%, respectively. And thus it is shown that, by increasing the number of measuring points, the average relative error decreases and the inversion accuracy improves.

4.2. The Impact of Initial Guess

The number of measuring points is , the standard deviation of measurement is , and three different initial guesses, , , and are used, respectively, for numerical test. The inversion result is shown in Figure 6.

, , and are used, respectively, by the least square method to calculate the value of and ; the average relative errors are shown in Table 3. When , the average relative errors of and are 0.02% and 4.4%; when , the average relative errors of and are 0.05% and 2.1%; when , the average relative errors of and are 0.06% and 3.0%. It can be seen that the initial guess has a certain effect on the result, but the satisfactory result can be obtained within a reasonable range.

4.3. The Impact of Measurement Error

The initial guess is , the number of measuring points is , and three groups of standard deviation, , , are used, respectively, for numerical test. The inversion result is shown in Figure 7.

The standard deviations of measurement are , , ; the values of , are calculated by the least square method; the average relative errors are shown in Table 4. When , the average relative errors of and are 0.05% and 2.1%; when , the average relative errors of and are 0.20% and 3.0%; when , the average relative errors of and are 1.40% and 36.1%. It can be seen that there are some measurement errors under the premise of large amount of measurement data, and the inversion results can still be satisfactory. However, the larger the standard deviation of measurement is, the more distorted the inversion results will be.

5. Conclusion

In this paper, the boundary element method and decentralized fuzzy inference algorithm are used to invert the thermal conductivity changing with temperature in two-dimensional unsteady-state heat transfer system. The effects of initial guess, different number of measuring points, and measurement errors on the results are discussed. It proves that if the initial guess is taken within a reasonable range and when there is some measurement error, inversion results can be satisfactory. Through the calculation and analysis of examples, the accuracy and stability of the thermal conductivity inversion algorithm are verified.

Data Availability

In order to better share our research results, we agree to distribute the data publicly. The data relating to this research is available upon request by contacting the corresponding author, Shoubin Wang ([email protected], +86-022-23773022).

Disclosure

The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; and in the decision to publish the results.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Shoubin Wang and Li Zhang proposed mathematical description of boundary element method; Xiaogang Sun proved the decentralized fuzzy inference method; Li Zhang and Huangchao Jia designed the experiments and analyzed the results. Li Zhang wrote part of code.

Acknowledgments

This work was financially supported by the National Key Foundation for Exploring Scientific Instrument of China (2013YQ470767), Tianjin Science and Technology Committee for Science and Technology Development Strategy Research Project (15ZLZLZF00350), Tianjin Science and Technology Commissioner Project (16JCTPJC53000), and the 13th Five-Year Plan (2016–2020) of Science Education Project in Tianjin City (HE1017).