Finite Element Solution of an Unsteady MHD Flow through Porous Medium between Two Parallel Flat Plates

Finite element solution of unsteady magnetohydrodynamics (MHD) flow of an electrically conducting, incompressible viscous fluid past through porous medium between two parallel plates is presented in the presence of a transverse magnetic field and Hall effect.The results obtained from some test cases are then compared with previous publishedwork using the finite differencemethod (FDM). Numerical examples show that the finite element method (FEM) gives more accurate results in comparison with the finite difference method (FDM).


Introduction
Theoretical study of magnetohydrodynamics (MHD) flow problems are frequently encountered in cooling systems of nuclear reactors, MHD generators, blood flow measurements, pumps, and accelerators.
Due to coupling of the equations for electrodynamics and fluid mechanics, exact solution is possible only for some simple situations.By using several numerical techniques, such as finite element method (FEM), finite volume method (FVM), and boundary element method (BEM), approximate solution for the MHD flow problems can be obtained.
Gapta and Singh [1] obtained the exact solutions for unsteady flow in some special cases.Ram and Mishra [2] investigated the unsteady flow through magnetohydrodynamic porous media.Singh and Lal [3] studied the FEM solution of time-dependent MHD flow equations.Ram and Jain [4] have discussed MHD free convective flow through a porous medium in a rotating fluid.Reddy and Bathaiah [5] have analyzed the Hall effects on MHD flow through a porous straight channel.Lee and Dulikravich [6] proposed FDM scheme for the 3-dimensional unsteady MHD flow together with temperature field.Sheu and Lin [7] presented a convection-diffusion-reaction model for solving the unsteady MHD flow using a FDM scheme.The stabilized FEM for solution of the 3-dimensional time-dependent MHD flow equations was given by Ben Salah and et al. [8].Chauhan and Rastogi [9] have studied the Hall effects on MHD slip flow and heat transfer through a porous medium over an accelerated plate in a rotating system.Saha and Chakrabarti [10] have investigated the impact of magnetic field strength on magnetic fluid flow through a channel.Moniem and Hassanin [11] have developed a solution of MHD flow past a vertical porous plate through a porous medium under oscillatory suction.Sa'adAldin and Qatanani [12] have studied the unsteady MHD flow through two parallel porous flat plates.Sivaiah and Srinivasa-Raju [13] have discussed the finite element solution of heat and mass transfer in MHD flow of a viscous fluid past a vertical plate under oscillatory suction velocity.Yuksel and Ingram [14] have investigated the numerical analysis of a finite element method, Crank-Nicolson discretization for MHD flows at small magnetic Reynold number.Bég et al. [15] have developed a finite element and network electrical simulation of rotating magnetofluid flow in nonlinear porous media with inclined magnetic field and Hall currents.Sa'ad Aldin and Qatanani [16] have studied the analytical and finite difference methods for solving unsteady MHD flow through porous medium between two parallel flat plates.
In this work, the finite element solution for the unsteady magnetohydrodynamics (MHD) flow of an electrically conducting, incompressible viscous fluid past through porous medium between two parallel plates in the presence of a transverse magnetic field and Hall effect is considered.A comparison study has been carried out between the finite difference and the finite element solutions.A case study is analyzed with both the finite element method (FEM) and the finite difference method (FDM), namely, the implicit scheme presented in [16].It was found that the finite element method (FEM) is more accurate for solving these type of problems.

Formulation of the Problem
We consider an unsteady flow of an electrically conducting, incompressible viscous fluid past through porous medium between two parallel plates with Hall effect.Let the -axis be taken along the plates and -axis be normal to the plates.The fluid is subjected to a constant transverse magnetic field of strength  0 in the  direction, with the flow being considered in the  direction, as illustrated in Figure 1.The governing equations for the unsteady, viscous incompressible flow of an electrically conducting fluid for the Brinkman-extended Darcy model are as follows [2]: Equation of continuity is Equation of motion is General Ohm's law is Gauss's law of magnetism is where q is the velocity vector,  is the fluid density, p is the pressure, J is the current density, B is the magnetic vector,  is the coefficient of viscosity,  is the electrical conductivity,  is the permeability of the medium,  is the electron frequency,  is the electron collision time,   is the electric charge,   is the number density of electron, P  is the electron pressure, and  is the electric field.
We assume  to be negligible and the magnetic Reynold's number is small so that magnetic induction effect is ignored.Moreover, in the absence of pressure gradient, the ion-slip effects and electron pressure gradient, we have Solving ( 6) and ( 8), we have As the plates are infinite, there is no  dependence.Consequently, ( 2) and (3) take the following form: where  is the axial velocity, ] is the kinematic viscosity, and  =  is the Hall parameter.The initial and boundary conditions are given by  = 0,  ≤ 0,  = 0,  = ±ℎ,  > 0.
Upon introducing the nondimensional quantities, where  is the Hartman number,  is the Darcy parameter, and  is the mean velocity of the fluid.Then, the partial differential equations (10) together with the initial and boundary conditions (11) become subject to the initial and boundary conditions: In virtue of ( 14), the pressure is independent of ; then it is a function of  only.In this case, we can take the pressure gradient as a constant quantity; that is, where  0 > 0; thus (13) becomes subject to the initial and boundary conditions: where  = ( 2 /(1 +  2 ) + 1/).

