Open access peer-reviewed chapter

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

Written By

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

Submitted: July 28th, 2018 Reviewed: August 24th, 2018 Published: November 19th, 2018

DOI: 10.5772/intechopen.81123

From the Edited Volume

## Pattern Formation and Stability in Magnetic Colloids

Edited by Nicola?s O. Rojas

Chapter metrics overview

View Full Metrics

## 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 h is 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=0 is maintained at temperature T1, while the cooler plate y=h is at temperature T2, T1>T2 and T0 is 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 fx is defined as

Fk=1k!dkfxdxkx=x0E9

where fx is the given function and Fk is the transformed function which is also known as the spectrum of fx. The inverse differential transformation of Fk is 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 n is 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 Fk and Θk are the transformed functions of fy and θ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, , Fn and Θ4, Θ5, Θn by setting k=0,1,2,. The values of Fj for j=4,5, and Θj for j=2,3, can now be evaluated in terms of F0, F1, F2, F3, Θ0 and Θ1. For convenience, the values of F0, F1, F2, F3, Θ0 and Θ1 are set as unknowns such as

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

The values of Fk and Θk for 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=1 is 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) kT02dTdy2 is entropy generation due to heat transfer; the next term μT0dudy2 is entropy generation due to viscous dissipation, 1a2d2udy22 and σB02T0u2 are 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 SGNs are the dimensional and dimensionless entropy generation rates.

The ratio of heat transfer entropy generation N1 to fluid friction entropy generation N2 is represented as

Φ=N2N1E20

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

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

Note that N1 represents heat transfer irreversibility, while N2 denotes 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 βvKn and Hartmann number H are 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ψ10 and the selected reference values of βvKn=0.05, ln=1.667. Furthermore, 0ξ5 and 0a2 with 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=1 but reduces towards y=0 in 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=0 is not significant; however, it is decelerated towards the microchannel plate y=1. Fluid temperature is considerably enhanced at plate y=1 for 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=1 is 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

## References

