Research Article  Open Access
Solution of the Boundary Layer Equation of the PowerLaw Pseudoplastic Fluid Using Differential Transform Method
Abstract
The boundary layer equation of the pseudoplastic fluid over a flat plate is considered. This equation is a boundary value problem (BVP) with the high nonlinearity and a boundary condition at infinity. To solve such problems, powerful numerical techniques are usually used. Here, through using differential transform method (DTM), the BVP is replaced by two initial value problems (IVP) and then multistep differential transform method (MDTM) is applied to solve them. The differential equation and its boundary conditions are transformed to a set of algebraic equations, and the Taylor series of solution is calculated in every sub domain. In this solution, there is no need for restrictive assumptions or linearization. Finally, DTM results are compared with the numerical solution of the problem, and a good accuracy of the proposed method is observed.
1. Introduction
Boundary layer equation of the Newtonian fluid over a flat plate is one of the classic issues of mechanical engineering. Because of high level of application of nonNewtonian fluids in industry, the boundary layer equation of the nonNewtonian fluids is being noted by engineers. On the other hand, nonlinearity of these equations attracted the interest of mathematicians to evaluate the power and accuracy of the approximate and numerical methods. One of the most significant non Newtonian models is the power law one. In this model, shear thinning or shear thickening property of the fluid is considered and the shear stress has a nonlinear relation with shear rate, while shear stress has a linear relation with shear rate for the Newtonian fluid.
As we pointed out above, boundary layer problems are among the nonlinear ones and most of them do not have an exact analytical solution. So, numerical and approximate methods are used by researchers to solve such equations. The most known numerical method used to solve boundary layer problems is the shooting one. Based on this method, boundary value problem is transformed to an initial value problem with unknown initial values. After that, the problem is replaced with a system of firstorder ordinary differential equations and usually is solved through RungeKutta method. Shooting method is appropriate to solve many boundary value problems, but this one is not useful for solving some BVPs because of the instability of the solution.
Approximate techniques like decomposition method (DM), homotopy analysis method (HAM), homotopy perturbation method (HPM), and variational iteration method (VIM) are good substitutes for numerical methods. During the recent years, boundary layer problems have been solved using some of these methods, such as HAM [1–9], HPM [10–14], VIM [15–17], and DM [18–21]. In most of the researches, some modifications were introduced to overcome the nonlinearity and the boundary condition at infinity.
Differential transform method (DTM) is also one of the other approximate methods to solve differential equations. Here, DTM and multistep DTM are used to solve the boundary layer equation of the pseudoplastic fluid. This method was introduced by Zhou [22] to solve initial value problems in analysis of the electrical circuits. After that, DTM is applied to differential algebraic equation [23, 24], partial differential equation [25–30], integral equation [31–33], ordinary differential equation [34–38], and fractional differential equation [39–42]. The method is an iterative technique to find the Taylor series solution of the problem. There is no need for the high calculation cost to determine the coefficients of Taylor series, which is the reverse of the standard Taylor series method.
This paper is organized as follows. In Section 2, DTM and multistep DTM are introduced. In Section 3, the boundary layer equation of the pseudoplastic fluid is extracted from continuity and momentum equations as presented in [43]. In Section 4, the problem is solved and the results are illustrated as some figures and tables.
2. Differential Transform Method
The differential transform is defined as follows: whereis an arbitrary function andis the transformed function. The inverse transformation is as follows: Substituting (1) into (2), we have The functionis usually considered as a series with limited terms and (2) can be rewritten as whererepresents the number of Taylor series’ components. Usually, through elevating this value, we can increase the accuracy of the solution.
Although the DTM series solution is a good approximate of the exact solution, but this series is diverged for greater areas. There are two ways to overcome this problem. One is the using pade approximate technique. Usually, Pade approximate gives us more exact information about the behavior of the solution (see [44]). Another one is to use multistep DTM. Based on this one, solution domain is divided to some sub domains.
To solve a differential equation in the domainusing multistep DTM, this domain is divided tosections. We suppose the sub domains are equal and length of each sub domain is. So, there is a separate function for every sub domain as follows: where. Multistep DTM for every sub domain is defined as The inverse multistep DTM is Some of the properties of DTM and multistep DTM are shown in Table 1. These properties are extracted from (1) and (6).

