Abstract
In this article, we explore the use of a nonperturbation linearization method to solve the coupled highly nonlinear system of equations due to flow over a vertical surface subject to a magnetic field. The linearization method is used in combination with an asymptotic expansion technique. The effects of Dufour, Soret and magnetic filed parameters are investigated. The velocity, temperature and concentration distributions as well as the skinfriction, heat and mass transfer coefficients have been obtained and discussed for various physical parametric values. The accuracy of the solutions has been tested using a local nonsimilarity method. The results show that the nonperturbation technique is an accurate numerical algorithm that converges rapidly and may serve as a viable alternative to finite difference and finite element methods for solving nonlinear boundary value problems.
Mathematics Subject Classification (2000): 76D05; 74S30; 76E06; 76M25.
Keywords:
crossdiffusion; free convection; linearization method; incompressible flow; magnetohydrodynamics (MHD)1 Introduction
Convection driven by density variations caused by two different components which have different rates of diffusion plays an important role in fluid dynamics since such flows occur naturally in many physical and engineering processes. Heat and salt in sea water provide perhaps the best known example of doublediffusive convection, Stern [1]. Other examples of double diffusive convection are encountered in diverse applications such as in chemical and petroleum industries, filtration processes, food processing, geophysics and in the modeling of solar ponds and magma chambers. A review of the literature in this subject can be found in Nield and Bejan [2].
One of the earliest studies of double diffusive convection was by Nield [3]. Baines and Gill [4] investigated linear stability boundaries, while Rudraiah et al. [5] used the nonlinear perturbation theory to investigate the onset of double diffusive convection in a horizontal porous layer. Poulikakos [6] presented the linear stability analysis of thermosolutal convection using the DarcyBrinkman model. Bejan and Khair [7] presented a multiple scale analysis of heat and mass transfer about a vertical plate embedded in a porous medium. They considered concentration gradients which aid or oppose thermal gradients. Related studies on double diffusive convection have been undertaken by, among others, Lai [8], Afify [9], and Makinde and Sibanda [10].
Investigations by, among others, Eckert and Drake [11] and Mortimer and Eyring [12] have provided examples of flows such as in the geosciences, where diffusionthermo and thermaldiffusion effects are quite significant. Anjalidevi and Devi [13] showed that diffusionthermo and thermaldiffusion effects are significant when density differences exist in the flow regime. In general, Diffusionthermo and thermaldiffusion effects have been found to be particularly important for intermediate molecular weight gases in binary systems that are often encountered in chemical engineering processes. Theoretical studies of the Soret and Dufour effects on double diffusive convection have been made by many researchers, among them, Kafoussias and Williams [14], Postelnicu [15], Mansour et al. [16], Narayana and Sibanda [17], and Awad et al. [18].
In this article, we investigate convective heat and mass transfer along a vertical flat plate in the presence of diffusionthermo, thermal diffusion effects, and an external magnetic field. The governing momentum, heat and mass transfer equations are, in general, strongly coupled and highly nonlinear. Apart from numerical methods, a number of semianalytical techniques have in recent years been proposed to find approximate solutions of nonlinear boundary value problems. Some of these techniques, such as the Adomian decomposition method, the variational iteration method, the homotopy analysis method, the homotopy perturbation method, the differential transform method, etc, are now well known and their strengths and weaknesses well understood. Some of these weaknesses include, for example, small regions of convergence and the use of artificially inserted parameters (Liao [19], Geng [20]). The recent spectral homotopy analysis method (see Motsa et el. [21,22]) sought to improve the accuracy and efficiency of the homotopy analysis method while retaining its essence. The use of spectral methods also provided greater flexibility in the choice of basis functions. Nonetheless, the challenge to find more accurate, robust and computationally efficient solution techniques for nonlinear problems in engineering and science still remains.
A recent method that remains to be generalized and whose robustness remains to be tested in the case of highly nonlinear equations with a strong coupling is the successive linearization method (Makukula et al. [23,24]). This method has been used in a limited number of studies by, for example, Awad et al. [25]) and Motsa et al. [26] to solve fluid flow problems. In this study the coupled set of differential equations that describe convective heat and mass transfer flow along a vertical flat plate in the presence of diffusionthermo, thermal diffusion effects and an external magnetic field are solved using the successive linearization method. A nonsimilarity technique is used to validate the linearization method.
2 Problem formulation
We consider the problem of double diffusive convection along a vertical plate with an external magnetic field imposed along the yaxis. The induced magnetic field is assumed to be negligible. The fluid temperature and solute concentration in the ambient fluid are T_{∞}, C_{∞ }and those at the surface are T_{w }and C_{w}, respectively. The coordinates system and the flow configuration are shown in Figure 1.
Figure 1. Physical model and coordinate system.
Under the usual boundary layer and Boussinesq approximations the governing equations for a viscous incompressible fluid may be written as (Kafoussias and Williams [14]);
subject to the boundary conditions
where u and v are the velocity components along the x and y axes, respectively T and C are the fluid temperature and solute concentration across the boundary layer, ν is the kinematic viscosity, ρ is the fluid density, σ is the electrical conductivity, B_{0 }is the uniform magnetic field, β_{T }and β_{C }are the coefficients of thermal and solutal expansions, D_{m }is the thermal diffusivity, k_{T }is the thermal diffusion ratio, c_{s }is the concentration susceptibility, c_{p }is the fluid specific heat capacity, T_{m }is the mean fluid temperature, U_{∞ }is the free stream velocity and g is the gravitational acceleration.
To satisfy the continuity equation (1), we define the stream function ψ in terms of the velocity by
and introduce the following dimensionless variables
where, without loss of generality, we take the constant length ℓ to be unity in the subsequent analysis. Using equations (7) in (2)(4), we get the transformed equations
with corresponding boundary conditions
The fluid and physical parameters in equations (8)(10) are the local thermal and solutal Grashof numbers Gr_{x }and Gc_{x}, the local magnetic field parameter Ha_{x}, the Prandtl number Pr, the Dufour number D_{f}, the Soret number S_{r}, and the Schmidt number Sc. These parameters are defined in equations (12), (13) below.
The parameters of engineering interest in heat and mass transport problems are the skin friction coefficient C_{fx}, the Nusselt number Nu_{x }and the Sherwood number Sh_{x}. These parameters characterize the surface drag, the wall heat and mass transfer rates, respectively, and are defined by
and
where Re_{x }= U_{∞}x/ν.
3 Method of solution
The main challenge in using the linearization method as described in Makukula et al. [23,24] is how to generalize the method so as to find solutions of partial differential equations of the form (8)(10). It is certainly not clear how the method may be applied directly to the terms on the right hand side of equations (8)(10). For this reason, equations (8)(10) are first simplified and reduced to sets of ordinary differential equations by assuming regular perturbation expansions for f, θ, and ϕ in powers of ξ (which is assumed to be small) as follows
where M_{1 }is the order of the approximate solution. Substituting (17) into equations (8)(10) and equating the coefficients of like powers of ξ, we obtain the zeroth order set of ordinary differential equations
with corresponding boundary conditions
The O(ξ^{1}) equations are
with boundary conditions
Finally, the O(ξ^{2}) equations are
These equations have to be solved subject to boundary conditions
The coupled system of equations (18)(20), (22)(24), and (26)(28) together with the associated boundary conditions (21), (25), and (29) may be solved independently pairwise one after another. These equations may now be solved using the successive linearization method in the manner described in Makukula et al. [23,24]). We begin by solving equations (18)(20) with boundary conditions (21).
The successive linearization method is a nonperturbation method requiring neither the presence of an embedded perturbation parameter nor the addition of an artificial parameter. The method is therefore free of the major limitations associated with other perturbation methods. In the SLM algorithm assumption is made that the functions f_{0}(η), θ_{0}(η), and ϕ_{0}(η) may be expressed as
where F_{i}, Θ_{i}, and Φ_{i }(i ≥ 1) are unknown functions and F_{m}, Θ_{m}, and Φ_{m }are successive approximations which are obtained by recursively solving the linear part of the system that is obtained from substituting equations (30) in (18)(20). In choosing the form of the expansions (30), prior knowledge of the general nature of the solutions, as is often the case with perturbation methods, is not necessary.
Suitable initial guesses F_{0}(η), Θ_{0}(η), and Φ_{0}(η) which are selected to satisfy the boundary conditions (21) are
The subsequent solutions F_{i}, Θ_{i}, and Φ_{i }are obtained by iteratively solving the linearized form of the equations that are obtained by substituting equation (30) in the governing equations (18)(20). The linearized equations to be solved are
subject to the boundary conditions
where the coefficients parameters a_{k, i1}, b_{k,i1}, c_{k,i1}, d_{k,i1}, and r_{k,i1 }are defined by
Once each solution F_{i}, Θ_{i}, and Φ_{i }(i ≥ 1), has been found from iteratively solving equations (32)(34), the approximate solutions for the system (18)(20) are obtained as
where M_{2 }is the order of the SLM approximations. In this study, we used the Chebyshev spectral collocation method to solve equations (32)(34). The physical region [0, ∞) is first transformed into the spectral domain [1,1] using the domain truncation technique in which the problem is solved on the interval [0, L], where L is a scaling parameter used to invoke the boundary condition at infinity. This is achieved by using the mapping
We discretise the spectral domain [1,1] using the GaussLobatto collocation points given by
where N is the number of collocation points used. The unknown functions F_{i}, Θ_{i }, and Φ_{i }are approximated at the collocation points by
where T_{k }is the kth Chebyshev polynomial defined as
The derivatives of the variables at the collocation points are represented as
where r is the order of differentiation and
Substituting equations (37)(41) into (32)(34) gives the following linear system of equations
subject to the boundary conditions
Here A_{i1 }is a 3(N + 1) × 3(N + 1) square matrix, while X_{i }and R_{i1 }are 3(N + 1) × 1 column vectors defined by
where
and
where [·] is a diagonal matrix of size (N + 1) × (N + 1) and a_{k,i1}, b_{k,i1}, c_{k,i1 }are diagonal matrices of size (N + 1) × (N + 1) and T is the transpose. After modifying the matrix system (43) to incorporate the boundary conditions (44), the solution is obtained as
Equation (53) gives a solution of (18)(20) for f_{0}, θ_{0 }and ϕ_{0}. The procedure is repeated to obtain the O(ξ^{1}) and O(ξ^{2}) solutions using equations (22)(24) and (26)(28), respectively.
4 Results and discussion
In generating the results in this article, we determined through numerical experimentation that L = 15, N = 40, and M_{2 }= 5 gave sufficient accuracy for the linearization method. The value of the Prandtl number used is Pr = 0.71 which physically corresponds to air. The Schmidt number used Sc = 0.22 is for hydrogen at approximately 25° and one atmospheric pressure (see Afify [27]). In concert with previous related studies, the Dufour and Soret numbers are chosen in such a way that their product is constant, provided the mean temperature T_{m }is also kept constant.
To determine the accuracy and validate the linearization method, equations (8)(10) were further solved using a local nonsimilarity method (LNSM) developed by Sparrow and Yu [28] and Sparrow et al. [29]. Previous studies have consistently used the Matlab bvp4c solver to evaluate the accuracy of the successive linearization method. However, as with other BVP solvers, the accuracy and convergence of the bvp4c algorithm depends on a good initial guess and works better for systems involving few equations (Shampine et al. [30]).
Tables 1, 2, and 3 show, firstly the effects of various parameters on the skinfriction and the local heat and mass transfer coefficients at different values of ξ and, second, give a sense of the accuracy and convergence rate of the linearization method. The results from the two methods are in excellent agreement with the second order SLM series giving accuracy of up to five significant figures.
Table 1. The effect of the magnetic field parameter Ha_{x }on f"(0), θ' (0) and ϕ'(0) when Gr_{x }= 0.5, Gc_{x }= 2, D_{f }= 0.2, and Sr = 0.3
Table 2. The effect of the Soret parameter Sr on f"(0), θ'(0) and ϕ'(0) when Gr_{x }= 0.5, Gc_{x }= 2, and Ha_{x }= 1
Table 3. Soret and Dufour effects of the skin friction coefficient C_{f}, Nusselt number Nu and Sherwood number Sh when Gr_{x }= 0.5, Gc_{x }= 2, Sc = 0.22, and Ha_{x }= 0.5
Table 1 shows the effect of increasing the magmatic field parameter Ha_{x }on the local skin friction, heat and mass transfer coefficients. We observe that increasing the magnetic field parameter reduces the local skin friction as well as the heat and mass transfer coefficients. In Table 2, we present the effect of the Soret parameter on f"(0), θ'(0), and ϕ'(0) which are, respectively, proportional to the local skin friction coefficient, the local Nusselt number and Sherwood number. We observe that f"(0) and θ'(0) increase with increases in Sr, while ϕ'(0) decreases as Sr increases. These results are confirmed in Table 3. Here the Nusselt number increases as the Soret number increases, while the opposite trend occurs as the Dufour number increases. The recent study by ElKabeir [31] shows that these results may be modified by injection, suction or the presence of a chemical reaction.
The influence of the various fluid and physical parameters on the fluid properties is given qualitatively in Figures 2, 3, 4, 5, 6 and 7. Figures 2 and 3 illustrate the effect of the magnetic filed parameter on the velocity f'(η), temperature θ and concentration ϕ profiles within the boundary layer. We observe that, as expected, strengthening the magnetic field slows down the fluid motion due to an increasing drag force which acts against the flow if the magnetic field is applied in the normal direction. We also observe that the magnetic field parameter enhances the temperature and concentration profiles. The effect of broadening both the temperature and concentration distributions is to reduce the wall temperature and concentration gradients thereby reducing the heat and mass transfer rates at the wall.
Figure 2. Effect of the magmatic field parameter Ha_{x }on f'(η) when Gr_{x }= 0.5, Gc_{x}= 0.1, D_{f }= 0.2, Sr = 0.3, and ξ = 0.01.
Figure 3. Effect of the magmatic field parameter Ha_{x }on θ(η) when Gr_{x }= 0.5, Gc_{x }= 0.1, D_{f }= 0.2, Sr = 0.3, and ξ = 0.01.
Figure 4. Effect of the magmatic field parameter Ha_{x }on ϕ(η) when Gr_{x }= 0.5, Gc_{x }= 0.1, D_{f }= 0.2, Sr = 0.3, and ξ = 0.01.
Figure 5. Effect of the Soret and Dufour parameters on the velocity f'(η) when Gr_{x }= 0.5, Gc_{x }= 2.5, Ha_{x }= 0.1, and ξ = 0.01.
Figure 6. Effect of Soret parameter S_{r }on θ(η) when Gr_{x }= 0.1, Gc_{x }= 2.5, Ha_{x }= 0.1, and ξ = 0.01.
Figure 7. Effect of Soret parameter Sr on ϕ(η) when Gr_{x }= 0.1, Gc_{x }= 2.5, Ha_{x }= 0.1, and ξ = 0.01.
Figure 5 shows the effect of increasing the Soret parameter (reducing the Dufour parameter) on the fluid velocity f'(η). The fluid velocity is found to increase with the Soret parameter.
The effect of Soret parameter on the temperature within the thermal boundary layer and the solute concentration is shown in Figures 6 and 7, respectively. An increase in the Soret effect reduces the temperature within the thermal boundary layer leading to an increase in the temperature gradient at the wall and an increase in heat transfer rate at the wall. On the other hand, increasing the Soret effect increases the concentration distribution which reduces the concentration gradient at the wall. These results are similar to the earlier findings by ElKaberir [31] and Alam and Rahman [32], although the latter studies were subject to injection/suction.
5 Conclusions
In this article, we have investigated MHD and crossdiffusion effects on doublediffusive convection from a vertical flat plate in a viscous incompressible fluid. Numerical approximations for the governing equations were found using a combination of a regular perturbation expansion and the successive linearization method. The solutions were validated by using a local similarity, nonsimilarity method. We determined the effects of various parameters on the fluid properties as well as on the skinfriction coefficient, the heat and the mass transfer rates. We have shown that the magnetic field parameter enhances the temperature and concentration distributions within the boundary layer. The effect of thermodiffusion is to reduce the temperature and enhance the velocity and the concentration profiles. The diffusionthermo effect enhances the velocity and temperature profiles while reducing the concentration distribution. The skinfriction, heat and mass transfer coefficients decrease with an increase in the magnetic field strength. The skinfriction and heat transfer coefficients increase whereas the mass transfer coefficient decreases with increasing Soret numbers.
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
The problem was formulated following discussions by all three authors. Numerical simulations were carried out by AAK and FGA while the discussion was authored principally by PS. All authors read and approved the final manuscript.
Acknowledgements
The authors would like to thank the referees for helpful comments and suggestions. This study was supported by the National Research Foundation (NRF). AAK was supported by the government of Sudan.
References