1. 1. Agboola OO, Opanuga AA, Okagbue HI, Bishop SA, Ogunniyi PO. Analysis of hall effects on the entropy generation of natural convection flow through a vertical microchannel. International Journal of Mechanical Engineering and Technology. 2008;9(8):712-721
2. 2. Arkilic EB, Schmidt MA, Breuer KS. Gaseous slip flow in long microchannels. Journal of Microelectromechanical Systems. 1997;6:167-178
3. 3. Sofonea V, Sekerka RF. Diffuse-reflection boundary conditions for a thermal lattice Boltzmann model in two dimensions: Evidence of temperature jump and slip velocity in microchannels. Physical Review E – Statistical Nonlinear and Soft Matter Physics. 2005;71:1-10 http://dx.doi.org/10.1103/PhysRevE.71.066709
4. 4. Lv Q, Liu X, Wang E, Wang S. Analytical solution to predicting gaseous mass flow rates of microchannels in a wide range of Knudsen numbers. Physical Review E. 2013;88:013007 http://dx.doi.org/10.1103/PhysRevE.88.013007
5. 5. Khadrawi AF, Al-Shyyab A. Slip flow and heat transfer in axially moving micro-concentric cylinders. International Communications in Heat and Mass Transfer. 2010;37:1149-1152 http://dx.doi.org/10.1016/j.icheatmasstransfer.2010.06.006
6. 6. Chen S, Tian Z. Simulation of thermal micro-flow using lattice Boltzmann method with Langmuir slip model. International Journal of Heat and Fluid Flow. 2010;31:227-235 http://dx.doi.org/10.1016/j.ijheatfluidflow.2009.12.006
7. 7. Larrode FE, Housiadas C, Drossinos Y. Slip-flow heat transfer in circular tubes. International Journal of Heat and Mass Transfer. 2000;43:2669-26680
8. 8. Adesanya SO. Free convective flow of heat generating fluid through a porous vertical channel with velocity slip and temperature jump. Ain Shams Engineering Journal. 2015;6:1045-1052. DOI: http://dx.doi.org/10.1016/j.asej.2014.12.008
9. 9. Aziz A, Niedbalski N. Thermally developing microtube gas flow with axial conduction and viscous dissipation. International Journal of Thermal Sciences. 2011;50:332-340 http://dx.doi.org/10.1016/j.ijthermalsci.2010.08.003
10. 10. Jha BK. Aina B. MHD natural convection flow in a vertical porous microchannel formed by nonconducting and conducting plates in the presence of induced magnetic field. Heat Transfer Research. 2017;48(15):1-24
11. 11. Jha BK, Aina B, Joseph SB. Natural convection flow in a vertical microchannel with suction/injection. Proceedings of the Institution of Mechanical Engineers, Part E: Journal of Process Mechanical Engineering. 2014;228(3):171-180. DOI: 10.1177/0954408913492719
12. 12. Jha BK, Aina B, Ajiya AT. MHD natural convection flow in a vertical parallel plate microchannel. Ain Shams Engineering Journal. 2015;6:289-295
13. 13. Bejan A. Entropy Generation through Heat and Fluid Flow. New York: Wiley; 1982
14. 14. Adesanya SO, Falade JA, Jangili S, Beg OA. Irreversibility analysis for reactive third-grade fluid flow and heat transfer with convective wall cooling. Alexandria Engineering Journal. 2017;56:153-160
15. 15. Adesanya SO, Kareem SO, Falade JA, Arekete SA. Entropy generation analysis for a reactive couple stress fluid flow through a channel saturated with porous material. Energy. 2015;93:1239-1245
16. 16. Jangili S, Adesanya SO, Falade JA, Gajjela N. Entropy generation analysis for a radiative micropolar fluid flow through a vertical channel saturated with non-Darcian porous medium. International Journal of Applied and Computational Mathematics. 3(4):3759-3782. DOI: 10.1007/s40819-017-0322-8
17. 17. Adesanya SO, Makinde OD. Irreversibility analysis in a couple stress film flow along an inclined heated plate with adiabatic free surface. Physica A: Statistical Mechanics and its Applications. 2015;432:222-229
18. 18. Adesanya SO, Makinde OD. Effects of couple stresses on entropy generation rate in a porous channel with convective heating. Computational & Applied Mathematics. 2015;34:293-307. DOI: 10.1007/s40314-014-0117-z
19. 19. Ajibade AO, Jha BK, Omame A. Entropy generation under the effect of suction/injection. Applied Mathematical Modelling. 2011;35:4630-4646
20. 20. Eegunjobi AS, Makinde OD. Effects of Navier slip on entropy generation in a porous channel with suction/injection. Journal of Thermal Science and Technology. 2002;7(4):522-535
21. 21. Eegunjobi AS, Makinde OD. Combined effect of buoyancy force and Navier slip on entropy generation in a vertical porous channel. Entropy. 2012;14:1028-1044. DOI: 10.3390/e14061028
22. 22. Das S, Jana RN. Entropy generation due to MHD flow in a porous channel with Navier slip. Ain Shams Engineering Journal. 2014;5:575-584
23. 23. Opanuga AA, Okagbue HI, Agboola OO, Imaga OF. Entropy generation analysis of buoyancy effect on hydromagnetic poiseuille flow with internal heat generation. Defect and Diffusion Forum. 2017;378:102-112
24. 24. Opanuga AA, Gbadeyan JA, Iyase SA. Second law analysis of hydromagnetic couple stress fluid embedded in a non-Darcian porous medium. IAENG International Journal of Applied Mathematics. 2017;47(3):287-294
25. 25. Opanuga AA, Okagbue HI, Agboola OO. Irreversibility analysis of a radiative MHD Poiseuille flow through porous medium with slip condition. Proceedings of The World Congress on Engineering 2017; July 5–7; London, U.K
26. 26. Opanuga AA, Bishop SA, Okagbue HI, Agboola OO. Hall current and joule heating effects on flow of couple stress fluid with entropy generation. Engineering, Technology & Applied Science Research. 2008;8(3):2923-2930
27. 27. Arikoglu A. Ozkol solution of difference equations by using differential transform method. Applied Mathematics and Computation. 2006;173(1):126-136. http://dx.doi.org/10.1016/j.amc.2005.06.013
28. 28. Biazar J, Eslami M. Differential transform method for quadratic Riccati differential equation. International Journal of Nonlinear Science. 2010;9(4):444-447
29. 29. Agboola OO, Gbadeyan JA, Opanuga AA, Agarana MC, Bishop SA, Oghonyon JG. Variational iteration method for natural frequencies of a cantilever beam with special attention to the higher modes. Proceedings of The World Congress on Engineering 2017; July 5–7 London; UK
30. 30. Agboola OO, Opanuga AA, Gbadeyan JA. Solution of third order ordinary differential equations using differential transform method. Global Journal of Pure and Applied Mathematics. 2015;11(4):2511-2517
31. 31. Odibat Z. Differential transform method for solving Volterra integral equation with separable kernels. Mathematical and Computer Modelling. 2008;48:144-149. http://dx.doi.org/10.1016/j.mcm.2007.12.022
32. 32. Opanuga AA, Agboola OO, Okagbue HI. Approximate solution of multipoint boundary value problems. Journal of Engineering and Applied Sciences. 2015;10(4):85-89
33. 33. Opanuga AA, Okagbue HI, Agboola OO, Bishop SA. Second law analysis of ion slip effect on MHD couple stress fluid. International Journal of Mechanics. 2018;12:96-101
34. 34. Chen CK, Weng HC. Natural convection in a vertical microchannel. Journal of Heat Transfer. 2005;127:1053-1056

Written By

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

Submitted: July 28th, 2018 Reviewed: August 24th, 2018 Published: November 19th, 2018