Advances in Mathematical Physics

Volume 2014 (2014), Article ID 809367, 13 pages

http://dx.doi.org/10.1155/2014/809367

Research Article

Analysis of Heat Transfer in Berman Flow of Nanofluids with Navier Slip, Viscous Dissipation, and Convective Cooling

O. D. Makinde,¹ S. Khamis,² M. S. Tshehla,¹ and O. Franks³

¹Faculty of Military Science, Stellenbosch University, Private Bag Box X2, Saldanha 7395, South Africa

²Mathematics and Computational Science and Engineering, Nelson Mandela African Institute of Science and Technology (NM-AIST), Arusha, Tanzania

³Faculty of Engineering Built Environment and Information Technology, Nelson Mandela Metropolitan University, P.O. Box 77000, Port Elizabeth 6031, South Africa

Received 10 February 2014; Accepted 26 February 2014; Published 31 March 2014

Academic Editor: R. N. Jana

Copyright © 2014 O. D. Makinde 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.

Abstract.

Heat transfer characteristics of a Berman flow of water based nanofluids containing copper (Cu) and alumina (Al₂O₃) as nanoparticles in a porous channel with Navier slip, viscous dissipation, and convective cooling are investigated. It is assumed that the exchange of heat with the ambient surrounding takes place at the channel walls following Newton’s law of cooling. The governing partial differential equations and boundary conditions are converted into a set of nonlinear ordinary differential equations using appropriate similarity transformations. These equations are solved analytically by regular perturbation methods with series improvement technique and numerically using an efficient Runge-Kutta Fehlberg integration technique coupled with shooting scheme. The effects of the governing parameters on the dimensionless velocity, temperature, skin friction, pressure drop, and Nusselt numbers are presented graphically and discussed quantitatively.

1. Introduction

The study of fluid flow and heat transfer between two porous boundaries has gained tremendous attention of researchers due to its wide applications in engineering and industrial processes. Some of the practical interests include problems dealing with transpiration cooling where the walls of a channel containing heated fluid are protected from overheating by passing cooler fluid over the exterior surface of the channel; fluid flow occurring during the separation of isotopes of Uranium-235 and Uranium-238 by gaseous diffusion in order to produce fuel for nuclear reactors; controlling boundary layer flow over aircraft wings by injection or suction of fluid out of or into the wing; lubrication of porous bearings; petroleum technology; ground water hydrology; seepage of water in river beds; purification and filtration processes; methods of decreasing rates of heat transfer in combustion chambers exhaust nozzles and porous walled flow reactors, and so forth. In a pioneering work, Berman [1] presented an exact solution of the Navier-Stokes equations that describes the steady two-dimensional flow of an incompressible viscous fluid along a channel with parallel rigid porous walls, under the action of uniform suction or injection of fluid at the surface. Sellars [2] extended Berman’s work to high suction Reynolds number. Yuan [3] considered the flow in a channel with porous walls. He obtained solutions for small suction and injection values and asymptotic solution valid at large injection values. Terrill [4, 5] gave an exact series solution for the fully developed laminar flow in a pipe of circular cross section with porous wall driven by a spatially variable and time independent suction or injection at the surface. Studies on developing flow in porous-walled ducts with suction and injection effects were carried out by Sorour et al. [6] and Zaturska et al. [7]. In view of the above interests, several researchers have also investigated the heat transfer problems between two permeable parallel walls under different physical situations [8–10].

Moreover, the applications of conventional heat transfer fluids such as water and glycol mixture in engineering flow processes are limited due to their low thermal properties. A potential solution to improve these thermal properties is to add nanoparticles into the conventional fluids, hence forming the so-called nanofluids as coined by Choi [11]. Nanofluids contain thermally conducting submicron solid particles and have great potential as a high-energy carrier. The nanoparticles such as copper, alumina, titania, and copper oxide, unlike larger-sized particles, can be suspended stably within the conventional fluids without settling out of suspension. Nanofluids are free from numerous problems such as abrasion, clogging, and high pressure loss and are considered to be the next-generation working fluids in modern heat transfer technologies [12]. Experimental results [13–15] have shown that, even with small solid volume fraction of nanoparticles (usually less than 5%), the thermal conductivity of heat transfer fluids can be enhanced by 10–50%. Several authors [16–18] have also theoretically investigated the heat transfer enhancement of nanofluids under different physical conditions.

