Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2013 / Article
Special Issue

New Developments in Fluid Mechanics and Its Engineering Applications

View this Special Issue

Research Article | Open Access

Volume 2013 |Article ID 765013 | 7 pages |

Spectral Local Linearisation Approach for Natural Convection Boundary Layer Flow

Academic Editor: Waqar Khan
Received27 Mar 2013
Accepted03 Aug 2013
Published04 Sep 2013


The present work introduces a spectral local linearisation method (SLLM) to solve a natural convection boundary layer flow problem with domain transformation. It is customary to find solutions of semi-infinite interval problems by first truncating the interval and subsequently applying a suitable numerical method. However, this gives rise to increased error terms in the numerical solution. Carrying out a transformation of the semi-infinite interval problems into singular problems posed on a finite interval can avoid the domain truncation error and enables the efficient application of collocation methods. The SLLM is based on linearising and decoupling nonlinear systems of equations into a sequence or subsystems of differential equations which are then solved using spectral collocation methods. A comparative study between the SLLM and existing results in the literature was carried out to validate the results. The method has shown to be a promising efficient tool for nonlinear boundary value problems as it gives converging results after very few iterations.

1. Introduction

Boundary layer problems involve a rapid change in the value of a physical variable over a limited region of time or space and therefore are governed by highly nonlinear equations. They are very common in physics and have numerous applications in extrusion, melt spinning, food processing, mechanical forming processes, and many other fields.

In recent years, a host of different methods have been developed with the aim of finding the most efficient method for solving the often complex nonlinear equation systems governing the fluid flows. Boundary layer problems are very difficult to solve because besides having very thin regions where there is rapid change of the fluid properties they are defined on unbounded domains. Analytical methods are being used to solve boundary layer problems but converge very slowly for some boundary layer problems, particularly those with very large parameters. Traditional numerical techniques such as finite differences, finite element methods, Keller-Box method [1], Runge-Kutta with shooting methods [24], and spectral collocation methods [57] are still being used in solving nonlinear BVPs arising in fluid mechanics and other engineering applications. However, numerical approaches also encounter problems in resolving the solution of the governing equations in the very thin regions and in cases where there are singularities or multiple solutions hence the need for the search for ways of improving numerical methods.

In this paper we present a new numerical technique based on linearising and decoupling nonlinear systems of equations into a sequence of or subsystems of differential equations which are then solved using spectral collocation methods. The spectral local linearisation method (SLLM) is used here to solve the problem of natural convection boundary layer flow that arises in glass-fibre production process. Free convection boundary layer problems arise in many applications including buoyancy-induced motions in the atmosphere, oceans, and quasi-solid bodies such as earth. The problem under study was earlier considered in the works of Kuiken [8, 9], Liao [10], and Motsa et al. [11]. In this work, however, the problem is first transformed from the semi-infinite domain into the finite interval to avoid using a scaling parameter to invoke the boundary condition at infinity. The skin friction and Nusselt number are generated by the SLLM and validated by comparing them with those of Kuiken [8].

2. Problem Statement

Consider the laminar free-convection flow of an incompressible viscous fluid about a flat plate parallel to the direction of the generating body force. Figure 1 presents the physical model of the flow. The basic governing boundary layer equations, following Kuiken [9], Ghotbi et al. [12], and Ali et al. [13], are as follows: with boundary conditions where and are the velocity components along the and -axes, respectively, is the kinematic viscosity of the fluid, is the acceleration due to gravity, is the coefficient of thermal expansion, is the fluid temperature, is the temperature of the surrounding fluid, and is the thermal diffusivity.

Applying the transformations below following [9, 12, 13] where is the stream function defined as and are similarity functions, (2) and (3) become along with the boundary conditions

Prime denotes differentiation with respect to the similarity variable and is the Prandtl number. It was reported in [9] that both and decay algebraically as . Equations (7) can be simplified by using the following transformations [10]: where is a spatial scale parameter. Substituting (9) in (7) and (8) results in with boundary conditions where the prime now denotes differentiation with respect to .

3. Domain Transformation

