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.