^{1}

^{2}

Mixed convective flow of a viscous incompressible electrically conducting micropolar fluid along a semi-infinite inclined permeable flat plate with viscous dissipation has been analyzed numerically. With appropriate transformations the boundary layer equations are transformed into a set of nonlinear ordinary differential equations. The local similarity solutions of the transformed dimensionless equations for the flow, microrotation and the heat transfer characteristics are evaluated using Nachtsheim-Swigert shooting iteration technique (guessing the missing value) together with sixth order Runge-Kutta-Butcher integration scheme. Numerical results are presented in the form of non-dimensional velocity, microrotation and temperature profiles within the boundary layer for different parameters entering into the analysis. The effects of pertinent parameters on the local skin friction coefficient (viscous drag), plate couple stress and rate of heat transfer (Nusselt number) are also displayed graphically.

Micropolar fluids are those which contain micro-constituents that can undergo rotation, the presence of which can affect the hydrodynamics of the flow so that it can be distinctly non Newtonian. These fluids are fluids with microstructure belonging to a class of complex fluids with nonsymmetrical stress tensor referred to as micromorphic fluids. It has many practical applications, for examples analyzing the behavior of exotic lubricants, the flow of colloidal suspensions or polymeric fluids, liquid crystals, additive suspensions, animal blood, body fluids, and turbulent shear flows.

Convective flows over a flat plate that is immersed in a micropolar fluid has attracted an increasing amount of attention since the early studies of Eringen [

The study of magnetohydrodynamic (MHD) flow of an electrically conducting fluid is of considerable interest in modern metallurgical and metal-working processes. There has been a great interest in the study of Magnetohydrodynamic flow and heat transfer in any medium due to the effect of magnetic field on the boundary layer flow control and on the performance of many systems using electrically conducting fluids. This type of flow has attracted the interest of many researchers due to its applications in many engineering problems such as MHD generators, plasma studies, nuclear reactors, geothermal energy extractions. By the application of magnetic field, hydromagnetic techniques are used for the purification of molten metals from non-metallic inclusions. So such type of problem that we are dealing with is very much useful to polymer technology and metallurgy. Rahman and Sattar [

The aim of the present work is to study the mixed convective flow of micropolar fluids past an inclined porous plate with viscous dissipation effect. The resulting governing equations are solved by Nachtsheim-Sweigert shooting iteration technique together with Runge-Kutta-Butcher initial value solver and reported graphically.

A two-dimensional steady MHD laminar mixed convective flow of a viscous incompressible electrically conducting micropolar fluid along a semi-infinite inclined permeable flat plate with an acute angle α to the vertical is considered. With x-axis measured along the plate, a magnetic field of uniform strength B_{0} is applied in the y-direction which is normal to the flow direction. Fluid injection or suction is imposed at the plate surface. The external flow takes place in the direction parallel to the inclined plate and has the uniform velocity U ∞ . The temperature of the surface is held uniform at T_{w} which is higher than the ambient temperature T ∞ . The flow configuration and coordinate system are shown in the following

We further assume that: 1) the fluid has constant kinematic viscosity and thermal diffusivity, and that the Boussinesq approximation may be adopted for steady laminar flow, 2) the magnetic Reynolds number is taken to be small enough so that induced magnetic field is assumed to be negligible in comparison with applied magnetic field so that B = ( 0 , B 0 , 0 ) , where b_{0} is the uniform magnetic field acting normal to the plate, 3) Joule heating is negligible. Within framework of such assumptions, the governing equations for this problem can be written as follows:

Continuity Equation:

∂ u ∂ x + ∂ v ∂ y = 0 (1)

Momentum Equation:

u ∂ u ∂ x + v ∂ u ∂ y = υ a ∂ 2 u ∂ y 2 + S ρ ∂ σ ∂ y + g 0 β ( T − T ∞ ) cos α − σ 0 B 0 2 ( u − U ∞ ) ρ (2)

Angular Momentum Equation:

u ∂ σ ∂ x + v ∂ σ ∂ y = υ s ρ j ∂ 2 σ ∂ y 2 − S ρ j ( 2 σ + ∂ u ∂ y ) (3)

Energy Equation:

u ∂ T ∂ x + v ∂ T ∂ y = κ ρ c p ∂ 2 T ∂ y 2 + υ a c p ( ∂ u ∂ y ) 2 (4)

