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

This chapter presents some new exact solutions corresponding to unsteady magnetohy- drodynamic (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.


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.

Governing equations
The equation of continuity and momentum of MHD flow passing through porous space is given by [7] 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 ¼ μ where S and pI represents the extra stress tensor and the indeterminate spherical stress, the dynamic viscosity is denoted by μ, A = L + L T 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 where σ represents electrical conductivity and β o the strength of magnetic field. For the Jeffrey fluid the Darcy's resistance satisfies the following equation 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 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 therefore from Eqs. (2), (5) and (6), it results that S xx = S yy = S yz = S zz = 0 and the relevant equations where τ 1 = S xy and τ 2 = S xz are the tangential stresses. In the absence of pressure gradient in the flow direction, the governing equation leads to ρ is the magnetic parameter, K ¼ ϕ κ is the porosity parameter and ν = μ/ρ is the kinematic viscosity.

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 t > 0, 0 < x < d and 0 < 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 which is the solution of the problem t > 0, 0 < x < d and 0 < 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

Calculation of the velocity field
Multiplying both sides of Eq. (13) where The Fourier transform F mn (t) have to satisfy the initial conditions We apply LT to Eq. (17) and using initial conditions (18) to get 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 We apply the discrete inverse LT to Eq. (20), to obtain Taking the inverse Fourier sine transform we get the analytic solution of the velocity field To obtain a more compact form of velocity field we write Eq. (22) in terms of Fox H-function, or where From Eq. (24), we obtain the velocity field due to cosine oscillations of the duct u x; y; t ð Þ¼ 16U cos wt ð Þ À wθ sin wt ð Þ ð Þ dh and the velocity field due to sine oscillations of the duct v x; y; t We use the following property of the Fox H-function [21] in the above equation

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 τ 2 x; y; t ð Þ¼τ 2c x; y; t ð Þþiτ 2s x; y; t ð Þ, in Eq. (7), we get We apply LT to Eqs. (29) and (30), to obtain Taking the inverse Fourier transform of Eq. (19) to get F x; y; s ð Þ and then by putting it into Eq. (31), we get or where We express Eq. (34) in series form in order to obtain a more suitable form of τ 1 , Using the inverse LT of the last equation, we obtain Lastly, we write the stress field in a more compact form by using Fox H-function From Eq. (37), we obtain the tangential tension due to cosine oscillations of the duct and the tangential tension corresponding to sine oscillations of the duct In the similar fashion we can find τ 2c (x, y, t) and τ 2s (x, y, t) from Eqs. (19) and (32).

Volume flux
The volume flux due to cosine oscillations is given by 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 Similarly, we obtain the volume flux of the rectangular duct due to the sine oscillations v x; y; t ð Þ ¼ 64U sin wt ð Þ À wθ cos wt ð Þ ð Þ dh

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.