MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 508431 10.1155/2013/508431 508431 Research Article Unsteady Hydromagnetic Rotating Flow through an Oscillating Porous Plate Embedded in a Porous Medium Khan I. 1,2 Khan A. 1 Farhad A. 1 Qasim M. 3 Sharidan S. 1 Chinyoka Tirivanhu 1 Department of Mathematical Sciences Faculty of Science Universiti Teknologi Malaysia 81310 Skudai Malaysia utm.my 2 College of Engineering Majmaah University P.O. Box 66, Majmaah 11952 Saudi Arabia utm.my 3 Department of Mathematics COMSATS Institute of Information Technology Park Road Chak Shahzad Islamabad 44000 Pakistan ciit.edu.pk 2013 17 12 2013 2013 23 07 2013 02 10 2013 09 10 2013 2013 Copyright © 2013 I. Khan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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.

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 . It admits an exact analytical solution. The Stokes’ problem not only is of fundamental theoretical interest but also occurs in many applied problems . The transient solutions in terms of tabulated functions for the flow of viscous fluid due to oscillating boundary have been expressed by Panton . Later, in 2000, Erdogan  studied the unsteady flow of Newtonian fluid past an oscillating horizontal plane wall and obtained the exact solutions. Fetecau et al.  presented new exact solutions corresponding to Stokes’ second problem for both small and large times. Erdoǧan and Imrak  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 , Ingham et al. , Ingham and Pop , Vafai , and Nield and Bejan  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. ). According to Hayat et al. , 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  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. , Sahoo et al. , Seth et al. , and Farhad et al.  studied the hydromagnetic flow of rotating fluids in porous media. Recently, Jana et al.  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.  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.  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=0  such 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 B0  applied in a direction normal to the plate by means of the plate oscillations in its plane (Hayat et al. ). 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 t0, 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  (1)ut-w0uz-2Ωv=ν2uz2-νϕku-σB02ρu,z,t>0,vt-w0vz+2Ωu=ν2vz2-νϕkv-σB02ρv,z,t>0, where  u  and v  are the velocity components in x- and y-directions,  w0  is the suction velocity at the plate, ν is the kinematic viscosity, t is the time, σ is the electrical conductivity, B0  is the applied magnetic field, ρ is the fluid density, ϕ is the porosity, and k is the permeability of the porous medium.

Geometry of the problem.

The corresponding initial and boundary conditions are (2)u=v=0;t0,(3)u=UH(t)cos(ω0t)orUsin(ω0t);v=0;atz=0;t>0,(4)u0;v0asz;t>0, in which U  is the constant plate velocity,  ω0  is the frequency of oscillation of the plate, and H(t)  is the Heaviside function.

Defining F=u+iv, the above system of equations reduces to (5)Ft-w0Fz+2iΩF=ν2Fz2-νϕkF-σB02ρF,z,t>0,F(z,0)=0;t0,F(0,t)=UH(t)cos(ω0t)orUsin(ω0t);t>0F0asz;t>0. Introducing the following dimensionless variables (see Jana et al. ): (6)G=FU,ξ=zUν,τ=U2tν,ω=ω0νU2, the dimensionless problem takes the following form: (7)Gτ-SGξ+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 S  are 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 q  is 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 c  and s  denote 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. ; 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 τc  and τs  at 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. . Further, the present solution given by (18) for suddenly moved plate (ω=0) reduces to the solution (19) obtained by Farhad et al.  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 number M, permeability parameter K, rotation parameter η, suction parameter S, and phase angle ωτ. Figures 26 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 of M.

Profiles of velocity for different values of K.

Profiles of velocity for different values of η.

Profiles of velocity for different values of S.

Profiles of velocity for different values of ωτ.

Figure 2 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 K  reduces 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.

Figure 5 shows the effect of suction parameter S  on 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. , 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.

