Temperature and Concentration Stratification Effects in Mixed Convection Flow of an Oldroyd-B Fluid with Thermal Radiation and Chemical Reaction.

This research addresses the mixed convection flow of an Oldroyd-B fluid in a doubly stratified surface. Both temperature and concentration stratification effects are considered. Thermal radiation and chemical reaction effects are accounted. The governing nonlinear boundary layer equations are converted to coupled nonlinear ordinary differential equations using appropriate transformations. Resulting nonlinear systems are solved for the convergent series solutions. Graphs are plotted to examine the impacts of physical parameters on the non-dimensional temperature and concentration distributions. The local Nusselt number and the local Sherwood number are computed and analyzed numerically.

Introduction

Analysis of non-Newtonian fluids has great importance due to its several industrial and engineering applications. In particular these fluids are encountered in the material processing, chemical and nuclear industries, bioengineering, oil reservoir engineering, polymeric liquids and foodstuffs. Several fluids like paints, paper pulp, shampoos, ketchup, apple sauce, slurries, certain oils and polymer solutions are the non-Newtonian fluids. The characteristics of all the non-Newtonian fluids cannot be explained via one constitutive relationship. Hence various fluid models are proposed in the literature for the properties of non-Newtonian fluids. Generally non-Newtonian materials are classified under three categories namely (i) differential type (ii) rate type and (iii) integral type. The Maxwell fluid model is the simplest subclass of rate type fluids. This model describes only the properties of relaxation time. The characteristics of retardation time cannot be predicted by the Maxwell fluid. An Oldroyd-B fluid model was developed to examine both the relaxation and retardation times characteristics. Instabilities in viscoelastic liquids were studied by Larson [1]. In this investigation, he discussed the instabilities in Taylor-Couette flows, instabilities in cone-and-plate and plate-and-plate flows, instabilities in parallel shear flows, instabilities in external and multi-dimensional flows. The instabilities in the flows occurring in the absence of inertial forces were investigated by Shaqfeh [2]. Laso and Ottinger [3] presented a study to examine the numerical simulation of viscoelastic liquids based on molecular models. Rajagopal and Bhatnagar [4] computed the asymptotically decaying solution of an Oldroyd-B fluid past an infinite porous plate. Thermodynamic properties of rate type non-Newtonian fluids were investigated by Rajagopal and Srinivasa [5]. Numerical solutions of Oldroyd-B and PTT-fluids with both the linear and exponential stress functions were developed by Alves et al. [6]. Some recent investigations on non-Newtonian fluids can be seen in the references [7-14].

Heat and mass transfer analysis in the boundary layer flow over a stretching surface has key role in the industrial and engineering applications, for example, manufacturing of plastic and rubber sheets, annealing and thinning of copper wires, drawing on stretching sheets through quiescent fluids, boundary layer along a liquid film condensation process, damage of crops due to freezing, desalination, refrigeration and air conditioning, compact heat exchangers, solar power collectors, human transpiration and many others (see refs. [15,16]). Heat and mass transfer effects in boundary layer flow of viscoelastic fluid with thermal slip condition were investigated by Turkyilmazoglu [17]. Hayat and Alsaedi [18] carried out a study to examine the heat and mass transfer phenomena in buoyancy driven flow of an Oldroyd-B fluid. Thermophoresis and Joule heating effects are further considered. Hayat et al. [19] presented the series solutions of magnetohydrodynamic (MHD) flow of Casson fluid with heat and mass transfer. Soret and Dufour effects are present in this investigation. Mixed convection flow of Jeffrey fluid in the presence of heat and mass transfer is investigated by Shehzad et al. [20]. Gupta et al. [21] discussed the effect of cadmium on growth and active constituents of bacopa monnieri. Induced magnetic field effect in mixed convection peristaltic flow of third order fluid with nanoparticles is discussed by Noreen [22]. Bachok et al. [23] studied the boundary layer flow of viscous fluid in presence of mixed convection and viscous dissipation. Su et al. [24] developed a lattice Boltzmann method coupled with the Oldroyd-B constitutive equation to stimulate flow of viscoelastic fluid. Here the numerical results of 2D channel flow agree well with the analytical and some experimental results reported in the previous studies. Slip effects in peristaltic flow of generalised Oldroyd-B fluids is explored by Tripathi et al. [25]. They computed the homotopic solutions of the modelled differential system.

