Mathematical Problems in EngineeringVolume 2011 (2011), Article ID 258734, 17 pagesdoi:10.1155/2011/258734

Research Article

Investigation of Rotating MHD Viscous Flow and Heat Transfer between Stretching and Porous Surfaces Using Analytical Method

M. Sheikholeslami,¹ H. R. Ashorynejad,¹ D. D. Ganji,¹ and A. Kolahdooz²

¹Faculty of Mechanical Engineering, Babol University of Technology, Babol, Iran²Faculty of Mechanical Engineering, Eslamic Azad University of Technology, Khomeninishahr Branch, Iran

Received 17 May 2011; Accepted 29 July 2011

Academic Editor: Ezzat G. Bakhoum

Copyright © 2011 M. Sheikholeslami 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.

Hydromagnetic flow between two horizontal plates in a rotating system, where the lower plate is a stretching sheet and the upper is a porous solid plate, is analyzed. Heat transfer in an electrically conducting fluid bonded by two parallel plates is studied in the presence of viscous dissipation. The equations of conservation of mass and momentum and energy are reduced to a nonlinear ordinary differential equations system. Homotopy perturbation method is used to get complete analytic solution for velocity and temperature profiles. Results show an acceptable agreement between this method results and numerical solution. Also the effects of different parameters are discussed through graphs.

1. Introduction

Flow of a viscose fluid over a stretching surface has important applications in polymer industries. For instance, a number of technical processes concerning polymers involve the cooling of continuous strips extruded from a die by drawing them through a quiescent fluid with controlled cooling system, and in the process of drawing, these strips are sometimes stretched.

Glass blowing, continuous casting of metals, and spinning of fibers also involve the flow over a stretching surface. In all these cases, the quality of the final product depends on the rate of heat transfer on the stretching surface.

Dutta et al. [1] studied the temperature field in the flow over a stretching surface subjected to uniform heat flux. Andersson et al. [2] investigated the unsteady two-dimensional non-Newtonian flow of a power-law fluid past a stretching surface. Bujurke et al. [3] and Dandapat and Gupta [4] examined the temperature distribution in the steady boundary layers of a second-grade fluid near a stretching surface. P. S. Gupta and A. S. Gupta [5] investigated the heat and mass transfer on a stretching sheet with suction or blowing. Sakiadis [6] firstly studied the boundary layer flow over a stretched surface moving with constant velocity. Erickson et al. [7] extended the work of Sakiadis to include blowing or suction at the stretched sheet surface on a continuous moving surface with constant speed and investigated its effects on the heat and mass transfer in the boundary layer.

In recent years, the effect of magnetic field in different engineering applications such as the cooling of reactors and many metallurgical processes involve the cooling of continuous tiles has been under more considerable attention. Several engineering processes, such as materials manufactured by extrusion processes and heat-treated materials traveling between a feed roll and a wind-up roll on convey belts possess the characteristics of a moving continuous surface, are just some examples of applications which involve the problem discussed above.

Chakrabarti and Gupta [8] studied the MHD flow of Newtonian fluids initially at rest, over a stretching sheet at a different uniform temperature. Vajravelu and Hadjinicolaou [9] made analysis to flows and heat transfer characteristics in an electrically conducting fluid near an isothermal sheet. Heat transfer analysis of MHD fluid over a uniformly stretching sheet was investigated by Chakrabarti and Gupta [10]; In 1983, Borkakoti and Bharali [11] studied the two-dimensional channel flow with heat transfer analysis of a hydromagnetic fluid where the lower plate was a stretching sheet. The flow between two rotating disks has important technical applications such as lubrication. Keeping this fact in mind, Vajravelu and Kumar [12] studied the effect of rotation on the two-dimensional channel flow. They solved the governing equations analytically and numerically. Most of engineering problems, especially some of heat transfer equations, are nonlinear; therefore, some of them are solved using numerical solution, and some are solved using the different analytic method, such as perturbation method, homotopy perturbation method, and variational iteration method introduced by He [13, 14].

Perturbation techniques are based on the existence of small or large parameters, the so-called perturbation quantity. Unfortunately, many nonlinear problems in science and engineering do not contain those kinds of perturbation quantities. Therefore, many different methods have recently introduced some ways to eliminate the small parameter. One of the semiexact methods which does not need small parameters is the homotopy perturbation method.

