Applied Mathematics
Vol.07 No.18(2016), Article ID:72925,20 pages

Combined Effects of Thermal Diffusion and Diffusion-Thermo Effects on Transient MHD Natural Convection and Mass Transfer Flow in a Vertical Channel with Thermal Radiation

Isah Bala Yabo1, Basant Kumar Jha2, Jeng-Eng Lin3

1Department of Mathematics, Usmanu Danfodiyo University, Sokoto, Nigeria

2Department of Mathematics, Ahmadu Bello University, Zaria, Nigeria

3Department of Mathematical Sciences, George Mason University, Virginia, USA

Copyright © 2016 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

Received: October 15, 2016; Accepted: December 19, 2016; Published: December 22, 2016


We study the nonlinear coupled evolution equations which model the transient MHD natural convection and mass transfer flow of viscous, incompressible and electrically conducting fluid between two infinite vertical plates in the presence of the transversal magnetic field, thermal radiation, thermal diffusion and diffusion-thermo effects. Both analytical and numerical methods are used for this study.


Transient MHD, Diffusion-Thermo, Thermal Diffusion

1. Introduction

Heat and mass transfer flow of an electrically conducting fluid between parallel plates in the presence of a magnetic field and thermal radiation is of a special significance both naturally and in many industrial applications such as magnetohydrodynamic power generators, pumps, cooling of nuclear reactors, geothermal systems, thermal insulators, nuclear waste disposal, petroleum and polymer technology, heat exchangers and others [1] . One of the earliest studies in this field was carried out by [2] . They analyzed the influence of the effects of a transverse uniform magnetic field on the flow of a viscous incompressible electrically conducting fluid exiting through parallel stationary plates that are insulated. Since then, several aspects of this problem have been investigated [3] .

It has been found that an energy flux can be generated not only by temperature gradients but also by concentration gradients. The heat transfer caused by concentration gradient is called the diffusion-thermo or Dufour effect while mass transfer caused by temperature gradients is called Soret or thermal diffusion effect. Thus, Soret effect is referred to species differentiation developing in an initial homogenous mixture submitted to a thermal gradient and the Dufour effect referred to the heat flux produced by a concentration gradient. Fouriers neglect these effects on the basis that they are of a smaller order of magnitude, than the effects described and Ficks laws in heat and mass transfer processes [4] . [5] found that the thermal diffusion effect assumes significance for instance in isotope separation and in mixtures between gases with very light molecular weight (H2, He) and for medium molecular weight (N2, air), the diffusion-thermo effects constitute a substantial effect such that it cannot be ignored. Hence several models with Soret and Dufour effect in different heat and mass transfer problems have been studied. Recently, [6] considered Dufour effect on the free-convection and mass transfer flow in a vertical channel when the boundaries are subjected to symmetric concentration and thermal input. [7] presented a theoretical treatment of unsteady hydromagnetic flow and heat and mass transfer of an incompressible electrically-conducting and radiating fluid in a vertical filled with porous medium taking into account the Soret number. [4] considered the effect of thermo-diffusion on MHD mixed convective heat and mass transfer flow of a viscous fluid through a porous medium with radiation, heat generation and chemical reaction.

The role of thermal radiation on MHD flow and heat transfer problem continues to have great deal of interest. For instance at high operating temperature, radiation effect can be quite significant. For instance, in many processes in engineering areas such as Nuclear power plants, gas turbines and the various propulsion devices for aircraft, missiles, satellites and space vehicles, knowledge of thermal radiation effect becomes very important for the design of the pertinent equipment [8] . [9] explained the behaviors of chemically reacted unsteady MHD free convection flow in the presence of thermal diffusion and diffusion thermo. [10] analyzed Soret and Dufour effects on unsteady MHD mixed convection flow past a radiative vertical porous plate embedded in a porous medium with chemical reaction. [11] discussed the effect of Soret and Dufour number on an unsteady magnetohydrodynamic free convective fluid flow past a vertical porous plate in the presence of suction or injection. [9] studied chemical reaction effect on an unsteady MHD free convection flow past an infinite vertical accelerated plate with constant heat flux, thermal diffusion and diffusion thermo. [12] examined the thermal radiation effect on the MHD three-dimensional flow of Eyring-power fluid. More recently, [13] carried out a parametric study on radiation, Soret and Dufour effects in MHD channel flow bounded by a long wavy wall and a uniformly moving parallel at wall.

