Picard Successive Approximation Method for Solving Differential Equations Arising in Fractal Heat Transfer with Local Fractional Derivative

Abstract

The Fourier law of one-dimensional heat conduction equation in fractal media is investigated in this paper. An approximate solution to one-dimensional local fractional Volterra integral equation of the second kind, which is derived from the transformation of Fourier flux equation in discontinuous media, is considered. The Picard successive approximation method is applied to solve the temperature field based on the given Mittag-Leffler-type Fourier flux distribution in fractal media. The nondifferential approximate solutions are given to show the efficiency of the present method.

(ProQuest: ... denotes non-US-ASCII text omitted.)

Ai-Min Yang 1,2 and Cheng Zhang 3 and Hossein Jafari 4 and Carlo Cattani 5 and Ying Jiao 6

Academic Editor:Ming Li

1, College of Science, Hebei United University, Tangshan, China
2, College of Mechanical Engineering, Yanshan University, Qinhuangdao, China
3, School of Civil Engineering and Architecture, Chongqing Jiaotong University, Chongqing 400074, China
4, Faculty of Basic Sciences, Department of Mathematics, Ayatollah Amoli Branch, Islamic Azad University, Amol 4615143358, Iran
5, Department of Mathematics, University of Salerno, Via Ponte don Melillo, Fisciano, 84084 Salerno, Italy
6, Qinggong College, Hebei United University, Tangshan 063000, China

Received 12 December 2013; Accepted 1 January 2014; 11 February 2014

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.

1. Introduction

Engineering problems can be mathematically described by differential equations. Many initial and boundary value problems associated with differential equations can be transformed into problems of solving some approximate integral equations. Heat transfer is described by theory of integral equations. Integral equation arising in heat transfer with smooth condition is valid for continuous media [1-4]. The common methods for solving the equations of heat transfer are purely mathematical are among them; the finite difference techniques (FDT) [5], the regression analysis (RA) [6], the Adomian decomposition method (ADM) [7], the combined Laplace-Adomian method (CLAM) [8], the homotopy analysis method (HAM) [9, 10], the differential transformation method (DTM) [11], the spline-wavelets techniques (SWT) [12], the boundary element method (BEM) [13], the heat-balance integral method (HBIM) [14, 15], the variational iteration method (VIM) [16], the local fractional variational iteration method (LFVIM) [17], and the Picard successive approximation method (PSAM) [18].

On the other hand, the nanoscale heat problem can be characterized as fractal behaviors. As usual, the materials are called the Cantor materials. Heat transfer in fractal media with nonsmooth conditions is a hot topic. For example, the heat transfer equations in a medium with fractal geometry [19] and fractal domains [20] were considered. The local fractional transient heat conduction equations based upon the Fourier law within local fractional derivative arising in heat transfer from discontinuous media were presented in [21-24].

Fractional calculus was successfully used to deal with the real world problems [25-30]. There is its limit that the operators do not deal with the local fractional continuous functions (nondifferential functions). Hence, the local fractional Fourier flux [21] is not handled by using some approaches from the classical and fractional operators. This paper focuses on analytical solution to local fractional Fourier flux in fractal media by using Picard's successive approximation method [18]. This paper is organized as follows. In Section 2, we give notations to local fractional derivative and integrals and investigate the heat transfer in fractal media. Section 3 is devoted to Picard's successive approximation method based upon local fractional integrals. Analysis solution is shown in Section 4. Conclusions are in Section 5.

2. Heat Transfer in Fractal Media with Local Fractional Derivative

In order to study the non-differential solution for the heat problem in fractal media with local fractional derivative, we here begin with the Fourier flux equation in discontinuous media.

The temperature field reads as [21] [figure omitted; refer to PDF] where f(x,y,z,τ) is local fractional continuous at fractal domain Ω .

For a given temperature field T , a local fractional temperature gradient [21] can be written as follows: [figure omitted; refer to PDF] where the local fractional partial derivative is defined by [21-24] [figure omitted; refer to PDF] where ^Δα (f(x,y)-f(_x0 ,y))[congruent with]Γ(1+α)Δ(f(x,y)-f(_x0 ,y)) .

