Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2011 / Article

Research Article | Open Access

Volume 2011 |Article ID 258734 | https://doi.org/10.1155/2011/258734

M. Sheikholeslami, H. R. Ashorynejad, D. D. Ganji, A. Kolahdooz, "Investigation of Rotating MHD Viscous Flow and Heat Transfer between Stretching and Porous Surfaces Using Analytical Method", Mathematical Problems in Engineering, vol. 2011, Article ID 258734, 17 pages, 2011. https://doi.org/10.1155/2011/258734

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

Academic Editor: Ezzat G. Bakhoum
Received17 May 2011
Accepted29 Jul 2011
Published13 Oct 2011

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 at𝑦=0and𝑦=ℎ. The lower plate is being stretched by two equal and opposite forces, so that the position of the point (0,0,0) remains unchanged. A uniform magnetic flux with density 𝐵0 is acting along 𝑦-axis about which the system is rotating. The upper plate is subjected to a constant flow injection with a velocity𝑣0. The governing equations of motion in a rotating frame of reference are𝜕𝑢+𝜕𝑥𝜕𝑣+𝜕𝑦𝜕𝑤𝑢𝜕𝑧=0,(2.1)𝜕𝑢𝜕𝑥+𝜈𝜕𝑢1𝜕𝑦+2Ω𝑤=−𝜌𝜕𝑝∗𝜕𝜕𝑥+𝜐2𝑢𝜕𝑥2+𝜕2𝑢𝜕𝑦2î‚¹âˆ’ğœŽğµ20𝜌𝑢𝑢,(2.2)𝜕𝑣1𝜕𝑦=−𝜌𝜕𝑝∗𝜕𝜕𝑦+𝜐2𝑣𝜕𝑥2+𝜕2𝑣𝜕𝑦2𝑢,(2.3)𝜕𝑤𝜕𝑥+𝜈𝜕𝑤𝜕𝜕𝑦−2Ω𝑤=𝜐2𝑤𝜕𝑥2+𝜕2𝑤𝜕𝑦2î‚¹âˆ’ğœŽğµ20𝜌𝑤,(2.4) 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.

The boundary conditions are𝑢=ğ‘Žğ‘¥,𝑣=0,𝑤=0at𝑦=0,𝑢=0,𝑣=−𝑣0,𝑤=0at𝑦=+ℎ.(2.5)

The following nondimensional variables are introduced:𝑦𝜂=ℎ,𝑢=ğ‘Žğ‘¥ğ‘“î…ž(𝜂),𝜈=âˆ’ğ‘Žâ„Žğ‘“(𝜂),𝑤=ğ‘Žğ‘¥ğ‘”(𝜂),(2.6) where a prime denotes differentiation with respect to 𝜂.

Substituting (2.6) in (2.1)–(2.4), we have−1ğœŒâ„Žğœ•ğ‘âˆ—ğœ•ğœ‚=ğ‘Ž2ğ‘¥î‚¸ğ‘“î…žâˆ’ğ‘“ğ‘“î…žî…žâˆ’ğ‘“î…žî…žî…žğ‘…+𝑀𝑅+2𝐾𝑟𝑅𝑔−1,(2.7)ğœŒâ„Žğœ•ğ‘âˆ—ğœ•ğœ‚=ğ‘Ž2â„Žî‚ƒğ‘“ğ‘“î…ž+1ğ‘…ğ‘“î…žî…žî‚„ğ‘”,(2.8)î…žî…žî€·ğ‘“âˆ’ğ‘…î…žğ‘”âˆ’ğ‘“ğ‘”î…žî€¸+2𝐾𝑟𝑓′−𝑀𝑔=0,(2.9) and the nondimensional quantities are defined, in which 𝑅 is the viscosity parameter, 𝑀 is the magnetic parameter, and 𝐾𝑟 is the rotation parameter𝑅=ğ‘Žâ„Ž2𝜐,𝑀=ğœŽğµ20ℎ2𝜌𝜐,𝐾𝑟=Ωℎ2𝜐.(2.10) Equation (2.7) with the help of (2.8) can be written as follows:ğ‘“î…žî…žî…žî‚ƒğ‘“âˆ’ğ‘…î…ž2âˆ’ğ‘“ğ‘“î…žî…žî‚„âˆ’2𝐾2𝑟𝑔−𝑀2ğ‘“î…ž=𝐴.(2.11) Differentiation of (2.11) with respect to 𝜂 givesğ‘“ğ‘–ğ‘£î€·âˆ’ğ‘…ğ‘“â€²ğ‘“î…žî…žâˆ’ğ‘“ğ‘“î…žî…žî€¸âˆ’2ğ¾ğ‘Ÿğ‘”â€²âˆ’ğ‘€ğ‘“î…žî…ž=0.(2.12)

