# Double Diffusive Convection in a Layer of Maxwell Viscoelastic Fluid in Porous Medium in the Presence of Soret and Dufour Effects

Authors: Ramesh Chand; G. C. Rana
Journal: Journal of Fluids
 (2014)

---

## Abstract

Double diffusive convection in a horizontal layer of Maxwell viscoelastic fluid in a porous medium in the presence of temperature gradient (Soret effects) and concentration gradient (Dufour effects) is investigated. For the porous medium Darcy model is considered. A linear stability analysis based upon normal mode technique is used to study the onset of instabilities of the Maxwell viscolastic fluid layer confined between two free-free boundaries. Rayleigh number on the onset of stationary and oscillatory convection has been derived and graphs have been plotted to study the effects of the Dufour parameter, Soret parameter, Lewis number, and solutal Rayleigh number on stationary convection.


---

## Body

## 1. Introduction

Bénard convection originated from the experimental works of Bénard [1] and theoretical analysis of Lord Rayleigh [2]. Lord Rayleigh studied the dynamic origins of convective cells and proposed his theory on the buoyancy driven convection. The detailed study of Bénard convection in Newtonian fluid in nonporous medium under varying assumptions of hydrodynamics and hydromagnetics had been given by Chandrasekhar [3]. Lapwood [4] had studied the stability of convective flow in hydromagnetics in a porous medium using Rayleigh’s procedure. The Rayleigh instability of a thermal boundary layer in flow through a porous medium had been considered by Wooding [5]. McDonnel [6] suggested the importance of porosity in the astrophysical context.Double-diffusive convection is referred to buoyancy-driven flows induced by combined temperature and concentration gradients. The onset of double diffusive convection in a fluid saturated in porous medium is regarded as a classical problem due to its wide range of applications in many engineering fields such as evaporative cooling of high temperature systems, agricultural product storage, soil sciences, enhanced oil recovery, packed-bed catalytic reactors, and the pollutant transport in underground. A detailed review of the literature concerning double diffusive convection in binary fluid in a porous medium was given by Nield and Bejan [7], Trevisan and Bejan [8], and Malashetty and Kollur [9]. Thermal convection in binary fluid driven by the Soret and Dufour effects had been investigated by Knobloch [10] and showed that the equations were identical to the thermosolutal problem except relation between the thermal and solutal Rayleigh numbers. The above literature dealt with Newtonian fluids.The study of natural convection of non-Newtonian fluids in a porous medium had gained much attention because of its engineering and industrial applications. These applications included design of chemical processing equipment, formation and dispersion of fog, distributions of temperature and moisture over agricultural fields and groves of fruit trees, and damage of crops due to freezing and pollution of the environment.The fluids that show distinct deviation from “Newtonian hypothesis” (stress on fluid is linearly proportional to strain rate of fluid) are called non-Newtonian fluids. Different models had been proposed to explain the behavior of non-Newtonian fluids. Maxwell model is one of them. These fluids help us to understand the wide variety of fluids that exist in the physical world and characterized by power-law model. The work on viscoelastic fluid appears to be that of Herbert on plane coquette flow heated from below. He found a finite elastic stress in the undistributed state to be required for the elasticity to affect the stability. Using a three constants rheological model due to Oldroyd [11], the author demonstrated, for finite rate of strain, that the elasticity has a destabilizing effect, which results solely from the change in apparent viscosity.The importance of the study of viscoelastic fluids in a porous medium has been increasing for the last few years. This is mainly due to their applications in petroleum drilling, manufacturing of foods and paper, and many others. The problem of convective instability of viscoelastic fluid heated from below was first studied by Green [12]. Vest and Arpaci [13] investigated problems of overstability in a horizontal layer of a viscoelastic fluid heated from below. Bhatia and Steiner [14] studied the problem of thermal instability of a Maxwellian viscoelastic fluid in the presence of rotation and found that rotation has a destabilizing influence in contrast to its stabilizing effect on a viscous Newtonian fluid. Bhatia and Steiner [15] have also studied the thermal instability of a Maxwellian viscoelastic fluid in hydromagnetic and found that magnetic field had stabilizing effect on the Maxwell fluid, just as in the case of Newtonian fluid. Sharma and Kumar [16] studied the Hall effect on thermosolutal instability in a Maxwellian viscoelastic fluid and found that Hall effect destabilizes the fluid layer while Kirti and Chand [17] studied the combined effect Hall Current, suspended particles, and variable gravity in a layer of Maxwell viscoelastic fluid. Chand and Kango [18], Chand [19–21], and Chand and Kumar [22] studied problems of thermal instability of Maxwell viscoelastic fluid in porous medium under various assumptions. Chand and Rana [23] investigated the Dufour and Soret effects in layer of elasticoviscous fluid in a porous medium and found that Dufour parameter destabilizes the fluid layer while Soret parameter has both the stabilizing and destabilizing effects on fluid layer depending upon certain conditions.In this paper an attempt has been made to study the Dufour and Soret effects on the onset of instability in a horizontal layer of Maxwell viscoelastic fluid in a porous medium.


