Open access peer-reviewed chapter

Convective Heat and Mass Transfer of Two Fluids in a Vertical Channel

Written By

Suresh Babu Baluguri and G. Srinivas

Submitted: 15 June 2020 Reviewed: 16 October 2020 Published: 12 December 2020

DOI: 10.5772/intechopen.94529

From the Edited Volume

Finite Element Methods and Their Applications

Edited by Mahboub Baccouch

Chapter metrics overview

775 Chapter Downloads

View Full Metrics

Abstract

A mathematical model for convective heat and mass transfer of two immiscible fluids in a vertical channel of variable width with thermo-diffusion, diffusion-thermal effects is presented. The governing boundary layer equations generated for momentum, angular momentum, energy and species concentration are solved with appropriate boundary conditions using Galeriken finite element method. The effects of the pertinent parameters are studied in detail. Furthermore, the rate of heat transfer, mass transfer and shear stress near both walls is analyzed.

Keywords

  • vertical channel
  • immiscible fluids
  • finite element method
  • heat and mass transfer

1. Introduction

The developments are carried out in the field of fluid dynamics which was initiated by Euler [1] by proving his famous equations of fluid flow for ideal (inviscid) fluids. Fluid dynamics is a subset of the science that looks at the materials in motion. Hydrodynamics deals with the fluids which are in motion. Fluid dynamics comes under science of fluid mechanics along with the other subcategories as fluid statics, which corresponds to fluids at rest, while fluid dynamics includes fluids in motion. Fluid is defined as a matter in a gas or liquid state. Fluid dynamics is governed by the regulations of preservation of mass, energy and linear momentum. These laws state that the total amount in a closed model remains unchangeable and the energy and mass cannot be formed or demolished. They can deform but will not disappear. Another governing law is the continuum hypothesis which defines that they are uninterrupted and their characteristics fluctuate all over. The history of fluid dynamics can be found in Rouse and Ince [2] and Tokaty [3].

1.1 Micropolar fluid

The subject of micro-polar fluids attained higher degree by many researchers because when the fluid is with the suspended particles we cannot analyze the properties of fluid flow by regular Newtonian fluid characteristics. Generally this fluid is defined as non-Newtonian comprising of small firm cylindrical matters, polymer liquids, liquefied suspensions, animal blood and such related components. The existence of dust or smoke especially, gas are characterized as micro-polar fluids in fluid dynamics. Eringen [4, 5] had taken the initiation in describing the subject of micropolar fluids. In his theory he considered the local impacts emerging from the micro-structure and the inherent movement of the fluid elements. Peddieson and McNitt [6], Ariman et al. [7] addressed many investigations and applications of micropolar fluid mechanics, which are also described in the works of Lukaszewicz [8] and Eringen [9]. A.J. Chamkha et al. [10] analyzed the wholly established free convection of micropolar fluid in a upright passage.

1.2 Magnetohydrodynamics (MHD)

MHD is the science concerned with the motions of electro fluids and their interactions with magnetic fields. It is a vital branch and comparatively new in the field of fluid dynamics. When a conducting fluid is moving along a magnetic field, it results in induction of an electric field current which in turn produces the body forces. According to Faraday’s principle, on passage of electric current in a magnetic area, it experiences a force making to direct it at right angles to the electric field. Similarly, if the conductor has electromagnetic forces of the same order as the hydrodynamical and inertial forces, these forces are taken in the equation of motion along with the other forces. The integration of Navier–Stokes relations of fluid dynamics and Maxwell’s expressions of electromagnetism describe magneto hydrodynamics which are to be solved simultaneously. There are many scientific & technical applications in the literature: heating and flow control in metal structures, power production from 2-phase models or seeded high temperature gases, magnetic constraints of extreme temperature plasma and dynamo that develop magnetic field in environmental matters.

The concept of MHD flow of the boundary layer in a vertical channel is greatly considered in present metallurgical and metal processing fields. Most of the metallic materials are manufactured from the molten state. It is significant to determine the heat transfer in metals, which are electric conductors. Therefore, a controlled cooling system is required, so that, it can be regulated through an external magnetic field.

1.3 Convective heat and mass transference

Convection is the movement of molecules within the fluids. It belongs to the fundamental means of heat and mass transference that is carried out by ways of diffusion and random Brownian movement of distinct liquid elements. In our context, convection refers to the totality of advective and diffusive transfer. However, it is taken for only advective phenomena. A mechanism of transfer of heat occurring due to bulk motion of fluids is regarded as convective heat transfer. Emphasis is given to heat that is being passed and distributed.

Extensive research has been done over convective heat and mass transference of the fluid flow in vertical channels and other geometries. The existence of temperature and concentration differences or gradients lead to the convective heat and mass transfer and it is regarded as an area of study for broad examination because it is applied in several engineering issues, which are common in atmospheric buoyancy induced actions, liquid and semi-solid bodies and so on. There are quite a large number of application in the heat and mass transfer flows like, rocket nozzles, nuclear power plants, air craft and its re-entry in atmosphere, chemical and process instruments, mist formation and dispersal, temperature and humidity circulation over cultivation farms, plants destruction because of freezing, etc. Packham [11] considered the steady co-current motion of two immiscible viscous fluids in a parallel tube, the fluid interface being ripple-free and plane. Shail [12] considered the Hartmann flow of a conducting fluid in the pass way between two parallel insulating sheets of unbounded length, there exists a sheet of non-conductive fluid between the conductive fluid and the upper passage layer and given the conclusion that considerable increase could be attained in the conductive fluid velocity for appropriate proportions of the depth and viscosity of the two fluids. Beckermann et al. [13] conducted a numerical and experimental study to analyze the fluid motion and heat transfer in a upright rectangular cover which is occupied partly with a vertical layer of a fluid-saturated absorbent structure, where it is determined that the fluid quantity entering the fluid area to the absorbent layer is dependent on the Darcy & Rayleigh numbers. Lohrasbi and Sahai [14] researched in 2-phase MHD flow and heat transfer with the 1-phase conductive fluid.

