Skip to main content

On a new numerical analysis of the Hall effect on MHD flow and heat transfer over an unsteady stretching permeable surface in the presence of thermal radiation and heat source/sink

Abstract

This paper employs a computational iterative approach known as the spectral local linearization method (SLLM) to analyze the Hall effect on MHD flow and heat transfer over an unsteady stretching permeable surface in the presence of thermal radiation and heat source/sink. To demonstrate the reliability of our proposed method, we made comparison with the Matlab bvp4c routine technique, and an excellent agreement was observed. The governing partial differential equations are transformed into a system of ordinary differential equations by using suitable similarity transformations. The results are obtained for velocity, temperature, skin friction and the Nusselt number.

MSC: 65PXX, 76-XX.

1 Introduction

Theoretical studies of magnetohydrodynamic flow and heat transfer over stretching surfaces have received great attention by virtue of their numerous applications in the fields of metallurgy and chemical engineering. Such applications include geothermal reservoirs, wire and fiber coating, food stuff processing, reactor fluidization, enhanced oil recovery, packed bed catalytic reactors and cooling of nuclear reactors. The primary aim in extrusion is to maintain the quality of the surface of the extricate. When the magnetic strength is strong, Hall currents cannot be neglected. In an ionized gas where the density is low and/or the magnetic field is very strong, the conductivity normal to the magnetic field is reduced due to the free spiraling of electrons and ions about the magnetic lines of force before suffering collisions, and a current is induced in a direction normal to both electric and magnetic fields. Due to Hall currents, the electrical conductivity of the fluid becomes anisotropic and this causes the secondary flow. The Hall effect is important when the Hall parameter, which is the ratio between the electron-cyclotron frequency and the electron-atom-collision frequency, is high. Hall currents are of great importance in many astrophysical problems, Hall accelerators, Hall sensors and flows of plasma in MHD power generators.

Examples of such studies include Sakiadis [1], [2] who did pioneering work on boundary layer flow on a continuously moving surface. Shateyi and Motsa [3] carried out a numerical analysis of the problem of magnetohydrodynamic boundary layer flow, heat and mass transfer rates on steady two-dimensional flow of an electrically conducting fluid over a stretching sheet embedded in a non-Darcy porous medium in the presence of thermal radiation and viscous dissipation. Shateyi [4] investigated thermal radiation and buoyancy effects on heat and mass transfer over a semi-infinite stretching surface with suction and blowing. Singh et al.[5] investigated two-dimensional unsteady flow of a viscous incompressible fluid about a stagnation point on a permeable stretching sheet. Shateyi and Motsa [6] numerically investigated the unsteady heat, mass and fluid transfer over a horizontal stretching sheet. More recently, Shateyi and Marewo [7] studied the magnetohydrodynamic boundary layer flow with heat and mass transfer of an UCM fluid over a stretching sheet in the presence of viscous dissipation and thermal radiation. Jhankal [8] considered the problem of unsteady boundary layer and heat transfer over a stretching surface in the presence of a transverse magnetic field.

Pal and Mondal [9] presented a model to study the effects of temperature-dependent viscosity and variable conductivity on mixed convective diffusion of species over a stretching sheet. When the conducting fluid is an ionized gas, and the strength of the applied magnetic field is large, the conductivity normal to the magnetic field is reduced due to the free spiraling of electrons and ions about the magnetic lines of force before collisions occur and a current is induced in a direction normal to both electric and magnetic fields. This phenomenon is known as the Hall effect or Hall current. The conductivity of the fluid is anisotropic and the effect of Hall current cannot be neglected when the medium is rarefied or if a strong magnetic field is present. The study of MHD viscous flows with Hall current has important applications in problems of Hall accelerators and sensors as well as flight magnetohydrodynamics. Pop and Watanabe [10] presented the problem of free convection flow of a conducting fluid which is permeated by a transverse magnetic field and the Hall effect is taken into account. Abd El-Aziz [11] investigated the effect of Hall currents on the flow and heat transfer of an electrically conducting fluid over an unsteady stretching surface in the presence of a strong magnetic field.