Stern, ME: The 'salt fountain' and thermohaline convection. Tellus. 12, 172–175 (1960). Publisher Full Text

Nield, DA, Bejan, A: Convection in Porous Media. Springer Verlag, New York (1999)

Nield, DA: Onset of thermohaline convection in a porous medium. Water Resour Res. 4, 553–560 (1968). Publisher Full Text

Baines, PG, Gill, AE: On thermohaline convection with linear gradients. J Fluid Mech. 37, 289–306 (1969). Publisher Full Text

Rudraiah, N, Srimani, PK, Friedrich, R: Finite amplitude convection in a two component fluid saturated porous layer. Int J Heat Mass Transfer. 25, 715–722 (1982). Publisher Full Text

Poulikakos, D: Double diffusive convection in a horizontally sparsely packed porous layer. Int Commun Heat Mass Transfer. 13, 587–598 (1986). Publisher Full Text

Bejan, A, Khair, KR: Heat and mass transfer by natural convection in porous medium. Int J Heat Mass Transfer. 28, 909–918 (1985). Publisher Full Text

Lai, FC: Coupled heat and mass transfer by natural convection from a horizontal line source in saturated porous medium. Int Commun Heat Mass Transfer. 17, 489–499 (1990). Publisher Full Text

Afify, AA: MHD free convective flow and mass transfer over a stretching sheet with chemical reaction. Heat Mass Transfer. 40, 495–500 (2004)