Effect of stratification is an important aspect in heat and mass transfer analyses. Stratification of fluids occurs due to temperature variations, concentration differences or the presence of different fluids of different densities. When the heat and mass transfer are present simultaneously then it is important to analyze the effect of double stratification on the convective flows. The analysis of mixed convection in a doubly stratified medium is an important problem. It is because of its occurrence in geophysical flows (see ref. [26]). Such flows involve in the rivers, lakes and seas, thermal energy storage systems and solar ponds etc. Chang and Lee [27] investigated the free convection flow by a vertical plate with uniform and constant heat flux in a thermally stratified micropolar fluid. Cheng [28] examined the combined heat and mass transfer effect in natural convection flow from a vertical wavy surface in a power-law fluid saturated porous medium. Both thermal and mass stratification effects were present. Srinivasacharya and Reddy [29] discussed the effect of double stratification in mixed convection flow of micropolar fluid. Effect of double stratification on MHD free convection flow of micropolar fluid is investigated by Srinivasacharya and Upendar [30]. Non-Darcy mixed convection flow in a doubly stratified medium under Soret and Dufour effects is studied by Srinivasacharya and Surender [31]. Srinivasacharya and Surender [32] addressed the effect of double stratification on mixed convection boundary layer flow of a nanofluid past a vertical plate in porous medium.

The basic theme of present study is to investigate the effects of thermal radiation, chemical reaction, thermal and solutal stratification in the mixed convection boundary layer flow of an Oldroyd-B fluid over a stretching surface. The studies available in the literature on this topic mostly dealt with the thermal stratification effect. Some recent aforementioned studies investigated the effects of both thermal and concentration stratification in viscous fluid flow. This is the first attempt to study such effects for non-Newtonian fluids. Mathematical modelling is developed under the consideration of thermal and concentration stratification effects. The series solutions to the resulting nonlinear differential systems are constructed via homotopy analysis method (HAM) [33-41]. The effects of various emerging parameters on the temperature and concentration fields are presented through plots and tables. The local Nusselt and the local Sherwood numbers are computed numerically and analyzed.

Mathematical Modeling

We consider the steady two-dimensional doubly stratified mixed convection flow of an incompressible Oldroyd-B fluid. The flow is caused by a linearly stretching surface at y = 0. The flow occupies the domain y > 0. Boundary layer flow is considered in the presence of thermal radiation and first order chemical reaction. The governing boundary layer equations for incompressible flow of an Oldroyd-B fluid with heat and mass transfer are given below (see Appendix for detailed derivation):

The appropriate boundary conditions are where u and v are the velocity components in the x– and y–directions respectively, λ 1 the relaxation time, ν = μ / ρ the kinematic viscosity, μ the dynamic viscosity, ρ the density of fluid, λ 2 the retardation time, g the gravitational acceleration, β the thermal expansion coefficient, T the temperature, β the concentration expansion coefficient, C the concentration, α = k / ρc the thermal diffusivity of the fluid, k the thermal conductivity, c the specific heat at constant pressure, q the radiative heat flux, D the diffusion coefficient, K 1 the reaction rate, T and T the temperatures of the surface and far away from the surface and C and C the concentrations at the surface and far away from the surface. The subscript w denotes the wall condition. This study assumes that the surface stretching velocity, wall temperature and wall concentration are where a, A 1, B 1, M 1, N 1, T ∞,0 and C ∞,0 are the positive constants. The radiative heat flux q via Rosseland's approximation can be expressed as follows: in which σ 1 is the Stefan-Boltzman constant and m is the mean absorption coefficient. We assume that the difference in temperature within the flow is such that T 4 can be written as a linear combination of temperature. By employing Taylor's series and neglecting higher order terms we have [18]:

Substituting Eq (10) in Eq (9) we get

Using Eq (11) in Eq (4) we have

The dimensionless variables can be defined as follows:

Incompressibility condition is now identically satisfied and Eqs (2)–(8) and (12) become

In above expressions β 1 and β 2 are the Deborah numbers in terms of relaxation and retardation times respectively, λ is the mixed convection parameter, Gr is the Grashof number, Re is the local Reynolds number, N is the buoyancy ratio parameter, Rd is the thermal radiation parameter, Pr is the Prandtl number, ε 1 is the thermal stratification parameter, Sc is the Schmidt number, γ is the chemical reaction parameter, ε 2 is the solutal stratification parameter and prime stands for differentiation with respect to η. Note that when β 2 = 0, this analysis reduced to the Maxwell fluid flow case. Eq (15) indicates that P is constant in the y–direction. The involved variables can be defined as follows:

The local Nusselt number Nu and the local Sherwood number Sh are given by

Series Solutions

The initial guesses and the linear operators are

The operators satisfy the following properties [33]: in which C (i = 1–7) are the arbitrary constants.

We can define the following zeroth-order deformation problems [36,40]:

Here p denotes the embedding parameter and the non-zero auxiliary parameters and N , N and N the nonlinear operators. Setting p = 0 and p = 1 we have [35,38]:

When p varies from 0 to 1 then and vary from the initial guesses f 0(η), θ 0(η) and ϕ 0(η) to the final solutions f(η), θ(η) and ϕ(η), respectively. Taylor series expansion gives [37,39]:

The convergence of above series strongly depends upon and Considering that and are chosen in such a manner that Eqs (37)–(39) converge at p = 1 then [33,34]:

The mth-order problems are [33]:

The mth-order deformation problems have the solutions in which and denote the special solutions.

Convergence Analysis

The series solutions (40)–(42) involve the auxiliary parameters and These parameters are useful in adjusting and controlling the convergence of the obtained series solutions. The proper values of these parameters are quite essential to construct the convergent solutions via homotopy analysis method. To choose the suitable values of and the curves are plotted at 13th order of approximations. Fig 1 clearly depicts that the convergence region lies within the domain and Furthermore the series solutions converge in the whole region of η when Table 1 shows that the 11th order of approximations are sufficient for the convergent series solutions.

Fig 1

The curves for the functions f(η), θ(η) and ϕ(η) when β 1 = β 2 = 0.2 = γ = Rd, λ = 0.1 = N, Pr = 1.0 = Sc and ε 1 = 0.3 = ε 2.

Table 1

Convergence of HAM solutions for different order of approximations when β 1 = β 2 = 0.2 = γ = Rd, λ = 0.1 = N, Pr = 1.0 = Sc, ε 1 = 0.3 = ε 2 and

Order of approximations	−f′′ (0)	−θ′ (0)	−ϕ′ (0)
1	0.91150	1.01333	1.15333
5	0.93405	0.98685	1.15531
11	0.93419	0.98680	1.15505
20	0.93419	0.98680	1.15505
35	0.93419	0.98680	1.15505
50	0.93419	0.98680	1.15505

Results and Discussion