1.4 Viscous dissipation

Viscous dissipation relates to the conversion of kinetic to internal energy (heating up the fluid) with respect to viscosity. It plays a significant part in normal convection in numerous units which hold huge deviations of gravitational force Gebhart [15]. Gebhart and Mollendorf [16] examined viscous dissipation in peripheral normal convection by taking in account of exponential deviation of wall temperature using resemblance relation. Fand and Brucker [17] stated that the impact of viscous dissipation is important in case of normal/natural convection in absorbent structure with respect to their investigational correlation for the heat transference in peripheral motions. Fand et al. [18] validated the comment for the Darcy method by experimental and analytical means when predicting the heat transfer coefficient from a parallel chamber implanted in an absorbent medium. Viscous dissipation performs as a heat source and heats the medium considerably. Nakayama and Pop [19] evaluated the influence of viscous dissipation on the Darcy’s free convection towards an arbitrary shaped non-isothermal matter placed in a permeable medium. Murthy and Singh [20] observed viscous dissipation on non-Darcy normal convection from an erect flat sheet in a permeable medium saturated with Newtonian fluid. It is deduced that heat transfer decreases significantly with the presence of viscous dissipation effect. El-Amin [21] analyzed the impact of viscous dissipation and Joule heating on magneto fluid dynamics forced convection jointly on a non-isothermal straight container fixed in a fluid saturated permeable membrane. Bejan [22] defined that the calculations are limitedin examining the dissipation effect by means of a stable, 1-D energy relation, based on the analogical form with viscous dissipation effect. Pantokratoras [23] evaluated the viscous dissipation effects in a normal convection using a warmed straight plate. Seddeek [24] investigated viscous dissipation effect and thermophoresis on Darcy Forchheimer mixed convection in a fluid saturated permeable medium. Duwairi et al. [25] studied the effects of viscous dissipation and Joule heating employing an isothermal cone in a saturated porous medium. Various non-Newtonian fluids have high viscosity because the irreparable criterion owing to viscous dissipation sometimes becomes vital. Hence it motivates investigators to analyze the effects of viscous dissipation in a non-Newtonian fluid saturated permeable medium. Cortell [26] analyzed viscous dissipation effect and thermal boundary layer radiation on a nonlinear wide plate. Kairi and Murthy [27] analyzed the viscous dissipation impact over normal convection heat and mass transference from an upright cone in a non-Newtonian fluid saturated non-Darcy absorbent structure. Cortell [28] analyzed the influences of suction, viscous dissipation and thermal radiation over heat transfer of a power-law fluid past a boundless permeable sheet.

1.5 Diffusion effects

The diffusion effects namely thermal-diffusion (Soret) and diffusion-thermo (Dufour) are highly important in fluid mechanics. Soret is the transfer of mass formed by temperature gradients, i.e. species diversity evolving in a primary homogeneous blend directed to a thermal gradient. Diffusion-thermo effect is the heat transfer or the heat flux formed by concentration gradient. The problems concerned to heat and mass transference and density variations with temperature and concentration lead to integrated buoyancy convected force. The diffusion impacts influence the flow field in boundary membrane on an upright channel.

Chapman and Cowling [29] developed the diffusion-thermo and thermal-diffused heat and mass effects. Eckert and Drake [30] suggested that Dufour effect has widened magnitude and so this effect should not be ignored. Kafoussias and Williams [31] included the boundary layer flows with Soret and Dufour effects for the combined forced-normal convection problem. Anghel et al. [32] analyzed the Dufour and Soret effects of a free convection boundary layer on a vertical field inserted in a permeable membrane. Postelnicu [33] evaluated the effects of thermal-diffusion and diffusion-thermo on combined heat and mass transference in natural convection boundary layer flow in a Darcian porous media under transverse magnetic effect. Alam and Rahman [34] analyzed the effects of thermal-diffusion and diffusion-thermo on combined and free convection heat and mass transference flow past an erect permeable flat sheet inserted in a porous membrane with or without flexible suction. In many studies, Dufour and Soret effects are ignored based on a minor magnitude order than Fourier’s and Fick’s laws effect. The effect of thermal-diffusion and diffusion-thermo influences over the motion area in mixed convection boundary-layer on an upright surface kept in a permeable medium and on mixed convection flow past a vertical permeable even sheet with varying suction. Chamkha and Ben-Nakhi [35] studied the combined convection flow in the existence of thermal radiation with an erect porous layer embedded in an absorbent media considering thermal-diffusion and diffusion-thermo effects. El-Aziz [36] examined the Dufour and Soret effects on MHD heat and mass transference on a porous widening layer in the presence of thermal radiation in a combined manner. Maleque [37] considered only the diffusion-thermo effect on convective heat and mass transference past a rotating permeable disk, in where the thermal-diffusion effect is ignored. Anwar Beg et al. [38] described the thermal-diffusion and diffusion-thermo impacts by numerically studying the free convection MHD heat and mass transfer over a stretching layer with saturated permeable structure. Pal and Chatterjee [39] study shows combined convected magneto hydrodynamic heat and mass transference past a stretching plate considering Ohmic dissipated thermal-diffusion and diffusion-thermo impacts with micro-polar fluid. MHD flow of a pair of immiscible and conducting fluids within isothermal and insulated moving sheets under an applied electric and inclined magnetic effect and with an induced magnetic field has been investigated by Stamenkovic et al. [40].

1.6 Two fluid flow

