International Journal of Engineering MathematicsVolume 2013 (2013), Article ID 903818, 9 pageshttp://dx.doi.org/10.1155/2013/903818

Research Article

Effects of Mass Transfer, Radiation, Joule Heating, and Viscous Dissipation on Steady MHD Marangoni Convection Flow over a Flat Surface with Suction and Injection

S. Mohammed Ibrahim

Department of Mathematics, Priyadarshini College of Engineering & Technology, Nellore 524004, India

Received 4 August 2013; Accepted 6 November 2013

Academic Editor: Giuseppe Carbone

Copyright © 2013 S. Mohammed Ibrahim. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The combined effects of radiation and mass transfer on a steady MHD two-dimensional Marangoni convection flow over a flat surface in presence of Joule heating and viscous dissipation under influence of suction and injection is studied numerically. The general governing partial differential equations are transformed into a set of nonlinear ordinary differential equations by using unique similarity transformation. Numerical solutions of the similarity equations are obtained using the Runge-Kutta method along with shooting technique. The effects of governing parameters on velocity, temperature, and concentration as well as interface velocity, the surface temperature gradient, and the surface concentration gradient were presented in graphical and tabular forms. Comparisons with previously published work are performed and the results are found to be in excellent agreement.

1. Introduction

In recent years, many researchers have investigated the Marangoni convection and applied such convection to their problems. This is mainly because it is important in crystal growth melts and greatly influences other industrial processes. Nishino and Kawamura [1] stated that Marangoni convection negatively affects the quality of silicon crystals for semiconductors and the convection also occurs in heat pipe for heat radiation devices of personal computers. The surface tension gradient variations along the interface may induce the Marangoni convection. In particular, the surface tension gradients that are responsible for Marangoni convection can be due to gradients of temperature (thermal convection) and/or concentration (solutal convection). A lot of analyses in Marangoni convection have been discovered in various geometries and conditions. Some of experimental works linked to Marangoni convection were discussed in several papers by Arafune and Hirata [2], Arafune et al. [3], Galazka and Wilke [4], Neumann et al. [5], Arendt and Eggers [6], Xu et al. [7], and Christopher and Wang [8].

On the other hand, the study of Magnetohydrodynamics (MHD) is important in the heat and mass transport process. The study of heat transfer is integral part of natural convection flow and belongs to the class of problems in boundary layer theory. The quality of heat transferred is highly dependent upon the fluid motion within the boundary layer. A large number of physical phenomena involve natural convection that was studied by Jaluria [9], which are enhanced and driven by internal heat generation. In such flows the buoyancy force is incremented due to heat generation resulting in modification of heat transfer characteristic. The effect of internal heat generation is especially pronounced for low Prandtl number fluid. MHD is the study of motion of electrically conducting incompressible fluid in the presence of magnetic field, that is, an electromagnetic field interacting with the velocity field of an electrically conducting fluid. Hydromagnetic flows have become important due to industrial applications; for instance it is used to deal with the problem of cooling of nuclear reactor by fluid having very low Prandtl number that was studied by Itaru et al. [10] and Fumizawa [11]. Moreover, Al-Mudhaf and Chamkha [12] investigated the similarity solution for MHD thermosolutal Marangoni convection over a flat surface in the presence of heat generation or absorption with fluid suction and injection.

It was realized that the studies of thermal radiation with heat and mass transfer are important in electrical power generation, astrophysical flows, solar power technology, and other industrial areas. Pathak and Maheshwari [13] have analyzed the influence of radiation on an unsteady free convection flow bounded by an oscillating plate with variable wall temperature. The effects of thermal radiation, buoyancy, and suction/blowing on natural convection heat and mass transfer over a semi-infinite stretching surface have been studied by Shateyi [14]. Moreover, Suneetha et al. [15] have investigated the radiation effects on the MHD free convection flow past an impulsively started vertical plate with variable surface temperature and concentration.

The effects of viscous dissipation and Joule heating are usually characterized by the Eckert number and the product of the Eckert number and magnetic parameter, respectively, and both have a very important part in geophysical flow and in nuclear engineering that was studied by Alim et al [16], the effects of suction or injection on boundary layer flow also have a huge influence over the engineering application and have been widely investigated by numerous researchers. With this understanding, many researchers studied the effects of suction or injection in various geometries, for example, Borisevish and Potanin [17] for heat transfer near a rotating disk. Duwairi [18] and Chen [19] for MHD convection flow in the presence of radiation.