## 2. Mathematical Formulations of the Problem

Consider an infinite horizontal layer of Maxwell viscoelastic fluid of thickness “d,” confined between the planes z=0 and z=d in a porous medium of porosity ε and medium permeability k1 and is acted upon by gravity g(0,0,-g). This layer of fluid is heated and soluted in such a way that a constant temperature and concentration distribution is prescribed at the boundaries of the fluid layer. The temperature (T) and concentration (C) are taken to be T0 and C0 at z=0 and T1 and C1 at z=d, (T0>T1,C0>C1). Let ΔT and ΔC be the difference in temperature and concentration across the boundaries.Letq(u,v,w), p,ρ,T,C,α,α′,μ, κ, and κ′ be the Darcy velocity vector, hydrostatic pressure, density, temperature, solute concentration, coefficient of thermal expansion, an analogous solvent coefficient of expansion, viscosity, thermal diffusivity, and solute diffusivity of fluid, respectively.


### 2.1. Assumptions

The mathematical equations describing the physical model are based upon the following assumptions.(i)
Thermophysical properties expect for density in the buoyancy force (Boussinesq hypothesis) are constant.(ii)
Darcy’s model with time derivative is employed for the momentum equation.(iii)
The porous medium is assumed to be isotropic and homogeneous.(iv)
No chemical reaction takes place in a layer of fluid.(v)
The fluid and solid matrix are in thermal equilibrium state.(vi)
Radiation heat transfer between the sides of the wall is negligible when compared with other modes of the heat transfer.


### 2.2. Governing Equations

According to the works of Bhatia and Steiner [14, 15], Sharma and Kumar [16], and Chand [19–21] the appropriate governing equations for Maxwell viscoelastic fluid in a porous medium are
(1)∇·q=0,0=(1+λ∂∂t)v×(-∇p+ρ(1-α(T-T0)-α′(C-C0))g)-μk1q,σ∂T∂t+q·∇T=κ∇2T+DTC∇2C,ε∂C∂t+q·∇C=κ′∇2C+DCT∇2T,
where DTC and DCT are the Dufour and Soret coefficients; σ=(ρcp)m/(ρcp)f is thermal capacity ratio, cp is specific heat, and the subscripts m and f refer to porous medium and fluid, respectively.We assume that temperature and concentration are constant at the boundaries of the fluid layer. Therefore, boundary conditions are(2)w=0,T=T0,C=C0atz=0,w=0,T=T1,C=C1atz=d.


### 2.3. Steady State and Its Solutions

