Open access peer-reviewed chapter

Thermal Radiation and Thermal Diffusion for Soret and Dufour’s Effects on MHD Flow over Rotating Infinite Disk

Written By

Gamal M. Abdel-Rahman and Faiza M.N. El-fayez

Submitted: 07 September 2018 Reviewed: 22 October 2018 Published: 22 January 2019

DOI: 10.5772/intechopen.82186

From the Edited Volume

Nanofluid Flow in Porous Media

Edited by Mohsen Sheikholeslami Kandelousi, Sadia Ameen, M. Shaheer Akhtar and Hyung-Shik Shin

Chapter metrics overview

1,081 Chapter Downloads

View Full Metrics

Abstract

In general, the thermal radiation and thermal diffusion effects over an electrically conducting, Newtonian fluid in a steady laminar magnetohydrodynamic convective flow over a porous rotating infinite disk with the consideration of heat and mass transfer in the presence of Soret and Dufour’s diffusion effects have been obtained and studied numerically. The governing continuity, momentum, energy and concentration equations are converted into a system of non-linear ordinary differential equations by means of similarity transformation. The resulting system of coupled non-linear ordinary differential equations is solved numerically. In this chapter, numerical results were presented for velocity (radial, axial and tangential), temperature, concentration and pressure profiles for different parameters of the problem Also, the effects of the pertinent parameters on the radial and tangential skin friction, the rate of heat and mass transfer are obtained and discussed numerically and illustrated graphically.

Keywords

  • MHD
  • heat and mass transfer
  • thermal radiation
  • slip flow
  • porous medium
  • rotating disk

1. Introduction

The flow due to rotating disks is of great interest in many practical and engineering aspects. Rotating disk flows of electrically conducting fluids have practical applications in many areas, such as rotating machinery, lubrication, oceanography, computer storage devices, viscometer and crystal growth processes etc. Also, the study is interesting from the mathematical point of view. During the last two decades, research on renewable energy sources, as for example, solar, wind energy or energy from hydro-power and the preparation of oxygenated additives to blend diesel fuel, has been intensified.

Pioneering study of fluid flow due to an infinite rotating disk was carried by authors [1, 2, 3]. Chemical reactions usually accompany a large amount of exothermic and endothermic reactions. These characteristics can be easily seen in a lot of industrial processes, it has been realized that it is not always permissible to neglect the convection effects in porous constructed chemical reactors [4]. The reaction produced in a porous medium was extraordinarily in common, such as the topic of PEM fuel cells modules and the polluted underground water because of discharging the toxic substance, etc.

Fourier’s law, for instance, described the relation between energy flux and temperature gradient. In other aspects, Fick’s law was determined by the correlation of mass flux and concentration gradient. Moreover, it was found that energy flux can also be generated by composition gradients, pressure gradients, or body forces. The energy flux caused by a composition gradient was discovered in 1873 by Dufour and was correspondingly referred to the Dufour effect.

It was also called the diffusion-thermo effect. On the other hand, mass flux can also be created by a temperature gradient, as was established by Soret. This is the thermal-diffusion effect. In general, the thermal-diffusion and the diffusion-thermo effects were of a smaller order of magnitude than the effects described by Fourier’s or Fick’s law and were often neglected in heat and mass transfer processes. There were still some exceptional conditions. The thermal-diffusion effect has been utilized for isotope separation and in mixtures between gases with very light molecular weight H2He and of medium molecular weight N2air, the diffusion-thermo effect was found to be of a magnitude such that it may not be neglected in certain conditions [5].

The first traceable interest in magnetohydrodynamics (MHD) flow was in 1907, when Northrop built an MHD pump prototype [6, 7]. Since then, analysis of the effects of both rotation and magnetic fields on fluid flows has been an active area of research. While technology expanded in many directions, the subject of MHD has developed in the use of magnetic fields and the range of fluid and thermal processes by [8, 9, 10, 11, 12, 13]. This study considers the effect of slip as a result of rarefied effect, a type of flow commonly encountered in many engineering tasks such as high altitude flight, micro-machines, vacuum technology, aerosol reactors, etc. In this study, the slip and no-slip regimes that lie in the range 0.1 > Kn > 0 are considered. A completely different extension of von Karman’s one-disk problem is the analysis of Sparrow et al. [14]. They considered the flow of a Newtonian fluid due to the rotation of a porous-surfaced disk and for that purpose replaced the conventional no-slip boundary conditions at the disk surface with a set of linear slip flow conditions. A substantial reduction in torque then occurred as a result of surface slip. Recently Frusteri and Osalusi [15] studied the effects of variable properties on MHD and slip flow over a porous rotating disk.