Meanwhile, advances in the manufacture of microdevices have enabled experimental investigation of fluid flow in nano- and microscale, and many experimental results have provided evidences to support the slip condition [19, 20]. In order to describe the slip characteristics of fluid on the solid surface, Navier [21] introduced a more general boundary condition, namely, the fluid velocity component tangential to the solid surface, relative to the solid surface, which is proportional to the shear stress on the fluid-solid interface. The proportionality is called the slip length which describes the slipperiness of the surface. However, from the literature survey, it is found that no study has been conducted on the heat transfer characteristics of Berman flow of nanofluids with Navier slip, viscous dissipation, or convective cooling at the walls. Hence, the present study is an attempt in this direction. The flow of water base nanofluids containing copper (Cu) and alumina (Al₂O₃) as nanoparticles in a uniformly porous wall channel with Navier slip, viscous dissipation, and convective heat exchange with the ambient surrounding is investigated. Sections 2–4 give in more details the model nonlinear governing equations together with the analytical and numerical solution techniques employed to tackle the problem. In Section 5, we present graphically and discuss the main features of the flow and heat transfer characteristics in a range of governing parameters. Final conclusions are drawn in Section 6. To the best of our knowledge, the results of this paper are new and they have not been published before.

2. Problem Formulation

Consider a two-dimensional steady flow of a viscous incompressible water base nanofluids containing copper (Cu) and alumina (Al₂O₃) as nanoparticles in a uniformly porous wall channel. The channel wall is subjected to Navier slip and convectively exchange heat with the ambient surrounding. We choose the Cartesian coordinates system in such a way that the -axis is taken along the channel and the -axis is normal to it as shown in Figure 1.

Figure 1: Schematic diagram of the physical system.

The governing equations which are those of conservation of mass, momentum, and energy are [1–5, 8–10]

where are the velocity components of the nanofluid in the directions, respectively, uniform wall suction velocity, is the pressure, is the nanofluid temperature, is the effective dynamic viscosity of the nanofluid [22], is the effective thermal conductivity of the nanofluid [23], and and are the nanofluid density and the heat capacitance of the nanofluid, respectively, which are given by [13–18].

where is the nanoparticles solid volume fraction, is the reference density of the fluid fraction, is the reference density of the solid fraction, is the viscosity of the fluid fraction, is the thermal conductivity of the fluid fraction, is the specific heat at constant pressure, and is the thermal conductivity of the solid volume fraction.

Due to the symmetric nature of the flow, the boundary conditions at the channel centerline and at the porous wall may be written as [1–4]

where is the ambient surrounding temperature, is the coefficient of heat transfer, and is the Navier slip coefficient. Introduce the stream function and vorticity into the governing equations (1)–(4) as follows:

After eliminating the pressure from (2) and (3), we obtain

The following dimensionless variables and parameters are introduced into (8), together with their corresponding boundary conditions:

and we obtain

with

where is the flow Reynolds number such that represents wall suction and represents wall injection; Ec is the Eckert number; Bi is the Biot number; Pr is the base fluid Prandtl number; , , , and can be easily determined from the thermophysical properties of the base fluid and the nanoparticles; and is the Navier slip parameter such that corresponds to no slip, while full lubrication is described in the limit . We seek a similarity form of solution due to Berman [1]; that is,

Equations (10)-(11) together with the boundary conditions in (12) then become

where . The dimensionless fluid axial pressure gradient is given as follows:

where

Other quantities of practical interest in this study are the local skin friction coefficient and the local Nusselt number Nu, which are defined as

where is the skin friction and is the heat flux at the channel walls which are given by

Using (9) and (13), we substitute (22) into (21) and obtain

In the following section, the boundary value problem (14)–(18) was solved analytically using regular perturbation method and numerically by the Runge-Kutta-Fehlberg method with shooting technique [24]. The results are utilised to compute the fluid pressure gradient, local skin friction, and local Nusselt number as highlighted in (20) and (23).

3. Perturbation Method

Due to the nonlinear nature of the model equations (14)–(18), it is convenient to form a power series expansion in the parameter ; that is,

Substituting the solution series in (24) into (14)–(18) and collecting the coefficients of like powers of , we obtain the following.

Zeroth Order. Consider

with

Higher Order . Consider

with

The equations are solved iteratively and the series solutions for the velocity and temperature fields are given as