Here, the local fractional derivative is defined on the fractal set like a Cantor set. For example, when we consider the Cantor set, we can find the local fractional derivative of discontinuous function T (however, T is a local fractional continuous function).

We consider the heat flux per unit fractal area q is proportional to the temperature gradient in fractal medium. Fourier law of heat conduction in fractal medium with local fractional derivative is expressed by [21] [figure omitted; refer to PDF] where ^K2α denotes the thermal conductivity of the fractal material, and it is related to fractal dimensions of materials. It is shown that the fractal dimensions of materials are an important characteristic value. Here, we consider the fractal Fourier flow, which is discontinuous; however, it is found that it is local fractional continuous. Like classical Fourier flow, its thermal conductivity is an approximate value for fractal one when α=1 .

Fourier law of one-dimensional heat conduction equation in fractal media reads as [21] [figure omitted; refer to PDF] where ^K2α denotes the thermal conductivity of the fractal materials.

When τ=_τ0 , from (5) we have [figure omitted; refer to PDF] at τ>_τ0 and in A , where ^K2α is the thermal conductivity of the fractal materials. Namely, T is a bi-Lipschitz mapping, and shows the fractal characteristic behavior [21].

Local fractional heat conduction equation with heat generation in fractal media can be written as [21] [figure omitted; refer to PDF] Local fractional heat conduction equation with no heat generation in fractal media is suggested as [21, 22] [figure omitted; refer to PDF] where ^∇2α is a local fractional Laplace operator [21].

3. The Method

In this section, we discuss the Picard successive approximation method. Meanwhile, we transfer the Fourier law of one-dimensional heat conduction equation in fractal media into the local fractional Volterra integral equation of the second kind.

3.1. Picard's Successive Approximation Method

This method is first proposed in [18]. Here we will give a short introduction to Picard's successive approximation method within the local fractional calculus.

In this method, we set [figure omitted; refer to PDF] We give the first approximation _u1 (x) by [figure omitted; refer to PDF] where the local fractional integral of f(x) of order α in the interval [a,b] is defined as follows [21-24]: [figure omitted; refer to PDF] with Δ_tj =_tj+1 -_tj , Δt=max...{Δ_t0 ,Δ_t1 ,Δ_t2 ,...} , and [_tj ,_tj+1 ] , j=0,...,N-1, _t0 =a , and _tN =b .

Here, we find that the equality _u1 (x) is a local fractional continuous function if f(x) , K(x,t) and _u0 (x) are local fractional continuous functions.

Continuing in this manner, we have the infinite sequences of functions [figure omitted; refer to PDF] such that the recurrence equations are given by [figure omitted; refer to PDF] where _u0 (x) is equivalent to any selected function, which is the local fractional continuous function.

Hence, we have successive approximation as follows: [figure omitted; refer to PDF] Thus, at the limit, the solution u(x) is written as [figure omitted; refer to PDF]

3.2. An Alternative Method from Local Fractional Derivative to Local Fractional Volterra Integral Equations

We directly observe that the local fractional differential equation of α order [figure omitted; refer to PDF] can be written immediately as the local fractional Volterra integral equation of thee second kind in the form [figure omitted; refer to PDF] where λ=1/^K2 .

The Mittag-Leffler type Fourier flux distribution in fractal media can be written as follows: [figure omitted; refer to PDF] Making use of (18), we can get the local fractional Volterra integral equation of the second kind in the form: [figure omitted; refer to PDF]

4. Approximate Solutions for Local Fractional Volterra Integral Equation of the Second Kind

Let us assume that the zeroth approximation is [figure omitted; refer to PDF] Then the first approximation can be written as follows: [figure omitted; refer to PDF] Here we obtain the second approximation, which reads as [figure omitted; refer to PDF] Proceeding in this manner, we have the third approximation in the following form: [figure omitted; refer to PDF] Hence, continuing in this manner, we obtain [figure omitted; refer to PDF] Taking the limit, we have [figure omitted; refer to PDF] where the term _Eα [^xα (1+^λα )] is a Mittag-Leffler type Fourier flux distribution in fractal media, which is related to the fractal coarse-grained mass function [21, 24]. When K=1 , we get u(x)=T(a)_Eα [2^xα ] . The nondifferentiable solution of (25) for parameters K=1 , α=ln...2/ln...3 , and T(a)=1 is shown in Figure 1; the non-differentiable solution of (25) for parameters K=1 , α=ln...2/ln...3 , and T(a)=2 is shown in Figure 2; the non-differentiable solution of (25) for parameters K=1 , α=ln...2/ln...3 , and T(a)=3 is shown in Figure 3; the non-differentiable solution of (25) for parameters K=1 , α=ln...2/ln...3 , and T(a)=4 is shown in Figure 4.

