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

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 𝐵₀ 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𝜕𝑢𝑢𝜕𝑢𝜕𝑥+𝜕𝜕𝑣𝑦+𝜕𝜕𝑤𝑧=0,(2.1)𝜕𝑥+𝜈𝜕𝜕𝑦𝑢+2Ω𝑤=−𝜌1𝜕𝜕𝑝𝑥^∗+𝜐𝜕²𝑢𝜕𝑥²+𝜕²𝑢𝑢𝜕𝑣𝜕𝑦²−𝜎𝐵𝜌²⁰𝑢,(2.2)𝜕𝑦=−𝜌1𝜕𝜕𝑝𝑦^∗+𝜐𝜕²𝑣𝜕𝑥²+𝜕²𝑣𝑢𝜕𝑤𝜕𝑦²,(2.3)𝜕𝑥+𝜈𝜕𝜕𝑤𝑦−2Ω𝑤=𝜐𝜕²𝑤𝜕𝑥²+𝜕²𝑤𝜕𝑦²−𝜎𝐵𝜌²⁰𝑤,(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.

Figure 1 

Geometry of the problem.

The boundary conditions are𝑢=𝑢=0,𝑎𝑥𝑣,=𝑣−=𝑣₀0,𝑤=0at𝑦=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𝜌ℎ𝜕𝜕𝑝𝜂^∗=𝑎²𝑥𝑓^−𝑓𝑓^−𝑓^−1𝑅+𝑀𝑅+2𝑅𝐾^𝑟𝑔,(2.7)𝜌ℎ𝜕𝜕𝑝𝜂^∗=𝑎²ℎ𝑓𝑓^+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.10) Equation (2.7) with the help of (2.8) can be written as follows:𝑓^−𝑅𝑓^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,𝑎ℎ.(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𝜕𝑢𝜕𝑧²+𝜇𝜑,𝜕𝑥²+𝜕𝑣𝜕𝑦²+𝜕𝑤𝜕𝑧²+𝜕𝜕𝑥𝑣+𝜕𝜕𝑦𝑢²+𝜕𝑤𝜕𝑦+𝜕𝜕𝑣𝑧²+𝜕𝑤𝜕𝑥+𝜕𝜕𝑧𝑢²−23𝜕𝜕𝑥𝑢+𝜕𝜕𝑣𝑦+𝜕𝜕𝑤𝑧².(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𝑅𝑓𝜃′+Ec4𝑓^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=𝑎²ℎ²/𝐶_𝑝(𝜃₀−𝜃_ℎ) is the Eckert number, Ec_𝑥=𝑎²𝑥²/𝐶_𝑝(𝜃₀−𝜃_ℎ) is the local Eckert number, and the nondimensional temperature is defined as𝜃(𝜂)=𝑇−𝑇^ℎ𝑇₀−𝑇_ℎ,(3.4) 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:𝐴(𝑢)−𝑓(𝑟)=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,(4.4) where𝜈(𝑟,𝑝)∶Ω×[0,1]⟶𝑅.(4.5)

In (2.5), 𝑝∈[0,1] 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:𝜈=𝜈₀+𝑝𝜈₁+𝑝²𝜈₂+⋯,(4.6) and the best approximation for solution is𝑢=lim_𝑝→1𝜈=𝜈₀+𝜈₁+𝜈₂+⋯.(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:𝐻𝐻((𝑓𝑔,,𝑝𝑝))==((11−−𝑝𝑝))𝑔𝑓^𝑖𝑣_−𝑓^0𝑖𝑣+𝑝𝑓^𝑖𝑣−𝑅𝑓^𝑓^−𝑓𝑓^−2𝐾_𝑟𝑔^−𝑀𝑓^=0,𝐻(𝜃,𝑝)=(1−𝑝)𝜃_−𝑔^0+𝑝𝑔^−𝑅𝑓^𝑔−𝑓𝑔^+2𝐾_𝑟𝑓^−𝑀𝑔=0,−𝜃_0+𝑝𝜃^+Pr𝑅𝑓𝜃^+Ec4𝑓^2+𝑔²+Ec^𝑥𝑓^2+𝑔^2=0.(5.1)

We consider 𝑓,𝑔,𝜃 as follows:𝑓(𝜂)=𝑓₀(𝜂)+𝑝𝑓₁(𝜂)+⋯=^𝑛_𝑖=0𝑔(𝜂)=𝑔₀𝑝^𝑖𝑓_𝑖(𝜂),(𝜂)+𝑝𝑔₁(𝜂)+⋯=^𝑛_𝑖=0𝜃(𝜂)=𝜃₀𝑝^𝑖𝑔_𝑖(𝜂),(𝜂)+𝑝𝜃₁(𝜂)+⋯=∑^𝑛_𝑖=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,An𝑓𝑔𝜃d′(((b000o)))u===nd111a,,,ry=co𝑔0𝑓𝜃n′.(d((111i)t))io===ns0𝜆0,a,,re𝑓′(1)=0,(5.3)

𝑝¹:−−0𝑔.^05−𝑓_00.5𝑓^0+0.5𝑓₀𝑓^0−0.5𝑓^0𝑓^0+𝑓^1𝑖𝑣=0,0.25𝑓_0𝑔₀−0.5𝑔₀+𝑔^1+0.5𝑓₀𝑔^0+𝑓^0=0,𝑓(²0)+=0.02,5A𝑔^0n𝑓𝑔d𝜃²′(((+b000o)))0u===.n2d5000a,,,r𝑔y₀c²o𝑔𝜃𝑓+n′(d((1𝜃11i)t))i^1o===n+0s00,0,,a.r2e5𝑓𝑓₀′𝜃(^11+)=𝑓0^0,²=0.(5.4)

𝑝²:0.5+𝑓0₁.5𝑓₀𝑓^1−𝑔^0+0.5𝑓₁𝑓^0+0.5𝑓_1𝑓^0−0.5𝑓^1+𝑓^2𝑖𝑣=0,0.25𝑓₀𝑔^0+0.5𝑓₀𝑔^1+𝑓^0−0.5𝑔₁−0.5𝑓^0𝑔₁+𝑔^2−0.5𝑓^1𝑔₀=0,𝜃^1+0.25𝑓₁𝜃𝑓^0(+0)0=.250A,𝑓nd𝑓^0𝑔𝜃′b²(((0o00+)u))n𝜃===d^2a000+r,,,y2c𝑓𝑔𝜃𝑓o′^0(n𝑓((111d^1)))i+t=i==o0n000.s.5,,𝑓ar𝑓^0e𝑓′(^11+)=0.50𝑔,₀𝑔₁+0.5𝑔^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.5𝜂²+𝜂,₀𝜃(𝜂)=0,₀𝑓(𝜂)=−𝜂+1,₁𝑔₁(𝜂)=−0.0007𝜂⁶+0.0083𝜂⁵−0.0416𝜂⁴+0.0611𝜂³−0.0270𝜂²,𝜃(𝜂)=0.1667𝜂³−0.5𝜂²+0.3333𝑥,₁𝑓(𝜂)=−0.0937𝜂⁴+0.3750𝜂³−0.6250𝜂²+0.3437𝜂,₂𝑔₂(𝜂)=0.0002𝜂⁸+0.0010𝜂⁷−0.0116𝜂⁶−0.0046𝜂⁵+0.0116𝜂⁴−0.0081𝜂³+0.0016𝜂²,𝜃(𝜂)=0.0041𝜂⁵−0.0222𝜂⁴+0.0368𝜂³−0.0188𝜂,₂(𝜂)=−0.0001𝜂⁸+0.0011𝜂⁷−0.0067𝜂⁶+0.0180𝜂⁵−0.0313𝜂⁴+0.0342𝜂³−0.0135𝜂²−0.0016𝜂.(5.6)

The solution of this equation, when 𝑝→1, will be as follows:𝑓𝑔((𝜂𝜂))==𝑓𝑔₀₀(𝜂)+𝑓₁(𝜂)+⋯+𝑓₇(𝜂)+𝑓₈(𝜂),𝜃(𝜂)=𝜃₀(𝜂)+𝑔₁(𝜂)+⋯+𝑔₇(𝜂)+𝑔₈(𝜂),(𝜂)+𝜃₁(𝜂)+⋯+𝜃₇(𝜂)+𝜃₈(𝜂).(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:𝑓𝑔((𝜂𝜂))==−−00..00000011𝜂𝜂⁸₈+0.0007𝜂⁷−0.0012𝜂⁶+0.0033𝜂⁵−0.0306𝜂⁴+0.0537𝜂³−0.5256𝜂²+𝜂,𝜃(𝜂)=0.0001𝜂₉+0.0002𝜂⁷+0.0009𝜂⁶+0.0019𝜂⁵−0.0203𝜂⁴+0.2015𝜂³−0.5𝜂²+0.3158𝜂,−0.0005𝜂⁸+0.0015𝜂⁷−0.0093𝜂⁶+0.0293𝜂⁵−0.1481𝜂⁴+0.4322𝜂³−0.6506𝜂²−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)

Table 1 

𝜃(𝜂) values in different steps of HPM solution at𝑅=0.5,𝐾_𝑟=0.5, 𝑀=0.5, Pr=0.5, 𝜆=0.5, and Ec=Ec_𝑥=0.5.

Table 2 

Comparison between numerical results and HPM solution for 𝑓, 𝑔, 𝜃 at 𝑅=0.5, 𝐾_𝑟=0.5, 𝑀=0.5, Pr=0.5, 𝜆=0.5, and Ec=Ec_𝑥=0.5.

Table 3 

Comparison between numerical results and HPM solution for 𝑓, 𝑔, 𝜃 at 𝑅=0.5, 𝐾_𝑟=0.5, 𝑀=0, Pr=0.5, 𝜆=0.5, and Ec=Ec_𝑥=0.5.

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 𝑅=2, 𝐾𝑟=0.5, Pr=1, 𝜆=0.5, and Ec=Ec𝑥=0.5.

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 𝑅=2, 𝐾𝑟=0.5, 𝑀=1, Pr=1, and Ec=Ec𝑥=0.5.

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) 𝑅=2, 𝑀=1, 𝜆=0.5 and (b) 𝐾𝑟=4, 𝑀=1, 𝜆=0.5, and Ec=Ec𝑥=0.5.

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.

(a)

(b)

(c)

Figure 5 

Temperature profile (𝜃) for variables 𝑅, Pr and 𝜆 at 𝐾𝑟=0.5, 𝑀=1, and E𝑐=Ec𝑥=0.5 and (a) 𝜆=0.5 and Pr=1, (b) 𝑅=2 and 𝜆=0.5, and (c) 𝑅=2, Pr=1.

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.

(a)

(b)

Figure 6 

Temperature profile for variables Ec andEc𝑥at 𝑅=2, 𝐾𝑟=0.5, 𝑀=1, Pr=1, and 𝜆=0.5 and (a) Ecx=0.5, (b) Ec=0.5.

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)

(b)

(c)

(d)

Figure 7 

The effect of active parameters on skin friction atEc=Ec𝑥=0.5 and (a) 𝑅=0.5, 𝑀=0.5, and 𝜆=1, (b) 𝑅=0.5, 𝐾𝑟=0.5, and 𝑀=0.5, (c) 𝐾𝑟=0.5, 𝑀=0.5, and 𝜆=1, (d) 𝑅=0.5, 𝐾𝑟=0.5, and 𝜆=1.

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.

(a)

(b)

(c)

(d)

Figure 8 

The effect of active parameters on Nusselt number (𝑁𝑢=𝜃(0)). (a) 𝐾𝑟=0.5, Pr=0.5, 𝜆=1.5 (b) 𝑅=0.5, Pr=0.5, 𝜆=1.5 (c) 𝑅=0.5, 𝐾𝑟=0.5, 𝜆=1.5 (d) 𝑅=0.5, 𝐾𝑟=0.5, Pr=0.5and 𝑀=0.5, Ec=Ec𝑥=0.5.

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.