American Journal of Computational Mathematics
Vol.06 No.01(2016), Article ID:65175,12 pages

MHD Stagnation Point Flow of a Casson Fluid over a Nonlinearly Stretching Sheet with Viscous Dissipation

Monica Medikare1, Sucharitha Joga1, Kishore Kumar Chidem2

1Department of Mathematics, Osmania University, Hyderabad, Telangana

2Department of Mathematics, Nizam College, Osmania University, Hyderabad, Telangana

Copyright © 2016 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

Received 6 February 2016; accepted 27 March 2016; published 30 March 2016


Study to analyze the MHD stagnation point flow of a Casson fluid over a nonlinearly stretching sheet with viscous dissipation was carried out. The partial differential equations governing this phenomenon were transformed into coupled nonlinear ordinary differential equations with suitable similarity transformations. These equations were then solved by finite difference technique known as Keller Box method. The various parameters such as Prandtl number (Pr), Eckert number (Ec), Magnetic parameter (M), Casson parameter (β) and non linear stretching parameter (n) determining the velocity and temperature distributions, the local Skin friction coefficient and the local Nusselt number governing such a flow were also analyzed. On analysis it was found that the Casson fluid parameter (β) decreased both the fluid velocity and temperature whereas an increase in (β) increased both the heat transfer rate and wall skin-friction coefficient.


Stagnation Point, Casson Fluid, MHD, Viscous Dissipation, Nonlinear Stretching Sheet

1. Introduction

The study of boundary layer flow over a stretching sheet has been a great challenge for researchers and it has immense applications in various industrial processes such as extraction of polymer sheet, paper production, wire drawing, glass-fiber production etc. The flow of incompressible fluid over a linearly stretching sheet was first investigated by Crane [1] . The Mass Transfer with suction or blowing over a stretching sheet was investigated by Gupta and Gupta [2] . The heat transfer over a stretching sheet of a hydromagnetic flow has been studied by Chakrabarthi and Gupta [3] . The flow over a stretching surface in three diemsional was discussed by Wang [4] . The MHD flow characteristics over a stretching sheet of a viscoelastic fluid was demonstrated by Andersson [5] . Later his work was extended by M. I. Char [6] with mass transfer.

The numerical analysis of magnetic field on Eyring-Powell fluid flow towards a stretching sheet has been discussed by N S Akbar [7] .

All the above investigations are restricted for the flow over a linearly stretching sheet but it is not necessary that the stretching sheet has to be linear. Vajravelu [8] analyzed the viscous flow over a nonlinearly stretching sheet. The heat transfer characteristics of a viscous fluid over a non linear stretching sheet was later introduced by Cortell [9] . A great number of studies for the boundary layer flow over a non linear stretching sheet under different aspects of heat and mass transfer, slip and convective boundary conditions etc., are presented in the literature [10] - [16] .

The study of stagnation point flow towards a stationary semi infinite wall was first introduced by Hiemenz [17] in a two dimensional case. In his work he reduced the Navier-Stokes equations into non linear ordinary differential equations with the help of similarity transformations. The same problem was extended by Homann [18] . Chiam [19] combined the works of Hiemenz and Crane i.e., the stagnation point flow towards a stretching sheet when the stretching rate of the plate is equal to the strain rate of the stagnation point flow and he found no boundary layer structure near the plate. El-Dab [20] obtained numerical solution of MHD flow of micropolar fluid with heat and mass transfer towards a stagnation point on a vertical plate.

In all these attempts, stagnation point flow due to the stretching sheet was analyzed. The boundary layer flow over a shrinking sheet was first investigated by Wang [21] . He dealt with stagnation point flow on a two dimensional shrinking sheet and axi symmetric stagnation point flow on a axisymmetric shrinking sheet. In contrast, Lok et al. [22] numerically studied non orthogonal stagnation point flow towards a stretching sheet. In his work, he determined that the obliqueness of a free stream line causes the shifting of the stagnation point towards the incoming flow.

Many fluids in industries resemble non-Newtonian behaviour. Non-Newtonian fluids are more appropriate than Newtonian fluids because of their varied industrial applications like petroleum drilling, polymer engineering, certain separation processes, food manufacturing etc.

For non-Newtonian fluids, the relationship between stress and the rate of strain is not linear and it is difficult to express all these properties in a single constitutive equation. Consequently, these fluid models [23] - [31] have been proposed depending on various physical characters.