Figure 1: The nondifferentiable solution for Mittag-Leffler type Fourier flux distribution for parameters K=1 , α=ln...2/ln...3 , and T(a)=1 .

[figure omitted; refer to PDF]

Figure 2: The nondifferentiable solution for Mittag-Leffler type Fourier flux distribution for parameters K=1 , α=ln...2/ln...3 , and T(a)=2 .

[figure omitted; refer to PDF]

Figure 3: The nondifferentiable solution for Mittag-Leffler type Fourier flux distribution for parameters K=1 , α=ln...2/ln...3 , and T(a)=3 .

[figure omitted; refer to PDF]

Figure 4: The nondifferentiable solution for Mittag-Leffler type Fourier flux distribution for parameters K=1 , α=ln...2/ln...3 , and T(a)=4 .

[figure omitted; refer to PDF]

5. Conclusions

This work studied the Fourier law of one-dimensional heat conduction equation in fractal media. Mittag-Leffler type Fourier flux distribution in fractal media with temperature field effect was considered. An approximation solution for the local fractional Volterra integral equation of the second kind derived from Fourier law of one-dimensional heat conduction equation for heat conduction in discontinuous media was studied by using Picard's successive approximation method. The non-differential approximate solutions were given to show the efficiency of the present method.

Acknowledgments

This work was supported by National Scientific and Technological Support Projects (no. 2012BAE09B00), the National Natural Science Foundation of China (no. 51274270), and the National Natural Science Foundation of Hebei Province (no. E2013209215).

Conflict of Interests

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

References

[1] J. Crank, P. Nicolson, "A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type," Mathematical Proceedings of the Cambridge Philosophical Society , vol. 43, pp. 50-67, 1947.

[2] T. M. Shih, "A literature survey on numerical heat transfer," Numerical Heat Transfer , vol. 5, no. 4, pp. 369-420, 1982.

[3] B. C. Choi, S. W. Churchill, "A technique for obtaining approximate solutions for a class of integral equations arising in radiative heat transfer," International Journal of Heat and Fluid Flow , vol. 6, no. 1, pp. 42-48, 1985.

[4] Y. Z. Povstenko, "Fractional heat conduction equation and associated thermal stress," Journal of Thermal Stresses , vol. 28, no. 1, pp. 83-102, 2004.

[5] M. Dehghan, "The one-dimensional heat equation subject to a boundary integral specification," Chaos, Solitons & Fractals , vol. 32, no. 2, pp. 661-675, 2007.

[6] Y. Ioannou, M. M. Fyrillas, C. Doumanidis, "Approximate solution to Fredholm integral equations using linear regression and applications to heat and mass transfer," Engineering Analysis with Boundary Elements , vol. 36, no. 8, pp. 1278-1283, 2012.

[7] S. Nadeem, N. S. Akbar, "Effects of heat transfer on the peristaltic transport of MHD Newtonian fluid with variable viscosity: application of adomian decomposition method," Communications in Nonlinear Science and Numerical Simulation , vol. 14, no. 11, pp. 3844-3855, 2009.

[8] A.-M. Wazwaz, M. S. Mehanna, "The combined Laplace-Adomian method for handling singular integral equation of heat transfer," International Journal of Nonlinear Science , vol. 10, no. 2, pp. 248-252, 2010.

[9] S. Abbasbandy, "The application of homotopy analysis method to nonlinear equations arising in heat transfer," Physics Letters A , vol. 360, no. 1, pp. 109-113, 2006.