Makinde, OD, Sibanda, P: MHD mixedconvective flow and heat and mass transfer past a vertical plate in a porous medium with constant wall suction. J Heat Transfer. 130, 1–8 (2008)

Eckert, ERG, Drake, RM: Analysis of heat and mass transfer. McGrawHill, New York (1972)

Mortimer, RG, Eyring, H: Elementary transition state theory of the Soret and Dufour effects. Proc Natl Acad Sci. 77, 1728–1731 (1980). PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Anjalidevi, SP, Devi, RU: Soret and Dufour effects on MHD slip flow with thermal radiation over a porous rotating infnite disk. Commun Nonlinear Sci. 16, 1917–930 (2011). Publisher Full Text

Kafoussias, NG, Williams, EW: Thermaldiffusion and diffusion thermo effects on mixed freeforced convective and mass transfer boundary layer flow with temperature dependent viscosity. Int J Eng Sci. 33, 1369–1384 (1995). Publisher Full Text

Postelnicu, A: Influence of a magnetic field on heat and mass transfer by natural convection from vertical surfaces in porous media considering Soret and Dufour effects. Int J Heat Mass Transfer. 47, 1467–1472 (2004). Publisher Full Text

Mansour, MA, ElAnssary, NF, Aly, AM: Effects of chemical reaction and thermal stratification on MHD free convective heat and mass transfer over a vertical stretching surface embedded in a porous media considering Soret and Dufour numbers. Chem Eng J. 145, 340–345 (2008). Publisher Full Text

