Abstract

The present analysis is concerned with the criteria for the onset of flow reversal of the fully developed mixed convection in a vertical channel under the effect of the chemical reaction. The governing equations and the critical values of the buoyancy force are solved and calculated numerically via MAPLE. Parameter zones for the occurrence of reversed flow are presented. The exothermic chemical reaction is found to enhance the flow reversal and made flow reversal possible for symmetrical walls temperature.

1. Introduction

The study of combined and free convection flow in a heated vertical parallel-plate channel has received considerable attention because of its wide range of applications such as geothermal reservoir, cooling of nuclear reactors, thermal insulation, energy storage, and conservation, chemical, food, and metallurgical industries, and petroleum reservoirs. A combined forced and natural convective flow puts on the situation where an externally driven flow joins with an native buoyant flow and when both are prominent. One of the earliest studies on mixed convection in a vertical channel with uniform wall temperatures had been studied analytically by Tao [1]. Then Habchi and Acharya [2] studied mixed convection in a vertical channel with asymmetric walls heating, where one plate is heated and the other is adiabatic. Fully developed flow was revisited by Aung and Worku [3] and Cheng et al. [4]. They analyzed effects of various boundary conditions and different constant temperature parameters on buoyancy-aided flow. Recently, Pop et al. [5] studied the mixed convection in a vertical channel with chemical reaction and found that dual solutions exist for both velocity and temperature.

The high buoyancy force due to a differentially heated wall combined with an upward flow leads to a high fluid flow adjacent to the walls that can precipitate a downward flow (i.e., reversal flow) emanating from the open top of the channel in order to augment the increased upward flow. Experiment on the appearance of flow reversal was conducted by Sparrow et al. [6] for free convection in a one-sided heated vertical channel. Aung and Worku [3] and Cheng et al. [4] studied analytically the occurrence flow reversal for the mixed convection. Criteria for the onset of flow reversal in terms of the Grashof to Reynolds number ratio were developed. Criteria for the occurrence of this flow, adjacent to the colder wall under effects of concentration, and magnetic field including internal heating, micropolar as well as thermophoretic, were presented by [79] and [10, 11], respectively.

Numerical studies by Pop et al. [5] found that reversed flow also occurs near the hot wall, where physically it is not possible [3]. They demonstrated that increasing the internal heating and constant temperature parameters could enhance the reversed flow but overlooked determining the criteria for the onset of flow reversal. Therefore, in this paper, we shall present the conditions for the onset of flow reversal of the fully developed flow mixed convection flow in a vertical channel under effect of the chemical reaction for various values of the constant temperature parameters studied by Pop et al. [5].

2. Mathematical Formulation

Consider the steady flow of a viscous and incompressible fluid between two vertical and parallel plane walls. The distance between the walls, that is, the channel width, is . A coordinate system is chosen such that the -axis is parallel to the gravitational acceleration vector g, but with the opposite direction. The -axis is orthogonal to the channel walls, and the origin of the axes is such that the positions of the channel walls are and , respectively. A sketch of the system and of the coordinate axes is reported in Figure 1. The walls at and are isothermal at given temperatures and , where we assume that . The fluid has a uniform vertical upward stream wise velocity distribution at the channel entrance. Thus, the basic equations for steady, fully developed flow of a viscous, incompressible fluid and the heat assumed supplied to the surrounding fluid by an exothermic surface reaction [5] are subject to the boundary conditions: Reaction takes place only on the surface which is governed by the single first-order Arrhenius kinetics. The closure system by mass flux conservation is . The dimensionless governing equations for the problem are [5] subject to the boundary conditions: And the global conservation of mass at any cross section in the channel is .

3. Solutions Method

Equations (3) and (4) subject to (2) have been solved numerically by Pop et al. [5] using the bvpp4c function from MATLAB. Pop et al. [5] found the solutions are possible only for the finite range of . The range depend on constant temperature parameter, , for example, at , the solution domain are in . We solved again (3) and (4) subject to (5) via dsolve from MAPLE. The nonlinearities of (4) were handled by applying continuation, optional in dsolve. The continuation is utilized to maintain Newton method (inside package dsolve) reaching the convergence.

4. Results and Discussion

The effects of various levels , , and are presented in Figure 2. These figures show that velocity increases as , , or increase. Flow reversal is found near the colder wall () for . Increasing the Frank-Kamenetskii number or constant temperature parameter enhances the flow reversal. The interesting behavior of the reversed flow occurs at in Figure 2(c), where reversal is found near both of the walls. This phenomenon due to heat generated by exothermic surface reaction and being conducted away into the surrounding fluid and this behavior progress continuously cause the rate of the increase in the temperature to rise sharply in middle region of the channel.

The result in Figure 2 also implies that there exists a critical value such that, for , flow reversal occurs at . By setting at and utilizing dsolve with continuation, the can be obtained within the solution domain. The numerical value is plotted in Figure 3.

It was observed that increasing the Frank-Kamenetskii number reduces (Figure 3(a)), while reducing the temperature difference ratio increases the (Figure 3(b)). When , we recovered the critical value which had been calculated by Aung and Worku [3], Cheng et al. [4], El-Din [7], Barletta [12], and Chamkha et al. [9]. We can not obtain the critical value of the ratio of Grashof number and Reynolds number () with symmetrical wall temperatures when there is no chemical reaction in the fluid. The exists only in the range of constant .

5. Conclusions

Parameter zones for the occurrences flow reversal by mixed convection under the effect of chemical reaction in a vertical parallel-plate channel are presented. We can conclude that flow reversal adjacent to the cold wall is found to exist within the channel as the ratio of Grashof number and Reynolds number is above a threshold value. The exothermic chemical reaction is found to enhance the flow reversal and made flow reversal possible for symmetrical walls temperatures. The aspect of the nonuniqueness of the solution as found by Pop et al. [5] will be the focus of our next research undertaking.

Nomenclature

:Concentration of reactant
:Axial pressure gradient,
:Grashof number
:Dimensionless parameter,
:Gravitational acceleration
:Activation energy
:Preexponential factor
:Frank-Kamenetskii number
:Channel width
:Exothermicity factor
:Universal gas constant
:Reynolds number
:Temperature difference ratio
:Temperature
:Dimensionless velocity components in the -direction
,  :Dimensionless space coordinates.
Greek Symbols
:Thermal diffusivity
:Thermal expansion coefficient
:Density
:Dimensionless temperature
:Dynamic viscosity
:Kinematic viscosity.
Subscript
0:Reference value
:Left wall
:Right wall
:Critical.