Open access peer-reviewed chapter

Unsteady Magnetohydrodynamic Flow of Jeffrey Fluid through a Porous Oscillating Rectangular Duct

By Amir Khan, Gul Zaman and Obaid Algahtani

Submitted: March 20th 2017Reviewed: September 8th 2017Published: April 26th 2018

DOI: 10.5772/intechopen.70891

Downloaded: 302

Abstract

This chapter presents some new exact solutions corresponding to unsteady magnetohydrodynamic (MHD) flow of Jeffrey fluid in a long porous rectangular duct oscillating parallel to its length. The exact solutions are established by means of the double finite Fourier sine transform (DFFST) and discrete Laplace transform (LT). The series solution of velocity field, associated shear stress and volume flow rate in terms of Fox H-functions, satisfying all imposed initial and boundary conditions, have been obtained. Also, the obtained results are analyzed graphically through various pertinent parameter.

Keywords

  • porous medium
  • Jeffrey fluid
  • oscillating rectangular duct
  • Fox H-function
  • MSC (2010): 76A05
  • 76A10

1. Introduction

Considerable progress has been made in studying flows of non-Newtonian fluids throughout the last few decades. Due to their viscoelastic nature non-Newtonian fluids, such as oils, paints, ketchup, liquid polymers and asphalt exhibit some remarkable phenomena. Amplifying interest of many researchers has shown that these flows are imperative in industry, manufacturing of food and paper, polymer processing and technology. Dissimilar to the Newtonian fluid, the flows of non-Newtonian fluids cannot be explained by a single constitutive model. In general the rheological properties of fluids are specified by their so-called constitutive equations. Exact recent solutions for constitutive equations of viscoelastic fluids are given by Rajagopal and Bhatnagar [1], Tan and Masuoka [2, 3], Khadrawi et al. [4] and Chen et al. [5] etc. Among non-Newtonian fluids the Jeffrey model is considered to be one of the simplest type of model which best explain the rheological effects of viscoelastic fluids. The Jeffrey model is a relatively simple linear model using the time derivatives instead of convected derivatives. Nadeem et al. [6] obtained analytic solutions for stagnation flow of Jeffrey fluid over a shrinking sheet. Khan [7] investigated partial slip effects on the oscillatory flows of fractional Jeffrey fluid in a porous medium. Hayat et al. [8] examined oscillatory rotating flows of a fractional Jeffrey fluid filling a porous medium. Khan et al. [9] discussed unsteady flows of Jeffrey fluid between two side walls over a plane wall.

Much attention has been given to the flows of rectangular duct because of its wide range applications in industries. Gardner and Gardner [10] discussed magnetohydrodynamic (MHD) duct flow of two-dimensional bi-cubic B-spline finite element. Fetecau and Fetecau [11] investigated the flows of Oldroyd-B fluid in a channel of rectangular cross-section. Nazar et al. [12] examined oscillating flow passing through rectangular duct for Maxwell fluid using integral transforms. Unsteady magnetohydrodynamic flow of Maxwell fluid passing through porous rectangular duct was studied by Sultan et al. [13]. Tsangaris and Vlachakis [14] discussed analytic solution of oscillating flow in a duct of Navier-Stokes equations.

In the last few decades the study of fluid motions through porous medium have received much attention due to its importance not only to the field of academic but also to the industry. Such motions have many applications in many industrial and biological processes such as food industry, irrigation problems, oil exploitation, motion of blood in the cardiovascular system, chemistry and bio-engineering, soap and cellulose solutions and in biophysical sciences where the human lungs are considered as a porous layer. Unsteady MHD flows of viscoelastic fluids passing through porous space are of considerable interest. In the last few years a lot of work has been done on MHD flow, see [15, 16, 17, 18, 19] and reference therein.

According to the authors information up to yet no study has been done on the MHD flow of Jeffrey fluid passing through a long porous rectangular duct oscillating parallel to its length. Hence, our main objective in this note is to make a contribution in this regard. The obtained solutions, expressed under series form in terms of Fox H-functions, are established by means of double finite Fourier sine transform (DFFST) and Laplace transform (LT). Finally, the obtained results are analyzed graphically through various pertinent parameter.