Shateyi and Motsa [12] numerically analyzed variable viscosity on magnetohydrodynamic fluid flow and heat transfer over an unsteady stretching surface with the Hall effect. Pal [13] studied the influence of Hall current and thermal radiation on flow and heat transfer characteristics in a viscous fluid over an unsteady stretching permeable surface. Zaman [14] considered the effects of Hall current on the flow of an incompressible, unsteady, viscous, MHD fluid with slip conditions. Lastly, Ali et al.[15] investigated heat and mass transfer of a steady flow of an incompressible electrically conducting fluid due to stretching plate under the influence of an applied uniform magnetic field and the effects of Hall current.

Governing equations modeling MHD flow and heat transfer over stretching surfaces are highly nonlinear, thereby exact solutions are impossible to obtain. Therefore, numerical solutions have always been developed and modified as a bid of getting more accurate and stable solutions. The current study seeks to investigate the Hall effect on MHD flow and heat transfer over an unsteady stretching permeable surface in the presence of thermal radiation and heat source/sink. We propose to numerically solve the present problem using a recently developed iterative method known as the spectral local linearization method (SLLM), Motsa [16]. The SLLM approach is based on transforming a nonlinear ordinary differential equation into an iterative scheme. The iterative scheme is then blended with the Chebyshev spectral method [17]. A similar approach to our current proposed method can be found in Motsa and Shateyi [18], [19] and Motsa et al.[20].

2 Mathematical formulation

We consider an unsteady two-dimensional laminar MHD boundary layer flow and heat transfer of an incompressible, viscous and electrically conducting fluid over a continuously moving stretching permeable surface. We take the x-axis along the stretching surface in the direction of the motion, and the y-axis is perpendicular to the sheet in the outward direction and the z-axis is transverse to the xy-plane. The flow is subjected to a transverse magnetic field of strength B 0 and the Hall current is taken into account in this study. The effects of the Hall current give rise to the Lorentz force in the z-direction, which induces a cross flow in this direction. We assume that there is no variation of flow and heat transfer quantities in the z-direction. The temperature at the sheet is maintained at a prescribed constant T w . The flow is generated due to a stretching sheet which is caused by simultaneously applying two equal but opposite forces. The velocity is assumed to be proportional to the distance from the origin. Under these assumptions, together with the usual boundary layer approximations, the dimensional governing equations of continuity, momentum and energy in the presence of a strong magnetic field, heat source and thermal radiation become as follows:

u x + v y =0,
(1)
u t +u u x +v u y =ν 2 u y 2 σ B 2 ρ ( 1 + m 2 ) (u+mw),
(2)
w t +u w x +v w y =ν 2 w y 2 σ B 2 ρ ( 1 + m 2 ) (muw),
(3)
T t +u T x +v T y =α 2 T y 2 1 ρ c p q r y + Q ρ c p (T T ).
(4)

The associated boundary conditions to the current problem are the following:

u = U w ( x , t ) , v = V w , T = T w ( x , t ) at  y = 0 , u 0 , w 0 , T T as  y ,
(5)

where u, v and w are the velocity components along x, y and z directions, respectively, and t is the time. T is the temperature within the fluid, c p is the specific heat at constant pressure, α is the thermal diffusivity, ν is the kinematic viscosity of the fluid density, T w (x,t) is the temperature on the stretching surface, T is the ambient temperature with T w > T . We have V w = ( ν U w / x ) 1 / 2 f(0) representing the mass transfer at the surface with V w >0 for injection and V w <0 for suction. We also have U w (x,t)=ax/(1ct), where a (stretching rate) and c are positive constants, with ct<1. It is noted that the stretching rate a/(1ct) increases with time since a>0. The surface temperature of the sheet varies with the distance x from the origin and time t and takes the form

T w (x,t)= T + b 2 x 2 ν ( 1 c t ) 3 / 2 ,
(6)

where b is a constant with b0.

We must remark that the particular forms of U w and T w (x,t) have been specifically chosen in order to come up with similarity transformations which then transform the governing partial differential equations into a set of nonlinear ordinary differential equations.

Following Ishak et al.[21], we introduce the following dimensionless functions f and θ, and the similarity variable η,

η = ( b ν ( 1 c t ) ) 1 2 y , ψ ( x , y , t ) = ( ν b 1 c t ) 1 2 x f ( η ) , T ( x , y , t ) = T + b 2 x 2 ν ( 1 c t ) 3 2 θ ( η ) , B 2 = B 0 2 ( 1 c t ) .
(7)