For many years, Scientists and Engineers have been showing interest in two phase flows, which arise in many industrial applications. The two-phase fluid flow phenomena are important in pipe flows, fluidized beds, sedimentation, gas purification, transport processes and shock waves. The study of dynamics of two phase fluid system is concerned with the motion of a liquid or gas containing immiscible, suspended stokesian solid particles. In the equations of motion of two phase fluid flows, which are the modified form of Navier-Stoke’s equations, the presence of dust adds an extra force term which represents the interaction between the dust and the fluid particles. The modified form of Navier–Stokes equations coupled with Euler equations of motion for perfect fluids are used as the equations of motion of fluid phase and particle phase, respectively. Practical application of these flows may be found in heat exchanges utilizing liquid metal or liquid sodium coolants in the area of thermal instability in boiling heat transfer studies. Malashetty and Leela [41, 42] studied two-fluid flow and heat transference in a parallel fluid passage in both conductive phases. Such investigations are beneficial to understand the slag layer effects over heat transfer features of a coal-fired magneto hydrodynamic generator. Vajravelu et al. [43] dealt with the hydromagnetic unstable motion of two immiscible conducting fluids between two porous media of different porosity. Malashetty and Umavathi [44] studied 2-phase magnetohydrodynamic flow and heat transference in a sloped passage, where 1-phase is conductive and the transport characteristics of the fluids are assumed to be unvarying. Srinivasan et al. [45] theoretically studied the two immiscible fluid models in a permeable membrane by considering the impacts of non-Darcian boundary and inertia. Malashetty et al. [46] explored the complexities of completely established 2-fluid magneto hydrodynamic flow including and excluding applied electric field in an slant pass way and described the solution of energy and momentum equations, using perturbation method for smaller value product of Prandtl and Eckert number in completely progressed free convection 2-fluid MHD flow of a tilted passage.

Advertisement

2. Mathematical formulation

The two infinite plates are kept at Y=h1 and Y=h2 initially as shown in Figure 1 and the two sheets are isothermal with dissimilar temperatures T1 and T2 respectively. The distance h1Y0 represents region-1 and the distance 0Yh2 represents region-2 where the first one occupies micropolar fluid and the other, viscous fluid. Here the buoyancy force determines the fluid flow.

Figure 1.

Schematic diagram of the problem.

The governing conditions for the problem are developed based on the assumptions stated below:

  1. The flow is assumed to be 1-D, steady, laminar, immiscible and incompressible.

  2. The transport characteristics of the two fluids are presumed to be constant.

  3. The fluid flow is fully developed.

  4. The temperature and heat flows are continuous at the interface.

  5. T1>T2,C1>C2.

Advertisement

3. Governing equations

The governing equations which are derived in Chapter-II under the above assumptions yields.

Region-1:

dU1dY=0Lawof conservation of massE1
ρ1=ρ01β1TT1T0β1CC1C0PhysicalstateE2
μ1+Kρ1d2U1dY2+Kρ1dndY+gβ1TT1T0+gβ1CC1C0σB02U1ρ1=0MomentumE3
γd2ndY2K2n+dU1dY=0Conservation of Angular MomentumE4

where γ=μ1+K2j

k1ρ1Cpd2T1dY2+1ρ1Cpμ1dU1dY2+ρ1D1KT1CS1d2C1dY2=0EnergyE5
D1d2C1dY2+D1KT1TMd2T1dY2=0DiffusionE6

Region-2:

dU2dY=0ContinuityE7
ρ2=ρ01β2TT2T0β2CC2C0StateE8
μ2ρ2d2U2dY2+gβ2TT2T0+gβ2CC2C0σB02U2ρ2=0MomentumE9
k2ρ2Cpd2T2dY2+1ρ2Cpμ2dU2dY2+ρ2D2KT2CS2d2C2dY2=0EnergyE10
D2d2C2dY2+D2KT2TMd2T2dY2=0DiffusionE11

The above equation models (1) to (11) are solved by the following boundary and interface parameters.

U1=0atY=h1, U2=0atY=h2,U10=U20,

T=T1atY=h1, T=T2atY=h2,T10=T20,

C=C1atY=h1, C=C2atY=h2,C10=C20,

n=0atY=h1,μ1+KdU1dY+Kn=μ2dU2dYatY=0,

dndY=0atY=0, k1dT1dY=k2dT2dYatY=0,D1dC1dY=D2dC2dYatY=0.

The following non dimensional variables form the equation systems (1) to (11) in to dimensionless form:

y=Yh1(region-1), y=Yh2 (region-2), u1=U1U0, u2=U2U0, θ1=T1T0ΔT,

θ2=T2T0ΔT, N=h1U0n, j=h2 (Characteristic length), K=Kμ1, c1=C1C0ΔC, c2=C2C0ΔT, Gr=gβ1TΔTh13ν12, Gc=gβ1CΔCh13ν12, R=U0h1v1, Sr=D1KT1ΔTTMΔCU0h1, Sc=v1D1, Du=D1KT1ΔCCPCS1ν1ΔT, M=σB02h12μ1, Pr=μ1Cpk1, Ec=U02CpΔT,

CS=CS1CS2, KT=KT1KT2, D=D1D2, h=h1h2, m=μ1μ2, α=k1k2, ρ=ρ1ρ2, b1=β1Tβ2T, b2=β1Cβ2C,ν=ν1ν2.

The dimensionless forms of governing equations thus obtained are:

Region-1:

d2Ndy22K2+K2N+du1dy=0E12
1+Kd2u1dy2+KdNdy+GrRθ1+GcRc1Mu1=0E13
1PrRd2θ1dy2+EcRdu1dy2+DuRd2c1dy2=0E14
1ScRd2c1dy2+Srd2θ1dy2=0E15

Region −2

d2u2dy2+mb1ρh2GrRθ2+mb2ρh2GcRc2mMh2u2=0E16
ρhα1PrRd2θ2dy2+ρhmEcRdu2dy2+cshDKTDuRd2c2dy2=0E17
hD1ScRd2c2dy2+hKTDSrd2θ2dy2=0E18

