Skip to main content

On cross-diffusion effects on flow over a vertical surface using a linearization method

Abstract

In this article, we explore the use of a non-perturbation 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 skin-friction, 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 non-similarity method. The results show that the non-perturbation 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.

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 double-diffusive 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 Darcy-Brinkman 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 diffusion-thermo and thermal-diffusion effects are quite significant. Anjalidevi and Devi [13] showed that diffusion-thermo and thermal-diffusion effects are significant when density differences exist in the flow regime. In general, Diffusion-thermo and thermal-diffusion 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 diffusion-thermo, 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 semi-analytical 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 diffusion-thermo, thermal diffusion effects and an external magnetic field are solved using the successive linearization method. A non-similarity 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 y-axis. 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
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]);

u x + v y = 0 ,
(1)
u u x + v u y = v 2 u y 2 - σ B 0 2 ρ u + g β T ( T - T ) + g β C ( C - C ) ,
(2)
u T x + v T y = α 2 T y 2 + D m k T c s c p 2 C y 2 ,
(3)
u C x + v C y = D m 2 C y 2 + D m k T T m 2 T y 2 ,
(4)

subject to the boundary conditions

u = 0 , v = 0 , T = T w , C = C w on y = 0 ,
(5)
u = U , T = T , C = C when y .
(6)

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, B0 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

u = ψ y and v = - ψ x ,

and introduce the following dimensionless variables

ξ = x , η = U ν ξ 1 / 2 y , ψ = ( U ν ξ ) 1 / 2 f ( ξ , η ) , θ ( ξ , η ) = T - T T w - T , ϕ ( ξ , η ) = C - C C w - C ,
(7)

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

f + 12 f f H a x f + G r x θ + G c x ϕ = ξ ( f f ξ f f ξ ) ,
(8)
1 P r θ + 1 2 f θ + D f ϕ = ξ f θ ξ - θ f ξ ,
(9)
1 S c ϕ + 1 2 f ϕ + S r θ = ξ f ϕ ξ - ϕ f ξ ,
(10)

with corresponding boundary conditions

f = 0 , f = 0 , θ = 1 , ϕ = 1 , on η = 0 , f = 1 , θ = 0 , ϕ = 0 , when η .
(11)

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.

G r x = g β T ( T w - T ) x ν 2 , G c x = g β C ( C w - C ) x ν 2 , H a x = σ B 0 2 x ρ U , P r = ν α ,
(12)
D f = D m K T ( C w - C ) α C s C p ( T w - T ) , S r = K T D m ( T w - T ) α T m ( C w - C ) , S c = ν D m .
(13)

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

C f x = μ ρ U 2 u y y = 0 = f ( ξ , 0 ) R e x ,
(14)
N u x = - x T w - T T y y = 0 = - R e x θ ( ξ , 0 ) ,
(15)

and

S h x = - x C w - C C y y = 0 = - R e x ϕ ( ξ , 0 ) ,
(16)

where Re x = Ux/ν.

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

f = f M 1 ( η , ξ ) = i = 0 M 1 ξ i f i ( η ) , θ = θ M 1 ( η , ξ ) = i = 0 M 1 ξ i θ i ( η ) , ϕ = ϕ M 1 ( η , ξ ) = i = 0 M 1 ξ i ϕ i ( η )
(17)

where M1 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

f 0 + 1 2 f 0 f 0 - H a x f 0 + G r x θ 0 + G c x ϕ 0 = 0 ,
(18)
1 P r θ 0 + 1 2 f 0 θ 0 + D f ϕ 0 = 0 ,
(19)
1 S c ϕ 0 + 1 2 f 0 ϕ 0 + S r θ 0 = 0 ,
(20)

with corresponding boundary conditions

f 0 = 0 , f 0 = 0 , θ 0 = 1 , ϕ 0 = 1 , at η = 0 f 0 = 1 , θ 0 = 0 , ϕ 0 = 0 , as η = .
(21)

The O(ξ1) equations are

