Research Article | Open Access

K. Ramesh, M. Devakar, "Effects of Heat and Mass Transfer on the Peristaltic Transport of MHD Couple Stress Fluid through Porous Medium in a Vertical Asymmetric Channel", *Journal of Fluids*, vol. 2015, Article ID 163832, 19 pages, 2015. https://doi.org/10.1155/2015/163832

# Effects of Heat and Mass Transfer on the Peristaltic Transport of MHD Couple Stress Fluid through Porous Medium in a Vertical Asymmetric Channel

**Academic Editor:**Jamshid M. Nouri

#### Abstract

The intrauterine fluid flow due to myometrial contractions is peristaltic type motion and the myometrial contractions may occur in both symmetric and asymmetric directions. The channel asymmetry is produced by choosing the peristaltic wave train on the walls to have different amplitude, and phase due to the variation of channel width, wave amplitudes and phase differences. In this paper, we study the effects of heat and mass transfer on the peristaltic transport of magnetohydrodynamic couple stress fluid through homogeneous porous medium in a vertical asymmetric channel. The flow is investigated in the wave frame of reference moving with constant velocity with the wave. The governing equations of couple stress fluid have been simplified under the long wave length approximation. The exact solutions of the resultant governing equations have been derived for the stream function, temperature, concentration, pressure gradient, and heat transfer coefficients. The pressure difference and frictional forces at both the walls are calculated using numerical integration. The influence of diverse flow parameters on the fluid velocity, pressure gradient, temperature, concentration, pressure difference, frictional forces, heat transfer coefficients, and trapping has been discussed. The graphical results are also discussed for four different wave shapes. It is noticed that increasing of couple stresses and heat generation parameter increases the size of the trapped bolus. The heat generation parameter increases the peristaltic pumping and temperature.

#### 1. Introduction

In recent years, the flow of non-Newtonian fluids has received much attention due to the increasing industrial, medical, and technological applications. Various researchers have attempted diverse flow problems related to several non-Newtonian fluids and couple stress fluid is one of them. The theory of couple stress fluids originated by Stokes [1] has many biomedical, industrial, and scientific applications and was used to model synthetic fluids, polymer thickened oils, liquid crystals, animal blood, and synovial fluid. Some earlier developments in couple stress fluid theory with some basic flows can be found in the book by Stokes [2]. Recently, few researchers have studied some couple stress fluid flows for different flow geometries [3â€“8].

Nowadays, peristaltic flows have gained much attention because of their applications in physiology and industry. Peristaltic transport is a form of fluid transport induced by a progressive wave of area contraction or expansion along the length of a distensible tube/channel and transporting the fluid in the direction of the wave propagation. This phenomenon is known as peristalsis. In physiology this plays an important role in various situations such as the food movement in the digestive tract, urine transport from kidney to bladder through ureter, movement of lymphatic fluids in lymphatic vessels, bile flow from the gall bladder into the duodenum, spermatozoa in the ductus efferentes of the male reproductive tract, ovum movement in the fallopian tube, blood circulation in the small blood vessels, the movement of the chyme in the gastrointestinal tract, intrauterine fluid motion, swallowing food bolus through esophagus, and transport of cilia. Many industrial and biological instruments such as roller pumps, finger pumps, heart-lung machines, blood pump machines, and dialysis machines are engineered based on the peristaltic mechanism [9]. The intrauterine fluid flow due to myometrial contractions is peristaltic in nature and these myometrial contractions occur in both symmetric and asymmetric directions and also when embryo enters the uterus for implantation there start the asymmetric contractions. The contractions inside the nonpregnant uterus are very complicated because they are composed of variable amplitudes and different wave lengths [10]. In view of this, Pandey and Chaube [11] have investigated the peristaltic transport of a couple stress fluid in a symmetrical channel using perturbation method in terms of small amplitude ratio. Ali and Hayat [12] have studied the peristaltic motion of micropolar fluid in an asymmetric channel. Naga Rani and Sarojamma [13] have analyzed the peristaltic transport of a Casson fluid in an asymmetric channel. Hayat et al. [14] have discussed the peristaltic flow of a Johnson-Segalman fluid in an asymmetric channel. Hayat and Javed [15] have studied the peristaltic transport of power-law fluid in asymmetric channel.