Casson fluid is one such type of such non-Newtonian fluid, which behaves like an elastic solid, with a yield shear stress existing in the constitutive equation. This fluid model has its origin in modelling of flow of many biological fluids especially blood. Examples of such fluids include foams, yoghurt, molten chocolate, cosmetics, nail polish, tomato puree etc. Casson [32] introduced this model to predict the flow behaviour of pigment oil suspensions of the printing ink type. Later on, several researchers studied Casson fluid pertaining to different flow situations. The unsteady boundary layer flow and heat transfer of a Casson fluid over a moving flat plate with a parallel free stream was studied by Mustafa et al. [33] . The exact solution for boundary layer flow of Casson fluid over a permeable stretching/shrinking sheet with and without external magnetic field was discussed by Bhattacharyya et al. [34] [35] . The Casson fluid has an infinite viscosity at zero rate of shear and a yield stress below which no flow occurs and a zero velocity at an infinite shear rate [36] [37] . An excellent collection of articles can be found in [38] - [40] .

Motivated by the above investigations on non-Newtonian fluids and its wide applications, the objective of the present study is to analyze MHD stagnation point flow of a Casson fluid over a non linear stretching sheet with viscous dissipation. We have extended the works of Vajravelu [8] by considering magnetic field, stagnation point and viscous dissipation effects to the flow regime for more physical implications. In addition to this, the roles of the other parameters are explained in brief. The results are displayed in figures and tables. A parametric study of the physical parameters is conducted and a representative set of numerical results for the velocity and temperature profiles as well as the skin friction coefficient and local Nusselt number is illustrated graphically to show interesting features of the solutions. Here the governing partial differential equations were transformed into coupled non linear ordinary differential equations using appropriate similarity tranformations. The transformed ODE’s were then solved by Keller Box method. However, to the best of our knowledge no work has been carried out yet to analyze the present result. The layout of this paper is as follows.

In Section 2 the problem is formulated and similarity transformation has been employed to transform the partial differential equation into a non linear ordinary differential equations. In Section 3 a numerical solution using Keller box method has been discussed. In Section 4 the numerical results have been discussed graphically. Finally, the main concluding observations are mentioned in Section 5.

2. Mathematical Formulation

Consider the steady two-dimensional stagnation point flow of an incompressible Casson fluid located at y = 0. The flow being confined in the region y > 0. We take the non-linear stretching sheet in the XOY plane (see Figure 1). Along the x-axis, two equal and opposite forces are being applied. The surface is stretched in the x-di- rection such that the x-component of the velocity varies non-linearly along it, i.e. where c > 0 is constant of proportionality and n is a power index. A magnetic field of uniform strength is applied perpendicular to the surface. The magnetic Reynolds number is taken to be small enough so that the induced magnetic field can be neglected in comparison to the applied magnetic field. It is also assumed that the ambient fluid is moved with a velocity, where a > 0 is a constant. The rheological equation of state for an isotropic and incompressible flow of the Casson fluid [Nakamura and Sawada [34] , Mustafa et al. [31] is given by

where and are the plastic dynamic viscosity, yield stress of the fluid respectively. Similarly is the product of the component of deformation rate with itself, , is the (i, j)-th component of the deformation rate and is a critical value of this product based on the non-Newtonian model.

The continuity, momentum and energy equations governing the fluid flow are given by

where u, v are the velocity components in x, y direction respectively, is the kinematic viscosity, is the viscosity, is the Casson fluid parameter, is the thermal diffusivity of the fluid, is thespecific heat and is the magnetic field.

The suitable boundary conditions are given by

Figure 1. Physical sketch of the given problem.

Here, c () is the surface stretching sheet related parameter. are the uniform temperature at the sheet, free stream temperature respectively and n is the power index related to the surface stretching speed.

With the help of following similarity transformations

The Equations (1), (2) and (3) are transformed into coupled non linear ordinary differential equations as follows.

and the boundary conditions are transformed into

where is the Prandtl Number, is the magnetic parameter,

is the Eckert number and is the velocity ratio parameter.

where prime denotes the differentiation with respect to.

The quantities of practical interest are the Skin friction coefficient and the Local Nusselt number which are defined as

where is the shear stress or skin friction along the stretching sheet and is the heat flux from the sheet.

Hence the dimensionless form of Skin friction and the Local Nusselt number are given by

where is the local Reynolds number.

3. Numerical Solution

The numerical solution for the above coupled ordinary differential equations for different values of governing parameters is obtained using finite difference scheme called Keller-box method. This method involves the four main steps which are as follows:

1) Reduce the equation or system of equations to a first order system.

2) Write the difference equations using central differences.

3) Linearize the resulting algebraic equations (if they are non-linear) by Newton’s method.