2. Governing equations

The equation of continuity and momentum of MHD flow passing through porous space is given by [7]

·V=0,ρdVdt=divT+J×B+R,E1

where velocity is represented by V, density by ρ, Cauchy stress tensor by T, magnetic body force by J × B, current density by J, magnetic field by B, and Darcy’s resistance in the porous medium by R.

For an incompressible and unsteady Jeffrey fluid the Cauchy stress tensor is defined as [9]

T=pI+S,S=μ1+λA+θAt+V.A,E2

where S and pI represents the extra stress tensor and the indeterminate spherical stress, the dynamic viscosity is denoted by μ, A = L + LT is the first Rivlin-Ericksen tensor, L is the velocity gradient, λ is relaxation time and θ is retardation time. The Lorentz force due to magnetic field is

J×B=σβo2V,E3

where σ represents electrical conductivity and βo the strength of magnetic field. For the Jeffrey fluid the Darcy’s resistance satisfies the following equation

R=μϕκ1+λ1+θtV,E4

where κ(>0) and ϕ(0 < ϕ < 1) are the permeability and the porosity of the porous medium.

In the following problem we consider a velocity field and extra stress of the form

V=00wxyt,S=SxytE5

where w is the velocity in the z-direction. The continuity equation for such flows is automatically satisfied. Also, at t = 0, the fluid being at rest is given by

Sxy0=0,E6

therefore from Eqs. (2), (5) and (6), it results that Sxx = Syy = Syz = Szz = 0 and the relevant equations

τ1=μ1+λ1+θtxwxyt,τ2=μ1+λ1+θtywxyt,E7

where τ1 = Sxy and τ2 = Sxz are the tangential stresses. In the absence of pressure gradient in the flow direction, the governing equation leads to

1+λtwxyt=ν1+θtx2+y2wxytνK1+θtwxytH1+λwxyt,E8

where H=σB02ρis the magnetic parameter, K=ϕκis the porosity parameter and ν = μ/ρ is the kinematic viscosity.

3. Statement of the problem

We take an incompressible flow of Jeffrey fluid in a porous rectangular duct under an imposed transverse magnetic field whose sides are at x = 0, x = d, y = 0, and y = h. At time t = 0+ the duct begins to oscillate along z-axis. Its velocity is of the form of Eq. (5) and the governing equation is given by Eq. (8). The associated initial and boundary conditions are

wxy0=twxy0=0,E9
w0yt=wx0t=wdyt=wxht=Ucoswt,E10

or

w0yt=wx0t=wdyt=wxht=Usinwt,E11
t>0,0<x<dand0<y<h.

The solutions of problems (8)(10) and (8), (9), (11) are denoted by u(x, y, t) and v(x, y, t) respectively. We define the complex velocity field

Fxyt=uxyt+ivxyt,E12

which is the solution of the problem

1+λtFxyt=ν1+θtx2+y2FxytνK1+θtFxytH1+λFxyt,E13
Fxy0=tFxy0=0,E14
F0yt=Fdyt=Fx0t=Fxht=Ueiwt,E15
t>0,0<x<dand0<y<h.

The solution of the problem (13)(15) will be obtained by means of the DFFST and LT.

The DFFST of function F(x, y, t) is denoted by

Fmnt=0d0hsinmπxdsinnπyhFxytdxdy,m,n=1,2,3,..E16

4. Calculation of the velocity field

Multiplying both sides of Eq. (13) by sinmπxdand sinnπyh, integrating w.r.t x and y over [0, d] × [0, h] and using Eq. (16), we get

(1+λ)Fmn(t)t+νλmn(1+θt)Fmn(t)+H(1+λ)Fmn(t)+νK(1+θt)Fmn(t)=νλmnU[1(1)m][1(1)n]ζmλn(1+iwθ)eiwt,E17

where

ζm=d,λn=handλmn=ζm2+λn2.

The Fourier transform Fmn(t) have to satisfy the initial conditions

Fmn0=tFmn0=0.E18

We apply LT to Eq. (17) and using initial conditions (18) to get

F¯mns=νλmnU11m1+iwθ11nζmλnsiw1+λs+H+ν1+θsλmn+K.E19