where . We remark here that if we set the parameters and in (29), we will recover the solution for the classical case of conventional fluid given in [1–4]. Using a computer symbolic algebra package (MAPLE), several terms of the above solution series in (29)–(31) are obtained. From (29)–(31) together with (20) and (23), we obtained the series solutions for the skin friction, Nusselt number, and axial pressure gradient as follows:

where the expression for , is given in the appendix. We are aware that these power series solutions are valid for very small parameter values of . However, using Hermite-Padé approximation technique (see Makinde [25]) that is based on the series summation and improvement method, the usability of the extended solution series is improved beyond small parameter values of .

4. Numerical Procedure

An efficient finite difference approach based on Runge-Kutta-Fehlberg method with shooting technique [24] has been employed to numerically solve the coupled nonlinear ordinary differential equations (14)–(16) subject to the boundary conditions (17)-(18) for different values of governing parameters. The boundary value problem is first transformed into an initial value problem (IVP). Let

Substituting (33) into (14)–(18), we obtain a system of first-order differential equations, respectively, as follows:

subject to the initial conditions:

By applying the shooting method with the Newton-Raphson algorithm to guess the unspecified conditions , , , and in (35), the resulting initial value problem is then integrated numerically until the boundary conditions at are achieved. The step size is taken as and the convergence criteria were set to 10⁻⁷.

5. Results and Discussion

The flow of water base nanofluids containing Cu and Al₂O₃ as nanoparticles and their heat transfer characteristics in a symmetrically porous channel with Navier slip and convective cooling at the surface are investigated. The governing partial differential equations and the corresponding boundary conditions are converted into a set of nonlinear ordinary differential equations and tackled both analytically using the perturbation method coupled with series improvement technique and numerically using Runge-Kutta Fehlberg integration technique coupled with shooting scheme. Thermophysical properties of base fluid and nanoparticles are presented in Table 1. For pure water, the momentum diffusivity is dominant and convection is very effective in transferring within the fluid in comparison to pure conduction. Following [16–18], we take in the numerical computation. Note that when , no nanoparticle is present in the based fluid (water). The solid volume fraction in the base fluid is taken as to 0.3 (i.e., ranging from 0 to 30 percent). In order to get a clear insight into the entire flow structure and thermal development, we have assigned numerical values to other parameters encountered in the problem. Numerical solutions are displayed in Tables 2 and 3 together with Figures 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, and 23. The numerical values of suction driven normal velocity profile are displayed in Table 2. In the absence of nanoparticles and Navier slip at the channel walls, the results agreed well with the one already in the literature (see [1–3]) for the flow of conventional fluid in a symmetrical porous wall channel. However, Cu-water nanofluid is affected more by the combined effects of suction and Navier slip at wall in comparison to Al₂O₃-water nanofluid. Table 3 shows the perfect agreement between the series solution and the numerical solution for the axial velocity profiles along the channel centerline with increasing concentration of nanoparticles for both Cu-water and Al₂O₃-water nanofluids.

Table 1: Thermophysical properties of the fluid phase (water) and nanoparticles [13–18].

Physical properties
Fluid phase (Water)
Cu
Al₂O₃

(J/kgK)
4179
385
765

(kg/m³)
997.1
8933
3970

(W/mK)
0.613
400
40

Table 2: Computation showing the normal velocity profiles for .

,  

Reference [1–3]

,  

Present

,  

Cu-water

,  

Al₂O₃-water

0.0
0.00000
0.000000
0.000000
0.000000

0.1
0.14874
0.148739
0.160680
0.160920

0.2
0.29455
0.294548
0.317732
0.318184

0.3
0.43449
0.434493
0.467527
0.468132

0.4
0.56564
0.565639
0.606428
0.607103

0.5
0.68504
0.685042
0.730781
0.731419

0.6
0.78974
0.789739
0.836887
0.837373

0.7
0.87672
0.876724
0.920958
0.921192

0.8
0.94292
0.942920
0.979045
0.978994

0.9
0.98514
0.985137
1.006937
1.006713

1.0
1.00000
1.000000
1.000000
1.000000

Table 3: Computation showing the axial velocity profiles for , , , and .

Cu-water

(series)

Cu-water

(numerical)

Al₂O₃-water

(series)

Al₂O₃-water

(numerical)

0.00
1.581405
1.581405
1.581405
1.581405

0.01
1.583733
1.583733
1.584084
1.584084

0.05
1.594688
1.594688
1.596218
1.596218