3. Mathematical Formulation
Twodimensional boundary layer equations for an incompressible fluid are whereis the fluid density,is the shear stress, andandare the velocities in direction ofand, respectively. For the nonNewtonian powerlaw fluid, the shear stress is calculated through the following relation: whereandare the flow consistency and powerlaw index, respectively.
Here, the following dimensionless parameters are used: whereis the plate length,is the far field velocity, andis generalized Reynolds number.
Combining (8), (10), and (11), the dimensionless momentum equation can be obtained: with the boundary conditions Similarity variable and stream function for the problem is defined as [45] Regarding (14), the velocities can be obtained as a function of the similarity variable: Substituting (15) in (12), we have Equation (16) is the boundary layer equation of the powerlaw fluid as a function of index. For pseudoplastic fluid () the value of theis more than zero [43, 45–48]. So, (16) can be simplified as
4. Solution of the Problem
In this section, we try to solve (17) using DTM and multistep DTM. The solution consists of two stages; first through mathematical relations and applying DTM, the value ofis found. After that, boundary layer equation of the pseudoplastic fluid is solved as an initial value problem (IVP) using multistep DTM.
4.1. Applying DTM (Obtaining)
The BVP (17) can be transformed to an initial value problem with the replacement of the unknown initial conditioninstead of the boundary condition at infinity: where By applying the DTM on (18) at, the recursive relation is derived: Regarding (19), Taking differential transform on (21), we have From (19), the value ofcan be obtained: Combining (20), (22), and (23), the coefficients of the Taylor series solution of (18) can be calculated: Now, consider a similar equation (18) in which: In the same manner, the coefficientsare Substituting (24) and (26) in (2), theandare obtained. These functions have a relation as follows: For, both values ofandare 1/3 for the Blasius equation. From (27), we have Regarding (28) and the boundary condition at infinity in (17), the value ofcan be obtained:
4.2. Applying MultiStep DTM (Solution of IVPs)
As was noted, multistep DTM usually is used for solving problems in which Taylor series is diverged in the solution domain. This technique was already used to solve engineering and computational problems [49–57].
To solve boundary layer problems, the domainis replaced by. Butshould be great enough that the solution is not to be dependent on. The solution domain should be divided toequal parts (). So, we have where Applying multistep DTM on (30) and (31), we have So And for The initial conditions of the problem are considered for first sub domain (). Regarding (6), (25), and (31), initial conditions transformed as follows: The boundary conditions of each sub domain are the continuity of the,, and. These boundary conditions can be obtained from (7): And we have for The value of thecan be calculated by differentiating from (7): The unknown parameter () is calculated from (29). Now, (18) is solved with a similar process like (25) using multistep DTM. The only difference is that the conditionis replaced by the condition.
In Figure 1 is illustrated for different values ofand compared with the numerical solution and results of Lemieux et al. [43]. The numerical solution of the problem is done by the MATLAB software. Figures 2 and 3 show the variation ofand, respectively. The most important step of this scheme is to choose the appreciate finite value of. Thus to estimate this value, we start with an initial guess forand solve the problem. The solution process is repeated with another larger value ofuntil the values ofanddiffer regarding the desired accuracy.
In Table 2 the values of theare compared with the numerical solution and results of Lemieux et al. [43]. The approximate solution of the problem is presented in Table 3 for a special case ().