The porous medium plays an important role in the study of transport process in biofluid mechanics, industrial mechanics, and engineering fields. The fluid transport through porous medium is widely applicable in the vascular beds, lungs, kidneys, tumorous vessels, bile duct, gall bladder with stones, and small blood vessels. In the pathological situations, the distribution of fatty cholesterol, artery clogging, blood clots in the lumen of coronary artery, transport of drugs and nutrients to brain cells, and functions of organs are modeled as porous medium [16]. Recently, Tripathi [17] studied the peristaltic hemodynamic flow of couple stress fluids through a porous medium. Tripathi and BÃ©g [18] have investigated the peristaltic flow of generalized Maxwell fluid through a porous medium using homotopy perturbation method. Abd elmaboud and Mekheimer [19] have discussed peristaltic transport of a second-order fluid through a porous medium using regular perturbation method. The magnetohydrodynamic flows also gained much attention due to the widespread applications in biofluid mechanics and industry. It is the fact that many fluids like blood are conductive in nature and gave a new direction for research. The indispensable role of biomagnetic fluid dynamics in medical science has been very helpful with many problems of physiology. It has wide range of applications, such as magnetic wound or cancer tumour treatment, bleeding reduction during surgeries, provocation of occlusion of feeding vessels of cancer tumor, cell separation, transport of drugs, blood pump machines, and magnetic resonance imaging to diagnose the disease and the influence of magnetic field which may be utilized as a blood pump in carrying out cardiac operations for the blood flow in arteries with arterial disease like arterial stenosis or arteriosclerosis. Specifically, the magnetohydrodynamic flows of non-Newtonian fluids are of great interest in magnetotherapy. The noninvasive radiological tests use the magnetic field to evaluate organs in abdomen [20]. Hayat et al. [21] have studied the peristaltic transport of magnetohydrodynamic Johnson-Segalman fluid for the case of a planar channel. Wang et al. [22] have investigated the peristaltic motion of a magnetohydrodynamic generalized second-order fluid in an asymmetric channel. Nadeem and Akram [23, 24] have discussed the peristaltic transport of a couple stress fluid and Williamson fluid in an asymmetric channel with the effect of the magnetic field.

Heat transfer plays a significant role in the cooling processes of industrial and medical applications. Such consideration is very important since heat transfer in the human body is currently considered as an important area of research. In view of the thermotherapy and the human thermoregulation system, the model of bioheat transfer in tissues has been attracted by the biomedical engineers. In fact the heat transfer in human tissues involves complicated processes such as heat conduction in tissues, heat transfer due to perfusion of the arterial-venous blood through the pores of the tissue, metabolic heat generation, and external interactions such as electromagnetic radiation emitted from cell phones [25]. Heat transfer also involves many complicated processes such as evaluating skin burns, destruction of undesirable cancer tissues, dilution technique in examining blood flow, paper making, food processing, vasodilation, and radiation between surface and its environment [26]. Mustafa et al. [27] have studied the peristaltic transport of nanofluid in a channel. The heat transfer characteristics of a couple stress fluid in an asymmetric channel have been analyzed by Abd elmaboud et al. [28]. Nadeem and Akbar [29] have discussed the influence of heat transfer and magnetic field on peristaltic flow of a Johnson-Segalman fluid in a vertical symmetric channel. Some more works regarding peristaltic flows with the effect of heat transfer and magnetic field can be seen in [30â€“33]. Srinivas et al. [34] have studied the effects of both wall slip conditions and heat transfer on peristaltic flow of MHD Newtonian fluid in a porous channel with elastic wall properties. Mass transfer is another important phenomenon in physiology and industry. This phenomenon has great applications such as nutrientsâ€™ diffusion out from the blood to neighboring tissues, membrane separation process, reverse osmosis, distillation process, combustion process, and diffusion of chemical impurities [35]. Recently, Noreen [36] studied the problem of mixed convection peristaltic flow of third-order nanofluid with an induced magnetic field. Saleem and Haider [37] have discussed the peristaltic transport of Maxwell fluid with heat and mass transfer in an asymmetric channel. Some more relevant works on the peristaltic transport with heat and mass transfer can be seen in [38â€“42].