By using the Rosseland approximation, the radiative heat flux is given by

q r = 4 σ 3 K T 4 y ,
(8)

where σ and K are, respectively, the Stephan-Boltzman constant and the mean absorption coefficient. Assume that the temperature differences within the flow are small such that T 4 can be expressed as a linear function. Expanding T 4 in a Taylor series about T and neglecting higher order terms, we get

T 4 4 T T3 T 4 .
(9)

By using the above transformations, the governing partial differential equations are transformed into a system of non-dimensional nonlinear and coupled ordinary differential equations as follows:

f +f f f 2 A ( f + 1 2 η f ) M 1 + m 2 ( f + m g ) =0,
(10)
g +f g f gA ( g + 1 2 η g ) M 1 + m 2 ( m f g ) =0,
(11)
( 1 + 4 3 R ) θ +Pr ( f θ 2 f θ ) Pr A 2 ( 3 θ + η θ ) +δθ=0.
(12)

Here M 2 =σ B 0 2 /ρa is the magnetic parameter, A=c/a is a parameter that measures the unsteadiness, Pr=ν/α is the Prandtl number, R=4 σ T 3 /k K s is the thermal radiation parameter. The boundary conditions are as follows:

f(0)= f w , f (0)=1,θ(0)=1,g(0)=0,
(13)
f 0,θ0,g0as η,
(14)

where f(0)= f w with f w <0 or f w >0 corresponding to injection or suction, respectively. The physical engineering quantities of interest in this problem are the skin friction coefficients in the x- and z-directions and the local Nusselt number, Nu x which are defined as follows:

C f x = 2 μ ( u / y ) y = 0 ρ U w 2 = 2 Re x 1 / 2 f ( 0 ) , C f z = 2 μ ( w / y ) y = 0 ρ U w 2 = 2 Re x 1 / 2 g ( 0 ) , Nu x = x q w κ ( T w T ) ,
(15)

where τ w =μ ( u y ) y = 0 is the wall shear stress, and q w =κ ( T y ) y = 0 is the surface heat flux, where μ and κ are the dynamic viscosity and thermal conductivity, respectively.

3 Method of solution

3.1 Basic idea of the spectral local linearization method (SLLM)

Consider the problem of finding Z=[ Z 1 (η), Z 2 (η),, Z m (η)] that satisfies the system

L i + N i = H i ,
(16)

of m differential equations, where i=1,2,,m, and each H i is a function of η[a,b]. Also, L i and N i are linear and nonlinear components, respectively.

Basically, the SLLM is an iterative method for solving differential equations such as (16) which begins with an initial approximation Z 0 of Z. Successive application of the SLLM generates approximations Z 1 , Z 2 , , where Z r =[ Z 1 , r ,, Z m , r ] for each r=0,1,2, . Upon linearizing nonlinear component N i , differential equation (16) shall be solved numerically using the Chebyshev spectral collocation method [17].

First, let

L i | r + 1 + N i | r + 1 = H i
(17)

denote the ith differential equation (16) after the first r+1 iterations of the SLLM. Let w r be an n-tuple of Z i , r and its derivatives. If we assume that N i is a function of w r only, then linearization of N i at w r is

N i | r + 1 N i | r + 1 + N i | r ( w r + 1 w r ),
(18)

which, upon inserting into equation (17) and re-arranging the result, gives

L | r + 1 + N i | r w r + 1 = H i + N i | r w r N i | r .
(19)

Equation (19) shall be solved using the Chebyshev spectral collocation method. For the sake of brevity, we leave out the details of this method. We urge the interested reader to see, for example, [22].

3.2 Solving current problem using the SLLM

If we let f =p, then equations (10)-(12) become

p +f p p 2 A ( p + η 2 p ) M 1 + m 2 (p+mg)=0,
(20)
g +f g pgA ( g + η 2 g ) + M 1 + m 2 (mpg)=0,
(21)
( 1 + 4 3 R ) θ +Pr ( f θ 2 p θ ) Pr A 2 ( 3 θ + η θ ) +δθ=0.
(22)

Equations (20)-(22) together with the change of variable f =p may be written as

