# Unsteady Hydromagnetic Rotating Flow through an Oscillating Porous Plate Embedded in a Porous Medium

Authors: I. Khan; A. Khan; A. Farhad; M. Qasim; S. Sharidan
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

This paper investigates unsteady hydromagnetic flow of a viscous fluid in a rotating frame. The fluid is bounded by an oscillating porous plate embedded in a porous medium. The Laplace transform and Fourier sine transform methods are employed to find the exact solutions. They satisfy all imposed initial and boundary conditions and as special cases are reduced to some published results from the literature. The graphical results are plotted for different values of pertinent parameters and some interesting conclusions are made.


---

## Body

## 1. Introduction

Stokes’ problem for the flow of an incompressible, viscous fluid over an oscillating plane serves as a benchmark in the literature of fluid dynamics [1]. It admits an exact analytical solution. The Stokes’ problem not only is of fundamental theoretical interest but also occurs in many applied problems [2]. The transient solutions in terms of tabulated functions for the flow of viscous fluid due to oscillating boundary have been expressed by Panton [3]. Later, in 2000, Erdogan [4] studied the unsteady flow of Newtonian fluid past an oscillating horizontal plane wall and obtained the exact solutions. Fetecau et al. [5] presented new exact solutions corresponding to Stokes’ second problem for both small and large times. Erdoǧan and Imrak [6] made comparative study for the solution of Stokes’ second problem by employing two different transform methods. They expressed transient solutions in terms of tabulated functions.On the other hand, flow through porous media is very prevalent in nature and, therefore, has become a principal interest of researchers in many scientific and engineering studies. For example, one can refer to the books of Pop and Ingham [7], Ingham et al. [8], Ingham and Pop [9], Vafai [10], and Nield and Bejan [11] for the detailed literature on this topic. Moreover, the rotating flow through porous media has been the subject of many studies in the last few decades due to their wide range of applications in cosmological and geophysical fluid dynamics, astrophysics, meteorology, petroleum, and hydrology to study the movement of underground water (Jana et al. [12]). According to Hayat et al. [13], the analysis of the effects of rotation and magnetic field in fluid flows has been an active area of research because of its geophysical and technological applications such as the MHD power generator and boundary layer flow control. Based on this motivation, Das et al. studied in [14] the simultaneous effects of rotation and magnetic field on the flow of a second grade fluid between two parallel plates. Bearing in mind the importance of MHD and porosity in a rotating fluid, Hayat et al. [15, 16], Abelman et al. [17], Sahoo et al. [18], Seth et al. [19], and Farhad et al. [20] studied the hydromagnetic flow of rotating fluids in porous media. Recently, Jana et al. [21] studied the unsteady flow of viscous fluid through a porous medium bounded by a porous plate in a rotating system.In order to further discuss the work of Jana et al. [21] and to make closer relations between this study and practical engineering, we study in this work hydromagnetic rotating flow of a viscous fluid bounded by an oscillating porous plate embedded in a porous medium. The governing mathematical problem is solved by using the Laplace transform and Fourier sine transform methods. The expressions for velocity and skin friction in case of cosine and sine oscillations of the plate are obtained. The analytical results are plotted and discussed. The results of Jana et al. [21] are recovered as a special case from our obtained solutions. The present study is useful in astronomy to study the rotating motion of astrophysical objects such as magnetic stars.


## 2. Mathematical Formulation of the Problem