The steady state is given by(3)u=v=w=0,p=p(z),T=Ts(z),C=Cs(z).
The solution of steady state is given as
(4)Ts=T0-ΔTdz,Cs=C0-ΔCdz,ps=p0-ρ0g(z+αΔT2dz2+α′ΔC2dz2),
where subscript 0 shows the value of the variable at boundary z=0.


### 2.4. Perturbation Solution

To study the stability of the system, we superimposed infinitesimal perturbations on the basic state, which are of the forms(5)q=0+q′,T=Ts+T′,C=Cs+C′,vvvvvvvvvvvvvvp=ps+p′,
where the prime denotes the perturbed quantities. Substituting (5) into (1) and neglecting higher order terms of the perturbed quantities, we get
(6)∇·q′=0,0=(1+λ∂∂t)(-∇p′+ρ0(αT′+α′C′)g)-μk1q′,σ∂T′∂t-w′ΔTd=κ∇2T′+DTC∇2C′,ε∂C′∂t-w′ΔCd=κ′∇2C′+DCT∇2T′.
Introduce the dimensionless variables as
(7)(x′′,y′′,z′′)=(x′,y′,z′d),(u′′,v′′,w′′,)=(u′,v′,w′κ)d,t′′=κσd2t,p′′=k1d2μκp′,T′′=T′ΔT,C′′=C′ΔC.
Thereafter drop the dashes (′′) for simplicity.Equation (6) in nondimensional form can be written as
(8)∇.q=0,0=(1+F∂∂t)(-∇p+RaT+RsC)-q,∂T∂t-w=∇2T+Df∇2C,εσ∂C∂t-w=1Le∇2C+Sr∇2T.
Here nondimensional parameters are as follows.Ra
=
g
ρ
0
α
k
1
Δ
T
d
/
μ
κ is the thermal Rayleigh number, Rs=gρ0α′k1ΔCd/μκ′ is the solutal Rayleigh number, Le=κ/κ′ is the Lewis number, F=(κ/σd2)λ is the stress relaxation parameter, Df=DTCΔC/κΔT is the Dufour parameter, and Sr=DCTΔT/κΔC is the Soret parameter.The nondimensional boundary conditions are(9)w=T=C=0atz=0,z=1.


## 2.1. Assumptions

The mathematical equations describing the physical model are based upon the following assumptions.(i)
Thermophysical properties expect for density in the buoyancy force (Boussinesq hypothesis) are constant.(ii)
Darcy’s model with time derivative is employed for the momentum equation.(iii)
The porous medium is assumed to be isotropic and homogeneous.(iv)
No chemical reaction takes place in a layer of fluid.(v)
The fluid and solid matrix are in thermal equilibrium state.(vi)
Radiation heat transfer between the sides of the wall is negligible when compared with other modes of the heat transfer.


## 2.2. Governing Equations

According to the works of Bhatia and Steiner [14, 15], Sharma and Kumar [16], and Chand [19–21] the appropriate governing equations for Maxwell viscoelastic fluid in a porous medium are
(1)∇·q=0,0=(1+λ∂∂t)v×(-∇p+ρ(1-α(T-T0)-α′(C-C0))g)-μk1q,σ∂T∂t+q·∇T=κ∇2T+DTC∇2C,ε∂C∂t+q·∇C=κ′∇2C+DCT∇2T,
where DTC and DCT are the Dufour and Soret coefficients; σ=(ρcp)m/(ρcp)f is thermal capacity ratio, cp is specific heat, and the subscripts m and f refer to porous medium and fluid, respectively.We assume that temperature and concentration are constant at the boundaries of the fluid layer. Therefore, boundary conditions are(2)w=0,T=T0,C=C0atz=0,w=0,T=T1,C=C1atz=d.


## 2.3. Steady State and Its Solutions

The steady state is given by(3)u=v=w=0,p=p(z),T=Ts(z),C=Cs(z).
The solution of steady state is given as
(4)Ts=T0-ΔTdz,Cs=C0-ΔCdz,ps=p0-ρ0g(z+αΔT2dz2+α′ΔC2dz2),
where subscript 0 shows the value of the variable at boundary z=0.