The homotopy perturbation method was proposed first by He in 1998 and was further developed and improved by He [15]. The method yields a very rapid convergence of the solution series in the most of cases. The HPM proved its capability to solve a large class of nonlinear problems efficiently, accurately, and easily with approximations convergency very rapidly to solution. Usually, few iterations lead to high-accuracy solution. Recently, this method is employed for many researches in engineering sciences. He’s homotopy perturbation method is applied to obtain approximate analytical solutions for the motion of a spherical particle in a plane couette flow Jalaal et al. [16]. Then Jalaal et al. [17] showed the effectiveness of HPM for unsteady motion of a spherical particle falling in a Newtonian fluid. Ghotbi et al. [18] used HPM to approximate the solution of the ratio-dependent predator-prey system with constant effort prey harvesting. Also homotopy perturbation method was used for solving nonlinear MHD Jeffery Hamel problem by Moghimi et al. [19]. Recently, Ganji et al. studied the steady-state flow of a Hagen-Poiseuille model in a circular pipe and entropy generation due to fluid friction and heat transfer using HPM [20].

In this study, the purpose is to solve nonlinear equations through the HPM. It can be seen that this method is strongly capable of solving a large class of coupled and nonlinear differential equations without tangible restriction of sensitivity to the degree of the nonlinear term.

2. Flow Analysis

2.1. Governing Equations

Consider the steady flow of an electrically conducting fluid between two horizontal parallel plates when the fluid and the plates rotate together around the -axis which is normal to the plates with an angular velocity.

A Cartesian coordinate system is considered as followes: the -axis is along the plate, the -axis is perpendicular to it, and the -axis is normal to the plane (see Figure 1). The origin is located on the lower plate, and the plates are located atand. The lower plate is being stretched by two equal and opposite forces, so that the position of the point remains unchanged. A uniform magnetic flux with density is acting along -axis about which the system is rotating. The upper plate is subjected to a constant flow injection with a velocity. The governing equations of motion in a rotating frame of reference are where , , and denote the fluid velocity components along the ,, and directions, is the kinematic coefficient of viscosity, is the fluid density, and is the modified fluid pressure. The absence of in (2.4) implies that there is a net cross-flow along the -axis.

Figure 1: Geometry of the problem.

The boundary conditions are

The following nondimensional variables are introduced: where a prime denotes differentiation with respect to .

Substituting (2.6) in (2.1)–(2.4), we have and the nondimensional quantities are defined, in which is the viscosity parameter, is the magnetic parameter, and is the rotation parameter Equation (2.7) with the help of (2.8) can be written as follows: Differentiation of (2.11) with respect to gives

Therefore, the governing equations and boundary conditions for this case in nondimensional form are given by subject to the following boundary conditions:

3. Heat Transfer Analysis

3.1. Energy Equation

The energy equation for the present problem with viscous dissipation in nondimensional form is given by

With replacing nondimensional variables and using similarity solution method, by neglecting the last term of viscous dissipation in the energy equation, we have the following energy equation: subject to the boundary conditions where is the Prandtl number, is the Eckert number, is the local Eckert number, and the nondimensional temperature is defined as where and are temperature at the lower and upper plates.

4. Analysis of the Homotopy Perturbation Method

To illustrate the basic ideas of this method, we consider the following equation: with the boundary condition of where is a general differential operator, is a boundary operator, is a known analytical function, and is the boundary of the domain .

can be divided into two parts which are and , where is linear and is nonlinear. Equation (4.1) can therefore be rewritten as follows:

Homotopy perturbation structure is shown as follows: where

In (2.5), is an embedding parameter, and is the first approximation that satisfies the boundary condition. We can assume that the solution of (4.5) can be written as a power series in , as follows: and the best approximation for solution is

5. Implementation of the Method

According to the so-called homotopy perturbation method (HPM), we construct a homotopy. Suppose the solution of (4.1) has the following form:

We consider as follows: with substituting from (5.2) into (5.1) and some simplification and rearranging based on powers of -terms, it can be obtained that

:

:

:

By solving (5.3)–(5.5) with boundary conditions for , , , , , and , it can be obtained that

The solution of this equation, when , will be as follows: where for , , , , , and , the following functions are obtained:

6. Results and Discussion

The objective of the present study was to apply homotopy perturbation method to obtain an explicit analytic solution of rotating MHD flow and heat transfer of viscous fluid over stretching and porous surface (Figure 1). As can be seen in Table 1, homotopy perturbation method is converged in step 8, and error has been minimized. There is an acceptable agreement between the results of numerical solution obtained by four-order Rung-kutte method and differential transformation method as shown in Tables 2 and 3. In those tables, error is introduced as follows:

Table 1: values in different steps of HPM solution at,, , , , and .

NM

% Error

% Error

% Error

% Error

