# Effects of Thermal Diffusion and Radiation on Unsteady MHD Free Convection Flow Past an Infinite Heated Vertical Plate in a Porous Medium

Authors: S. Sivaiah; K. Anitha; S. Venkataramana
Journal: ISRN Thermodynamics
 (2012)

---

## Abstract

We have explored the influence of thermal diffusion and radiation on unsteady magnetohydrodynamic free convection flow past an infinite heated vertical plate in a porous medium. The governing boundary layer equations are written into a dimensionless form by similarity transformations. The transformed nonlinear differential equations are solved numerically with finite element methods. Numerical calculations are carried out for different values of dimensionless parameters. The results are presented graphically for velocity, temperature, and concentration profiles and show that the flow field and other quantities of physical interest are significantly influenced by these parameters.


---

## Body

## 1. Introduction

Flow through porous medium past infinite vertical plate is common in nature and has many applications in engineering and science. A number of workers have investigated such flows and excellent literature on the properties and phenomenon may be found in the literature [1–4]. For example, Soundalgekar [3] investigated the effects of free convection currents on the oscillatory type boundary layer flow past an infinite vertical plate with constant suction where the plate temperature differs from the free steam temperature.Recently, attention has been on the effects of transversely applied magnetic field and thermal perturbation on the flow of electrically conducting viscous fluids such as plasma. Various properties associated with the interplay of magnetic fields and thermal perturbation in porous medium past vertical plate find useful applications in astrophysics, geophysical fluid dynamics, and engineering. Researches in these fields have been conducted by many investigators [5–9]. For example, Soundalgekar [10] investigated a two-dimensional steady free convection flow of an incompressible, viscous, and electrically conducting fluid past an infinite vertical porous plate with constant suction and plate temperature when the difference between the plate temperature and free stream is moderately large to cause free convection currents. In another study Israel-Cookey and Sigalo [11] investigated the problem of unsteady MHD past a semi-infinite vertical plate in an optically thin environment with simultaneous effects of radiation, free convection parameters, and time-dependent suction. Chamka [12] investigated the unsteady convective heat and mass transfer past an infinite permeable moving plate with heat absorption where it was found that increase in Solutal Grashoff number enhanced the concentration buoyancy effects leading to an increase in the velocity. Anand Rao and Sivaiah [13] studied the chemical reaction effects on an unsteady MHD free convective flow past an infinite vertical porous plate with constant suction and heat source. Anand Rao et al. [14, 15] studied finite element solution of MHD transient flow past an impulsively started infinite horizontal porous plate in a rotating fluid with hall current and also studied heat and mass transfer in MHD flow of a viscous fluid past a vertical plate under oscillatory suction velocity.The objective of the present chapter is to examine the effects of thermal diffusion and radiation on unsteady MHD free convection flow past an infinite heated vertical plate in a porous medium. The equations of continuity, linear momentum, energy, and diffusion, which govern the flow field, are solved by using Galerkin finite element method. Similarity solutions are then obtained numerically for various parameters, entering into the problem, and discussed from the physical point of view.


## 2. Formulation of the Problem

