Numerical solution of the imprecisely defined inverse heat conduction problem
Tapaswini Smita, Chakraverty S.†, Behera Diptiranjan
Department of Mathematics, National Institute of Technology Rourkela, Odisha 769 008, India

Corresponding author. E-mail: sne_chak@yahoo.com

Abstract

This paper investigates the numerical solution of the uncertain inverse heat conduction problem. Uncertainties present in the system parameters are modelled through triangular convex normalized fuzzy sets. In the solution process, double parametric forms of fuzzy numbers are used with the variational iteration method (VIM). This problem first computes the uncertain temperature distribution in the domain. Next, when the uncertain temperature measurements in the domain are known, the functions describing the uncertain temperature and heat flux on the boundary are reconstructed. Related example problems are solved using the present procedure. We have also compared the present results with those in [ Inf. Sci. (2008)178 1917] along with homotopy perturbation method (HPM) and [ Int. Commun. Heat Mass Transfer (2012)39 30] in the special cases to demonstrate the validity and applicability.

Keyword: 02.60.Cb; 02.60.–x; 02.90.+p; 07.05.Mh; triangular fuzzy number; double parametric form of fuzzy numbers; uncertain inverse heat conduction; variational iteration method (VIM)
1. Introduction

In heat conduction problems, when the heat flux and/or temperature histories at the surface of a solid body are known as functions of time, then we find the temperature distribution. This is termed as a direct problem. However, in many heat transfer situations, the surface heat flux and temperature histories must be determined from temperature measurements at one or more interior locations. This is an inverse problem. Inverse heat conduction problems (IHCPs) play a very important role in many branches of science and technology. The problem consists in determining the temperature and heat flux. The IHCP is much more difficult to solve than the direct heat conduction problem in which the initial and boundary conditions are given and the temperatures are to be determined. This IHCP equation is analyzed in the literature.[114] Hetmaniok et al.[1] applied the homotopy perturbation method for the solution of inverse heat conduction problem. Moreover, Slota[2] implemented homotopy perturbation method for the solution of one phase. A hybrid inverse scheme was applied by Chen and Wu[3] for estimation of surface conditions for nonlinear IHCPs. Cialkowski et al.[4] used Trefftz non-continuous method for the solution of a stationary inverse heat conduction problem. Kalman filter-enhanced Bayesian back propagation was successfully applied by Deng and Hwang[5] to obtain the solution of IHCPs. Frackowiak et al.[6] developed the fitting algorithm for the solution of inverse problems of heat conduction. Different finite element methods were used by Grvsa and Lesniewska[7] for the solution of IHCPs. Grzymkowski and Slota[8] solved one-phase inverse Stefan problems by the Adomian decomposition method. The Laplace transform was used by Monde et al.[9] to obtain the solution of two-dimensional inverse heat conduction problems. Wroblewska et al.[10] implemented the modified method of elementary balances to obtain the numerical solution of a direct and inverse heat transfer problem. Woodfielda et al.[11] presented the solution of two-dimensional inverse heat conduction problem. An analytical solution based on Greens function has been obtained by Fernandes et al.[12] to reduce computational time in inverse heat conduction problems. The lattice Boltzmann method was used by Kameli and Kowsary[13] to obtain the solution of inverse heat conduction problem. Johansson et al.[14] proposed a new method to obtain the solutions for the radially symmetric inverse heat conduction problem.

