MHD Free Convection Flow past an Inclined Stretching Sheet with Considering Viscous Dissipation and Radiation

Abstract

The present study concentrates on the analysis of MHD free convection flow past an inclined stretching sheet. The viscous dissipation and radiation effects are assumed in the heat equation. Approximation solutions have been derived for velocity, temperature, concentration, Nusselt number, skin friction and Sherwood number using Nachtsheim-Swigert shooting iteration technique along with the six-order Runge-Kutta iteration scheme. Graphs are plotted to find out the characteristics of different physical parameters. The variations of physical parameters on skin friction coefficient, Nusselt number and Sherwood number are displayed via table.

Share and Cite:

Hasan, M. , Karim, E. and Samad, A. (2017) MHD Free Convection Flow past an Inclined Stretching Sheet with Considering Viscous Dissipation and Radiation. Open Journal of Fluid Dynamics, 7, 152-168. doi: 10.4236/ojfd.2017.72010.

1. Introduction

In recent years, considerable interest has been shown in investigating radiation interaction with natural convection flow commonly known as free convection for heat transfer in fluid. This is due to the significant role of thermal radiation in the surface heat transfer when convection heat transfer is small particularly in free convection problems. Again the boundary layer flow on continuous surfaces is an important type of flow which occurs in a number of technical processes. Examples are paper production, crystal growing and glass blowing, aerodynamics extrusion of plastic sheets and fibers. Thus, the study of heat transfer has become important industrially for determining the quality of the final product. Laminar natural convection flow and heat transfer of fluid with and without heat source in channels with constant wall temperature have been extensively studied by Ostrach [1] . Hossain and Takhar [2] analyzed the effect of radiation effects on free convection flow of a gas past a semi infinite flat plate using the Cogley-Vincnitine-Giles equilibrium model. Ali et al. [3] studied the radiation effect on natural convection flow over a vertical surface in a gray gas. Following Ali et al., Mansour [4] studied the interaction of mixed convection with thermal radiation in laminar boundary layer flow over a horizontal, continuous moving sheet with suction/injection. Alabraba et al. [5] studied the same problem considering magnetic effect taking into account the chemical reaction and soret-doufour effects. Sattar and Kalim [6] made a study of the combined unsteady free convection dynamic boundary layer and thermal radiation boundary layer on a semi-infinite vertical plate by using the Rosseland diffusion approximation. Chen [7] investigated natural convection flow over a permeable surface with variable wall temperature and concentration.

The viscous dissipation effect plays an important role in natural convection. Natural convection flow is often encountered in the cooling of nuclear reactors. Viscous dissipation effects on non-linear MHD flow in a porous medium over a stretching porous surface have been studied by S.P. Anjali Devi and B. Ganga [8] . Jha and Ajibade [9] studied the effect of viscous dissipation on natural convection flow between vertical parallel plates with time-periodic boundary conditions. Ferdows et al. [10] described that in the presence of uniform magnetic field with viscous dissipation at the wall, the thermophoretic parameter is one of the most useful parameters to control the boundary layer of the fluid.

The analysis and discussion of natural convection flow, the viscous dissipation effect is generally ignored but here considered the combined effect of viscous dissipation and radiation on free convection flow an inclined stretching sheet.

2. Mathematical Analysis

A steady-state two-dimensional heat and mass transfer flow of an electrically conducting viscous incompressible fluid along an isothermal stretching permeable inclined sheet with an angle to the vertical embedded in a porous medium with heat generation/absorption is considered. A strong magnetic field is applied in the y-axis direction. Here the effect of the induced magnetic field is neglected in comparison to the applied magnetic field. The electrical current flowing in the fluid gives rise to an induced magnetic field if the fluid were an electrical insulator, but here we have taken the fluid to be electrically conducting. Hence, only the applied magnetic field of strength B0 plays a role which gives rise

to magnetic forces in x-direction, where is the electrical con-

ductivity assumed to be directly proportional to the x- translational velocity (u) of the fluid found by Helmy [11] and is the density of the fluid. Two equal and opposite forces are introduced along the x- axis so that the sheet is stretched keeping the origin fixed as shown in Figure 1.