We consider the unsteady flow of an incompressible, viscous, and radiating hydromagnetic fluid past an infinite porous heated vertical plate with time-dependent suction in an optically thin environment. The physical model and the coordinate system are shown in Figure1. The x′-axis is taken along the vertical infinite porous plate in the upward direction and the y′-axis normal to the plate.Figure 1
The physical model and coordinate system of the problem.At timet′=0, the plate is maintained at a temperature Tw′, which is high enough to initiate radiative heat transfer. A constant magnetic field H0′2 is maintained in the y′ direction and the plate moves uniformly along the positive x′ direction with velocity U0. Under Boussinesq approximation the flow is governed by the following equations:
(1)∂w′∂y′=0,(2)∂u′∂t′+w′∂u′∂y′=υ∂2u′∂y′2+∂U′∂t′-(μ2σCH0′2ρ+υK)(u′-U′)+gβ(T′-T∞′)+gβ*(C′-C∞′),(3)∂T′∂t′+w′∂T′∂y′=kρCP(∂2T′∂y′2-∇qz′′)+μρCP(∂u′∂y′)2-QoρCP(T′-T∞′),(4)∂2qz′′∂y′2-3α2qz′′-16ασT∞3∂T′∂y′=0,(5)∂C′∂t′+w′∂C′∂y′=D∂2C′∂y′2+DmkTTm∂2T′∂y′2.The boundary conditions are(6)u′=0,T′=Tw′,C′=Cw′ony′=0,u′=0,T′=T∞′,C′=C∞′asy′→∞,
where (u′,v′,w′) are dimensional velocity components; (x′,y′) is dimensional cartesian coordinates; H0′2 is constant transverse magnetic field; w0′ is dimensional suction velocity; Tw′ is wall temperature; T∞′ is reference temperature; U′ is dimensional free stream velocity; t′ is dimensional time; Cw′ is concentration near the plate; C∞′ is concentration in the fluid for away from the plate the plate; ρ is density; g is acceleration due to gravity; K is dimensional porosity parameter; Cp is specific heat capacity; k is thermal conductivity; qz′ is radiative heat flux; U0 is mean velocity of U′(t′); σC is electrical conductivity; Q0 is dimensional heat absorption coefficient; υ is kinematic coefficient; μ is permeability; ε is time corrective parameter; β is coefficients of volume expansion due to temperature; β* is coefficient of volume expansion due concentration; D is chemical diffusivity; Dm is molecular diffusivity; kT is mean absorption coefficient; Tm is mean fluid temperature; α2 is absorption coefficient; ω′ is dimensional free stream frequency of oscillation.Since the medium is optically thin with relatively low density andα≪1, the radiative heat flux given by (4) in the spirit of Bestman [7] becomes
(7)∂qz′′∂y′=4α2(T′-T∞′),
where
(8)α2=∫0∞δλ∂B∂T′.
Here B is the Planck’s constant, δ is the radiation absorption, and λ is the frequency.                                                     Further, from (1) it is clear that the suction velocity w′ at the plate is either a constant or a function of time only. So we assume it in the form
(9)w′=-w0′(1+εAeiω′t′),
where A is small positive parameter and ε is small such that εA≪1, and the negative sign indicates that the suction velocity is towards the plate.In view of (4), (8), and (9), (2), (3), and (5) become
(10)14∂u∂t-(1+εAeiωt)∂u∂y=14∂U∂t+∂2u∂y2-(M2+χ2)(u-U)+Grθ+GcC,(11)14Pr∂θ∂t-Pr(1+εAeiωt)∂θ∂y=(∂2∂y2-R2)θ+PrEc(∂u∂y)2-PrSθ,(12)14Sc∂C∂t-Sc(1+εAeiωt)∂C∂y=∂2C∂y2+(Sr)(Sc)(∂2θ∂y2).In order to write the governing equations and the boundary conditions in a dimensionless form, the following nondimensional quantities are introduced:(13)t=w0′2t′4υ,y=w0′y′υ,u=u′U0,ω=4υω′w0′2,U=U′U0,θ=T′-T∞′Tw′-T∞′,Pr=μcPk,Gr=gβυ(Tw′-T∞′)U0w0′2,Gc=gβ*υ(Cw′-C∞′)U0w0′2,Ec=U02CP(Tw′-T∞′),R2=4α2ρCPkw0′2(Tw′-T∞′),S=υQ0ρcPw0′2,M2=μ2σCH0′2υρw0′2,Sr=DmkT(Tw′-T∞′)υTm(Cw′-C∞′),Sc=υD,C=(C′-C∞′)(Cw′-C∞′),χ2=υ2Kw0′2,
where M2 is Hartmann number; Pr is Prandtl number; Sc is Schmidt number; Ec is Eckert number; θ is nondimensional temperature; C is nondimensional species concentration; Gr is Grashof number; Gc is modified Grashof number; χ2 is Darcy number; Sr is Soret number; R2 is radiation parameter; S is heat source.Equations (10), (11), and (12) are now subject to the boundary conditions
(14)u=0,θ=1,C=1ony=0,u→0,θ→0,C→0asy→∞.
The mathematical statement of the problem is now complete and embodies the solution of (10), (11), and (12) subject to boundary conditions (14).