where u, v are the velocity components along x, y co-ordinates respectively, v is the kinematic viscosity υ a = μ + S ρ is the apparent kinematic viscosity, μ is the

coefficient of dynamic viscosity, S is the microrotation coupling coefficient (also known as the coefficient of gyro-viscosity or as the vortex viscosity), ρ is the mass density of the fluid, σ is the microrotation component normal to the

xy-plane, b_{0} is the uniform magnetic field acting normal to the plate, υ s = ( μ + S 2 ) j

is the microrotation viscosity or spin-gradient viscosity, σ 0 is the magnetic permeability, j is the micro-inertia density, T is the temperature of the fluid within the boundary layer, T ∞ is the temperature of the fluid outside the boundary layer, U ∞ is the velocity of the fluid far away from the plate, c p is the specific heat of the fluid at constant pressure, k is the thermal conductivity, g_{0} is the acceleration due to gravity, β is the volumetric coefficient of thermal expansion.

The appropriate boundary conditions of the model are:

u = 0 , v = v 0 ( x ) , σ = − n ∂ u ∂ y , T = T w at y = 0 u = U ∞ , σ = 0 , T = T ∞ as y → ∞ } (5)

where the subscripts w and ∞ refer to the wall and boundary layer edge, respectively. The value v 0 = 0 corresponds to an impermeable plate. When microrotation parameter n = 0 we obtain σ = 0 which represents no-spin condition i.e. the microelements in a concentrated particle flow-close to the wall are not able to rotate as stated by Jena and Mathur [

In order to obtain similarity solution of the problem we introduce the following non-dimensional variables:

η = y U ∞ 2 υ a x , ψ = 2 υ a U ∞ x f ( η ) , σ = U ∞ 3 2 υ a x g , θ = T − T ∞ T w − T ∞ (6)

where ψ is the stream function, η is the dimensionless distance normal to the sheet, f is the dimensionless stream function, θ is the dimensionless fluid temperature

Since u = ∂ ψ ∂ y and v = − ∂ ψ ∂ x we have the velocity components

u = U ∞ f ′ , v = − υ a U ∞ 2 x ( f − η f ′ ) (7)

where prime denotes the derivative with respect to η .

Now substituting Equations (6) and (7) into Equations (2)-(4) we obtain

f ‴ + f f ″ + K g ′ + 2 γ θ cos α − 2 M ( f ′ − 1 ) = 0 (8)

G 2 g ″ − 2 G 1 ( 2 g + f ″ ) + f ′ g + g ′ f = 0 (9)

θ ″ + P r f θ ′ + P r E c ( f ″ ) 2 = 0 (10)

In the above equations we have used the following non-dimensional parameters:

K = S ρ υ a is the coupling parameter,

γ = G r x Re x 2 is the Richardson parameter,

G r x = g 0 β ( T w − T ∞ ) x 3 υ a 2 is the local Grashof number,

Re x = U ∞ x υ a is the local Reynolds number,

G 1 = S x ρ j U ∞ is the vortex viscosity parameter,

G 2 = υ s ρ j υ a is the spin gradient viscosity parameter,

M = σ 0 B 0 2 x ρ U ∞ is the local magnetic field parameter,

Pr = ρ υ a c p k is the Prandtl number.

The corresponding boundary conditions become

f = f w , f ′ = 1 , g = − n f ″ , θ = 1 at η = 0 f ' = 1 , g = 0 , θ = 0 as η → ∞ } (11)

where f w = ∓ v 0 2 x υ a U ∞ is the suction/injection velocity at the plate for v 0 > 0 and v 0 < 0 .

The parameters f W γ , M , G 1 and f w correspond to local effects i.e. pertaining to specific values of x . Therefore (8)-(10) are ordinary differential equations that are locally similar together with the boundary conditions (11). The above-noted systems have been solved numerically for various values of the parameters entering into the problem.

The set of nonlinear ordinary differential Equations (8)-(10) with boundary conditions (11) are non-linear and coupled. It is difficult to solve them analytically. Hence we adopt a procedure to obtain the solution numerically. Here we use the standard initial-value solver shooting method namely Nachtsheim-Swigert (1965) shooting iteration technique (guessing the missing value) together with sixth order Runge-Kutta-Butcher initial value solver.

Important Physical ParametersThe quantities of chief physical interest are the skin friction coefficient (viscous drag), plate couple stress and the Nusselt number i.e., rate of heat transfer.

Skin friction coefficient:

The equation defining the wall shear stress is

τ w = ( μ + S ) ( ∂ u ∂ y ) y = 0 + S ( σ ) y = 0 (12)

The local skin friction coefficient is defined as

C f = τ w 1 2 ρ U ∞ 2 = ( 2 Re x − 1 ) 1 / 2 [ f ″ ( 0 ) + K g ( 0 ) ] = ( 2 Re x − 1 ) 1 2 ( 1 − n K ) f ″ ( 0 ) (13)

Thus from Equation (13) we see that the local values of the skin friction coefficient C f is proportional to f ″ ( 0 ) .

Plate couple stress:

The equation defying the plate couple stress is

M w = υ s ( ∂ σ ∂ y ) y = 0 (14)

The dimensionless couple stress is defined by

M x = M w 1 2 ρ υ a U ∞ = G 2 K G 1 g ′ ( 0 ) (15)

Thus the local plate couple stress in the boundary layer is proportional to g ′ ( 0 ) .

Nusselt number:

The local heat flux may be written by Fourier’s law as

q w ( x ) = − k ( ∂ T ∂ y ) y = 0 = − k ( T w − T ∞ ) ( ρ U ∞ 2 μ x ) 1 2 θ ′ ( 0 ) (17)

The local heat transfer coefficient is given by

h ( x ) = q w ( x ) ( T w − T ∞ ) (17)

The local Nusselt number may be written as

N u x = x h ( x ) k = − ( 1 2 Re x ) 1 2 θ ′ ( 0 ) (18)

Thus from Equation (18) we see that the local Nusselt number N u x is proportional to – θ ′ ( 0 ) .

The velocity, temperature and microrotation profiles of convective flow of a viscous incompressible electrically conducting micropolar fluid past an inclined porous flat plate considering viscous dissipation effect has been investigated numerically by using Nachtsheim-Swigert shooting iteration technique. The values of the parameters γ , K , G 1 , G 2 are taken to be fixed as 1, 0.2, 0.2 and 2.0 respectively, whereas the effects of the other important parameters, namely magnetic field parameter M , Prandtl number Pr , Eckert number E c , suction parameter f w over these profiles have been studied. The value of the Richardson parameter (buoyancy parameter) γ = 1 corresponds to mixed convection.

Figures 2(a)-(c), respectively, show the velocity, temperature and microrotation profiles for different values of Eckert number. From

Figures 3(a)-(c), respectively, illustrate the velocity, temperature and microrotation profiles for different values of suction parameter f w . From

Figures 4(a)-(c), respectively, describe the velocity, temperature and microrotation profiles for different values of magnetic field parameter M. From

Figures 5(a)-(c), respectively, are drawn to show the velocity, temperature and microrotation profiles for different values of Prandtl number Pr. From

In Figures 6(a)-(c), respectively, we have presented the physical parameters f ″ ( 0 ) , g ′ ( 0 ) and – θ ′ ( 0 ) which are respectively proportional to skin friction coefficient, plate couple stress and Nusselt number for different values of angle of inclination α and magnetic field parameter M. From _{f}) and plate couple stress ( M x ) decreases with the increase of α and increases with the increase of M.

Figures 7(a)-(c), respectively, show that skin friction coefficient, plate couple stress and Nusselt number for different values of Eckert number Ec and suction parameter f w . From

A numerical study has been performed to investigate the effect of various parameters such as suction parameter f w , Prandtl number Pr, Eckert number Ec, magnetic field parameter M on the velocity, temperature and microrotation profiles and discussed them from the physical point of view.

Skin friction coefficient increases with the increase of magnetic field parameter M, suction parameter f w , Eckert number Ec. This coefficient decreases with the increase of angle of inclination α and Prandtl number Pr.

Nusselt number increases with the increase of suction parameter f w , angle of inclination α, and Prandtl number Pr. Nusselt number decreases with the increase of magnetic field parameter M, Eckert number Ec.

Plate couple stress increases with the increase of magnetic field parameter M, suction parameter f w , Eckert number Ec. This coefficient decreases with the increase of angle of inclination α and Prandtl number Pr.

Sheikh, N. and Hasan, M. (2017) Mixed Convective Flow of Micropolar Fluids past an Inclined Porous Flat Plate. Open Journal of Fluid Dynamics, 7, 642-656. https://doi.org/10.4236/ojfd.2017.74042