International Journal of Chemical Engineering

Volume 2013 (2013), Article ID 310273, 4 pages

http://dx.doi.org/10.1155/2013/310273

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

^{1}School of Mathematical Sciences, Universiti Kebangsaan Malaysia, 43600 Bangi, Selangor, Malaysia^{2}Solar Energy Research Institute, Universiti Kebangsaan Malaysia, 43600 Bangi, Selangor, Malaysia^{3}Mathematics Sciences, UIN Sultan Syarif Kasim Riau, Panam, Pekanbaru 28293, Indonesia

Received 12 March 2013; Accepted 21 April 2013

Academic Editor: Donald L. Feke

Copyright © 2013 Habibis Saleh et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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

#### References

- L. N. Tao, “On combined free and forced convection in channels,”
*Journal of Heat Transfer*, vol. 82, pp. 233–238, 1960. View at Google Scholar - S. Habchi and S. Acharya, “Laminar mixed convection in a symmetrically or asymmetrically heated vertical channel,”
*Numerical Heat Transfer A*, vol. 9, no. 5, pp. 605–618, 1986. View at Google Scholar · View at Scopus - W. Aung and G. Worku, “Theory of fully developed, combined convection including flow reversal,”
*Journal of Heat Transfer*, vol. 108, pp. 299–304, 1986. View at Google Scholar - C. H. Cheng, H. S. Kou, and W. H. Huang, “Flow reversal and heat transfer of fully developed mixed convection in vertical channels,”
*Journal of Thermophysics and Heat Transfer*, vol. 4, no. 3, pp. 375–383, 1990. View at Google Scholar · View at Scopus - I. Pop, T. Grosan, and R. Cornelia, “Effect of heat generated by an exothermic reaction on the fully developed mixed convection flow in a vertical channel,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 15, no. 3, pp. 471–474, 2010. View at Publisher · View at Google Scholar · View at Scopus - E. M. Sparrow, G. M. Chrysler, and L. F. Azevedo, “Observed flow reversals and measuredpredicted nusselt numbers for natural convection in a one-sided heated vertical channel,”
*Journal of Heat Transfer*, vol. 106, no. 2, pp. 325–332, 1984. View at Google Scholar · View at Scopus - M. M. S. El-Din, “Fully developed forced convection in a vertical channel with combined buoyancy forces,”
*International Communications in Heat and Mass Transfer*, vol. 19, no. 2, pp. 239–248, 1992. View at Google Scholar · View at Scopus - A. J. Chamkha, “On laminar hydromagnetic mixed convection flow in a vertical channel with symmetric and asymmetric wall heating conditions,”
*International Journal of Heat and Mass Transfer*, vol. 45, no. 12, pp. 2509–2525, 2002. View at Publisher · View at Google Scholar · View at Scopus - A. J. Chamkha, T. Groşan, and I. Pop, “Fully developed mixed convection of a micropolar fluid in a vertical channel,”
*International Journal of Fluid Mechanics Research*, vol. 30, no. 3, pp. 251–263, 2003. View at Publisher · View at Google Scholar · View at Scopus - T. Grosan, R. Pop, and I. Pop, “Thermophoretic deposition of particles in fully developed mixed convection flow in a parallel-plate vertical channel,”
*Heat and Mass Transfer*, vol. 45, no. 4, pp. 503–509, 2009. View at Publisher · View at Google Scholar · View at Scopus - E. Magyari, “Thermophoretic deposition of particles in fully developed mixed convection flow in a parallel-plate vertical channel: the full analytical solution,”
*Heat and Mass Transfer*, vol. 45, no. 11, pp. 1473–1482, 2009. View at Publisher · View at Google Scholar · View at Scopus - A. Barletta, “Laminar mixed convection with viscous dissipation in a vertical channel,”
*International Journal of Heat and Mass Transfer*, vol. 41, no. 22, pp. 3501–3513, 1998. View at Publisher · View at Google Scholar · View at Scopus