In this section we transform the problem from the semi-infinite domain into the finite domain . The transformation procedure proposed by Hammerling et al. [14] transforms a general BVP of the form into the domain by change of the independent variable. A transformation of the form transforms the semi-infinite interval into . For the general case, the function and its derivatives are then expressed as

The boundary conditions are transformed in a similar manner with an exception of a problem with Neumann boundary conditions at the free end. In that case the problem has to be reduced to have admissible boundary conditions [14]. Previous studies carried out by Budd et al. [15] and Kitzhofer et al. [16] showed that transforming a problem into a finite domain before it is being solved facilitates the efficient use of collocation methods, and classical convergence orders can be observed.

Using (13), (10) together with boundary conditions (11) is transformed into

4. Spectral Local Linearisation Method (SLLM)

The transformed system of (14) and (15) is solved in this section using the spectral local linearisation method. The SLLM is based on developing a decoupled iterative scheme that is then chronologically solved using spectral methods. Local linearisation is applied to the equations before generating the iterative scheme in a manner similar to the Gauss-Seidel approach of decoupling linear algebraic systems. Consider a system of nonlinear ordinary differential equations in unknown functions , written as a sum of its linear and nonlinear terms as follows:

Let the subscripts and represent the previous and current iteration, respectively. Local linearisation of at the previous iteration about is then carried out as follows: so that, at the current iteration, (16) becomes where represents . Using the Gauss-Seidel approach, whereby the updated solutions , obtained at previous equations, are used to find a solution at the current iteration level , the local linearisation iteration scheme can be generated as follows: where at the th iteration represents . Hence, starting from an initial approximation , the iterative scheme (19) is solved until convergence is reached for all unknowns.

Following the guidelines above, we generate the iterative scheme for (14) with the corresponding boundary conditions as follows:

Since the coefficient parameters and the right hand side of (20) are known (from previous iterations), the equation system (20) can easily be solved using any numerical method such as finite differences, finite elements, Runge-Kutta-based shooting methods, or collocation methods. In this work, the equations are solved using the Chebyshev spectral collocation method. This method is based on approximating the unknown functions by the Chebyshev interpolating polynomials in such a way that they are collocated at the Gauss-Lobatto points defined as where is the number of collocation points used (see, e.g., [17, 18]). In order to implement the method, the physical region is transformed into the region using the following mapping: The unknown functions and are approximated at the collocation points by where is the number of collocation points and is the th Chebyshev polynomial defined as

The derivatives of the variables at the collocation points are represented as where with being the Chebyshev spectral differentiation matrix (see, e.g., [17, 18]). Substituting (22)–(25) in (20) we obtain the following system: where

In the above equations and are the values of functions and , respectively, and is a zero vector.

Equations (26) are then solved chronologically starting from the initial approximate solutions which are chosen to satisfy the boundary conditions (15). These initial guesses were suggested in [10] and are presented here after being transformed into the finite interval. The parameter appearing in (28) is the so-called auxiliary parameter which is used to control convergence of the solution.

The convergence rate of the SLLM algorithm can be significantly improved by applying the successive overrelaxation (SOR) technique on (20). Under the SOR framework, a convergence controlling relaxation parameter is introduced and the SLLM scheme for finding, say , is modified to

The results in the following section show that, for , applying the SOR method improves the efficiency and accuracy of the SLLM.

5. Results and Discussion

Numerical computations for the solution of the governing equations (14) and (15) were carried out using the spectral local linearization method. Values of the skin friction and the Nusselt number were generated at different Prandtl numbers. Results for the typical velocity and temperature profiles were also presented. In order to measure the accuracy and convergence of the SLLM, the present results were verified by comparing them with those obtained by Kuiken [8]. The plots were generated using 50 collocation points unless specified.

As pointed out in Section 4 that is a convergence controlling parameter, a proper selection of is necessary. To choose proper values of we use the norm of residual function of () and that of () where and are plotted as functions of :

The valid region of where the solution converges becomes the horizontal part of the residual curve. Numerical iterations are carried out until a representative solution residual falls below a set tolerance. In this work the tolerance level was set to be . Illustrative residual curves are shown in Figures 2 and 3. For example, in Figure 2, the valid values of that will give converging results for are and about for .