LakshmiNarayana, PA, Sibanda, P: Soret and Dufour effects on free convection along a vertical wavy surface in a fluid saturated Darcy porous medium. Int J Heat Mass Transfer. 53, 3030–3034 (2010). Publisher Full Text

Awad, FG, Sibanda, P, Motsa, SS: On the linear stability analysis of a Maxwell fluid with doublediffusive convection. Appl Math Model. 34, 3509–3517 (2010). Publisher Full Text

Liao, SJ: Beyond Perturbation: Introduction to the Homotopy Analysis Method. CRC Press, Boca Raton (2003)

Geng, F: Modified variational iteration method for second order initial value problems. Appl Appl Math. 1, 73–81 (2010)

Motsa, SS, Sibanda, P, Shateyi, S: A new spectralhomotopy analysis method for solving a nonlinear second order BVP. Commun Nonlinear Sci Numer Simulat. 15, 2293–2302 (2010). Publisher Full Text

Motsa, SS, Sibanda, P, Awad, FG, Shateyi, S: A new spectralhomotopy analysis method for the MJD JefferyHamel problem. Comput. Fluids. 39, 1219–1225 (2010). Publisher Full Text

Makukula, ZG, Sibanda, P, Motsa, SS: A note on the solution of the Von Karman equations using series and chebyshev spectral methods.p. 17. Boundary Value Problems (2010)