Consider the unsteady flow of an incompressible viscous fluid occupying the upper porous half-space of the(x,y) plane.The fluid is bounded by a porous plate at z=0such that the positive z-axis is taken normal to the plate and the x-axis is taken parallel to the plate. We consider the hydromagnetic flow induced in the fluid in the presence of a uniform magnetic field of strength B0applied in a direction normal to the plate by means of the plate oscillations in its plane (Hayat et al. [13]). The electric field due to polarization of charges is neglected. Both of the fluid and plate are in a state of rigid body rotation with constant angular velocity Ω=Ωk^; k^ is a unit vector in the z-direction. The geometry of the problem is shown in Figure 1.Initially, for t≤0, both of the fluid and plate are at rest. At time t=0+, the lower plate suddenly starts to move in its own plane with oscillating velocity in the flow direction. Under these assumptions, the equations of momentum in a rotating frame along the x- and y-directions are [18–20]
(1)∂u∂t-w0∂u∂z-2Ωv=ν∂2u∂z2-νϕku-σB02ρu,z,t>0,∂v∂t-w0∂v∂z+2Ωu=ν∂2v∂z2-νϕkv-σB02ρv,z,t>0,
whereuand vare the velocity components in x- and y-directions,w0is the suction velocity at the plate, ν is the kinematic viscosity, t is the time, σ is the electrical conductivity, B0is the applied magnetic field, ρ is the fluid density, ϕ is the porosity, and k is the permeability of the porous medium.Figure 1
Geometry of the problem.The corresponding initial and boundary conditions are(2)u=v=0;t≤0,(3)u=UH(t)cos(ω0t)orUsin(ω0t);v=0;atz=0;t>0,(4)u⟶0;v⟶0asz⟶∞;t>0,
in which Uis the constant plate velocity,ω0is the frequency of oscillation of the plate, and H(t)is the Heaviside function.DefiningF=u+iv, the above system of equations reduces to
(5)∂F∂t-w0∂F∂z+2iΩF=ν∂2F∂z2-νϕkF-σB02ρF,z,t>0,F(z,0)=0;t≤0,F(0,t)=UH(t)cos(ω0t)orUsin(ω0t);t>0F⟶0asz⟶∞;t>0.
Introducing the following dimensionless variables (see Jana et al. [21]):
(6)G=FU,ξ=zUν,τ=U2tν,ω=ω0νU2,
the dimensionless problem takes the following form:
(7)∂G∂τ-S∂G∂ξ+2iηG=∂2G∂ξ2-M2G-1KG,(8)G(ξ,0)=0;τ≤0,(9)G(0,τ)=H(τ)cos(ωτ)orsin(ωτ);τ>0,(10)G(∞,τ)=0;τ>0,(11)M2=σB02νρU2,1K=ν2ϕkU2,η=ΩνU2,S=w0U,
where M,K,η, and Sare the magnetic, permeability, rotation, and suction parameters, respectively.


## 3. Solution of the Problem by Using Laplace Transform

Taking the Laplace transform of (7)  and using the initial condition (8), we get the following transformed differential equation:
(12)d2G¯(ξ,q)dξ2+SdG¯(ξ,q)dξ-[q+M2+1K+2iη]G¯(ξ,q)=0,
where qis the transform parameter. In view of the boundary conditions (9) and (10), the Laplace transform G¯(ξ,q)of G(ξ,τ)  has to satisfy the following conditions:
(13)G¯(0,q)=qq2+ω2orG¯(0,q)=ωq2+ω2,(14)G¯(0,q)=0.
Solution of (12) under the boundary conditions (13) and (14) yields
(15)G¯c(ξ,q)=12e-ξS/2{1q+iωe-ξq+a+1q-iωe-ξq+a},
respectively,
(16)G¯s(ξ,q)=12ie-ξS/2{-1q+iωe-ξq+a+1q-iωe-ξq+a},
where the subscripts cand sdenote the cosine and sine oscillations of the plate and
(17)a=S24+1K+M2+2iη.
Now, in order to find the inverse Laplace transform of (15) and (16), we use formula (A4) from the Appendix of Farhad et al. [20]; hence, we get the following solutions:
(18)Gc(ξ,τ)=H(τ)e-ξS/24×[e-iωτ{e-ξa-iωerfc(ξ2τ-(a-iω)τ)e-iωτ+eξa-iωerfc(ξ2τ+(a-iω)τ)}]+H(τ)e-ξS/24×[eiωτ{e-ξa+iωerfc(ξ2τ-(a+iω)τ)+e-ξa+iωerfc(ξ2τ+(a+iω)τ)}],
respectively,
(19)Gs(ξ,τ)=e-ξS/24i×[-e-iωτ{e-ξa-iωerfc(ξ2τ-(a-iω)τ)+eξa-iωerfc(ξ2τ+(a-iω)τ)}]+e-ξS/24i[eiωτ{e-ξa+iωerfc(ξ2τ-(a+iω)τ)+eξa+iωerfc(ξ2τ+(a+iω)τ)}].


## 4. Solution of the Problem by Using Fourier Sine Transform