In all these studies Soret and Dufour effects were assumed to be negligible. Such effects are significant when density differences exist in the flow regime. For example, when species are introduced at a surface in fluid domain, with different (lower) density than the surrounding fluid, both Soret (thermo-diffusion) and Dufour (diffusion-thermal) effects can be influential. An analytical study of convection along a horizontal cylinder for a Helium-air system was reported subsequently by Sparrow et al. [16]. In view of the importance of above mentioned effects, Maleque [17] studied Soret effect on convective heat and mass transfer past a rotating porous disk and he neglected the Dufour effect. Ahmed [18] investigated the Dufour and Soret effects on free convective heat and mass transfer over a stretching surface considering suction or injection. Recently, numerical study of free convection magnetohydrodynamic heat and mass transfer due to a stretching surface under saturated porous medium with Soret and Dufour effects was also discussed by Anwar Beg et al. [19].

In these two papers [15, 20] they have studied the effect of the magnetic field on the equations of motion that I found them used in the cylindrical coordinate, although the magnetic field parameters are in the Cartesian coordinate, which is wrong.

Advertisement

2. Mathematical formulation

Consider the steady, axially symmetric, incompressible flow of an electrically conducting fluid with heat and mass transfer flow due to a rotating porous disk in the presence of radiation has been considered. Assume that the fluid is infinite in extent in the positive z-direction. Let rφz be the set of cylindrical polar coordinates and let the disk rotate with constant angular velocity Ω and is placed at z=0. The components of the flow velocity are uvw in the directions of increasing rφz, respectively. p, TandC are the pressure, the temperature and the concentration distribution, respectively. The surface of the rotating disk is maintained at a uniform temperature Tw and uniform concentration Cw. Far away from the surface, the free stream is kept at a constant temperature T, at a constant concentration C and at a constant pressure p. The fluid is assumed to be gray, emitting and absorbing heat, but not scattering medium and is assumed to be Newtonian. The physical model and geometrical coordinates are shown in Figure 1.

Figure 1.

Schematic diagram of the problem.

The MHD body forces J¯×B¯ the Maxwell’s equations:

divB¯=0,CurlB¯=μmJ¯anddivE¯=0

where J¯ is the electric current density, B¯=B+b is the total magnetic field, μm is the magnetic permeability and b is the induced magnetic field. The external uniform magnetic field B is imposed in the direction normal to the surface of the disk which is assumed unchanged by taking small magnetic Reynolds number, so that the flow induction distortion of the applied magnetic field can be neglected as in the case with most of conducting fluids. In addition, a uniform suction is applied at the surface of the disk for the entire range.

The magnetic body force J¯×B¯ takes the form σV¯×B¯×B¯, therefore, σV¯×B¯×B¯=σB2V¯/r, where σ is the electrical conductivity of fluid and V¯ is velocity vector V¯=uvw and B¯=00B. The Lorentz force (MHD body force) has two components:

Fr=σB2u/r,Fθ=σB2v/r.

Under these assumptions, the governing equations for the continuity, momentum, energy and concentration in laminar incompressible flow can be written as follows:

ur+ur+wz=0,E1
uurv2r+wuz=1ρpr+ν2ur2+1rurur2+2uz2σB2ρruνK1u+gβTTT+gβCCC,E2
uvr+uvr+wvz=ν2vr2+1rvrvr2+2vz2σB2ρrvνK1v,E3
uwr+wwz=1ρpz+ν2wr2+1rwr+2wz2,E4
uTr+wTz=kρcp2Tr2+1rTr+2Tz21ρcpqrz+DkTcscp2Cr2+1rCr+2Cz2+σB2ρcpru2+QρcpTT,E5
uCr+wCz=D2Cr2+1rCr+2Cz2+DkTTm2Tr2+1rTr+2Tz2k1CCE6

When the mean free path of the fluid particle is comparable to the characteristic dimensions of the flow field domain, Naiver-Stokes equations break down since the assumption of continuum media fails. In the range of 0.1<kn<10 Knudsen Number, the higher order continuum equation, for example, Burnett equation should be used. For the range of 0.001kn0.1, no slip boundary conditions cannot be used and should be replaced with the following expression (Gad-el-Hak [20]):