In all of the above cited papers, we may observe that the authors have considered the parameters and the initial/boundary conditions as crisp or defined exactly. In addition, many other authors[1532] have also investigated some application problems with crisp variables and parameters. However, in actual practice, rather than the particular value, only uncertain or vague estimates about the variables and parameters may be known. This is because these are found in general by some observations, experiments, or experiences. Hence, to handle these uncertainties and vagueness, one may use interval or fuzzy parameters and variables in the governing differential equations. In this context, some papers related to fuzzy differential equations are given in Refs.  [33]– [37]. The concept of fuzzy derivative was first introduced by Chang and Zadeh.[33] It was followed by Dubois and Prade.[34] The fuzzy differential equations and fuzzy initial value problems were studied by Kaleva[35] and Seikkala.[36] The analytical method is developed by Bede[37] using the Hukuhara derivative to obtain the solution of fuzzy differential equations. Various numerical methods for solving fuzzy differential equations are also introduced in Refs.  [38]– [46]. Ma et al.[38] developed a scheme based on the classical Euler method to solve the fuzzy ordinary differential equations. A two-dimensional differential transform method to solve fuzzy partial differential equations (FPDEs) has been studied by Mikaeilvand and Khakrangin.[39] An algorithm based on α - cut of a fuzzy set has been developed by Akin et al.[40] for the solution of second-order fuzzy initial value problems. Variation of constant formula is handled by Khastan et al.[41] to solve first-order fuzzy differential equations. The concept of generalized H-differentiability is studied by Chalco-Cano and Roman Flores[42] to solve fuzzy differential equations. Behera and Chakraverty[43] obtained the uncertain impulse response of an imprecisely defined half-order mechanical system. Also, the fuzzy hyperbolic reaction-diffusion equation was solved by Tapaswini and Chakraverty[44] to analyze the uncertain forest fire. A new double parametric form of fuzzy number has been developed by Tapaswini and Chakraverty, [45] and then they applied the Adomian decomposition method to obtain the solution of uncertain vibration equation for very large membranes.

Recently, the variational iteration method (VIM) is found to be a powerful tool for the analysis of linear and non-linear physical problems. The VIM was first developed by He[46, 47] and was successfully applied to solve various linear and non-linear differential equations of scientific and engineering problems. Convergence of the VIM is proved in Ref.  [48]. Very recently, VIM has been applied to a wide class of other physical problems.[4955]

In the present analysis, VIM is applied for the numerical solution of the uncertain inverse heat conduction problem. Uncertainty in the initial condition is defined in term of fuzzy numbers. Our literature review reveals that fuzzy differential equations are always converted to two crisp differential equations in general to obtain the solution bounds. However, in the proposed methodology the fuzzy differential equation has been converted to a single crisp differential equation using double parametric form of fuzzy numbers[45] and then the corresponding crisp differential equation is solved to obtain the final fuzzy solution by substituting the parametric values.

This paper is organized as follows. In Section 2, some basic preliminaries related to the present investigation are given. In Section 3, the basic idea of the VIM is discussed. The VIM is applied with the proposed technique in Section 4 to solve uncertain inverse heat conduction problems. In Section 5, fuzzy and interval solutions for different types of initial conditions are given.

2. Preliminaries

In this section, we present some notations, definitions, and preliminaries, which are used further in this paper.[5660]

Definition 1. Fuzzy number A fuzzy number Ũ is a convex normalized fuzzy set Ũ of the real line R such that

where μ Ũ is called the membership function of the fuzzy set and it is piecewise continuous.

Definition 2. Triangle fuzzy number A fuzzy number Ũ is said to be triangular if

(i) there exists exactly one x0R with μ Ũ (x0) = 1 (x0 is called the mean value of Ũ ), where μ Ũ is called the membership function of the fuzzy set; (ii) μ Ũ (x) is piecewise continuous.

Let us consider an arbitrary triangular fuzzy number Ũ = (a, b, c). The membership function μ Ũ of Ũ will be defined as follows:

Definition 3. Singular parametric form of fuzzy number The triangular fuzzy number Ũ = (a, b, c) in single parametric form can be represented as

where r ∈ [0, 1]. It may be noted that the lower and upper bounds of the fuzzy numbers satisfy the following requirements:

(i) u(r) is a bounded left continuous non-decreasing function over [0, 1],

(ii) ū (r) is a bounded right continuous non-increasing function over [0, 1], and

(iii) u(r) ≤ ū (r), 0 ≤ r ≤ 1.

Definition 4. Fuzzy arithmetric For any two arbitrary fuzzy numbers and = [y(r), ȳ (r)] and scalar k, the fuzzy arithmetic is defined as follows:

(i) if and only if and ,

(ii),

(iii), where ,

(iv)

Definition 5. Double parametric form of fuzzy number Using the parametric form as discussed in Definition 3 we have Ũ = [u(r), ū (r)]. Now one may write this in double parametric form as Ũ (r, β ) = β (ū (r) − u(r)) + u(r), where r and β ∈ [0, 1].