[10] B. Raftari, K. Vajravelu, "Homotopy analysis method for MHD viscoelastic fluid flow and heat transfer in a channel with a stretching wall," Communications in Nonlinear Science and Numerical Simulation , vol. 17, no. 11, pp. 4149-4162, 2012.

[11] A. A. Joneidi, D. D. Ganji, M. Babaelahi, "Differential transformation method to determine fin efficiency of convective straight fins with temperature dependent thermal conductivity," International Communications in Heat and Mass Transfer , vol. 36, no. 7, pp. 757-762, 2009.

[12] C. Cattani, E. Laserra, "Spline-wavelets techniques for heat propagation," Journal of Information & Optimization Sciences , vol. 24, no. 3, pp. 485-496, 2003.

[13] N. Simões, A. Tadeu, J. António, W. Mansur, "Transient heat conduction under nonzero initial conditions: a solution using the boundary element method in the frequency domain," Engineering Analysis with Boundary Elements , vol. 36, no. 4, pp. 562-567, 2012.

[14] J. Hristov, "Approximate solutions to fractional sub-diffusion equations: the heat-balance integral method," The European Physical Journal , vol. 193, no. 1, pp. 229-243, 2011.

[15] J. Hristov, "Heat-balance integral to fractional (half-time) heat diffusion sub-model," Thermal Science , vol. 14, no. 2, pp. 291-316, 2010.

[16] D. D. Ganji, H. Sajjadi, "New analytical solution for natural convection of Darcian fluid in porous media prescribed surface heat flux," Thermal Science , vol. 15, no. 2, pp. 221-227, 2011.

[17] X. J. Yang, D. Baleanu, "Fractal heat conduction problem solved by local fractional variation iteration method," Thermal Science , vol. 17, no. 2, pp. 625-628, 2013.

[18] X.-J. Yang, "Picard's approximation method for solving a class of local fractional Volterra integral equations," Advances in Intelligent Transportation Systems , vol. 1, no. 3, pp. 67-70, 2012.

[19] R. R. Nigmatullin, "The realization of the generalized transfer equation in a medium with fractal geometry," Physica Status Solidi B , vol. 133, no. 1, pp. 425-430, 1986.

[20] K. Davey, R. Prosser, "Analytical solutions for heat transfer on fractal and pre-fractal domains," Applied Mathematical Modelling , vol. 37, no. 1-2, pp. 554-569, 2013.

[21] X. J. Yang Advanced Local Fractional Calculus and Its Applications , World Science, New York, NY, USA, 2012.

[22] Y. Z. Zhang, A. M. Yang, X.-J. Yang, "1-D heat conduction in a fractal medium: a solution by the local fractional Fourier series method," Thermal Science , vol. 17, no. 3, pp. 953-956, 2013.

[23] M.-S. Hu, D. Baleanu, X.-J. Yang, "One-phase problems for discontinuous heat transfer in fractal media," Mathematical Problems in Engineering , vol. 2013, 2013.

[24] X. J. Yang Local Fractional Functional Analysis and Its Applications , Asian Academic, Hong Kong, 2011.

[25] M. Li, W. Zhao, "On bandlimitedness and lag-limitedness of fractional Gaussian noise," Physica A , vol. 392, no. 9, pp. 1955-1961, 2013.

[26] M. Li, "Approximating ideal filters by systems of fractional order," Computational and Mathematical Methods in Medicine , vol. 2012, 2012.

[27] M. Li, W. Zhao, "On 1/f noise," Mathematical Problems in Engineering , vol. 2012, 2012.

[28] M. H. Heydari, M. R. Hooshmandasl, F. M. Maalek Ghaini, M. Li, "Chebyshev wavelets method for solution of nonlinear fractional integro-differential equations in a large interval," Advances in Mathematical Physics , vol. 2013, 2013.

[29] S. S. Ray, R. K. Bera, "Analytical solution of a fractional diffusion equation by Adomian decomposition method," Applied Mathematics and Computation , vol. 174, no. 1, pp. 329-336, 2006.

[30] Q. Wang, "Numerical solutions for fractional KdV-Burgers equation by Adomian decomposition method," Applied Mathematics and Computation , vol. 182, no. 2, pp. 1048-1055, 2006.