In the present article, we consider the thermal diffusion and diffusion-thermo effect on MHD free convection and mass transfer flow which, will be analysed by means of non-linear Rosseland diffusion approximation.

2. Mathematical Analysis

Consider a two-dimensional transient combined free convective and mass transfer flow of a viscous, incompressible and electrically conduction fluid between two infinite vertical parallel plates. A uniform transverse magnetic field of magnitude is applied in the presence of an incident radiation flux of intensity, which absorbed by the plate and transferred to the fluid as shown in Figure 1. At time, both the fluid and plates are assumed to be at rest at constant temperature and constant concentration respectively. At time the temperature of the plate situated at and concentration rise to and while the other plate at a distance from it, is fixed and maintained at temperature. The stream wise coordinate is denoted by taken vertically upward direction and that normal to it is denoted by. The flow is assumed laminar and fully developed means that the axial (direction) velocity depends only on transverse coordinate,. Therefore from the continuity equation, the transverse velocity, the variable is equated to zero. Since the plates are of infinite length, the velocity, temperature and concentration are function of and alone. Using the Boussinesq’s approximation, the governing equations for the present physical situation in the dimensional form are:




where is the initial temperature of the fluid and porous plates, is the dimensional temperature of the fluid, is the thermal diffusivity, is the thermal

Figure 1. Schematic diagram of the problem.

conductivity, is the density of the fluid, is the coefficient of the thermal expansion, is the fluid electrical conductivity, is the gravitational acceleration and is the strength of applied magnetic field. The flow is assumed laminar and fully developed. Meaning that the axial (direction) velocity depends only on transverse coordinate,.


The quantity appearing on the right hand side of Equation (2) represents the radiative heat flux in the direction where the radiative heat flux in the direction is considered insignificant in comparison with that in the direction. The radiative heat flux term in the problem is simplified by using the Rosseland diffusion approximation for an optically thick fluid according to [14] , where is Stefan-Boltzmann constant and the mean absorption coefficient.

This approximation is valid for intensive absorption, that is, for an optically thick boundary layer. Despite these shortcomings, the Rosseland approximation has been used with success in a variety of problems ranging from the transport of radiation through gases at low density to the study of the effects of radiation on blast waves by nuclear explosion [14] .

The required initial and boundary conditions to be satisfied are


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


Using Equations ((4) and (6)) in Equations (1)-(3), we obtain the following dimensionless equations for velocity, temperature and concentration respectively.





3. Analytical Solution

The governing equations presented in the previous section are highly nonlinear and exhibited no exact solutions. In general such solution can be very useful in validating computer routines of complicated time dependent two or three-dimensional free convective and radiating conducting fluid and comparison with experimental data. It is therefore of interest to reduce the governing equations of the present problem to the form that can be solved analytically. A special case of the present problem that exhibit analytical solution is the problem of steady state MHD natural convection flow in a vertical channel in a presence of thermal radiation, thermal diffusion and diffusion-thermo effects. The resulting steady state equations and boundary conditions for this special case can be written as




the boundary conditions are


To construct analytical solutions of Equations ((11) to (13)) subject to (14), it is assumed that the radiation parameter is small and taking a power series expansion in the radiation parameter employs a regular perturbation method.


where is the radiation parameter. The second and higher order of terms of gives correction to, and account for thermal radiation effects. Substituting Equation (15) in Equations ((11), (13)) and equating like powers of, one obtains the boundary value problem for and as







The relevant boundary conditions to be satisfied are:


The required solutions of the governing steady state energy, concentration and momentum Equations ((16)-(21)) subject to boundary condition (22) are:




Using (23), we write the steady state rate of heat (Nusselt number) transfer on the boundary:


The steady state rate of mass (Sherwood number) transfer on the boundary plate from (24) as:


Also, using (25), we write the skin friction coefficient as:


where and the constant and are remove in order to reduce the size of the work.

4. Numerical Procedure