The dimensionless boundary and interface conditions thus formed are:

u1=0aty=1,u2=0aty=1,u10=u20,θ1=1aty=1,θ2=0aty=1,θ10=θ20,c1=1aty=1,c2=0aty=1,c10=c20,N=0aty=1,du1dy+K1+KN=1mh1+Kdu2dyaty=0,dNdy=0aty=0,dθ1dy=1hαdθ2dyaty=0,dc1dy=1hDdc2dyaty=0.E19
Advertisement

4. Solution of the problem

The finite element method as described in Chapter-II is applied in solving the dimensionless coupled differential equations generated by the fluid flows. For the problem discussed here, it is considered that each region is classified into 100 linear elements and each element is 3 nodded.

The element equations associated with Eqs. (12) to (18) is as follows:

yiyi+1d2Ndy22K2+K2N+du1dyηkdy=0E20
yiyi+11+Kd2u1dy2+KdNdy+GrRθ1+GcRc1Mu1ηkdy=0E21
yiyi+11PrRd2θ1dy2+EcRdu1dy2+DuRd2c1dy2ηkdy=0E22
yiyi+11ScRd2c1dy2+Srd2θ1dy2ηkdy=0E23
yiyi+1d2u2dy2+mb1ρh2GrRθ2+mb2ρh2GcRc2mMh2u2χkdy=0E24
yiyi+1ρhα1PrRd2θ2dy2+ρhmEcRdu2dy2+cshDKTDuRd2c2dy2χkdy=0E25
yiyi+1hD1ScRd2c2dy2+hKTDSrd2θ2dy2χkdy=0E26

Where ηk and χk denotes the shape functions of a typical element yiyi+1 in the region 1 and 2 correspondingly.

On integrating the above equations and by replacing the finite element Galerkin calculations,

u1i=j=13ujiηji, c1i=j=13cjiηji, Ni=j=13Njiηji, θ1i=j=13θjiηji, u2i=j=13ujiχji, c2i=j=13cjiχji, θ2i=j=13θjiχji.

From Eq. (20) we get

yiyi+1dNidydηkdydy+2K2+Kyiyi+12Niηkdydηkdyu1idy=ηkdNidy+ηku1iyiyi+1

The stiffness matrix equation corresponding to the above is

akjiNki+bkjiuki=Q1jiE27

where akji=yiyi+1dηkdydηjidydy+2K2+Kyiyi+12ηkηjidy

bkji=2K2+Kyiyi+1dηjidyηkdy
Q1ji=ηkdNidy+ηku1iyiyi+1

From Eq. (21) we get

yiyi+11+Kdηkdydu1idydy+yiyi+1KdηkdyNidyGrRyiyi+1ηkθ1idyGcRyiyi+1ηkc1idy+Myiyi+1ηku1idy=1+Kηkdu1idyKηkNiyiyi+1

The stiffness matrix equation corresponding to the above is

ckjiuki+dkjiNki+kjiθki+fkjicki=Q2jiE28

where ckji=yiyi+11+Kdηjidydηkdydy+Myiyi+1ηjiηkdy, dkji=Kyiyi+1dηjidyηkdy

kji=GrRyiyi+1ηjiηkdy,fkji=GcRyiyi+1ηjiηkdy
Q2ji=1+Kηkdu1idyKηkNiyiyi+1

From Eq. (22) we get

1PrRyiyi+1dηkdydθ1idydyEcRyiyi+1ηkdu1idy2dyDuRyiyi+1dηkdydc1idydy=1PrRηkdθ1idyDuRηkdc1idyyiyi+1

The stiffness matrix equation corresponding to the above is

gkjiθki+ukiThkjiuki+mkjicki=Q3jiE29

where

gkji=1PrRyiyi+1dηkdydηjidydy,hkji=EcRyiyi+1dη1idy2dη1idydη2idydη1idydη3idydη2idydη1idydη2idy2dη2idydη3idydη3idydη1idydη3idydη2idydη3idy2ηkdy.
mkji=DuRyiyi+1dηkdydηjidydy,Q3ji=1PrRηkdθ1idyDuRηkdc1idyyiyi+1.

From Eq. (23) we get

1ScRyiyi+1dηkdydc1idydy+Sryiyi+1dηkdydθ1idydy=1ScRηkdc1idy+Srηkdθ1idyyiyi+1

The stiffness matrix equation corresponding to the above is

nkjicki+pkjiθki=Q4jiE30

where nkji=1ScRyiyi+1dηkdydηjidydy, pkji=Sryiyi+1dηkdydηjidydy

Q4ji=1ScRηkdc1idy+Srηkdθ1idyyiyi+1

From Eq. (24) we get

yiyi+1dχkdydu2idydymb1ρh2GrRyiyi+1χkθ2idymb2ρh2GcRyiyi+1χkc2idy+mMh2yiyi+1χku2idy=χkdu2idyyiyi+1

The stiffness matrix equation corresponding to the above is

Ckjiuki+Dkjiθki+Ekjicki=Q5jiE31

where Ckji=yiyi+1dχjidydχkdydy+mMh2yiyi+1χjiχkdy, Dkji=mb1ρh2GrRyiyi+1χjiχkdy, Ekji=mb2ρh2GcRyiyi+1χjiχkdy, Q5ji=χkdu2idyyiyi+1.

From Eq. (25) we get

ρhα1PrRyiyi+1dχkdydθ2idydyρhmEcRyiyi+1χkdu2idy2dycshDKTDuRyiyi+1dχkdydc2idydy=ρhα1PrRχkdθ2idycshDKTDuRχkdc2idyyiyi+1

The stiffness matrix equation corresponding to the above is

Fkjiθki+ukiTGkjiuki+Hkjicki=Q6jiE32