## 2.4. Perturbation Solution

To study the stability of the system, we superimposed infinitesimal perturbations on the basic state, which are of the forms(5)q=0+q′,T=Ts+T′,C=Cs+C′,vvvvvvvvvvvvvvp=ps+p′,
where the prime denotes the perturbed quantities. Substituting (5) into (1) and neglecting higher order terms of the perturbed quantities, we get
(6)∇·q′=0,0=(1+λ∂∂t)(-∇p′+ρ0(αT′+α′C′)g)-μk1q′,σ∂T′∂t-w′ΔTd=κ∇2T′+DTC∇2C′,ε∂C′∂t-w′ΔCd=κ′∇2C′+DCT∇2T′.
Introduce the dimensionless variables as
(7)(x′′,y′′,z′′)=(x′,y′,z′d),(u′′,v′′,w′′,)=(u′,v′,w′κ)d,t′′=κσd2t,p′′=k1d2μκp′,T′′=T′ΔT,C′′=C′ΔC.
Thereafter drop the dashes (′′) for simplicity.Equation (6) in nondimensional form can be written as
(8)∇.q=0,0=(1+F∂∂t)(-∇p+RaT+RsC)-q,∂T∂t-w=∇2T+Df∇2C,εσ∂C∂t-w=1Le∇2C+Sr∇2T.
Here nondimensional parameters are as follows.Ra
=
g
ρ
0
α
k
1
Δ
T
d
/
μ
κ is the thermal Rayleigh number, Rs=gρ0α′k1ΔCd/μκ′ is the solutal Rayleigh number, Le=κ/κ′ is the Lewis number, F=(κ/σd2)λ is the stress relaxation parameter, Df=DTCΔC/κΔT is the Dufour parameter, and Sr=DCTΔT/κΔC is the Soret parameter.The nondimensional boundary conditions are(9)w=T=C=0atz=0,z=1.


## 3.  Normal Modes and Stability Analysis

Analyze the disturbances into the normal modes and assume that the perturbed quantities are of the form(10)[w,T,C]=[W(z),Θ(z),Γ(z)]exp⁡(ikxx+ikyy+nt),
where kx,ky are wave numbers along x and y directions, respectively, and n is growth rate of disturbances.Using (10), (8) becomes
(11)(D2-a2)W+(1+Fn)(a2RaΘ+a2RsΓ)=0,W+(D2-a2-n)Θ+Df(D2-a2)Γ=0,W+Sr(D2-a2)Θ+(1Le(D2-a2)-εσn)Γ=0,
where D=d/dz and a=kx2+ky2 is the dimensionless resultant wave number.The boundary conditions are(12)W=0,D2W=0,Θ=0,Γ=0atz=0,W=0,D2W=0,Θ=0,Γ=0atz=1.
We assume the solution to W, Θ, and Γ is of the form
(13)W=W0sinπz,Θ=Θ0sin⁡πz,Γ=Γ0sin⁡πz,
which satisfy boundary conditions (12).Substituting solution (13) in (11), integrating each equation from z=0 to z=1 by parts, we obtain the following matrix equation as
(14)[J-a2(1+Fn)Ra-a2(1+Fn)Rs-1(J+n)DfJ-1SrJ(JLe+εnσ)][W0Θ0Γ0]=[000],
where J=π2+a2.The nontrivial solution of the above matrix requires that(15)Ra=(J(J+n)((J/Le)+(εn/σ))-SrDfJ2)a2(1+Fn)J((1/Le)-Df)+(εn/σ)+SrJ-(J+n)J((1/Le)-Df)+(εn/σ)Rs.
For neutral instability n=iω, (where ω is real and dimensionless frequency of oscillation) and equating real and imaginary parts of (15), we have
(16)J((J2Le-ω2εσ)-SrDfJ2)+a2Ra(J(Df-1Le)+ω2εFσ)-a2Rs(J-SrJ-ω2F)=0,J2(1Le+εσ)+a2Ra(FJ(Df-1Le)-εσ)-a2Rs(JF-SrJF+1)=0.
For stationary convection ω = 0 (n=0), we have
(17)Ra=(π2+a2)2a2(DfSrLe-1DfLe-1)+(1-Sr)LeDfLe-1Rs.
It is clear from (17) that stationary Rayleigh number Ra is a function of dimensionless wave number a, Dufour parameter Df, Soret parameter Sr, Lewis number Le and solutal Rayleigh number Rs, and independent of stress relaxation parameter F. Thus for stationary convection the Maxwell viscoelastic fluid behaves like an ordinary Newtonian fluid. This result is the same as obtained by Motsa [24] and Chand and Rana [23].The critical cell size at the onset of instability is obtained from the condition(∂Ra/∂a)a=ac=0, which gives ac=π.This result is the same as obtained by Lapwood [4] for Newtonian fluid.The corresponding critical Rayleigh numberRac for steady onset is
(18)Rac=4π2(DfSrLe-1DfLe-1)+(Sr-1)Le1-DfLeRs.
This result is the same as obtained by Motsa [24] and Chand and Rana [23].IfRs=Df=Sr=0 then
(19)Rac=4π2.
This is exactly the same result which was obtained by Nield [25].