f 1 + 1 2 f 0 f 1 - ( H a x + f 0 ) f 1 + 3 2 f 0 f 1 + G r x θ 1 + G c x ϕ 1 = 0 ,
(22)
1 P r θ 1 + 1 2 f 0 θ 1 - f 0 θ 1 + 3 2 θ 0 f 1 + D f ϕ 1 = 0 ,
(23)
1 S c ϕ 1 + 1 2 f 0 ϕ 1 - f 0 ϕ 1 + 3 2 ϕ 0 f 1 + S r θ 1 = 0 ,
(24)

with boundary conditions

f 1 = 0 , f 0 = 0 , θ 1 = 0 , ϕ 1 = 0 , at η = 0 f 1 = 0 , θ 1 = 0 , ϕ 1 = 0 , as η = .
(25)

Finally, the O(ξ2) equations are

f 2 + 1 2 f 0 f 2 - ( H a x + 2 f 0 ) f 2 + 5 2 f 0 f 2 + G r x θ 2 + G c x ϕ 2 = f 1 f 1 - 3 2 f 1 f 1 ,
(26)
1 P r θ 2 + 1 2 f 0 θ 2 - 2 f 0 θ 2 + 5 2 θ 0 f 2 + D f ϕ 2 = f 0 θ 1 - 3 2 f 1 θ 1 ,
(27)
1 S c ϕ 2 + 1 2 f 0 ϕ 2 - 2 f 0 ϕ 2 + 5 2 ϕ 0 f 2 + S r θ 2 = f 1 ϕ 1 - 3 2 f 1 ϕ 1 .
(28)

These equations have to be solved subject to boundary conditions

f 2 = 0 , f 2 = 0 , θ 2 = 0 , ϕ 2 = 0 , at η = 0 f 2 = 0 , θ 2 = 0 , ϕ 2 = 0 , as η = .
(29)

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 non-perturbation 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 f0(η), θ0(η), and ϕ0(η) may be expressed as

f 0 ( η ) = F i ( η ) + m = 0 i - 1 F m ( η ) , θ 0 ( η ) = Θ i ( η ) + m = 0 i - 1 Θ m ( η ) , ϕ ( η ) = Φ i ( η ) + m = 0 i - 1 Φ m ( η ) ,
(30)

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 F0(η), Θ0(η), and Φ0(η) which are selected to satisfy the boundary conditions (21) are

F 0 ( η ) = η + e - η - 1 , Θ 0 ( η ) = e - η , Φ 0 ( η ) = e - η ,
(31)

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

F i + a 1 , i - 1 F i - H a x F i + a 2 , i - 1 F i + G r x Θ i + G c x Φ i = r 1 , i - 1 ,
(32)
1 P r Θ i + b 1 , i - 1 Θ i + b 2 , i - 1 F i + D f Φ i = r 2 , i - 1 ,
(33)
1 S c Φ i + c 1 , i - 1 Φ i + c 2 , i - 1 F i + S r Θ i = r 3 , i - 1 ,
(34)

subject to the boundary conditions

F i ( 0 ) = F i ( 0 ) = F i ( ) = Θ i ( 0 ) = Θ i ( ) = Φ i ( 0 ) = Φ i ( ) ,
(35)

where the coefficients parameters ak, i-1, bk,i-1, ck,i-1, dk,i-1, and rk,i- 1are defined by

a 1 , i - 1 = 1 2 m = 0 i - 1 F m , a 2 , i - 1 = 1 2 m = 0 i - 1 F m , b 1 , i - 1 = 1 2 m = 0 i - 1 F m b 2 , i - 1 = 1 2 m = 0 i - 1 Θ m , c 1 , i - 1 = 1 2 m = 0 i - 1 F m c 2 , i - 1 = 1 2 m = 0 i - 1 Φ m , r 1 , i - 1 = H a x m = 0 i - 1 F m - 1 2 m = 0 i - 1 F m m = 0 i - 1 F m - m = 0 i - 1 F m - G r x m = 0 i - 1 Θ m - G c x m = 0 i - 1 Φ m , r 2 , i - 1 = - 1 P r m = 0 i - 1 Θ m - 1 2 m = 0 i - 1 F m m = 0 i - 1 Φ m - D f m = 0 i - 1 Φ m , r 3 , i - 1 = - 1 S c m = 0 i - 1 Φ m - 1 2 m = 0 i - 1 F m m = 0 i - 1 Φ m - S r m = 0 i - 1 Φ m ,

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