where Fkji=ρhα1PrRyiyi+1dχkdydχjidydy,

Gkji=ρhmEcRyiyi+1dχ1idy2dχ1idydχ2idydχ1idydχ3idydχ2idydχ1idydχ2idy2dχ2idydχ3idydχ3idydχ1idydχ3idydχ2idydχ3idy2χkdy Hkji=cshDKTDuRyiyi+1dχkdydχjidydy, Q6ji=ρhα1PrRχkdθ2idycshDKTDuRχkdc2idyyiyi+1.

From Eq. (26) we get

hD1ScRyiyi+1dχkdydc2idydy+hKTDSryiyi+1dχkdydθ2idydy=hD1ScRχkdc2idy+hKTDSrχkdθ2idyyiyi+1

The stiffness matrix equation corresponding to the above is

Lkjicki+Mkjiθki=Q7jiE33

where Lkji=hD1ScRyiyi+1dχkdydχjidydy, Mkji=hKTDSryiyi+1dχkdydχjidydy

Q7ji=hD1ScRχkdc2idy+hKTDSrχkdθ2idyyiyi+1

The Langrange’s interpolation polynomials are used as the shape functions at each of the nodes are considered as follows:

ηi1=y2i101100y2i1001002i1021002i1011002i1021002i100100,ηi2=y2i102100y2i1001002i1011002i1021002i1011002i100100,ηi3=y2i101100y2i1021002i1001002i1021002i1001002i101100.

and similarly for χi1,χi2,χi3.

The shear stress values, heat (Nusselt number) and mass transfer rate (Sherwood number) are calculated at both walls as per the following relations:

St1=u1yy=1,St2=u2yy=1,Nu1=θ1yy=1,Nu2=θ2yy=1,Sh1=c1yy=1,Sh2=c2yy=1.
Advertisement

5. Results and discussion

The numerical solution of the system of equations is analyzed for several values of the governing factors and its corresponding graphical representations are resulted. Thermal Grashof (Gr), Molecular Grashof (Gc) and Reynolds numbers (R), Magnetic field (M) and Material parameters (K′) and Dufour (Du), Schmidt (Sc), Soret (Sr) and Eckert numbers (Ec) are fixed as Gr = 5, Gc = 5, R = 3, M = 3, K′ = 0.1, Du = 0.08, Sr. = 0.1, Sc = 0.66, Sr. = 0.001 for all the profiles excepting the varying parameter.

The profiles of all the governing parameters are depicted from Figures 210. The flow in micropolar region is found to be more than the flow in viscous region. The variations of linear momentum and angular momentum are clear for each and every governing parameter. The variation of temperature and diffusion are very narrow except for the parameters R, Du, Sr., and Sc. The temperature and diffusion are uniform across the channel and are found to be significant at the mid region of the channel. The diffusion is slightly effected at the interface due to two fluids. Hence the two fluid flow model has much importance in the real time systems. All our results are compared with earlier studies and they are validated.

Figure 2.

(a) Represents behavior of u. (b) Represents behavior of N for Gr. (c) Represents behavior of θ. (d) Represents behavior of c for Gr.

Figure 3.

(a) Represents behavior of u. (b) Represents behavior of N for Gc. (c) Represents behavior of θ. (d) Represents behavior of c for Gc.

Figure 4.

(a) Represents behavior of u. (b) Represents behavior of N for R. (c) Represents behavior of θ. (d) Represents behavior of c for R.

Figure 5.

(a) Represents behavior of u. (b) Represents behavior of N for M. (c) Represents behavior of θ. (d) Represents behavior of c for M.

Figure 6.

(a) Represents behavior of u. (b) Represents behavior of N for K′. (c) Represents behavior of θ. (d) Concentration profiles for K′.

Figure 7.

(a) Represents behavior of u. (b) Represents behavior of N for Du. (c) Represents behavior of θ. (d) Represents behavior of c for Du.

Figure 8.

(a) Represents behavior of u. (b) Represents behavior of N for Sr. (c) Represents behavior of θ. (d) Represents behavior of c for Sr.

Figure 9.

(a) Represents behavior of u. (b) Represents behavior of N for Sc. (c) Represents behavior of θ. (d) Represents behavior of c for Sc.

Figure 10.

(a) Represents behavior of u. (b) Represents behavior of N for Ec. (c) Represents behavior of θ. (d) Represents behavior of c for Ec.

Figure 2(a)(d) illustrate the effect of Grashof numbers on velocity, angular velocity, temperature and diffusion. As Gr increases the velocity and angular velocity increases substantially. The buoyancy enhances the flow in both regions i.e. thermal buoyancy force dominates the viscous force in both regions of the channel and it is found to be more in micropolar region. The lowest velocity corresponds to Gr = 2. Higher Gr values boost up the flow in both regions. As Gr increases the minute enhancement of temperature and diffusion are observed. Similar observations are noticed with all the variations of Gc which are plotted in Figure 3(a)(d).

Figure 4(a)(d) describe the Reynolds number (R) impact on velocity, angular velocity, temperature and diffusion. The reduction of velocity is found with increase of Reynolds number due to domination of inertial force on viscous force in both regions of the channel and found more drastic in viscous region. Also reduces the micro rotation with increase of Reynolds number. The effect of inertial forces enhances the temperature and reduction of the diffusion is shown with increase of R.

Figure 5(a)(d) describe the magnetic field (M) effect on velocity, angular velocity, temperature and diffusion. They portray that there could be seen reduction in velocity and angular velocity as M increases. It shows that magnetic field has a tendency to retard fluid velocity and angular velocity due to the formation of resistive Lorentz force, where when magnetic effect is applied to the fluid, it tends to retard the fluid motion. The magnetic field parametric impact on temperature and diffusion is minute.