Ut=2ξξλUtn,E7

Where Ut is the tangential velocity, n is the normal direction to the wall, ξ is the tangent momentum accommodation coefficient and λ is the mean free path. For kn<0.001, the no-slip boundary condition is valid, therefore, the velocity at the surface is equal to zero. In this study the slip and the no-slip regimes of the Knudsen number that lies in the range 0.1>kn>0 is considered. By using Eq. (7), the appropriate boundary conditions for the flow induced by an infinite disk z=0 which rotates with constant angular velocity Ω subjected to uniform suction w0 through the disk are given by

z=0:u=2ξξλuz,v=Ωr+2ξξλvz,w=w0,T=Tw,C=Cw,z:u0,v0,TT,CC,pp.E8

where k is the thermal conductivity, ν=μ/ρ is the kinematic viscosity of the ambient fluid, σ is the electrical conductivity, K1 is the permeability of the porous medium, g is the gravitational acceleration, βTandβC are the expansion coefficients of temperature and concentration, ρ,μandcp are the density, dynamic viscosity and the specific heat at constant pressure, respectively, Q is the volumetric heat generation/absorption rate, D is the molecular diffusion coefficient, kT,cs,Tmandqr are the thermal-diffusion rate, concentration susceptibility, the mean fluid temperature and the radiative heat flux. Using the Rosseland approximation (Rashed [21]), the radiative heat flux qr could be expressed by

qr=4σ3kT4zE9

Where the σ represents the Stefan-Boltzman constant and k is the Rosseland mean absorption coefficient.

Assuming that the temperature difference within the flow is sufficiently small such that T4 could be approached as the linear function of temperature

T44T3T3T4E10

For the flow under study, it is relevant to assume that the applied magnetic field Br has the form cobble [22] B=B0r, where B0 is constant magnetic flux density.

To obtain the non-dimensional form of the above equations, the following dimensionless variables are introduced.

R¯=rL,Z¯=zL,U¯=uΩL,V¯=vΩL,λ¯=λL,W¯=wΩL,P¯=ppρΩ2L2,ν¯=νΩL2,k1¯=k1L2,T¯=TTTwT,C¯=CCCwC,k1=k¯L2E11

Substituting Eqs. (9)(11) in Eqs. (1)(6), we obtain the following dimensionless equations:

U¯R¯+U¯R¯+W¯Z¯=0,E12
U¯U¯R¯V2¯R¯+W¯U¯Z¯=P¯R¯+ν¯2U¯R2¯+1R¯U¯R¯U¯R2¯+2U¯Z2¯σB02ρΩU¯ν¯k¯1U¯+gβTTwTΩ2LT¯+gβcCwCΩ2LC¯,E13
U¯V¯R¯U¯V¯R¯+W¯V¯Z¯=ν¯2V¯R2¯+1R¯V¯R¯V¯R2¯+2V¯Z2¯σB02ρΩV¯ν¯k¯1V¯,E14
U¯W¯R¯+W¯W¯Z¯=P¯Z¯+ν¯2W¯R2¯+1R¯W¯R¯+2W¯Z2¯,E15
U¯T¯R¯+W¯T¯Z¯=kρcpΩL22T¯R2¯+1R¯T¯R¯+2T¯Z2¯+16σT33ρcpΩkL22T¯Z2¯+DkTCwCcscpΩL2TwT2C¯R2¯+1R¯C¯R¯+2C¯Z2¯+σB02ΩLρcpTwTU2¯R¯+QρcpΩT¯,E16
U¯C¯R¯+W¯C¯Z¯=DΩL22C¯R2¯+1R¯C¯R¯+2C¯Z2¯+DkTTwTTmΩL2CwC2T¯R2¯+1R¯T¯R¯+2T¯Z2¯k¯ΩL2C¯.E17

The boundary conditions (8) are reduced to

Z¯=0:U¯=2ξξλ¯U¯Z¯,V¯=R¯+2ξξλ¯V¯Z¯,W¯=w0ΩL,T¯=1,C¯=1,Z¯:U¯0,V¯0,T¯0,C¯0,P¯0.E18

The governing equations are obtained by introducing a dimensionless normal distance from the disk, η=Z¯/ν¯ along with the von-Karman transformations,

U¯=R¯Fη,V¯=R¯Gη,W¯=ν¯Hη.T¯=θη,C¯=φη,P¯=ν¯P¯ηE19