Now, applying Fourier sine transform to (12) and using boundary conditions (13) and (14), we obtain
(20)G~c(ζ,q)=2(ζ+ζS)q(ζ2+ζ2S+q+a)(q2+ω2),
respectively,
(21)G~s(ζ,q)=2(ζ+ζS)ω(ζ2+ζ2S+q+a)(q2+ω2),
where G~(ζ,q)denotes Fourier sine transform of G¯(ξ,q).Now, we take the inverse Fourier sine transform of (20), which yields
(22)Gc(ξ,q)=22qπ(q2+ω2)∫0∞‍(ζ+ζS)(ζ2+ζ2S+q+a)sin(ζξ)dζ.Taking the inverse Laplace transform of (22), we obtain
(23)Gc(ξ,τ)=2π∫0∞‍ζ(ζ2+ζ2S+a)e-(ζ2+ζ2S+a)τ((ζ2+ζ2S+a)2+ω2)×sin(ζξ)dζ-2π∫0∞‍ζS(ζ2+ζ2S+a)e-(ζ2+ζ2S+a)τ((ζ2+ζ2S+a)2+ω2)×sin(ζξ)dζ+2π∫0∞‍ζ(ζ2+ζ2S+a)((ζ2+ζ2S+a)2+ω2)×cos(ωτ)sin(ζξ)dζ+2π∫0∞‍ζω((ζ2+ζ2S+a)2+ω2)×sin(ωτ)sin(ζξ)dζ+2π∫0∞‍ζS(ζ2+ζ2S+a)((ζ2+ζ2S+a)2+ω2)×cos(ωτ)sin(ζξ)dζ+2π∫0∞‍ζSω((ζ2+ζ2S+a)2+ω2)×sin(ωτ)sin(ζξ)dζ.
Similarly for the sine oscillation of the boundary we get from (21) the following:
(24)Gs(ξ,τ)=2π∫0∞‍ζωe-(ζ2+ζ2S+a)τ((ζ2+ζ2S+a)2+ω2)sin(ζξ)dζ-2π∫0∞‍ζω((ζ2+ζ2S+a)2+ω2)cos(ωτ)sin(ζξ)dζ+2π∫0∞‍ζ(ζ2+ζ2S+a)ω((ζ2+ζ2S+a)2+ω2)sin(ωτ)sin(ζξ)dζ+2π∫0∞‍ζSωe-(ζ2+ζ2S+a)τ((ζ2+ζ2S+a)2+ω2)sin(ζξ)dζ-2π∫0∞‍ζSω((ζ2+ζ2S+a)2+ω2)cos(ωτ)sin(ζξ)dζ+2π∫0∞‍ζS(ζ2+ζ2S+a)((ζ2+ζ2S+a)2+ω2)sin(ωτ)sin(ζξ)dζ.
The nondimensional shear stresses τcand τsat the plate ξ=0, due to cosine and sine oscillations of the plate, are obtained from (18) and (19) as follows:
(25)τc=-(τx+iτy)=-H(τ)e-at-iωτ4πτ×[(-4eiωτ-Sπτeaτ-Sπτeaτ+2iωτ2(a+iω)πτeaτ+2iωτerfc((a+iω)τ))-2eaτ(a-iω)πτ-2eaτ+2iωτ(a+iω)πτ+2eaτ(a-iω)πτerfc((a-iω)τ)+2(a+iω)πτeaτ+2iωτerfc((a+iω)τ))],τs=-(τx+iτy)=ie-at-iωτ4πτ×[(-4eiωτ-Sπτeaτ-Sπτeaτ+2iωτ2(a+iω)πτeaτ+2iωτerfc((a+iω)τ)-2eaτ(a-iω)πτ-2eaτ+2iωτ(a+iω)πτ+2eaτ(a-iω)πτerfc((a-iω)τ)+2(a+iω)πτeaτ+2iωτerfc((a+iω)τ))].


## 5. Special Cases