f 0 ( η ) i = 0 M 2 F i ( η ) , θ 0 ( η ) i = 0 M 2 Θ i ( η ) , ϕ 0 ( η ) i = 0 M 2 Φ i ( η ) ,
(36)

where M2 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

η L = ς + 1 2 , - 1 ς 1 ,
(37)

We discretise the spectral domain [-1,1] using the Gauss-Lobatto collocation points given by

ς j = cos π j N , j = 0 , 1 , , N ,
(38)

where N is the number of collocation points used. The unknown functions F i , Θ i , and Φ i are approximated at the collocation points by

F i ( ς ) k = 0 N F i ( ς k ) T k ( ς j ) , Θ i ( ς ) k = 0 N Θ i ( ς k ) T k ( ς j ) Φ i ( ς ) k = 0 N Φ i ( ς k ) T k ( ς j ) , j = 0 , 1 , , N
(39)

where T k is the k th Chebyshev polynomial defined as

T k ( ς ) = cos [ k cos - 1 ( ς ) ] .
(40)

The derivatives of the variables at the collocation points are represented as

F i ( r ) = k = 0 N D k j r F i ( ς k ) , Θ i ( r ) = k = 0 N D k j r Θ i ( ς k ) , Φ i ( r ) = k = 0 N D k j r Φ i ( ς k ) , j = 0 , 1 , , N
(41)

where r is the order of differentiation and D = 2 L D with D being the Chebyshev spectral differentiation matrix whose entries are defined as;

D j k = c j ( - 1 ) j + 1 c k ς j - ς k j k ; j , k = 0 , 1 , , N , D k k = - ς k 2 ( 1 - ς k 2 ) k = 1 , 2 , , N - 1 , D 00 = 2 N 2 + 1 6 = - D N N .
(42)

Substituting equations (37)-(41) into (32)-(34) gives the following linear system of equations

A i - 1 X i = R i - 1
(43)

subject to the boundary conditions

F i ( ς 0 ) = k = 0 N D 0 k F i ( ς k ) = g i ( ς 0 ) = Θ i ( ς 0 ) = Φ i ( ς 0 ) = F i ( ς N ) = g i ( ς N ) = Θ i ( ς N ) = Φ i ( ς N ) = 0 .
(44)

Here Ai- 1is a 3(N + 1) × 3(N + 1) square matrix, while X i and Ri- 1are 3(N + 1) × 1 column vectors defined by

A i - 1 = A 11 A 12 A 13 A 21 A 22 A 23 A 31 A 32 A 33 , X i = F ̃ i Θ ̃ i Φ ̃ i , R i - 1 = r 1 , i - 1 r 2 , i - 1 r 3 , i - 1
(45)

where

F ̃ i = [ F i ( ς 0 ) , F i ( ς 1 ) , , F i ( ς N - 1 ) , F i ( ς N ) ] T ,
(46)
Θ ̃ i = [ Θ i ( ς 0 ) , Θ i ( ς 1 ) , , Θ i ( ς N - 1 ) , Θ i ( ς N ) ] T ,
(47)
Φ ̃ i = [ Φ i ( ς 0 ) , Φ i ( ς 1 ) , , Φ i ( ς N - 1 ) , Φ i ( ς N ) ] T ,
(48)
r 1 , i - 1 = [ r 1 , i - 1 ( ς 0 ) , r 1 , i - 1 ( ς 1 ) , , r 1 , i - 1 ( ς N - 1 ) , r 1 , i - 1 ( ς N ) ] T
(49)
r 2 , i - 1 = [ r 2 , i - 1 ( ς 0 ) , r 2 , i - 1 ( ς 1 ) , , r 2 , i - 1 ( ς N - 1 ) , r 2 , i - 1 ( ς N ) ] T
(50)
r 3 , i - 1 = [ r 3 , i - 1 ( ς 0 ) , r 3 , i - 1 ( ς 1 ) , , r 3 , i - 1 ( ς N - 1 ) , r 3 , i - 1 ( ς N ) ] T
(51)