Where F,G,H,θ,φandP¯ are non-dimensionless functions in terms of vertical co-ordinate η. Substituting these transformations into Eqs. (12)(17) gives the nonlinear ordinary differential equations, expressed as

2F+H/=0,E20
F//HF/F2+G2M+SF+αθ+=0,E21
G//HG/2FGM+SG=0,E22
H//HH/P/¯=0,E23
1Pr1+43Rdθ//Hθ/+Duφ//+JF2+δθ=0,E24
1Scφ//Hφ/+S0θ//βφ=0.E25

With the appropriate boundary conditions:

η=0;F0=γF/0,G0=1+γG/0,H0=Ws,θ0=1,φ0=1,η;F=0,G=0,θ=0,φ=0,P¯=0.E26

Where γ=2ξλΩ/ξν the slip is factor and Ws=w0/νΩ represents a uniform suction Ws<0 at the disk surface. The boundary conditions given by Eq. (26) imply that the radial F, the tangential G components of velocity, temperature and concentration vanish sufficiently far away from the disk, whereas the axial velocity component H is anticipated to approach a yet unknown asymptotic limit for sufficiently large η-values.

Advertisement

3. Skin-friction coefficient, Nusselt number and Sherwood number

The parameters of engineering interest for the present problem are the local skin-friction coefficients and the local rates of heat and mass transfer to the surface are calculated. The radial shear stress and tangential shear stress are given by:

τφz=μvz+1rwφz=0,
τzr=μuz+wrz=0,

The tangential and radial skin-friction coefficient are, respectively, given by, Eq. (27):

Cf1=τφz/ρΩ2L2=ReG/0,Cf2=τzr/ρΩ2L2=ReF/0,E27

Now the heat flux qw and the mass flux Jw at the surface are given by.

qw=kTzz=0,andMw=DCzz=0.

Hence the Nusselt number, Nu and Sherwood number, Sh are obtained as Eqs. (28) and (29), respectively:

Nu=LqwkTwTν¯1/2=θ/0,E28

And

Sh=LMwDCwCν¯1/2=φ/0.E29

Where Re=R2/¯ν¯ is the rotational Reynolds number.

Advertisement

4. Numerical results and discussion

The system of non-linear ordinary differential Eqs. (20)(25) together with the boundary conditions (26) are locally similar and solved numerically by using the Keller-box method. In order to gain physical insight, the velocity (radial, axial and tangential), temperature, concentration and pressure profiles have been discussed by assigning numerical values to the parameter, encountered in the problem which the numerical results are tabulated and displayed with the graphical illustrations.

In order to verify the accuracy of our present method, we have compared our results with those of Frusteri and Osalusi [15], Osalusi and Sibanda [23] and Maleque and Sattar [24]. Table 1 shows the values of F/0,G/0andθ/0 for various values of Ws. The comparisons in all above cases are found to be excellent and agreed, so it is good.

WsF|(0)–G|(0)–θ|(0)
Frusteri et al. [15]Osalusi et al. [23]Maleque et al. [24]Present studyFrusteri et al. [15]Osalusi et al. [23]Maleque et al. [24]Present studyFrusteri et al. [15]Osalusi et al. [23]Maleque et al. [24]Present study
γ = 0, ε = 0
0.00.42410.42410.51020.46270.65140.65140.61600.62410.53870.53870.32580.3171
−2.00.23240.23240.24250.24352.06872.06872.03912.04521.51961.51961.44211.4746
−4.00.12460.12460.12480.12454.00654.00654.00544.00522.85202.85202.84472.8432
−5.00.09990.09990.09990.09985.00295.00295.00315.00323.55413.55413.55413.5539
−10.00.04990.04990.05060.050210.000310.000310.001710.00197.10027.10027.10207.1018

Table 1.

Comparison of the current and recent numerical values of the radial and tangential skin-friction coefficients and the rate of heat transfer coefficient for various of Ws with S=0.0,α=0.0,N=0.0,J=0.0,δ=0.0,Pr=0.71,β=0.0,Rd=109,Du=0.0,S0=0.0,Sc=1.0andM=0.