3. Formulation of uncertain inverse heat conduction problem

We consider the uncertain heat conduction equation in the domain D as

where a is the thermal diffusivity, ũ is the uncertain temperature, and t and x refer to time and spatial location, respectively.

Let the domain be defined in R2 as D = {(x, t) : x ∈ [0, b], t ∈ [0, t* )} and t* denote the end of the process. On the boundary of the domain, three components are distributed as[1]

where the initial and boundary conditions are given. On the boundary Γ 0, the uncertain initial condition is defined as

On the boundary Γ 1 the uncertain Dirichlet boundary condition is assumed as

In the uncertain inverse problem, the temperature distribution ũ (x, t) in region D along with uncertain temperature and uncertain heat flux ) on boundary Γ 2 are to be determined. On Γ 2 we define the Dirichlet and the Neumann boundary condition as

4. Variational iteration method[46, 47]

To illustrate the basic idea of the technique, we consider the following general nonlinear system:

where L is a linear operator, N is a nonlinear operator, and g(t) is a given continuous function.

The basic character of the method is to construct a correction functional for Eq.  (6) as follows:

where λ is a general Lagrangian multiplier which can be identified via variational theory, un is the n-th approximate solution, and û n denotes a restricted variation, i.e., δ û n = 0. The initial approximation u0 can be freely chosen if it satisfies the initial and boundary conditions of the problem. However, the success of the method depends on the proper selection of the initial approximation u0 . We approximate the solution as

5. Double parametric based uncertain inverse heat conduction problem

First we convert the uncertain heat conduction equation to an interval based uncertain heat conduction equation using single parametric form of fuzzy numbers. Then, the interval based differential equation is reduced to crisp differential equation by using double parametric form[45] of fuzzy numbers. Next, we apply VIM to obtain the solution in double parametric form.

As per the single parametric form, we may write the above fuzzy heat conduction equation as

which subjects to the fuzzy initial condition

Using the double parametric form (as discussed in Definition 5), equation  (8) can be expressed as

subjects to the fuzzy initial condition

where r, β ∈ [0, 1].

Let us now denote

Substituting these into Eq.  (9), we obtain

with the initial condition

Solving Eq.  (10), one may obtain the uncertain temperature as ũ (x, t; r, β ). To obtain the lower and upper bounds of the solution in single parametric form, one may substitute β = 0 and 1, respectively, and thus the solution bounds may be represented as

6. Solution of the fuzzy inverse heat conduction problem using VIM

First we have applied the VIM to solve Eq.  (10). According to VIM, we may construct a correction functional as follows:

By making the above correction functional (i.e., Eq.  (13)) stationary, and noticing that we obtain

Thus, we obtain the Euler– Lagrange equation as

with a natural boundary

Hence, the Lagrange multiplier can be easily identified as follows:

By substituting the identified Lagrange multiplier into Eq.  (13), the following variational iteration formula can be obtained as:

Knowing the uncertain temperature distribution ũ (x, t; r, β ) or its approximation ũ n(x, t; r, β ), one may easily obtain the missing fuzzy boundary conditions as

7. Solution bounds for different fuzzy initial conditions

In this section we have illustrated the solution procedure by solving some example problems.

Example 1 In this example we consider b = 2, a = 1, and k = 1. The double parametric form of triangular fuzzy initial condition is assumed as

We start with an initial approximation ũ 0 = ũ (x, 0; r, β ) given by Eq.  (17) and by the variational iterational formula (16), and have

for j ≥ 4.

Therefore, the exact uncertain temperature distribution in the whole region under consideration can be written as

Knowing the uncertain temperature distribution in the whole area, one can easily determine the functions describing the fuzzy boundary conditions as

To obtain the solution bound, i.e., uncertain temperature distribution in single parametric form, we may put β = 0, 1 for lower and upper bounds of the solution, respectively. Hence, we obtain

We have the function describing the boundary condition in term of bounds as

and

One may note that in the special case, when r = 1, the results (crisp) obtained by the proposed method are exactly the same as that of the solution obtained by Hetmaniok et al.[1] in crisp case.

Example 2 Next, let us take b = 1.2, , k = 2, and the triangular double parametric form of the fuzzy initial condition as