Makukula, ZG, Sibanda, P, Motsa, SS: A novel numerical technique for twodimensional laminar flow between two moving porous walls. Math Problems Eng. 2010, 15 (2010) Article ID 528956

Awad, FG, Sibanda, P, Motsa, SS, Makinde, OD: Convection from an inverted cone in a porous medium with crossdiffusion effects. Comput Math Appl. 61, 1431–1441 (2011). Publisher Full Text

Motsa, SS, Sibanda, P, Shateyi, S: On a new quasilinearization method for systems of nonlinear boundary value problems. Math Methods Appl Sci. 34, 1406–1413 (2011). Publisher Full Text

Afify, AA: Similarity solution MHD: Effects of thermal and diffusion thermo on free convective heat and mass transfer over a stretching surface considering suction or injection. Commun. Nonlinear Sci Numer Simulat. 14, 2202–2214 (2009). Publisher Full Text

Sparrow, EM, Yu, HS: Local nonsimilarity thermal boundary layer solutions. J Heat Transfer. 93, 328–334 (1971). Publisher Full Text

Sparrow, EM, Quack, H, Boerner, CJ: Local nonsimilarity boundary layer solutions. AIAA. 8, 1936–1942 (1970). Publisher Full Text

Shampine, LF, Ketzscher, R, Forth, SA: Using AD to solve BVPs in Matlab. ACM Trans Math Softw. 31(1), 79–94 (2005). Publisher Full Text

ElKabeir, SMM: Soret and Dufour effects on heat and mass transfer by mixed convection over a vertical surface saturated porous medium with temperature dependent viscosity. Int J Numer Meth Fluids (2011)

Alam, MS, Rahman, MM: Dufour and Soret effects on mixed convection flow past a vertical porous flat plate with variable suction. Nonlinear Anal: Model Control. 11, 3–12 (2006)