Figures 2, 3 and 11ad display the velocity (radial, axial and tangential), temperature, concentration and pressure profiles under the magnetic field parameter, porosity parameter and Schmidt number. The (radial, axial and tangential) components of the velocity and pressure profile decrease with increase of magnetic field due to the inhibiting influence of the Lorentz force and increasing of all porosity parameter and Schmidt number, while the temperature and the concentration profiles increase with increasing of all magnetic field parameter, porosity parameter and Schmidt number. In Figures 4 and 5ad, it is clear that the (radial and axial) components of the velocity and pressure profiles increase with increase of temperature buoyancy parameter and increasing concentration buoyancy parameter, while the tangential component of the velocity, the temperature and the concentration profiles decrease with increasing temperature buoyancy parameter and concentration buoyancy parameter.

Figure 2.

Effect of magnetic field parameter on (a) the velocity (radial, axial and tangential) profile, (b) the temperature profile, (c) the concentration and (d) the pressure profile.

Figure 3.

Effect of porosity parameter on (a) the velocity (radial, axial and tangential) profile, (b) the temperature profile, (c) the concentration and (d) the pressure profile.

Figure 4.

Effect of temperature buoyancy parameter on (a) the velocity (radial, axial and tangential) profile, (b) the temperature profile, (c) the concentration and (d) the pressure profile.

Figure 5.

Effect of concentration buoyancy parameter on (a) the velocity (radial, axial and tangential) profile, (b) the temperature profile, (c) the concentration and (d) the pressure profile.

The effects of Pr on the (radial and axial) components of the velocity, the temperature, concentration and pressure profiles are shown in Figure 6ad, respectively. It is observed that both the (radial and axial) components of the velocity, the temperature and pressure profiles decrease with the increase of Prandtl number. While, the concentration profile increase. Physically, it means that thermal boundary layer thickness gets decreased. In fact, it is well known that the thermal boundary layer thickness is inversely proportional to the square root of Prandtl number. Hence, the decrease of temperature profile with increasing Pr is straightforward.

Figure 6.

Effect of Prandtl number on (a) the velocity (radial and axial) profile, (b) the temperature profile, (c) the concentration and (d) the pressure profile.

Figures 7 and 10ad show the effects of radiation parameter and heat source parameter on the velocity (radial and axial), temperature, concentration and pressure profiles, also, we found that the (radial and axial) components of the velocity, temperature and pressure profiles increase with the increase of radiation parameter and increasing heat source parameter. While the concentration profile decrease. The effects of DuandS0 on the velocity (radial and axial), temperature, concentration and pressure profiles, are shown in Figure 8ad, respectively. It is observed that the (radial and axial) components of the velocity, temperature, concentration and pressure profiles increase with the decreasing the Dufour’s number and increasing Soret number.

Figure 7.

Effect of radiation parameter on (a) the velocity (radial and axial) profile, (b) the temperature profile, (c) the concentration and (d) the pressure profile.

Figure 8.

Effect of Soret and Dufour’s number on (a) the velocity (radial and axial) profile, (b) the temperature profile, (c) the concentration and (d) the pressure profile.

Figure 9.

Effect of Joule heating parameter on (a) the velocity (radial and axial) profile, and (b) the pressure profile.

Figure 10.

Effect of heat source parameter on (a) the velocity (radial and axial) profile, (b) the temperature profile, (c) the concentration and (d) the pressure profile.

The effects of J on the velocity (radial and axial) and pressure profiles, are shown in Figure 9a and b, respectively. We found that the (radial and axial) components of the velocity and pressure profiles decrease with the increase Joule heating parameter. In Figure 12ac, it is clear that the (radial and axial) components of the velocity, concentration and pressure profiles decrease with increase of chemical reaction parameter. Figure 13ad displays the velocity (radial, axial and tangential), temperature, concentration and pressure profiles under the suction parameter. The (radial, axial and tangential) components of the velocity, temperature, concentration and pressure profiles decrease with increase of suction parameter. The effects of γ on the velocity (radial, axial and tangential) and pressure profiles are shown in Figure 14a and b, respectively. It is observed that the (radial and axial) components of the velocity, and pressure profiles increase with the increasing slip parameter. While the tangential component of the velocity profile decrease.

Figure 11.

Effect of Schmidt number on (a) the velocity (radial, axial and tangential) profile, (b) the temperature profile, (c) the concentration and (d) the pressure profile.

Figure 12.

Effect of chemical reaction parameter on (a) the velocity (radial and axial) profile, (b) the concentration and (c) the pressure profile.

Figure 13.

Effect of suction parameter on (a) the velocity (radial, axial and tangential) profile, (b) the temperature profile, (c) the concentration and (d) the pressure profile.

Figure 14.

Effect of slip parameter on (a) the velocity (radial, axial and tangential) profile, and (b) the pressure profile.