and

A 11 = D 3 + a 1 , i - 1 D 2 - H a x D + [ a 3 , i - 1 ] , A 12 = [ G r x ] , A 13 = [ G c x ] , A 21 = [ b 2 , i - 1 ] , A 22 = P r - 1 D 2 + b 1 , i - 1 D , A 23 = D f D 2 A 31 = [ c 2 , i - 1 ] , A 32 = S r D 2 , A 33 = S c - 1 D 2 + c 1 , i - 1 D ,
(52)

where [·] is a diagonal matrix of size (N + 1) × (N + 1) and ak,i- 1, bk,i- 1, ck,i- 1are 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

X i = A i - 1 - 1 R i - 1
(53)

Equation (53) gives a solution of (18)-(20) for f0, θ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 M2 = 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 non-similarity 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 skin-friction 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 El-Kabeir [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
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
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
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
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
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
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 El-Kaberir [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 cross-diffusion effects on double-diffusive 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, non-similarity method. We determined the effects of various parameters on the fluid properties as well as on the skin-friction 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 thermo-diffusion is to reduce the temperature and enhance the velocity and the concentration profiles. The diffusion-thermo effect enhances the velocity and temperature profiles while reducing the concentration distribution. The skin-friction, heat and mass transfer coefficients decrease with an increase in the magnetic field strength. The skin-friction and heat transfer coefficients increase whereas the mass transfer coefficient decreases with increasing Soret numbers.

References

  1. Stern ME: The 'salt fountain' and thermohaline convection. Tellus 1960, 12: 172-175. 10.1111/j.2153-3490.1960.tb01295.x

    Article  Google Scholar 

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

    Chapter  Google Scholar 

  3. Nield DA: Onset of thermohaline convection in a porous medium. Water Resour Res 1968, 4: 553-560. 10.1029/WR004i003p00553

    Article  Google Scholar 

  4. Baines PG, Gill AE: On thermohaline convection with linear gradients. J Fluid Mech 1969, 37: 289-306. 10.1017/S0022112069000553

    Article  Google Scholar 

  5. Rudraiah N, Srimani PK, Friedrich R: Finite amplitude convection in a two component fluid saturated porous layer. Int J Heat Mass Transfer 1982, 25: 715-722. 10.1016/0017-9310(82)90177-6

    Article  MATH  Google Scholar 

  6. Poulikakos D: Double diffusive convection in a horizontally sparsely packed porous layer. Int Commun Heat Mass Transfer 1986, 13: 587-598. 10.1016/0735-1933(86)90035-7

    Article  Google Scholar 

  7. Bejan A, Khair KR: Heat and mass transfer by natural convection in porous medium. Int J Heat Mass Transfer 1985, 28: 909-918. 10.1016/0017-9310(85)90272-8

    Article  MATH  Google Scholar 

  8. Lai FC: Coupled heat and mass transfer by natural convection from a horizontal line source in saturated porous medium. Int Commun Heat Mass Transfer 1990, 17: 489-499. 10.1016/0735-1933(90)90067-T

    Article  Google Scholar 

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

    Article  Google Scholar 

  10. Makinde OD, Sibanda P: MHD mixed-convective flow and heat and mass transfer past a vertical plate in a porous medium with constant wall suction. J Heat Transfer 2008, 130: 1-8.

    Article  Google Scholar 

  11. Eckert ERG, Drake RM: Analysis of heat and mass transfer. McGraw-Hill, New York; 1972.

    Google Scholar 

  12. Mortimer RG, Eyring H: Elementary transition state theory of the Soret and Dufour effects. Proc Natl Acad Sci 1980, 77: 1728-1731. 10.1073/pnas.77.4.1728

    Article  Google Scholar 

  13. Anjalidevi SP, Devi RU: Soret and Dufour effects on MHD slip flow with thermal radiation over a porous rotating infnite disk. Commun Nonlinear Sci 2011, 16: 1917-930. 10.1016/j.cnsns.2010.08.020

    Article  Google Scholar 

  14. Kafoussias NG, Williams EW: Thermal-diffusion and diffusion thermo effects on mixed free-forced convective and mass transfer boundary layer flow with temperature dependent viscosity. Int J Eng Sci 1995, 33: 1369-1384. 10.1016/0020-7225(94)00132-4

    Article  MATH  Google Scholar 

  15. 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 2004, 47: 1467-1472. 10.1016/j.ijheatmasstransfer.2003.09.017

    Article  MATH  Google Scholar 

  16. Mansour MA, El-Anssary 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 2008, 145: 340-345. 10.1016/j.cej.2008.08.016

    Article  Google Scholar 

  17. Lakshmi-Narayana 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 2010, 53: 3030-3034. 10.1016/j.ijheatmasstransfer.2010.03.025

    Article  Google Scholar 

  18. Awad FG, Sibanda P, Motsa SS: On the linear stability analysis of a Maxwell fluid with double-diffusive convection. Appl Math Model 2010, 34: 3509-3517. 10.1016/j.apm.2010.02.038

    Article  MATH  MathSciNet  Google Scholar 

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

    Chapter  Google Scholar 

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

    Google Scholar 

  21. Motsa SS, Sibanda P, Shateyi S: A new spectral-homotopy analysis method for solving a nonlinear second order BVP. Commun Nonlinear Sci Numer Simulat 2010, 15: 2293-2302. 10.1016/j.cnsns.2009.09.019

    Article  MATH  MathSciNet  Google Scholar 

  22. Motsa SS, Sibanda P, Awad FG, Shateyi S: A new spectral-homotopy analysis method for the MJD Jeffery-Hamel problem. Comput. Fluids 2010, 39: 1219-1225. 10.1016/j.compfluid.2010.03.004

    MATH  MathSciNet  Google Scholar 

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

    Google Scholar 

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

    Article  MathSciNet  Google Scholar 

  25. Awad FG, Sibanda P, Motsa SS, Makinde OD: Convection from an inverted cone in a porous medium with cross-diffusion effects. Comput Math Appl 2011, 61: 1431-1441. 10.1016/j.camwa.2011.01.015

    Article  MATH  MathSciNet  Google Scholar 

  26. Motsa SS, Sibanda P, Shateyi S: On a new quasi-linearization method for systems of nonlinear boundary value problems. Math Methods Appl Sci 2011, 34: 1406-1413. 10.1002/mma.1449

    Article  MATH  MathSciNet  Google Scholar 

  27. 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 2009, 14: 2202-2214. 10.1016/j.cnsns.2008.07.001

    Article  MATH  Google Scholar 

  28. Sparrow EM, Yu HS: Local non-similarity thermal boundary layer solutions. J Heat Transfer 1971, 93: 328-334. 10.1115/1.3449827

    Article  Google Scholar 

  29. Sparrow EM, Quack H, Boerner CJ: Local non-similarity boundary layer solutions. AIAA 1970, 8: 1936-1942. 10.2514/3.6029

    Article  MATH  Google Scholar 

  30. Shampine LF, Ketzscher R, Forth SA: Using AD to solve BVPs in Matlab. ACM Trans Math Softw 2005, 31(1):79-94. 10.1145/1055531.1055535

    Article  MATH  MathSciNet  Google Scholar 

  31. El-Kabeir 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.

    Google Scholar 

  32. 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 2006, 11: 3-12.

    MATH  MathSciNet  Google Scholar 

Download references

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.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Precious Sibanda.

Additional information

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.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Sibanda, P., Khidir, A.A. & Awad, F.G. On cross-diffusion effects on flow over a vertical surface using a linearization method. Bound Value Probl 2012, 25 (2012). https://doi.org/10.1186/1687-2770-2012-25

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2012-25

Keywords