We will apply the discrete inverse LT technique [20] to obtain analytic solution for the velocity fields and to avoid difficult calculations of residues and contour integrals, but first we express Eq. (19) in series form as

F¯mns=νλmnU11m1+iwθ11nζmλnsiwp=0q=0r=0s=0l=0×νp+1λsθqKprHlλmnr+1ΓqpΓrpΓs+p+1Γl+p+11p+q+r+s+lq!r!s!l!ΓpΓpΓ1+pΓ1+pslq+p+1.E20

We apply the discrete inverse LT to Eq. (20), to obtain

Fmnt=eiwtU11m11nνλmn1+iwθζmλnp=0q=0r=0s=0l=0×νp+1λsθqKprHlλmnr+1ΓqpΓrpΓs+p+1Γl+p+1tlq+p1p+q+r+s+lq!r!s!l!ΓpΓpΓ1+pΓ1+pΓlq+p+1.E21

Taking the inverse Fourier sine transform we get the analytic solution of the velocity field

Fxyt=4dhm=1n=1sinζmxsinλnyFmnxyt=4eiwtU1+iwθdhm=1n=111m11nsinζxsinλnyζmλn×p=0q=0r=0s=0l=0νp+1λsθqKprHlλmnr+1tlq+p1p+q+r+s+lq!r!s!l!×ΓqpΓrpΓs+p+1Γl+p+1ΓpΓpΓ1+pΓ1+pΓlq+p+1.E22

To obtain a more compact form of velocity field we write Eq. (22) in terms of Fox H-function,

Fxyt=4eiwtU1+iwθdhm=1n=1sinζmxsinλny11m11nζmλn×p=0q=0r=0s=01p+q+r+sνp+1λsθqKprλmnr+1tq+pq!r!s!×H4,61,4[Ht1q+p0,1r+p0,sp0,p101,1p0,1p0,p0,p0,qp1].E23

or

Fxyt=16eiwtU1+iwθdhc=0e=0sinζcxsinλeyζcλe×p=0q=0r=0s=01p+q+r+sνp+1λsθqKprλcer+1tq+pq!r!s!×H4,61,4[Ht1q+p0,1r+p0,1sp0,p101,1p0,1p0,p0,p0,qp1]E24

where

ζc=2m+1πd,λe=2n+1πh,c=2m+1,e=2n+1.

From Eq. (24), we obtain the velocity field due to cosine oscillations of the duct

uxyt=16Ucoswtsinwtdhc=0e=0sinζcxsinλeyζcλe×p=0q=0r=0s=01p+q+r+sνp+1λsθqKprλcer+1tq+pq!r!s!×H4,61,4[Ht1q+p0,1r+p0,1sp0,p101,1p0,1p0,p0,p0,qp1]E25

and the velocity field due to sine oscillations of the duct

vxyt=16Usinwtcoswtdhc=0e=0sinζcxsinλeyζcλe×p=0q=0r=0s=01p+q+r+sνp+1λsθqKprλcer+1tq+pq!r!s!×H4,61,4[Ht1q+p0,1r+p0,1sp0,p101,1p0,1p0,p0,p0,qp1].E26

We use the following property of the Fox H-function [21] in the above equation

Hp,q+11,p[χ1a1A1,1a2A2,,1apAp10,1b1B1,,1bqBq]=k=0Γa1+A1kΓap+Apkk!Γb1+B1kΓbq+Bqkχk.

5. Calculation of the shear stress

We denote the tangential tensions for the cosine oscillations of the duct by τ1c(x, y, t), τ2c(x, y, t) and for sine oscillations by τ1s(x, y, t), τ2s(x, y, t).

If we introduce

τ1xyt=τ1cxyt+iτ1sxyt,E27
τ2xyt=τ2cxyt+iτ2sxyt,E28

in Eq. (7), we get

τ1xyt=μ1+λ1+θtxFxyt,E29
τ2xyt=μ1+λ1+θtyFxyt.E30

We apply LT to Eqs. (29) and (30), to obtain

τ¯1xys=μ1+θs1+λxF¯xys,E31
τ¯2xys=μ1+θs1+λyF¯xys.E32

Taking the inverse Fourier transform of Eq. (19) to get F¯xysand then by putting it into Eq. (31), we get

