# Effects of Mass Transfer, Radiation, Joule Heating, and Viscous

Authors: S. Mohammed Ibrahim
Journal: International Journal of Engineering Mathematics
 (2013)

---

## 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.


---

## Body

## 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. Thex-axis is taken along the surface and y-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 B0 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(1)∂u∂x+∂u∂y=0.
Momentum equation
(2)u∂u∂x+v∂u∂y=ν∂2u∂y2-ΔB02ρu.
Energy equation
(3)u∂T∂x+v∂T∂y=kρcp∂2T∂y2-1ρcp∂qr∂y+ΔB02ρcpu2+νcp(∂u∂y)2.
Species equation
(4)u∂C∂x+v∂C∂y=D∂2C∂y2.
The boundary conditions for the velocity, temperature, and concentration fields are
(5)v(x,0)=-vw,T(x,0)=T∞+Ax2,C(x,0)=C∞+A*x2,μ∂u∂y=-dσdT∂T∂x-dσdC∂C∂xaty=0(6)u(x,∞)=0,T(x,∞)=T∞,C(x,∞)=C∞asy⟶∞,
where u and v are the velocity components in the x andy directions, respectively. ν is the kinematic viscosity, Δ is the electric conductivity, σ is the surface tension, B0 is the uniform magnetic field strength, ρ is the density of the fluid, cp is the specific heat at constant pressure, T is the fluid temperature, T∞ is the fluid temperature far from the surface, C is the species concentration, C∞ is the species concentration far from the surface, A is the temperature gradient coefficient, A* is the concentration gradient coefficient, μ is the dynamic viscosity, k is the thermal conductivity of the fluid, qr is the radiative heat flux, and vw is constant suction (vw>0) or injection (vw<0) 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]
(7)qr=-4σ*3k*∂T4∂y,
where σ* is the Stefan-Boltzmann constant and k* is the mean absorption coefficient. It is assumed that the temperature differences within the flow are small, so that the term T4 may be expressed as a linear function of temperature. Hence by expanding T4 in a Taylor’s series about T∞ and neglecting higher-order terms,
(8)T4≈4T∞3T-3T∞4.
Using (7) and (8), (3) reduces [22]:
(9)u∂T∂x+v∂T∂y=α(1+Nr)∂2T∂y2+ΔB02ρcpu2+νcp(∂u∂y)2,
where α=k/ρcp is the thermal diffusivity and Nr=16σ*T∞3/3kk* 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 u=∂ψ/∂y and v=-∂ψ/∂x to obtain the similarity solution of the problems. The similarity transformations are given by
(10)η=C1y,f(η)=C2x-1ψ(x,y),θ(η)=[T(x,y)-T∞]Ax-2,ϕ(η)=[C(x,y)-C∞]A*x-2,
where
(11)A=ΔTL2,A*=ΔCL2,C1=ρA(dσ/dT)μ23,C2=ρ2μA(dσ/dT)3
with L being the length of the surface and ΔT and ΔC are the constant characteristics temperature and concentration, respectively. Substituting (10) into (2), (4), and (9), we obtain the following ordinary differential equations:
(12)f′′′+ff′′-(f′)2-M2f′=0,(1+Nr)Prθ′′-2f′θ+fθ′+Ec[M2(f′)2+(f′′)2]=0ϕ′′+Sc(fϕ′-2f′ϕ)=0,
where primes denote differentiation with respect to η, Pr=ν/α is the Prandtl number, M2=ΔB02C2/ρC1 is the magnetic field parameter, Ec=C12/AcpC22 is the Eckert number, Sc=ν/D is the Schmidt number, fw(>0) is the constant suction parameter, and fw(<0) is the constant injection parameter. It should be mentioned here again that the viscous dissipation effect is examined using the Eckert number Ec, while the product of the Eckert number and the magnetic field parameter M gives the joule heating.The dimensionless form of the boundary conditions become(13)f(0)=fw,f′′(0)=-2,θ(0)=1,ϕ(0)=1,f′(∞)=0,θ(∞)=0,ϕ(∞)=0.