0
1
1
0
1
0
1
0
1
0

0.1
0.928387
0.928223
0.017678
0.928444
0.006163
0.928449
0.006698635
0.928449
0.00670429

0.2
0.846138
0.84596
0.021041
0.846272
0.01581
0.84628
0.016736093
0.84628
0.016746103

0.3
0.755385
0.755242
0.018984
0.755583
0.026195
0.755592
0.027417617
0.755592
0.027431147

0.4
0.657944
0.657829
0.017493
0.658177
0.035331
0.658186
0.036816201
0.658187
0.036832934

0.5
0.555344
0.555229
0.020556
0.555575
0.041731
0.555585
0.043510198
0.555585
0.043530456

0.6
0.44884
0.448702
0.030858
0.449039
0.044167
0.449048
0.046342095
0.449048
0.046367093

0.7
0.339433
0.339262
0.050375
0.339574
0.041441
0.339583
0.044186822
0.339583
0.044218728

0.8
0.227869
0.227685
0.080907
0.227943
0.03217
0.227951
0.035722981
0.227951
0.035765242

0.9
0.114645
0.114502
0.124586
0.114662
0.014558
0.114667
0.019207699
0.114667
0.019265704

1
0
0
0
0
0
0
0
0
0

Table 2: Comparison between numerical results and HPM solution for , , at , , , , , and .

% Error

% Error

% Error

0
0
0
0
0
0
0
1
1
0

0.1
0.094742
0.094794
0.054323
0.026758
0.026778
0.073964
0.928387
0.928449
0.006224175

0.2
0.179187
0.179354
0.093254
0.044707
0.044738
0.069201
0.846138
0.84628
0.014169522

0.3
0.253603
0.253898
0.116124
0.054988
0.055018
0.054436
0.755385
0.755592
0.020721081

0.4
0.318182
0.318574
0.123244
0.058698
0.058716
0.030566
0.657944
0.658187
0.024234016

0.5
0.373035
0.373467
0.115934
0.056892
0.056891
0.000929
0.555344
0.555585
0.02417436

0.6
0.418197
0.418601
0.096615
0.050587
0.050568
0.037868
0.44884
0.449048
0.020811422

0.7
0.453627
0.45394
0.068972
0.040777
0.040746
0.077163
0.339433
0.339583
0.015009305

0.8
0.479208
0.479391
0.038258
0.028434
0.028402
0.114496
0.227869
0.227951
0.008149803

0.9
0.494751
0.49481
0.011815
0.014521
0.0145
0.143769
0.114645
0.114667
0.002208717

1
0.5
0.5
0
0
0
0
0
0
0

Table 3: Comparison between numerical results and HPM solution for , , at , , , , , and .

% Error

% Error

% Error

0
0
0
0
0
0
0
1
1
0

0.1
0.094908
0.094961
0.055277
0.027788
0.02781
0.076047
0.9286
0.928667
0.007196834

0.2
0.179709
0.179881
0.095531
0.046609
0.046643
0.071459
0.846601
0.846752
0.017836439

0.3
0.254498
0.254803
0.119875
0.057507
0.057539
0.056888
0.75604
0.756261
0.029298263

0.4
0.319345
0.319755
0.128301
0.061531
0.061552
0.033084
0.658692
0.658953
0.039612508

0.5
0.374292
0.374748
0.121774
0.059736
0.059736
0.001376
0.556078
0.556341
0.047236899

0.6
0.419351
0.419781
0.102421
0.05317
0.053151
0.036168
0.449468
0.449697
0.050848494

0.7
0.454508
0.454843
0.073797
0.04288
0.042847
0.076493
0.339887
0.340054
0.049140518

0.8
0.479718
0.479916
0.041304
0.029902
0.029868
0.11517
0.228122
0.228215
0.040617504

0.9
0.494912
0.494976
0.012863
0.015268
0.015246
0.14576
0.114722
0.114749
0.023395295

1
0.5
0.5
0
0
0
0
0
0
0

Figure 2 shows the magnetic field effect on nondimensional velocity component (,, and ). The decrease of curve is observed by applying higher magnetic field intensity, and values increase near stretching sheet and decrease under porous sheet, while at the middle point, these values are constant.

(a)

(b)

(c)

Figure 2: Velocity components profile (a) , (b) , and (c) for variable at , , , , and .

At low Reynolds numbers, the velocity profile exhibits center line symmetry indicating a Poiseuille flow for non-Newtonian fluids. At higher Reynolds numbers, the maximum velocity point is shifted to the streiching wall where shear stress becomes larger as the Reynolds number grows.