L 1 + N 1 = H 1 ,
(23)
L 2 + N 2 = H 2 ,
(24)
L 3 + N 3 = H 3 ,
(25)
L 4 + N 4 = H 4 ,
(26)

where

L 1 = f p, N 1 =0, H 1 =0,
(27)
L 2 = p ( A + M 1 + m 2 ) p A η 2 p M m 1 + m 2 g , N 2 ( p , p ) = f p p 2 , H 2 = 0 ,
(28)
L 3 = g ( A + M 1 + m 2 ) g A η 2 g + M m 1 + m 2 p , N 3 ( g , g ) = f g p g , H 3 = 0 ,
(29)
L 4 = ( 1 + 4 3 R ) θ Pr A η 2 θ + ( δ 3 2 Pr A ) θ , N 4 ( θ , θ ) = Pr ( f θ 2 p θ ) , H 4 = 0 .
(30)

As an illustration we derive equation (19) when i=2. In this case, equation (19) is

L 2 | r + 1 + N 2 p | r + 1 p r + 1 + N 2 p | r + 1 p r + 1 = H 2 + N 2 p | r + 1 p r + N 2 p | r + 1 p r N 2 | r ,
(31)

which, upon substituting (28) and simplifying, becomes

p r + 1 ( A + M 1 + m 2 ) p r + 1 A η 2 p r + 1 2 p r p r + 1 + f r p r + 1 M m 1 + m 2 g r + 1 = p r 2 .
(32)

A similar approach for each i=1,3,4 yields

f r + 1 = p r + 1 ,
(33)
g r + 1 A η 2 g r + 1 ( A + M 1 + m 2 ) g r + 1 p r g r + 1 + f r g r + 1 + M m 1 + m 2 p r + 1 =0,
(34)
( 1 + 4 3 R ) θ r + 1 PrA η 2 θ r + 1 + ( δ 3 2 Pr A ) θ r + 1 2Pr p r θ r + 1 +Pr f r θ r + 1 =0.
(35)

Equations (32)-(35) are subject to the boundary conditions

p r + 1 (0)=1, p r + 1 ()=0,
(36)
f r + 1 (0)= f w ,
(37)
g r + 1 (0)=0, g r + 1 ()=0,
(38)
θ r + 1 (0)=1, θ r + 1 ()=0,
(39)

respectively. Chebyshev differentiation [17] replaces equations (32) through (35) with the discrete form

( D I O O O A 11 A 12 O O A 21 A 22 O O O O A 3 )( f r + 1 p r + 1 g r + 1 Θ r + 1 )=( 0 B 1 B 2 B 3 )
(40)

or equivalently

D f r + 1 = p r + 1 ,
(41)
( A 11 A 12 A 21 A 22 )( p r + 1 g r + 1 )=( B 1 B 2 ),
(42)
A 3 Θ r + 1 = B 3
(43)

subject to the boundary conditions

f r + 1 ( τ N ) = f w , p r + 1 ( τ N ) = 1 , p r + 1 ( τ 0 ) = 0 , g r + 1 ( τ N ) = 0 , g r + 1 ( τ 0 ) = 0 , θ r + 1 ( τ N ) = 1 , θ r + 1 ( τ 0 ) = 0 ,

where

A 11 = D 2 + diag { f r A η 2 } D 2 diag { p r } ( A + M 1 + m 2 ) I , A 12 = M m 1 + m 2 I , B 1 = p r 2 , A 22 = D 2 + diag { f r A η 2 } D diag { p r } ( A + M 1 + m 2 ) I , A 21 = M m 1 + m 2 I , B 2 = 0 , A 3 = ( 1 + 4 3 R ) D 2 + Pr diag { f r A η 2 } D 2 Pr diag { p r } + ( δ 3 2 Pr A ) I , B 3 = 0

and O, 0 are zero matrices of order (N+1)×(N+1), (N+1)×1, respectively. Before we solve linear system (40), we need to apply boundary conditions. This is easier to illustrate using the equivalent form, equations (41)-(43), as follows:

( D 0 0 1 ) ( f r + 1 ( τ 0 ) f r + 1 ( τ N ) ) = ( p r + 1 f w ) , ( 1 0 0 0 0 0 A 11 A 12 0 0 1 0 0 0 0 0 0 1 0 0 A 21 A 22 0 0 0 0 0 1 ) ( p r + 1 ( τ 0 ) p r + 1 ( τ N ) g r + 1 ( τ 0 ) g r + 1 ( τ N ) ) = ( 0 B 1 1 0 B 2 0 ) , ( 1 0 0 A 3 0 0 1 ) ( θ r + 1 ( τ 0 ) θ r + 1 ( τ N ) ) = ( 0 B 3 1 ) .

We choose initial approximations

f 0 ( η ) = f w + 1 e η , p 0 ( η ) = θ 0 ( η ) = e η , g 0 ( η ) = η e η

that satisfy boundary conditions (36)-(39) so that the SLLM generates subsequent approximations f r , p r , g r , θ r for each r=1,2, .

4 Results and discussion

The numerical results iteratively generated by the SLLM for the main parameters that have significant effects on the flow properties are presented in this section. All the SLLM results presented in this work were obtained using N=50 collocation points, and we are glad to highlight that convergence was achieved in as few as six iterations. We take the infinity value η to be 40. The magnetic field is taken quite strong by assigning large values of M to ensure the generation of Hall currents. In order to validate our numerical method, it was compared to MATLAB routine bvp4c which is an adaptive Lobatto quadrature iterative scheme.

In Table 1 we present a comparison between the SLLM approximate results and the bvp4c results for selected default values of the stretching parameter A. It can be clearly seen from this table that there is an excellent agreement between the results from the two methods. Also Table 1 shows that an increase in the unsteadiness parameter leads to an increase in the skin-friction coefficients in both directions. Also the heat transfer gradient increases as the values of the unsteadiness parameter increase. The negative values of f (0) mean that the solid surface exerts a drag force on the fluid. This is due to the development of the velocity boundary layer which in the current study is caused solely by the stretching sheet.

Table 1 Comparison of the SLLM results of f (0), g (0), θ (0) with those obtained by bvp4c for different values of the unsteadiness parameter

In Table 2 we display the effect of the magnetic parameter on the skin friction coefficients and the Nusselt number. The magnetic parameter M represents the significance of the magnetic field on the flow properties. As the magnetic strength increases, the dragging effect is clearly seen by the significant increments in the skin friction. We also observe that increasing the values of the Hartman number leads to the lowering of the values of the Nusselt number. Application of a strong magnetic field reduces the velocity, which in turn increases heat diffusion within the fluid flow. This physically explains why heat transfer at the wall is reduced as M is increased.

Table 2 Comparison of the SLLM results of f (0), g (0), θ (0) with those obtained by bvp4c for different values of the magnetic parameter

Table 3 displays the influence of the Hall current on the skin friction coefficients as well as the Nusselt number. The skin friction coefficient is reduced as the values of the Hall current parameter increase. This explains why the skin friction coefficient in the axial direction decreases. However, in the transverse direction the skin friction increases as the Hall current increases. There is a small effect of the Hall current on heat transfer rate on the stretching surface.

Table 3 Comparison of the SLLM results of f (0), g (0), θ (0) with those obtained by bvp4c for different values of the Hall parameter

The effect of the unsteadiness parameter A on the axial velocity f (η) is presented in Figure 1. It can be observed from this figure that the axial velocity profiles decrease as A increases. Increasing the unsteadiness parameter A causes the velocity boundary layer thickness to decrease. This in turn may cause delays in transition from laminar to turbulent flows.

Figure 1
figure 1

Graph of the SLLM solutions for the horizontal velocity for different values of A , with M=1 , Pr=0.71 , δ=0.5 , R=1 , m=0.5 , f w =1 .

The influence of the suction/injection parameter f w on the axial velocity is depicted in Figure 2. The axial velocity is significantly influenced by this parameter. The velocity boundary layer is greatly enhanced when fluid is injected ( f w <0) into the flow system thereby increasing the velocity profiles. However, removing fluid from the flow system through suction, as expected, drastically reduces the velocity profiles as can be clearly seen in Figure 2.

Figure 2
figure 2

Graph of the SLLM solutions for the horizontal velocity for different values of f w , with M=1 , Pr=0.71 , δ=0.5 , R=1 , m=0.5 , A=0.5 .

The application of a magnetic field perpendicular to the flow produces a drag force known as the Lorentz force. This force reduces the axial velocity as can be observed in Figure 3. Thus increasing the magnetic strength parameter M reduces the velocity boundary layer and thereby reduces the velocity profiles.