Figure 6(a)(d) explain the effect of material parameter (K′) on velocity, angular velocity, temperature and diffusion. The effect of this parameter is very significant in both velocity and angular velocity, as K′ increases the velocity decreases significantly and this is reversed with respect to angular velocity. Minute effect is observed for both temperature and diffusion.

Figure 7(a)(d) explain the effect of Dufour number (Du) on velocity, angular velocity, temperature and diffusion. Figure 7(a) depict the fact that as Du increases, i.e. when molecular diffusivity increases, the velocity is reduced and leads to reduction of micro rotation also from Figure 7(b) it clearly indicates the influence of the concentration gradients to the thermal energy flux in the flow. Figure 7(c) specifies that the temperature reduces with increase of molecular diffusivity over the thermal diffusivity. It is clear that the diffusion profiles increase with increase of dufour number as observed in Figure 7(d).

Figure 8(a)(d) relates the Soret number (Sr) impact on velocity, angular velocity, temperature and diffusion. Figure 8(a) shows that as Sr. increases i.e. thermal diffusivity increases the decrease in velocity is found and micro rotation also decreases with increase of Sr. from Figure 8(b). Soret number states the impact of temperature gradients stimulating considerable mass diffusion effects. Here, as Soret number increases it leads to rise in temperature and shows the decay in the fluid concentration from Figure 8(c) and (d).

Figure (a)–(d) specify the Schmidt number (Sc) effect on velocity, angular velocity, temperature and diffusion. Figure 9(a) and (b) show that as Sc increases, velocity and angular velocity decreases significantly. From Figure 9(c) and (d) it is found that the temperature increases with increase of Sc and fluid concentration reduces with increase in Schmidt number.

Figure 10(a)(d) define the effect of Eckert number (Ec) on velocity, angular velocity, temperature and diffusion. From all the Figures it is concluded that the enthalpy is not having much influence over the flow for small variation of enthalpy. Figure 10(a) and (b) show that the velocity and angular velocity increase when Eckert number increases. So it is observed that momentum and angular momentum are inversely proportional to enthalpy. From Figure 10(c) the temperature increases with increase of kinetic energy. The kinetic energy reduces the concentration of the fluid as shown from Figure 10(d).

Table 1 shows the Shear stress and Nusselt and Sherwood numbers values with all the effects of all governing functions. From this table, it is observed that the absolute Shear stress enhances with increase in Gr on both the boundaries y=1 and y=1 because of buoyancy forces and similar nature is observed for Gc also. For increase of Reynolds number, magnetic field and Material parameter and Dufour, Soret and Schmidt numbers, the stress reduces on both the boundaries. This case is reversed for dissipation effect. The Nusselt number i.e. rate of heat transfer decreases on the boundary at y=1 and increases on the other boundary y=1 for the parameters Gr, Gc, R, Sr., Sc and Ec. The rise in convection is leading to reduction of heat transfer rate on the plate bounding the region −1, the reverse effect is observed for boundary of the region-2. Drastic heat transfer rate is observed for the variations of the Reynolds Number. The increase in the Reynolds number decreases the heat transfer rate on the left plate and enhances on the right plate. For the other parameters M, K′, Du the effect is reversal. The Sherwood number i.e. rate of mass transfer increase on the boundary at y=1 and decrease at the boundary y=1 for the parameters Gr, Gc, R, Sr., Sc, Ec. This is because the rise in convection and inertial forces leading to enhance the concentration. For the other parameters M, K′, Du the effect is reversal i.e. mass transfer increases at the left boundary and decreases at the right.

GrSt-ISt-IINu-INu-IISh-ISh-II
2−0.8307280.2387730.5018020.4981510.6744870.325508
5−1.213080.3679720.5017050.4981960.6745120.325473
8−1.595440.4971880.5015710.498260.6745470.325425
10−1.850370.5833440.501460.4983120.6745760.325385
GcSt-ISt-IINu-INu-IISh-ISh-II
2−0.8675710.2763850.501790.4981610.6744910.3255
5−1.213080.3679720.5017050.4981960.6745120.325473
8−1.558580.4595530.5015930.4982430.6745410.325438
10−1.788910.5206030.5015040.4982790.6745640.32541
RSt-ISt-IINu-INu-IISh-ISh-II
1−3.754111.22110.511290.4877490.5562650.443686
2−1.848760.5817030.5071410.4926270.6143050.385672
3−1.213080.3679720.5017050.4981960.6745120.325473
5−0.7032330.1956590.4898270.510140.801780.198213
MSt-ISt-IINu-INu-IISh-ISh-II
1−1.563810.6320440.501520.4983220.6745630.325375
2−1.354980.4688030.501640.4982380.674530.325441
3−1.213080.3679720.5017050.4981960.6745120.325473
5−1.027720.2522610.5017720.4981590.6744940.325503
K′St-ISt-IINu-INu-IISh-ISh-II
0.05−1.250540.3697860.5016970.4981980.6745140.325472
0.1−1.213080.3679720.5017050.4981960.6745120.325473
0.2−1.14610.3641310.5017190.4981930.6745090.325476
0.3−1.087730.3601460.5017320.498190.6745060.325478
DuSt-ISt-IINu-INu-IISh-ISh-II
0.08−1.213080.3679720.5017050.4981960.6745120.325473
0.1−1.2130.367890.5022070.4976950.6742360.325748
0.3−1.212050.3669210.508120.4918040.6709890.328973
0.5−1.210730.3655830.5162950.4837050.6664990.333392
SrSt-ISt-IINu-INu-IISh-ISh-II
0.05−1.241930.3974330.5101330.4897620.5849410.415051
0.08−1.224810.3799480.5051310.4947680.6381040.361884
0.1−1.213080.3679720.5017050.4981960.6745120.325473
0.13−1.194990.3495120.4964230.5034810.7306380.269345
ScSt-ISt-IINu-INu-IISh-ISh-II
0.22−1.251210.4069090.5128440.4870490.5561310.443864
0.66−1.213080.3679720.5017050.4981960.6745120.325473
1.02−1.179750.3339450.4919690.5079380.7779670.222014
1.5−1.131990.2851910.4780180.5218970.9261940.073786
EcSt-ISt-IINu-INu-IISh-ISh-II
0.001−1.213080.3679720.5017050.4981960.6745120.325473
0.002−1.213090.3679790.5015270.4982760.6745590.325413
0.02−1.213310.3680930.498320.4997030.6753910.324325
0.05−1.213670.3682850.4929680.5020860.6767790.322509