## 4. Result and Discussion

The onset of double diffusive convection in a horizontal layer of Maxwell viscoelastic fluid in the presence of Soret and Dufour in a porous medium is investigated analytically and graphically. The expressions for both the stationary and oscillatory Rayleigh numbers, which characterize the stability of the system, are obtained analytically. The stationary critical Rayleigh number is found to be independent of the viscoelastic parameterF; thus Maxwell viscoelastic binary fluid behaves like ordinary Newtonian binary fluid. The stationary critical Rayleigh number and critical wave number are independent of viscoelastic parameter because of the absence of base flow in the present case. The computations are carried out for different values solutal Rayleigh number Rs, Soret parameter Sr, Dufour parameter Df, and Lewis number Le. The parameters considered are in the range of 102≤Ra≤105 (thermal Rayleigh number), 102≤Rs≤103 (solutal Rayleigh number), 0≤Sr≤1 (Soret parameter), 0≤Df≤1 (Dufour parameter), [24, 26], and 10-2≤Le≤1 (Lewis number) [27, 28], [29]. We choose the values of Dufour parameter Df and Lewis number Le in such a way that DfLe≠1.In order to investigate effects of the Dufour parameterDf, Soret parameter Sr, Lewis number Le, and solutal Rayleigh number Rs on stationary convection, we examine the behavior of ∂Ra/∂Df,∂Ra/∂Sr,∂Ra/∂Le,and∂Ra/∂Rs analytically.From (17), we have(i)
(20)∂Ra∂Df=(1-Sr)Le(DfLe-1)2[(π2+a2)2a2Df-LeRs],∂Ra∂Df>0ifLeRs>(π2+a2)2a2Df,<0ifLeRs<(π2+a2)2a2Df,=0 at two specific values of wave number, which depend upon Lewis number Le, solutal Rayleigh number Rs, and Dufour parameter Dfand is independent of Soret parameter Sr.Thus for the stationary convection Dufour parameterDf has a stabilizing effect if LeRs>((π2+a2)2/a2)Df and destabilizing effect if LeRs<((π2+a2)2/a2)Df:(ii)
(21)∂Ra∂Sr=LeDfLe-1[(π2+a2)2a2Df-Rs],∂Ra∂Sr>0if(π2+a2)2a2Df<Rs,<0if(π2+a2)2a2Df>Rs,=0 at two specific values of wave number, which depend upon solutal Rayleigh number Rs and Dufour parameter Df and is independent of Lewis number Le and Soret parameter Sr.Thus for stationary convection Soret parameterSr, therefore, has a stabilizing effect if ((π2+a2)2/a2)Df<Rs and destabilizing effect if ((π2+a2)2/a2)Df>Rs:(iii)
(22)∂Ra∂Le=(1-Sr)(DfLe-1)[(π2+a2)2a2Df-Rs],∂Ra∂Le>0if(π2+a2)2a2Df<Rs,<0if(π2+a2)2a2Df>Rs,=0 at two specific values of wave number, which depend upon solutal Rayleigh number Rs and Dufour parameter Df and is independent of Lewis number Le and Soret parameter Sr.Thus for stationary convection Lewis number Le has a stabilizing effect if((π2+a2)2/a2)Df<Rs and destabilizing effect if ((π2+a2)2/a2)Df>Rs:(iv)
(23)∂Ra∂Rs=(1-Sr)LeDfLe-1<0.Thus for stationary convection solutal Rayleigh number Rs has a destabilizing effect.Now we discussed the effects of various parameters on the onset of double diffusive convection of Maxwell viscoelastic fluid in a porous medium for stationary convection graphically. The convection curves for solutal Rayleigh number Rs, Soret parameterSr, Dufour parameter Df, and Lewis number Le in the (Ra, a) plane are shown in Figures 1–4.Figure 1
Variation of Rayleigh number Ra with wave numbera for different values of Dufour parameter Df.Figure 2
Variation of Rayleigh numberRa with wave number a for different values of Soret parameter Sr.Figure 3
Variation of Rayleigh number Ra with wave numbera for different values of Lewis number Le.Figure 4
Variation of Rayleigh number Ra with wave numbera for different values of solutal Rayleigh number Rs.Figure1 shows the variation of stationary Rayleigh number with wave number for different values of Dufour parameter and it is found that the Rayleigh number first increases then decreases and finally increases with increase in the value of Dufour parameter; thus for stationary convection Dufour parameter has both the stabilizing and destabilizing effects depending upon certain conditions.Figure2 shows the variation of stationary Rayleigh number with wave number for different values of Soret parameter and it is found that the Rayleigh number first decreases then increases and finally decreases with increase in the value of Soret parameter; thus for stationary convection Soret parameter has both the stabilizing and destabilizing effects depending upon certain conditions.Figure3 shows the variation of stationary Rayleigh number with wave number for different values of Lewis number and it is found that the Rayleigh number first increases then decreases and finally increases with increase in the value of Lewis number; thus for stationary convection Lewis number has both the stabilizing and destabilizing effects depending upon certain conditions.Figure4 shows the variation of Rayleigh number with wave number for different value of the solutal Rayleigh number Rs and it is found that the Rayleigh number decreases with increase in the value of solutal Rayleigh number Rs; thus solutal Rayleigh number Rs has destabilizing effect on the stationary convection.Curves in Figures1–3 intersect two specific values of wave number a which is independent of Lewis number Le and Soret parameter Sr as explained in the analytical part of the results and discussion.


## 5. Conclusions

A linear stability analysis of double diffusive convection in a horizontal layer of Maxwell viscoelastic fluid in the presence of Soret and Dufour in a porous medium is investigated analytically and graphically. The expressions for both the stationary and oscillatory Rayleigh numbers, which characterize the stability of the system, are obtained.The main conclusions are as follows.(i)
In stationary convection Maxwell viscoelastic fluid behaves like ordinary Newtonian fluid.(ii)
Dufour parameter, Soret parameter, and Lewis parameter have both stabilizing and destabilizing effects on the stationary convection.(iii)
Solutal Rayleigh number destabilizes the stationary convection.(iv)
In limiting case when Rs =Sr = Df = 0 the critical thermal Rayleigh number obtained is the same as reported by Nield [25].



---

*Source: 479107-2014-07-17.xml*