## 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 h=Δη=0.01 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 Figures1–18 which illustrates the influence of physical parameters namely, the magnetic parameter M, Prandtl number Pr, radiation parameter Nr, Eckert number Ec, Schmidt number Sc,Suction, and Injection fw on the velocity f′(η), 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 Pr=0.72, M=1.0, Nr=0.5, Ec=0.1, Sc=0.6, and fw=1.0, unless otherwise indicated.Figure 1
Effects ofM on velocity profiles with various values of fw.Figure 2
Effects ofM on temperature profiles with various values of fw.Figure 3
Effects ofM on concentration profiles with various values of fw.Figure 4
Effects of Pr on temperature profiles with various values offw.Figure 5
Variation of radiation parameterNr with interface heat transfer -θ′(0) for different values of fw.Figure 6
Temperature profiles for different values ofNr when fw=1.Figure 7
Temperature profiles for different values ofNr when fw=0.Figure 8
Temperature profiles for different values ofNr when fw=-1.Figure 9
Variation of Eckert number Ec with interface heat transfer-θ′(0) for different values of fw.Figure 10
Temperature profiles for different values of Ec whenfw=1.Figure 11
Temperature profiles for different values of Ec whenfw=0.Figure 12
Temperature profiles for different values of Ec whenfw=-1.Figure 13
Variation of heat Schmidt number Sc with interface mass transfer-ϕ′(0) for different values of fw.Figure 14
Concentration profiles for different values ofSc.Figure 15
Velocity profiles for different values offw.Figure 16
Temperature profiles for different values offw.Figure 17
Concentration profiles for different values offw.Figure 18
Variation of the surface temperature gradient-θ′(0) with the Eckert number Ec for different values of Pr.Figures1, 2, and 3 show the effect of magnetic parameter M 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 M 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 M increases. These behaviors are depicted in Figures 1, 2, and 3.Figure4 illustrates the effect of Prandtl number Pr 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 parameterNr on temperature of the fluid with the influence of the suction or injection parameter fw is presented in Figures 5, 6, 7, and 8. Figure 5 displays the variation with Nr of the reduced temperature gradient, -θ′(0), with different values of fw. One can see that the reduced temperature gradient, -θ′(0), 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 Nr increases. Thus, radiation can be used to control the thermal boundary layers quite effectively.In Figure9 we can see the combined effects of Joule heating and viscous dissipation on heat transfer rate. When Ec=0, there are no effects of Joule and viscous heating. It can be observed that -θ′(0) reduced sharply with the increasing of Ec.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 Ec>3. 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 Ec increase the temperature.The influence of the Schmidt numberSc on the concentration of the fluid in presence of suction or injection parameter fw is presented in Figures 13 and 14. Figure 13 presents the effect of Schmidt number Sc on the reduced species gradient, -ϕ′(0), under the presence of different values of fw.It is noted that fluid suction or injection has the tendency to increase the interface mass transfer. The effect of Schmidt number Sc on the concentration is presented in Figure 14. It is evident that the concentration decreases with an increase in Sc.Figures15, 16, and 17 illustrate the influence of the suction or injection parameter fw on the velocity, temperature, and concentration, respectively. Physically speaking, imposition of fluid suction (fw>0) 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 (fw>0) produces the opposite effect, namely, an increase in the fluid velocity, temperature, and concentration.Figure18 illustrates the variation of the surface temperature gradient -θ′(0) with the Eckert number Ec for different values of the Prandtl number, that is, Pr = 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.Tables1 and 2 show the comparison of f′(0) and -θ′(0) 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 off′(0) and -θ′(0) with different parameter M.
M
Al-Mudhaf and Chamkha [12]
Present
f
′
(
0
)
-
θ
′
(
0
)
f
′
(
0
)
-
θ
′
(
0
)
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.27452Table 2
Comparison values off′(0) and -θ′(0) with different parameter fw.
f
w
Al-Mudhaf and Chamkha [12]
Present
f
′
(
0
)
-
θ
′
(
0
)
f
′
(
0
)
-
θ
′
(
0
)
−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.01785Further, Table3 shows the effects of magnetic parameter M, Prandtl number Pr, radiation parameter Nr, viscous dissipation Ec, Schmidt number Sc and suction fw on the physical parameters surface velocity gradient f′(0), surface temperature gradient -θ′(0) and surface concentration gradient -ϕ′(0), respectively. It can be seen that all of f′(0), -θ′(0), and -ϕ′(0) decrease as magnetic field parameter M increases. -θ′(0) increases as Prandtl number Pr increases. -θ′(0) decreases as radiation Nr and viscous dissipation Ec increase. -ϕ′(0) decreases as Schmidt number Sc increase and as suction parameter fwincreases f′(0) decreases while -θ′(0) and -ϕ′(0) increase.Table 3
Computation showingf′(0),-θ′(0) and -ϕ′(0) for different embedded flow parameter values.
M
Pr
Nr
Ec
Sc
f
w
f
′
(
0
)
-
θ
′
(
0
)
-
ϕ
′
(
0
)
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 numberEc or radiation parameter Nr increases, while temperature decreases as Prandtl number Pr or magnetic field Mincreases, whereas concentration decreases as Schmidt number Sc 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.



---

*Source: 903818-2013-12-14.xml*