5. Conclusion
In this paper, the boundary layer equation of the pseudoplastic fluid over a flat plate was extracted from the boundary layer theory as presented in [43]. This equation is a boundary value problem with high nonlinearity and boundary condition at infinity. Using DTM, the BVP transformed to a pair of initial value problems and multistep DTM is used to solve them.
Using the method, the differential equation and boundary conditions are transformed into a recurrence set of equations. Finally, the coefficients of power series are obtained based on the solution of this set of equations. The results validated with the numerical solution, and a good accuracy is observed. The proposed method overcame nonlinearity without using linearization or restrictive assumptions.
References
 S.J. Liao, “A general approach to get series solution of nonsimilarity boundarylayer flows,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 5, pp. 2144–2159, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. R. Ghotbi, H. Bararnia, G. Domairry, and A. Barari, “Investigation of a powerful analytical method into natural convection boundary layer flow,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 5, pp. 2222–2228, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Z. Ziabakhsh and G. Domairry, “Analytic solution of natural convection flow of a nonNewtonian fluid between two vertical flat plates using homotopy analysis method,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 5, pp. 1868–1880, 2009. View at: Publisher Site  Google Scholar
 S. Abbasbandy and T. Hayat, “Solution of the MHD FalknerSkan flow by homotopy analysis method,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 910, pp. 3591–3598, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. Cheng, S. Liao, R. N. Mohapatra, and K. Vajravelu, “Series solutions of nano boundary layer flows by means of the homotopy analysis method,” Journal of Mathematical Analysis and Applications, vol. 343, no. 1, pp. 233–245, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 S.J. Liao and I. Pop, “Explicit analytic solution for similarity boundary layer equations,” International Journal of Heat and Mass Transfer, vol. 47, no. 1, pp. 75–85, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 S. Liao, “A new branch of solutions of boundarylayer flows over an impermeable stretched plate,” International Journal of Heat and Mass Transfer, vol. 48, no. 12, pp. 2529–2539, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 T. Hayat and M. Sajid, “Homotopy analysis of MHD boundary layer flow of an upperconvected Maxwell fluid,” International Journal of Engineering Science, vol. 45, no. 2–8, pp. 393–401, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. R. Ghotbi, “Homotopy analysis method for solving the MHD flow over a nonlinear stretching sheet,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 6, pp. 2653–2663, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 B. Raftari and A. Yildirim, “The application of homotopy perturbation method for MHD flows of UCM fluids above porous stretching sheets,” Computers and Mathematics with Applications, vol. 59, no. 10, pp. 3328–3337, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. Esmaeilpour and D. D. Ganji, “Application of He's homotopy perturbation method to boundary layer flow and convection heat transfer over a flat plate,” Physics Letters A, vol. 372, no. 1, pp. 33–38, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 L. Xu, “He's homotopy perturbation method for a boundary layer equation in unbounded domain,” Computers and Mathematics with Applications, vol. 54, no. 78, pp. 1067–1070, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. Fathizadeh and F. Rashidi, “Boundary layer convective heat transfer with pressure gradient using Homotopy Perturbation Method (HPM) over a flat plate,” Chaos, Solitons and Fractals, vol. 42, no. 4, pp. 2413–2419, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 H. Bararnia, E. Ghasemi, S. Soleimani, A. R. Ghotbi, and D. D. Ganji, “Solution of the FalknerSkan wedge flow by HPMPade' method,” Advances in Engineering Software, vol. 43, no. 1, pp. 44–52, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 H. Mirgolbabaei, D. D. Ganji, M. M. Etghani, and A. Sobati, “Adapted variational iteration method and axisymmetric flow over a stretching sheet,” World Journal of Modelling and Simulation, vol. 5, no. 4, pp. 307–314, 2009. View at: Google Scholar
 S. T. MohyudDin, A. Yildirim, S. Anl Sezer, and M. Usman, “Modified variational iteration method for freeconvective boundarylayer equation using padé approximation,” Mathematical Problems in Engineering, vol. 2010, Article ID 318298, 11 pages, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A.M. Wazwaz, “The variational iteration method for solving two forms of Blasius equation on a halfinfinite domain,” Applied Mathematics and Computation, vol. 188, no. 1, pp. 485–491, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A.M. Wazwaz, “The modified decomposition method and Padé approximants for a boundary layer equation in unbounded domain,” Applied Mathematics and Computation, vol. 177, no. 2, pp. 737–744, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 E. Alizadeh, K. Sedighi, M. Farhadi, and H. R. EbrahimiKebria, “Analytical approximate solution of the cooling problem by Adomian decomposition method,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 2, pp. 462–472, 2009. View at: Publisher Site  Google Scholar
 S. A. Kechil and I. Hashim, “Nonperturbative solution of freeconvective boundarylayer equation by Adomian decomposition method,” Physics Letters A, vol. 363, no. 12, pp. 110–114, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 E. Alizadeh, M. Farhadi, K. Sedighi, H. R. EbrahimiKebria, and A. Ghafourian, “Solution of the FalknerSkan equation for wedge by Adomian Decomposition Method,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 3, pp. 724–733, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. K. Zhou, Differential Transformation and Its Applications for Electrical CircuIts, Huazhong University Press, Wuhan, China, 1986, (Chinese).
 F. Ayaz, “Applications of differential transform method to differentialalgebraic equations,” Applied Mathematics and Computation, vol. 152, no. 3, pp. 649–657, 2004. View at: Publisher Site  Google Scholar
 H. Liu and Y. Song, “Differential transform method applied to high index differentialalgebraic equations,” Applied Mathematics and Computation, vol. 184, no. 2, pp. 748–753, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 F. Ayaz, “Solutions of the system of differential equations by differential transform method,” Applied Mathematics and Computation, vol. 147, no. 2, pp. 547–567, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. S. V. Ravi Kanth and K. Aruna, “Differential transform method for solving linear and nonlinear systems of partial differential equations,” Physics Letters A, vol. 372, no. 46, pp. 6896–6898, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 F. Ayaz, “On the twodimensional differential transform method,” Applied Mathematics and Computation, vol. 143, no. 23, pp. 361–374, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 X. Yang, Y. Liu, and S. Bai, “A numerical solution of secondorder linear partial differential equations by differential transform,” Applied Mathematics and Computation, vol. 173, no. 2, pp. 792–802, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 S.H. Chang and I.L. Chang, “A new algorithm for calculating twodimensional differential transform of nonlinear functions,” Applied Mathematics and Computation, vol. 215, no. 7, pp. 2486–2494, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 B. Jang, “Solving linear and nonlinear initial value problems by the projected differential transform method,” Computer Physics Communications, vol. 181, no. 5, pp. 848–854, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Z. M. Odibat, “Differential transform method for solving Volterra integral equation with separable kernels,” Mathematical and Computer Modelling, vol. 48, no. 78, pp. 1144–1149, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. Arikoglu and I. Ozkol, “Solutions of integral and integrodifferential equation systems by using differential transform method,” Computers and Mathematics with Applications, vol. 56, no. 9, pp. 2411–2417, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. Arikoglu and I. Ozkol, “Solution of boundary value problems for integrodifferential equations by using differential transform method,” Applied Mathematics and Computation, vol. 168, no. 2, pp. 1145–1158, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 C.L. Chen and Y.C. Liu, “Differential transformation technique for steady nonlinear heat conduction problems,” Applied Mathematics and Computation, vol. 95, no. 23, pp. 155–164, 1998. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 S. Mosayebidorcheh and T. Mosayebidorcheh, “Series solution of convective radiative conduction equation of the nonlinear fin with temperature dependent thermal conductivity,” In press. View at: Publisher Site  Google Scholar
 S.U. Islam, S. Haq, and J. Ali, “Numerical solution of special 12thorder boundary value problems using differential transform method,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 4, pp. 1132–1138, 2009. View at: Publisher Site  Google Scholar
 M. Torabi and A. Aziz, “Convectiveradiative fins with simultaneous variation of thermal conductivity, heat transfer coefficient, and surface emissivity with temperature,” Heat Transfer, vol. 41, no. 2, pp. 99–113, 2012. View at: Publisher Site  Google Scholar
 A. A. Joneidi, D. D. Ganji, and 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. View at: Publisher Site  Google Scholar
 D. Nazari and S. Shahmorad, “Application of the fractional differential transform method to fractionalorder integrodifferential equations with nonlocal boundary conditions,” Journal of Computational and Applied Mathematics, vol. 234, no. 3, pp. 883–891, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Z. Odibat, S. Momani, and V. S. Erturk, “Generalized differential transform method: application to differential equations of fractional order,” Applied Mathematics and Computation, vol. 197, no. 2, pp. 467–477, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 V. S. Erturk, S. Momani, and Z. Odibat, “Application of generalized differential transform method to multiorder fractional differential equations,” Communications in Nonlinear Science and Numerical Simulation, vol. 13, no. 8, pp. 1642–1654, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. Arikoglu and I. Ozkol, “Solution of fractional differential equations by using differential transform method,” Chaos, Solitons and Fractals, vol. 34, no. 5, pp. 1473–1481, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 P. F. Lemieux, R. N. Dubey, and T. E. Unny, “Variational method for a pseudoplastic fluid in a laminar boundary layer over a flat plate,” Journal of Applied Mechanics, Transactions ASME, pp. 345–349, 1971. View at: Google Scholar
 G. A. Baker, Essentials of Pade Approximants, Academic Press, New York, NY, USA, 1975.
 T. G. Myers, “An approximate solution method for boundary layer flow of a power law fluid over a flat plate,” International Journal of Heat and Mass Transfer, vol. 53, no. 1112, pp. 2337–2346, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. Benlahsen, M. Guedda, and R. Kersner, “The Generalized Blasius equation revisited,” Mathematical and Computer Modelling, vol. 47, no. 910, pp. 1063–1076, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. P. Denier and P. P. Dabrowski, “On the boundarylayer equations for powerlaw fluids,” Proceedings of the Royal Society A, vol. 460, no. 2051, pp. 3143–3158, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 C.C. Hsu, “A simple solution for boundary layer flow of power law fluids past a semiinfinite flat plate,” AIChE Journal, vol. 15, no. 3, pp. 367–370, 1969. View at: Publisher Site  Google Scholar
 H. Yaghoobi and M. Torabi, “The application of differential transformation method to nonlinear equations arising in heat transfer,” International Communications in Heat and Mass Transfer, vol. 38, no. 6, pp. 815–820, 2011. View at: Publisher Site  Google Scholar
 M. M. Rashidi, A. J. Chamkha, and M. Keimanesh, “Application of multistep differential transform method on flow of a secondgrade fluid over a stretching or shrinking sheet,” American Journal of Computational Mathematics, vol. 6, pp. 119–128, 2011. View at: Google Scholar
 M. Keimanesh, M. M. Rashidi, A. J. Chamkha, and R. Jafari, “Study of a third grade nonNewtonian fluid flow between two parallel plates using the multistep differential transform method,” Computers and Mathematics with Applications, vol. 62, no. 8, pp. 2871–2891, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 L.T. Yu and C.K. Chen, “The solution of the Blasius equation by the differential transformation method,” Mathematical and Computer Modelling, vol. 28, no. 1, pp. 101–111, 1998. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 B.L. Kuo, “Application of the differential transformation method to the solutions of FalknerSkan wedge flow,” Acta Mechanica, vol. 164, no. 34, pp. 161–174, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Z. M. Odibat, C. Bertelle, M. A. AzizAlaoui, and G. H. E. Duchamp, “A multistep differential transform method and application to nonchaotic or chaotic systems,” Computers and Mathematics with Applications, vol. 59, no. 4, pp. 1462–1472, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 H.S. Peng and C.L. Chen, “Hybrid differential transformation and finite difference method to annular fin with temperaturedependent thermal conductivity,” International Journal of Heat and Mass Transfer, vol. 54, no. 1112, pp. 2427–2433, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 Y.L. Yeh, C. C. Wang, and M.J. Jang, “Using finite difference and differential transformation method to analyze of large deflections of orthotropic rectangular plate problem,” Applied Mathematics and Computation, vol. 190, no. 2, pp. 1146–1156, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 H. P. Chu and L. Ch. Chen, “Hybrid differential transform and finite difference method to solve the nonlinear heat conduction problem,” Communications in Nonlinear Science and Numerical Simulation, vol. 13, pp. 1605–1614, 2008. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 Sobhan Mosayebidorcheh. 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.