In order to verify the correctness of our obtained solutions, it is important to note that (18) and (19) satisfy the governing equation and the imposed initial and boundary conditions. Further, these solutions are more general comparing to the existing solutions in the literature. In this section, we derive some special cases of these solutions. The solution given by (18) for hydrodynamic fluid (M=0) over an impulsively moved plate (ω=0) reduces to the following form:
(26)G(ξ,τ)=H(τ)e-ξS/22×[e-ξS2/4+1/K+2iη×erfc(ξ2τ+(S24+1K+2iη)τ),e-ξS24+1K+2iη×erfc(ξ2τ-(S24+1K+2iη)τ)+eξS2/4+1/K+2iη×erfc(ξ2τ+(S24+1K+2iη)τ)e-ξS24+1K+2iη],
which is identical to the solution given by (13),  obtained by Jana et al. [21]. Further, the present solution given by (18) for suddenly moved plate (ω=0) reduces to the solution (19) obtained by Farhad et al. [20] in the absence of Hall current and slip boundary condition. Hence this provides a useful mathematical check to our calculi.


## 6. Results and Discussion

The exact solutions for the unsteady hydromagnetic flow of viscous fluid bounded by a porous plate in a porous medium are obtained. The analytical results are displayed for various values of emerging parameters such as Hartmann numberM, permeability parameter K, rotation parameter η, suction parameter S, and phase angle ωτ. Figures 2–6 are plotted for the cosine velocity given by (18). In these figures, (a) and (b) show the real and imaginary parts of velocity.Profiles of velocity for different values ofM.
(a)
(b)Profiles of velocity for different values ofK.
(a)
(b)Profiles of velocity for different values ofη.
(a)
(b)Profiles of velocity for different values ofS.
(a)
(b)Profiles of velocity for different values ofωτ.
(a)
(b)Figure2 is prepared to show the variation of velocity for different values of Hartmann number M. It is found that the real part of velocity and boundary layer thickness decreases with increasing values of Hartmann number. It is due to the fact that the application of transverse magnetic field results a resistive type force (called Lorentz force) similar to drag force and upon increasing the values of M, the drag force increases which leads to the deceleration of the flow. However, it is observed that this behavior is quite opposite for the imaginary part of velocity. Figure 3 reveals that the presence of porous medium increases the resistance to flow thus reducing its velocity. It may also be expected due to the fact that increasing values of Kreduces the friction forces which assists the fluid considerably to move fast and increases the boundary layer thickness. The graphs for the rotation parameter are plotted in Figure 4. It is observed that the real part of velocity increases when ηis increased. However, when observed carefully, it is found from Figure 4(b) that the behavior of imaginary part of velocity is more oscillatory in nature. The velocity increases first and then decreases when ηis increased.Figure5 shows the effect of suction parameter Son the flow through a porous medium in a rotating frame. Clearly, the absolute value of velocity and boundary layer thickness for both real and imaginary parts of velocity decrease with an increase in suction parameter. The variation of velocity for different values of phase angle ωτis shown in Figure 6. It is found that the real part of velocity decreases with increasing phase angle ωτ. Two different values of the phase angle ωτ=0π and π/2 are chosen. It is interesting to note that, when ωτ=0π, the real part of velocity is 1 which corresponds to the impulsive motion of the plate. Moreover, for ωτ=π/2, the real part of velocity is 0. The absolute value of the imaginary part of velocity is increasing with increasing phase angle ωτ. Furthermore, for large values of ξ, that is, ξ→∞, both of the real and imaginary parts of velocity are approaching zero. The graphical behavior in this figure is in accordance with the imposed boundary conditions (3) and (4). Moreover, it is interesting to note from Figure 6 that for a hydrodynamic fluid and impulsive motion of the plate, that is, when M=0 and ω=0, the graph for velocity matches with that of Jana et al. [21], which ensures the accuracy of the obtained analytical results.


## 7. Conclusions

The unsteady hydromagnetic rotating flow of viscous fluid bounded by an oscillating porous plate embedded in a porous medium is studied. The Laplace transform and Fourier sine transform methods are used for finding the solutions of the problem. The analytical results for nondimensional velocity and skin friction are obtained. The graphical results are plotted. The results show that velocity increases with increasing rotation parameter and permeability parameter whereas it decreases with increasing Hartmann number, suction parameter, and phase angle. Moreover, as the permeability of the medium increases, the velocity field increases in the boundary layer. Thus we can control the velocity field by introducing porous medium in a rotating system.



---

*Source: 508431-2013-12-17.xml*