0.10
1.612850
1.612850
1.615303
1.615303

0.15
1.637660
1.637660
1.640343
1.640343

0.20
1.672037
1.672037
1.674097
1.674097

0.25
1.721196
1.721196
1.721395
1.721395

0.30
1.795420
1.795420
1.791485
1.791485

Figure 2: Nanofluids axial velocity profiles.

Figure 3: Axial velocity profiles with increasing .

Figure 4: Axial velocity profiles with increasing .

Figure 5: Axial velocity profiles with increasing .

Figure 6: Nanofluids normal velocity profiles.

Figure 7: Normal velocity profiles with increasing .

Figure 8: Normal velocity profiles with increasing .

Figure 9: Normal velocity profiles with increasing .

Figure 10: Nanofluids temperature profiles.

Figure 11: Temperature profiles with increasing .

Figure 12: Temperature profiles with increasing .

Figure 13: Temperature profiles with increasing .

Figure 14: Temperature profiles with increasing .

Figure 15: Temperature profiles with increasing .

Figure 16: Temperature profiles with increasing .

Figure 17: Skin friction with increasing .

Figure 18: Skin friction with increasing and .

Figure 19: Axial pressure gradient with increasing .

Figure 20: Axial pressure gradient with increasing and .

Figure 21: Nusselt number with increasing .

Figure 22: Nusselt number with increasing and .

Figure 23: Nusselt number with increasing , , and .

5.1. Velocity Profiles with Parameter Variation

The effects of parameter variation on both the axial and normal velocity components are displayed in Figures 2–9. Generally, it is interesting to note that the effect of Navier slip is to cause flow reversal at the channel walls. In Figure 2, it is observed that Cu-water moves faster along the channel centerline region and subjected to higher flow reversal at the wall as compared to Al₂O₃-water in the presence of suction. With Cu-water as the working nanofluid, it is observed that increasing nanoparticles volume fraction concentration from 0 to 30% increases both the axial velocity along the centerline region and the flow reversal at the channel walls as shown in Figure 3. Similar trend is observed in Figures 4-5 with a growing suction parameter and Navier slip parameter. Figures 6–9 show the nanofluids normal velocity profiles. With suction and Navier slip, the Cu-water moves faster to the wall as compared to Al₂O₃-water as illustrated in Figure 6. A further increase in normal velocity towards the walls is observed with a growing in suction parameter, Navier slip parameter, and nanoparticles volume fraction as shown in Figures 7–9.

5.2. Temperature Profiles with Parameter Variation

Figures 10–16 illustrate the nanofluids temperature profiles across the channel with different parameter variation. Generally, a decrease in the fluid temperature near the channel walls is observed due to convective heat loss to ambient surrounding. It is noteworthy that the temperature of Cu-water nanofluid is generally higher than that of Al₂O₃-water nanofluid under the same flow condition as shown in Figure 10. In Figure 11, it is observed that the nanofluid temperature decreases with growing in nanoparticles volume fraction. Similar effect of a decrease in nanofluid temperature is observed in Figure 12 with Cu-water as working nanofluid as the Biot number increases. This is expected, since an increase in Biot number indicates a rise in convective cooling due to heat loss to the ambient surrounding from the channel walls. Meanwhile, a combine increase in the suction, Navier slip, and viscous dissipation as shown in Figures 13–15 causes an increase in the nanofluid temperature. This may be attributed to the fact that , , and Ec increase the internal heat generation within the fluid due velocity gradient increases, leading to a rise in temperature. Figure 16 elucidates the temperature profiles with increasing axial distance along the channel. The nanofluid temperature decreases within the channel centerline region and increases near the wall region with increasing axial distance. Moreover, it is interesting to note that at point within the channel, the nanofluid temperature is independent of the axial distance.

5.3. Skin Friction, Pressure Gradient, and Nusselt Number

