Open access peer-reviewed chapter

# Convection Flow of MHD Couple Stress Fluid in Vertical Microchannel with Entropy Generation

By Abiodun A. Opanuga, Olasunmbo O. Agboola, Hilary I. Okagbue and Sheila A. Bishop

Submitted: July 9th 2018Reviewed: August 24th 2018Published: November 19th 2018

DOI: 10.5772/intechopen.81123

## Abstract

Entropy generation of fully developed steady, viscous, incompressible couple stress fluid in a vertical micro-porous-channel in the presence of horizontal magnetic field is analysed in this work. The governing equations for the flow are derived, and nondimensionalised and the resulting nonlinear ordinary differential equations are solved via a rapidly convergent technique developed by Zhou. The solution of the velocity and temperature profiles are utilised to obtain the flow irreversibility and Bejan number. The effects of couple stresses, fluid wall interaction parameter (FSIP), effective temperature ratio (ETR), rarefaction and magnetic parameter on the velocity profile, temperature profile, entropy generation and Bejan number are presented and discussed graphically.

### Keywords

• microchannel
• entropy generation
• MHD
• natural convection
• differential transform method (DTM)

## 1. Introduction

In the last decades, the study of microchannel flows has become an important subject for researchers because of the reduction in the size of such devices which increases the dissipated heat per unit area. The effective performance of these devices is dependent on the temperature; as a result a comprehensive knowledge of such flow behaviours is required for accurate prediction of performance during the design process. These microfluidics have characteristic lengths of 1100μm and are categorised by the dimensionless quantity called the Knudsen number Kn. Researchers have shown that microchannel flows are influenced by several parameters of which velocity slip and temperature jump occurring at the solid-fluid interface in small-scale systems are the most important [2, 3, 4]. Velocity slip and temperature jump become more significant at higher Knudsen number. The latter boundary conditions are assumed at Knudsen number greater than 0.01 since below this value the classical Navier-Stokes equations is no longer valid.

The influence of velocity slip and temperature jump on microchannel flows has been extensively studied. Khadrawi and Al-Shyyab [5] obtained a close form solution of the effect of velocity slip and temperature jump on heat and fluid flow for axially moving micro-concentric cylinders. Chen and Tian [6] applied lattice Boltzmann numerical technique with Langmuir model for the velocity slip and temperature jump to investigate fluid flow and heat transfer between two horizontal parallel plates. The study suggested the application of Langmuir-slip model as an alternative for the Maxwell-slip model. Earlier on, Larrode et al. [7] in his work, slip-flow heat transfer in circular tubes, has stated the significance of fluid-wall interaction. Moreover, Adesanya [8] has analytically studied the effects of velocity slip and temperature jump on the unsteady free convective flow of heat-generating/heat-absorbing fluid with buoyancy force. The study concluded that an increase in the slip parameter enhanced fluid motion, while an increase in the temperature jump parameter increased fluid temperature. Aziz and Niedbalski [9] applied finite difference technique to investigate thermally developing microtube gas flow with respect to both radial and axial coordinates. The result indicated that increase in Knudsen number reduced local Nusselt number (Nu). Several investigations regarding microchannel flows have also been carried out by Jha and collaborators [10, 11, 12].

From the perspective of energy management, it has been established that thermal processes such as the microscale fluid flow and heat transfer modelling are irreversible. This implies that entropy generation which destroys the available energy leading to inefficiency of thermal designs exists. Therefore, the aim of this research endeavour is the minimisation of irreversibility associated with microchannel flow of couple stress fluid by applying the robust approach proposed by Bejan [13] and applied by numerous researchers such as Adesanya and collaborators [14, 15, 16], Adesanya and Makinde [17, 18], Ajibade and Jha [19], Eegunjobi and Makinde [20, 21], Das and Jana [22] and more recently Opanuga and his collaborators [23, 24, 25, 26].

In this work, an efficient technique introduced by Zhou has been employed to construct the solutions of the velocity and temperature profiles. Due to the accuracy of this technique in handling numerous linear and nonlinear models of both ordinary and partial equations, it has gained wide applications by investigators over the last decades. Arikoglu and Ozkol [27] applied it to obtain the solution of difference equations. Biazar and Eslami [28] solved quadratic Riccati differential equation using this method. Agboola et al. [29, 30] applied it to third-order ordinary differential equations and natural frequencies of a cantilever beam. Solutions of Volterra integral equation via this technique was obtained by Odibat [31], while Opanuga et al. [32] compared the method with Adomian decomposition method to find the solution of multipoint boundary-value problem and more recently in couple stress fluid model [33].

## 2. Problem formulation