Blowing velocity parameter has a noticable effect of nondimensional velocity component as shown in Figure 3, which by increasing profile of and becomes nonlinear, and the maximum amount of increases, and velocity component in direction increases severely.

(a)

(b)

(c)

Figure 3: Velocity components profile (a) , (b) , and (c) for variable at , , , , and .

Also it shows that increasing the blowing velocity parameter leads to increase, which shows that blowing velocity parameter and magnetic field effects on are in opposite.

Figure 4 shows that by increasing rotating parameter , values of transverse velocity component between two sheets increase, and the location of maximum amount of approaches stretching sheet. Coriolis force has inverse effect on in comparison with Lorenz force which means that with increasing the rotating parameter transverse , velocity component between two plates increases as shown in Figure 4. And it shows that the viscosity parameter affects profile similarly to magnetic field; however, with less changes in intensity, also with increasing , the location of maximum amount of approaching stretching sheet, hat indicates decreasing of boundary layer thickness near stretching plate.

(a)

(b)

Figure 4: Velocity components profile for variables and at (a) , , and (b) , , , and .

As can be seen in Figure 5, increasing viscosity parameter leads to increasing the curve of temperature profile and the decreasing of values, and it can be shown that increasing in presence of viscous dissipation leads to increasing temperature between two plates. Also increasing temperature between two plates observed, which is caused by increasing this effect, is more sensible near the stretching plate.

(a)

(b)

(c)

Figure 5: Temperature profile for variables , and at , , and and (a) and , (b) and , and (c) , .

The effects of viscous dissipation for which the Eckert number and the local Eckert number are responsible are shown in Figure 6. It is obvious from the graphs that by increasing and , the temperature near the stretching wall increases. This is due to the fact that heat energy is stored in the fluid due to the frictional heating.

(a)

(b)

Figure 6: Temperature profile for variables andat , , , , and and (a) , (b) .

In Figure 7, a coefficient of skin friction in stretching plate and porous plate is discussed with hanging effective parameters. In stretching plate, with the increase of viscosity parameter, skin friction increases, and increasing rotating parameter leads to a similar effect on skin friction. Applying higher magnetic field intensity leads to skin friction reduction. Increasing blowing velocity parameter leads to skin friction increasing. For porous plate, with the increase of , skin friction decreases, while with and increase, the reduction in skin friction observed; is also with reduction, skin friction increases.

(a)

(b)

(c)

(d)

Figure 7: The effect of active parameters on skin friction at and (a) , , and , (b) , , and , (c) , , and , (d) , , and .

A coefficient of Nusselt number consulted changes of effective parameters. Increasing the or leads to Nu decreasing, while by increasing , , and , the Nusselt number increases, as shown in Figure 8.

(a)

(b)

(c)

(d)

Figure 8: The effect of active parameters on Nusselt number . (a) , , (b) , , (c) , , (d) , , and , .

7. Conclusion

In this paper, hydromagnetic flow problem between two horizontal plates in a rotating system, where the lower plat is a stretching sheet and the upper is a porous solid plate, has been solved via a sort of analytical method, homotopy perturbation method. Also this problem is solved by a numerical method (the Runge-Kutta method of order 4), and some of the conclusions were summarized as follows:(a)homotopy perturbation method is a powerful approach for solving nonlinear differential equation such as the discussed problem, and it can be observed that there is a good agreement between the present and numerical results; (b)presence of magnetic field leads to creating a Lorentz force which causes transverse velocity component reduction between two plates although this force does not have a noticeable effect on temperature profile;(c)increasing in presence of viscous dissipation leads to temperature increasing between two plates, while in absence of viscous dissipation, the changes are inverse; (d)increasing temperature between two plates is due to increasing viscosity parameter or increasing viscous dissipation, whose effect is more sensible near stretching plate;(e)increasing magnetic parameter or viscosity parameter leads to decreasing Nu, while with increasing the rotation parameter, blowing velocity parameter, and , the Nusselt number increases.

Nomenclature

: | Magnetic field (wb·m⁻²)
: | Magnetic parameter
: | Constant function
: | Viscosity parameter
: | Rotation parameter
: | Modified fluid (pressure)
: | Prandtl number
: | Eckert number
: | Injection velocity
: | Maximum value of velocity
: | Velocity components along , , and axes, respectively.

Greek Symbols

: | Kinematic viscosity
: | Angle of the channel
: | Any angle
: | Dimensionless angle
: | Fluid density.

References