Schlichting H. Grenzchicht Theorie 1982 8th ed Karlsruhe, Germany Braun Tokuda N. On the impulsive motion of a flat plate in a viscous fluid Journal of Fluid Mechanics 1968 33 4 657 672 10.1017/S0022112068001606 Panton R. The transient solution for Stokes's oscillating plane: a solution in terms of tabulated functions Journal of Fluid Mechanics 1968 31 4 819 825 10.1017/S0022112068000509 Erdogan M. E. Note on an unsteady flow of a viscous fluid due to an oscillating plane wall International Journal of Non-Linear Mechanics 2000 35 1 1 6 2-s2.0-0033687788 10.1016/S0020-7462(99)00019-0 MR1723758 ZBL1006.76028 Fetecau C. Vieru D. Fetecau C. A note on the second problem of Stokes for Newtonian fluids International Journal of Non-Linear Mechanics 2008 43 5 451 457 2-s2.0-43249131215 10.1016/j.ijnonlinmec.2007.12.022 Erdoǧan M. E. Imrak C. E. On the comparison of the solutions obtained by using two different transform methods for the second problem of Stokes for Newtonian fluids International Journal of Non-Linear Mechanics 2009 44 1 27 30 2-s2.0-56949093137 10.1016/j.ijnonlinmec.2008.08.004 Pop I. Ingham D. B. Convective Heat Transfer: Mathematical and Computational Modeling of Viscous Fluids and Porous Media 2001 Oxford, UK Pergamon Ingham D. B. Bejan A. Mamut E. Pop I. Emerging Technologies and Techniques in Porous Media 2004 Dodrecht, The Netherlands Kluwer Academic Ingham D. B. Pop I. Transport Phenomena in Porous Media 2005 Oxford, UK Pergamon Vafai K. Handbook of Porous Media 2005 New York, NY, USA Taylor & Francis Nield D. A. Bejan A. Convection in Porous Media 2006 3rd New York, NY, USA Springer MR1656781 Jana M. Das S. Jana R. N. Unsteady Couette flow through a porous medium in a rotating system Open Journal of Fluid Dynamics 2012 2 149 158 10.4236/ojfd.2012.24016 Hayat T. Nadeem S. Siddiqui A. M. Asghar S. An oscillating hydromagnetic non-Newtonian flow in a rotating system Applied Mathematics Letters 2004 17 6 609 614 2-s2.0-3242729440 10.1016/S0893-9659(04)90134-6 MR2057358 ZBL1060.76130 Das S. Maji S. L. Guria M. Jana R. N. Unsteady MHD Couette flow in a rotating system Mathematical and Computer Modelling 2009 50 7-8 1211 1217 2-s2.0-69249220121 10.1016/j.mcm.2009.05.036 MR2568111 ZBL1185.76915 Hayat T. Fetecau C. Sajid M. Analytic solution for MHD Transient rotating flow of a second grade fluid in a porous space Nonlinear Analysis: Real World Applications 2008 9 4 1619 1627 2-s2.0-43649092634 10.1016/j.nonrwa.2007.04.006 MR2422567 ZBL1154.76391 Hayat T. Fetecau C. Sajid M. On MHD transient flow of a Maxwell fluid in a porous medium and rotating frame Physics Letters A 2008 372 10 1639 1644 2-s2.0-38949198352 10.1016/j.physleta.2007.10.036 Abelman S. Momoniat E. Hayat T. Steady MHD flow of a third grade fluid in a rotating frame and porous space Nonlinear Analysis: Real World Applications 2009 10 6 3322 3328 2-s2.0-68749085867 10.1016/j.nonrwa.2008.10.067 MR2561342 ZBL1269.76004 Sahoo S. N. Panda J. P. Dash G. C. Hydromagnetic oscillatory flow and heat transfer of a viscous liquid past a vertical porous plate in a rotating medium Indian Journal of Science and Technology 2010 3 817 821 Seth G. S. Ansari M. S. Nandkeolyar R. Unsteady hydromagnetic Couette flow within porous plates in a rotating system Advances in Applied Mathematics and Mechanics 2010 2 3 286 302 MR2610416 Farhad A. Norzieha M. Sharidan S. Khan I. Samiulhaq Hydromagnetic rotating flow in a porous medium with slip condition and Hall current International Journal of Physical Sciences 2012 7 10 1540 1548 10.5897/IJPS11.1767 Jana M. Maji S. L. Das S. Jana R. N. Unsteady flowof viscous fluid through a porous medium bounded bya porous plate in a rotating system Journal of Porous Media 2010 13 7 645 653 2-s2.0-77958558347 10.1615/JPorMedia.v13.i7.60