# Flow Reversal of Fully Developed Mixed Convection in a Vertical Channel with Chemical Reaction

Authors: Habibis Saleh; Ishak Hashim; Sri Basriati
Journal: International Journal of Chemical Engineering
 (2013)

---

## 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.


---

## Body

## 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 [7–9] 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 x-axis is parallel to the gravitational acceleration vector g, but with the opposite direction. The y-axis is orthogonal to the channel walls, and the origin of the axes is such that the positions of the channel walls are y=0 and y=ℓ, respectively. A sketch of the system and of the coordinate axes is reported in Figure 1. The walls at y=0 and y=ℓ are isothermal at given temperatures T1 and T2, where we assume that T1≥T2. The fluid has a uniform vertical upward stream wise velocity distribution U0 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
(1)gβ(T-T0)-1ρdpdx+νd2udy2=0,αd2Tdy2+Qk0ae-E/RT=0,
subject to the boundary conditions:
(2)u(0)=0,u(ℓ)=0,T(0)=T1,T(ℓ)=T2,(T1≥T2).
Reaction takes place only on the surface which is governed by the single first-order Arrhenius kinetics. The closure system by mass flux conservation M is ∫0L‍UdY=M. The dimensionless governing equations for the problem are [5]
(3)d2UdY2=-GRΘ+A,(4)d2ΘdY2=-KFeΘ,
subject to the boundary conditions:
(5)U(0)=0,U(1)=0,Θ(0)=RT,Θ(1)=-RT.
And the global conservation of mass at any cross section in the channel is ∫01‍UdY=1.Figure 1
Schematic representation of the model.


## 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 KF. The range depend on constant temperature parameter, RT, for example, at RT=1, the solution domain are in 0≤KF<3.155. 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 levelsGR, KF, and RT are presented in Figure 2. These figures show that velocity increases as GR, KF, or RT increase. Flow reversal is found near the colder wall (Y=1) for GR≥50. Increasing the Frank-Kamenetskii number or constant temperature parameter enhances the flow reversal. The interesting behavior of the reversed flow occurs at RT=0.1 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.Plots ofU versus Y for different values of (a) GR where KF=1.5,RT=1 (b) KF where GR=300, RT=1, and (c) RT where GR=300, KF=3.1.
(a)
(b)
(c)The result in Figure2 also implies that there exists a critical value GRc such that, for GR>GRc, flow reversal occurs at Y=1. By setting dU/dY=0 at Y=1 and utilizing dsolve with continuation, the GRc can be obtained within the solution domain. The numerical value GRc is plotted in Figure 3.The critical valueGRc for the onset of a flow reversal with (a) increasing KF and (b) reducing RT.
(a)
(b)It was observed that increasing the Frank-Kamenetskii number reducesGRc (Figure 3(a)), while reducing the temperature difference ratio increases the GRc (Figure 3(b)). When KF=0, 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 (GR) with symmetrical wall temperatures RT=0 when there is no chemical reaction in the fluid. The GRc exists only in the range of constant RT≤0.05.


## 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.



---

*Source: 310273-2013-05-12.xml*