Variational Formulations and Galerkin Approximation.
The dimensionless partial differential equation ( 17) subject to the initial and boundary conditions (18) is solved by weighted residual Galerkin finite element method.The standard approach to deriving a Galerkin scheme is to multiply both sides of ( 17) by a test function  ∈  1 0 [−1, 1] and integrate over the domain where Integrating by parts, we obtain where ⟨⋅, ⋅⟩ denotes the  2 -inner product and ⟨, ⟩ = 0 for  ≤ 0.
Here  is the Stiffness matrix and  is the Mas matrix defined, respectively, as with   (  ) =   being the usual finite element basis corresponding to the partition  ℎ .Thus, to compute the entries of the Stiffness matrix , next, we need to determine    (). Here and then where ℎ  =   −  −1 .

Time
Stepping.In this section, we consider the semidiscretization in time.We first discretize the time interval (0, ) into a uniform grid with size  = /.Approximating the derivative in (24) at time level   by the Crank-Nicolson scheme with  0 = 0, we have Then, we rewrite the Crank-Nicolson method as Thus, we have the full discretization which is simply a combination of discretization in space and time:

Numerical Results and Discussion
To show the efficiency of the FEM described in the previous parts and to draw a comparison with the FDM, we present   Table 1 compares the exact values for the velocity (, ) with both the FEM and the FDM values (for more details on the exact and FDM solutions, see [16]).A further comparison between the exact, FEM, and FDM values for the velocity can be observed in Figure 2. A plot of the absolute error that resulted from the FEM can be seen in Figure 3. Figure 4   compares the absolute errors obtained from the FEM and the FDM solutions.

Numerical Example 2.
As for another test case, we take  = 2,  = 1,  = 0.3,  0 = 1,  = 0.5 fixed, and  ∈ [0, 1]. Figure 5 compares the exact, FEM, and FDM values for the velocity.Figure 6 presents a plot of the absolute error that resulted from the FEM.A comparison between the absolute errors obtained from the FEM and the FDM solutions can be seen in Figure 7.

Conclusions
MHD flow problems, which have a very important place in physics and engineering, are usually hard to solve analytically.Therefore, it is required to obtain approximate solutions using computational methods.In this work, the problem of unsteady MHD flow through porous medium in the presence of magnetic field between two parallel flat plates has been investigated and solved using the FEM.A comparison between FFM and FDM has been carried out.The exact results and the numerical results using the FEM have shown to be in closed agreement.This can clearly be seen in Figures 2, 3, 5, and 6.The results of the numerical examples indicate that the FFM is more accurate than the FDM (see Figures 4 and 7).This asserts the ability and reliability of the FEM for solving these types of problems.

Figure 1 :
Figure 1: Schematic diagram of the system.

Figure 2 :
Figure 2: The exact, FEM, and FDM values for the velocity.

Figure 3 :Numerical Example 1 .
Figure 3: The absolute error resulted from the FEM approximation.

Figure 4 :
Figure 4: The absolute error resulted from the FEM and FDM approximations.

Figure 5 : 1 E
Figure 5: The exact, FEM, and FDM values for velocity case 2.

Figure 6 : 1 E
Figure 6: The absolute error resulted from FEM approximation case 2.

Figure 7 :
Figure 7: The absolute error that resulted from FEM and FDM approximations case 2.

Table 1 :
The exact, (FEM) and (FDM) solutions of the velocity . Exact solution   FEM solution  FE |  −  FE | by FEM FDM solution  FD |  −  FD | by FDM