Unsteady electromagnetic free convection flows of two-dimensional micropolar fluid through in a porous medium parallel to a vertical porous plate have been investigated numerically. Similarity analysis has been used to transform the governing equations into its non-dimensional form by using the explicit finite difference method to obtain numerical solutions. Estimated results have been gained for various values of Prandtl number, Grashof number, material parameters, micropolar parameter, electric conductivity, electric permeability, thermal relaxation time and the permeability of the porous medium. The effect
s
of pertinent parameters on the velocity, electric induction, magnetic induction, microrotation and temperature distributions have been investigated briefly and illustrate
d
graphically.
Micropolar Fluid Free Convection Porous Medium Explicit Finite Difference1. Introduction
Fluids with microstructure are micropolar fluids which are randomly oriented or composed of spherical particles that are rigid with their rotation and also ceased in a viscous medium. It has been known that Navier-Stokes equations are unable to explain the phenomena at micro and nanoscales; on the other hand, MFD can express the physical phenomena at micro and nanoscales owing to its additional degree of freedom for circulation. Physical examples of micropolar fluids may present in the non-Newtonian fluids, blood flows, polymer fluids and liquid crystals and all of them containing intrinsic polarities. The presence of dust or fumes in a gas can be especially modeled using micropolar fluid dynamics. The porous media heat transfer problems have various practical uses in engineering applications such as geothermal systems, crude oil extraction and groundwater pollution. Eringen first proposed [1] and [2] the general theory of micropolar fluids which illustrate certain microscopic effects arising from the microstructure and micro motions of the fluid flow. The interaction of natural convection with thermal radiation in laminar boundary layer flow over an isothermal, horizontal flat plate is studied by Ali et al. [3]. Harutha and Devasena [4] investigated the steady mixed convection flow of a viscous incompressible micropolar fluid through a porous medium towards a stagnation point over a vertical surface when the buoyancy forces assist. Hudimoto and Tokuoka [5] have devised the two-dimensional parallel shear flow of a linear micropolar fluid. They analyzed and compared it with the colloidal suspensions. Rees and Pop [6] expressed the steady micropolar free convection fluid flow from a vertical isothermal flat plate. Elbarbary [7] discussed a new Chebyshev finite difference method is proposed for solving the governing equations of the boundary layer flow. Nandhini and Ramya [8] analyzed the heat and mass transfer of the free convection flow in a micropolar fluid past an inclined stretching sheet. Hassanien and Glora [9] analyzed the heat transfer on a non-isothermal stretching sheet to a micropolar fluid. Kartini Ahmad et al. [10] described a micropolar fluid flow and heat transfer past a non-linearly stretching plate. Khonsari and Brewe [11] investigated and compared the parameters of micropolar fluids with finite length lubricated that resulted significantly higher load carrying capacity than Newtonian fluids. Effects of free convection currents with one relaxation time on the flow of a viscoelastic conduction fluid through a porous medium, which is bounded by a vertical plane surface, have studied by Ezzatand Abd-Ellal [12]. Edlabe and Mohammed [13] determined the heat and mass transfer occurring in the hydromagnetic flow of the non-Newtonian fluid on a linearly accelerating surface with temperature dependent heat source subject to suction or blowing. Edlabe and Ouaf [14] are obtained the heat and mass transfer in a hydro magnetic flow of a micropolar fluid past a stretching surface with Ohmic heating and viscous dissipation. Aydin and Pop [15] analyzed the two-dimensional steady laminar natural convective flow and heat transfer of micropolar fluids in a square enclosure. Muthu et al. [16] investigated the oscillatory flow of micropolar fluid in an annular region with constriction, provided by variation of the outer tube radius. Glora [17] presented an unsteady combined convection of a micropolar fluid among a vertical plate. Hsu and Wang [18] presented a numerical study of the laminar mixed convection of micropolar fluids in a square cavity with localized heat source Lok et al. [19] studied a microplar mixed convection boundary layer fluid flow near the region of the stagnation point of on a double-infinite vertical flat plate. Zakaria [20] investigated the influence of a transverse magnetic field on the motion of an electrically conducting micropolar fluid through a porous medium in one-dimensional and used the Laplace transformation with ε -algorithm technique to find its solution in the Laplace transformation domain numerically.
In the present work, our aim is to study that the numerical investigation on unsteady two-dimensional electromagnetic free convection micropolar fluid flows through a porous medium along a vertical porous plate. The obtained equations are non-linear coupled partial differential equations, which are solved by using explicit finite difference method and the results are shown graphically and also discussed its behavior in detail for the velocity, induced magnetic field, induced electric field, micro rotation and temperature distribution with respect to its pertinent parameters.
2. Problem Formulation
Considered unsteady MHD micropolar fluid flow embedded in a porous medium along a vertical porous plate. The velocity at the wall is zero and also outside of the boundary layer is zero. The temperature of the plate is raised from T w to T ∞ , where T w and T ∞ is the temperature at the plate and outside of the boundary layer respectively. The magnetic Reynolds number is taken large enough so that the induced magnetic field equation is considerable for our assumption. The Physical model of the system is shown in the following Figure 1.
The flow is governed by the equation of continuity, the momentum equation, induction magnetic field equation, the electric field equation, angular momentum equation and the energy equation are as follows:
∂ u + ∂ x + + ∂ v + ∂ y = 0 (1)
∂ u + ∂ t + + u + ∂ u + ∂ x + + v + ∂ u + ∂ y = g β ( T + − T ∞ + ) + ( μ + μ * ρ ) ∂ 2 u + ∂ y + 2 + μ * ρ ∂ N + ∂ y + + α 2 H 0 ( ∂ H x + ∂ y + + ε 0 + ∂ E + ∂ t + ) − μ ρ κ u + (2)
∂ H x + ∂ t + + u + ∂ H x + ∂ x + + v + ∂ H x + ∂ y + = ν m ∂ 2 H x + ∂ y + 2 − ε 0 + σ 0 + ∂ 2 H x + ∂ t + 2 + H 0 ∂ u + ∂ y + (3)
∂ E + ∂ t + + u + ∂ E + ∂ x + + v + ∂ E + ∂ y + = − σ 0 + ε 0 + E + − 1 ε 0 + ∂ H x + ∂ y + − H 0 σ 0 + μ 0 ε 0 + u + (4)
∂ N + ∂ t + + u + ∂ N + ∂ x + + v + ∂ N + ∂ y + = γ ρ j ∂ 2 N + ∂ y + 2 − μ * ρ j ( ∂ u + ∂ y + + 2 N + ) (5)
∂ T + ∂ t + + u + ∂ T + ∂ x + + v + ∂ T + ∂ y + = k ρ c p ∂ 2 T + ∂ y + 2 − τ 0 + ∂ 2 T + ∂ t + 2 (6)
with boundary conditions are as follows:
y + = 0 : u + = 0 , N + = 0 , T + = T w + , H x + = 0 y + → ∞ : u + → 0 , N + → 0 , T + → T ∞ + , H x + → 0 (7)
Similarity Analysis
Now introducing the following non-dimensional quantities as
x = α ρ μ x + , y = α ρ μ y + , u = u + α , w = w + α , H = H x + H 0 , E = E + μ 0 H 0 α ,
t = α 2 ρ μ t + , θ = T + − T ∞ + T w + − T ∞ + , N = μ α 2 ρ N + , ε 0 = μ 0 α 2 ε 0 + , σ 0 = μ μ 0 ρ σ 0 +
and τ 0 = α 2 ρ μ τ 0 + .
Using these quantities into the above Equations (1)-(7), we obtain the following dimensionless form of the given equations:
∂ u ∂ x + ∂ v ∂ y = 0 (8)
∂ u ∂ t + u ∂ u ∂ x + v ∂ u ∂ y = G r θ + ( 1 + R ) ∂ 2 u ∂ y 2 + R ∂ N ∂ y + ∂ H ∂ y + ε 0 ∂ E ∂ t − u K (9)
∂ H ∂ t + u ∂ H ∂ x + v ∂ H ∂ y = 1 b ∂ 2 H ∂ y 2 − ε 1 b ∂ 2 H ∂ t 2 + ∂ u ∂ y (10)
∂ E ∂ t + u ∂ E ∂ x + v ∂ E ∂ y = − σ 0 ε 0 E − σ 0 ε 0 u − 1 ε 0 ∂ H ∂ y (11)
∂ N ∂ t + u ∂ N ∂ x + v ∂ N ∂ y = λ ∂ 2 N ∂ y 2 − σ ( ∂ u ∂ y + 2 N ) (12)
∂ θ ∂ t + u ∂ θ ∂ x + v ∂ θ ∂ y = 1 P r ∂ 2 θ ∂ y 2 − τ 0 ∂ 2 θ ∂ t 2 (13)
with the corresponding boundary conditions:
y = 0 : u = 0 , N = 0 , θ = 1 , H = 0 y → ∞ : u → 0 , N → 0 , θ → 0 , H → 0
where G r = β g μ ( T w + − T ∞ + ) / ( ρ α 3 ) is the Grashof number, P r = c p μ / k is the Prandtl number, R = μ * / μ is the micropolar parameter, K = ρ 2 α 2 κ / μ 2 is the Permeability parameter. Also σ = μ μ * / ( α 2 j ρ 2 ) , b = μ / ρ ν m , ε 1 = b ε 0 / σ 0 and λ = γ / j μ are the dimensionless material parameter.
3. Method of Solution
The explicit finite difference method has been used to solve the governing non-linear coupled dimensionless partial differential Equations (8) to (13) together with its boundary conditions. The finite difference schemes with respect to t, x and y are as follows:
∂ u ∂ t = U i , j k + 1 − U i , j k Δ t ; ∂ u ∂ x = U i , j k − U i − 1 , j k Δ x ;
∂ u ∂ y = U i , j k − U i , j − 1 k Δ y ; ∂ 2 H ∂ t 2 = h i , j k + 2 − 2 h i , j k + 1 + h i , j k Δ t 2 ;
∂ 2 u ∂ y 2 = U i , j + 1 k − 2 U i , j k + U i , j − 1 k Δ y 2
Here, the subscript i and j refer to x and y and the superscript k refers to time t. Finite difference Schemes for the other variables have been written in the same way. The graphical representations of this problem have been illustrated by using Compaq visual FORTRAN 6.6 a tools.
4. Results and Discussion
The behavior of the velocity (u), induced magnetic field (H), induced electric field (E), microrotation (N) and temperature ( θ ) distributions have been analyzed for the different values of Prandtl number ( P r ), Grashof number ( G r ), permeability of porous medium (K), micropolar parameter (R) and thermal relaxation time ( τ 0 ) with the values of time t = 1 . The flow characteristics have been shown graphically from Figures 2-25.
4.1. Time and Mesh Sensitivity Test
To get the steady-state solution, the computations are carried out for different time t = 20 , 25 , 29 and 30 with time increment Δ t = 0.001 for the velocity distribution, which have shown in Figure 2. It is found that after t = 30 , there are very negligible changes. On the other hand, to choose the appropriate mesh, a solutions find out for different pairs of meshes like as ( m , n ) = ( 70 , 70 ) ; ( m , n ) = ( 80 , 80 ) and ( m , n ) = ( 90 , 90 ) on the velocity distributions, which have shown in Figure 3. For those three chosen different values of meshes, the profiles are likely unchanged. There are same situations for the other distributions. Therefore our estimated steady-state time is at t = 30 with time increment Δ t = 0.001 and mesh pair is ( m , n ) = ( 80 , 80 ) with the fixed values of G r = 1.0 ; P r = 7.0 ; R = 7.0 ; τ 0 = 0.02 ; σ = 0.1 ; σ 0 = 1.0 ; λ = 0.2 ; ε 0 = 1.0 ; K = 1.0 and b = 0.2 .
4.2. Comparison with Previous Results
Zakaria [20] investigated the influence of the Grashof number on the velocity u, which is shown in Figure 4(a). Here the velocity decreases with the increase of Grashof number Gr. But in Figure 4(b), it is found that the velocity increases with the same increasing values of Grashof number. In this case the maximum time has taken t = 1 .
4.3. Primary Velocity Distributions
Figure 5 depicts that the velocity u is increased with the increase of K. In Figure 6, it is observed that the velocity u is decreased with the increasing values of Pr. But in Figure 7, showed a cross-flow for the velocity, here velocity distribution is decreased within the interval 0 < y < 6 (approx.) and thereafter it has very minor increasing effect with the increase of R. Figure 8 represented that the velocity has an increasing effect with the increase of t0.
4.4. Induced Magnetic Field Distributions
Figure 9 and Figure 10 illustrate that the induced magnetic field distribution H has a cross-flow for the different values of Gr. It is obvious that near the plate, H has a minor increasing effect and thereafter found a large decreasing effect for increasing values of Gr and K. Figure 11 represents that H has a very minor increasing effect near the plate and thereafter a decreasing effect with the increase of Pr. But from Figure 12, it is observed that H has an increasing effect with the rising values of R.
4.5. Induced Electric Field Distributions
Profiles in Figure 13 and Figure 14, represented that the induced electric field E is decreased with the increase of Gr and K respectively. But E has an increasing effect with the rising values of Pr which is shown in Figure 15.
4.6. Microrotation Distributions
The microrotation N has a cross-flow depicts in Figure 16. It has a decreasing effect within 0 < y < 2.9 (approx.) and thereafter it has an increasing effect with the increase of K. But Figure 17 noticed that the increasing values of R, the micropolar rotation N has a decreasing effect within 0 < y < 2.2 (approx.) and then it has an increasing effect.
4.7. Temperature Distributions
Figure 18 displays the effect of the Prandtl number Pr on the temperature ( θ ). As shown, temperature is decreasing with the increasing of Pr.
4.8. Skin-Friction, Current Density and Rate of Heat Transfer
The effects of various parameters on local and average shear stress from the velocity profile have been investigated. The non-dimensional form of the local shear stress and average shear stress in x-direction is given by the relations τ L = μ ∂ u ∂ y | y = 0 and τ A = 1 L ∫ 0 L μ ∂ u ∂ y | y = 0 d x respectively. From the temperature profile, the effects of various parameters on local and average Nusselt numbers have been calculated. The local Nusselt number and the average Nusselt number are given by N u L = − μ ∂ θ ∂ y | y = 0 and N u A = − 1 L ∫ 0 L μ ∂ θ ∂ y | y = 0 d x respectively. Similarly, analyze the effects of various parameters on the local and average Sherwood numbers from the concentration field. The rate of mass transfer at the plate is defined as the Sherwood number; the local Sherwood number and the average Sherwood number is defined by S h L = − μ ∂ φ ∂ y | y = 0 and S h A = − 1 L ∫ 0 L μ ∂ φ ∂ y | y = 0 d x respectively.
Figures (a) and (b) of Figures 19-22 are displayed the variations of the local shear stress and average shear stress respectively. It is obtained from Figure 19 and Figure 20 that the local (or average) shear stress increases with the increase of Gr and K. The same effects of velocity are also represented. But from Figure 21 and Figure 22, show that both of the local and average shear stress decreases with the increase of Pr and R respectively.
Again, figures (a) and (b) of Figure 23, Figure 24 showed the variations of the local and average current density respectively. It is presented that JwL (or Jw) has both decreasing effect with the increasing value of Gr and K. But Figure 25(a), Figure 25(b) have shown the variation of local (or average) Nusselt number for different values of the Prandtl number (Pr), it depicts from the figure that Nusselt number increases with the increase of Pr.
5. Conclusions
In the present study, the influence of various values of Prandtl number, Grashof number, permeability parameter, micropolar parameter, electric conductivity, electric permeability and thermal relaxation time has been investigated. The non-linear coupled governing equations have been solved numerically and the main findings can be summarized as follows:
1) The velocity u increases with the increase of Gr, K and τ 0 , while it decreases with the increase of Pr and R.
2) Induced magnetic field H has cross-flow near the plate. But in major space, it has been increasing effect with the increase of Pr and R, while it decreases with the increase of Gr and K.
3) Induced electric field E increases with the increase of Pr, while it decreases with the increase of Gr and K.
4) Microrotation N has cross-flow for all the different values of all the parameters. First portion near the plate N has increasing effect with the increase of R. Thereafter, it has a decreasing effect. But for Gr and K, it has reverse effect.
5) Temperature θ decreases with the increase of Pr.
6) Local (or average) Shear stress increases of Gr and K, while it decreases with the increase of Pr and R.
7) Local (or average) Current density decreases with the increase of Gr and K.
8) Local (or average) Nusselt number increases with the increase Pr.
The accuracy of this work is qualitatively good in case of all the flow parameters.
Conflicts of Interest
The authors declare no conflicts of interest regarding the publication of this paper.
Cite this paper
Islam, M.R., Nasrin, S. and Alam, M.M. (2020) Unsteady Electromagnetic Free Convection Micropolar Fluid Flow through a Porous Medium along a Vertical Porous Plate. Open Journal of Applied Sciences, 10, 701-718. https://doi.org/10.4236/ojapps.2020.1011049
NomenclatureReferencesEringen, A.C. (1966) Theory of Micropolar Fluids. Journal of Mathematics and Mechanics, 16, 1-18. https://doi.org/10.1512/iumj.1967.16.16001Eringen, A.C. (1972) Theory of Thermo Micropolar Fluids. Journal of Mathematical Analysis and Applications, 38, 480-496. https://doi.org/10.1016/0022-247X(72)90106-0Ali, M.M., Chen, T.S. and Armaly, B.F. (1984) Natural Convection Radiation Interaction in Boundary Layer Flow over Horizontal Surface. AIAA Journal, 22, 1797-1803. https://doi.org/10.2514/3.8854Harutha, A. and Devasena, Y. (2016) MHD Mixed Convection Flow of a Micropolar Fluids through Porous Medium towards a Stagnation, Point on a Vertical Porous Surface. IOSR Journal of Mathematics, 12, 32-37. https://doi.org/10.9790/5728-1204043237Hudimoto, B. and Tokuoka, T. (1969) Two Dimensional Shears Flows of Linear Micro Polar Fluids. International Journal of Engineering Science, 7, 515-522. https://doi.org/10.1016/0020-7225(69)90036-6Rees, D.A.S. and Pop, I. (1998) Free Convection Boundary Layer Flow of a Micropolar Fluid from a Vertical Flat Plate. IMA Journal of Applied Mathematics, 61, 179-197. https://doi.org/10.1093/imamat/61.2.179Elbarbary, E.M.E. (2005) Chebyshev Finite Difference Method for the Solution of Boundary Layer. Applied Mathematics and Computation, 160, 487-498. https://doi.org/10.1016/j.amc.2003.11.016Nandhini, E. and Ramya, M. (2018) MHD Free Convection Flow in a Micropolar Fluid past an Inclined Stretching Sheet with Considering Viscous Dissipation and Radiation. International Journal of Scientific Research in Science, Engineering and Technology, 4, Issue1.Hassanien, I.A. and Glora, R.S.R. (1990) Heat Transfer to a Micropolar Fluid from a Non-Isothermal Stretching Sheet with Suction and Blowing. Acta Mechanica, 84, 191-199. https://doi.org/10.1007/BF01176097Ahmad, K., Ishak, A. and Nazar, R. (2013) Micropolar Fluid Flow and Heat Transfer over a Nonlinearly Stretching Plate with Viscous Dissipation. Mathematical Problems in Engineering, 2013, Article ID: 257161. https://doi.org/10.1155/2013/257161Khonsari, M.M. and Brewe, D. (1989) On the Performance of Finite Journal Bearings Lubricated with Micro Polar Fluids. Tribiology Transactions, 32, 155-160. https://doi.org/10.1080/10402008908981874Ezzat, M.A. and Abd-Ellal, M.Z. (1997) Free Convection Effects on a Viscoelastic Boundary Layer Flow with One Relaxation Time through a Porous Medium. Journal of Franklin Institute, 334, 685-706. https://doi.org/10.1016/S0016-0032(96)00095-6Edlabe, N.T.M. and Mohammed, M.A.A. (2002) Heat and Mass Transfer in Hydrodynamic Flow of the Non-Newtonian Fluid with Heat Source over an Accelerating Surface through a Porous Medium. Chaos, Solitons and Fractals, 13, 907-917. https://doi.org/10.1016/S0960-0779(01)00066-2Edlabe, N.T.M. and Ouaf, M.E.M. (2006) Chebyshev Finite Difference Method for Heat and Mass Transfer in a Hydrodynamic Flow of a Micropolar Fluid past a Stretching Surface with Ohmic Heating and Viscous Dissipation. Applied Mathematics and Computation, 177, 561-571. https://doi.org/10.1016/j.amc.2005.07.071Aydin, O. and Pop, L. (2007) Natural Convection in a Differentially Heated Enclosure Filled with a Micropolar Fluid. International Journal of Thermal Science, 46, 963-969. https://doi.org/10.1016/j.ijthermalsci.2006.11.018Muthu, P., Rathish Kumar, B.V. and Chandra, P. (2008) A Steady of Micropolar Fluid in an Annular Tube with Application to Blood Flow. Journal of Mechanics in Medicine and Biology, 8, 561-576. https://doi.org/10.1142/S0219519408002541Glora, R.S.R. (1995) Unsteady Mixed Convection in Micropolar Boundary Layer Flow on a Vertical Plate. Fluid Dynamics Research, 15, 237-250. https://doi.org/10.1016/0169-5983(95)94957-UHsu, T.H. and Wang, S.G. (2000) Mixed Convection of Micropolar Fluids in a Cavity. International Journal of Heat and Mass Transfer, 43, 1563-1572. https://doi.org/10.1016/S0017-9310(99)00242-2Lok, Y.Y., Amin, N. and Pop, I. (2006) Unsteady Mixed Convection Flow of a Micropolar Fluid near the Stagnation Point on a Vertical Surface. International Journal of Thermal Science, 45, 1149-1157. https://doi.org/10.1016/j.ijthermalsci.2006.01.015Zakaria, M. (2004) Problem in Electromagnetic Free Convection Flow of a Micropolar Fluid with Relaxation Time through a Porous Medium. Applied Mathematics and Computation, 151, 601-613. https://doi.org/10.1016/S0096-3003(03)00365-5