This section presents the impacts of various emerging parameters including Deborah number in terms of relaxation time β 1, Deborah number in terms of retardation time β 2, Prandtl number Pr, thermal radiation parameter Rd, Schmidt number Sc, chemical reaction parameter γ, thermal stratification parameter ε 1 and solutal stratification parameter ε 2 on the dimensionless temperature profile θ(η) and concentration profile ϕ(η). This purpose is achieved through the plots in the Figs 2–9. Fig 2 is plotted to examine the effects of Deborah number β 1 on the temperature profile θ(η) and concentration profile ϕ(η) when β 1 = 0.0, 0.25, 0.50 and β 2 = 0.2 = γ = Rd, λ = 0.1 = N, Pr = 1.0 = Sc, ε 1 = 0.3 = ε 2. Fig 2 examined that the temperature profile θ(η) and concentration profile ϕ(η) are enhanced when we use the larger values of Deborah number β 1. Since Deborah number β 1 has dependence on the relaxation time [4]. Larger values of Deborah number β 1 implies to higher relaxation time. It is well known fact that the larger relaxation time fluids have higher temperature and concentration and smaller relaxation time fluids possess lower temperature and concentration. In view of this argument, both temperature profile θ(η) and concentration profile ϕ(η) are enhanced via larger Deborah number β 1. The influence of Deborah number β 2 on the dimensionless temperature and concentration fields when β 2 = 0.0, 0.25, 0.50 and β 1 = 0.2 = γ = Rd, λ = 0.1 = N, Pr = 1.0 = Sc, ε 1 = 0.3 = ε 2 is studied in Fig 3. Fig 3 clearly depicts that the temperature θ(η) and concentration ϕ(η) are decreasing functions of Deborah number β 2 [11]. Here the Deborah number β 2 is dependent on the retardation time. When we increase the values of Deborah number β 2, the retardation time is increased. Such increase in retardation time is responsible for the reduction in the temperature θ(η) and concentration ϕ(η). Here it is interesting to mention that β 1 = 0 = β 2 correspond to viscous fluid case and β 2 = 0 shows the Maxwellian fluid flow situation. From experimental point of view, it is quite obvious that the values of β 2 are not much than the values of β 1. Influence of thermal stratification parameter ε 1 on the temperature θ(η) and concentration ϕ(η) is shown in Fig 4 when ε 1 = 0.0, 0.1, 0.2 and β 1 = β 2 = 0.2 = γ = Rd, λ = 0.1 = N, Pr = 1.0 = Sc, ε 2 = 0.3. Here the temperature and thermal boundary layer thickness are decreased while concentration and its related boundary layer thickness are increased when we increase in thermal stratification parameter. When the thermal stratification effect is taken into account, the effective temperature difference between the surface and the ambient fluid is decreased while opposite behavior is observed for concentration profile [28]. Influence of solutal stratification parameter ε 2 on the temperature profile θ(η) and concentration profile ϕ(η) is shown in Fig 5 when ε 2 = 0.0, 0.1, 0.2 and β 1 = β 2 = 0.2 = γ = Rd, λ = 0.1 = N, Pr = 1.0 = Sc and ε 1 = 0.3. The temperature profile is enhanced while the concentration profile is reduced with an increase in solutal stratification parameter [29]. Influence of Prandtl number on the temperature profile is shown in Fig 6 when Pr = 0.5, 0.75, 1.0, 1.25 and β 1 = β 2 = 0.2 = γ = Rd, λ = 0.1 = N, Sc = 1.0, ε 1 = 0.3 = ε 2. The temperature and thermal layer thickness are reduced for the increasing values of Prandtl number. Physically larger Prandtl fluids possess lower thermal diffusivity and smaller Prandtl fluids have higher thermal diffusivity. This change in thermal diffusivity causes a reduction in the temperature and thermal boundary layer thickness. Basically Prandtl number is the ratio of momentum diffusivity to thermal diffusivity. In heat transfer, Prandtl number is used to control the thicknesses of momentum and thermal boundary layers. Fig 7 is plotted to examine the change in temperature profile when Rd = 0.0, 0.3, 0.6, 1.0 and β 1 = β 2 = 0.2 = γ, λ = 0.1 = N, Pr = 1.0 = Sc, ε 1 = 0.3 = ε 2. Fig 7 describes that the temperature and thermal boundary layer thickness are enhanced with an increase in the thermal radiation parameter. Larger values of thermal radiation parameter provide more heat to working fluid that shows an enhancement in the temperature and thermal boundary layer thickness [20]. Influence of Schmidt number on the concentration field is shown in Fig 8 when Sc = 0.5, 0.75, 1.0, 1.25 and β 1 = β 2 = 0.2 = γ = Rd, λ = 0.1 = N, Pr = 1.0, ε 1 = 0.3 = ε 2. It is clearly observed that the concentration and its related boundary layer thickness are decreasing functions of Schmidt number. Schmidt number is inversely proportional to the diffusion coefficient. Hence an increase in Schmidt number corresponds to a smaller diffusion coefficient. Such smaller diffusion coefficient creates a reduction in the concentration field. Fig 9 is plotted to investigate the effects of chemical reaction parameter when γ = 0.0, 0.3, 0.6, 1.0 and β 1 = β 2 = 0.2 = Rd, λ = 0.1 = N, Pr = 1.0 = Sc, ε 1 = 0.3 = ε 2. It is noticed from Fig 9 that the concentration and its associated boundary layer thickness are decreasing functions of chemical reaction parameter. Chemical reaction increases the rate of interfacial mass transfer. The reaction reduces the local concentration, thus increasing the concentration gradient and its flux. As a result, concentration of the chemical species in the boundary layer decreases with an increase in chemical reaction parameter. Tables 2 and 3 show the numerical values of the local Nusselt and the local Sherwood numbers for different values of β 1, β 2, λ, N, Pr, Sc, Rd, γ, ε 1 and ε 2. The values of local Nusselt and the local Sherwood numbers are decreased by increasing ε 1 and ε 2 while these values are increased for the larger λ and N. Table 4 is computed to validate the present results with the previous published results in a limiting sense. Here we compared our results for a Maxwell fluid case. From this Table, we examined that the present series solutions have good agreement with the numerical solutions of Megahed [42] in limiting sense.

Fig 2

Temperature distribution function θ(η) and concentration distribution function ϕ(η) when β 1 = 0.0, 0.25, 0.50 and β 2 = 0.2 = γ = Rd, λ = 0.1 = N, Pr = 1.0 = Sc, ε 1 = 0.3 = ε 2.