The radial and tangential skin frictions and the heat and mass transfer coefficients are tabulated in Table 2 for various values of M, S, α, N, Pr, Rd, Du&S0, J, δ, Sc, β, Ws and γ. We observed that increase for all magnetic field parameter M and porosity parameter S leads to an decrease in the all tangential skin friction G/0, heat transfer rate θ/0 and mass transfer rate φ/0, while an increase in the radial skin friction F/0, the increase for all radial skin friction, tangential skin friction, heat transfer rate and mass transfer rate, with increasing of the temperature buoyancy parameter α and the concentration buoyancy parameter N. We found that the radial skin friction, tangential skin friction, mass transfer rate decreases while heat transfer rate increase with increasing of Prandtl number, Dufour number decreases and Soret number increases. It can that be seen that the radial skin friction, tangential skin friction and mass transfer rate increase while heat transfer rate decrease with increasing values of Rd and δ. It is observed that an increase in the Joule heating parameter, results in a decrease in the tangential Skin-friction coefficient, Nusselt number and Sherwood number. The tangential skin friction and heat transfer rate decrease but the radial skin friction and mass transfer rate increase with increasing the Schmidt number. It also can be seen from this table that increasing the chemical reaction parameter to decrease in the radial skin friction, tangential skin friction and heat transfer rate while increase the mass transfer rate. We found also the tangential skin friction increase but the radial skin friction, Nusselt number and Sherwood number decrease with increasing the suction parameter. Finally, the radial skin friction and the tangential skin friction decrease but Nusselt number and Sherwood number increase with increasing slip parameter.

MSαNPrRdDu&S0JδScβWsγRe1l2Cf1Re1l2Cf2NuSh
0.12.00.50.51.05.00.25 & 0.20.50.00.640.5−1.00.20.38571.62570.31441.0486
2.00.33231.72750.31161.0367
3.00.31121.77320.31051.0319
0.10.20.55341.39780.32201.0870
1.00.46291.50710.31851.0663
2.00.38571.62570.31441.0486
0.12.00.00.1970−1.60050.30391.0056
0.50.38571.62570.31441.0486
0.80.49531.63970.31981.0725
0.12.00.50.00.25101.60960.30931.0247
0.50.38571.62570.31441.0486
0.90.48941.63770.31801.0667
0.12.00.50.50.710.38791.62610.29441.0522
1.00.38571.62570.31441.0486
2.30.37571.62410.41551.0326
0.12.00.50.51.01.00.36831.62290.50221.0194
5.00.38571.62570.31441.0486
10.00.38921.62630.28281.0541
0.12.00.50.51.05.01.0 & 0.050.39101.62660.23241.0591
0.25 & 0.20.38571.62570.31441.0486
0.2 & 0.250.38561.62560.31951.0454
0.12.00.50.51.05.00.25 & 0.20.10.38580.31571.0487
0.50.38570.31441.0486
1.50.38540.31111.0485
0.12.00.50.51.05.00.25 & 0.20.50.50.37951.62470.39381.0370
0.00.38571.62570.31441.0486
0.50.39311.62700.22391.0625
0.12.00.50.51.05.00.25 & 0.20.50.00.220.60831.41210.34510.5590
0.640.38571.62570.31441.0486
0.780.37771.62780.30841.2094
0.12.00.50.51.05.00.25 & 0.20.50.00.640.20.39271.62680.31930.9130
0.50.38571.62570.31441.0486
1.00.37701.62440.30751.2368
0.12.00.50.51.05.00.25 & 0.20.50.00.640.5−1.00.38571.62570.31441.0486
−1.50.37371.76840.34451.2854
2.00.35901.90950.37641.5404
0.12.00.50.51.05.00.25 & 0.20.50.00.640.5−1.00.20.38571.62570.31441.0486
0.50.24031.09680.31561.0573
0.80.17570.82640.31621.0613

Table 2.

Numerical of the values skin-friction coefficient Re1l2Cf1Re1l2Cf2, Nusselt number Nu and Sherwood number Sh at the surface with M, S, α, N, Pr, Rd, Du&S0, J, δ, Sc, β, Wsand γ.

Advertisement

5. Conclusions