The aim of the present study is to investigate the influence of heat and mass transfer on the peristaltic flow of magnetohydrodynamic couple stress fluid through homogeneous porous medium in a vertical asymmetric channel. This paper is arranged as follows. Section 2 presents the mathematical formulation for the problem. The solution of the problem is obtained in Section 3. The four different wave forms are presented in Section 4 while the computational results are discussed in Section 5. The last section, Section 6, presents the conclusions of the present study.

#### 2. Formulation of the Problem

Let us consider magnetohydrodynamic couple stress fluid in a vertical asymmetric channel through the porous medium with the width of . The surfaces and of the asymmetric channel are maintained at constant temperatures and and the constant concentrations and , respectively (see Figure 1). The porous medium is assumed to be homogeneous. The motion is induced by sinusoidal wave trains propagating with constant speed along the channel walls as defined by the following:where and are the wave amplitudes, is the wave length, is the channel width, is the velocity of propagation, is the time, and is the direction of wave propagation. The phase difference varies in the range , in which corresponds to symmetric channel with waves out of phase and corresponds to waves in phase, and further , , , , and meet the following relation .

The continuity, momentum, energy, and concentration equations for an MHD incompressible couple stress fluid, in the absence of body couples, are [8, 16]in which represents the material derivative and is Darcy's resistance in the porous medium which are given bywhere is the velocity vector, is the density, is the pressure, is the viscosity, is material constant associated with couple stress, is the electric current density, is the total magnetic field, is the acceleration due to the gravity, is the coefficient of thermal expansion, is the coefficient of expansion with concentration, is the specific heat at constant pressure, is the temperature, is the mass concentration, is the thermal conductivity, is the constant heat generation parameter, is the coefficient of mass diffusivity, is the thermal diffusion ratio, is the mean temperature, and is the permeability parameter. The viscous dissipation is neglected in the energy equation.

In the fixed frame, governing equations for the peristaltic motion of an incompressible magnetohydrodynamic couple stress fluid through homogeneous porous medium in the two-dimensional vertical channel arein which and are the respective velocity components, is the pressure, is the temperature and is the concentration in the reference to fixed frame system, is the electrical conductivity of the fluid, and is the applied magnetic field. In the above, the induced magnetic field is neglected since the magnetic Reynolds number is assumed to be small.

The coordinates, velocities, pressure, temperature, and concentration in the fixed frame and wave frame are related by the following expressions:in which , , , , and are velocity components, pressure, temperature, and concentration in the wave frame, respectively.

Using (5), the governing equations in the wave frame are given as follows:We introduce the following dimensionless parameters:where is the dimensionless wave number, is the Reynolds number, is the Hartmann number, is the Darcy number, is the couple stress parameter, is the local temperature Grashof number, is the local concentration Grashof number, is the Prandtl number, is the dimensionless temperature, is the dimensionless concentration, is the heat generation parameter, is the Schmidt number, and is the Soret number.

In terms of these nondimensional variables, the governing equations (6), after dropping the bars, becomeThe dimensionless velocity components in terms of stream function are related by the following relations:Using (13), the governing equations (9)â€“(12) reduced toThe dimensionless boundary conditions can be put in the formsThe dimensionless mean flow rate in fixed frame is related to the nondimensional mean flow rate in wave frame byin whichWe note that and represent the dimensionless forms of the peristaltic wallswhere , , , and satisfy the relation .

#### 3. Solution of the Problem

Assuming that the wave length of the peristaltic wave is very large compared to the width of the channel, the wave number becomes very small. This assumption is known as long wave length approximation. Since is very small, all the higher powers of are also very small. Therefore, neglecting terms containing and its higher powers from (14), we getElimination of pressure form from (19) and (20) yieldsSolving (21) and (22) with the boundary conditions (15), the temperature and concentration are obtained aswhereInserting (24) in (23), with the help of boundary conditions (15), we obtain the stream function asin whichUsing (24) and (26) in (19), the pressure gradient is given bywhereThe nondimensional expressions for the pressure difference for one wave length , the frictional forces at both walls at and at , and the heat transfer coefficients and at the right and left walls are defined as follows [32]:Using (30), the heat transfer coefficients at the right and left walls, respectively, are obtained asThe expressions for pressure rise and frictional forces at both walls at and at involve the integration of . Due to the complexity of , the analytical integration of integrals of (30) is not possible. In view of this, a numerical integration scheme is used for the evaluation of the integrals.