However, the study of the combined effects of mass transfer and radiation on steady MHD Marangoni convection flow of a dissipative fluid has received a little attention. Hence, the object of the present paper is to analyze the combined effects of radiation and mass transfer on steady MHD laminar Marangoni convection boundary layer flow of an electrically conducting fluid past a flat surface, by taking Joule heating and viscous dissipation under influence of suction or injection. The governing partial differential equations are reduced to a system of self-similar equations using the similarity transformations. The resultant equations are then solved numerically using the Runge-Kutta fourth order technique along with shooting method. The effects of governing physical parameters on the velocity, temperature, and concentration as well as surface velocity, surface temperature, and surface concentration gradient are computed and presented in graphical and tabular forms. To verify the obtained results, we have compared the present numerical results with previous work by Al-Mudhaf and Chamkha [12]. The comparison results show a good agreement and we are confident that our present numerical results are accurate.

2. Mathematical Analysis

A steady, two-dimensional, laminar boundary layer flow of a viscous incompressible electrically conducting and radiating fluid over a flat surface is considered. The fluid is assumed to be gray, absorbing but nonscattering. The -axis is taken along the surface and -axis normal to it. The surface is assumed to be in the presence of the surface tension due to temperature and concentration gradients at the wall. Further, a uniform magnetic field of strength is applied normal to the surface which then produces the magnetic forces along the surface. The fluid is assumed to be slightly conducting, and hence the magnetic Reynolds number is much less than unity and the induced magnetic field is negligible in comparison with the applied magnetic field. Under the above assumptions and with the usual boundary layer approximations, the governing equations are [20] as follows.

Continuity equation Momentum equation Energy equation Species equation The boundary conditions for the velocity, temperature, and concentration fields are where and are the velocity components in the and directions, respectively. is the kinematic viscosity, is the electric conductivity, is the surface tension, is the uniform magnetic field strength, is the density of the fluid, is the specific heat at constant pressure, is the fluid temperature, is the fluid temperature far from the surface, is the species concentration, is the species concentration far from the surface, is the temperature gradient coefficient, is the concentration gradient coefficient, is the dynamic viscosity, is the thermal conductivity of the fluid, is the radiative heat flux, and is constant suction or injection velocity. The second, third, and fourth terms on the right hand side of (3) represent the radiation, Joule heating, and viscous dissipation, respectively. The fourth condition of (5) represents the Marangonic coupling condition at the interface.

When we use the Rosseland approximation for radiation, the radiative heat flux can be simplified as [21] where is the Stefan-Boltzmann constant and is the mean absorption coefficient. It is assumed that the temperature differences within the flow are small, so that the term may be expressed as a linear function of temperature. Hence by expanding in a Taylor’s series about and neglecting higher-order terms, Using (7) and (8), (3) reduces [22]: where is the thermal diffusivity and is the radiation parameter.

Further, we use the similarity transformation by Al-Mudhaf and Chamkha [12] and the standard definition of the stream function such that and to obtain the similarity solution of the problems. The similarity transformations are given by where with being the length of the surface and and are the constant characteristics temperature and concentration, respectively. Substituting (10) into (2), (4), and (9), we obtain the following ordinary differential equations: where primes denote differentiation with respect to , is the Prandtl number, is the magnetic field parameter, is the Eckert number, is the Schmidt number, (>0) is the constant suction parameter, and (<0) is the constant injection parameter. It should be mentioned here again that the viscous dissipation effect is examined using the Eckert number , while the product of the Eckert number and the magnetic field parameter gives the joule heating.

The dimensionless form of the boundary conditions become

3. Solution of the Problem

The nonlinear ordinary differential equations (12) subjects to the boundary conditions (13) are solved numerically using the Runge-Kutta fourth order technique along with shooting method. In this method, it is most important to choose the appropriate finite value of the edge of boundary layer, (say ), that is between 4 to 8, which is in accordance with the standard practice in the boundary layer analysis. First of all, higher-order nonlinear differential equations (12) are converted into simultaneous linear differential equations and they are further transformed into initial value problem by applying the shooting method. The resultant initial value problem is solved by employing Runge-Kutta fourth order technique. The initial step size is used to obtain the numerical solution with five decimal place accuracy as the criterion of convergence.

4. Results and Discussion

In order to get a physical insight into the problem, a representative set of numerical results are shown in Figures 1–18 which illustrates the influence of physical parameters namely, the magnetic parameter , Prandtl number , radiation parameter Nr, Eckert number , Schmidt number , Suction, and Injection on the velocity , temperature and concentration , profiles as well as the reduced velocity at the interface, reduced heat transfer, and reduced mass transfer at the interface. Throughout the calculations, the parametric values are fixed to be , , , , , and , unless otherwise indicated.