Therefore, the governing equations and boundary conditions for this case in nondimensional form are given byğ‘“ğ‘–ğ‘£î€·âˆ’ğ‘…ğ‘“â€²ğ‘“î…žî…žâˆ’ğ‘“ğ‘“î…žî…žî€¸âˆ’2ğ¾ğ‘Ÿğ‘”â€²âˆ’ğ‘€ğ‘“î…žî…žğ‘”=0,î…žî…žâˆ’ğ‘…(𝑓′𝑔−𝑓𝑔′)+2𝐾𝑟𝑓′−𝑀𝑔=0,(2.13) subject to the following boundary conditions:𝑓=0,ğ‘“î…ž=1,𝑔=0at𝜂=0,𝑓=𝜆,ğ‘“î…žğ‘£=0,𝑔=0at𝜂=1,𝜆=0.ğ‘Žâ„Ž(2.14)

3. Heat Transfer Analysis

3.1. Energy Equation

The energy equation for the present problem with viscous dissipation in nondimensional form is given by𝑢𝜕𝑇𝜕𝑥+𝑣𝜕𝑇𝜕𝑦+𝑤𝜕𝑇=𝜕𝑧𝑘𝜌𝑐𝑝𝜕2𝑇𝜕𝑥2+𝜕2𝑇𝜕𝑦2+𝜕2𝑇𝜕𝑧2+𝜇𝜑,𝜙=2𝜕𝑢𝜕𝑥2+𝜕𝑣𝜕𝑦2+𝜕𝑤𝜕𝑧2+𝜕𝑣+𝜕𝑥𝜕𝑢𝜕𝑦2+𝜕𝑤+𝜕𝑦𝜕𝑣𝜕𝑧2+𝜕𝑤+𝜕𝑥𝜕𝑢𝜕𝑧2−23𝜕𝑢+𝜕𝑥𝜕𝑣+𝜕𝑦𝜕𝑤𝜕𝑧2.(3.1)

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:ğœƒî…žî…žî€ºî€·+Pr𝑅𝑓𝜃′+Ec4ğ‘“î…ž2+𝑔2+Ecğ‘¥î€·ğ‘“î…žî…ž2+ğ‘”î…ž2=0,(3.2) subject to the boundary conditions𝜃(0)=1,𝜃(1)=0,(3.3) where Pr=𝜇𝐶𝑝/𝑘 is the Prandtl number, Ec=ğ‘Ž2ℎ2/𝐶𝑝(𝜃0âˆ’ğœƒâ„Ž) is the Eckert number, Ec𝑥=ğ‘Ž2𝑥2/𝐶𝑝(𝜃0âˆ’ğœƒâ„Ž) is the local Eckert number, and the nondimensional temperature is defined as𝜃(𝜂)=ğ‘‡âˆ’ğ‘‡â„Žğ‘‡0âˆ’ğ‘‡â„Ž,(3.4) where 𝑇0 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:𝐴(𝑢)−𝑓(𝑟)=0𝑟∈Ω,(4.1) with the boundary condition of𝐵𝑢,𝜕𝑢𝜕𝑛=0𝑟∈Γ,(4.2) 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:𝐿(𝑢)+𝑁(𝑢)−𝑓(𝑟)=0𝑟∈Ω.(4.3)

Homotopy perturbation structure is shown as follows:𝐻𝐿𝑢(𝜈,𝑝)=(1−𝑝)(𝜈)−𝐿0[𝐴]+𝑝(𝜈)−𝑓(𝑟)=0,(4.4) where[]𝜈(𝑟,𝑝)∶Ω×0,1⟶𝑅.(4.5)

In (2.5), 𝑝∈[0,1] is an embedding parameter, and 𝑢0 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:𝜈=𝜈0+𝑝𝜈1+𝑝2𝜈2+⋯,(4.6) and the best approximation for solution is𝑢=lim𝑝→1𝜈=𝜈0+𝜈1+𝜈2+⋯.(4.7)

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:𝐻𝑓(𝑓,𝑝)=(1−𝑝)𝑖𝑣−𝑓0𝑖𝑣𝑓+ğ‘ğ‘–ğ‘£î€·ğ‘“âˆ’ğ‘…î…žğ‘“î…žî…žâˆ’ğ‘“ğ‘“î…žî…žî€¸âˆ’2ğ¾ğ‘Ÿğ‘”î…žâˆ’ğ‘€ğ‘“î…žî…žî€¸î€·ğ‘”=0,𝐻(𝑔,𝑝)=(1−𝑝)î…žî…žâˆ’ğ‘”0î…žî…žî€¸î€·ğ‘”+ğ‘î…žî…žî€·ğ‘“âˆ’ğ‘…î…žğ‘”âˆ’ğ‘“ğ‘”î…žî€¸+2ğ¾ğ‘Ÿğ‘“î…žî€¸î€·ğœƒâˆ’ğ‘€ğ‘”=0,𝐻(𝜃,𝑝)=(1−𝑝)î…žî…žâˆ’ğœƒ0î…žî…žî€¸î‚€ğœƒ+ğ‘î…žî…žî‚ƒ+Prğ‘…ğ‘“ğœƒî…žî‚€+Ec4ğ‘“î…ž2+𝑔2+Ec𝑥𝑓2+ğ‘”î…ž2=0.(5.1)

We consider 𝑓,𝑔,𝜃 as follows:𝑓(𝜂)=𝑓0(𝜂)+𝑝𝑓1(𝜂)+⋯=𝑛𝑖=0𝑝𝑖𝑓𝑖(𝜂),𝑔(𝜂)=𝑔0(𝜂)+𝑝𝑔1(𝜂)+⋯=𝑛𝑖=0𝑝𝑖𝑔𝑖(𝜂),𝜃(𝜂)=𝜃0(𝜂)+𝑝𝜃1(𝜂)+⋯=𝑛∑𝑖=0𝑝𝑖𝜃𝑖.(𝜂)(5.2) with substituting 𝑓,𝑔,𝜃 from (5.2) into (5.1) and some simplification and rearranging based on powers of 𝑝-terms, it can be obtained that

𝑝0:𝑓0𝑖𝑣𝑔=0,0î…žî…žğœƒ=0,0=0.Andboundaryconditionsare𝑓(0)=0,𝑓′(0)=1,𝑓(1)=𝜆,𝑓′(1)=0,𝑔(0)=1,𝑔′(1)=0,𝜃(0)=1,𝜃(1)=0,(5.3)

𝑝1:âˆ’ğ‘”î…ž0−0.5𝑓0+0.5𝑓0𝑓0−0.5ğ‘“î…ž0𝑓0+𝑓1𝑖𝑣=0,−0.5ğ‘“î…ž0𝑔0−0.5𝑔0+𝑔1+0.5𝑓0ğ‘”î…ž0+ğ‘“î…ž0𝑓=0,0.2502𝑔+0.2502𝑔+0.2502+𝜃1+0.25𝑓0ğœƒî…ž1+î€·ğ‘“î…ž02=0.Andboundaryconditionsare𝑓(0)=0,𝑓′(0)=0,𝑓(1)=0,𝑓′(1)=0,𝑔(0)=0,𝑔′(1)=0,𝜃(0)=0,𝜃(1)=0,(5.4)

𝑝2:+0.5𝑓0𝑓1î…žî…žâˆ’ğ‘”î…ž0+0.5𝑓1𝑓0+0.5𝑓1î…žğ‘“0−0.5𝑓1+𝑓2𝑖𝑣=0,0.5𝑓1ğ‘”î…ž0+0.5𝑓0ğ‘”î…ž1+ğ‘“î…ž0−0.5𝑔1−0.5ğ‘“î…ž0𝑔1+𝑔2−0.5ğ‘“î…ž1𝑔0=0,0.25𝑓0ğœƒî…ž1+0.25𝑓1ğœƒî…ž0𝑓+0.2502+𝜃2+2ğ‘“î…ž0ğ‘“î…ž1+0.5𝑓0î…žî…žğ‘“1+0.5𝑔0𝑔1+0.5ğ‘”î…ž0ğ‘”î…ž1=0.Andboundaryconditionsare𝑓(0)=0,𝑓′(0)=0,𝑓(1)=0,𝑓′(1)=0,𝑔(0)=0,𝑔′(1)=0,𝜃(0)=0,𝜃(1)=0.(5.5)

By solving (5.3)–(5.5) with boundary conditions for 𝑅=0.5, 𝐾𝑟=0.5, 𝑀=0.5, Pr=0.5, 𝜆=0.5, and Ec=Ec𝑥=0.5, it can be obtained that𝑓0(𝜂)=−0.5𝜂2𝑔+𝜂,0(𝜃𝜂)=0,0𝑓(𝜂)=−𝜂+1,1(𝜂)=−0.0007𝜂6+0.0083𝜂5−0.0416𝜂4+0.0611𝜂3−0.0270𝜂2,𝑔1(𝜂)=0.1667𝜂3−0.5𝜂2𝜃+0.3333𝑥,1(𝜂)=−0.0937𝜂4+0.3750𝜂3−0.6250𝜂2𝑓+0.3437𝜂,2(𝜂)=0.0002𝜂8+0.0010𝜂7−0.0116𝜂6−0.0046𝜂5+0.0116𝜂4−0.0081𝜂3+0.0016𝜂2,𝑔2(𝜂)=0.0041𝜂5−0.0222𝜂4+0.0368𝜂3𝜃−0.0188𝜂,2(𝜂)=−0.0001𝜂8+0.0011𝜂7−0.0067𝜂6+0.0180𝜂5−0.0313𝜂4+0.0342𝜂3−0.0135𝜂2−0.0016𝜂.(5.6)

The solution of this equation, when 𝑝→1, will be as follows:𝑓(𝜂)=𝑓0(𝜂)+𝑓1(𝜂)+⋯+𝑓7(𝜂)+𝑓8𝑔(𝜂),(𝜂)=𝑔0(𝜂)+𝑔1(𝜂)+⋯+𝑔7(𝜂)+𝑔8(𝜂),𝜃(𝜂)=𝜃0(𝜂)+𝜃1(𝜂)+⋯+𝜃7(𝜂)+𝜃8(𝜂).(5.7) where for 𝑅=0.5, 𝐾𝑟=0.5, 𝑀=0.5, Pr=0.5, 𝜆=0.5, and Ec=Ec𝑥=0.5, the following functions are obtained:𝑓(𝜂)=−0.0001𝜂8+0.0007𝜂7−0.0012𝜂6+0.0033𝜂5−0.0306𝜂4+0.0537𝜂3−0.5256𝜂2+𝜂,𝑔(𝜂)=−0.0001𝜂8+0.0002𝜂7+0.0009𝜂6+0.0019𝜂5−0.0203𝜂4+0.2015𝜂3−0.5𝜂2+0.3158𝜂,𝜃(𝜂)=0.0001𝜂9−0.0005𝜂8+0.0015𝜂7−0.0093𝜂6+0.0293𝜂5−0.1481𝜂4+0.4322𝜂3−0.6506𝜂2−0.6546𝜂+1.(5.8)

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: ||||%Error=𝑓(𝜂)NM−𝑓(𝜂)HPM𝑓(𝜂)NM||||×100.(6.1)


𝜂 NM 𝑛 = 2 % Error 𝑛 = 4 % Error 𝑛 = 6 % Error 𝑛 = 8 % Error

0110101010
0.10.9283870.9282230.0176780.9284440.0061630.9284490.0066986350.9284490.00670429
0.20.8461380.845960.0210410.8462720.015810.846280.0167360930.846280.016746103
0.30.7553850.7552420.0189840.7555830.0261950.7555920.0274176170.7555920.027431147
0.40.6579440.6578290.0174930.6581770.0353310.6581860.0368162010.6581870.036832934
0.50.5553440.5552290.0205560.5555750.0417310.5555850.0435101980.5555850.043530456
0.60.448840.4487020.0308580.4490390.0441670.4490480.0463420950.4490480.046367093
0.70.3394330.3392620.0503750.3395740.0414410.3395830.0441868220.3395830.044218728
0.80.2278690.2276850.0809070.2279430.032170.2279510.0357229810.2279510.035765242
0.90.1146450.1145020.1245860.1146620.0145580.1146670.0192076990.1146670.019265704
1000000000


𝜂 𝑓 𝑔 𝜃
N M H P M % Error N M H P M % Error N M H P M % Error

0000000110
0.10.0947420.0947940.0543230.0267580.0267780.0739640.9283870.9284490.006224175
0.20.1791870.1793540.0932540.0447070.0447380.0692010.8461380.846280.014169522
0.30.2536030.2538980.1161240.0549880.0550180.0544360.7553850.7555920.020721081
0.40.3181820.3185740.1232440.0586980.0587160.0305660.6579440.6581870.024234016
0.50.3730350.3734670.1159340.0568920.0568910.0009290.5553440.5555850.02417436
0.60.4181970.4186010.0966150.0505870.0505680.0378680.448840.4490480.020811422
0.70.4536270.453940.0689720.0407770.0407460.0771630.3394330.3395830.015009305
0.80.4792080.4793910.0382580.0284340.0284020.1144960.2278690.2279510.008149803
0.90.4947510.494810.0118150.0145210.01450.1437690.1146450.1146670.002208717
10.50.50000000


𝜂 𝑓 𝑔 𝜃
N M H P M % Error N M H P M % Error N M H P M % Error

0000000110
0.10.0949080.0949610.0552770.0277880.027810.0760470.92860.9286670.007196834
0.20.1797090.1798810.0955310.0466090.0466430.0714590.8466010.8467520.017836439
0.30.2544980.2548030.1198750.0575070.0575390.0568880.756040.7562610.029298263
0.40.3193450.3197550.1283010.0615310.0615520.0330840.6586920.6589530.039612508
0.50.3742920.3747480.1217740.0597360.0597360.0013760.5560780.5563410.047236899
0.60.4193510.4197810.1024210.053170.0531510.0361680.4494680.4496970.050848494
0.70.4545080.4548430.0737970.042880.0428470.0764930.3398870.3400540.049140518
0.80.4797180.4799160.0413040.0299020.0298680.115170.2281220.2282150.040617504
0.90.4949120.4949760.0128630.0152680.0152460.145760.1147220.1147490.023395295
10.50.50000000

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.

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.

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.

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 Pr 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.

The effects of viscous dissipation for which the Eckert number (Ec) and the local Eckert number (Ec𝑥) are responsible are shown in Figure 6. It is obvious from the graphs that by increasing Ec and Ec𝑥, 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.

In Figure 7, a coefficient of skin friction in stretching plate (ğ‘“î…žî…ž(0)) and porous plate (ğ‘“î…žî…ž(1)) 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 coefficient of Nusselt number (ğœƒî…ž(0)) consulted changes of effective parameters. Increasing the 𝑀 or 𝑅 leads to Nu decreasing, while by increasing 𝐾𝑟, Pr, and 𝜆, the Nusselt number increases, as shown in Figure 8.

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 Pr 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 Pr, the Nusselt number increases.

Nomenclature

𝐵0:Magnetic field (wb·m−2)
𝑀:Magnetic parameter
𝐶𝑖:Constant function
𝑅:Viscosity parameter
𝐾𝑟:Rotation parameter
𝑝∗:Modified fluid (pressure)
Pr:Prandtl number
Ec:Eckert number
𝑣0:Injection velocity
𝑈max: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

  1. 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. View at: Google Scholar
  2. 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. View at: Publisher Site | Google Scholar
  3. 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. View at: Publisher Site | Google Scholar
  4. 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. View at: Google Scholar
  5. 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. View at: Google Scholar
  6. B. C. Sakiadis, “Boundary layer behavior on continuous solid flat surfaces,” Aiche Journal, vol. 7, pp. 26–28, 1961. View at: Google Scholar
  7. 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. View at: Google Scholar
  8. 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. View at: Google Scholar
  9. 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. View at: Google Scholar
  10. 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. View at: Google Scholar
  11. 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. View at: Google Scholar
  12. 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. View at: Publisher Site | Google Scholar
  13. A. H. Nayfeh, Introduction to Perturbation Techniques, Wiley, New York, NY, USA, 1979.
  14. 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  15. J. H. He, “Some asymptotic methods for strongly nonlinear equations,” International Journal of Modern Physics B, vol. 20, no. 10, pp. 1141–1199, 2006. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  16. 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. View at: Publisher Site | Google Scholar
  17. 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. View at: Publisher Site | Google Scholar
  18. 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  19. 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  20. 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. View at: Publisher Site | Google Scholar

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.


More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views1094
Downloads939
Citations

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.