Solution of the boundary layer flow of various nanofluids over a moving semi-infinite plate using HPM/Ivairiu nanoskysciu ribinio tekejimo virs judancios pusiau begalines ploksteles tyrimas homotopiniu perturbaciju metodu.

Dalir, N. ; Nourazar, S.S.

1. Introduction

The performance of thermal systems is subjected to a primary limitation due to the low thermal conductivity of conventional heat transfer fluids. In order to improve the thermal conductivity, nanoscale particles being dispersed in a base fluid, known as nanofluid, are used. Thus nanofluids are suspensions of nanoparticles in fluids. Nanofluids offer considerable advantages over conventional heat transfer fluids. Compared to pure fluids, they enhance thermal and transport properties considerably. Wherever heat transfer enhancement is crucial such as in nuclear reactors, transportation and electronics, nanofluids are usable. They can enhance thermal conductivity of the base fluid enormously. Nanofluids are also very stable and have no additional problems such as non-Newtonian behavior, the reason being the tiny size of nanoparticles.

The research on nanofluids began over a decade ago, by focusing on measuring and modeling their effective thermal conductivity and viscosity. Choi et al. [1] added a small amount of nanoparticles to conventional heat transfer fluids and observed the increase of thermal conductivity. Das [2] presented a numerical investigation on the convective heat transfer performance of nanofluids over a permeable stretching surface in the presence of partial slip, thermal buoyancy and temperature dependent internal heat generation or absorption. Makinde and Aziz [3] studied numerically the boundary layer flow induced in a nanofluid due to a linearly stretching sheet with a convective boundary condition at the sheet surface. Kandasamy et al. [4] solved numerically the problem of laminar fluid flow which results from the stretching of a vertical surface with variable stream conditions in a nanofluid. They used a model for the nanofluid which incorporates the effects of Brownian motion and thermo-phoresis in the presence of magnetic field. Anwar et al. [5] investigated theoretically the problem of free convection boundary layer flow of nanofluids over a non-linear stretching sheet, incorporating the effects of buoyancy parameter, the solutal buoyancy parameter and the power law velocity parameter. Boundary-layer flow problem over a moving or fixed flat plate is a classical problem, which has been investigated by many researchers, for example Bachok et al. [6] studied the steady-state boundary-layer flow of a nanofluid over a moving semi-infinite flat plate in a uniform free stream, and found that dual solutions exist when the plate and the free stream move in the opposite directions. Bachok et al. [7] investigated the problem of a uniform free stream of nanofluid parallel to a fixed or moving flat plate. They solved the problem using the shooting method. Ahmad et al. [8] solved the Blasius and Sakiadis problems in nanofluids and concluded that the inclusion of nanoparticles into the base fluid had resulted in an increase of the skin friction and heat transfer coefficients.

In the present paper, the two-dimensional steadystate boundary layer flow of nanofluids over an impermeable semi-infinite moving flat horizontal plate embedded in the water-based nanofluid is studied. It is assumed that the flat plate moves with a constant velocity. The governing equations, i.e. mass and momentum conservation equations, are transformed using the similarity transformations to a nonlinear ordinary differential equation (ODE), and then the resulting ODE is solved using the homotopy perturbation method (HPM). Six types of nanoparticles, i.e., copper (Cu), alumina ([Al.sub.2][O.sub.3]), titania (Ti[O.sub.2]), copper oxide (CuO), silver (Ag), and silicon (Si[O.sub.2]) in the water based fluid with Pr = 6.2 are considered. The velocity and stream function profiles are plotted for various nanoparticles and for various values of the nanoparticle volume fraction. The effect of the nanoparticle volume fraction on the flow characteristics, and mainly on the local skin friction coefficient, is investigated.

2. Mathematical formulation

The steady-state two-dimensional laminar boundary layer flow over a continuously moving flat horizontal plate embedded in a water-based nanofluid is considered. The nanofluid can contain each of six types of nanoparticles including Cu, [Al.sub.2][O.sub.3], Ti[O.sub.2], CuO, Ag, and Si[O.sub.2]. It is assumed that the plate has a constant velocity. A uniform spherical size and shape is assumed for the nanoparticles. It is also assumed that the base fluid and the nanoparticles are in the thermal equilibrium, and no velocity slip occurs between the base fluid and the nanoparticles [9]. Considering these assumptions, the laminar boundary layer equations of mass and momentum conservation are as follows:

[partial derivative]u/[partial derivative]x + [partial derivative]v/[partial derivative]y = 0;(1

u [partial derivative]u/[partial derivative]x + v [partial derivative]u/[partial derivative]y = [[[mu].sub.nf] / [[rho].sub.nf]] [[[partial derivative].sup.2]u/ [[partial derivative]y.sup.2]]. (2)

The boundary conditions for the fluid velocity are as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (3)

in which [U.sub.w] is the plate velocity which is constant, and u and v are the velocity components in x- and y-directions, respectively. [[rho].sub.nf] is the density of the nanofluid, and [[mu].sub.nf] is the viscosity of the nanofluid, which are given by the following relations [10]:

[[rho].sub.nf] = (1-[phi]) [[rho].sub.f] + [[phi] [[rho].sub.s], [[mu].sub.nf] = [[mu].sub.f] /[(1-[phi]).sup.2.5], (4)

where [phi] is the nanoparticle volume fraction, and [[rho].sub.f] and [[rho].sub.s] are the densities of fluid and solid fractions, respectively.

The dimensionless similarity variable and the dimensionless stream-function used to transform the governing equations to an ordinary differential equation are defined as:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (5)

where [Re.sub.x] = [U.sub.w]x/[v.sub.f] is the local Reynolds number, in which [v.sub.f] is the kinematic viscosity of the base fluid (water). [psi] (x,y) is the stream function which identically satisfies Eq. (1) and is defined as u = d[psi] / dy, v = - d[psi]/ dx. By the use of the similarity parameters (5), the boundary layer momentum Eq. (2) and the boundary conditions (3) transform to the following forms:

f'" + 1/2 [(1 - [phi]).sup.2.5](1 - [phi] + [phi] [[rho].sub.s]/ [[rho].sub.f]) ff" = 0; (6)

F (0)= 0, f' (0) = 1, f' ([infinity])=0 . (7)

In Eqs. (6) and (7), prime denotes differentiation with respect to [eta]. The significant quantity is the local skin friction coefficient [C.sub.f,x] defined as [C.sub.f,x] = [[tau].sub.w]/ [[rho].sub.f][U.sub.w.sup.2], in which the plate surface shear stress is given as [[tau].sub.w] = [[mu].sub.nf] [(du/dy).sub.y=0]. Use of the similarity parameters (5) gives [11]:

[C.sub.f,x] [Re.sup.0.5.sub.x] = f" (0)/[(1-[phi].sup.2.5] (8)

3. Solution by homotopy perturbation method (HPM)

Using HPM [12], the original nonlinear ODE (which cannot be solved easily) is divided into some linear ODEs (which are solved easily in a recursive manner by mathematical symbolic software such as Mathematica or Maple).

At first, the governing ODE (6) and the boundary conditions (7) are written as:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]; (9)

u (0) = 0, u' (0) = 1, u' ([infinity]) = 0. (10)

Then, a homotopy is constructed in the following form:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (11)

According to HPM, the following serious in terms of powers of p is substituted in Eq. (11):

u = [u.sub.0] + [pu.sup.1] + [p.sup.2] [u.sub.2] + ... (12)

After some algebraic manipulation, equating the identical powers of p to zero gives:

[p.sup.0:] [u'".sub.0] - [[alpha].sup.2] [u'.sub.0] = 0; [u.sub.0] (0) = 0; [u'.sub.0] (0) = 1; [u'.sub.0] ([infinity]) = 0; (13)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (14)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (15)

Eq. (13) for [p.sup.0] has the following solution:

[u.sub.0] ([eta]) = 1/[alpha] (1-exp(- [alpha][eta])). (16)

Here [alpha] is a constant which is further to be determined. If solution (16) for [u.sub.0] is substituted in the equation for [p.sup.1], Eq. (14), it will become as:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (17)

Eq. (17) for [u.sub.1] can be solved in an unbounded domain under the boundary conditions [u.sub.1] (0) = 0, [u'.sub.1] (0) = 0, [u'.sub.1] ([infinity]) = 0 (as it is shown in the Appendix) [13], which gives [u.sub.1] as:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (18)

in which [alpha] = [([OMEGA] / 2).sup.0.5] and [OMEGA] = [(1 - [phi]).sup.2.5] [1 - [phi] + [phi] ([[rho].sub.s] /[[rho].sub.f])]. It should be noted that [alpha] can be [alpha] = [+ or -] [([OMEGA] / 2).sup.0.5], but here as [alpha] is demanded to be positive ([alpha] > 0), therefore [alpha] = [([OMEGA] / 2).sup.0.5]. Thus the first-order approximate semianalytical solution f([eta]) = u([eta]) = [u.sub.0]( [eta])+[u.sub.1]( [eta]) becomes as:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (19)

According to Eq. (19), the dimensionless plate surface shear stress is as:

f" (0) = [3/2] [alpha] + 5 [OMEGA]/12 [alpha] (20)

4. Results and discussion

A small computer code in the symbolic software Mathematica is written and HPM solutions to the governing ordinary differential Eq. (6) with the boundary conditions (7) are obtained. The solutions are shown in tables and diagrams. In Table 1 the density of water and nanoparticles used in the present study are given.

Table 2 gives the HPM solution values of the dimensionless fluid velocity gradient at the plate surface [alpha] = f" (0) for Cu-water, [Al.sub.2][O.sub.3]-water, Ti[O.sub.2]-water, CuO-water, Ag-water and SiO2-water working fluids for various values of the nanoparticle volume fraction [phi]. It can be seen that the values of f" (0) are equal for Cu-water, [Al.sub.2][O.sub.3]-water, Ti[O.sub.2]-water, CuO-water, Ag-water and Si[O.sub.2]-water working fluids in zero nanoparticle volume fraction (i.e., [phi] = 0). This is logical because the governing Eq. (6) reduces to f'" + 0.5f" = 0 in [phi] = 0, which is the governing equation of boundary layer flow of a pure fluid. Therefore the semi-analytical results for f" (0) in [phi] = 0 are not changed by the type of nanoparticle used. It can also be seen that when the solid nanoparticle volume fraction [phi] increases, the magnitude of f" (0) increases slightly in Cu-water, CuO-water and Ag-water working fluids, but it decreases slightly in [Al.sub.2][O.sub.3]-water, Ti[O.sub.2]-water and Si[O.sub.2]-water working fluids.

Table 3 compares HPM solution and numerical solution [7] values of the local skin friction coefficient ([C.sub.fx] [Re.sub.x.sup.0.5] = f" (0) / [(1 - [phi]).sup.2.5]) for Cu-water, [Al.sub.2][O.sub.3]-water, Ti[O.sub.2]-water, CuO-water, Ag-water and Si[O.sub.2]-water working fluids for various values of the nanoparticle volume fraction [phi] (0 [less than or equal to] [phi] [less than or equal to] 0.2). It can be seen that the HPM solutions agree within 2% error with the numerical solutions obtained using a shooting method. HPM results for Cu-water working fluid are also compared with the experimental data [8] of the local skin friction coefficient in table 3, where a good agreement within 1% error is observed. It is also seen that when [phi] increases, the local skin friction coefficient magnitude increases. It is also observed that when [phi] = 0 the local skin friction coefficient ([C.sub.f,x] [Re.sub.x.sup.0.5]) values are equal for all the working fluids. The reason is that when [phi] = 0 the nanofluid boundary layer flow problem reduces to the regular fluid boundary layer problem, and thus the nanoparticle type does not alter the values of skin friction coefficient.

Fig. 1 presents the variations of f" (0) with [phi] for various nanoparticles (i.e., Cu, [Al.sub.2][O.sub.3], Ti[O.sub.2], CuO, Ag and Si[O.sub.2]) using HPM solution from Table 2. It is seen that with the increase of [phi] the magnitude of f "(0) increases for Ag-water, Cu-water and CuO-water working fluids, but the magnitude of f" (0) decreases for Ti[O.sub.2]-water, [Al.sub.2][O.sub.3]-water and Si[O.sub.2]-water working fluids with the increase of q. Comparison of Fig. 1 with the nanoparticles densities in table 1 makes it clear that the nanoparticles with higher density result in higher magnitudes of f" (0) and the nanoparticles with lower density result in lower f" (0) magnitudes.

[FIGURE 1 OMITTED]

Fig. 2 demonstrates the variations of dimensionless skin friction group ([C.sub.f,x] [Re.sub.x.sup.0.5]) with the nanoparticle volume fraction q for various nanoparticles using HPM solution from table 3. It can be seen that for all types of nanoparticles (Ag, Cu, CuO, Ti[O.sub.2], [Al.sub.2] [O.sub.3] and Si[O.sub.2]) the dimensionless skin friction group magnitude at the plate surface [C.sub.f,x] [Re.sub.x.sup.0.5] increases when q increases. Thus it can be said that the addition of any type of nanoparticle to a regular fluid enhances the skin friction. It can also be observed that a higher nanoparticle volume fraction results in a higher dimensionless skin friction group. Thus the addition of more and more amounts of nanoparticles of any type to a fluid (up to [phi] [less than or equal to] 0.2) causes the skin friction boost. Nevertheless, as it is clear from Fig. 2, the amount of increase in [C.sub.f,x] [Re.sub.x.sup.0.5] by the addition of nanoparticles to the regular fluid is not the same for all types of nanoparticles. For instance, for Ag nanoparticles the increase of [C.sub.f,x] [Re.sub.x.sup.0.5] is higher compared to all the other nanoparticle types, and for the Si[O.sub.2] nanoparticles it is lower compared to the other nanoparticle types. Here, similar to Fig. 2, the trend of [C.sub.f,x] [Re.sub.x.sup.0.5] increase is proportional to the density of nanoparticles.

The variations of local skin friction coefficient ([C.sub.f,x]) with the local Reynolds number ([Re.sub.x]) for various values of [phi] for Cu-water working fluid is plotted in Fig. 3 using HPM solution from Table 3. The horizontal axis gives the Re number values in the laminar boundary layer flow range ([Re.sub.x] [less than or equal to] [10.sup.5]). It is seen that the [C.sub.f,x] magnitude decreases with the increase of the Re number, and lower [phi] values result in lower [C.sub.f,x]'s. Thus when the situation favors the use of nanofluid along with lower skin friction coefficient, lower nanoparticle volume fractions with higher Reynolds numbers are ideal.

Fig. 4 is the curve for the local skin friction coefficient ([C.sub.f,x]) as a function of the Reynolds number for various nanoparticles. As it can be seen, the higher the Reynolds number, the lower the skin friction coefficient values. It can also be seen that the Ag nanoparticles give the highest Cf/s and the SiO2 nanoparticles give the lowest values of the [C.sub.f,x].

[FIGURE 2 OMITTED]

[FIGURE 3 OMITTED]

Fig. 5 depicts the velocity profiles f'([eta]) for some values of [phi] ([phi] = 0, 0.1, 0.2) for Cu-water working fluid using HPM solution. It can be observed that the velocity profiles are steeper for the nanofluid cases (i.e., [phi] = 0.1 and 0.2). Thus the velocity boundary layer is considerably thinner for the nanofluid cases (i.e., [phi] = 0.1 and 0.2) compared to regular fluid case ([phi] = 0).

[FIGURE 4 OMITTED]

[FIGURE 5 OMITTED]

In Fig. 6 the velocity profiles f' ([eta]) for various nanoparticles (Cu, [Al.sub.2][O.sub.3], Ti[O.sub.2], CuO, Ag and Si[O.sub.2]) in [phi] = 0.1 are demonstrated using HPM solution of the present paper. It can be seen that the boundary layer velocity profiles are affected by the types of nanoparticles used in nanofluids. It is also observable that the steepest velocity profile is for Ag nanoparticles, and therefore the velocity boundary layer has the lowest thickness for Ag nanoparticles. With regard to the steepness of the velocity profile, the nanoparticle types Cu, CuO, Ti[O.sub.2], [Al.sub.2][O.sub.3] and Si[O.sub.2] stand on the next steps. Thus Si[O.sub.2] nanoparticles generate the thickest velocity boundary layer. It is worth mentioning that the nanofluid boundary layer thickness is inversely proportional to the density of nanoparticles used in the working fluid.

Fig. 7 shows the stream-function profiles f ([eta]) for some values of the nanoparticle volume fraction [phi] for Cu-water working fluid using HPM. Here again, thicker boundary layer can be observed for the regular fluid ([phi] = 0) compared to the nanofluid cases ([phi] = 0.1 and 0.2). It can also be seen from the HPM results that the stream-function values decrease when the nanoparticle volume fraction [phi] increases.

[FIGURE 6 OMITTED]

[FIGURE 7 OMITTED]

[FIGURE 8 OMITTED]

Finally, the stream-function profiles f ([eta]) for various nanoparticles in [phi] = 0.1 using the HPM solution are shown in Fig. 8. It is seen that the highest stream-function values are for the Ag nanoparticles and the lowest stream-function values are for the Si[O.sub.2] nanoparticles. With regard to the stream-function values, the other types of nanoparticles fall with the order of Cu, CuO, Ti[O.sub.2] and [Al.sub.2][O.sub.3] between Ag and Si[O.sub.2] nanoparticles. It should be said that this order of variation for values of the stream-function is seen to be inversely proportional to the densities of relevant nanoparticles.

5. Conclusions

The two-dimensional boundary layer flow of nanofluids over an impermeable consciously moving horizontal plate is studied. The continuity and momentum conservation equations are transformed by the similarity method to a nonlinear ordinary differential equation which is solved using the homotopy perturbation method (HPM) for various types of nanoparticles including copper (Cu), alumina ([Al.sub.2][O.sub.3]), titania (Ti[O.sub.2]), copper oxide (CuO), silver (Ag) and silicon (Si[O.sub.2]) in the water based fluid. The results show that the present HPM solution with only two terms agrees within 2% error with the previous numerical solutions and within 1% error with the experimental data for the local skin friction coefficient. The investigation shows that the inclusion of nanoparticles in the base fluid causes an increase in the local skin friction coefficient, which also increases with the boost in the nanoparticle volume fraction. The results also show that the increase of the local skin friction coefficient depends highly on the type of nanoparticles, such that Ag nanoparticles result in the highest values of the local skin friction coefficient.

Appendix

The equation for u1 (Eq. (17)) is solved using the symbolic software Mathematica under the boundary conditions [u.sub.1] (0) = 0, [u'.sub.1] (0) = 0, i.e.:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (A.1)

which gives the following solution:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.2)

Here C (1) is integration constant and [OMEGA] = [(1 - [phi]).sup.2.5] x [1 - [phi] + [phi] ([[rho].sub.s]/ [[rho].sub.f])]. Applying the boundary condition [u'.sub.1] ([infinity]) = 0 gives:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.3)

If C(1) = 0 is substituted in u1(t) of (A.2), it gives:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.4)

By checking, it is seen that for Eq. (A.4), [u.sub.1] (0) = 0, [u'.sub.1] (0) = 0. If the third boundary condition [u'.sub.1] ([infinity]) = 0 is applied to Eq. (A.4), it gives the value of [alpha] = [+ or -] [([OMEGA] / 2).sup.0.5]. The obtained value of [alpha] removes the secular term from the ordinary differential equation (ODE) for [u.sub.1]. If [alpha] is substituted in the last term of Eq. (A.4), it gives:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.5)

crossref http://dx.doi.org/10.5755/j01.mech.20.1.3406

References

[1.] Choi, S.U.S.; Zhang, Z.G.; Yu, W.; Lockwood, F.E.; Grulke, E.A. 2001. Anomalously thermal conductivity enhancement in nanotube suspensions, Appl. Phys. Lett. 79: 2252-2254.

http://dx.doi.org/10.1063/L1408272.

[2.] Das, K. 2012. Slip flow and convective heat transfer of nanofluids over a permeable stretching surface, Computers & Fluids 64: 34-42.

http://dx.doi.org/10.1016/j.compfluid.2012.04.026.

[3.] Makinde, O.D.; Aziz, A. 2011. Boundary layer flow of a nanofluid past a stretching sheet with a convective boundary condition, Int. J. Thermal Sci. 50: 1326-1332. http://dx.doi.org/10.1016/j.ijthermalsci.2011.02.019.

[4.] Kandasamya, R.; Loganathanb, P.; Puvi-Arasu, P. 2011. Scaling group transformation for MHD boundary-layer flow of a nanofluid past a vertical stretching surface in the presence of suction/injection, Nuclear Engineering and Design 241: 2053-2059. http://dx.doi.org/10.1016/j.nucengdes.2011.04.011.

[5.] Anwar, M.I.; Khan, I.; Sharidan, S.; Salleh, M.Z. 2012. Conjugate effects of heat and mass transfer of nanofluids over a nonlinear stretching sheet, Int. J. Physical Sciences 7: 4081-4092. http://dx.doi.org/10.5897/IJPS12.358.

[6.] Bachok, N.; Ishak, A.; Pop, I. 2010. Boundary-layer flow of nanofluids over a moving surface in a flowing fluid, Int. J. Thermal Sci. 49: 1663-1668. http://dx.doi.org/10.1016/j.ijthermalsci.2010.01.026.

[7.] Bachok, N.; Ishak, A.; Pop, I. 2012. Flow and heat transfer characteristics on a moving plate in a nanofluid, Int. J. Heat Mass Transfer 55: 642-648. http://dx.doi.org/10.1016/j .ijheatmasstransfer.2011.10.047.

[8.] Ahmad, S.; Rohni, A.M.; Pop, I. 2011. Blasius and Sakiadis problems in nanofluids, Acta Mech. 218: 195-204. http://dx.doi.org/10.1007/s00707-010-0414-6.

[9.] Wang, X.Q.; Mujumdar, A.S. 2008. A review on nanofluids-Part I: theoretical and numerical investigations, Brazilian J. Chem. Engng. 25: 613-630. http://dx.doi.org/10.1590/S0104-66322008000400001.

[10.] Wang, X.Q.; Mujumdar, A.S. 2008. A review on nanofluids - Part II: experiments and applications, Brazilian J. Chem. Engng. 25: 631-648. http://dx.doi.org/10.1590/S0104-66322008000400002.

[11.] Noghrehabadi, A.; Pourrajab, R.; Ghalambaz, M. 2012. Effect of partial slip boundary condition on the flow and heat transfer of nanofluids past stretching sheet prescribed constant wall temperature, Int. J. Thermal Sciences 54: 253-261. http://dx.doi.org/10.1016/j.ijthermalsci.2011.11.017.

[12.] He, J.H. 1999. Homotopy perturbation technique, Comput. Methods Appl. Mech. Engrg. 178: 257-262. http://dx.doi.org/10.1016/S0045-7825(99)00018-3.

[13.] Raftari, B.; Mohyud-din, S.T.; Yildirim, A. 2011. Solution to the MHD flow over a non-linear stretching sheet by homotopy perturbation method, Science China Physics Mechanics & Astronomy 54: 342-345. http://dx.doi.org/10.1007/s11433-010-4180-1.

Received January 30, 2013

Accepted January 07, 2014

N. Dalir *, S.S. Nourazar **

* Department of Mechanical Engineering, Amirkabir University of Technology, Tehran, Iran, E-mail: dalir@aut.ac.ir

** Department of Mechanical Engineering, Amirkabir University of Technology, Tehran, Iran, E-mail: icp@aut.ac.ir

Table 1
Density of water and nanoparticles [10]

     Density          Fluid Phase (water)    Cu    [Al.sub.2][O.sub.3]

[rho], kg/[m.sup.3]          997.1          8933          3970

     Density          Ti[O.sub.2]   CuO     Ag     Si[O.sub.2]

[rho], kg/[m.sup.3]      4250       6500   10500      2670

Table 2
Values of [alpha] = f" (0) for various working fluids using HPM

[phi]   Cu-water working fluid    Al2O3 -water working fluid

0.0            -0.461410                  -0.461410
0.1            -0.538417                  -0.455451
0.2            -0.557023                  -0.444905

[phi]   TiO2-water working fluid   CuO-water working fluid

0.0            -0.461410                  -0.461410
0.1            -0.460196                  -0.514793
0.2            -0.447039                  -0.517323

[phi]   Ag-water working fluid     SiO2-water working fluid

0.0           -0.461410                   -0.461410
0.1           -0.577509                   -0.446562
0.2           -0.608019                   -0.412185

Table 3
Values of the local skin friction coefficient ([C.sub.f,x]
[Re.sub.x.sup.0.5]) for various working fluids using HPM,
numerical solution and experimental data

             Cu-water working fluid        Al2O3-water working fluid

         Present    Numerical   Experi-     Present     Numerical
[phi]    results    solution     mental     results      solution
                       [7]      data [8]                   [7]

 0.0    -0.461410    -0.4438     -0.455    -0.461410     -0.4438
 0.1    -0.700670    -0.6784     -0.691    -0.592702     -0.5767
 0.2    -0.973080    -0.9442     -0.961    -0.777218     -0.7410

        TiO2-water working fluid       CuO-water working fluid

         Present     Numerical      Present    Numerical   Experi-
[phi]    results      solution      results    solution     mental
                        [7]                       [7]      data [8]

 0.0    -0.461410     -0.4438      -0.461410    -0.4438     -0.455
 0.1    -0.598875     -0.5830      -0.669925
 0.2    -0.780946     -0.7540      -0.903726

        Ag-water working fluid   SiO2-water working fluid

         Present    Numerical      Present    Numerical
[phi]    results    solution       results     solution
                       [7]                       [7]

 0.0    -0.461410    -0.4438      -0.461410    -0.4438
 0.1    -0.751541                 -0.581133
 0.2    -1.062170                 -0.720057