Figure 1: Effects of on velocity profiles with various values of .

Figure 2: Effects of on temperature profiles with various values of .

Figure 3: Effects of on concentration profiles with various values of .

Figure 4: Effects of Pr on temperature profiles with various values of .

Figure 5: Variation of radiation parameter with interface heat transfer for different values of .

Figure 6: Temperature profiles for different values of when .

Figure 7: Temperature profiles for different values of when .

Figure 8: Temperature profiles for different values of when .

Figure 9: Variation of Eckert number Ec with interface heat transfer for different values of .

Figure 10: Temperature profiles for different values of Ec when .

Figure 11: Temperature profiles for different values of Ec when .

Figure 12: Temperature profiles for different values of Ec when .

Figure 13: Variation of heat Schmidt number Sc with interface mass transfer for different values of .

Figure 14: Concentration profiles for different values of .

Figure 15: Velocity profiles for different values of .

Figure 16: Temperature profiles for different values of .

Figure 17: Concentration profiles for different values of .

Figure 18: Variation of the surface temperature gradient with the Eckert number for different values of .

Figures 1, 2, and 3 show the effect of magnetic parameter in the influence of suction and injection on velocity, temperature, and concentration profiles. Application of a transverse magnetic field results in a drag-like force called the Lorentz force. This force tends to slow down the movement of the fluid along surface and to increase its temperature and the concentration species. This is evident in the decreases in the velocities and increases in the temperature and concentration as increases. In addition, as the strength of magnetic field increases, the hydrodynamic velocity boundary layer decreases while the thermal and solutal (concentration) boundary layer. It is also noticed from Figure 1 that the wall velocity is nonzero due to the Marangoni or surface tension effect and it decreases as increases. These behaviors are depicted in Figures 1, 2, and 3.

Figure 4 illustrates the effect of Prandtl number in the influence of suction and injection on temperature profile. The result indicates that the increasing of Prandtl number decreases the temperature profile.

Effect of radiation parameter on temperature of the fluid with the influence of the suction or injection parameter is presented in Figures 5, 6, 7, and 8. Figure 5 displays the variation with of the reduced temperature gradient, , with different values of . One can see that the reduced temperature gradient, , decreases as Nr increases. Thus, the heat transfer rate at the surface decreases in the presence of radiation. This result qualitatively agrees with expectations, since the effect of radiation is to decrease the rate of energy transport to the fluid, thereby decreasing the temperature of the fluid. From Figures 6–8, we observe that the temperature profiles increase as radiation parameter increases. Thus, radiation can be used to control the thermal boundary layers quite effectively.

In Figure 9 we can see the combined effects of Joule heating and viscous dissipation on heat transfer rate. When , there are no effects of Joule and viscous heating. It can be observed that reduced sharply with the increasing of . Further, imposition of the fluid suction has the tendency to increase the interface heat transfer. On other hand, fluid injection tends to decrease the heat transfer rate at the wall but rapidly increase the rate as . Meanwhile in Figures 10, 11, and 12 present the effect of the viscous dissipation on temperature profiles in the presence of suction and injection. It is noticed that increases as increase the temperature.

The influence of the Schmidt number on the concentration of the fluid in presence of suction or injection parameter is presented in Figures 13 and 14. Figure 13 presents the effect of Schmidt number on the reduced species gradient, , under the presence of different values of . It is noted that fluid suction or injection has the tendency to increase the interface mass transfer. The effect of Schmidt number on the concentration is presented in Figure 14. It is evident that the concentration decreases with an increase in .

Figures 15, 16, and 17 illustrate the influence of the suction or injection parameter on the velocity, temperature, and concentration, respectively. Physically speaking, imposition of fluid suction at the wall has the tendency to decrease the fluid velocity and thickness of the hydrodynamic boundary layer. As a result, the fluid temperature and concentration boundary layer decrease as well. However, fluid injection produces the opposite effect, namely, an increase in the fluid velocity, temperature, and concentration.

Figure 18 illustrates the variation of the surface temperature gradient with the Eckert number Ec for different values of the Prandtl number, that is, = 0.72, 7.0 and 10. From Figure 18, it is observed that as Ec increases the surface temperature gradient is found to be decreased in the presence of thermal radiation Nr. Furthermore, if we consider higher Prandtl number then we notice further reduction in the surface temperature gradient.

Tables 1 and 2 show the comparison of and with those reported by Al-Mudhaf and Chamkha [12], which show an excellent agreement and we are confident that our present numerical results are accurate.