τ¯1xys=4μ1+θsdh1+λm=1n=1cosζmxsinλny11m11n1+λs+H+ν1+θsλmn+K×λmn1+iwθλnsiw,E33

or

τ¯1xys=16μ1+θsdh1+λc=0e=0cosζcxsinλeyλce1+iwθλesiw1+λs+H+ν1+θsλce+K,E34

where

ζc=2m+1πd,λe=2n+1πh,c=2m+1,e=2n+1.

We express Eq. (34) in series form in order to obtain a more suitable form of τ1,

τ¯1xys=16μU1+iwθdhsiwc=0e=0cosζcxsinλeyλep=0q=0r=0s=0l=0×νp+1λsθqKprHlλcer+1Γqp1ΓrpΓs+p+2Γl+p+11p+q+r+s+lq!r!s!l!ΓpΓp+1Γ2+pΓ1+pslq+p+1.E35

Using the inverse LT of the last equation, we obtain

τ1xyt=16μeiwtU1+iwθdhc=0e=0cosζcxsinλeyλep=0q=0r=0s=0×l=0νp+1λsθqKprHlλcer+1tlq+pΓqp1Γs+p+21p+q+r+s+lq!r!s!l!ΓpΓp+1Γ2+pΓ1+p×Γl+p+1Γlq+p+1.E36

Lastly, we write the stress field in a more compact form by using Fox H-function

τ1xyt=16μeiwtU1+iwθdhc=0e=0cosζcxsinλeyλep=0×q=0r=0s=0νp+1λsθqKprλmnr+1tq+p1p+q+r+sq!r!s!×H4,61,4[Ht2q+p0,1r+p0,1sp0,p101,p0,1p0,p0,p0,qp1].E37

From Eq. (37), we obtain the tangential tension due to cosine oscillations of the duct

τ1cxyt=16coswtsinwtdhc=0e=0cosζcxsinλeyλe×p=0q=0r=0s=01p+q+r+sνp+1λsθqKprλmnr+1tq+pq!r!s!×H4,61,4[Ht2q+p0,1r+p0,1sp0,p101,p0,1p0,p0,p0,qp1]E38

and the tangential tension corresponding to sine oscillations of the duct

τ1sxyt=16sinwtcoswtdhc=0e=0cosζcxsinλeyλep=0×q=0r=0s=01p+q+r+sνp+1λsθqKprλmnr+1tq+pq!r!s!×H4,61,4[Ht2q+p0,1r+p0,1sp0,p101,p0,1p0,p0,p0,qp1].E39

In the similar fashion we can find τ2c(x, y, t) and τ2s(x, y, t) from Eqs. (19) and (32).

6. Volume flux

The volume flux due to cosine oscillations is given by

Qcxyt=0d0huxytdxdy,E40

putting u(x, y, t) from Eq. (25) into the above equation, we obtain the volume flux of the rectangular duct due to cosine oscillations

uxyt=64Ucoswtsinwtdhc=0e=01ζcλe2×p=0q=0r=0s=0νp+1λsθqKprλcer+1tq+p1p+q+r+sq!r!s!×H4,61,4[Ht1q+p0,1r+p0,1sp0,p101,1p0,1p0,p0,p0,qp1].E41

Similarly, we obtain the volume flux of the rectangular duct due to the sine oscillations

vxyt=64Usinwtcoswtdhc=0e=01ζcλe2×p=0q=0r=0s=0νp+1λsθqKprλcer+1tq+p1p+q+r+sq!r!s!×H4,61,4[Ht1q+p0,1r+p0,1sp0,p101,1p0,1p0,p0,p0,qp1].E42

7. Numerical results and discussion

We have presented flow problem of MHD Jeffrey fluid passing through a porous rectangular duct. Exact analytical solutions are established for such flow problem using DFFST and LT technique. The obtained solutions are expressed in series form using Fox H-functions. Several graphs are presented here for the analysis of some important physical aspects of the obtained solutions. The numerical results show the profiles of velocity and the adequate shear stress for the flow. We analyze these results by variating different parameters of interest.