The fluid is considered to be gray, absorbing emitting radiation but non-scat- tering medium and the Rosseland approximation is used to describe the radia-

Figure 1. Physical model and coordinate system.

tion heat flux in the energy equation. The radiative heat flux in the -direction is negligible to the flux in the y-direction. The plate temperature and concentration are initially raised to and respectively which are thereafter maintained constant. The ambient temperature of the flow is and the concentration of the uniform flow is.

Under the usual boundary layer and Boussinesq approximations and using the Darcy-Forchhemier model, the flow and heat transfer in the presence of radiation are governed by the following equations.

Continuity equation

(1)

Momentum equation

(2)

Energy equation

(3)

Concentration equation

(4)

where u and v are the velocity components in the x-direction and y-direction respectively, is the kinematic viscosity, is the acceleration due to gravity, is the volumetric coefficient of thermal expansion, is the angle of inclination, k is the Darcy permeability constant, and are the fluid temperature within the boundary layer and in the free-stream respectively, while C is the concentration of the fluid within the boundary layer, is the electric conductivity, is the uniform magnetic field strength (magnetic induction), is the density of the fluid, is the thermal conductivity of the fluid, is the specific heat at constant pressure, is the volumetric rate of heat generation/absorption and is the chemical molecular diffusivity.

The corresponding boundary conditions are

(5)

where is a constant, is a velocity component at the wall having positive value to indicate suction, is the uniform sheet temperature and is the concentration of the fluid at the sheet.

By using Rosseland approximation, takes the form

(6)

where is the Stefan-Boltzmann constant and is the mean absorption coefficient. It is assumed that the temperature difference within the flow are sufficiently small such that may be expressed as a linear function of temperature.

Thus (7)

Using the Equations (6) and (7) in Equation (3), we get

(8)

2.1. Similarity Analysis

In order to obtain similarity solution for the problem under consideration, we may take the following suitable similarity variables

(9)

where is the stream function, is the dimensionless distance normal to the sheet, is the dimensionless stream function, is the dimensionless fluid temperature and is the dimensionless concentration.

Since and we have the velocity components from Equa-

tion (6) given by

(10)

where prime denotes the derivative with respect to.

Now introducing the similarity variables from Equation (9) and using Equation (10), Equations (2), (8) and (4) are reduced to the dimensionless equations given by

(11)

(12)

(13)

where is the buoyancy parameter, is the magnetic field parameter, is the local Darcy number, is the Reynolds number, is the Forchhemier number, is the Prandtl number, is the heat source /sink parameter, is the Eckert number, is the radiation parameter, and is the Schmidt number.

The transformed boundary conditions are

(14)

where is the suction parameter for.

The nonlinear ordinary differential Equations (11), (12) and (13) under the boundary conditions (14) are solved numerically for various values of the parameters entering into the problems.

2.2. Skin Friction, Rate of Heat and Mass Transfer

The parameters of engineering interest for the present problem are the skin friction coefficient, local Nusselt number and Sherwood number which indicate physically the wall shear stress, the rate of heat transfer and the local surface mass flux respectively. From the following definitions

(15)

where is the viscosity, is the thermal conductivity and is the mass diffusivity. The dimensionless local wall shear stress, local surface heat flux and the local surface mass flux for an impulsively started plate are respectively obtained as

(16)

(17)

(18)

where is the Reynolds number. Hence the values proportional

to the skin-friction coefficient, Nusselt number and Sherwood number are, and respectively.

2.3. Numerical Computation

The numerical solutions of the non-linear differential Equations (11)-(13) under the boundary conditions (14) have been performed by applying a shooting method namely Nachtsheim and Swigert [12] iteration technique (guessing the missing values) along with sixth order Runge-Kutta iteration scheme. We have chosen a step size to satisfy the convergence criterion of in all cases. The value of has been found to each iteration loop by. The maximum value of to each group of parameters and has been determined when the values of the unknown boundary conditions at not change to successful loop with error less than.

In order to verify the effects of the step size , we have run the code for our model with three different step sizes as and in each case we have found excellent agreement among them shown in Figures 2-4.

3. Results and Discussion

For the purpose of discussing the results of the flow field represented in the Figure 1, the numerical calculations are presented in the form of non-dimensional velocity, temperature and concentration profiles. The value of buoyancy parameter is taken to be positive to represent cooling of the plate. The parame-