In this work, thermal radiation and thermal diffusion effects over an electrically conducting, Newtonian fluid in a steady laminar magnetohydrodynamic convective flow over a porous rotating infinite disk with the consideration of heat and mass transfer in the presence of Soret and Dufour’s diffusion effects have been obtained and studied numerically. Magnetic field parameter, porosity parameter, temperature buoyancy parameter, concentration buoyancy parameter, Prandtl number, radiation parameter, Soret and Dufour’s number, Joule heating parameter, heat source parameter, Schmidt number, chemical reaction parameter, suction parameter, slip parameter effects were considered in the separate cases. The subsequent outcome may be drawn as:

  1. The components (radial, axial and tangential) of the velocity, temperature, concentration and pressure profiles under the porosity parameter. The (radial, axial and tangential) components of the velocity and pressure profile decrease with increasing porosity parameter, while the temperature and the concentration profiles increase with increasing porosity parameter. And also, that increase porosity parameter S leads to a decrease in the all tangential skin friction G/0, heat transfer rate θ/0 and mass transfer rate φ/0, while an increase in the radial skin friction F/0.

  2. The (radial and axial) components of the velocity and pressure increase with increasing of temperature buoyancy parameter and concentration buoyancy parameter, while the tangential component of the velocity, the temperature and the concentration profiles decrease with increasing temperature buoyancy parameter and concentration buoyancy parameter. We found that the increase for all radial skin friction, tangential skin friction, heat transfer rate and mass transfer rate, with increasing of the temperature buoyancy parameter α and the concentration buoyancy parameter N.

  3. The (radial and axial) components of the velocity, temperature and pressure profiles increase with the increase of heat source parameter. While the concentration profile decrease. And also, the radial skin friction, tangential skin friction and mass transfer rate increase while heat transfer rate decrease with increasing value of δ.

  4. The (radial and axial) components of the velocity, concentration and pressure profiles decrease with increase of chemical reaction parameter. And in Table 2, an increase the chemical reaction parameter, results in an decrease in the (radial and tangential) Skin-friction coefficient and Nusselt number, while Sherwood number increases.

Advertisement

Nomenclature

B

external uniform magnetic field

B0

constant magnetic flux density

b

induced magnetic field

C

concentration distribution

Cw

uniform concentration

C

constant concentration

cp

specific heat at constant pressure

cs

concentration susceptibility

Cf1

tangential skin-friction coefficient

Cf2

radial skin-friction coefficient

D

molecular diffusion coefficient

Du

Dufour number =DkTCwC/νcscpTwT

g

gravitational acceleration, ms2

k

thermal conductivity

k1

rate of chemical reaction

KT

thermal-diffusion rate

K1

permeability of the porous medium

k

Rosseland mean absorption coefficient

F,G,H

radial F, tangential G and axial H components of dimensionless velocity

M

magnetic field parameter =σB02/ρΩ

N

concentration buoyancy parameter =gβcCwC/LR¯Ω2

Nu

Nusselt number

n

normal direction to the wall

J

electric current density

J

Joule heating parameter =σB02ΩL/ρcpTwT

Pr

Prandtl number =νρcp/k

p

pressure distribution, Nm2

p

constant pressure

Q

volumetric heat generation/absorption rate, Wm3

qr

radiative heat flux, Wm2

Rd

radiation parameter =kk/4σT3

S

porosity parameter =ν¯/k1¯

S0

Soret number =DkTTwT/νTmCwC

Sc

Schmidt number =ν/D

Sh

Sherwood number

T

temperature distribution K

Tm

mean fluid temperature K

Tw

uniform temperature K

T

constant temperature K

uvw

components of the flow velocity are in the directions of increasing rφz, respectively, ms1

Ut

tangential velocity, ms1

w0

uniform suction

Greek symbolsμm

magnetic permeability, kgm1s1

ξ

tangent momentum accommodation coefficient

λ

mean free path

ρ

fluid density, kgm3

μ

dynamic viscosity, kgm1s1

ν

kinematic viscosity of the ambient fluid, m2s1

δ

heat source parameter =Q/Ωρcp

α

temperature buoyancy parameter =gβTTwT/LR¯Ω2

γ

slip factor

η

similarity variable

θ

dimensionless temperature distribution

φ

dimensionless concentration distribution

σ

electrical conductivity

σ

Stephan-Boltzmann constant

Ω

constant angular velocity

β

chemical reaction parameter =k¯/ΩL2

βT

coefficient of temperature

βC

coefficient of concentration

Subscriptsw,

conditions at the surface and in the free stream