## 3. Solution of the Problem

By applying Galerkin finite element method for (10) over the element (e) (yj≤y≤yk) is
(15)∫yjyk{N(e)T[4∂2u(e)∂y2-∂u(e)∂t+4B∂u(e)∂y-Du(e)+P]}dy=0,
where P=(∂U/∂t)+4(Gr)θ+4(Gc)C+DU, B=1+εAeiωt, and D=4(M2+χ2).Integrating the first term in (15) by parts, one obtains
(16)4N(e)T{∂u(e)∂y}yjyk-∫yjyk{4∂N(e)T∂y∂u(e)∂y+N(e)T(∂u(e)∂t-4B∂u(e)∂y+Du(e)-P)∂N(e)T∂y}dy=0.
Neglecting the first term in (16), one gets
(17)∫yjyk{4∂N(e)T∂y∂u(e)∂y+N(e)T(∂u(e)∂t-4B∂u(e)∂y+Du(e)-P)∂N(e)T∂y}dy=0.
Let u(e)=N(e)ϕ(e) be the linear piecewise approximation solution over the element (e) (yj≤y≤yk), where N(e)=[NjNk], ϕ(e)=[ujuk]T, Nj=(yk-y)/(yk-yj), and Nk=(y-yj)/(yk-yj) are the basis functions. One obtains
(18)4∫yjyk{[Nj′Nj′Nj′Nk′Nj′Nk′Nk′Nk′][ujuk]}dy+∫yjyk{[NjNjNjNkNjNkNkNk][u˙ju˙k]}dy-4B∫yjyk{[NjNj′NjNk′Nj′NkNk′Nk][ujuk]}dy+D∫yjyk{[NjNjNjNkNjNkNkNk][ujuk]}dy=P∫yjyk[NjNk]dy.
Simplifying, we get
(19)4l(e)2[1-1-11][ujuk]+16[2112][u˙ju˙k]-4B2l(e)[-11-11][ujuk]+D6[2112][ujuk]=P2[11],
where prime and dot denote differentiation with respect to y and time t, respectively. Assembling the element equations for two consecutive elements (yi-1≤y≤yi) and (yi≤y≤yi+1), one obtains
(20)4l(e)2[1-10-12-10-11][ui-1uiui+1]+16[210141012][u˙i-1u˙iu˙i+1]-4B2l(e)[-110-1010-11][ui-1uiui+1]+D6[210141012][ui-1uiui+1]=P2[121].Now put row corresponding to the nodei to zero, from (20) the difference schemes with l(e)=h are
(21)4h2[-ui-1+2ui-ui+1]+16[u˙i-1+4u˙i+u˙i+1]-4B2h[-ui-1+ui+1]+D6[ui-1+4ui+ui+1]=P.
Applying the trapezoidal rule, the following system of equations in Crank-Nicholson method is obtained
(22)A1ui-1n+1+A2uin+1+A3ui+1n+1=A4ui-1n+A5uin+A6ui+1n+P*.
Now from (11) and (12) the following equations are obtained:
(23)B1θi-1n+1+B2θin+1+B3θi+1n+1=B4θi-1n+B5θin+B6θi+1n+Q*,C1Ci-1n+1+C2Cin+1+C3Ci+1n+1=C4Ci-1n+C5Cin+C6Ci+1n+R**,
where
(24)A1=2+12Brh+Dk-24r,A2=8+4Dk+48r,A3=2-12Brh+Dk-24r,A4=2-12Brh-Dk+24r,A5=8-4Dk-48r,A6=2+12Brh+Dk+24r,B1=2(Pr)+12(Pr)Brh+4R2k-24r+4S(Pr)k,B2=8(Pr)+48r+16R2k+16S(Pr)k,B3=2(Pr)-12(Pr)Brh+4R2k-24r+4S(Pr)k,B4=2(Pr)-12(Pr)Brh-4R2k+24r-4S(Pr)k,B5=8(Pr)-48r-16R2k-16S(Pr)k,B6=2(Pr)+12(Pr)Brh-4R2k+24r-4S(Pr)k,C1=2(Sc)+12(Pr)Brh-24r,C2=8(Sc)+48r,C3=2(Sc)-12(Sc)Brh-24r,C4=2(Sc)-12(Sc)Brh+24r,C5=8(Sc)-48r,C6=2(Sc)+12(Sc)Brh+24r,P*=12Pk=12k((∂U/∂t)+4(Gr)θ+4(Gc)C+DU),Q*=12Qk=48(Pr)k(Ec)(∂u/∂y)2,R**=12R*k=48(Sc)(Sr)k(∂2θ/∂y2).Herer=k/h2 and h, k are mesh sizes along y-direction and time-direction, respectively. Index i refers to space and j refers to the time. In (22) and (23), taking i=1(1)n and using boundary conditions (14), then the following system of equations is obtained
(25)AiXi=Bii=1(1)3,
where Ai′s are matrices of order n and Xi, Bi′s are column matrices having n-components. The solutions of the above system of equations are obtained by using Thomas algorithm for velocity, temperature, and concentration. Also, numerical solutions for these equations are obtained by C-programme. In order to prove the convergence and stability of Galerkin finite element method, the same C-programme was run with smaller values of h and k; no significant change was observed in the values of u, θ, and C. Hence the Galerkin finite element method is stable and convergent.