The nonlinear partial differential Equations ((7)-(9)) are solved numerically using semi-implicit finite difference scheme. We used forward difference formulas for all time derivatives and approximate both the second and first spatial derivatives with second order central differences. The semi implicit finite difference equation corresponding to Equations ((7)-(9)) is as follows:




Using the known values of and at grid point and reducing the solution to tri-diagonal matrix we obtained temperature and concentration fields at time using the known values of the previous time for all. Then the velocity field is evaluated using the already known value of temperature and concentration fields obtained at. These processes are repeated till the required solution of, and are gained at convergence criteria.


The iterative system does not restrict time step and the technique is always convergent and unconditionally stable.

5. Results and Discussion

This paper, examined the transient natural convective flow through a vertical channel in the presence of a thermal diffusion, diffusion-thermo and radiative heat transfer under the influence of a uniform magnetic field. The system of governing Equations (7)- (9) with the boundary conditions (10) is solved employing the finite difference method. The effects of the flow governing parameters on velocity, temperature, skin friction, Nusselt number and Sherwood number, have been shown using line graph as demonstrated in Figure 2 through Figure 25. These result show the variations in the velocity field, temperature distribution, skin friction, and heat transfer rate, influenced by the material parameters of the flow problem, that is, the magnetic parameter, radiation parameter, temperature difference parameter, sustention parameter, dimensionless time, Prandtl number, Schmidt number, thermal diffusion parameter and diffusion-thermo parameter. The radiation parameter is in the range of because terms associated with behave as strong heat sources and large value of leads to finite time temperature blow up (see [15] ). Prandtl number 0.71, 1.0, 3.0 and 7.0, which physically represent air, electrolyte solution, saturated liquid freon and water respectively. In air the diffusing chemical species of common interest have Schmidt number 0.1 1.0 [16] : the work consider for, for and for respectively. Moreover time is selected in the range, in order to capture the steady state situation for air and water. Besides all other parameter is choosing arbitrary.

Figures 2-4 display the velocity, temperature and concentration profiles for different values of dimensionless time when other controlling parameters are fixed. It is observed that velocity, temperature as well as concentration increases as dimensionless time increases and attains steady state.

Figure 5 indicates that high values of sustention parameter substantially increase the velocity.

Figures 6-8 depict the effect of Prandtl number on velocity, temperature and concentration. It is seen that an increase in the values of leads to the decrease in fluid velocity, temperature and concentration. This may be attributed to the fact that at higher Prandtl number, convection currents become weak and qualitatively decreases the temperature, which leads to decrease in velocity and concentration.

Figure 9 illustrates the effect of magnetic parameter on velocity. It is seen that an increase in the values of leads to the decrease in fluid velocity. This is true since magnetic parameter produce resistive force, which acts opposite direction to the fluid

Figure 2. Velocity profile when, , , and.

Figure 3. Temperature profile when, , , , and.

Figure 4. Concentration profile when, , , and.

Figure 5. Velocity profile when, , , , , , and.

Figure 6. Velocity profile when, , , , , , and.

Figure 7. Temperature profile when, , , , and.

Figure 8. Concentration profile when, , , , and.

Figure 9. Velocity profile when, , , , , , and.

motion. This is an important controlling mechanism inflow and heat transfer processes so that the finished product meet the desire quality specification.

Figure 10 show the influence of Schmidt number on velocity. The velocity is observed to descend with Schmidt number when other numerical values are treated constant.

Figure 11 and Figure 12 illustrate the disparity in both velocity and temperature due to a change in thermal-diffusion parameter when other parameters are fixed. From these Figures it is noted that the effect of is to increase both velocity and temperature respectively.

Figures 13-15 explain the effect of increasing diffusion-thermo parameter. With all other parameters constant, the velocity, temperature and concentration increases with increase in.

Figure 16 and Figure 17 illustrate the effect of temperature difference parameter on velocity and temperature profiles, respectively. It is observed that as the value of increases there is corresponding increase in the temperature of the fluid, and hence cause increase in velocity.

Figure 18 and Figure 19 inspect the influence of the radiation parameter on the velocity and temperature, respectively. It is observed that the impact of has insignificant effect on both velocity and temperature.