The effects of relaxation time λ of the model are important for us to be discuss. In Figure 1 we depict the profiles of velocity and shear stress for three different values of λ. It is observed from these figures that the flow velocity as well as the shear stress decreases with increasing λ, which corresponds to the shear thickening phenomenon. Figure 2 are sketched to show the velocity and the shear stress profiles at different values of retardation time θ. It is noticeable that velocity as well as the shear stress decreases by increasing θ. In order to study the effect of frequency of oscillation ω, we have plotted Figure 3, where it appears that the velocity is also a strong function of ω of the Jeffrey fluid. The effect of frequency of oscillation on the velocity profile for cosine oscillation is same as that of the retardation time θ. The effect of magnetic parameter H of the model is important for us to be discussed. In Figure 4, we depict the profiles of velocity and shear stress for three different values of H. It is observed from these figures that the flow velocity as well as the shear stress decreases with increasing H, which corresponds to the shear thickening phenomenon. Figure 5 is sketched to show the velocity and the shear stress profiles at different values of K. It is noticeable that velocity as well as the shear stress increases by increasing K. In order to study the effects of t, we have plotted Figure 6, where it appears that the velocity is also a strong function of t of the Jeffrey fluid. It can be observed that the increase of t acts as an increase of the magnitude of velocity components near the plate, and this corresponds to the shear-thinning behavior of the examined non-Newtonian fluid. Figure 7 presents the velocity field and the shear stress profiles at different values of y. It is noticeable that velocity and shear stress both decreases by increasing y.

Figure 1.

Velocity and shear stress profiles corresponding to the cosine oscillations of the duct for different values of λ. Other parameters are taken as x = 0.5, y = 0.3, U = 0.2, H = 0.5, K = 0.6, d = 1, h = 2, θ = 0.6, ω = 0.5 and ν = 0.1.

Figure 2.

Velocity and shear stress profiles corresponding to the cosine oscillations of the duct for different values of θ. Other parameters are taken as x = 0.5, y = 0.3, U = 0.2, H = 0.5, K = 0.6, d = 1, h = 2, λ = 1.4, ω = 0.5 and ν = 0.1.

Figure 3.

Velocity and shear stress profiles corresponding to the cosine oscillations of the duct for different values of ω. Other parameters are taken as x = 0.5, y = 0.3, U = 0.2, H = 0.5, K = 0.6, d = 1, h = 2, θ = 0.6, λ = 1.4 and ν = 0.1.

Figure 4.

Velocity and shear stress profiles corresponding to the cosine oscillations of the duct for different values of H. Other parameters are taken as x = 0.5, y = 0.3, U = 0.2, λ = 1.4, K = 0.6, d = 1, h = 2, θ = 0.6, ω = 0.5 and ν = 0.1.

Figure 5.

Velocity and shear stress profiles corresponding to the cosine oscillations of the duct for different values of K. Other parameters are taken as x = 0.5, y = 0.3, U = 0.2, H = 0.5, λ = 1.4, d = 1, h = 2, θ = 0.6, ω = 0.5 and ν = 0.1.

Figure 6.

Velocity and shear stress profiles corresponding to the cosine oscillations of the duct for different values of t. Other parameters are taken as x = 0.5, λ = 1.4, U = 0.2, H = 0.5, K = 0.6, d = 1, h = 2, θ = 0.6, ω = 0.5 and ν = 0.1.

Figure 7.

Velocity and shear stress profiles corresponding to the cosine oscillations of the duct for different values of y. Other parameters are taken as x = 0.5, λ = 1.4, U = 0.2, H = 0.5, K = 0.6, d = 1, h = 2, θ = 0.6, ω = 0.5 and ν = 0.1.

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Amir Khan, Gul Zaman and Obaid Algahtani (April 26th 2018). Unsteady Magnetohydrodynamic Flow of Jeffrey Fluid through a Porous Oscillating Rectangular Duct, Porosity - Process, Technologies and Applications, Taher Hcine Ghrib, IntechOpen, DOI: 10.5772/intechopen.70891. Available from:

chapter statistics

302total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Porous Structures in Heat Pipes

By Patrik Nemec

Related Book

First chapter

Study of CrAlN Multilayred Thin Films

By Tlili Ibrahim and Taher Ghrib

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.

More About Us