Figure 2. Velocity profiles for different step sizes.

Figure 3. Temperature profiles for different step sizes.

Figure 4. Concentration profiles for different step sizes.

ters are chosen arbitrarily where corresponds physically to air at 20˚C. corresponds to electrolyte solution such as salt water and corresponds to water, and and 1.0 correspond to hydrogen, water vapor and methanol respectively at approximate 25˚C and 1 atmosphere.

Due to free convection problem positive large values of is chosen. The value of Re is kept 100 and Da equal to 0.20. However, numerical computations have been carried out for different values of the suction parameter (Fw), magnetic field parameter (M), angle of inclination (), Prandtl number (Pr), heat source parameter (Q), radiation parameter (N) and Eckert number Ec. The numerical results for the velocity, temperature and concentration profiles are displayed in Figures 5-25.

The effect of the angle of inclination of the sheet on the velocity field is shown in the Figure 5. From this figure, we see that the velocity decreases with the increase of swiftly up to. After, the velocity increases because the buoyancy force decreases. Figure 6 shows that temperatures rise with the grow of. Finally, we observe that the angle of inclination affects the concentration very slowly near the plate surface. Away from the plate, however, the effect on the concentration profile is significant.

Figure 5. Velocity profiles for different values of α.

Figure 6. Temperature profiles for different values of α.

Figure 7. Concentration profiles for different values of α.

Figure 8. Velocity profiles for different values of Ec.

Figures 8-10 are drawn to discuss the influence of Eckert number Ec on velocity, temperature and concentration profiles. Figure 8 shows that the velocity profiles increase with the increase of Ec upto. After the velocity profiles reduce. Again Figure 9 shows quick increasing effect on temperature profiles. On the other hand, has significant decreasing effect on concentration profiles observed in Figure 10.

Figure 11 to Figure 13 demonstrate the impact of the suction parameter Fw on the velocity, temperature and concentration profiles. It is observed that when Fw increases, the velocity, temperature and concentration decrease monotoni-

Figure 9. Temperature profiles for different values of Ec.

Figure 10. Concentration profiles for different values of Ec.

Figure 11. Velocity profiles for different values of Fw.

Figure 12. Temperature profiles for different values of Fw.

Figure 13. Concentration profiles for different values of Fw.

Figure 14. Velocity profiles for different values of M.

Figure 15. Temperature profiles for different values of M.

Figure 16. Concentration profiles for different values of M.

Figure 17. Velocity profiles for different values of N.

Figure 18. Temperature profiles for different values of N.

Figure 19. Concentration profiles for different values of N.

Figure 20. Velocity profiles for different values of Pr.

Figure 21. Temperature profiles for different values of Pr.

Figure 22. Concentration profiles for different values of Pr.

Figure 23. Velocity profiles for different values of Q.

Figure 24. Temperature profiles for different values of Q.

Figure 25. Concentration profiles for different values of Q.

cally. These figures indicate that since the boundary layer thickness getting smaller but the temperature and concentration gradient at the stretching sheet getting steeper, the increasing suction enhance the heat and mass transfer coefficient.

The influence of magnetic field parameter M on velocity, temperature and concentration profiles are plotted in Figures 14-16. Here from Figure 14, we see that first the velocity decreases gradually up to as M increases. On the other hand, the temperature as well as concentration profiles increase with the increase of M.

Figures 17-19 describe the dimensionless velocity, temperature and concentration profiles for different values of radiation parameter N. A strong decline in the velocity and temperature profiles are caused by increasing N have found in Figure 17 and Figure 18. Again the concentration profiles are significantly increased for the increasing value of N shown in Figure 19.

The behavior of Prandtl number Pr on velocity, temperature and concentration distributions are displayed in Figures 20-22. From the Figure 20 and Figure 21, we see that the velocity profiles and the temperature profiles decrease rapidly with the increase of Prandtl number Pr. Again Figure 22 shows that the concentration profiles increase sharply as Prandtl number increases.