Figures 20-22 present variation of skin friction at for different values of, , and. The results show that the values of skin friction increases with increase of, and respectively. In addition the result show that skin friction increases with dimensionless time until steady state is attain.

Figures 23-25 present the variation of rate of heat transfer for different values of, , and at. It is observed in Figure 23 that the rate of heat transfer increases with increase in. The physical fact is that an increase in radiation adds more

Figure 10. Velocity profile when, , , , , and.

Figure 11. Velocity profile when, , , , , , and.

Figure 12. Temperature profile when, , , , and.

Figure 13. Velocity profile when, , , , , , , and.

Figure 14. Temperature profile when, , , , and.

Figure 15. Concentration profile when, , , , and.

Figure 16. Velocity profile when, , , , , , and.

Figure 17. Temperature profile when, , , , and.

Figure 18. Velocity profile when, , , , , , , and.

Figure 19. Temperature profile when, , , , , and.

Figure 20. Skin friction against at.

Figure 21. Skin friction against at.

Figure 22. Skin friction against at.

Figure 23. Nusselt number against at.

Figure 24. Nusselt number against at.

Figure 25. Nusselt number against at.

heat to the fluid and lead to an increase temperature, which correspond to the increase in rate of heat transfer. Figure 24 and Figure 25 show that the effect of and has insignificant effect to the values of rate of heat transfer. Furthermore as dimension time increases the rate of heat transfer decreases (see Figures 23-25). This may be attributed that the heat is been driven away from plate through fluid to the other plate.

6. Result Validation

In order to verify the accuracy of our results, the analytical solutions (steadystate) obtained from Equation (25) by Perturbation method was compared with that of numerical solution derived by semi-implicit finite difference from Equation (29). These computed results are tabulated in Table 1. It is interesting to observe from this Table 1 that the transient and steady state solutions results (under some limiting conditions) are in very good agreement at large value of time, which clearly shows the correctness of our numerical (computed) scheme.

7. Conclusions

Transient MHD natural convection flow in a vertical channel in a presence of thermal radiation, thermal diffusion and diffusion-thermo effects confined between two infinite vertical parallel plates has been investigated. The model representing the present physical situation is highly nonlinear due to the presence of thermal radiation effect. The non-linear differential equations under appropriate initial and boundary condition are numerically solved using implicit finite difference method. Also analytical solutions are derived by regular perturbation method for steady state situation. The impact of the essential dimensionless controlling parameter on velocity, temperature, concentration, skin-friction and Nusselt number is presented with the aid of line graphs and discussed.

Table 1. Comparison between numerical values of transient velocity obtained using finite difference method and steady state velocity obtained by Perturbation method.

From the indicated results of the problem, the following observations were concluded:

1) As the radiation, temperature difference, sustention parameter, thermal-diffusion, diffusion-thermo, and non-dimensional time parameters increase, both the velocity and temperature increase while velocity decreases with Schmidt number and magnetic parameter.

2) Increasing Prandtl number reduces values of velocity, temperature and concentration.

3) Skin friction, Nusselt and concentration increase with time.

4) The values of Nusselt number are higher in case of radiation and contrast in case of thermal-diffusion and diffusion-thermo at.

5) Good agreement between steady state and transient solution was found after some sufficiently large time.

6) During numerical computation, it is found that the time required reaching steady state velocity and temperature is directly proportional to the Prandtl number of the working fluid. That is the time required to reach steady state velocity and temperature is approximately 0.71 (dimensionless time) for air for fixed values of other controlling parameters.


We thank the Editor and the referee for their comments. This support is greatly appreciated.

Cite this paper