A fully developed laminar, viscous, incompressible and electrically conducting couple stress fluid in a vertical parallel microchannel of width his considered. The x-axis is such that it is vertically upward along the plates while the y-axis is taken normal to it. There is an asymmetric heating of the plates such that the hotter plate y=0is maintained at temperature T1, while the cooler plate y=his at temperature T2, T1>T2and T0is the reference frame. Velocity slip and temperature jump are incorporated. Furthermore, induced magnetic field effect arising due to the motion of an electrically conducting fluid is taken into consideration (Figure 1). Following Cheng and Weng [34], the velocity slip and temperature jump are given as

us=2σvσvλdudyy=h,TsTw=2σtσt2γγ+1λPrdTdyy=hE1

While the governing equations for the flow are stated as [10]

ρνdudy=dpdx+μd2udy2ηd4udy4+ρgβTT0σB02uE2
ρcpv0dTdy=kd2Tdy2+μdudy2+ηd2udy22+σB02u2E3

The boundary conditions are

u=2fvfvλdudy,d2udy2=0,T=T2+2ftft2γγ+1λPrdTdy,y=0u=2fvfvλdudy,d2udy2=0,T=T12ftft2γγ+1λPrdTdy,y=hE4

The dimensionless parameters for this flow are

y=yh,u=uν,θ=TT0T1T0,Re=ρv0hμ,G=h2μν0dpdx,a2=μh2η,Pr=νρcpk,Br=μν0kT1T0,Bi1=γ1hk,Bi2=γ2hk,ν=μρ,Ω=T1T0T0,H2=σB02h2μ,ξT2T0=T1T0,βv=2fvfv,βt=2fvfv2γsγs+11Pr,γs=cpcv,Kn=λh,ψ=βtβv,Ns=SGT0h2kT1T0E5

Using (5) in Eqs. (2)(4) yields the boundary-value problems:

Redudy=G+d2udy21a2d4udy4H2u+Grθ,E6
d2θdy2=Re prdyBrdudy2+1a2d2udy22+H2u2,E7
u=βvKndudy,d2udy2=0,θ=ξ+βvKnψdy,y=0u=βvKndudy,d2udy2=0,θ=1βvKnψdy,y=1E8

## 3. Method of solution

DTM is applied in this work to obtain the solution of the velocity and temperature profiles. The results are used to calculate the entropy generation and irreversibility ratio.

### 3.1 Differential transformation method (DTM)

Consider a function, fx. The differential transformation of the function fxis defined as

Fk=1k!dkfxdxkx=x0E9

where fxis the given function and Fkis the transformed function which is also known as the spectrum of fx. The inverse differential transformation of Fkis given by

fx=k=0xkFkE10

In practise, the function stated in Eq. (10) is usually represented by a finite series of the form:

fx=k=0nxkFkE11

where nis the size of the series (Arikoglu and Ozkol [27], Biazar and Eslami [28]).

To apply DTM to the problem in view, the basic properties of DTM, which are outlined in Table 1, are invoked in Eqs. (6)(8). Doing this, one obtains the following recurrence relations:

Fk+4=1a2k+4!k+1k+2Fk+2Rek+1Fk+1H2Fk+GE12
Θk+2=1k+2![RePrk+1Θk+1Brr=0kr+1Fr+1kr+1Fkr+1Bra2r=0kr+1r+2Fr+2kr+1kr+2Fkr+2BrH2r=0kFrFkrE13
Original functionTransformed function
fy=uy±wyFk=Uk±Wk
fy=dnuydynFk=k+n!k!Uk+n
fy=u2Fk=r=0kUrUkr
fy=duydy2Fk=r=0kr+1kr+1Ur+1Ukr+1
fy=d2uydy22Fk=r=0kr+1r+2kr+1kr+2Ur+2Ukr+2

### Table 1.

Operations and properties of differential transform method.

where Fkand Θkare the transformed functions of fyand θy, respectively. These are given by

fy=k=0ykFk,θy=k=0ykΘk.E14

Using the recursive relations of Eqs. (12) and (13), we can obtain the differential coefficients, F4, F5, , Fnand Θ4, Θ5, Θnby setting k=0,1,2,. The values of Fjfor j=4,5,and Θjfor j=2,3,can now be evaluated in terms of F0, F1, F2, F3, Θ0and Θ1. For convenience, the values of F0, F1, F2, F3, Θ0and Θ1are set as unknowns such as

F0=a1,F1=a2,Θ0=b1,Θ1=b2.F2=a3,F3=a4,E15

The values of Fkand Θkfor k=0,1,are now substituted back into Eq. (14) to obtain the series solutions in the form:

fy=k=0nykFk,θy=k=0nykΘk.E16

We next invoke the transformed form of boundary condition (8) on (16) to determine the values of all the unknown coefficients stated in (15). Coding Eqs. (12)(15) in symbolic Maple software yields the approximate solution. The results are presented in Figures 26.

To verify the accuracy of the results, the exact solution of the velocity profile (6) subject to the boundary conditions (8) at βvKn=0.05, Gr=0, a=1, Re=0.1, H=1is obtained as.

uy=(3.063639e0.866146y+2.057901Cos0.470472y+1.e1.732291yCos0.528272y+1.808594Sin0.4704729y0.093418e1.732291ySin0.528272y)E17

The above solution is compared with DTM solution as displayed in Table 2.

βvKn=0.05, Gr=0, a=1, Re=0.1, H=1
yExactDTM
00.0018728118137910.001872811814000
0.10.0055463525402520.005546352541866
0.20.0088196333394400.008819633342380
0.30.0113792862670870.011379286271194
0.40.0130037218376600.013003721842694
0.50.0135607413095680.013560741315200
0.60.0130061534362520.013006153442060
0.70.0113833862157980.011383386221257
0.80.0088240853939640.008824085398432
0.90.0055496925550640.005549692557769
10.0018739964848020.001873996484819

### Table 2.

Comparison of the exact solution with the values of velocity (u).

### 3.2 Analysis of entropy generation

The local entropy generation for the flow is given as Bejan [13]:

SG=kT02dTdy2+μT0dudy2+1a2d2udy22+σB02T0u2.E18

The first term in Eq. (18) kT02dTdy2is entropy generation due to heat transfer; the next term μT0dudy2is entropy generation due to viscous dissipation, 1a2d2udy22and σB02T0u2are couple stress and magnetic entropy generation.

Using (4) in Eq. (18), the dimensionless form of entropy generation is written as

Ns=dy2+BrΩdudy2+BrΩa2d2udy22+BrM2Ωu2,E19

where SGNsare the dimensional and dimensionless entropy generation rates.

The ratio of heat transfer entropy generation N1to fluid friction entropy generation N2is represented as

Φ=N2N1E20

Alternatively, Bejan number gives the entropy generation distribution ratio parameter; it represents the ratio of heat transfer entropy generation N1to the total entropy generation Nsdue to heat transfer and fluid friction; it is defined as

Be=N1Ns=11+Φ,E21
Be=0,N2N10.5,N1=N21,N2N1E22

Note that N1represents heat transfer irreversibility, while N2denotes irreversibility due to viscous dissipation, couple stresses and magnetic field.

## 4. Results and discussion

In this work, investigation has been conducted on fully developed, steady, viscous and incompressible flow of couple stress fluid in a vertical micro-porous-channel in the presence of magnetic field. Effects of couple stress parameter a, fluid wall interaction parameter ψ, effective temperature ratio (ETR) ξ, rarefaction βvKnand Hartmann number Hare presented in this section. Reasonable intervals for the above parameters as used by Chen and Weng [32] are adopted in this investigation: 0H10, 0vKn0.1, 0ψ10and the selected reference values of βvKn=0.05, ln=1.667. Furthermore, 0ξ5and 0a2with reference values of a=1, ξ=0.5.

### 4.1 Influence of couple stress parameter

Figure 2A illustrates the influence of couple stress parameter on fluid velocity. The plot shows a significant reduction in fluid velocity as the values of couple stress parameter increase. This observation is due to the increase in the dynamic viscosity of the fluid. In Figure 2B, it is illustrated that fluid temperature drops as couple stress parameter increases. It is depicted in Figure 2C that fluid entropy generation decreases at the microchannel walls as couple stress parameter increases. However, the effect is opposite near the middle of the microchannel. The same scenario is observed in Figure 2D; it is shown that Bejan number reduces in value at microchannel right wall which is an indication that fluid friction irreversibility is the major contributor to entropy generation as couple stress parameter increases.

### 4.2 Influence of fluid-structure interaction parameter

Figure 3A presents the effect of fluid-structure interaction parameter on fluid velocity. It is noted that FSIP does not have any significant effect on the slip velocity at the walls as well as the microchannel. However, Figure 3B reveals that fluid temperature is enhanced as FSIP increases. Response to the enhancement in fluid temperature in Figure 3B is the rise in fluid entropy generation at the hotter wall of the microchannel, while entropy generation is reduced at the cooler wall as depicted in Figure 3C. In Figure 3D, it is noticed that Bejan number increases at the hotter wall of the microchannel, while it reduces at the cooler wall. The implication of the latter is that heat transfer irreversibility is the cause of entropy generation at the hotter wall, while on the other hand, fluid friction irreversibility is the major contributor at the cooler wall of the microchannel.

### 4.3 Influence of effective temperature ratio

Next is the response of fluid velocity, fluid temperature and entropy generation to variation in effective temperature ratio. In Figure 4, it is observed that velocity is accelerated slightly except towards the microchannel centre and plate y=0. In Figure 4B, fluid temperature is significantly enhanced at increasing values of effective temperature ratio. The effect of this is noticed in Figure 4C with an increase in entropy production approaching the microchannel plate y=1. Furthermore, Bejan number rises approaching the plate y=1but reduces towards y=0in Figure 4D. It is then concluded that fluid friction irreversibility is dominant at plate y=0, while heat transfer irreversibility is dominant at plate y=1.

### 4.4 Influence of rarefaction

In Figure 5A, the effect of rarefaction parameter on fluid velocity is presented. The plot indicates that rarefaction parameter increases and fluid velocity at plate y=0is not significant; however, it is decelerated towards the microchannel plate y=1. Fluid temperature is considerably enhanced at plate y=1for different values of rarefaction as displayed in Figure 5B. Figures 5C and 5D presents similar results at plate y=1. Both entropy generation and Bejan number reduce as the values of rarefaction parameter increase. The implication of this observation is that entropy generation at microchannel wall y=1is a consequence of viscous dissipation.

### 4.5 Influence of Hartmann number

Finally, the response of fluid velocity, temperature, entropy generation and Bejan number to variation in Hartmann number is presented in Figure 6. In Figure 6A, fluid velocity is found to have decelerated within the microchannel region. The observed reduction in the motion of fluid is attributed to the presence of applied magnetic field which usually induces a resistive type of force known as Lorentz force. Also the presence of Ohmic heating in the flow significantly enhanced fluid temperature as depicted in Figure 6B. Figures 6C and 6D displays increase in fluid entropy generation and Bejan number at the microchannel walls as Hartmann number increases from 1 to 5. It is deduced that fluid entropy generation is induced by heat transfer irreversibility.

## 5. Conclusions

Entropy generation of fully developed steady, viscous, incompressible couple stress fluid in a vertical micro-porous-channel in the presence of magnetic field is analysed in this work. The equations governing the fluid flow are solved via an efficient technique proposed by Zhou. Then fluid entropy generation and Bejan number are calculated by the results obtained. This work reduces to Chen and Weng [34] when Hartmann number, couple stress parameter, entropy generation and Bejan number are neglected H0a0Ns0Be0. Furthermore, it agrees with Jha and Aina [10] in the absence of couple stress parameter, entropy generation and Bejan number a0Ns0Be0. The present study is significant in the cooling of microchannel devices and conservation of useful energy. The following conclusions are made based on the results above:

1. Couple stress parameter reduces fluid velocity, velocity slip, temperature and entropy generation.

2. Increase in fluid-structure interaction parameter increases fluid temperature. However, entropy generation enhances at the hotter wall and reduces at the cooler region of the microchannel. Furthermore, fluid irreversibility is enhanced at the hot wall and increases at the cold wall.

3. Effective temperature ratio slightly enhances fluid velocity and velocity slip, raising fluid temperature significantly.

4. An increase in the values of rarefaction mainly reduces fluid velocity and velocity slip, increases fluid temperature and reduces entropy generation.

5. Hartmann number decreases fluid velocity and velocity slip, while the temperature is enhanced considerably. Entropy generation and Bejan number are enhanced at microchannel walls.

## Acknowledgments

Authors appreciate the funding provided by Covenant University, Ota.

## Nomenclature

 B0 uniform magnetic field u fluid velocity h channel width ft,fv thermal and tangential momentum accommodation coefficients, respectively Cp specific heat capacity Re Reynolds number a couple stress parameter k thermal conductivity Kn Knudsen number H Hartmann number Pr Prandtl number T temperature of fluid T0 reference temperature Br Brinkman number EG local volumetric entropy generation rate Be Bejan number Cv specific heats at constant volume Ns dimensionless entropy generation parameter Greek letters ρ fluid density βt,βv dimensionless variables γs ratio of specific heat μ dynamic viscosity ξ effective temperature ratio σ electrical conductivity Ω temperature difference η fluid particle size effect due to couple stresses ψ fluid wall interaction parameter

## How to cite and reference

### Cite this chapter Copy to clipboard

Abiodun A. Opanuga, Olasunmbo O. Agboola, Hilary I. Okagbue and Sheila A. Bishop (November 19th 2018). Convection Flow of MHD Couple Stress Fluid in Vertical Microchannel with Entropy Generation, Pattern Formation and Stability in Magnetic Colloids, Nicola?s O. Rojas, IntechOpen, DOI: 10.5772/intechopen.81123. Available from:

### chapter statistics

1Crossref citations

### Related Content

#### Pattern Formation and Stability in Magnetic Colloids

Edited by Nicolás O. Rojas

Next chapter

#### Free Convection in a MHD Open Cavity with a Linearly Heated Wall Using LBM

By Raoudha Chaabane

#### Swirling Flows and Flames

Edited by Toufik Boushaki

First chapter

#### Introductory Chapter: Swirling Flows and Flames

By Toufik Boushaki

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

View all Books