Figures 17-18 depict the skin friction profiles for both Cu-water and Al₂O₃-water nanofluids at the channel walls. The skin friction generally increases with an increase in nanoparticles volume fraction; however, it is noticed that the skin friction produced by Cu-water is more intense than the one produced by Al₂O₃-water as shown in Figure 17. This is expected since the velocity gradient of Cu-water near the channel walls is higher than that of Al₂O₃-water. In Figure 18, it is observed that the skin friction generally increases with an increase in Navier slip. Meanwhile, a growing in suction increases the skin friction, while a growing in injection decreases the skin friction. The pressure drop along the channel is illustrated in Figures 19 and 20. For both Cu-water and Al₂O₃-water nanofluids, the pressure drop increases with increasing nanoparticles volume fraction. Interestingly, the pressure drop produced by Al₂O₃-water is slightly higher than that of Cu-water as shown in Figure 19. Figure 20 shows a general increase in pressure drop with a rise in Navier slip. A growing in suction decreases the pressure drop, while a growing in injection increases the pressure drop along the channel. Figures 21–23 elucidate the heat transfer rate at the channel walls with different parameter variation. It is observed that the wall heat transfer rate (Nu) decreases with an increase in the nanoparticles volume fraction as shown in Figure 21. A slight increase in Al₂O₃-water Nusselt number is noticed as compared to Cu-water Nusselt number. In Figure 22, it is seen that the increase in the Navier slip parameter results in increase in Nusselt number. This may be due to a rise in the nanofluid temperature gradient at the channel walls. Meanwhile, as the suction increases, the heat flux at the wall increases, while a decrease in wall heat is observed with a rise in injection. The strength of the wall heat flux is enhanced with increasing axial distance, viscous dissipation, and convective cooling as illustrated in Figure 23. This can be attributed to a rise in the temperature gradient due to convective heat exchange with the ambient along the channel walls.

6. Conclusion

The combined effects of viscous dissipation, Navier slip, and convective cooling on Berman flow and heat transfer of water base nanofluids containing Cu and Al₂O₃ as nanoparticles are investigated. The nonlinear model problem is tackled both analytically using perturbation series method and numerically using Runge-Kutta Fehlberg integration technique coupled with shooting scheme. We summarize below some of the essential features of physical interest from the above analysis.

(i)Cu-water nanofluid moves faster with enhanced flow reversal at the walls as compared to Al₂O₃-water nanofluid.(ii)Nanofluids velocity and flow reversal at the walls increase with suction, , and .(iii)Cu-water produced higher temperature as compared to Al₂O₃-water. The nanofluids temperature increases with suction, , and Ec, but it decreases with Bi and .(iv)The skin friction produced by Cu-water is more intense than that of Al₂O₃-water. The skin friction increases with suction , , and but decreases with injection .(v)The pressure drop produced by Al₂O₃-water in more than that of Cu-water. The pressure drop is enhanced by injection, , and but decreases by suction.(vi)The Nusselt number increases with suction, , Bi, Ec, and , but it decreases with injection and . A slight increase in Nu for Al₂O₃-water is noticed as compared to Cu-water.

Appendix

Nomenclature

: | Velocity components
: | Coordinates
: | Nanofluid thermal conductivity
: | Prandtl number
: | local Biot number
: | Ambient temperature
: | Dimensionless stream function
: | Wall suction velocity
: | Temperature
: | Dimensionless temperature
: | Reynolds number
: | Specific heat at constant pressure
: | Solid fraction thermal conductivity
: | Base fluid thermal conductivity
: | Dimensionless axial velocity
: | Eckert number
: | Dimensionless axial coordinate
: | Axial pressure gradient coefficient
: | Channel half width
: | Nusselt number
: | Skin friction coefficient.

Greek Symbols

: | Stream function
: | Dimensionless temperature
: | Nanofluid dynamic viscosity
: | Nanofluid thermal diffusivity
: | Dimensionless normal coordinate
: | Nanofluid density
: | Solid fraction density
: | Base fluid kinematic viscosity
: | Base fluid dynamic viscosity
: | Solid volume fraction parameter
: | Base fluid density
: | Slip coefficient
: | Navier slip parameter
: | Vorticity
: | Dimensionless stream function
: | Dimensionless vorticity
: | Dimensionless temperature.

Conflict of Interests

The authors declare that there is no conflict of interests.

References

A. S. Berman, “Laminar flow in channels with porous walls,” Journal of Applied Physics, vol. 24, pp. 1232–1235, 1953.

J. R. Sellars, “Laminar flow in channels with porous walls at high suction Reynolds numbers,” Journal of Applied Physics, vol. 26, no. 4, pp. 489–490, 1955.

S. W. Yaun, “Further investigation of laminar flow in channels with porous walls,” Journal of Applied Physics, vol. 27, no. 3, pp. 267–269, 1956.