## 4. Results and Discussion

In the previous sections, we have formulated and solved the problem of an unsteady MHD free convection flow past an infinite heated vertical plate in a porous medium with radiation. By invoking, the optically thin differential approximation for the radiative heat flux in the energy equation. In the numerical computation, the Prandtl number (Pr=0.71) which corresponds to air and various values of the material parameters are used. In addition, the boundary condition y→∞ is approximated by ymax=2, which is sufficiently large for the velocity to approach the relevant stream velocity. The temperature and the species concentration are coupled to the velocity via Grashof number (Gr) and modified Grashof number (Gc) as seen in (9). For various values of Grashof number and modified Grashof number, the velocity profiles u are plotted in Figures 2 and 3. The Grashof number (Gr) signifies the relative effect of the thermal buoyancy force to the viscous hydrodynamic force in the boundary layer. As expected, it is observed that there is a rise in the velocity due to the enhancement of thermal buoyancy force. Also, as (Gr) increases, the peak value of the velocity increases rapidly near the porous plate and then decays smoothly to the free stream velocity. The modified Grashof number (Gc) defines the ratio of the species buoyancy force to the viscous hydrodynamic force. As expected, the fluid velocity increases and the peak value is more distinctive due to the increase in the species buoyancy force. The velocity distribution attains a distinctive maximum value in the vicinity of the plate and then decreases properly to approach the free stream value. It is noticed that the velocity increases with the increasing values of modified Grashof number (Gc). Figure 4 illustrates the velocity profiles for different values of Prandtl number Pr. The numerical results show that the effect of increasing values of Prandtl number results in decreasing velocity. The nature of velocity profiles in the presence of foreign species such as hydrogen (Sc = 0.22), helium (Sc = 0.30), water vapour (Sc = 0.60), and oxygen (Sc = 0.66) is shown in Figure 5. The flow field suffers a decrease in primary velocity at all points in the presence of heavier diffusing species. The effect of the magnetic field parameter M is shown in Figure 6 in case of cooling of the plate. It is observed that the velocity of the fluid decreases with the increase of the magnetic field parameter values. The decrease in the velocity as the Hartmann number M increases is because the presence of a magnetic field in an electrically conducting fluid introduces a force called the Lorentz force, which acts against the flow if the magnetic field is applied in the normal direction, as in the present study. This resistive force slows down the fluid velocity component as shown in Figure 6. The influence of the viscous dissipation parameter, that is, the Eckert number Ec on the velocity, and temperature are shown in Figures 7 and 13, respectively. The Eckert number Ec expresses the relationship between the kinetic energy in the flow and the enthalpy. It embodies the conversion of kinetic energy into internal energy by work done against the viscous fluid stresses. Greater viscous dissipative heat causes a rise in the temperature as well as the velocity.Figure 2
Velocity profiles for different values ofGr.Figure 3
Velocity profiles for different values ofGc.Figure 4
Velocity profiles for different values ofPr.Figure 5
Velocity profiles for different values ofSc.Figure 6
Velocity profiles for different values ofM.Figure 7
Velocity profiles for different values ofEc.The effect of the thermal radiation parameterR on the primary velocity and temperature profiles in the boundary layer is illustrated in Figures 8 and 15, respectively. Increasing the thermal radiation parameter R produces significant increase in the thermal condition of the fluid and its thermal boundary layer. This increase in the fluid temperature induces more flow in the boundary layer causing the velocity of the fluid there to increase. Figures 9 and 14 have been plotted to depict the variation of velocity and temperature profiles against y for different values of heat source parameter S by fixing other physical parameters. From this graph we observe that velocity and temperature decrease with increase in the heat source parameter S because when heat is absorbed, the buoyancy force decreases the temperature profiles. Figure 10 shows the effects of Darcy number χ on the velocity profiles for cooling as well as heating the plate. For a cooling plate, fluid velocity increases, whereas for a heating plate it decreases with increase of χ. Darcy number is the measurement of the porosity of the medium. As the porosity of the medium increases, the value of χ increases. For large porosity of the medium, fluid gets more space to flow, as a consequence its velocity increases. The variations of tangential velocity distribution with y for different values of the Soret number (Sr) are shown in Figure 11. It can be clearly seen that the velocity distribution in the boundary layer increases with the Soret number. Figure 12 illustrates the temperature profiles for different values of Prandtl number Pr. It is observed that the temperature decreases as the Prandtl number increases. The reason is that smaller values of Pr are equivalent to increase in the thermal conductivity of the fluid and therefore heat is able to diffuse away from the heated surface more rapidly for higher values of Pr. Hence in the case of smaller Prandtl number the thermal boundary layer is thicker and the rate of heat transfer is reduced. The effects of Schmidt number (Sc) and Soret number (Sr) on the concentration field are presented in Figures 16 and 17. Figure 16 shows that the concentration field is due to variation in Schmidt number (Sc) for the gasses hydrogen, helium, water vapour, oxygen, and ammonia. It is observed that the concentration field is steady for hydrogen and falls rapidly for oxygen and ammonia in comparison to water vapour. Thus hydrogen can be used for maintaining effective concentration field and water vapour can be used for maintaining normal concentration field. In Figure 17, it is observed that an increase in the Soret number (Sr) leads to increase in the concentration field.Figure 8
Velocity profiles for different values ofR.Figure 9
Velocity profiles for different values ofS.Figure 10
Velocity profiles for different values ofχ.Figure 11
Velocity profiles for different values ofSr.Figure 12
Temperature profiles for different values ofPr.Figure 13
Temperature profiles for different values ofEc.Figure 14
Temperature profiles for different values ofS.Figure 15
Temperature profiles for different values ofR.Figure 16
Concentration profiles for different values ofSc.Figure 17
Concentration profiles for different values ofSr.


## 5. Conclusions

In conclusion, therefore, the flow of an unsteady MHD free convection past an infinite heated vertical plate in a porous medium under the simultaneous effects of thermal diffusion and radiation is affected by the material parameters. The governing equations are approximated to a system of linear partial differential equations by using Galerkin finite element method. The results are presented graphically and we can conclude that the flow field and the quantities of physical interest are significantly influenced by these parameters.(1)
The velocity increases as Grashof numberGr, modified Grashof number Gc, Eckert number, thermal radiation parameter R, Darcy parameter χ, Ec, and Soret number Sr increase. However, the velocity was found to be decreasing as the Hartmann number M, Prandtl number Pr, Schmidt number Sc, and heat source parameter S are increasing.(2)
The fluid temperature was found to be decreasing as the heat source parameterS and Prandtl number Pr are increasing and found to be increasing as Eckert number Ec, and thermal radiation parameter R are increasing.(3)
The fluid concentration was found to be decreasing as the Schmidt numberSc and increasing as the Soret number Sr is increasing.



---

*Source: 670396-2012-11-28.xml*