4) Write them in matrix-vector form and use the block-tridiagonal-elimination technique to solve the linear system.

This method has been widely used in laminar and turbulent boundary layer flows. It seems to be much faster, easier to program, more efficient and flexible to use than other methods.

4. Results and Discussion

The numerical computations are carried out with the help of Keller box method. The effects of different parameters like Stagnation point, MHD, viscous dissipation on velocity and temperature profiles has been clearly analyzed including with its physical quantities of significance. The effect of Casson fluid parameter on skin friction and local Nusselt number has been examined. Apart from those the various quantities of non linear stretching parameter were also analyzed.

To verify the accuracy of the present results, comparison has been made with the previous results of Vajravelu and Cortell (Table 1).

Table 1 clearly shows the values of skin friction coefficient and heat transfer rate at the surface for different values of non linear stretching parameter n which are found to be in excellent agreement.

The effects of magnetic parameter M and velocity ratio parameter λ on the flow field velocity is displayed in Figure 2. It shows that the velocity of the fluid is not affected by various values of M for λ = 1 i.e., when the velocity of ambient fluid and the velocity of the sheet are same, where as for λ = 0.2 an increase in M decreases the velocity boundary layer thickness. However, with an increase in M the flow field velocity increases and boundary layer thickness decreases for λ = 2. Due to the Lorentz force induced by the action of magnetic field ,the thickness of boundary layer decreases for both λ = 0.2 and λ = 2.

Table 1. Values, of with Ec = 0, M = 0.

Figure 2. Velocity profile for different values of M when λ = 0.2, 1, 2 and β = 2, n = 10, Pr = 0.71.

Figure 3 depicts the effect of the velocity ratio parameter on the flow field velocity. It can be observed that for λ < 1 i.e., when the stretching velocity of the sheet exceeds the free stream velocity, the velocity of the fluid and the boundary layer thickness increases with an increase in λ. Moreover, for λ > 1 i.e., when the free stream exceeds the stretching velocity, the flow velocity increases and the boundary layer thickness decreases with an increase in λ. If λ = 1 i.e., when the stretching and free stream velocities are equal, then there is no boundary layer of Casson fluid flow near the sheet.

Figure 4 and Figure 5 indicates the effect of Casson parameter β on velocity and temperature profiles for a non linear stretching parameter n = 10. It can be observed from Figure 4 that the fluid velocity decreases when β increases. Due to the increase of Casson parameter β, the yield stress reduces and hence the momentum boundary layer thickness decreases.

Figure 3. Velocity profile for various values of λ when M = 1, β = 2, n = 10, Pr = 0.71 and Ec = 1.

Figure 4. Velocity profile for various values of β.

Figure 5 shows that the temperature decreases with increase in β implying that the thermal boundary layer decreases. In the vicinity of the sheet the temperature profile curves for larger values of β are lower than for those with smaller values of β. However, as we move far from the sheet within the thermal boundary layer, viscous dissipation vanishes and the curves for larger values of β tend to override for smaller values of β. Hence the temperature profiles intersect each other at a point and the behaviour goes other way round at this point.

Figure 6 and Figure 7 indicates the velocity and the temperature profiles for various values of nonlinear

Figure 5. Temperature profile for various values of β.

Figure 6. Velocity profile for various values of n.

Figure 7. Temperature profile for various values of n.

stretching parameter n. Figure 6 shows that the velocity of the fluid decreases with an increase in power index n, while the temperature increases (Figure 7). From the above two figures no specific variations were observed

when n is large. This is because of the term which approximately approaches to 2 when n reaches infin-

ity, as mentioned in the subfigures. Therefore the observations for the large values of n is not a study of interest.

Figure 8 indicates the effect of magnetic parameter M on velocity and temperature profiles. The velocity of the fluid decreases with an increase in M. This causes retarding effect on the flow field leading to the prominent reduction in velocity due to Lorentz force effect. Therefore, the Lorentz force increases the opposition of the flow of fluid reducing the velocity of the flow.

Figure 9 shows the temperature profile for various values of Pr. It is clear that the dimensionless parameter θ

Figure 8. Velocity profile for various values of M.

Figure 9. Temperature profile for various values of Pr.

decreases with the increase in Prandtl number. Since the Prandtl number is the ratio of momentum diffusivity to thermal diffusivity; it reduces the thermal boundary layer thickness. In general the Prandtl number is used in heat transfer problems to reduce the relative thickening of the momentum and the thermal boundary layers.