Table 1: Comparison values of and with different parameter M.

Al-Mudhaf and Chamkha [12]
Present

0
1.587671
1.442203
1.58155
1.440125

1
1.315181
1.206468
1.31379
1.20448

2
0.903945
0.759604
0.90321
0.75752

3
0.644888
0.442240
0.64402
0.44196

4
0.493358
0.272847
0.49248
0.27452

Table 2: Comparison values of and with different parameter .

Al-Mudhaf and Chamkha [12]
Present

−2
2.383451
1.251341
2.38963
1.25986

−1
2.000379
1.336441
2.00753
1.32874

0
1.587671
1.442203
1.58125
1.43987

1
1.179708
1.634990
1.17324
1.63123

2
0.848026
2.020945
0.83924
2.01785

Further, Table 3 shows the effects of magnetic parameter , Prandtl number , radiation parameter , viscous dissipation , Schmidt number and suction on the physical parameters surface velocity gradient , surface temperature gradient and surface concentration gradient , respectively. It can be seen that all of , , and decrease as magnetic field parameter increases. increases as Prandtl number increases. decreases as radiation and viscous dissipation increase. decreases as Schmidt number increase and as suction parameter increases decreases while and increase.

Table 3: Computation showing and for different embedded flow parameter values.

Pr
Nr
Ec
Sc

1.0
0.72
0.5
0.1
0.6
1.0
0.999612
0.941925
1.22659

2.0
0.72
0.5
0.1
0.6
1.0
0.732048
0.777927
1.03547

3.0
0.72
0.5
0.1
0.6
1.0
0.550868
0.683809
0.915823

1.0
1
0.5
0.1
0.6
1.0
0.999612
1.21094
1.22659

1.0
2.0
0.5
0.1
0.6
1.0
0.999612
2.11956
1.22659

1.0
0.72
1.0
0.1
0.6
1.0
0.999612
0.78197
1.22659

1.0
0.72
2.0
0.1
0.6
1.0
0.999612
0.625854
1.22659

1.0
0.72
0.5
0.2
0.6
1.0
0.999612
0.887262
1.22659

1.0
0.72
0.5
0.3
0.6
1.0
0.999612
0.832598
1.22659

1.0
0.72
0.5
0.1
1.0
1.0
0.999612
0.941925
1.81808

1.0
0.72
0.5
0.1
1.5
1.0
0.999612
0.941925
2.51464

1.0
0.72
0.5
0.1
0.6
2.0
0.752099
1.19244
1.54709

1.0
0.72
0.5
0.1
0.6
3.0
0.578583
1.55153
1.99824

5. Conclusions

The combined effects of mass transfer and thermal radiation on a steady laminar MHD Marangoni convection boundary layer flow over a flat surface by taking the Joule heating and viscous dissipation in the influence of fluid suction and injection. The governing partial differential equations are reduced to a system of self-similar equations using the similarity transformations. The resultant equations are then solved numerically using the Runge-Kutta method along with shooting technique. Comparison with previously published work was performed and the results were found to be in excellent agreement. The effects of governing physical parameters on the velocity, temperature, and concentration as well as surface temperature gradient and surface concentration are computed and presented in graphical and tabular forms. It was found that the temperature increases as Eckert number or radiation parameter increases, while temperature decreases as Prandtl number or magnetic field increases, whereas concentration decreases as Schmidt number increases. It can be drawn from the present results that when the radiation parameter increases, the heat transfer rate at the surface decreases. Meanwhile, the imposition of suction is to decrease the fluid velocity, temperature, and concentration profiles, whereas injection shows the opposite effects. It should be noted that results obtained in this work can used for the analysis of Marangoni flow and heat and mass transfer for flow over curved surfaces provided that the curvature is much greater than the boundary layer thickness.

Conflict of Interests

Dr. S. Mohammed Ibrahim, the author of this paper declares that there is no conflict of interests regarding the publication of this paper.

References

K. Nishino and H. Kawamura, Chaos, Turbulence and Its Transition Process in Marangoni Convection, KIBO Japanese Experiment Module, Japan Aerospace Exporation Agency, 2007.

K. Arafune and A. Hirata, “Thermal and solutal marangoni convection in In-Ga-Sb system,” Journal of Crystal Growth, vol. 197, no. 4, pp. 811–817, 1999.

K. Arafune, K. Yamamoto, and A. Hirata, “Interactive thermal and solutal Marangoni convection during compound semiconductor growth in a rectangular open boat,” International Journal of Heat and Mass Transfer, vol. 44, no. 13, pp. 2404–2411, 2001.