Figure 3
figure 3

Graph of the SLLM solutions for the horizontal velocity for different values of the magnetic parameter M , with A=0.5 , Pr=0.71 , δ=0.5 , R=1 , m=0.5 , f w =1 .

The effect of the Hall current parameter m on the axial velocity is shown in Figure 4. The velocity is enhanced as the values of m increase. However, the axial velocity profiles approach their classical values when the Hall current parameter m becomes large (m>1.5) in our current study. Any further increase of the Hall current would make the magnetic effect insignificant.

Figure 4
figure 4

Graph of the SLLM solutions for the horizontal velocity for different values of the Hall parameter m , with M=1 , Pr=0.71 , δ=0.5 , R=1 , A=0.5 , f w =1 .

Figure 5 displays the effect of the unsteadiness parameter on the transverse velocity profiles. This parameter has very significant influence on the transverse velocity. It can clearly be seen from this figure that the effect of increasing the unsteadiness parameter A is to greatly reduce the transverse velocity g(η) near the plate. In Figure 6, we display the influence of the injection/suction parameter on the transverse velocity. Injecting/blowing fluid into the flow system causes the transverse velocity profiles to significantly increase, whereas suction, as expected, reduces the transverse velocity distribution.

Figure 5
figure 5

Graph of the SLLM solutions for the transverse velocity for different values of the unsteadiness parameter A , with M=1 , Pr=0.71 , δ=0.5 , R=1 , m=0.5 , f w =1 .

Figure 6
figure 6

Graph of the SLLM solutions for the transverse velocity for different values of f w , with M=1 , Pr=0.71 , δ=0.5 , R=1 , m=0.5 , A=0.5 .

The influence of the magnetic parameter M on the transverse velocity is depicted in Figure 7. One can clearly observe in this figure that the magnetic parameter M has a very significant effect on the transverse velocity profiles. An increase in the values of M leads to a huge increase in the values of the transverse velocity due to the Lorentz force which acts in the direction hence tends to accelerate the fluid flow in the traverse direction.

Figure 7
figure 7

Graph of the SLLM solutions for the transverse velocity for different values of the magnetic parameter M , with A=0.5 , Pr=0.71 , δ=0.5 , R=1 , m=0.5 , f w =1 .

The effect of the Hall current parameter on the transverse velocity is displayed in Figure 8. Increasing the values of m causes the transverse velocity to rapidly increase. Figure 9 shows the effect of the unsteadiness parameter A on the temperature profiles. Increasing the values of A reduces the thermal boundary layer thickness thus reducing the fluid temperature distribution.

Figure 8
figure 8

Graph of the SLLM solutions for the transverse velocity for different values of the Hall parameter m , with M=1 , Pr=0.71 , δ=0.5 , R=1 , A=0.5 , f w =1 .

Figure 9
figure 9

Graph of the SLLM solutions of the temperature profiles for different values of A , with M=1 , Pr=0.71 , δ=0.5 , R=1 , m=0.5 , f w =1 .

In Figure 10 we have the effect of the heat source/sink parameter δ on the temperature profiles. As expected, it is observed in this figure that the temperature in the boundary layer increases with increasing values of δ. The heat absorption due to a uniform sink (δ<0) leads to the reduction of the thermal boundary layer thickness, whereas this layer increases significantly with increases in δ>0.

Figure 10
figure 10

Graph of the SLLM solutions of the temperature profiles for different values of the heat source/sink parameter δ , with M=1 , Pr=0.71 , A=0.5 , R=1 , m=0.5 , f w =1 .

Figure 11 shows the effect of injection/suction on the temperature profiles. Injecting fluid into the flow system causes enhancement of the thermal boundary layer thickness. This in turn causes the temperature distribution in the flow system to greatly increase. On the other hand, suction reduces the thermal boundary layer thickness.

Figure 11
figure 11

Graph of the SLLM solutions of the temperature profiles for different values of the suction/injection parameter f w , with M=1 , Pr=0.71 , δ=0.5 , R=1 , m=0.5 , A=0.5 .