The influence of Eckert number Ec for linearly/non linearly stretching parameters i.e., n = 1, n = 10 is depicted in Figure 10. It illustrates that the temperature increases with an increase in Ec. The increase in temperature profile is higher for non linear stretching parameter than that for linear stretching parameter. The viscous dissipation produces heat due to drag between the fluid particles and this extra heat causes an increase of the initial fluid temperature, and these variations are clearly mentioned for linearly (n = 1) and nonlinearly stretching cases (n = 10). Both Figure 10 and Figure 11 reveal that the viscous dissipation may not only cause thermal reversal but also increase the thermal boundary layer whether or.

Figure 10. Temperature profile for various values of Ec, λ = 0.2.

Figure 11. Temperature profile for various values of Ec, λ = 1.5.

Figure 12 and Figure 13 indicates the skin friction coefficient and Local Nusselt number with variations of Casson parameter β respectively. Figure 12 depicts the variation of the skin-friction coefficient with β for two values of non linear stretching parameter n = 1, n = 10. It is found that as β increases the skin-friction coefficient increases. From this figure, it can be observed that the shear stress at the wall is negative. Physically, a negative sign implies that a drag force has been exerted by the surface on the fluid, and a positive sign implies the opposite.

Figure 13 indicates the variation of heat transfer coefficient with Casson parameter β. As β increases, the heat transfer coefficient gradually increases. Wall temperature gradient decreases with increase in n.

Figure 12. Skin friction coefficient for various values of Casson parameter β.

Figure 13. Local Nusselt number for various values of Casson parameter β.

5. Conclusion

In the present study, MHD stagnation point flow of a Casson fluid over a non linearly stretching sheet are investigated with viscous dissipation. The numerical solution is obtained by Keller box technique. The effects of various governing parameters on heat flow characteristics were analyzed. Briefly the above discussion can be summarized as follows.

・ The velocity boundary layer thickness reduces for magnetic parameter M.

・ An increase in Casson parameter β decreases the velocity of the fluid as well as the thermal boundary layer thickness.

・ The velocity of the fluid decreases and temperature increases with an increase in nonlinear stretching parameter..

・ The Eckert number increases the thermal boundary layer thickness where as the Prandtl number decreases it.

・ Both the skin friction coefficient and heat transfer coefficient increases with Casson parameter β.

Cite this paper