Fig 9

Concentration distribution function ϕ(η) when γ = 0.0, 0.3, 0.6, 1.0 and β 1 = β 2 = 0.2 = Rd, λ = 0.1 = N, Pr = 1.0 = Sc, ε 1 = 0.3 = ε 2.

Fig 3

Temperature distribution function θ(η) and concentration distribution function ϕ(η) when β 2 = 0.0, 0.25, 0.50 and β 1 = 0.2 = γ = Rd, λ = 0.1 = N, Pr = 1.0 = Sc, ε 1 = 0.3 = ε 2.

Fig 4

Temperature distribution function θ(η) and concentration distribution function ϕ(η) when ε 1 = 0.0, 0.1, 0.2 and β 1 = β 2 = 0.2 = γ = Rd, λ = 0.1 = N, Pr = 1.0 = Sc, ε 2 = 0.3.

Fig 5

Temperature distribution function θ(η) and concentration distribution function ϕ(η) when ε 2 = 0.0, 0.1, 0.2 and β 1 = β 2 = 0.2 = γ = Rd, λ = 0.1 = N, Pr = 1.0 = Sc and ε 1 = 0.3.

Fig 6

Temperature distribution function θ(η) when Pr = 0.5, 0.75, 1.0, 1.25 and β 1 = β 2 = 0.2 = γ = Rd, λ = 0.1 = N, Sc = 1.0, ε 1 = 0.3 = ε 2.

Fig 7

Temperature distribution function θ(η) when Rd = 0.0, 0.3, 0.6, 1.0 and β 1 = β 2 = 0.2 = γ, λ = 0.1 = N, Pr = 1.0 = Sc, ε 1 = 0.3 = ε 2.

Fig 8

Concentration distribution function ϕ(η) when Sc = 0.5, 0.75, 1.0, 1.25 and β 1 = β 2 = 0.2 = γ = Rd, λ = 0.1 = N, Pr = 1.0, ε 1 = 0.3 = ε 2.

Table 2

Values of local Nusselt number (Re)−1/2 Nu and local Sherwood number (Re)−1/2 Sh for different values of β 1, β 2, λ and N when Pr = 1.0 = Sc, Rd = 0.2 = γ and ε 1 = 0.3 = ε 2.

β 1	β 2	λ	N	−(1 + Rd)θ′ (0)	−ϕ′ (0)
0.0	0.2	0.1	0.1	0.97961	0.97656
0.2				0.95671	0.96073
0.4				0.93563	0.94621
0.2	0.0	0.1	0.1	0.92966	0.94168
	0.2			0.95671	0.96073
	0.4			0.97899	0.97648
0.2	0.2	0.0	0.1	0.95032	0.95590
		0.2		0.96267	0.96528
		0.5		0.97862	0.97761
0.2	0.2	0.1	0.0	0.95636	0.96045
			0.2	0.95706	0.96101
			0.5	0.95809	0.96186

Table 3

Values of local Nusselt number (Re)−1/2 Nu and local Sherwood number (Re)−1/2 Sh for different values of Pr, Sc, Rd, γ, ε 1 and ε 2 when β 1 = 0.2 = β 2 and λ = 0.1 = N.

Pr	Sc	Rd	γ	ε 1	ε 2	−(1 + Rd)θ′ (0)	−ϕ′ (0)
0.5	1.0	0.2	0.2	0.3	0.3	0.60210	0.96913
1.0						0.95671	0.96073
1.5						1.24219	0.95807
1.0	0.5	0.2	0.2	0.3	0.3	0.95754	0.61970
	1.0					0.95671	0.96073
	1.5					0.95647	1.22937
1.0	1.0	0.0	0.2	0.3	0.3	0.89755	0.95936
		0.3				0.98338	0.96143
		0.5				1.03215	0.96284
1.0	1.0	0.2	0.0	0.3	0.3	0.95674	0.89936
			0.2			0.95671	0.96073
			0.5			0.95667	1.04369
1.0	1.0	0.2	0.2	0.0	0.3	1.08817	0.96706
				0.5		0.86751	0.95648
				1.0		0.63912	0.94576
1.0	1.0	0.2	0.2	0.3	0.0	0.95736	1.12364
					0.5	0.95628	0.85206
					1.0	0.95519	0.58011

Table 4