Figure 12 shows that the thermal boundary layer becomes thicker when the values of the magnetic parameter are increased. Application of the magnetic field gives rise to a resistive-type force which slows down the motion of the fluid and in turn increases the thermal boundary layer, hence increasing the temperature of the fluid.

Figure 12
figure 12

Graph of the SLLM solutions of the temperature profiles for different values of the magnetic parameter M , with A=0.5 , Pr=0.71 , δ=0.5 , R=1 , m=0.5 , f w =1 .

Figure 13 is plotted to depict the influence of the thermal radiation parameter R on the temperature profiles. We clearly observe that the temperature in the boundary layer increases with increasing values of the thermal radiation parameter. This is due to the fact that the divergence of the radiative heat flux increases as the Rosseland radiative absorption K decreases, which in turn increases the rate of radiative heat transfer to the fluid. Thus the presence of thermal radiation enhances thermal state of the fluid causing its temperature to significantly increase.

Figure 13
figure 13

Graph of the SLLM solutions of the temperature profiles for different values of the thermal radiation parameter R , with M=1 , Pr=0.71 , δ=0.5 , A=0.5 , m=0.5 , f w =1 .

5 Conclusion

The present work analyzed MHD unsteady flow and heat transfer of an electrically conducting fluid over a stretching sheet in the presence of thermal radiation and the Hall effect. Much attention was given to trying to investigate how the velocity field, skin friction, temperature distribution and heat transfer are influenced by the parameters of importance in this study. The governing partial differential equations are transformed into a system of nonlinear ordinary differential equations by using suitable similarity variables. The resultant system of nonlinear ordinary differential equations is solved numerically by a recently developed technique known as the spectral local linearization method. The accuracy of the SLLM is validated against the MATLAB in-built bvp4c routine for solving boundary value problems. The following conclusions were drawn in our investigation.

An excellent agreement was observed between our results and those obtained using the bvp4c routine technique giving confidence to our present results.

The unsteadiness parameter A has significant effects on the velocity components and temperature profiles. The maximum axial velocity, transverse velocity and temperature profiles are attained when the flow is steady (A=0).

Increasing the values of the magnetic field strength decreases the momentum boundary layer thickness while increasing the thermal boundary layer thickness.

The velocity components are enhanced as the Hall parameter increases.

The fluid temperature increases with increasing values of thermal radiation as well as a heat source.

The heat transfer rate and the skin friction coefficient in the x-direction are increased while the skin friction in the z-direction decreases as the unsteadiness parameter increases.

The skin friction coefficients are enhanced while the heat transfer rate is depressed by increasing the values of magnetic strengths.