#### 4. Expressions for Wave Shapes

The nondimensional expressions for the four considered wave forms are given in the following.(1)Sinusoidal wave: (2)Triangular wave: (3)Square wave: (4)Trapezoidal wave:

#### 5. Results and Discussion

This section is dedicated to discussion and analysis of the velocity distribution, pumping characteristics, heat and mass characteristics, and trapping phenomena for different flow parameters.

##### 5.1. Flow Characteristics

Figures 2(a)â€“2(f) illustrate the influence of Hartmann number , Schmidt number , couple stress parameter , Darcy number , heat generation parameter , and Grashof number on axial velocity profile across the channel. From these figures, it is observed that the maximum velocities are always located near the centre of the channel and the velocity profiles are nearly parabolic in all cases. It is noted from Figure 2(a) that as Hartmann number increases, the velocity decreases near the centre of the channel and it is increased in the neighborhood of the walls. This seems realistic because the magnetic field acts in the transverse direction to the flow and magnetic force resists the flow. The similar behavior is observed in [25]. The same behavior can be seen with increasing of Schmidt number and couple stress parameter (see Figures 2(b) and 2(c)). It is observed from Figure 2(d) that increasing of Darcy number increases the velocity near the centre of the channel and decreases the velocity of the fluid near the peristaltic walls. The same trend is followed with the increasing of heat generation parameter (see Figure 2(e)). It is noticed from Figure 2(f) that, with increasing of Grashof number , the velocity at the left wall increases while a reverse trend is seen at the right wall.

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

##### 5.2. Pumping Characteristics

Figure 3 illustrates the variation of pressure gradient over one wave length . The effects of , , and on pressure gradient are displayed in Figures 3(a)â€“3(c). It can be seen from Figure 3(a) that increasing of Hartmann number increases the pressure gradient. It shows that when strong magnetic field is applied to the flow field then higher pressure gradient is needed to pass the flow. This result suggests that the fluid pressure can be controlled by the application of suitable magnetic field strength. This phenomenon is useful during surgery and critical operation to control excessive bleeding. It is also observed that increasing of and increases the pressure gradient. From Figures 3(d)â€“3(f), it is noted that with the increasing of , , and the pressure gradient decreases. It is noticed that, in the wider part of the channels and , the pressure gradient is small, so the flow can be easily passed without the imposition of large pressure gradient. However, in the narrow part of the channel the pressure gradient is large; that is, much larger pressure gradient is needed to maintain the same given volume flow rate. Figure 4 is prepared to see the behaviour of pressure gradient for different four wave forms. It is observed from Figures 4(a)â€“4(d) that, in all the wave forms, increase in decreases pressure gradient.

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

**(a) Sinusoidal wave**

**(b) Triangular wave**

**(c) Square wave**

**(d) Trapezoidal wave**

The dimensionless pressure difference per unit wave length versus time mean flow rate has been plotted in Figure 5. We split the whole region into four segments as follows: peristaltic pumping region where and and augmented pumping region when and . There is retrograde pumping region when and . Free pumping region corresponds to . The region where and is known as peristaltic pumping region. In this region, the positive value of is entirely due to the peristalsis after overcoming the pressure difference. The region where and is known as copumping or augmented pumping region. In this region, a negative pressure difference assists the flow due to the peristalsis of the walls. The region where and is called retrograde pumping region. In this region, the flow is opposite to the direction of the peristaltic motion. In the free pumping region, the flow is caused purely by the peristalsis of the walls. It is evident from Figure 5 that there is an inversely linear relation between and . From Figure 5(a), it is clear that with the increasing of , in the augmented pumping and free pumping regions, the pumping decreases, in the peristaltic pumping region, the pumping increases up to a critical value of and decreases after the critical value, and in the retrograde pumping region the pumping increases. It is observed from Figure 5(b) that, with the increasing of , the behaviour is quite opposite with . It is noticed from Figure 5(c) that in the augmented pumping region the pumping increases and in the peristaltic pumping and retrograde pumping regions the pumping decreases. Figure 5(d) depicts that, in all the pumping regions, the pumping decreases by increasing . It is noted from Figures 5(e)-5(f) that the behaviour is quite opposite with while increasing and .

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**