Important dimensionless parameters of interest.
This paper presents an essential study of scale analysis and double diffusive free convection boundary layer laminar flow of low Prandtl fluids over an inclined wall kept at uniform surface temperature. Buoyancy effect (N) was considered for an assisting flow when N ≥ 0, which implies that the thermal and solutal forces are consolidating each other to help drive the fluid flow in the same direction. Scale analysis and similarity transformation methods are used to obtain the governing equations, and the resulting system of coupled ordinary differential equations (ODEs) is solved with the differential transform method (DTM). Results for the distributions of velocity, temperature, and concentration boundary layer of the fluid adjacent to the wall are presented. The study includes the effects of the ratio of solutal buoyancy to thermal buoyancy and important dimensionless parameters used in this work with varying angles of inclination of the wall on fluid flow and heat transfer.
- scale analysis
- free convection
- double diffusion
- low Prandtl flow
- boundary layer
The impact of temperature and species concentration distribution on heat and mass transfer of fluid flow has received renewed interest to researchers and the academic community due to its multiple application areas notably in physical and chemical processes, food and manufacturing industries, geophysics, oceanography, and photosynthesis. These occurrences of thermo-solutal convection not only involve temperature variation but also concentration variation. Free convection problems driven only by temperature difference have been studied extensively by many investigators notably Akter et al. , Schlichting , Venkateswara , Gebhart and Pera , Khair and Bejan , and Mongruel et al. . Species concentration variation sometimes plays a major role in creating the buoyancy needed in driving flow and influencing rate of heat transfer. Double diffusion has been studied by few investigators over the years. The pioneer of this area of research is the work of Gebhart and Pera in 1971  where they investigated the combined buoyancy effects of thermal and mass diffusion on natural convection flow. Also, Bejan and Khair  carried out some analysis on heat and mass transfer by natural convection in porous media. Furthermore, the Schmidt number is the appropriate number in the concentration equation for
Some other research studies carried out were by Angirasa and Peterson , who considered free convection due to combined buoyancy forces for
However, research conducted to critically analyze fluid behavior with the effect of species concentration and thermal diffusion on heat and mass transfer particularly for low Prandtl flows past an inclined wall is very rare. This gap has been captured in this study. The objective of this research is to investigate the effect of combined heat and species concentration involving a low Prandtl number fluid flow over an inclined wall using the method of scale analysis in formulation of the model along with the similarity transformation technique to convert partial differential equations to ordinary differential equation. The resulting dimensionless coupled and non-linear equations are solved using differential transform method. The numerics of the computation are discussed for different values of dimensionless parameters and are graphically presented.
2. Problem formulation and scale analysis
The problem of combined heat and mass transfer over a heated semi-infinite inclined solid wall is considered. The fluid is assumed to be steady, Newtonian, viscous, and incompressible. It is assumed that the wall is maintained at uniform surface temperature and concentration and it is immersed in fluid reservoir at rest which is kept at uniform ambient temperature and concentration such that > and >. Boundary layer flow over an inclined wall driven by both thermal gradient and concentration gradient, respectively, are thereby set up due to the difference between wall values and quiescent fluid values. Hence, it is called combined heat and mass transfer phenomenon over an inclined wall (Figure 1).
This problem is governed by the non-linear and coupled conservation equations. Using the Boussinesq approximation and boundary layer simplifications, we have the following:
Here, is the velocity along x-axis, is the velocity along y-axis or along the vertical wall,
Following the procedures as outlined by Khair and Bejan , it can be clearly shown that for low Prandtl number flows, the velocity, concentration, and temperature boundary layer scales (for N = 0) are:
While vertical velocity (
The similarity variable for velocity layer scales as
And corresponding stream function scales as
where is the velocity function.
Boussinesq approximation is used in Eq. (2), and the PDEs are reduced to a set of coupled ODEs using similarity variable. It can easily be shown that for the inner layer of low Prandtl number flows, the dimensionless momentum, energy and concentration equations give:
the similarity variable for concentration layer is defined as , where is the thin viscous layer closest to the wall and is the concentration boundary layer. The attendant stream function is obtained as
The ordinary differential equations governing the momentum, energy, and species concentration become:
In the outer layer, where there is inertia-buoyancy balance, which is the similarity variable for thermal layer is defined as and the associated stream function obtained as .
The resulting dimensionless equations for low Prandtl number flow are:
In further works, these equations will be solved asymptotically as
2.1 Method of solution
Such that the governing equations of motion become:
Due to limitation of convergence of the classical DTM which is only valid near the multi-step transformation is used. Carrying out multi-step differential transformations, we have:
The transformation of the associated boundary conditions follows as:
, , , , ,
|The local Nusselt number: |
The Local Sherwood number:
The shearing stress on the plate:
The coefficient of skin friction:
3. Results and discussion
Table 1 shows the expressions for the dimensionless parameters that are of interest in this work. The solutions to Eqs. (13)–(15) subjected to Eq. (16) solved using the multi-step DTM method are presented graphically in the figures below. The results shown are for Prandtl numbers of 0.01, 0.1, 0.5, and 0.72, respectively.
Figure 2 shows the profiles of local skin friction against Prandtl number for various angles of inclination at a constant Lewis number. It could be observed from the plots that the shearing stress decreases as the Prandtl number increases for all the plate angles considered. More importantly, it is illustrated by the graphs that the local skin friction also decreases with increase in the buoyancy ratio and also with respect to the angle of inclination of the wall.
Figure 3 shows the plots of local Nusselt number against Prandtl number for various angles of inclination at a constant Lewis number. It is clearly seen from the results that the rate of heat transfer (Nusselt number) increases as the Pr increases for all the wall inclination angles. Also to note is the fact that Nusselt number increases as the buoyancy ratio is increased.
Figure 4 shows the results of how the local Sherwood number changes as the Lewis number is increased for various angles of inclination at a constant Prandtl number. It is noted from the graphs that the local Sherwood number increases as the Lewis number increases for all angles of inclination. Also, it is observed from the figures that the rate of increase of Sherwood number is dependent on N as well as the angle of inclination of the wall.
In Figure 5, the similarity profiles of the effect of the buoyancy ratio and the angles of inclination on the dimensionless velocities at fixed Prandtl and Lewis numbers is presented. It could be interpreted from the results for N = 0 that a maximum velocity is obtained for a vertical wall while at an angle of 60°, the minimum velocity. Also, it is clearly seen from the plots that aside from the vertical wall possessing the maximum velocity, its vertical velocity value for N = 1 is higher when compared to the case when N = 0. When the buoyancy ratio N is further increased to 1.5, the velocity for the vertical wall also increases. The trend in the figure also shows that increasing the inclination angle increases the velocity boundary layer thickness for all buoyancy ratio.
In Figure 6, the similarity profiles of dimensionless temperature for the wall angles of inclinations considered in this study at constant Lewis number are presented. The temperature profiles show that despite the varying buoyancy ratio, the thermal boundary layer thickness increases as the inclination angle of the wall increases.
Figure 7 presents the similarity profiles of dimensionless species concentration distributions for the angles of inclination of the wall at a constant Lewis number. It is clearly seen that as the buoyancy ratio is increased, the concentration boundary layer thickness has a decreasing trend under this same flow configuration but as the angle of inclination increases, increases.
Figure 8 shows the plots of the coupled similarity profiles of dimensionless temperature, species concentration, and velocity for wall inclination angle of 60° under the constant Lewis number of 10. It can be seen from the results that increasing N has negligible effect on the trio of velocity, concentration and temperature boundary layer thicknesses for a fixed wall angle. However, the effect of N is very noticeable in the vertical velocity values which is clearly higher when N = 1.5 compared to when N = 0.
The problem of double diffusive convection and its associated boundary layer flow is of tremendous interest in academic research and various manufacturing and process industries because of its implications in energy and mass transfer efficiency in engineering and scientific applications. In this study, scale analysis and double diffusive free convection of low Prandtl fluid flow over an inclined wall is investigated in the presence of species concentration and thermal diffusion. The governing boundary layer equations obtained by scale analysis are numerically solved using differential transform method (DTM). The conclusions reached as a result of the parametric study conducted are presented below:
The velocity boundary layer thickness is maximum when the plate is in the vertical position, while it is minimum when the plate is at the inclined angle of 60° to the vertical for any value of Lewis number but within a certain buoyancy ratio.
The thermal boundary thickness is maximum when the wall is inclined at 60° to the vertical and minimum when in a vertical position while keeping N constant.
Thermal boundary thickness decreases with increase in Prandtl number for all angles of inclination while keeping both the Lewis number and buoyancy ratio constant.
The concentration decreases with the increase in Lewis number for all range of values considered for N and all angles of inclination which is in agreement with Akter et al. . Also, as N increases, both concentration and temperature increase while velocity decreases with increase in plate angles.
The velocity, concentration, and thermal boundary layer thicknesses increase with an increase in the angle of inclination of the wall.
chemical species concentration
chemical species diffusivity
thermal Rayleigh number
velocity component in x-direction
velocity component in y-direction
coefficient of thermal expansion
coefficient of specie expansion
thermal diffusivity of fluid
wall derivative of dimensionless concentration
velocity boundary layer thickness
thermal boundary layer thickness
concentration boundary layer thickness
wall derivative of dimensionless temperature
density of fluid
constant wall dimensionless heat flux
condition at infinity
condition at the wall
Akter F, Islam MM, Islam A, Khan MS, Hossain MS. Chemical reaction and thermal diffusion effects on mass transfer flow through an inclined plate. Open Journal of Fluid Dynamics. 2016; 6:62. DOI: 10.4236/ojfd.2016.61006
Schlichting H. Boundary-Layer Theory. 7th ed. New York: McGraw-Hill; 1968
Rajua KV, Reddyb PBA, Suneethac S. Thermophoresis effect on a radiating inclined permeable moving plate in the presence of chemical reaction and heat absorption. International Journal of Dynamics of Fluids. 2017; 13:89-112
Gebhart B, Pera L. The nature of vertical natural convection flows resulting from the combined buoyancy effects of thermal and mass diffusion. International Journal of Heat and Mass Transfer. 1971; 14:2025-2050. DOI: 10.1016/0017-9310(71)90026-3
Khair KR, Bejan A. Mass transfer to natural convection boundary layer flow driven by heat transfer. Journal of Heat Transfer. 1985; 107:979-981. DOI: 10.1115/1.3247535
Mongruel A, Cloitre M, Allain C. Scaling of boundary-layer flows driven by double-diffusive convection. International Journal of Heat and Mass Transfer. 1996; 39:3899-3910. DOI: 10.1016/0017-9310(96)00054-3
Bejan A, Khair KR. Heat and mass transfer by natural convection in a porous medium. International Journal of Heat and Mass Transfer. 1985; 28:909-918. DOI: 10.1016/0017-9310(85)90272-8
Allain C, Cloitre M, Mongruel A. Scaling in flows driven by heat and mass convection in a porous medium. EPL (Europhysics Letters). 1992; 20:313. DOI: 10.1209/0295-5075/20/4/005
Angirasa D, Srinivasan J. Natural convection flows due to the combined buoyancy of heat and mass diffusion in a thermally stratified medium. Journal of Heat Transfer. 1989; 111:657-663. DOI: 10.1115/1.3250733
Angirasa D, Peterson GP, Pop I. Combined heat and mass transfer by natural convection with opposing buoyancy effects in a fluid saturated porous medium. International Journal of Heat and Mass Transfer. 1997; 40:2755-2773. DOI: 10.1016/S0017-9310(96)00354-7
Murthy SK, Kumar BR, Chandra P, Sangwan V, Nigam M. A study of double diffusive free convection from a corrugated vertical surface in a Darcy porous medium under Soret and Dufour effects. Journal of Heat Transfer. 2011; 133:092601. DOI: 10.1115/1.4003813
Ferdows M, Khan MS, Bég OA, Azad MAK, Alam MM. Numerical study of transient magnetohydrodynamic radiative free convection nanofluid flow from a stretching permeable surface. Proceedings of the Institution of Mechanical Engineers, Part E: Journal of Process Mechanical Engineering. 2014; 228:181-196. DOI: 10.1177/0954408913493406
Khan MS, Karim I, Ali LE, Islam A. Unsteady MHD free convection boundary-layer flow of a nanofluid along a stretching sheet with thermal radiation and viscous dissipation effects. International Nano Letters. 2012; 2:24. DOI: 10.1186/2228-5326-2-24
Umemura A, Law CK. Natural-convection boundary-layer flow over a heated plate with arbitrary inclination. Journal of Fluid Mechanics. 1990; 219:571-584. DOI: 10.1017/S0022112090003081
Chamka A, Khaled ARA. Simultaneously heat and mass transfer in free convection. Industrial Engineering Chemical. 2001; 49:961-968. DOI: 10.1021/ie50570a025
Reddy MG, Reddy NB. Mass transfer and heat generation effects on MHD free convection flow past an inclined vertical surface in a porous medium. Journal of Applied Fluid Mechanics. 2011; 4:7-11
Singh PK. Heat and mass transfer in MHD boundary layer flow past an inclined plate with viscous dissipation in porous medium. International Journal of Scientific and Engineering Research. 2012; 3:1-11
Ali LE, Islam A, Islam N. Investigate micropolar fluid behavior on MHD free convection and mass transfer flow with constant heat and mass fluxes by finite difference method. American Journal of Applied Mathematics. 2015; 3:157-168. DOI: 10.11648/j.ajam.20150303.23
Islam M, Akter F, Islam A. Mass transfer flow through an inclined plate with porous medium. American Journal of Applied Mathematics. 2015; 3:215-220. DOI: 10.11648/j.ajam.20150305.12