References

  1. 1. Von Karman T. Uber laminare und turbulente Reibung. Zeitschrift für Angewandte Mathematik und Mechanik. 1921;1:1233-1255
  2. 2. Cochran WG. The flow due to a rotating disk. Proceedings of the Cambridge Philosophical Society. 1934;30:365-375
  3. 3. Benton ER. On the flow due to a rotating disk. Journal of Fluid Mechanics. 1966;24:781-800
  4. 4. Nield DA, Bejan A. Convection in Porous Media. 2nd ed. New York: Springer-Verlag; 1999
  5. 5. Eckert ERG, Drake RM. Analysis of Heat and Mass Transfer. New York: McGraw-Hill; 1972
  6. 6. Northrup EF. Some newly observed manifestations of forces in the interior of an electrical conductor. Physics Review. 1907;24:474-497
  7. 7. Tillack MS, Morley NB. Magnetohydrodynamics. Standard Handbook for Electrical Engineer. McGraw Hill; 1998
  8. 8. Hayat T, Nadeem S, Siddiqui A, Asghar S. An oscillating hydromagnetic non-Newtonian flow in a rotating system. Applied Mathematics Letters. 2004;17:609-614
  9. 9. Ahn CH, Allen MG. Fluid micropumps based on rotary magnetic actuators. Proceedings of MEMS. 1995;95:408-412
  10. 10. Mensing G, Pearce T, Bebe DJ. An in-plane active magnetic mixer for microfluidic applications. In: IEEEEMB Conference. 2002. pp. 531-534
  11. 11. Mahmud S, Tasnim SH, Mamun MAH. Thermodynamic analysis of mixed convection in a channel with transverse hydromagnetic effect. International Journal of Thermal Sciences. 2003;42:731-740
  12. 12. Miklavcic M. The flow due to a rough rotating disk. Zeitschrift für Angewandte Mathematik und Physik. 2004;55:235-246
  13. 13. Mizzi S, Emerson DR, Stefanov SK, Barber RW, Reese JM. Effects of rarefaction on cavity flow in the slip regime. Journal of Computational and Theoretical Nanoscience. 2007;4:817-822
  14. 14. Sparrow EM, Beavers GS, Hung LY. Flow about a porous surface rotating disk. International Journal of Heat and Mass Transfer. 1971;14:993-996
  15. 15. Frusteri F, Osalusi E. On MHD and slip flow over a rotating porous disk with variable properties. International Communications in Heat and Mass Transfer. 2007;14:34, 492-501
  16. 16. Sparrow EM, Minkowycz WJ, Eckert ERG. Transpiration induced buoyancy and thermal diffusion-diffusion thermo in helium-air free convection boundary layer. ASME Journal of Heat Transfer. 1964;86:508-520
  17. 17. Abdul Maleque K. Magnetohydrodynamic convective heat and mass transfer flow due to a rotating disk with thermal diffusion effect. Journal of Heat Transfer. 2009;131:1-8
  18. 18. Afify Ahmed A. Effects of thermal diffusion and diffusion thermo on free convective heat and mass transfer over a stretching surface considering suction or injection. Communications in Nonlinear Science and Numerical Simulation. 2009;14:2202-2214
  19. 19. Beg A, Bakier AY, Prasad VR. Numerical study of free convection magnetohydrodynamic heat and mass transfer from a stretching surface to a saturated porous medium with Soret and Dufour effects. Computational Materials Science. 2009;46:57-65
  20. 20. Gad-el-Hak M. The fluid mechanics of microdevices. The freeman scholar lecture. Journal of Fluids Engineering-Transactions of the ASME. 1999;121:5-33
  21. 21. Rashed GMA. Chemical entropy generation and MHD effects on the unsteady heat and fluid flow through a porous medium. Journal of Applied Mathematics. 2016;2016:9 pages, Article ID 1748312
  22. 22. Cobble MH. Magnetohydrodynamic flow with a pressure gradient and fluid injection. Journal of Engineering Mathematics. 1977;11:249-256
  23. 23. Osalusi E, Sibanda P. On variable laminar convective flow properties due to a porous rotating disk in a magnetic field. Romanian Journal of Physics. 2006;51(9–10):933-944
  24. 24. Maleque AK, Sattar AM. Steady laminar convective flow with variable properties due to a porous rotating disk. Journal of Heat Transfer. 2005;127:1406-1409

Written By

Gamal M. Abdel-Rahman and Faiza M.N. El-fayez

Submitted: 07 September 2018 Reviewed: 22 October 2018 Published: 22 January 2019