Comparative values of –f″(0) for various values of β 1 when β 2 = λ = N = 0.

β 1	Megahed [42]	Present work
0.0	0.999978	1.0000
0.2	1.051945	1.0519
0.4	1.101848	1.1019
0.6	1.150160	1.1501
0.8	1.196690	1.1967
1.2	1.285253	1.2853
1.6	1.368641	1.3686
2.0	1.447616	1.4476

Conclusions

Influence of double stratification in mixed convection flow of an Oldroyd-B fluid with thermal radiation and chemical reaction are examined. This analysis reduces to the Maxwell fluid flow case when β 2 = 0. The main findings of this research are given below.

Temperature and concentration fields are increased when we increase the values of β 1. Here the relaxation time is enhanced when we give rise to Deborah number that leads to the higher temperature and concentration fields. This observation is similar that obtained in [18].

Temperature and thermal boundary layer thickness are decreased when Prandtl number increases. Prandtl number is used to control the heat transfer rate in industrial process [29]. The proper value of Prandtl number is quite essential to control the heat transfer rate in industrial and engineering processes.

Temperature and thermal boundary layer thickness are increasing functions of thermal radiation parameter. An increase in thermal radiation parameter provides more heat to fluid due to which larger temperature and thicker thermal boundary layer thickness are achieved.

Influence of thermal stratification parameter on the temperature and concentration fields are opposite [30]. Here temperature is decreased while concentration is enhanced for the higher values of thermal stratification parameter.

Temperature and concentration fields show opposite behavior for solutal stratification parameter. It is also observed that the thermal and concentration boundary layer thicknesses are reverse for the larger solutal stratification parameter.

Concentration field and associated boundary layer thickness are reduced when we increase the values of chemical reaction parameter. For γ = 0 our analysis reduces to the case when there is no chemical reaction.

Table 4 shows that our solutions have an excellent agreement with the previous published numerical results in limiting sense.

The used technique for the solutions development and analysis has advantages over the other in the sense of following points:

It is independent of small/large physical parameters.

It provides a simple way to ensure the convergence of series solutions.

It provides a large freedom to choose the base functions and related auxiliary linear operators.

Besides this the presented analysis is capable of describing relaxation and retardation times feature which many polymers show. Such analysis is particularly useful in polymer extrusion coating process, blood related viscoelastic effects in hemodynamic [43-45] etc. Such analysis provides a stimulus for future investigations on the topic in regimes of magnetohydrodynamics and convective conditions of heat transfer at the surface. It should be pointed out that presented analysis is not able to describe the rheological fluid properties in terms of normal stress effects, shear thinning and shear thickening features. Further the present analysis just deals with the hydrodynamic flow situation due to which the effects of viscous dissipation and Joule heating are ignored. In future one can discuss the hydromagnetic flow case in the presence of Joule heating and viscous dissipation.

Appendix

Here we include the derivation of the governing equation. The constitutive equation for an Oldroyd-B fluid is expressed as follows: where τ represents the components of Cauchy stress tensor, p the pressure, δ the components of identity tensor and the components of an extra stress tensor S is defined as follows: where μ denotes the dynamic viscosity, λ 1 and λ 2 are the relaxation and retardation times respectively, the components of first Rivlin-Ericksen tensor and the contravariant convective derivative which can be written below in the forms

The above expressions represent a contravariant vector and a contravariant tensor having rank 1 and 2 respectively, (where v denote the components of velocity and “,” represents the covariant derivative). In case of Cartesian coordinates the covariant derivative reduces to the usual partial derivative. The equations governing the flow are where ρ denotes the fluid density, g the gravitational field and the definition of acceleration ‘a’ is

By applying the operator Eq (A6) gives or

Assuming to derive the above equation that and Eq (A2), we get the following equation:

For an incompressible steady and two-dimensional flow one can write

For x-component of momentum equation we have

Using Eqs (A12–A18) and Eq (A11) becomes

Similarly y-component of momentum equation can be written as follows:

Here u and v show the velocities parallel to the x- and y-axes respectively and ν the kinematic viscosity.

Using the boundary layer approximations, i.e.,

Eqs (A20) and (A21) reduce to the form

Using boundary condition (6), Eq (A23) becomes

Putting Eq (A25) in Eq (A23) we get

Taylor's series expansion about ρ gives

Neglecting square and higher terms of (T – T ) and (C – C ) we have

This is Boussinesq approximation with

Substituting Eq (A28) into Eq (A26) we get