B. K. Dutta, P. Roy, and A. S. Gupta, “Temperature field in flow over a stretching sheet with uniform heat flux,” International Communications in Heat and Mass Transfer, vol. 12, no. 1, pp. 89–94, 1985.

H. I. Andersson, J. B. Aarseth, N. Braud, and B. S. Dandapat, “Flow of a power-law fluid film on an unsteady stretching surface,” Journal of Non-Newtonian Fluid Mechanics, vol. 62, no. 1, pp. 1–8, 1996.

N. M. Bujurke, S. N. Biradar, and P. S. Hiremath, “Second-order fluid flow past a stretching sheet with heat transfer,” Zeitschrift für Angewandte Mathematik und Physik, vol. 38, no. 4, pp. 653–657, 1987.

B. S. Dandapat and A. S. Gupta, “Flow and heat transfer in a viscoelastic fluid over a stretching sheet,” International Journal of Non-Linear Mechanics, vol. 24, no. 3, pp. 215–219, 1989.

P. S. Gupta and A. S. Gupta, “Heat and mass transfer on a stretching sheet with suction and blowing,” Journal of Chemical and Engineering Data, vol. 55, pp. 744–745, 1977.

B. C. Sakiadis, “Boundary layer behavior on continuous solid flat surfaces,” Aiche Journal, vol. 7, pp. 26–28, 1961.

L. E. Erickson, L. T. Fan, and V. G. Fox, “Heat and mass transfer on a moving continuous flat plate with suction or injection,” Industrial and Engineering Chemistry Fundamentals, vol. 5, no. 1, pp. 19–25, 1966.

A. Chakrabarti and A. S. Gupta, “Hydromagnetic flow and heat transfer over a stretching sheet,” Quarterly of Applied Mathematics, vol. 37, no. 1, pp. 73–78, 1979.

K. Vajravelu and A. Hadjinicolaou, “Convective heat transfer in an electrically conducting fluid at a stretching surface with uniform free stream,” International Journal of Engineering Science, vol. 35, no. 13-12, pp. 1237–1244, 1997.

M. Subhas Abel and N. Mahesha, “Heat transfer in MHD viscoelastic fluid flow over a stretching sheet with variable thermal conductivity, non-uniform heat source and radiation,” Applied Mathematical Modelling, vol. 32, no. 10, pp. 1965–1983, 2008.

A. K. Borkakoti and A. Bharali, “Hydromagnetic flow and heat transfer between two horizontal plates, the lower plate being a stretching sheet,” Quarterly of Applied Mathematics, vol. 41, no. 4, pp. 461–467, 1983.

K. Vajravelu and B. V. R. Kumar, “Analytical and numerical solutions of a coupled non-linear system arising in a three-dimensional rotating flow,” International Journal of Non-Linear Mechanics, vol. 39, no. 1, pp. 13–24, 2004.

A. H. Nayfeh, Introduction to Perturbation Techniques, Wiley, New York, NY, USA, 1979.

J.-H. He, “A coupling method of a homotopy technique and a perturbation technique for non-linear problems,” International Journal of Non-Linear Mechanics, vol. 35, no. 1, pp. 37–43, 2000.

J. H. He, “Some asymptotic methods for strongly nonlinear equations,” International Journal of Modern Physics B, vol. 20, no. 10, pp. 1141–1199, 2006.

M. Jalaal, M. G. Nejad, P. Jalili et al., “Homotopy perturbation method for motion of a spherical solid particle in plane couette fluid flow,” Computers and Mathematics with Applications, vol. 61, no. 8, pp. 2267–2270, 2011.

M. Jalaal, D. D. Ganji, and G. Ahmadi, “Analytical investigation on acceleration motion of a vertically falling spherical particle in incompressible Newtonian media,” Advanced Powder Technology, vol. 21, no. 3, pp. 298–304, 2010.

A. R. Ghotbi, A. Barari, and D. D. Ganji, “Solving ratio-dependent predator-prey system with constant effort harvesting using homotopy perturbation method,” Mathematical Problems in Engineering, Article ID 945420, 8 pages, 2008.

S. M. Moghimi, D. D. Ganji, H. Bararnia, M. Hosseini, and M. Jalaal, “Homotopy perturbation method for nonlinear MHD Jeffery-Hamel problem,” Computers and Mathematics with Applications, vol. 61, no. 8, pp. 2213–2216, 2011.

D. D. Ganji, H. R. Ashory Nezhad, and A. Hasanpour, “Effect of variable viscosity and viscous dissipation on the Hagen-Poiseuille flow and entropy generation,” Numerical Methods for Partial Differential Equations, vol. 27, no. 3, pp. 529–540, 2011.