Effects of Hall Current on Flow of Unsteady MHD Axisymmetric Second-Grade Fluid with Suction and Blowing over an Exponentially Stretching Sheet ()
1. Introduction
The theoretical study of boundary layer flows induced by a stretching sheet [1] is of considerable interest because their applications in fibers spinning, manufacturing of plastic and rubber sheets, the aerodynamic extrusion of plastic sheets, hot rolling and cooling of an infinite metallic plate in a cooling bath. Partha et al. [2] studied the effect of viscous dissipation on the mixed convection heat transfer from an exponentially stretching surface. Sajid and Hayat [3] extended this problem by investigating the radiation effects on the flow over an exponentially stretching sheet, and solved the problem analytically using the homotopy analysis method. The numerical solution for the same problem was then given by Bidin and Nazar [4] . Recently Ishak [5] investigated the thermal radiation effects on hydro-magnetic flow due to an exponentially stretching sheet. Magyari and Keller [6] investigated the steady boundary layers on an exponentially stretching continuous surface with an exponential temperature distribution. In recent years the theoretical study of Effects of Hall current on MHD flows has been a subject of great interest due to its widely spread applications in power generators and pumps, Hall accelerators, refrigeration coils, electric transformers, in flight MHD, solar physics involved in the sunspot development, the solar cycle, the structure of magnetic stars, electronic system cooling, cool combustors, fiber and granular insulation, oil extraction, thermal energy storage and flow through filtering devices and porous material regenerative heat exchangers. Usually Hall term representing the Hall current was ignored in applying Ohm’s law, because it has no remarkable effect for small and moderate values of the magnetic field. The effects of Hall current are very important if the strong magnetic field is applied [7] , because for strong magnetic field electromagnetic force is noticeable. The recent investigation for the applications of MHD is towards a strong magnetic field, due to which study of Hall current is important. In presence of a strong magnetic field in an ionized gas of low density, the conductivity normal to the magnetic field is decreased by free spiraling of electrons and ions about the magnetic lines of force before suffering collisions. A current induced in a direction normal to the electric and magnetic fields is called Hall current [8] . Some interesting old and new studies regarding the effects of Hall current on MHD flow are done by [9] -[17] .
The present investigation is to analyze the effects of Hall current on flow induced by an exponentially stretching sheet of unsteady MHD axisymmetric second-grade fluid with suction and blowing. Here, we assume that sheet is stretching exponentially with the radius. The existing literature on the topic shows that effects of Hall current on an axisymmetric flow over an exponentially stretching sheet in cylindrical coordinates have not been investigated so far. The arising non-linear problem is solved by the homotopy analysis method (HAM), which is a novel technique and has been used by many researchers [18] -[26] .
2. Mathematical Formulation of the Problem and Analytic Solution
We consider the unsteady flow of an electrically conducting incompressible second-grade fluid over a porous exponentially stretching sheet placed in the plane in the presence of transverse magnetic field. The cylindrical coordinates are used and it is assumed that the flow takes place in the upper half plane. We assume rotational symmetry of the flow so all physical quantities are independent of θ i.e., the azimuthal component of velocity vanishes identically. It is assumed that the sheet is being stretched with an exponential velocity, where is the reference velocity and is the reference radius. The Cauchy stress tensor for a second-grade fluid is given as [27]
(1)
where is the scalar pressure, is the identity tensor, is the coefficient of viscosity, are the material parameters of second-grade fluid, and are the first two Rivlin-Ericksen tensors defined by [28]
(2)
(3)
It is assumed that the flow meets the Clausius-Duhem inequality and that the specific Helmholtz free energy of the fluid is minimum at equilibrium [29] when
(4)
For detailed analysis about the signs of these normal stress moduli one may see [30] . The equations governing the magnetohydrodynamic flow with Hall effects are:
Velocity field:
(5)
Continuity equation:
(6)
Equation of motion:
(7)
Equations for the stream function:
(8)
Maxwell equations:
(9)
Generalized Ohm’s law:
(10)
where is the radial velocity and is the axial velocity, is time, is the total magnetic field, is the applied magnetic field, is the induced magnetic field, is the current density, is the electrical conductivity of the fluid, is the electric field, is the magnetic permeability, is the fluid density, is the stream function, and are the cyclotron frequency and collision time of the electrons respectively. We assume that, the quantities and are constants throughout the flow field, the magnetic field is normal to the velocity vector and the induced magnetic field is neglected compared with the imposed magnetic field, so that the magnetic Reynolds number is small [31] . We also suppose that and (where and are the cyclotron frequency and collision time for ions respectively). The radial velocity is zero far from the sheet and the pressure is uniform, so we neglect the pressure gradient term [32] . Under aforementioned assumptions and using the boundary layer approximations [33] , the unsteady problem with Hall currents become
(11)
(12)
The boundary conditions applicable to the flow are
(13)
where is the kinematic viscosity, is the Hall parameter and is the suction or blowing velocity in the, the positive sign denotes suction while negative sign denotes blowing. In order to nondimensionalize the problem let us introduce the similarity transformations
(14)
where is the dimensionless stream function. Equation (11) is identically satisfied and Equations (12) and (13) become [16]
(15)
(16)
where prime denotes differentiation with respect to, is the dimensionless modified Hartmann number [31] , is a dimensionless coordinate, is the dimensionless second grade fluid parameter and is the dimensionless suction or blowing parameter,
corresponds to suction and corresponds to blowing. The local skin friction coefficient or fractional drag coefficient on the surface of the exponentially stretching sheet is
(17)
Now using Equations (1), (2), (3) and (14), Equation (17) can be written in dimensionless variables as
(18)
where is the local Reynolds number and is a dimensionless constant. Note that Equation (18) is obtained by taking, for it is too difficult to find it.
To start with the homotopy analysis method, due to the boundary conditions (16) it is reasonable to choose the initial guess approximation and the auxiliary linear operator
(19)
(20)
where satisfy
(21)
where, and are arbitrary constants. Following the HAM and trying higher iterations with the unique and proper assignment of the results converge to the exact solution:
(22)
Using the symbolic computation software such as MATLAB, MAPLE, MATHEMATICA to successively obtain
(23)
similarly, , and so on are calculated. The obtained values of, , , - - - -leads us to take
(24)
After very lengthy calculations we arrive at the total complex analytic solution in compact form as
(25)
where from initial guess in Equation (19) we obtain
(26)
all other unknown constants can be determined by utilizing first four given in Equation (26) by using the recurrence relations, which we calculated but it is not possible to write here due to their length.
3. Graphs, Tables and Discussion
Here we discuss the convergence of the analytic solution and effect of emerging parameters on the radial and axial velocities of the fluid. The auxiliary parameter, gives the convergence region and rate of approximation for the homotopy analysis method. The is plotted for real part of axial velocity and from Figure 1 we observe that. Our calculations depict that the series of the dimensionless stream function in Equation (26) converges in the whole region of and for. Figure 2 indicates the variation of the real part of the axial velocity with η for suction. This figure shows that in case of suction for fixed values of, N, , , and s with increase in dimensionless time ξ real part of the axial velocity decreases. In Figure 2 boundary layer structure is observed and the boundary layer thickness increases with increasing time, which results in thickening of the boundary layer, which is used in recognition of the
fluid. Figure 3 elucidates that in case of blowing with the increase in Hall current imaginary part of the axial velocity decreases in magnitude and boundary layer thickness increases with increasing Hall current. Figure 4 illustrates that with the increase in second-grade parameter real part of the axial velocity increases for blowing. Figure 5 describes that for suction with the increase in exponential stretching parameter imaginary part of the axial velocity decreases in magnitude and boundary layer thickness increases with increasing and shear thickening effects are seen. Figure 6 represents that in case of suction with increase in dimensionless time real part of the radial velocity decreases. Clearly the viscosity induces in the fluid with the passage of time and after some time the steady state is achieved which proves the uniform validity of the solution for all time in the whole spatial region. Figure 7 shows that with the increase in second-grade parameter real part of the radial velocity increases for blowing. Figure 8 describes that for blowing with the increase in exponential stretching parameter real part of the radial velocity increases.
Tables 1 and 2 are prepared for the variation of the absolute values of the shear stress at the wall. It is observed from Table 1 that for suction and blowing for fixed values of N, , , and with increase in second-grade fluid parameter absolute values of the shear stress at the wall decreases. It is also observed from Table 1 that with increase in Hall current absolute values of the decreases for suction and blowing. When magnetic field is applied normal to the fluid velocity then it gives rise to a drag-like or resistive force which slow down or suppress the motion of the fluid on the stretching surface. This leads to a re-
duction in the velocity of the fluid and flow rates. With the increase in the strength of the magnetic field the motion of the particulate suspension on the surface reduces due to which shear stress at the wall reduces with increase in and, as observed in Table 1. From Table 1 we note that with increase in dimensionless time the absolute values of the increases.
Table 2 shows that for fixed values of, , , and with increase in exponential stretching absolute value of the decreases for suction and blowing. Also with increase in Hall parameter the absolute values of the decreases for suction and blowing.
4. Conclusion
Effects of Hall current on flow of unsteady MHD axisymmetric second-grade fluid with suction and blowing over an exponentially stretching sheet are seen first time. The present complex explicit analytic solution is uniformly valid for all dimensionless time in the entire flow regime. Convergence of the solution is appropriately discussed. Graphical and tabular results for Hall parameter, second-grade parameter, exponentially stretching parameter and dimensionless time reveal that Hall parameter, second-grade parameter, exponential stretching and dimensionless time have a significant influence on the radial and axial components of velocity.