For different values of heat source parameter Q, the velocity, temperature and concentration profiles are illustrated in Figures 23-25. Here we have plotted non-dimensional velocity, temperature and concentration profiles against for some representative values of the heat source parameter. The positive value of Q represents source i.e., heat generation in the fluid. For heat generation, the peak velocity occurs near the surface of the stretching plate. This is corroborated by Figure 24 where it is seen that the temperatures do indeed rapid increase as Q increases. Figure 25 shows that the concentration profiles decrease with the increase of heat source parameter.

4. Conclusions

The main goal of this study was the mathematical and numerical study of the viscous dissipation and radiation effect on MHD free convection flow past an inclined stretching sheet. The numerical solutions of the governing differential equations were obtained by using the shooting method. We observed the behavior of the physical parameters and also commented the numerical results from their plots.

Finally, the effects of various parameters on the skin friction, local Nusselt number and local Sherwood number for are shown in the Table 1.

Table 1. and for different values

We can make the following conclusions from the present study:

1. The effect of heat generation is remarkable. An increase in heat generation results in increasing velocity and temperature within the boundary layer.

2. Eckert number has effects on velocity as well as temperature profiles.

Radiation has significant effects on temperature profiles. So we can control the temperature field by using this parameter.

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1] Ostrach, S. (1952) Laminar Natural Convection Flow and Heat Transfer of Fluid with and without Heat Source in Channels with Constant Wall Temperature. NACA TN 2863.
[2] Hossain, M.A. and Takhar, H.S. (1996) Radiation Effects on Mixed Convection Along a Vertical Plate with Uniform Surface Temperature. Heat and Mass Transfer, 31, 243-248.
https://doi.org/10.1007/BF02328616
[3] Ali, M.M., Chen, T.S. and B.F. Armaly (1984) Natural Convection Radiation Interaction in Boundary Layer Flow over Horizontal Surface. AIAA Journal, 22, 1797-1803.
https://doi.org/10.2514/3.8854
[4] Mamsour, M.A. (1990) Radiation and Free Convection Effects on the Oscillatory Flow Past a Vertical Plate. Astrophysics and Space Science, 166, 269-275.
https://doi.org/10.1007/BF01094898
[5] Alabrara, M.A., Bestman, A.R. and Ogulu, A. (1992) Laminar Convection in Binary Mixed of Hydromagnetic Flow with Radiation Heat Transfer. Astrophysics and Space Science, 195, 431-439. https://doi.org/10.1007/BF00646774
[6] Sattar, M.A. and Kalim, H. (1996) Unsteady Free Convection Interaction with Thermal Radiation in a Boundary Layer Flow past a Vertical Porous Plate. Journal of Mathematical and Physical Sciences, 30, 25-37.
[7] Chen, C.-H. (2004) Heat and Mass Transfer in MHD Flow by Natural Convection from a Permeable, Inclined Surface with Variable Wall Temperature and Concentration. Acta Mechanica, 172, 219-235. https://doi.org/10.1007/s00707-004-0155-5
[8] Anajali Devi, S.P. and Ganga, B. (2009) Viscous Dissipation Effects on Non-Linear MHD Flow in a Porous Medium over a Stretching Surface. International Journal of Applied Mathematics and Mechanics, 5, 45-59.
[9] Jha, B.K. and Ajibade, A.O. (2011) Effect on Viscous Dissipation on Natural Convection Flow between Vertical Parallel Plates with Time-Periodic Boundary Conditions. Thermophysics and Aeromechanics, 18, 561-571. https://doi.org/10.1134/S0869864311040056
[10] Ferdows, M., Nazmul and Ota, M. (2011) Thermophoresis and Chemical Reaction Effects on MHD Natural Convective Heat and Mass Transfer Flow in a Rotating Fluid Considering Heat and Mass Fluxes. Canadian Journal on Science and Engineering Mathematics, 2, 114-139.
[11] Hemly, K.A. (1995) MHD Boundary Layer Equations for Power Law Fluids with Variable Electric Conductivity. Meccanica, 30, 187-200. https://doi.org/10.1007/BF00990456
[12] Nachtsheim, P.R. and Swigert, P. (1965) Satisfaction of the Asymptotic Boundary Conditions in Numerical Solution of the System of Non-Linear Equations of Boundary Layer Type. NASA TND-3004.

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.