Table 1.

Shear stress, Nusselt number, Sherwood numbers.

Advertisement

6. Conclusions

  • Significant effect of Soret number on momentum, diffusion are observed

  • Diffusion effects are reducing the momentum and are more pronounced in concentration.

  • For enhancement of inertial force the momentum and diffusion reduces to a higher extent in case of micropolar region than viscous region.

  • The diffusion parameters are reducing the magnitude of the shear stress on both the boundaries.

  • For dufour heat transfer rate effect is enhancing on the hot plate and reducing on the cold plate but reverse effect is observed for rate on mass transfer. Exactly opposite is observed for Soret number.

  • Reduction of heat transference rate on hot plate and enhancement on the cold plate is observed due to viscous dissipation.

References

  1. 1. L. Euler, General principles of the motion of fluids, Physica. D, Vol. 237(14), 2008, pp. 1825–1839
  2. 2. H. Rouse, S. Ince, History of Hydraulics, Iowa Inst, Hydraulic Research, Vol. 269, 1957
  3. 3. G.A.Tokaty, A History and Philosophy of Fluid Mechanics, GT Foulis & Co. Ltd., Oxfordshire, 1971
  4. 4. A.C.Eringen, Simple Microfluids, International Journal of Engineering Science,Vol 2(2), 1964, pp. 205–217
  5. 5. A.C.Eringen, Theory Of Thermomicrofluids, Journal of Mathematical Analysis and Applications, Vol. 38(2), 1972, pp. 480–496
  6. 6. J. Peddieson and R.P. McNitt, Boundary layer theory for a micropolar fluid, Adv. Eng. Sci.Vol.5,1970, pp. 405
  7. 7. T.Ariman, M.A.Turk, N.D.Sylvester, Application of microcontinum fluid mechanics, International Journal of Engineering Science, Vol.12, 1974, pp. 273–293
  8. 8. G. Lukaszewicz, Micropolar fluids: theory and application. Birkhäuser, Basel. 1999
  9. 9. A.C. Eringen, Microcontinuum field theories II: fluent media. Springer, NewYork, 2001
  10. 10. A.J. Chamkha, T.Groşan and I. Pop, Fully developed free convection of a micropolar fluid in a vertical channel, International Communications in Heat and Mass Transfer, Vol 29(8), 2002, pp. 1119–1127
  11. 11. Packham, B. A. and R. Shall, Stratified laminar flow of two immiscible fluids, Mathematical Proceedings of the Cambridge Philosophical Society. Vol. 69(3), Cambridge University Press, 1971
  12. 12. R. Shail, On laminar two-phase flows in magnetohydrodynamics, International Journal of Engineering Science, Vol.11 (10), 1973, pp. 1103–1108
  13. 13. C. Beckermann, S.Ramadhyani, R. Viskanta, Natural convection flow and heat transfer between a fluid layer and a porous layer inside a rectangular enclosure, ASME Journal of Heat Transfer Vol.109(2), 1987, pp. 363–370
  14. 14. J.Lohrasbi and V. Sahai, Magnetohydrodynamic heat transfer in two-phase flow between parallel plates, Applied Scientific Research, Vol. 45(1), 1988, pp. 53–66
  15. 15. B.Gebhart, Effects of viscous dissipation in natural convection, Journal of Fluid Mechanics, Vol. 14(02), 1962, pp. 225–232
  16. 16. B. Gebhart and J. Mollendorf, Viscous dissipation in external natural convection flows, Journal of fluid Mechanics, Vol.38(01), 1969, pp. 97–107
  17. 17. M.Fand and James Brucker, A correlation for heat transfer by natural convection from horizontal cylinders that account for viscous dissipation, International Journal of Heat and Mass Transfer, Vol. 26(5), 1983, pp. 709–716
  18. 18. M.Fand, T.E.Steinberger and P.Cheng, Natural convection heat transfer from a horizontal cylinder embedded in a porous medium, International Journal of Heat and Mass Transfer, Vol. 29(1), 1986, pp.119–133
  19. 19. A. Nakayama and I. Pop, Free convection over a non isothermal body in a porous medium with viscous dissipation, International communications in heat and mass transfer, Vol.16(2),1989, pp.173–180
  20. 20. P.V.S.N.Murthy and P. Singh, Effect of viscous dissipation on a non-Darcy natural convection regime, International journal of heat and mass transfer, Vol 40(6), 1997, pp. 1251–1260
  21. 21. El-Amin. M. F., Combined effect of viscous dissipation and Joule heating on MHD forced convection over a non-isothermal horizontal cylinder embedded in a fluid saturated porous medium, Journal of Magnetism and Magnetic materials, Vol. 263(3), 2003, pp.337–343
  22. 22. Bejan.A and S.Lorente, The constructal law and the thermodynamics of flow systems with configuration, International journal of heat and mass transfer, Vol. 47(14), 2004, pp. 3203–3214
  23. 23. A.Pantokratoras, Effect of viscous dissipation in natural convection along a heated vertical plate, Applied Mathematical Modelling, Vol. 29(6), 2005, pp. 553–564
  24. 24. Seddeek. M. A. and M. S. Abdelmeguid, Effects of radiation and thermal diffusivity on heat transfer over a stretching surface with variable heat flux, Physics Letters A, Vol.348 (3), 2006, pp. 172–179
  25. 25. Duwairi.H.M, Osama Abu-Zeid and A. Damesh Rebhi, Viscous and Joule heating effects over an isothermal cone in saturated porous media, Jordan Journal of Mechanical and Industrial Engineering, Vol.1(2), 2007, pp.113–118
  26. 26. Cortell Rafael, Effects of viscous dissipation and radiation on the thermal boundary layer over a nonlinearly stretching sheet, Physics Letters A, Vol. 372(5), 2008, pp. 631–636
  27. 27. Kairi.R. R. and P. V. S. N. Murthy, Effect of viscous dissipation on natural convection heat and mass transfer from vertical cone in a non-Newtonian fluid saturated non-Darcy porous medium, Applied mathematics and computation, Vol. 217(20), 2011, pp. 8100–8114
  28. 28. Cortell Rafael, Heat and fluid flow due to non-linearly stretching surfaces, Applied Mathematics and Computation, Vol. 217(19), 2011, pp. 7564–7572
  29. 29. Chapman Sydney and Thomas George Cowling, The Mathematical Theory of Non-uniform Gases: An Account of the Kinetic Theory of Viscosity, Thermal Conduction and Diffusion of Gases, Notes Added in 1951. Cambridge university press, 1952
  30. 30. Eckert RG and Robert M Drake, Analysis of Heat Transfer, McGraw-Hill, 1972
  31. 31. Kafoussias N.G and E.W.Williams, Thermal-diffusion and diffusion-thermo effects on mixed free-forced convective and mass transfer boundary layer flow with temperature dependent viscosity, International Journal of Engineering Science, Vol. 33(9), 1995, pp. 1369–1384
  32. 32. Anghel M, H. S. Takhar and I. Pop. Dufour and Soret effects on free convection boundary layer over a vertical surface embedded in a porous medium. Studia Universitatis Babes-Bolyai, Mathematica, Vol. 45(4), 2000, pp.11–21
  33. 33. Postelnicu Adrian, Influence of a magnetic field on heat and mass transfer by natural convection from vertical surfaces in porous media considering Soret and Dufour effects, International Journal of Heat and Mass Transfer, Vol. 47(6), 2004, pp. 1467–1472
  34. 34. Alam M. S and M. M. Rahman, Dufour and Soret effects on mixed convection flow past a vertical porous flat plate with variable suction, Nonlinear Analysis: Modelling and Control, Vol.11(1),2006, pp.3–12
  35. 35. A.J. Chamkha and Abdullatif Ben-Nakhi, MHD mixed convection–radiation interaction along a permeable surface immersed in a porous medium in the presence of Soret and Dufour’s effects, Heat and Mass transfer, Vol. 44(7), 2008, pp. 845
  36. 36. El-Aziz, Thermal-diffusion and diffusion-thermo effects on combined heat and mass transfer by hydromagnetic three-dimensional free convection over a permeable stretching surface with radiation, Physics Letters A, Vol . 372(3), 2008, pp.263–272
  37. 37. Maleque Abdul, Effects of combined temperature-and depth-dependent viscosity and Hall current on an unsteady MHD laminar convective flow due to a rotating disk, Chemical Engineering Communications, Vol. 197(4), 2009, pp. 506–521
  38. 38. O.Anwar Beg, A .Y. Bakier and V. R. Prasad, Numerical study of free convection magnetohydrodynamic heat and mass transfer from a stretching surface to a saturated porous medium with Soret and Dufour effects, Computational Materials Science, Vol. 46(1), 2009, pp. 57–65
  39. 39. Pal Dulal and Sewli Chatterjee, Mixed convection magnetohydrodynamic heat and mass transfer past a stretching surface in a micropolar fluid-saturated porous medium under the influence of Ohmic heating, Soret and Dufour effects, Communications in Nonlinear Science and Numerical Simulation, Vol. 16 (3), 2011, pp. 1329–1346
  40. 40. Z.M.Stamenkovic, D.D.Nikodijevic, M.M. Kocic and J.D. Nikodijević, Mhd flow and heat transfer of two immiscible fluids with induced magnetic field effects, Thermal Science, 16(suppl. 2),2012, pp.323–336
  41. 41. M.S. Malashetty and V. Leela, Magnetohydrodynamic heat transfer in two fluid flow, Proc. of National Heat Transfer conferences sponsored by AIChE and ASME– HTD, Phase Change Heat Transfer, Vol. 159, 1991, pp.171–175
  42. 42. M.S. Malashetty and V. Leela, Magnetohydrodynamic heat transfer in two phase flow, International Journal of Engineering Science, Vol. 30, 1992, pp.371–377
  43. 43. K .Vajravelu, P.V.Arunachalam and S.Sreenadh, Unsteady flow of two immiscible conducting fluids between two permeable beds, Journal of mathematical analysis and applications, Vol.196(3), 1995, pp. 1105–1116
  44. 44. M. S. Malashetty, J. C. Umavathi, Two-phase magneto hydrodynamicflow and heat transfer in an inclined channel. Int J Multiphase Flow, Vol. 22, 1997, pp. 545–560
  45. 45. V. Srinivasan and K.Vafai, Analysis of Linear Encroachment in Two-Immiscible Fluid Systems, ASME Journal of Fluid Engineering, Vol. 116,1994, pp. 135–139
  46. 46. M. S. Malashetty, J. C. Umavathi and J. Prathap Kumar, Convective magnetohydrodynamic two fluid flow and heat transfer in an inclined channel, Heat and Mass Transfer Vol. 37(2), 2001, pp. 259–264

Written By

Suresh Babu Baluguri and G. Srinivas

Submitted: 15 June 2020 Reviewed: 16 October 2020 Published: 12 December 2020