Z. Galazka and H. Wilke, “Influence of Marangoni convection on the flow pattern in the melt during growth of Y₃Al₅O₁₂ single crystals by the Czochralski method,” Journal of Crystal Growth, vol. 216, no. 1, pp. 389–398, 2000.

H. Neumann, Y. Plevachuk, and F. Allenstein, “Investigation of Marangoni convection in monotectic melts by resistance measurements,” Materials Science and Engineering A, vol. 361, no. 1-2, pp. 155–164, 2003.

B. Arendt and R. Eggers, “Interaction of Marangoni convection with mass transfer effects at droplets,” International Journal of Heat and Mass Transfer, vol. 50, no. 13-14, pp. 2805–2815, 2007.

Y. L. Xu, Z. B. Dong, Y. H. Wei, and C. L. Yang, “Marangoni convection and weld shape variation in A-TIG welding process,” Theoretical and Applied Fracture Mechanics, vol. 48, no. 2, pp. 178–186, 2007.

D. M. Christopher and B. Wang, “Prandtl number effects for Marangoni convection over a flat surface,” International Journal of Thermal Sciences, vol. 40, no. 6, pp. 564–570, 2001.

Y. Jaluria, Natural Convection, Heat and Mass Transfer, Pergamon Press, New York, NY, USA, 1980.

M. Itaru, T. Osamu, and S. Akimi, “Natural convection heat transfer from a horizontal cylinder to mercury under magnetic field,” International Journal of Heat and Mass Transfer, vol. 19, no. 9, pp. 1021–1029, 1976.

M.-O. Fumizawa, “Natural convection experiment with liquid nak under transverse magnetic field,” Journal of Nuclear Science and Technology, vol. 17, no. 2, pp. 98–105, 1980.

A. Al-Mudhaf and A. J. Chamkha, “Similarity solutions for MHD thermosolutal Marangoni convection over a flat surface in the presence of heat generation or absorption effects,” Heat and Mass Transfer/Waerme- und Stoffuebertragung, vol. 42, no. 2, pp. 112–121, 2005.

G. Pathak and C. H. Maheshwari, “Effect of radiation on unsteady free convection flow bounded by an oscillating plate with variable wall temperature,” International Journal of Applied Mechanics and Engineering, vol. 11, no. 2, pp. 371–382, 2006.

S. Shateyi, “Thermal radiation and buoyancy effects on heat and mass transfer over a semi-infinite stretching surface with suction and blowing,” Journal of Applied Mathematics, vol. 2008, Article ID 414830, 12 pages, 2008.

S. Suneetha, N. B. Reddy, and V. R. Prasad, “Thermal radiation effects on MHD free convection flow past an impulsively started vertical plate with variable surface temperature and concentration,” Journal of Naval Architecture and Marine Engineering, vol. 2, pp. 57–70, 2008.

M. A. Alim, M. D. Alam, and A. Mamun, “Joule heating effect on the coupling of conduction with magnetohydrodynamic free convection flow from a vertical plate,” Nonlinear Analysis Modelling and Control, vol. 12, no. 3, pp. 307–316, 2007.

V. D. Borisevich and E. P. Potanin, “Effects of viscous dissipating and Joule heat on heat transfer near a rotating disk in the presence of intensive suction,” Journal of Engineering Physics, vol. 55, no. 5, pp. 1220–1223, 1988.

H. M. Duwairi, “Viscous and Joule heating effects on forced convection flow from radiate isothermal porous surfaces,” International Journal of Numerical Methods for Heat and Fluid Flow, vol. 15, no. 5, pp. 429–440, 2005.

C.-H. Chen, “Combined effects of Joule heating and viscous dissipation on Magnetohydrodynamic flow past a permeable, stretching surface with free convection and radiative heat transfer,” Journal of Heat Transfer, vol. 132, no. 6, pp. 1–5, 2010.

T. S. Khaleque and M. A. Samad, “Effects of radiation, heat generation and viscous dissipation on MHD free convection flow along a stretching sheet,” Research Journal of Applied Sciences, Engineering and Technology, vol. 2, no. 4, pp. 368–377, 2010.

M. Q. Brewster, Thermal Radiative Transfer Properties, John Wiley & Sons, New York, NY, USA, 1992.

A. Ishak, “Radiation effects on the flow and heat transfer over a moving plate in a parallel stream,” Chinese Physics Letters, vol. 26, no. 3, pp. 1–4, 2009.