Abstract

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 published work using the finite difference method (FDM). Numerical examples show that the finite element method (FEM) gives more accurate results in comparison with the finite difference method (FDM).

1. 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.

2. 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 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 iswhere is the velocity vector, is the fluid density, is the pressure, is the current density, 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, 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

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) becomesubject 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 ; thus (13) becomessubject to the initial and boundary conditions:where .

3. Finite Element Method

3.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 and integrate over the domainwhere

Integrating by parts, we obtainwhere denotes the -inner product and for .

We shall approximate the solution of (21) by assuming that and lie in finite dimensional subspace of for each . Let for and assume that the set is linearly independent. Further, let be a partition of the interval into subintervals . Now we define the finite dimensional space spanned by asTo this end, the approximate solution isInserting (23) into (21) and selecting as trial function the basis function of , we obtain a system of ODEs:where and . Here is the Stiffness matrix and is the Mas matrix defined, respectively, aswith being the usual finite element basis corresponding to the partition . Thus, to compute the entries of the Stiffness matrix , next, we need to determine .

Hereand thenwhere .

Using (26) and (27), with a uniform mesh , we get

Scheme (24) is called semidiscretization, since is still a continuous function of [17].

3.2. Time Stepping

In this section, we consider the semidiscretization in time. We first discretize the time interval into a uniform grid with size . Approximating the derivative in (24) at time level by the Crank-Nicolson scheme with , 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:

4. 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 some examples. These tests are chosen such that there exist analytical solutions for them to give an obvious overview of the methods presented in this work.

Numerical Example  1. As an application to the FEM and FDM, we consider the following test case with the values , , , , fixed, and . Next the computed matrices , and given in (28), (29), and (30), respectively, are used in (32) to obtain the velocity .
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 ,  ,  ,  ,   fixed, and . 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.

5. 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.

Conflicts of Interest

The authors declare that they have no conflicts of interest.