Ojinga Gideon Omiunu,1 1,Monica Medikare,Sucharitha Joga,Kishore Kumar Chidem, (2016) MHD Stagnation Point Flow of a Casson Fluid over a Nonlinearly Stretching Sheet with Viscous Dissipation. American Journal of Computational Mathematics,06,37-48. doi: 10.4236/ajcm.2016.61005


  1. 1. Hauser, E. (2003) The WTO Food Security and Poverty Reduction—A Contradiction in Terms? 2nd Edition, Agriculture and Rural Development, Swaziland Printing and Publishing Company Limited, 8-11.

  2. 2. Akande, T. (2003) The Rice Sector in Nigeria. United Nation Crop Project (UNCP) Country Agricultural Project on Trade Liberalization in Agricultural Sector and the Environment, Geneva, 10. In: Reference .

  3. 3. Nigerian Population Census (2006) National Population Commission, Federal Republic of Nigeria.

  4. 4. Umeh, J.C. and Ataborh, E.M. (2006) Efficiency of Rice Farmers in Nigeria: Potentials for Food Security and Poverty Alleviation.

  5. 5. Ogundari, K. (2006) Determinants of profit Efficiency among Small-Scale Rice Farmers in Nigeria: A Profit Function Approach. Paper Presented at International Association of Agricultural Economists Conference, Gold Coast, 12-18 August 2006, 2.

  6. 6. Ayoola, J.B., Kudi, T.M., Dangbegnon, C., Daudu, C.K., Mando, A., Amapu, I.Y., Adeosun, J.O. and Ezui, K.S. (2012) Gender Perspectives of Action Research for Improved Rice Value Chain in Northern Guinea Savanna, Nigeria. Journal of Agricultural Science, 4, 211-218.

  7. 7. Ologbon, O.A.C., Ikheloa, E.E. and Akerele, E.O. (2012) Adoption of “Ofada” Rice Variety and Technical Efficiency of Rice-Based Production Systems in Ogun State, Nigeria. World Journal of Agricultural Sciences, 8, 624-631.

  8. 8. National Economic Empowerment and Development Strategy (NEEDS) (2007) The NEEDS Secretariat, National Planning Commission. Federal Secretariat, Government Press, Abuja, Nigeria.

  9. 9. Balogun, O.L., Akinlade, R.J. and Campbell, O.A. (2012) Impact of Microfinance on Rural Households. Cooperative Finance in Developing Economies, 2012 International Year of Cooperatives. Soma Prints Ltd., Lagos.

  10. 10. Longtau, S.R. (2000) Multi-Agency Partnership in West African Agriculture: A Review and Description of Rice Production System in Nigeria. Monograph, Ecosystem Development Organization, Jos, 47 p.

  11. 11. Rahman, S.A., Gabriel, J. and Marcus, N.D. (2004) Gender Differentials in Labour Contribution and Productivity in Farm Production. Empirical Evidence from Kaduna State of Nigeria. The National Conference on Family, Makurdi, 1-5 March 2004.

  12. 12. Ojo, C.O. (2012) Technical Efficiency of Rural Women Farmers in Borno State, Nigeria. Developing Country Studies, 2, 61-67.

  13. 13. Mahmood, R. and Hanafi, N. (2012) Entrepreneurial Orientation and Business Performance of Women-Owned Small and Medium Enterprises in Malaysia: Competitive Advantage as a Mediator. International Journal of Business and Social Science, 4, 82-90.

  14. 14. Rahman, S.A. (2008) Women’s Involvement in Agriculture in Northern and Southern Kaduna State, Nigeria. Journal of Gender Studies, 17, 17-26.

  15. 15. Ajayi (1998) Small Scale Farming in Nigeria: Problems and Prospect. The National Workshop on Strategies for Increasing the Agricultural Productivity of Nigerian Women, National Productivity Centre (NPC), 21-23 July.

  16. 16. Watson, J. (2003) SME Performance: Does Gender Matter? The Small Enterprise Association of Australia and New Zealand 16th Annual Conference, Ballarat, 28 September-1 October 2003, Paper No. 32, 1.

  17. 17. Ayoola1, J.B., Dangbegnon, C., Daudu, C.K., Mando, A., Kudi, T.M., Amapu, I.Y., Adeosun, J.O. and Ezui, K.S. (2011) Socio-Economic Factors Influencing Rice Production among Male and Female Farmers in Northern Guinea Savanna Nigeria: Lessons for Promoting Gender Equity in Action Research. Agriculture and Biology Journal of North America, 2, 1010.

  18. 18. Olowa, O.W. and Olowa, O.A. (2011) Issues, Problems and Policies in Agricultural Credit: A Review of Agricultural Credit in Nigeria. Bangladesh e-Journal of Sociology, 8, 87-108.

  19. 19. Nsien, G.F. (2000) A Survey on the Use of Micro-Credit in the Empowerment of Nigerian Women. University of Nigeria Research Publications.

  20. 20. FAO (1984) Grasroots Training for Swamp Rice Development. Manual on Swamp Construction. TCP/SIL/2304(T), FAO, Rome.

  21. 21. Fonjong, L.N. and Athanasia, M.F. (2007) The Fortunes and Misfortunes of Women Rice Producers in Ndop, Cameroon and the Implications for Gender Roles. Journal of International Women’s Studies, 8, 133-147.

  22. 22. Sen, A. (1990) Gender and Cooperative Conflict. In: Tinker, I., Ed., Persistent Inequality, Oxford University Press Oxford, 123-148.

  23. 23. Agawal, B. (1997) “Bargaining” and Gender Relations: Within and Beyond the Household. Feminist Economics, 3, 1-51.

  24. 24. Adebayo, O.O. and Adeola, R.G. (2008) Sources and Uses of Agricultural Credit by Small-Scale Farmers in Surulere Local Government Area of Oyo State. Anthropologist, 10, 313-314.

  25. 25. WARDA (2003) Rice in West Africa.

  26. 26. Amaza, P.S., Kwagbe, P.V. and Amos, A.A. (1999) Analysis of Women’s Participation in Agricultural Cooperatives: Case Study of Borno State, Nigeria. Annals of Borno, 15-16, 187-196.

  27. 27. Olagunju, F.I., Fakayode, S.B., Babatunde, R.O. and Ogunwole-Olapade, F. (2013) Gender Analysis of Sweet Potato Production in Osun State, Nigeria. Asian Journal of Agricultural Extension, Economics & Sociology, 2, 1-13.

  28. 28. Dayo, P., Ephraim, N., John, P. and Ayoola, O.O. (2009) Constraints to Increasing Agricultural Productivity in Nigeria: A Review. Nigeria Strategy Support Program (NSSP) Background Paper No. NSSP 006, September 2009.

  29. 29. Ndanitsa, M.A., Musa, S.A. and Umar, I.S. (2011) Impact of Microfinance Institutions’ Micro-Credit Program on Poverty Alleviation through Agricultural Finance Intermediation to Maize Farmers in Niger State, Nigeria. Savannah Journal of Agriculture, 6, 82-92.

  30. 30. Kazeem, M.A. (1991) Economics of Rice Production in Lagos State: A Case Study of Itoikin Irrigation Project, Epe, Lagos State. Unpublished M.Sc. Thesis, Department of Agricultural Economics, University of Agriculture, Abeokuta.

  31. 31. Alfred, S.D. (2002) Socio-Economic Factors Affecting Rice Production by Women Farmers in Ondo State. The Ogun Journal of Agricultural Sciences, 151-165. In: Ologbon, O.A.C., Ikheloa, E.E. and Akerele, E.O. (2012) Adoption of “Ofada” Rice Variety and Technical Efficiency of Rice-Based Production Systems in Ogun State, Nigeria. World Journal of Agricultural Sciences, 8, 624-631.

  32. 32. Rajasekhar, D. (2004) Institutional Finance for Agriculture Analysis at Micro and Macro Levels. Ph.D. Thesis, University of Mysore, Mysore.

  33. 33. Devi, R.U. and Govt, S.R.K. (2012) The Role of Credit Co-Operatives in the Agricultural Development of Andhra Pradesh, India. International Journal of Cooperative Studies, 1, 55-64.

  34. 34. Bratton, M. (1986) Financing Smallholder Production: A Comparison of Individual and Group Credit Scheme in Zimbabwe. Public Administration and Development, 6, 115-132.

  35. 35. Ellinger, P.N. and Barry, P.J. (2004) A Farmer’s Guide to Agricultural Credit; A Publication of the Centre for Farm and Rural Business Finance: A finance-Focused Initiative Providing Information to Farm and Rural Businesses and their Knowledge Providers. University of Illinois, Urbana Champaign.

  36. 36. Adegbite, D.A., Oloruntoba, A.O. and Olaoye, O.J. (2009) Performance Assessment of Ogun State Agricultural and Multi-Purpose Credit Agency (OSAMCA) in Credit Delivery and Operation (2004-2006). Journal of Sustainable Development in Africa, 10, 127-153.

  37. 37. Badiru, I.O. (2010) Review of Small farmer Access to Agricultural Credit in Nigeria. Nigeria Strategy Support Program. Policy Note 25. International Food Policy Research Institute (IFPRI).

  38. 38. Anyiro, C.O. and Oriaku, B.N. (2011) Access to and Investment of Formal Micro Credit by Small Holder Farmers in Abia State, Nigeria. A Case Study of Absu Micro Finance Bank, Uturu. Journal of Agricultural Sciences, 6, 69-76.

  39. 39. Odebiyi, O.C. and Olaoye, O.J. (2012) Small and Medium Scale Aquaculture Enterprises (SMES) Development in Ogun State, Nigeria: The Role of Microfinance Banks. Libyan Agriculture Research Center Journal International, 3, 1-6.

  40. 40. Ashaolu, O.F., Momoh, S., Phillip, B.B. and Tijani, I.A. (2011) Microcredit Effect on Agricultural Productivity: A Comparative Analysis of Rural Farmers in Ogun State, Nigeria. International Journal of Applied Agricultural and Apicultural Research, 7, 23-35.

  41. 41. Mgbada, J.U. (2000) Production of Staple Crops by Rural Women in Enugu and Ebonyi States. Lessons for Enhancing Poverty Alleviation Programmes. Proceedings of the 6th Annual National Conference of the Agricultural Extension Society of Nigeria, 10-12 April 2000, Ibadan, 11.

  42. 42. Ayoola, J.B. (1999) Economic Analysis of Women Entrepreneurship in Low-Income Households of Benue State, Nigeria. Ph.D. Thesis, Department of Agricultural Economics, Abubarkar Tafawa Balewa University, Bauchi.

  43. 43. Rahman, S.A. and Usman, J.I. (2004) Comparative Analysis of Women’s Participation in Agricultural Production in Northern and Southern Kaduna State, Nigeria. Mobilizing Investors for Sustainable Agricultural Research, Development and Production in Nigeria. Proceedings of the 38th Annual Conference of the Agricultural Society of Nigeria, 17-21 October 2004, 105-108.

  44. 44. Nigerian Institute of Social and Economic Research (NISER) (2002) Assessment of the Economic, Social and Environmental Impact of Rice Production in Nigeria within the Trade Libralisation Framework, A Research Paper, 93 p.

  45. 45. Crane, LJ. (1970) Flow Past a Stretching Plate. Zeitschrift für Angewandte Mathematik und Physik, 21, 645-647.

  46. 46. Gupta, P.S. and Gupta, A.S. (1977) Heat and Mass Transfer on a Stretching Sheet with Suction or Blowing. The Canadian Journal of Chemical Engineering, 55, 744-746.

  47. 47. Chakrabarti, A. and Gupta, A.S. (1979) Hydromagnetic Flow and Heat Transfer over a Stretching Sheet. Quarterly of Applied Mathematics, 37, 73-78.

  48. 48. Wang, C.Y. (1984) The Three Dimensional Flow Due to a Stretching Flat Surface. Physics of Fluids, 27, 1915-1917.

  49. 49. Andersson, H.I. (1992) MHD Flow of a Viscoelastic Fluid Past a Stretching Surface. Acta Mechanica, 95, 227-230.

  50. 50. Char, M.I. (1994) Heat and Mass Transfer in a Hydromagnetic Flow of the Viscoelastic Fluid over a Stretching Sheet. Journal of Mathematical Analysis and Applications, 186, 674-689,

  51. 51. Akbar, N.S. and Ebaid, A. (2015) Numerical Analysis of Magnetic Field on Eyring-Powell Fluid Flow towards a Stretching Sheet. Journal of Magnetism and Magnetic Materials, 382, 355-358.

  52. 52. Vajravelu, K. (2001) Viscous Flow over a Nonlinearly Stretching Sheet. Applied Mathematics and Computation, 124, 281-288

  53. 53. Cortell, R. (2007) Viscous Flow and Heat Transfer over a Nonlinearly Stretching Sheet. Applied Mathematics and Computation, 184, 864-873.

  54. 54. Abbas, Z. and Hayat, T. (2011) Stagnation Slip Flow and Heat Transfer over a Nonlinear Stretching Sheet. Numerical Methods for Partial Differential Equations, 27, 302-314.

  55. 55. Hayat, T., Javed, T. and Abbas, Z. (2009) MHD Flow of a Micropolar Fluid near a Stagnation-Point Towards a Non-Linear Stretching Surface. Nonlinear Analysis: Real World Applications, 10, 1514-1526

  56. 56. Mabood, F., Khan, W.A. and Ismail, A.I.M. (2015) MHD Boundary Layer Flow and Heat Transfer of Nanofluids over a Nonlinear Stretching Sheet: A Numerical Study. Journal of Magnetism and Magnetic Materials, 374, 569-576.

  57. 57. Prasad, K.V., Vajravelu, K. and Dattri, P.S. (2010) Mixed Convection Heat Transfer over a Non-Linear Stretching surface with Variable Fluid Properties. International Journal of Non-Linear Mechanics, 45, 320-330.

  58. 58. Rana, P. and Bhargava, R. (2012) Flow and Heat Transfer of a Nanofluid over a Nonlinearly Stretching Sheet: A Numerical Study. Communications in Nonlinear Science and Numerical Simulation, 17, 212-226.

  59. 59. Ghotbi, A.R. (2009) Homotopy Analysis Method for Solving the MHD flow over a Non-Linear Stretching Sheet. Communications in Nonlinear Science and Numerical Simulation, 14, 2653-2663.

  60. 60. Fanga, T. (2014) Magneto-Hydrodynamic Viscous Flow over a Nonlinearly Moving Surface: Closed-Form Solutions. The European Physical Journal Plus, 129, 92.

  61. 61. Hiemenz, K. (1911) Die Grenzschicht an einem in den gleichformigen Flüssigkeitsstrom eingetauchten geraden Kreiszylinder. Dinglers Polytech Journal, 326, 321-410.

  62. 62. Homann, F. (1936) Der Einfluss grosset Zahigkeit bei der Stromung um den Zylinder and um die Kugel. Zeitschrift für Angewandte Mathematik und Mechanik, 16, 153-164.

  63. 63. Chiam, T.C. (1994) Stagnation-Point Flow towards a Stretching Plate. Journal of the Physical Society of Japan, 63, 2443-2444.

  64. 64. El-Dab, N.T., Ghaly, A.Y., Rizkallah, R.R., Ewis, K.M. and Al-Bareda, A.S. (2015) Numerical Solution of MHD Flow of Micropolar Fluid with Heat and Mass Transfer towards a Stagnation Point on a Vertical Plate. American Journal of Computational Mathematics, 5, 158-174.

  65. 65. Wang, C.Y. (2008) Stagnation Point Flow towards a Shrinking Sheet. International Journal of Non-Linear Mechanics, 43, 377-382.

  66. 66. Lok, Y.Y., Amin, N. and Pop, I. (2006) Non-Orthogonal Stagnation Point Flow towards a Stretching Sheet. International Journal of Non-Linear Mechanics, 41, 622-627.

  67. 67. Fox, V.G., Erickson, L.E. and Fan, L.T. (1969) The Laminar Boundary Layer on a Moving Continuous Flat Sheet Immersed in a Non-Newtonian Fluid. AIChE Journal, 15, 327-333.

  68. 68. Wilkinson, W. (1970) The Drainage of a Maxwell Liquid down a Vertical Plate. Chemical Engineering Journal, 1, 255-257.

  69. 69. Djukic, D.S. (1974) Hiemenz Magnetic Flow of Power-Law Fluids. Journal of Applied Mechanics, 41, 822-823.

  70. 70. Rajagopal, K.R. (1980) Viscometric Flows of Third Grade Fluids. Mechanics Research Communications, 7, 21-25.

  71. 71. Rajagopal, K.R. and Gupta, A.S. (1981) On a Class of Exact Solutions to the Equations of Motion of a Second Grade Fluid. International Journal of Engineering Science, 19, 1009-1014.

  72. 72. Dorier, C. and Tichy, J. (1992) Behavior of a Bingham-Like Viscous Fluid in Lubrication Flows. Journal of Non-Newtonian Fluid Mechanics, 45, 291-310.

  73. 73. Zhou, X.-F. and Gao, L. (2007) Effect of Multipolar Interaction on the Effective Thermal Conductivity of Nanofluids. Chinese Physics B, 16, 2028-2032.

  74. 74. Cui, Z.-W., Liu, J.-X., Yao, G.-J. and Wang, K.-X. (2010) Borehole Guided Waves in a Non-Newtonian (Maxwell) Fluid-Saturated Porous Medium. Chinese Physics B, 19, Article ID: 084301.

  75. 75. Bhattacharyya, K. and Layek, G.C. (2011) MHD Boundary Layer Flow of Dilatant Fluid in a Divergent Channel with Suction or Blowing. Chinese Physics Letters, 28, Article ID: 084705.

  76. 76. Casson, N. (1959) A Flow Equation for Pigment Oil Suspensions of the Printing Ink Type. In: Mill, C.C., Ed., Rheology of Disperse Systems, Pergamon Press, Oxford, 84-102.

  77. 77. Mustafa, M., Hayat, T., Pop, I. and Aziz, A. (2011) Unsteady Boundary Layer Flow of a Casson Fluid Due to an Impulsively Started Moving Flat Plate. Heat Transfer, 40, 563-576.

  78. 78. Bhattacharyya, K., Hayat, T. and Alsaedi, A. (2014) Exact Solution for Boundary Layer Low of Casson Fluid over a Permeable Stretching/Shrinking Sheet. Zeitschrit fur Angewandte Mathematik und Mechanik, 94, 522-528.

  79. 79. Bhattacharyya, K., Hayat, T. and Alsaedi, A. (2013) Analytic Solution for Magneto Hydrodynamic Boundary Layer Low of Casson Fluid over a Stretching/Shrinking Sheet with Wall Mass Transfer. Chinese Physics B, 22, Article ID: 024702.

  80. 80. Nakamura, M. and Sawada, T. (1988) Numerical Study on the Flow of a Non-Newtonian Fluid through an Axisymmetric Stenosis. Journal of Biomechanical Engineering, 110, 137-143.

  81. 81. Bird, R.B., Dai, G.C. and Yarusso, B.J. (1983) The Rheology and Flow of Viscoplastic Materials. Reviews in Chemical Engineering, 1, 1-70.

  82. 82. Akbar, N.S. (2015) Influence of Magnetic Field on Peristaltic Flow of a Casson Fluid in an Asymmetric Channel: Application in Crude Oil Refinement. Journal of Magnetism and Magnetic Materials, 378, 320-326.

  83. 83. Nadeem, S. (2014) MHD Three Dimensional Boundary Layer Flow of Casson Nanofluid Past a Linearly Stretching Sheet with Convective Boundary Condition. IEEE Transactions on Nanotechnology, 13, 109-115.

  84. 84. Hussain, T., Shehzad, S.A., Alsaedi, A., Hayat, T. and Ramzan, M. (2015) Flow of Casson Nanofluid with Viscous Dissipation and Convective Conditions: A Mathematical Model. Journal of Central South University, 22, 1132-1140.