R. M. Terrill, “An exact solution for flow in a porous pipe,” Journal of Applied Mathematics and Physics, vol. 33, no. 4, pp. 547–552, 1982.

R. M. Terrill, “Laminar flow in a porous tube,” Journal of Fluids Engineering, vol. 105, pp. 303–306, 1983.

M. M. Sorour, M. A. Hassab, and S. Estafanous, “Developing laminar flow in a semiporous two-dimensional channel with nonuniform transpiration,” International Journal of Heat and Fluid Flow, vol. 8, pp. 44–54, 1987.

M. B. Zaturska, P. G. Drazin, and W. H. H. Banks, “On the flow of a viscous fluid driven along a channel by suction at porous walls,” Fluid Dynamics Research, vol. 4, no. 3, pp. 151–178, 1988.

R. J. Pederson and R. B. Kinney, “Entrance-region heat transfer for laminar flow in porous tubes,” International Journal of Heat and Mass Transfer, vol. 14, pp. 159–161, 1971.

O. D. Makinde, “Extending the utility of perturbation series in problems of laminar flow in a porous pipe and a diverging channel,” Australian Mathematical Society B. Applied Mathematics, vol. 41, no. 1, pp. 118–128, 1999.

G. D. Raithby, “Laminar heat transfer in the thermal entrance region of circular tubes and two-dimensional rectangular ducts with wall suction and injection,” International Journal of Heat and Mass Transfer, vol. 14, pp. 224–243, 1971.

S. U. S. Choi, “Enhancing thermal conductivity of fluids with nanoparticles,” in Proceedings of the ASME International Mechanical Engineering Congress and Exposition, pp. 99–105, ASME, San Francisco, Calif, USA, 1995.

B. C. Pak and Y. I. Cho, “Hydrodynamic and heat transfer study of dispersed fluids with submicron metallic oxide particles,” Experimental Heat Transfer, vol. 11, no. 2, pp. 151–170, 1998.

D. Wen and Y. Ding, “Experimental investigation into the pool boiling heat transfer of aqueous based c-alumina nanofluids,” Journal of Nanoparticle Research, vol. 7, no. 2, pp. 265–274, 2005.

J. Wang, J. Zhu, X. Zhang, and Y. Chen, “Heat transfer and pressure drop of nanofluids containing carbon nanotubes in laminar flows.,” Experimental Thermal and Fluid Science, vol. 44, pp. 716–721, 2013.

Y. Xuan and Q. Li, “Investigation on convective heat transfer and flow features of nanofluids,” Journal of Heat Transfer, vol. 125, no. 1, pp. 151–155, 2003.

H. F. Oztop and E. Abu-Nada, “Numerical study of natural convection in partially heated rectangular enclosures filled with nanofluids,” International Journal of Heat and Fluid Flow, vol. 29, pp. 1326–1336, 2008.

O. D. Makinde, “Computational modelling of nanofluids flow over a convectively heated unsteady stretching sheet,” Current Nanoscience, vol. 9, pp. 673–678, 2013.

S. Kakac and A. Pramuanjaroenkij, “Review of convective heat transfer enhancement with nanofluids,” International Journal of Heat and Mass Transfer, vol. 52, no. 13/14, pp. 3187–3196, 2009.

P. Huang and K. S. Breuer, “Direct measurement of slip length in electrolyte solutions,” Physics of Fluids, vol. 19, Article ID 028104, 2007.

M. J. Martin and D. I. Boyd, “Momentum and heat transfer in a laminar boundary layer with slip flow,” Journal of Thermophysics and Heat Transfer, vol. 20, no. 4, pp. 710–719, 2006.

C. L. M. H. Navier, “Mémoire sur les lois du mouvement des fluides,” Mémoires de l'Académie Royale des Sciences de l'Institut de France, vol. 6, pp. 389–440, 1823.

H. C. Brinkman, “The viscosity of concentrated suspensions and solution,” The Journal of Chemical Physics, vol. 20, pp. 571–581, 1952.

J. C. Maxwell, Electricity and Magnetism, Clarendon, Oxford, UK, 3rd edition, 1904.

T. Y. Na, Computational Methods in Engineering Boundary Value Problems, Academic Press, New York, NY, USA, 1979.

O. D. Makinde, “Hermite-Padé approach to thermal radiation effect on inherent irreversibility in a variable viscosity channel flow,” Computers & Mathematics with Applications, vol. 58, no. 11-12, pp. 2330–2338, 2009.