Again, by applying the procedure discussed previously, we obtain the uncertain temperature distribution as

and so on. In the similar manner, higher components of the iteration formula  (16) can be obtained. Finally, the approximate solution may be written as

One can now try to make a higher number of iterations, which result in the convergence of ũ n = (x, t; r, β ) to the exact solution

The corresponding boundary conditions are found to be

By putting β = 0, 1 into ũ (x, t; r, β ), we may obtain the approximate solution of lower and upper bounds, respectively, as

and

Similarly, we obtain the function describing the boundary condition in terms of bounds as

It is interesting to note again that the exact solution obtained by the proposed method for r = 1 agrees completely with that of the exact solution (crisp result) given by Hetmaniok et al.[1] in the crisp case.

8. Numerical results and discussion

Numerical results for uncertain inverse heat conduction problems with different fuzzy initial conditions are computed. Obtained results of the present analysis are compared with the existing solution of Hetmaniok et al.[1] in special case (r = 1). The computed results are shown in figures and tables.

Fuzzy solutions for both Examples  1 and 2 are depicted in Figs.  1 and 2 by varying t from 0 to 5 and keeping x = 1 as constant, respectively, with α = 0.4, 0.6. Interval solution plots for both the examples have been given in Figs.  3– 6. It is also interesting to note from Figs.  3 and 4 that the left and right bounds of the uncertain temperature, i.e., ũ (x, t), (with particular values of x = 1 and t = 1) gradually increase with the increase in t and x, respectively. Solution bounds (viz. interval solution) for uncertain temperature and heat flux on boundary Γ 2 for r = 0.9 are shown in Figs.  7 and 8 and Figs.  9 and 10, respectively, for both of the examples.

 Figure Option Fig.  1. Fuzzy solution of Example 1 for x = 1.

 Figure Option Fig.  2. Fuzzy solution of Example 2 for x = 1.

 Figure Option Fig.  3. Interval solutions for Example 1.

 Figure Option Fig.  4. Interval solutions for Example 1.

 Figure Option Fig.  5. Interval solutions for Example 2.

 Figure Option Fig.  6. Interval solutions for Example 1 at t = 1.

 Figure Option Fig.  7. Interval solutions of uncertain temperature on boundary Γ 2 for Example 1 at r = 0.9.

 Figure Option Fig.  8. Interval solutions of uncertain temperature on boundary Γ 2 for Example 2 at r = 0.9.

 Figure Option Fig.  9. Interval solutions of uncertain heat flux on boundary Γ 2 for Example 1 at r = 0.9.

 Figure Option Fig.  10. Interval solutions of uncertain heat flux on boundary Γ 2 for Example 2 at r = 0.9.

The computed results in digital form for both the examples are also shown in Tables  1 and 3, respectively, by varying x from 0 to 5 and taking t = 1 with r = 0.3. Similar results are tabulated in Tables  2 and 4, respectively, by varying t from 0 to 5 and taking x = 1 with r = 0.3. Also Tables  5 and 6 contains the numerical results for uncertain temperature on boundary Γ 2 for both the examples for r = 0.6. Similarly, Tables  7 and 8 tabulate the obtained results of uncertain heat flux for r = 0.8. It is interesting to note that in the special case for r = 1, the results exactly match with Hetmaniok et al.[1] and Bede[37] with HPM. However, in the Bede method, [37] it requires more computation time because we may have to convert the fuzzy differential equations to a system of differential equations depending on the sign of the coefficient.

Table 1. Comparison of present, Bede[37] with HPM, and Hetmaniok et al.[1] solutions of Example 1 with t = 1.
Table 2. Comparison of present, Bede[37] with HPM and Hetmaniok et al.[1] solutions of Example 1 for x = 1 and r = 0.3.
Table 3. Comparison of present, Bede[37] with HPM, and Hetmaniok et al.[1] solutions of Example 2 for t = 1 and r = 0.5.
Table 4. Comparison of present, Bede[37] with HPM, and Hetmaniok et al.[1] solutions of Example 2 for x = 1 and r = 0.5.
Table 5. Comparison results for uncertain temperature of present, Bede[37] with HPM, and Hetmaniok et al.[1] of Example 1 for r = 0.6.
Table 6. Comparison results for uncertain temperature of present, Bede[37] with HPM, and Hetmaniok et al.[1] of Example 2 for r = 0.6.
Table 7. Comparison results for uncertain heat flux of present, Bede[37] with HPM, and Hetmaniok et al.[1] of Example 1 for r = 0.8.
Table 8. Comparison results for uncertain heat flux of present, Bede[37] with HPM, and Hetmaniok et al.[1] of Example 2 for r = 0.8.
9. Conclusions