Figures 4 and 5 show the velocity and temperature profiles, respectively, when the Prandtl number is varied. The dimensionless velocity profiles are observed to decrease with increase in in Figure 4. Because of the nature of the boundary conditions (slip boundary condition not considered), the velocity profile is seen to increase to reach a particular maximum point and then reduces to zero as . Similarly, as increases, the temperature profiles also decrease as can be seen from Figure 5. For larger values of the Prandtl number, the thermal boundary layer thickness tends to be shorter because it is when the effective thermal diffusivity is low for a fixed kinematic viscosity. The results are consistent with results in the literature (see, e.g., [10, 12]).

Tables 1 and 2 present a comparison of the SLLM solutions of the skin friction and the Nusselt number for the flow at different values of the Prandtl number against those of Kuiken [8]. The effect of the Prandtl number on the skin friction coefficient and the Nusselt number is further shown in Figures 6 and 7. Convergence to the numerical solutions obtained by Kuiken [8] is achieved for up to eight decimal places. Values of the skin friction are accurate for all values of while the Nusselt number is accurate for up to six decimal places in some cases. The skin friction decreased with increase in the Prandtl numbers while the heat transfer rates increased with increased Prandtl numbers. The tables also reveal the advantage of relaxation on the method. Convergence rates are improved remarkably for each case. Worth noting also is the advantage of transforming the domain into the finite interval. There was no need to search for a suitable scaling parameter to invoke the boundary conditions at infinity while a minimal number of collocation points were used to generate the results.

Pr Basic SLLM SLLM with SOR ( ) Reference [8]
Iter. Iter.

0.001400.1391.12313813 0.85 19 1.12313813 1.12313813
0.01500.143 1.063380860.90 19 1.06338086 1.06338086
0.1 300.1370.92408304 0.85 19 0.92408304 0.92408304
1 100.1250.69321163 0.90 13 0.69321163 0.69321163
10 51360.44711652 0.8513 0.44711652 0.44711652
100 51460.264523540.85 14 0.26452354 0.26452354
1000 51490.151290200.85 15 0.15129020 0.15129020
100005150 0.08554085 0.85 160.08554085 0.08554085

Pr Basic SLLM SLLM with SOR ( ) Reference [8]
Iter. Iter.

0.001 300.1 36 0.04680746 0.85 24 0.04680746 0.04680746
0.01 300.1 41 0.13576074 0.85 19 0.135760740.13576074
0.1 300.5 37 0.35005967 0.90 18 0.35005967 0.35005967
1 5 1 27 0.76986120 0.90 13 0.76986119 0.76986120
10 5 1 36 1.49709921 0.85 13 1.49709921 1.49709921
100 5 1 46 2.74688550 0.85 14 2.74688549 2.74688550
1000 5 1 49 4.93494763 0.85 15 4.93494756 4.93494763
10000 51 50 8.80444927 0.85 16 8.80444960 8.80444927

6. Conclusion