References

  1. Sakiadis BC: Boundary-layer behaviour on continuous solid surfaces: I. Boundary-layer equations for two-dimensional and axisymmetric flow. AIChE J. 1961, 7: 26-28. 10.1002/aic.690070108

    Article  Google Scholar 

  2. Sakiadis BC: Boundary-layer behaviour on continuous solid surfaces: II. Boundary-layer equations for two-dimensional and axisymmetric flow. AIChE J. 1961, 7: 221-225. 10.1002/aic.690070211

    Article  Google Scholar 

  3. Shateyi S, Motsa SS: Hydromagnetic non-Darcy flow, heat and mass transfer over a stretching sheet in the presence of thermal radiation and ohmic dissipation. Can. J. Chem. Eng. 2011, 89: 1388-1400. 10.1002/cjce.20499

    Article  Google Scholar 

  4. Shateyi S: Thermal radiation and buoyancy effects on heat and mass transfer over a semi-infinite stretching surface with suction and blowing. J. Appl. Math. 2008., 2008:

    Google Scholar 

  5. Singh P, Jangid A, Tomer NS, Sinha D: Effects of thermal radiation and magnetic field on unsteady stretching permeable sheet in presence of free stream velocity. Int. J. Inf. Math. Sci. 2010, 6(3):160-166.

    Google Scholar 

  6. Shateyi S, Motsa SS: Thermal radiation effects on heat and mass transfer over an unsteady stretching surface. Math. Probl. Eng. 2009., 2009: 10.1155/2009/965603

    Google Scholar 

  7. Shateyi S, Marewo GT: A new numerical approach of MHD flow, heat and mass transfer for the UCM fluid over a stretching surface in the presence of thermal radiation. Math. Probl. Eng. 2013., 2013:

    Google Scholar 

  8. Jhankal AK: Unsteady MHD boundary layer flow and heat transfer over a stretching surface. Int. J. Phys. Math. Sci. 2013, 3(2):110-114.

    Google Scholar 

  9. Pal D, Mondal H: Effects of temperature-dependent viscosity and variable thermal conductivity on MHD non-Darcy mixed convective diffusion of species over a stretching sheet. J. Egypt. Math. Soc. 2013.

    Google Scholar 

  10. Pop I, Watanabe T: Hall effect on magnetohydrodynamic free convection about a semi-infinite vertical plate. Int. J. Eng. Sci. 1994, 32: 1903-1911. 10.1016/0020-7225(94)90087-6

    Article  Google Scholar 

  11. Abd El-Aziz M: Flow and heat transfer over an unsteady stretching surface with Hall effect. Meccanica 2010, 45: 97-109. 10.1007/s11012-009-9227-x

    Article  MathSciNet  Google Scholar 

  12. Shateyi S, Motsa SS: Variable viscosity on magnetohydrodynamic fluid flow and heat transfer over an unsteady stretching surface with Hall effect. Bound. Value Probl. 2010., 2010: 10.1155/2010/257568

    Google Scholar 

  13. Pal D: Hall current and MHD effects on heat transfer over an unsteady stretching permeable surface with thermal radiation. Comput. Math. Appl. 2013, 63(7):1161-1180. 10.1016/j.camwa.2013.07.010

    Article  Google Scholar 

  14. Zaman H: Hall effects on the unsteady incompressible MHD fluid flow with slip conditions and porous walls. Appl. Math. Phys. 2013, 1(3):31-38.

    Google Scholar 

  15. Ali M, Alam MS, Asaduzzaman : Hall effects on steady MHD boundary layer flow and heat transfer due to stretching plate in the presence of heat source or sink. Int. J. Sci. Eng. Res. 2013, 4(8):2178-2183.

    Google Scholar 

  16. Motsa SS: A new spectral local linearization method for nonlinear boundary layer flow problems. J. Appl. Math. 2013., 2013:

    Google Scholar 

  17. Trefethen LN: Spectral Methods in MATLAB. SIAM, Philadelphia; 2000.

    Book  Google Scholar 

  18. Motsa SS, Shateyi S: A successive linearization method approach to solving Lane-Emden type of equations. Math. Probl. Eng. 2012., 2012:

    Google Scholar 

  19. Motsa SS, Shateyi S: Successive linearization analysis of the effects of partial slip, thermal diffusion, and diffusion-thermo on steady MHD convective flow due to a rotating disk. Math. Probl. Eng. 2012., 2012:

    Google Scholar 

  20. Motsa SS, Makinde OD, Shateyi S: Application of successive linearisation method to squeezing flow with bifurcation. Adv. Math. Phys. 2014., 2014:

    Google Scholar 

  21. Ishak A, Nazar R, Pop I: Heat transfer over an unsteady stretching permeable surface with prescribed wall temperature. Nonlinear Anal., Real World Appl. 2009, 10: 2909-2913. 10.1016/j.nonrwa.2008.09.010

    Article  MathSciNet  Google Scholar 

  22. Shateyi S, Marewo GT: On a new numerical approach for the laminar boundary layer flow and heat transfer along a stretching cylinder embedded in a porous medium with variable thermal conductivity. J. Appl. Math. 2013., 2013: 10.1155/2013/576453

    Google Scholar 

Download references

Acknowledgements

The authors wish to acknowledge financial support from the University of Venda and NRF. The authors are grateful to the reviewers for their constructive suggestions.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Stanford Shateyi.

Additional information

Competing interests

The authors declare that they do not have any conflict of interests in their submitted paper.

Authors’ contributions

SS wrote the introduction, formulated the problem, generated and discussed the results and GTM wrote the code, generated and also discussed the results and proof read the manuscript. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0), which permits use, duplication, adaptation, distribution, and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Shateyi, S., Marewo, G.T. On a new numerical analysis of the Hall effect on MHD flow and heat transfer over an unsteady stretching permeable surface in the presence of thermal radiation and heat source/sink. Bound Value Probl 2014, 170 (2014). https://doi.org/10.1186/s13661-014-0170-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13661-014-0170-y

Keywords