In this paper the double parametric form of fuzzy numbers has successfully been applied to the solution of an uncertain inverse heat conduction problem using the variational iteration method (VIM). The double parametric form approach is found to be easy and straight forward. The double parametric form approach is easy and straightforward. This does not change the order of original fuzzy differential equations. However, Bede[37] in his method converted the fuzzy differential equations to two crisp differential equations or system of differential equations depending on the sign of the coefficients. Hence, it takes more time than the present method, which is a double parametric form of fuzzy number. Here, the performance of the method has been shown by using the initial condition as the triangular fuzzy number. Firstly, the uncertain temperature distribution in the domain is computed. Then, reconstruction of functions describing the uncertain temperature and heat flux on the boundary has been done when uncertain temperature measurement is known. It is interesting to note for r = 1 in both the examples that the lower and upper uncertain temperature distributions are equal. From the results, one may conclude that the HPM and VIM can be used as alternative and equivalent methods for obtaining analytical and approximate solutions of forward and inverse heat conduction problems because the solution takes the form of a convergent series with easily computable components.

Acknowledgment

The first author would like to thank the UGC, Government of India, for financial support under the Rajiv Gandhi National Fellowship (RGNF).

Reference
 1 Hetmaniok E, Nowak I, Slota D and Witula R 2012 Int. Commun. Heat Mass Transfer 39 30 [Cited within:7] [JCR: 2.208] 2 Slota D 2011 Numer. Heat Transfer Part A 59 755 DOI:10.1080/10407782.2011.572763 [Cited within:1] [JCR: 1.803] 3 Chen H T and Wu X Y 2007 Numer. Heat Transfer Part B 51 159 [Cited within:1] [JCR: 1.955] 4 Cialkowski M J, Frackowiak A and Grysa K 2007 Int. J. Heat Mass Transfer 50 217 [Cited within:1] [JCR: 2.315] 5 Deng S and Hwang Y 2007 Int. J. Heat Mass Transfer 50 2089 DOI:10.1016/j.ijheatmasstransfer.2006.11.019 [Cited within:1] [JCR: 2.315] 6 Frackowiak A, Botkin N D, Cialkowski M and Hoffmann K H 2010 Int. J. Heat Mass Transfer 53 2123 DOI:10.1016/j.ijheatmasstransfer.2009.12.039 [Cited within:1] [JCR: 2.315] 7 Grvsa K and Lesniewska R 2010 Inverse Problems Sci. Eng. 18 3 DOI:10.1080/17415970903233556 [Cited within:1] 8 Grzymkowski R and Slota D 2006 Comput. Math. Appl. 51 3 [Cited within:1] [JCR: 0.75] 9 Monde M, Arima H, Liu W, Mitsutake Y and Hammad J 2003 Heat Transfer-Asian Research 32 618 DOI:10.1002/(ISSN)1523-1496 [Cited within:1] 10 Wroblewska A, Frackowiak A and Cialkowski M J 2008 Arch. Thermodyn. 1 1 [Cited within:1] 11 Woodfielda P L, Mondeb M and Mitsutakeb Y 2006 Int. J. Heat Mass Transfer 49 187 DOI:10.1016/j.ijheatmasstransfer.2005.06.035 [Cited within:1] [JCR: 2.315] 12 Fernand es A P, Sousa P F B, Borges V L and Guimaraes G 2010 Appl. Math. Model 34 4040 DOI:10.1016/j.apm.2010.04.006 [Cited within:1] [JCR: 1.706] 13 Kameli H and Kowsar F 2012 Int. Commun. Heat Mass Transfer 39 1410 DOI:10.1016/j.icheatmasstransfer.2012.07.032 [Cited within:1] 14 Johansson B T, Lesnic D and Reeve T 2012 Int. Commun. Heat Mass Transfer 39 887 DOI:10.1016/j.icheatmasstransfer.2012.05.011 [Cited within:2] [JCR: 2.208] 15 Jiang M and Wang M Chin. Phys. Lett. 30 010201 DOI:10.1088/0256-307X/30/1/0102012013 [Cited within:1] [JCR: 0.811] [CJCR: 0.4541] 16 Xia L 2011 Chin. Phys. Lett. 28 120202 DOI:10.1088/0256-307X/28/12/120202 [Cited within:1] [JCR: 0.811] [CJCR: 0.4541] 17 Chen J and Zhang Y 2010 Chin. Phys. Lett. 27 026502 DOI:10.1088/0256-307X/27/2/026502 [Cited within:1] [JCR: 0.811] [CJCR: 0.4541] 18 WangQ H, Li Z H, Lai J C and He A Z 2007 Chin. Phys. Lett. 24 107 [Cited within:1] [JCR: 0.811] [CJCR: 0.4541] 19 Cui K and Yang W 2005 Chin. Phys. Lett. 22 2738 [Cited within:1] [JCR: 0.811] [CJCR: 0.4541] 20 Ding L, Liu P G, He J G, Zakaria A and LoVetri J 2014 Acta Phys. Sin. 63 044102(in Chinese) [Cited within:1] [JCR: 1.016] [CJCR: 1.691] 21 MaY M, Lü B, Liu Y, Qia L and Pei 2013 Acta Phys. Sin. 62 204202(in Chinese) [Cited within:1] [JCR: 1.016] [CJCR: 1.691] 22 MaY M, Lü B, Cao X Q, Song J Q, Zhang W M, Zhao Y L and Liu B N Acta Phys. Sin. 62 170504 (in Chinese) DOI:10.7498/aps.62.1705042013 [Cited within:1] 23 Zhang P and Zhang Acta Phys. Sin. 62 1642011(in Chinese) DOI:10.7498/aps.62.1642012013 [Cited within:1] [JCR: 1.016] [CJCR: 1.691] 24 Gu T Rao H and Zhang Q Acta Phys. Sin. 62 169501(in Chinese) DOI:10.7498/aps.62.1695012013 [Cited within:1] [JCR: 1.016] [CJCR: 1.691] 25 Wang Q and Wang M 2014 Chin. Phys. B 23 028703 DOI:10.1088/1674-1056/23/2/028703 [Cited within:1] [JCR: 1.148] [CJCR: 1.2429] 26 Cui C, Liu X and Duan S 2013 Chin. Phys. 22 104501 DOI:10.1088/1674-1056/22/10/104501 [Cited within:1] [JCR: 0.811] [CJCR: 0.4541] 27 Xu W L, Zhang J C, Li X, Xu B Q and Tao G S 2013 Chin. Phys. B 22 010203 DOI:10.1088/1674-1056/22/1/010203 [Cited within:1] [JCR: 1.148] [CJCR: 1.2429] 28 Wang F and Zhang R 2012 Chin. Phys. B 21 050204 DOI:10.1088/1674-1056/21/5/050204 [Cited within:1] [JCR: 1.148] [CJCR: 1.2429] 29 Li L and Zhang Q 2011 Chin. Phys. B 20 100201 DOI:10.1088/1674-1056/20/10/100201 [Cited within:1] [JCR: 1.148] [CJCR: 1.2429] 30 Huang C H and Tsai C C 1998 Heat Mass Transfer 34 47 DOI:10.1007/s002310050230 [Cited within:1] [JCR: 0.84] 31 Fan C L, Sun F R and Yang L 2009 Inverse Prob. Sci. Eng. 17 885 DOI:10.1080/17415970902916813 [Cited within:1] 32 Ozisik M N and Oriand e H R B 2000 Inverse Heat Transfer: Fundamentals and Applications New York Taylor & Francis [Cited within:1] 33 Chang S L and Zadeh L A 1972 IEEE Trans. Syst. Man. Cyber. 2 30 DOI:10.1109/3477.678652 [Cited within:2] 34 Dubois D and Prade H 1982 Fuzzy Sets Syst. 8 225 DOI:10.1016/S0165-0114(82)80001-8 [Cited within:1] [JCR: 1.749] 35 Kaleva O 1990 Fuzzy Sets Syst. 35 389 DOI:10.1016/0165-0114(90)90010-4 [Cited within:1] [JCR: 1.749] 36 Seikkala S 1987 Fuzzy Sets Syst. 24 319 DOI:10.1016/0165-0114(87)90030-3 [Cited within:1] [JCR: 1.749] 37 Bede B 2008 Inf. Sci. 178 1917 DOI:10.1016/j.ins.2007.11.016 [Cited within:5] [CJCR: 1.112] 38 Ma M, Friedman M and Kand el A 1999 Fuzzy Sets Sys. 105 13 DOI:10.1016/S0165-0114(97)00227-3 [Cited within:2] 39 Mikaeilvand N and Khakrangin S 2012 Neural Comput. Applic. 21 S307 DOI:10.1007/s00521-012-0901-x [Cited within:1] [JCR: 1.168] 40 Akin Ö, Khaniyev T, Orucc Ö and Türksen I B 2013 Expert Syst. Appl. 40 953 DOI:10.1016/j.eswa.2012.05.052 [Cited within:1] [JCR: 1.854] 41 Khastan A, Nieto J J and Rodriguez-Lopez R 2011 Fuzzy Sets Syst. 17 20 [Cited within:1] [JCR: 1.749] 42 Chalco-Cano Y and Roman-Flores H 2008 Chaos Soliton. Frac. 38 112119 [Cited within:1] 43 Behera D and Chakraverty S 2014 Annals Fuzzy Math. Inf. 7 401 [Cited within:1] 44 Tapaswini S and Chakraverty S 2014 Fire Safety J. 66 8 DOI:10.1016/j.firesaf.2014.04.002 [Cited within:1] [JCR: 1.222] 45 Tapaswini S and Chakraverty S 2014 The Scientific World J. 2014 1 [Cited within:3] 46 He J H 1999 Int. J. Nonlinear Mech. 34 699 DOI:10.1016/S0020-7462(98)00048-1 [Cited within:3] 47 He J H 2000 Appl. Math. Comput. 34 115 [Cited within:2] [JCR: 0.75] 48 Abbaoui K and Cherruault Y 1995 Comput. Math. Appl. 29 103 [Cited within:1] [JCR: 0.75] 49 Maidi A and Corriou J P 2013 Appl. Math. Comput. 219 8632 DOI:10.1016/j.amc.2013.02.075 [Cited within:1] [JCR: 0.75] 50 Wazwaz A M 2009 Appl. Math. Comput. 212 120 DOI:10.1016/j.amc.2009.02.003 [Cited within:1] [JCR: 0.75] 51 Huang Y J and Liu H K 2013 Appl. Math. Model. 378 118 [Cited within:1] [JCR: 1.706] 52 Mosayebidorcheh S and Mosayebidorche T 2012 Int. J. Heat Mass Transfer 55 6589 [Cited within:1] [JCR: 2.315] 53 Miansari M, Ganji D D and Miansari M 2008 Phys. Lett. A 372 779 DOI:10.1016/j.physleta.2007.08.065 [Cited within:1] [JCR: 1.11] 54 Ganji D D, Afrouzi G A and Talarposhti R A 2007 Phys. Lett. A 368 450 [Cited within:1] [JCR: 1.11] 55 Khaleghi H, Ganji D D and Sadighi A 2007 Num. Heat Transfer Part A 52 25 [Cited within:1] 56 Hanss M 2005 Applied Fuzzy Arithmetic: An Introduction with Engineering Applications Berlin Springer-Verlag [Cited within:1] 57 Jaulin L, Kieffer M, Didrit O and Walter E 2001 Applied Interval Analysis France Springer [Cited within:1] 58 Ross T J 2004 Fuzzy Logic with Engineering Applications New York John Wiley & Sons [Cited within:1] 59 Zimmermann H J 2001 Fuzzy Set Theory and its Application London: Kluwer Academic Publishers [Cited within:1] 60 Chakraverty S 2014 Mathematics of Uncertainty Modeling in the Analysis of Engineering and Science Problems USA IGI Global Publication [Cited within:1]