Yabo, I.B., Jha, B.K. and Lin, J.-E. (2016) Combined Effects of Thermal Diffusion and Diffusion-Thermo Effects on Transient MHD Natural Convection and Mass Transfer Flow in a Vertical Channel with Thermal Radiation. Applied Mathematics, 7, 2354-2373.


  1. 1. Moreau, R. (1990) Magnetohydrodynamics. Academic Publishers, Dordrecht.

  2. 2. Hartman, J. and Lakatos, F. (1937) Kongelige danke videnskabernes selskab. Matematisk-Fysiske Meddelelser, 15, 6-7.

  3. 3. Makinde, O.D. and Sibanda, P. (2008) Magnetohydrodynamic Mixed Convective Flow and Heat and Mass Transfer past a Vertical Plate in a Porous Medium with Constant Wall Suction. ASME-Journal of Heat Transfer, 130, 11260.

  4. 4. Kesavaraya, D.C., Satyanarayana, P.V. and Venkataramana, S.V. (2013) Radiation and Thermo-Diffusion Effects on Mixed Convective Heat and Mass Transfer Flow of a Viscous Dissipated Fluid over a Vertical Surface in the Presence of Chemical Reaction with Heat Source. International Journal of Scientific Engineering and Technology, 2, 56-72.

  5. 5. Eckert, E.R.G. and Dranke, R.M. (1972) Analysis of Heat and Mass Transfer. McGrow-Hill, New York.

  6. 6. Jha, B.K., Samaila, A.K. and Ajibade, A.O. (2011) Diffusion-Thermo Effects on Free Convective Heat and Mass Transfer Flow in a Vertical Channel with Symmetric Boundary Conditions. Journal of Heat Transfer, 133, 1-8.

  7. 7. Kumar, R. (2012) Thermo Diffusion and Radiative Heat Effect on Unsteady Oscillatory Mixed Convection MHD Flow through a Planer Channel Filled with Porous Medium. International Journal of Physics and Mathematical Sciences, 3, 13 p.

  8. 8. Ibrahim, F.S.M.E.A. and Bakr, A.A. (2006) Effect of the Chemical Reaction and Radiation Absorption on the Unsteady MHD Free Convection Flow past a Semi Infinite Vertical Permeable Moving Plate with Heat Source and Suction. Communication in Nonlinear Science and Numerical Simulation, 13, 1056-1066.

  9. 9. Sudhakar, K., Raju, R.S. and Rangamma, M. (2012) Chemical Reaction Effect on an Unsteady MHD Free Convection Flow past an Infinite Vertical Accelerated Plate with Constant Heat Flux, Thermal Diffusion and Diffusion Thermo. International Journal of Modern and Engineering Research, 2, 3329-3339.

  10. 10. Sharma, B.K., Yazad, K., Mishra, N.K. and Chaudhary, R.C. (2012) Soret and Dufour Effects on Unsteady MHD Mixed Convection Flow past a Radiative Vertical Porous Plate Embedded in a Porous Medium with Chemical Reaction. Applied Mathematics, 3, 717-723.

  11. 11. Sarada, K. and Shankar, B. (2013) The Effects of Soret and Dufour on an Unsteady MHD Free Convection Flow past a Vertical Porous Plate in the Presence of Suction or Injection. International Journal of Engineering Science, 2, 13-25.

  12. 12. Hayat, T., Awais, M. and Asghar, S. (2013) Radiation Effects in a Three-Dimensional Flow of MHD Eyring-Powell Fluid. Journal of the Egyptian Mathematical Society, 21, 379-384.

  13. 13. Ahmed, N., Bhaattacharya, D.J. and Barua, D.P. (2014) Radiation, Soret and Dufour Effects in MHD Channel Flow Bounded by a Long Wavy Wall and a Uniformly Moving Parallel at Wall. Journal of Applied Mathematics and Fluid Mechanics, 6, 1-19.

  14. 14. Ali Agha, H., Bouaziz, M.N. and Hanini, S. (2014) Free Convection Boundary Layer Flow from a Vertical Flat Plate Embedded in a Darcy Medium Filled with a Nanofluid: Effects of Magnetic Field and Thermal Radiation. Arabian Journal for Science and Engineering, 39, 8331-8340.

  15. 15. Makinde, O.D. and Chinyoka, T. (2010) Numerical Investigation of Transient Heat Transfer to Hydromagnetic Channel Flow with Radiative Heat and Convective Cooling. Communication in Nonlinear Science and Numerical Simulation, 15, 3919-3930.

  16. 16. Jha, B.K. and Ajibade, A.O. (2009) Free Convective Flow of Heat Generating/Absorbing Fluid between Vertical Porous Plates with Periodic Heat Input. International Communications in Heat and Mass Transfer, 36, 624-631.