The problem of natural convection boundary layer flow with heat transfer has been considered in the current study. The coupled nonlinear system of differential equations was transformed into a finite domain and solved using the spectral local linearisation method (SLLM). The convergence rate of the SLLM algorithm was improved by applying the successive overrelaxation (SOR) technique. Solving the problem in a finite domain avoided the errors that come with approximating the scaling parameter for the boundary conditions at infinity leading to more accurate results. In this work the efficiency of domain transformation on spectral methods was revealed as accurate results were achieved while using relatively small numbers of collocation points. The flow velocity, temperature profiles, skin friction, and rate of surface heat transfer were successfully computed. Both the velocity and temperature of the flow were seen to be reduced by increased Prandtl numbers. The effect of the Prandtl number was seen to reduce the skin friction while increasing the rate of heat transfer. The SLLM results were also found to be consistent with results available in the literature. The method gave converging results for very large Prandtl numbers. While quite simple and easy to implement, the method has produced accurate results at quite remarkable convergence rates, thus displaying its potential as an improved numerical method for solving nonlinear boundary value problems.


  1. T. Cebeci and P. Bradshaw, Physical and Computational Aspects of Convective Heat Transfer, Springer, New York, NY, USA, 1984. View at: Zentralblatt MATH | MathSciNet
  2. M. S. Alam, M. M. Rahman, and M. A. Sattar, “Similarity solutions for hydromagnetic free convective heat and mass transfer flow along a semi-infinite permeable inclined flat plate with heat generation and thermophoresis,” Nonlinear Analysis: Modelling and Control, vol. 12, no. 4, pp. 433–445, 2007. View at: Google Scholar | Zentralblatt MATH | MathSciNet
  3. K. V. Prasad, M. S. Abel, and S. K. Khan, “Momentum and heat transfer in visco-elastic fluid flow in a porous medium over a non-isothermal stretching sheet,” International Journal of Numerical Methods for Heat and Fluid Flow, vol. 10, no. 8, pp. 786–801, 2000. View at: Publisher Site | Google Scholar
  4. M. Yürüsoy, “Similarity solutions of boundary layer equations for a special non-Newtonian fluid in a special coordinate system,” Applied Mathematics and Mechanics, vol. 25, no. 5, pp. 587–594, 2004. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  5. E. M. E. Elbarbary and N. S. Elgazery, “Chebyshev finite difference method for the effect of variable viscosity on magneto-micropolar fluid flow with radiation,” International Communications in Heat and Mass Transfer, vol. 31, no. 3, pp. 409–419, 2004. View at: Publisher Site | Google Scholar
  6. E. M. E. Elbarbary and N. S. Elgazery, “Chebyshev finite difference method for the effects of variable viscosity and variable thermal conductivity on heat transfer from moving surfaces with radiation,” International Journal of Thermal Sciences, vol. 43, no. 9, pp. 889–899, 2004. View at: Publisher Site | Google Scholar
  7. N. S. Elgazery, “The effects of chemical reaction, Hall and ion-slip currents on MHD flow with temperature dependent viscosity and thermal diffusivity,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 4, pp. 1267–1283, 2009. View at: Publisher Site | Google Scholar
  8. H. K. Kuiken, “On boundary layers in fluid mechanics that decay algebraically along stretches of wall that are not vanishingly small,” IMA Journal of Applied Mathematics, vol. 27, no. 4, pp. 387–405, 1981. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  9. H. K. Kuiken, “A “backward” free-convective boundary layer,” View at: Google Scholar
  10. S. Liao, Beyond Perturbation: Introduction to Homotopy Analysis Method, CRC Series: Modern Mechanics and Mathematics, Chapman & Hall/CRC Press, Boca Raton, Fla, USA, 2003. View at: MathSciNet
  11. S. S. Motsa, S. Shateyi, and Z. Makukula, “Homotopy analysis of free convection boundary layer flow with heat and mass transfer,” Chemical Engineering Communications, vol. 198, no. 6, pp. 783–795, 2011. View at: Publisher Site | Google Scholar
  12. 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
  13. F. M. Ali, R. Nazar, and N. M. Arifin, “Numerical investigation of free convective boundary layer in a viscous fluid,” The American Journal of Scientific Research, no. 5, pp. 13–19, 2009. View at: Google Scholar
  14. R. Hammerling, O. Koch, Ch. Simon, and E. Weinmuller, “Numerical solution of singular eigenvalue problems for ODEs with a focus on problems posed on semi-infinite intervals,” 2012, View at: Google Scholar
  15. C. Budd, O. Koch, and E. Weinmuller, “Computation of self-similar solution profiles for the nonlinear Schrodinger equation,” 2012, View at: Google Scholar
  16. G. Kitzhofer, O. Koch, G. Pulverer, C. Simon, and E. B. Weinmüller, “BVPSUITE, a new MATLAB solver for singular implicit boundary value problems,” ASC Report No. 35, 2009. View at: Google Scholar
  17. C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang, Spectral Methods in Fluid Dynamics, Springer Series in Computational Physics, Springer, New York, NY, USA, 1988. View at: MathSciNet
  18. L. N. Trefethen, Spectral Methods in MATLAB, vol. 10 of Software, Environments, and Tools, SIAM, Philadelphia, Pa, USA, 2000. View at: Publisher Site | MathSciNet

Copyright © 2013 S. S. Motsa 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.

776 Views | 480 Downloads | 11 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder