Investigation of Nanofluid MHD Flow and Heat Transfer in a Channel

Similar documents
ABSTRACT. Keywords: Laminar viscous flow; magnetic field; optimal homotopy asymptotic method; porous channel ABSTRAK

Analysis of Natural Convection Flow in a Trapezoidal Cavity Containing a Rectangular Heated Body in Presence of External Oriented Magnetic Field

Parash Moni Thakur. Gopal Ch. Hazarika

USING MULTI-WALL CARBON NANOTUBE (MWCNT) BASED NANOFLUID IN THE HEAT PIPE TO GET BETTER THERMAL PERFORMANCE *

Application of Homotopy Perturbation Method (HPM) for Nonlinear Heat Conduction Equation in Cylindrical Coordinates

Application of Reconstruction of Variational Iteration Method on the Laminar Flow in a Porous Cylinder with Regressing Walls

OPTIMIZATION OF THE LONGITUDINAL FINS WITH DIFFERENT GEOMETRIES FOR INCREASING THE HEAT TRANSFER

THERMAL RADIATION EFFECTS ON MAGNETOHYDRODYNAMIC FLOW AND HEAT TRANSFER IN A CHANNEL WITH POROUS WALLS OF DIFFERENT PERMEABILITY

Numerical Analysis of MHD Flow of Fluid with One Porous Bounding Wall

G. C. Hazarika 2 Department of Mathematics Dibrugarh University, Dibrugarh

Homotopy Perturbation Method for Solving MHD Nanofluid Flow with Heat and Mass Transfer Considering Chemical Reaction Effect

NUMERICAL SOLUTION OF MHD FLOW OVER A MOVING VERTICAL POROUS PLATE WITH HEAT AND MASS TRANSFER

MIXED CONVECTION OF NEWTONIAN FLUID BETWEEN VERTICAL PARALLEL PLATES CHANNEL WITH MHD EFFECT AND VARIATION IN BRINKMAN NUMBER

FALLING FILM FLOW ALONG VERTICAL PLATE WITH TEMPERATURE DEPENDENT PROPERTIES

THE UNSTEADY FREE CONVECTION FLOW OF ROTATING MHD SECOND GRADE FLUID IN POROUS MEDIUM WITH EFFECT OF RAMPED WALL TEMPERATURE

ISSN Article

CONVECTIVE HEAT AND MASS TRANSFER IN A NON-NEWTONIAN FLOW FORMATION IN COUETTE MOTION IN MAGNETOHYDRODYNAMICS WITH TIME-VARING SUCTION

Application of He s homotopy perturbation method to boundary layer flow and convection heat transfer over a flat plate

Thermal Analysis of Convective-Radiative Circular Porous Fins with Sinusoidal Cross-Section Geometries using Galerkin Method

International Journal of Pure and Applied Mathematics

*Corresponding Author: Surajit Dutta, Department of Mathematics, C N B College, Bokakhat, Golaghat, Assam, India

Laplace Technique on Magnetohydrodynamic Radiating and Chemically Reacting Fluid over an Infinite Vertical Surface

Natural Convection in an Enclosure: Effect of Magnetic Field Dependent Thermal Conductivity

Chemical reaction Soret and Dufour Effect on Micropolar Fluid

Porous Plate In The Presence of Magnetic Field

Journal of Engineering Science and Technology Review 2 (1) (2009) Research Article

Riyadh 11451, Saudi Arabia. ( a b,c Abstract

Heat and Mass Transfer

Table of Contents. Foreword... xiii. Preface... xv

International Journal of Thermal Sciences

The Effect Of MHD On Laminar Mixed Convection Of Newtonian Fluid Between Vertical Parallel Plates Channel

UNSTEADY MAGNETOHYDRODYNAMICS THIN FILM FLOW OF A THIRD GRADE FLUID OVER AN OSCILLATING INCLINED BELT EMBEDDED IN A POROUS MEDIUM

Numerical Analysis of Laminar flow of Viscous Fluid Between Two Porous Bounding walls

Joule Heating Effect on the Coupling of Conduction with Magnetohydrodynamic Free Convection Flow from a Vertical Flat Plate

Boundary Layer Flow and Heat Transfer due to an Exponentially Shrinking Sheet with Variable Magnetic Field

Heat source/sink and thermal conductivity effects on micropolar nanofluid flow over a MHD radiative stretching surface

EffectofVariableThermalConductivityHeatSourceSinkNearaStagnationPointonaLinearlyStretchingSheetusingHPM

Problem 4.3. Problem 4.4

Finite Difference Solution of Unsteady Free Convection Heat and Mass Transfer Flow past a Vertical Plate

Effects of Viscous Dissipation on Unsteady Free Convection in a Fluid past a Vertical Plate Immersed in a Porous Medium

FORCED CONVECTION BOUNDARY LAYER MAGNETOHYDRODYNAMIC FLOW OF NANOFLUID OVER A PERMEABLE STRETCHING PLATE WITH VISCOUS DISSIPATION

Buoyancy Driven Heat Transfer of Water-Based CuO Nanofluids in a Tilted Enclosure with a Heat Conducting Solid Cylinder on its Center

Radiation Effect on MHD Casson Fluid Flow over a Power-Law Stretching Sheet with Chemical Reaction

Numerical solution of non-newtonian nanofluid flow over a stretching sheet

Analysis of the flow and heat transfer characteristics for MHD free convection in an enclosure with a heated obstacle

Research Article Innovation: International Journal of Applied Research; ISSN: (Volume-2, Issue-2) ISSN: (Volume-1, Issue-1)

On the Analysis of Laminar Flow of Viscous Fluid through a Porous Channel with Suction/Injection at Slowly Expanding or Contracting Walls

Hydromagnetic oscillatory flow through a porous medium bounded by two vertical porous plates with heat source and soret effect

Introduction. Statement of Problem. The governing equations for porous materials with Darcy s law can be written in dimensionless form as:

MHD Flow and Heat Transfer over an. Exponentially Stretching Sheet with Viscous. Dissipation and Radiation Effects

SOLUTION TO BERMAN S MODEL OF VISCOUS FLOW IN POROUS CHANNEL BY OPTIMAL HOMOTOPY ASYMPTOTIC METHOD

MHD Mixed Convection in Double Lid- Driven Differentially Heated Trapezoidal Cavity

COMBINED EFFECTS OF RADIATION AND JOULE HEATING WITH VISCOUS DISSIPATION ON MAGNETOHYDRODYNAMIC FREE CONVECTION FLOW AROUND A SPHERE

Numerical Solution of Mass Transfer Effects on Unsteady Flow Past an Accelerated Vertical Porous Plate with Suction

Thermophysical characteristics of ZnO nanofluid in L-shape enclosure.

Computers and Mathematics with Applications

Analysis of a Double Pipe Heat Exchanger Performance by Use of Porous Baffles and Nanofluids

MHD Non-Newtonian Power Law Fluid Flow and Heat Transfer Past a Non-Linear Stretching Surface with Thermal Radiation and Viscous Dissipation

EFFECTS OF CIRCULAR CORNERS AND ASPECT-RATIO ON ENTROPY GENERATION DUE TO NATURAL CONVECTION OF NANOFLUID FLOWS IN RECTANGULAR CAVITIES

A new approach for local similarity solutions of an unsteady hydromagnetic free convective heat transfer flow along a permeable flat surface

Effect of radiation on MHD nanofluid flow considering effective thermal conductivity and viscosity due to Brownian motion

Studies on flow through and around a porous permeable sphere: II. Heat Transfer

Flow of a micropolar fluid in channel with heat and mass transfer

Kabita Nath Department of Mathematics Dibrugarh University Dibrugarh, Assam, India

Mixed convection boundary layers in the stagnation-point flow toward a stretching vertical sheet

MIXED CONVECTION SLIP FLOW WITH TEMPERATURE JUMP ALONG A MOVING PLATE IN PRESENCE OF FREE STREAM

MAGNETOHYDRODYNAMIC GO-WATER NANOFLUID FLOW AND HEAT TRANSFER BETWEEN TWO PARALLEL MOVING DISKS

Nonlinear Analysis: Modelling and Control, 2008, Vol. 13, No. 4,

Conceptual Study of the Effect of Radiation on Free Convective Flow of Mass and Heat Transfer over a Vertical Plate

Unsteady MHD Couette Flow with Heat Transfer in the Presence of Uniform Suction and Injection

Numerical Study of Steady MHD Plane Poiseuille Flow and Heat Transfer in an Inclined Channel

UNSTEADY FREE CONVECTION BOUNDARY-LAYER FLOW PAST AN IMPULSIVELY STARTED VERTICAL SURFACE WITH NEWTONIAN HEATING

An assessment of a semi analytical AG method for solving two-dimension nonlinear viscous flow

Flow and heat transfer over a longitudinal circular cylinder moving in parallel or reversely to a free stream

Principles of Convection

The Effect of Suction and Injection on the Unsteady Flow Between two Parallel Plates with Variable Properties

Steady MHD Natural Convection Flow with Variable Electrical Conductivity and Heat Generation along an Isothermal Vertical Plate

Natural Convection and Entropy Generation in a Porous Enclosure with Sinusoidal Temperature Variation on the Side Walls

FREE CONVECTION OF HEAT TRANSFER IN FLOW PAST A SEMI-INFINITE FLAT PLATE IN TRANSVERSE MAGNETIC FIELD WITH HEAT FLUX

Unsteady Laminar Free Convection from a Vertical Cone with Uniform Surface Heat Flux

THE EFFECT OF VARIATION OF BASE FLUID ON NATURAL CONVECTION IN CAVITY FILLED WITH NAOFLUID IN THE PRESENCE OF MAGNETIC FIELD

MHD Boundary Layer Stagnation Point Flow and Heat Generation/ Absorption of a Micropolar Fluid with Uniform Suction / Injection

Unsteady MHD Convective Heat and Mass Transfer of a Casson Fluid Past a Semi-infinite Vertical Permeable Moving Plate with Heat Source/Sink

INFLUENCE OF VARIABLE PERMEABILITY ON FREE CONVECTION OVER VERTICAL FLAT PLATE EMBEDDED IN A POROUS MEDIUM

MHD free convection heat and mass transfer flow over a vertical porous plate in a rotating system with hall current, heat source and suction

FINITE ELEMENT ANALYSIS OF MIXED CONVECTION HEAT TRANSFER ENHANCEMENT OF A HEATED SQUARE HOLLOW CYLINDER IN A LID-DRIVEN RECTANGULAR ENCLOSURE

MHD flow and heat transfer near the stagnation point of a micropolar fluid over a stretching surface with heat generation/absorption

CONVECTIVE HEAT TRANSFER

Numerical Study of Natural Convection in. an Inclined L-shaped Porous Enclosure

Unsteady MHD Free Convection Flow past an Accelerated Vertical Plate with Chemical Reaction and Ohmic Heating

SIMULATION OF NATURAL CONVECTION HEAT TRANSFER USING NANOFLUID IN A CONCENTRIC ANNULUS

Rotating Flow of Magnetite-Water Nanofluid over a Stretching Surface Inspired By Non-Linear Thermal Radiation and Mass Transfer

Combined Effect of Buoyancy Force and Navier Slip on Entropy Generation in a Vertical Porous Channel

Radiative Mhd Stagnation Point Flow Over A Chemical Reacting Porous Stretching Surface With Convective Thermal Boundary Condition


Oscillatory MHD Mixed Convection Boundary Layer Flow of Finite Dimension with Induced Pressure Gradient

Numerical study of entropy generation and melting heat transfer on MHD generalised non-newtonian fluid (GNF): Application to optimal energy

A Semi-Analytical Solution for a Porous Channel Flow of a Non-Newtonian Fluid

MAGNETOHYDRODYNAMIC FLOW OF NANOFLUID OVER PERMEABLE STRETCHING SHEET WITH CONVECTIVE BOUNDARY CONDITIONS

Investigation of blood flow as third order non-newtonian fluid inside a porous artery in the presence of a magnetic field by an analytical method

Transcription:

Copyright 2014 by American Scientific Publishers All rights reserved. Printed in the United States of America Journal of Advanced Physics Vol. 3, pp. 1 11, 2014 www.aspbs.com/jap Investigation of Nanofluid MHD Flow and Heat Transfer in a Channel M. Fakou 1,,D.D.Ganji 2, and M. Abbasi 1 1 Department of Mechanical Engineering, Sari Branch, Islamic Azad University, Sari, Iran 2 Department of Mechanical Engineering, Babol University of Technology, Babol, Iran In this study, laminar fluid flow and heat transfer in channel with permeable walls in the presence of a transverse magnetic field is investigated. Homotopy perturbation method HPM has been employed for computing approximate solutions of nonlinear differential equations governing the problem. We have tried to show reliability and performance of the present method compared with the numerical method Runge-Kutta fourth-rate to solve this problem. The influence of the four dimensionless numbers: the Hartmann number, Reynolds number, Prandtl number and Eckert number on the non-dimensional velocity and temperature profiles have been investigated. The results show Homotopy perturbation method is very close to numerically method. In general, increasing the Reynolds number is increases the nanofluid flow velocity in the channel and the maximum amount of temperature increase and by applied magnetic field Ha, velocity and the maximum amount of temperature in the channel is reduces. Increasing Prandtl number and Eckert, respectively, the maximum value of theta increases and decreases. KEYWORDS: Homotopy Perturbation Method HPM, Laminar Viscous Flow, Heat Transfer, Uniform Magnetic Field. 1. INTRODUCTION Flow problem in a porous tube or channel received much attention in recent years due to its various applications in medical engineering, for example, in the dialysis of blood in artificial kidney, 1 in the flow of blood in the capillaries, 2 in the flow in blood oxygenators, 3 as well as in many other engineering areas such as the design of filters, 4 in transpiration cooling boundary layer control 5 and gaseous diffusion. 6 In 1953, Berman 7 described an exact solution of the Navier-Stokes equation for steady two-dimensional laminar flow of a viscous, incompressible fluid in a channel with parallel, rigid, porous walls driven by uniform, steady suction or injection at the walls. This mass transfer is paramount in some industrial processes. More recently, Chandran and Sacheti 8 analyzed the effects of a magnetic field on the thermodynamic flow past a continuously moving porous plate. Slow viscous flow problem in a semi-porous channel in the presence of transverse magnetic field is investigated by Sheikholeslami et al. 9 They show that the method is asymptotically optimal Homotopy a powerful method for solving nonlinear differential equations, such as the problem. Soleimani et al. 10 studied of natural convection heat Author to whom correspondence should be addressed. Emails: mehdi_fakour@yahoo.com, mehdi_fakoor8@yahoo.com Received: 13 October 2014 Accepted: 23 October 2014 transfer in an enclosure filled mid-loop with nanofluid using the Control Volume based Finite Element Method. They founded turn angle has a significant effect on flow lines, isotherms and local Nusselt number is maximum or minimum values. Sheikholeslami et al. 11 investigated the flow of nanofluid and heat transfer characteristics between two horizontal plates in a rotating system. Their results showed that for suction and injection, the heat transfer rate increases with the nanoparticle volume fraction, Reynolds number, and parameter injection/suction increases and then decreases with the strength of the spin parameter. Steady magneto hydrodynamic free convection boundary layer flow past a vertical semi-infinite flat plate embedded in water filled with a nanofluid has been theoretically studied by Hamad et al. 12 They found that copper and silver nanoparticles have proved that the highest cooling performance for this problem. Natural convection of a non- Newtonian copper-water nanofluid between two infinite parallel vertical flat plates was investigated by Domairry et al. 13 They concluded that as the size of nanoparticles increases, the boundary layer thickness increases, the thermal boundary layer thickness decreases. Sheikholeslami et al. 14 studied the natural convection in a concentric annulus between a cold outer square and heated inner circular cylinders in the presence of static radial magnetic field. Sheikholeslami et al. 15 performed a numerical analysis for natural convection heat transfer of Cu-water nanofluid in a cold outer circular enclosure containing J. Adv. Phys. 2014, Vol. 3, No. 4 2168-1996/2014/3/001/011 doi:10.1166/jap.2014.1153 1

Investigation of Nanofluid MHD Flow and Heat Transfer in a Channel Fakou et al. a hot inner sinusoidal circular cylinder in the presence of horizontal magnetic field using the Control Volume based Finite Element Method. They concluded that in the absence of a magnetic field, increasing the Rayleigh number increases, while the opposite trend was observed in the presence of a magnetic field decreases. Sheikholeslami et al. 16 studied the effects of magnetic field and nanoparticle on the Jeffery Hamel flow using Adomian decomposition method. They showed that increasing Hartmann number will lead to back flow reduction. Recently several authors investigated about nanofluid flow and heat 17 25 32 transfer. The main goal of this paper is to examine the laminar nanofluid flow in channel with permeable walls in the presence of transverse magnetic field using Homotopy perturbation method. Effective volume fraction nanofluid, Hartmann number, Reynolds number, Prandtl number, and Eckert number on the velocity and temperature considered. In addition to speed and temperature for different structures nanofluid copper and silver nanoparticles in water or ethylene glycol depicted. Fig. 2. Schematic of the problem nanofluid in a porous Channel and magnetic field. u + x y = 1 nf c p u T T + x y P nf y + nf 2 x + 2 2 y 2 2 T = k nf x + 2 T 2 y 2 u 2 + nf 4 y The suitable boundary conditions for the velocity are: 3 2. PROBLEM DESCRIPTION Steady two-dimensional laminar flow of an incompressible viscous electrically conducting fluid in a channel with permeable walls with a long rectangular plate with uniform translation in x,l x over an infinite porous plate and made to consider. The distance between the two plates is h. We are observing a normal velocity q on the porous walls. A uniform magnetic field B is assumed to be applied towards direction y Figs. 1 2. In the case of a short circuit to neglect the electrical field and perturbations to the basic normal field and without any gravity forces, the governing equations are: 9 u x + v u u u + x y = 1 P nf x + nf y = 0 1 2 u x + 2 u 2 y 2 u nf B 2 nf 2 y = 0 u = u 0 = q T = T 0 5 y = h u = u 0 = q T = T h 6 Calculating a mean velocity U by the relation: U h = h 0 u dy = L x q 7 We consider the following transformations: x = x L x y = y h u u= U = q P y = P q 2 = T T 0 T h T 0 Then, we can consider four dimensionless numbers: the Hartmann number Ha for the description of magnetic forces 11 and the Reynolds number Re for dynamic forces, Prandtl number and Eckert number: f Ha = Bh 9 f f 8 Re = hq nf 10 Table I. Different models for simulation of Dynamic viscosity. Model Electrical conductivity Dynamic viscosity I nf 3 s / f 1 = 1 + f s / f + 2 s / f 1 nf = f 1 2 5 Fig. 1. Three dimensional Schematic of the problem. II nf 3 s / f 1 = 1 + f s / f + 2 s f 1 nf = f 1 + 7 3 + 123 2 2 J. Adv. Phys., 3, 1 11, 2014

Fakou et al. Investigation of Nanofluid MHD Flow and Heat Transfer in a Channel Table II. Thermo physical properties of nanofluids and nanoparticles. Density Electrical Specific heat Thermal Material conductivity capacity C P conductivity K Silver 10500 6.30 10 7 230 418 Copper 8933 5.96 10 7 385 401 Ethylene 1113.2 1.07 10 4 2410 0.252 glycol Drinking 997.1 0.05 4179 0.613 water Ec = Pr = nf C P k nf 11 U 2 C p T h T 0 12 where the effective density nf and specific heat capacity C P are defined as: 12 nf = f 1 + s C P nf = C p f 1 + C p s 13 The effective thermal conductivity of the nanofluid can be approximated by the Maxwell Garnett s MG model as: 13 k nf = k s + 2k f 2 k f k s 14 k f k s + 2k f + k f k s where is the nanoparticle volume fraction. Different models for simulating dynamic viscosity of the nanofluids are shown in Table I. In the first model, the effective thermal conductivity and viscosity of nanofluid are calculated by the Maxwell 13 and Brinkman models, Table I. The thermo physical properties of the nanofluid are given in Table II. Introducing Eqs. 7 12 into Eqs. 1 4 leads to the dimensionless equations: u x + y = 0 15 u u u + P x y = 2 y x + nf hp u Ha2 Re 2 2 u x + 2 u 2 y 2 B A 16 Fig. 3. Comparison of HPM and numerical results for dimensionless velocities and temperature a U y,bv y and c y. J. Adv. Phys., 3, 1 11, 2014 3

Investigation of Nanofluid MHD Flow and Heat Transfer in a Channel Fakou et al. u + x C D Pr y = p y y + nf hq [ Re u + x y 2 = 2 x + 2 2 y 2 2 2 x + 2 2 y 2 u 2 ] y Ec where A, B, C and D are constant parameters: A = 1 + s f 3 B s / f 1 = 1 + s / f + 2 s / f 1 C = k s + 2k f 2 k f k s k s + 2k f + k f k s D = 1 + C P s C P f 17 18 19 20 Quantity of is defined as the aspect ratio between distance h and a characteristic length L x of the slider. This ratio is normally small. Berman s similarity transformation is used to be free from the aspect ratio of : = V y u= u U = u 0 U y + x dv dy 21 Introducing Eq. 21 in the second momentum Eqs. 19, 20 shows that quantity P y / y does not depend on the longitudinal variable x. With the first momentum equation, we also observe that 2 P y / x 2 is independent of x. Table IV. Comparison between y results from applied method for Re = Ha = Pr = Ec = 1and = 0.04. Y HPM NUM Error 0.00 0 00000000 0 00000000 0 000000000 0.05 0 03926615 0 044244071 0 004617456 0.0 0 081805159 0 090339094 0 008533934 0.15 0 126616048 0 138245080 0 011629032 0.20 0 174041620 0 187923007 0 013881387 0.25 0 223960615 0 239318270 0 015357655 0.30 0 276148632 0 292338734 0 016190102 0.35 0 330284539 0 346832006 0 016547467 0.40 0 385962387 0 402566078 0 016603691 0.45 0 442708104 0 459216991 0 016508887 0.50 0 499999999 0 516366327 0 016366328 0.55 0 557291894 0 573510120 0 016218226 0.60 0 614037612 0 630079324 0 016041712 0.65 0 669715459 0 685470318 0 015754859 0.70 0 723851367 0 739082677 0 015231310 0.75 0 776039384 0 790360529 0 014321145 0.80 0 825958379 0 838833283 0 012874904 0.85 0 873383951 0 884151530 0 010767579 0.90 0 918194840 0 926114515 0 007919675 0.95 0 960373383 0 964686543 0 00431316 1.00 1 000000000 1 000000000 0 00000000 We omit asterisks for simplicity. Then a separation of variables leads to: 11 UV VU = 1 1 Re A 1 2 5 U Ha 2 B 1 2 5 U 22 V IV = Ha 2 B 1 2 5 V + ReA 1 2 5 V V VV 23 Table III. Comparison between V y and U y results from applied method for Re = Ha = 1and = 0.04. Y HPM U HPM V NUM U NUM V Error U Error 0.00 1.000000000 1 00000000 1.000000000 1 000000000 0 000000 0.000000 0.05 0.970000437 0 985502286 0.970010437 0 9855016891 0 000010 0.000006 0.10 0.938748271 0 943969728 0.938750271 0 9439678628 0 000002 0.000002 0.15 0.907737512 0 878387967 0.907734512 0 8783847961 0 000003 0.000003 0.20 0.878269638 0 791779213 0.878265638 0 7917750853 0 000004 0.000004 0.25 0.851462274 0 687179379 0.851458274 0 6871748223 0 000004 0.000004 0.30 0.828255114 0 567622050 0.828252114 0 5676176222 0 000003 0.000004 0.35 0.809413786 0 436128623 0.809409786 0 4361248265 0 000004 0.000004 0.40 0.795532216 0 295703859 0.795527216 0 2957010956 0 000005 0.000003 0.45 0.787033986 0 149335935 0.787030986 0 1493344855 0 000002 0.000002 0.50 0.784173037 0 0000000000 0.784171037 0 0000000000 0 000001 0.000000 0.55 0.787033986 0 1493359392 0.787028986 0 1493344855 0 000005 0.000002 0.60 0.795532217 0 2957038599 0.795531217 0 2957010956 0 000002 0.000003 0.65 0.809413786 0 4361286244 0.809411786 0 4361248265 0 000001 0.000004 0.70 0.828255114 0 5676220507 0.828251114 0 5676176222 0 000004 0.000004 0.75 0.851462273 0 6871793806 0.851460273 0 6871748223 0 000002 0.000004 0.80 0.878269638 0 7917792140 0.878265638 0 7917750853 0 000004 0.000004 0.85 0.907737512 0 8783879694 0.907734512 0 8783847961 0 000003 0.000003 0.90 0.938748271 0 9439697288 0.938746271 0 9439678628 0 000002 0.000002 0.95 0.970000437 0 9855022883 0.970010437 0 9855016891 0 000001 0.000006 1.00 1.000000000 1 0000000000 1.000000000 1 0000000000 0 000000 0.000000 4 J. Adv. Phys., 3, 1 11, 2014

Fakou et al. Investigation of Nanofluid MHD Flow and Heat Transfer in a Channel = C D 1 2 5 Pr Re 1 2 5 V +Ec u 0 U 2 24 Where primes denote differentiation with respect to y and asterisks have been omitted for simplicity. The dynamic boundary conditions are: y = 0 U= 1 V = 1 V = 0 = 0 y = 1 U= 1 V = 1 V = 0 = 1 25 where A u is defined as follows: A u = L u + N u 28 Homotopy-perturbation structure is shown as: H p = L L u 0 + p L u 0 + p N f r = 0 29 or 3. DESCRIBE HOMOTOPY PERTURBATION METHOD AND APPLIED TO THE PROBLEM 3.1. Describe Homotopy Perturbation Method Consider the following function: A u f r = 0 26 with the boundary condition of: B = u u = 0 27 n H p = 1 p L L u 0 + p A f r = 0 30 where r p 0 1 R 31 Obviously, considering Eqs. 29 and 30, we have: H 0 = L L u 0 = 0 32 H 1 = A f r = 0 Where p 0 1 is an embedding parameter and u 0 is the first approximation that satisfies the boundary condition. Fig. 4. Effect of nanoparticle volume fraction,, onav y,bu y and c y, for water with copper nanoparticles when Re = 1andHa= 1. J. Adv. Phys., 3, 1 11, 2014 5

Investigation of Nanofluid MHD Flow and Heat Transfer in a Channel Fakou et al. The process of the changes in p from zero to unity is that of v r p changing from u 0 to u r. We consider v as: = 0 + p 1 + p 2 2 + 33 and the best approximation is: u = lim p 1 = 0 + 1 + 2 + 34 the above convergence is discussed in Refs. [28 30]. 3.2. The HPM Applied to the Problem A homotopy perturbation method can be constructed as follows: H U p = 1 p 1 Re 1 U A 1 2 5 + p UV VU 1 Re A 1 2 5 U Ha 2 B 1 2 5 U 1 35 H V p = 1 p V IV + p V IV Ha 2 B 1 2 5 V + ReA 1 2 5 V V VV 36 H p = 1 p + p + C 1 2 5 D Pr Re 1 2 5 V + Ec u 0 U 2 37 One can now try to obtain a solution of Eqs. 35 37 in the form of: U y = U 0 y + pu 1 y + p 2 U 2 y + 38 V y = V 0 y + pv 1 y + p 2 V 2 y + 39 y = 0 y + p 1 y + p 2 2 y + 40 where v i Y i = 1 2 3 are functions yet to be determined. According to Eqs. 35 37 the initial approximation to satisfy boundary condition is: U 0 Y = 1 41 V 0 Y = 4Y 3 + 6Y 2 1 42 0 Y = Y 43 Fig. 5. Effect of Hartmann number Ha on dimensionless velocities and temperature for water with copper nanoparticles, = 0.04, a V y, Re = 1, b U y,re= 1, c y, Re= 1, d V y,re= 5, e U y,re= 5, f y, Re= 5. 6 J. Adv. Phys., 3, 1 11, 2014

Fakou et al. Investigation of Nanofluid MHD Flow and Heat Transfer in a Channel Substituting Eqs. 38 40 into Eqs. 35 37 yields: 0 841 0 828 d2 dy U 1 Y 12Y 2 + 12Y = 0 44 2 28 98 4Y 3 + 6Y 2 1 1 015 24Y + 12 1 207 12Y 2 + 12Y 24Y + 12 + d4 dy 4 V 1 Y = 0 45 1 226 4y 3 + 6y 2 1 + d2 = 0 46 d 1 2 The solutions of Eqs. 44 46 may be written as follow: U 1 Y = 0 507Y 2 1 207Y 4 + 2 415Y 3 1 707Y 47 V 1 Y = 0 276Y 7 0 966Y 6 1 023Y 5 1 207Y 4 + 0 507Y 4 + 0 206Y 3 + 0 002Y 2 48 1 Y = 0 661Y 2 + 0 1Y 5 0 25Y 4 0 245Y 49 In the same manner, the rest of components were obtained by using the Maple package. According to the HPM, we can conclude: U Y = 1 0 527Y 2 1 427Y 4 + 2 328Y 3 0 596Y 0 062Y 8 + 0 249Y 7 0 455Y 6 + 0 490Y 5 50 V Y = 0 655Y 4 + 0 030Y 8 + 0 067Y 9 + 1 083Y 5 +5 995Y 2 0 703Y 6 + 0 071Y 7 1 0 044Y 10 +0 008Y 11 3 852Y 3 51 Y = 0 424Y 5 + 0 764Y + 0 384Y 7 + 0 045Y 9 +0 285Y 6 0 349Y 8 0 471Y 4 +0 254Y 3 + 0 463Y 2 52 and so on. In the same manner the rest of the components of the iteration formula can be obtained. 4. RESULTS AND DISCUSSION In the present study HPM method is applied to obtain an explicit analytic solution of the laminar nanofluid flow and Fig. 6. Effect of Reynolds number Re on dimensionless velocities and temperature for water with copper nanoparticles, = 0.04, a V y, Ha = 1, b U y,ha=1, c y, Ha= 1, d V y,ha= 10, e U y,ha= 10, f y, Ha= 10. J. Adv. Phys., 3, 1 11, 2014 7

Investigation of Nanofluid MHD Flow and Heat Transfer in a Channel Fakou et al. Fig. 7. Effect of Prandtl and Eckert number on dimensionless temperature for water with copper nanoparticles, = 0.04, a y, Ha= Re = Ec = 1, b y, Ha= Re = Ec = 10, c y, Ha= Re = Pr = 1. Table V. U y variations for the different types of nanofluids and nanoparticles, Re = Ha = 1and = 0 04. Y Water-copper Water-silver Ethylene glycol-copper Ethylene glycol-silver 0.00 1.0000000000 1.0000000000 1.00000000000 1.00000000000 0.05 0.9498232203 0.94904876890 0.95030075741 0.94958595388 0.10 0.9007322187 0.89906408666 0.90175676951 0.0022196494 0.15 0.85411387721 0.85149514560 0.85571750591 0.85331377222 0.20 0.81119606695 0.80762715725 0.81337655710 0.81010659565 0.25 0.77304840015 0.76858150005 0.77577269854 0.77168570082 0.30 0.74058140741 0.73531420441 0.74378943819 0.73897536682 0.35 0.71454486429 0.70861357206 0.71815372724 0.71273700140 0.40 0.69552553611 0.68909722667 0.69943408574 0.69356668020 0.45 0.68394474261 0.67720904675 0.68803851565 0.68189252520 0.50 0.68005622911 0.67321653057 0.68421265137 0.67797242697 0.55 0.68394474260 0.67720904674 0.68803851567 0.68189252521 0.60 0.69552553609 0.68909722665 0.69943408579 0.69356668022 0.65 0.71454486427 0.70861357203 0.71815372730 0.71273700143 0.70 0.74058140739 0.73531420438 0.74378943826 0.73897536685 0.75 0.77304840012 0.76858150001 0.77577269862 0.77168570086 0.80 0.81119606692 0.80762715721 0.81337655719 0.81010659569 0.85 0.85411387719 0.85149514557 0.85571750600 0.85331377226 0.90 0.90073221874 0.89906408664 0.90175676957 0.90022196497 0.95 0.94982322029 0.94904876889 0.95030075745 0.94958595390 1.00 1.0000000000 1.0000000000 1.00000000000 1.00000000000 8 J. Adv. Phys., 3, 1 11, 2014

Fakou et al. Investigation of Nanofluid MHD Flow and Heat Transfer in a Channel Table VI. V y variations for the different types of nanofluids and nanoparticles, Re = Ha = 1and = 0 04. Y Water-copper Water-silver Ethylene glycol-copper Ethylene glycol-silver 0.00 1 0000000000 1 0000000000 1 0000000000 1 0000000000 0.05 0 9854928636 0 9855016891 0 9854875921 0 9854955316 0.10 0 9439412149 0 9439678628 0 9439253032 0 9439492696 0.15 0 8783404086 0 8783847961 0 8783139116 0 8783538238 0.20 0 7917180520 0 7917750853 0 7916840134 0 7917352875 0.25 0 6871123986 0 6871748223 0 6870751489 0 6871312617 0.30 0 5675573543 0 5676176222 0 5675213954 0 5675755650 0.35 0 4360734267 0 4361248266 0 4360427614 0 4360889573 0.40 0 2956638427 0 2957010956 0 2956416186 0 2956750985 0.45 0 1493149602 0 1493344855 0 1493033123 0 1493208597 0.50 0 0000000000 0 0000000000 0 00000000000 0 00000000000 0.55 0 14931496023 0 14933448559 0 14930331238 0 14932085970 0.60 0 29566384271 0 29570109567 0 29564161866 0 29567509852 0.65 0 43607342672 0 43612482655 0 43604276149 0 43608895731 0.70 0 56755735436 0 56761762221 0 56752139548 0 56757556509 0.75 0 68711239862 0 68717482237 0 68707514897 0 68713126177 0.80 0 79171805207 0 79177508530 0 79168401349 0 79173528756 0.85 0 87834040866 0 87838479611 0 87831391166 0 87835382383 0.90 0 94394121497 0 94396786284 0 94392530329 0 94394926968 0.95 0 98549286362 0 98550168912 0 98548759210 0 98549553164 1.00 1 0000000000 1 0000000000 1 00000000000 1 00000000000 heat transfer in a channel with permeable walls in the presence of uniform magnetic field Fig 1. First, a comparison between the applied method, HPM and Numerical method is investigated. For this aim Eqs. 22 24 are solved for different nanofluid structures see Table II and comparison between the described methods is demonstrated in Figure 3. Other comparisons are presented in Tables III and IV. It can be concluded from table values that HPM is very close to numerically method. Figure 4 shows the effect of nanoparticle volume fraction on U y, V y and y for water with copper nanoparticles when Re = Ha = Pr = Ec = 1. Nanoparticle volume fraction has no significant effect on velocity and temperature profiles. Effect of Hartmann number Ha on dimensionless velocities and temperature for water with copper is shown in Figure 5. Applied magnetic field, the channels will not have much impact on the V y and y profiles. But generally increasing Hartmann number decreases rapidly, especially in the middle of the channel and this decreases rapidly is Table VII. y variations for the different types of nanofluids and nanoparticles, Re = Ha = Pr = Ec = 1and = 0.04. Y Water-copper Water-silver Ethylene glycol-copper Ethylene glycol-silver 0.00 0 00000000000 0 00000000000 0.000000000000 0.00000000000 0.05 0 05417352416 0 044244071846 0.045212682600 0.04467504873 0.10 0 10791227135 0 090339094430 0.092021986721 0.09108912824 0.15 0 16140594703 0 13824508022 0.140402578954 0.13920828022 0.20 0 21490232217 0 18792300769 0.190333492174 0.18900126077 0.25 0 26867040127 0 23931827096 0.241780641748 0.24042250035 0.30 0 32295613003 0 29233873463 0.294675146740 0.29339020952 0.35 0 37793572995 0 34683200620 0.348891448808 0.34776396742 0.40 0 43367207592 0 40256607872 0.404228821405 0.40332569755 0.45 0 49007966534 0 45921699182 0.460399412729 0.45976745571 0.50 0 54690317049 0 51636632725 0.517025246847 0.51668867320 0.55 0 60371317931 0 57351012053 0.573645558505 0.57360434716 0.60 0 65992073776 0 63007932481 0.629734614965 0.62996432379 0.65 0 71480991850 0 68547031849 0.684728783926 0.68518228727 0.70 0 76758533037 0 73908267738 0.738060515464 0.73867187539 0.75 0 81742979732 0 79036052912 0.789196112463 0.78988648748 0.80 0 86356627659 0 83883328324 0.837673673345 0.83835883932 0.85 0 90531772608 0 88415153093 0.883137534341 0.88373629275 0.90 0 94215924433 0 92611451569 0.925365998457 0.92580852790 0.95 0 97375826447 0 96468654335 0.964289923574 0.96452501175 1.00 1 0000000000 1 00000000000 1.000000000000 1.00000000000 J. Adv. Phys., 3, 1 11, 2014 9

Investigation of Nanofluid MHD Flow and Heat Transfer in a Channel Fakou et al. more sensible in high Reynolds number. The maximum value of U y in the middle of the channel when the Hartmann and Reynolds number is small, it happens. Effect of Reynolds number Re on dimensionless velocities and temperature is shown in Figure 6. When the Hartmann number is small, the effect of Reynolds number on the velocities and temperature profiles are more sensible. Generally decreases rapidly with increasing Reynolds number and the maximum value of U y in the middle of the channel when the Hartmann and Reynolds number is very small, it happens. when Re = 0 the temperature profile is linear and with increasing Reynolds number the temperature profile is parabolic and with increasing Reynolds number, decreases temperature and heat transfer will be done later. Effect of Prandtl and Eckert number on dimensionless temperature is shown in Figure 7. When Hartmann and Reynolds number is small, the effect of Prandtl number on the temperature profile is more sensible. When the Pr = 0, the temperature profile is linear and with the Prandtl number increases profile is parabolic and temperature decreases. Eckert number variation has little impact on the temperature profile. The effects of the nanoparticle and liquid phase material on velocities and temperature profiles are shown in Tables V VII for U y, V y and y, respectively. These tables reveal that when nanofluid includes copper as nanoparticles or ethylene glycol as fluid phase in its structure, the U y, V y and y values are greater than the other structures. 5. CONCLUSION In this study, Homotopy Perturbation Method is applied to solve the problem of laminar nanofluid flow and heat transfer in a channel with permeable walls in the presence of uniform magnetic field. First a comparison between the applied method, HPM and Numerical method is investigated. The results show Homotopy perturbation method is very close to numerically method. In general, increasing the Reynolds number is increases the nanofluid flow velocity in the channel and the maximum amount of temperature increase and by applied magnetic field Ha, velocity and the maximum amount of temperature in the channel is reduces. Increasing Prandtl number and Eckert, respectively, the maximum value of theta increases and decreases. Also nanoparticle volume fraction has no significant effect on velocity and temperature profiles but the effect in the velocity profile is minimal when volume fraction of nanoparticle is decreased. NOMENCLATURE HPM Homotopy perturbation method NUM Numerical method P Pressure q Mass transfer parameter Re Reynolds number U Dimensionless velocity in the x direction V Dimensionless velocity in the y direction h Suspension height Ec Eckert number Ha Hartmann number Pr Prandtl number Length of the slider L x C p Specific heat k Thermal conductivity T = T h T 0 Difference temperature between the plates Fluid thermal diffusivity u 0 x velocity of the pad u Dimensionless x-component velocity v Dimensionless y-component velocity Velocity component in the x direction Velocity component in the y direction x Dimensionless horizontal coordinate y Dimensionless vertical coordinate x Distance in the x direction parallel to the plates y Distance in the y direction parallel to the plates. u v Greek Symbols Fluid density Dimensionless temperature Kinematic viscosity Electrical conductivity Aspect ratio h/lx. References and Notes 1. V. Wernert, O. Schäf, H. Ghobarkar, and R. Denoyel, Micro Porous Mesoporous Mater 83, 101 2005. 2. A. Jafari, P. Zamankhan, S. M. Mousavi, and P. Kolari, Commun. Nonlinear Sci. Numer. Simul. 14, 1396 2009. 3. A. R. Goerke, J. Leung, and S. R. Wickramasinghe, Chem. Eng. Sci. 57, 2035 2002. 4. S. S. Mneina and G. O. Martens, AEU Int. J. Electron. Commun 63, 83 2009. 5. Y. H. Andoh and B. Lips, Appl. Therm. Eng. 23, 1947 2003. 6. A. Runstedtler, Chem.Eng.Sci.61, 5021 2006. 7. A. S. Berman, J. Appl. Phys. 24, 1232 1953. 8. A. K. S. P. Chandran and N. C. Sacheti, Int. Commun. Heat Mass Transf. 23, 889 1996. 9. M. Sheikholeslami, D. Domairry, H. R. Ashorynejad, and I. Hashim, Sains Malaysiana 10, 1177 2012. 10. S. Soleimani, M. Sheikholeslami, D. D. Ganji, and M. Gorji-Bandpy, Int. Commun. Heat Mass Transf. 39, 565 2012. 11. M. Sheikholeslami, H. R. Ashorynejad, G. Domairry, and I. Hashim, Hindawi Publ. Corp. J. Appl. Math. 2012, 19 2012. 12. M. A. A. Hamad, I. Pop, and A. I. Md Ismail, Nonlinear Anal. Real World Appl. 12, 1338 2011. 13. D. Domairry, H. R. Sheikholeslami, M. Ashorynejad, R. Subba, R. Gorla, and M. Khani, Nanoeng. Nanosyst. 225, 115 2012. 14. M. Sheikholeslami, M. Gorji-Bandpy, and D. D. Ganji, Int. Commun. Heat Mass Transf. 39, 978 2012. 10 J. Adv. Phys., 3, 1 11, 2014

Fakou et al. Investigation of Nanofluid MHD Flow and Heat Transfer in a Channel 15. M. Sheikholeslami, M. Gorji-Bandpy, D. D. Ganji, S. Soleimani, and S. M. Seyyedi, Int. Commun. Heat Mass Transf. 39, 1435 2012. 16. M. Sheikholeslami, D. D. Ganji, H. R. Ashorynejad, and H. B. Rokni, Appl. Math. Mech. English Ed. 1, 25 2012. 17. M. Hatami and D. D. Ganji, J. Mol. Liq. 188, 155 2013. 18. M. Hatami, R. Nouri, and D. D. Ganji, J. Mol. Liq. 187, 294 2013. 19. M. Hatami and D. D. Ganji, J. Mol. Liq. 190, 159 2014. 20. G. Domairry and M. Hatami, J. Mol. Liq. 193, 37 2014. 21. M. Hatami and D. D. Ganji, Case Studies in Thermal Engineering 2, 14 2014. 22. M. Hatami and D. D. Ganji, Energy Convers. Manage. 78, 347 2014. 23. M. Hatami, M. Sheikholeslami, and D. D. Ganji, J. Mol. Liq. 195, 230 2014. 24. M. Hatami, M. Sheikholeslami, M. Hosseini, and D. D. Ganji, J. Mol. Liq. 194, 251 2014. 25. M. Sheikholeslami, M. Hatami, and D. D. Ganji, J. Mol. Liq. 194, 30 2014. 26. H. R. Ashorynejad, A. A. Mohamad, and M. Sheikholeslami, Int. J. Therm. Sci. 64, 240 2013. 27. A. Desseaux, Int. J. Eng. Sci. 37, 1781 1999. 28. J. C. Maxwell, A Treatise on Electricity and Magnetism, Second edn,, Oxford University Press, Cambridge 1904, pp. 435 441. 29. M. Fakour, A. Vahabzadeh, D. D. Ganji, and I. Rahimi Petroudi, Cent. Eur. J. Eng. 4, 341 2014. 30. D. D. Ganji and A. Sadighi, J. Comput. Appl. Math. 207, 24 2007. 31. M. Fakour, A. Vahabzadeh, and D. D. Ganji, International Journal of Case Studies in Thermal Engineering 2014, http://dx.doi.org/10.1016/j.csite.2014.05.003. 32. M. Sheikholeslami, M. Hatami, and D. D. Ganji, Powder Technol. 246, 327 2013. 33. M. Fakour, D. D. Ganji, A. Vahabzadeh, and S. H. H. Kachapi, Walailak Journal of Science and Technology 11, 593 2014. J. Adv. Phys., 3, 1 11, 2014 11