`ISRN ThermodynamicsVolume 2012, Article ID 794741, 10 pageshttp://dx.doi.org/10.5402/2012/794741`
Research Article

## Unsteady/Steady Free Convective Couette Flow of Reactive Viscous Fluid in a Vertical Channel Formed by Two Vertical Porous Plates

1Department of Mathematics, Ahmadu Bello University, Zaria, Nigeria
2Department of Mathematics, Usmanu Danfodiyo University, P.M.B +234-2346, Sokoto, Nigeria

Received 28 March 2012; Accepted 15 May 2012

Academic Editors: M. Appell, G. L. Aranovich, A. Ghoufi, and B. Merinov

Copyright © 2012 Basant K. Jha 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

This paper presents unsteady as well as steady-state free convection Couette flow of reactive viscous fluid in a vertical channel formed by two infinite vertical parallel porous plates. The motion of the fluid is induced due to free convection caused by the reactive nature of viscous fluid as well as the impulsive motion of one of the porous plates. The Boussinesq assumption is applied, and the nonlinear governing equations of motion and energy are developed. The time-dependent problem is solved using implicit finite difference method, and steady-state problem is solved by applying regular perturbation technique. During the course of computation, an excellent agreement was found between the well-known steady-state solutions and transient solutions at large value of time.

#### 1. Introduction

The study of convection processes in porous channels is a well-developed field of investigation because of its importance to a variety of situations. For example, thermal insulations, geothermal system, surface catalysis of chemical reactions, solid matrix heat exchanger, microelectronic heat transfer equipment, porous flat collectors, coal and grain storage, petroleum industries, dispersion of chemical contaminants in various processes, nuclear waste material, and catalytic beds. Some of the examples mentioned above involve two or more fluids, multidimensional and unsteady flows .

The literature on the topic of unsteady as well as steady-state free convection flow problems is well surveyed by Singh et al. , Florio and Harnoy , Jha and Ajibade , Paul et al. , and Langellotto et al. . Muhuri  presented natural convection Couette flow between two vertical porous plates when one of the plates moves with uniform acceleration. Singh  investigated the motion of the fluid induced by the impulsive motion of one of the plates in the presence of convection currents due to asymmetric thermal condition on the boundaries. Abdulaziz and Hashim  considered free convection flow between porous vertical plates with asymmetric wall temperatures and concentrations and used homotopy analysis method to solve boundary-value problems. Fang [10, 11] analysed unsteady velocity and temperature profiles for Couette flow and pressure-driven flow in a channel with porous walls. Unsteady flow investigations with porous boundaries include the works of Wang et al. , Oxarango et al. , and Makinde and Ogulu . Makinde and Maserumule  presented thermal criticality and entropy analysis for a variable viscosity Couette flow.

In all previous studies the working fluid is considered as nonreactive viscous incompressible fluid. However, unsteady as well as steady-state free convection Couette flow of a reactive viscous fluid in a vertical channel formed by two vertical porous plates can be of importance for the design of equipment used in several types of engineering systems. There are many chemical reactions with important practical applications. A common configuration for such reactions is for the reactants to be made to flow over solid catalyst, with the reaction taking place on the surface of the catalyst. A full discussion of catalysis and description of many of its practical applications is given by Chaudhary and Merkin , Merkin and Chaudhary , and Chaudhary et al. . These authors studied flow configuration and assume that the reaction takes place only on the catalytic surface and can be represented schematically by the single first-order Arrhenius kinetics. According to these authors there is a three-way coupling between fluid flow, fluid/surface temperature, and reactant species concentration. Ayeni  and Dainton  studied a realistic mathematical description of thermal explosion which include the effects of Arrhenius temperature dependence with variable preexponential factor. Recently Jha et al.  investigated transient natural convection flow of reactive viscous flow in a vertical channel. Most recently Hazem Attia  studied the effect of suction and injection on unsteady Couette flow.

The objective of the present work is to analyze the effect of suction/injection on time dependent unsteady as well as steady-state free convection Couette flow of viscous reactive fluid in a vertical channel formed by two infinite vertical parallel porous plates.

#### 2. Mathematical Analysis

Consider the transient free convective Couette flow of viscous reactive fluid in a vertical channel formed by two infinite vertical parallel porous plates. The system under consideration is sketched in Figure 1. The flow is assumed transient and fully developed. The -axis is taken along the direction of the flow and parallel to the infinite vertical porous plates and -direction perpendicular to the flow direction. At time , the fluid and the porous plates are assumed to be same temperature , and there is no fluid motion. At , one of the porous plates starts moving with constant velocity and fluid generate heat due to its reactive nature causes to set up fluid motion inside the channel. In addition, at the same time the flow is subjected to suction of the fluid from one porous plate and at the same rate fluid is being injected through the other porous plate.

The fluid is assumed to be Newtonian and obeys the Boussinesq’s approximation. Under the previous assumptions the energy and momentum equations in dimensional form are where the initial fluid and wall temperature, dimensional temperature of fluid, the heat of reaction, the rate constant, the activation energy, the universal gas constant, the initial concentration of the reactant species, kinematic viscosity, and gap between the porous plates.

The initial and boundary conditions to be satisfied are Equations (1) and (2) can be nondimensionalized using the variables The physical quantities used in (3) are defined in the nomenclature.

Using the expressions (3) and (1) the dimensionless momentum and energy equations are while the initial and boundary conditions in dimensionless form are

#### 3. Analytical Solutions

The governing equations (4) presented in the previous section are highly nonlinear and exhibit no analytical solutions. The importance of analytical solutions which refers to steady free convection Couette flow of viscous reactive fluid relies on the chance to obtain nontrivial benchmarks to test the reliability of numerical codes developed for more complex physical situations. Analytical solutions are often an opportunity to inspect the internal consistency of mathematical models and of the approximations adopted, as well as to develop new theoretical results. The mathematical model representing the steady-state free convection Couette flow can be obtained by setting and in (4) to get The boundary conditions are In order to construct an approximate solution to (6) subject to (7), we employed a regular perturbation method by taking a power series expansion in the reactant consumption parameter : Substituting (8) into (6) and collecting the coefficients of like powers of , the solution of the governing equations is obtained as Steady-state skin frictions on the boundary plates is The steady-state rate of heat transfer on the boundary plates is The constants are defined in the appendix section.

#### 4. Numerical Solutions

To solve the time-dependent equation (4), the differential equations have been transformed into the corresponding finite difference equation. The procedure involves discretization of the momentum and energy equations into the finite difference equations at the grid point in which the time derivatives are approximated by the backward difference while the spatial derivatives are replaced by the central difference formula. The above equations are solved by Thomas algorithm by manipulating into a system of linear algebraic equations in the tridiagonal form.

In each time step, firstly the temperature field has been solved and then the velocity field is evaluated using the already known values of the temperature field. The process of computation is advanced until a solution is approached by satisfying the following convergence criterion: with respect to temperature and velocity fields. Here stands for the velocity or temperature fields, is the number of interior grid points, and is the maximum absolute value of .

#### 5. Results and Discussion

Figure 5: Temperature profile .
Figure 6: Temperature profile .
Figure 7: Temperature profile for .
Figure 8: Temperature profile .
Figure 9: Variation of skin friction at .
Figure 10: Variation of skin friction at .
Figure 11: Variation of skin friction at .
Figure 12: Variation of Nusselt number at .
Figure 13: Variation of Nusselt number at .
Figure 14: Variation of Nusselt number at .
Figure 15: Variation of Nusselt number at .
Figure 16: Variation of Nusselt number at .
Figure 17: Variation of Nusselt number at .

#### 6. Conclusion

The problem of unsteady as well as steady-state natural convection Couette flow of reactive viscous fluid in a vertical channel formed by two vertical porous plates has been presented. The velocity field and temperature field are obtained analytically by perturbation series method for steady free convective Couette flow of viscous reactive fluid in a vertical channel formed by two vertical porous plates and numerically by implicit finite difference technique for unsteady free convective Couette flow of viscous reactive fluid in a vertical channel formed by two vertical porous plates. Graphical results for the velocity, temperature, skin friction, and the Nusselt number variations were presented and discussed for various physical parametric values.

The main findings are as follows.(i)Skin friction is always higher in the case of air than water .(ii)It is also seen that increase in increases heat transfer on the porous plates.(iii)The heat transfer is higher at left porous plate where injection takes place in comparison to right plate where suction takes place .(iv)The introduction of suction/injection has distorted the symmetric nature of the flow.(v)During the course of computation it is observed that at steady state the effect of on the fluid flow is to increase the velocity and temperature as increases when injection is considered at while fluid velocity and temperature decreases as increases in the presence of suction at , which is not true in transient case (i.e., velocity and temperature decreases as increases).

Consider

#### Nomenclature

 𝐶𝑝: Specific heat of the fluid at constant pressure 𝐸: Activation energy 𝑔: Acceleration due to gravity 𝐻: Gap between the channels 𝑁𝑢0: Nusselt number at 𝑦=0 𝑁𝑢1: Nusselt number at 𝑦=1 Pr: Prandtl number 𝑄: Heat reaction parameter 𝑅: Universal gas constant 𝑠: Suction/injection parameter 𝑡′: Dimensional time 𝑡: Dimensionless time 𝑇: Dimensional temperature of the fluid 𝑇0: Initial temperature of the fluid 𝑢: Dimensional velocity of the fluid 𝑢: Dimensionless velocity of the fluid 𝑥: Dimensional coordinate parallel to channel 𝑦: Dimensional coordinate perpendicular to channel 𝑦: Dimensionless coordinate.
Greek Letters
 𝛽: Volumetric coefficient of thermal expansion 𝜌: Density of the fluid 𝜏0: Dimensionless skinfriction at 𝑦=0 𝜏1: Dimensionless skin friction at 𝑦=1 𝛼: Thermal diffusivity 𝜃: Dimensionless temperature 𝜆: Reactant consumption parameter 𝑣: Kinematic viscosity 𝜀: Activation energy parameter 𝛿: Grashof number.

#### Acknowledgment

The author A. K. Samaila is thankful to Usmanu Danfodiyo University, Sokoto, for financial support.

#### References

1. S. Mahmud and R. Andrew Fraser, “Flow, thermal, and entropy generation characteristics inside a porous channel with viscous dissipation,” International Journal of Thermal Sciences, vol. 44, no. 1, pp. 21–32, 2005.
2. A. K. Singh, H. R. Gholami, and V. M. Soundalgekar, “Transient free convection flow between two vertical parallel plates,” Heat and Mass Transfer/Waerme- und Stoffuebertragung, vol. 31, no. 5, pp. 329–331, 1996.
3. L. A. Florio and A. Harnoy, “Augmenting natural convection in a vertical flow path through transverse vibrations of an adiabatic wall,” Numerical Heat Transfer A, vol. 52, no. 6, pp. 497–530, 2007.
4. B. K. Jha and A. O. Ajibade, “Free convective flow between vertical porous plates with periodic heat input,” ZAMM Zeitschrift fur Angewandte Mathematik und Mechanik, vol. 90, no. 3, pp. 185–193, 2010.
5. T. Paul, B. K. Jha, and A. K. Singh, “Transient natural convection in a vertical channel,” International Journal of Applied Mechanics and Engineering, vol. 6, no. 4, pp. 913–922, 2001.
6. L. Langellotto, O. Manca, and S. Nardini, “Numerical investigation of transient natural convection in air in a convergent vertical channel symmetrically heated at uniform heat flux,” Numerical Heat Transfer A, vol. 51, no. 11, pp. 1065–1086, 2007.
7. P. K. Muhuri, “Flow formation in Couette motion in magnetohydrodynamics with suction,” Journal of the Physical Society of Japan, vol. 18, no. 11, pp. 1671–1675, 1963.
8. A. K. Singh, “Natural convection in unsteady Couette motion,” Defence Science Journal, vol. 38, no. 1, pp. 35–41, 1988.
9. O. Abdulaziz and I. Hashim, “Fully developed free convection heat and mass transfer of a micropolar fluid between porous vertical plates,” Numerical Heat Transfer A, vol. 55, no. 3, pp. 270–288, 2009.
10. T. Fang, “Further discussion on the incompressible pressure-driven flow in a channel with porous walls,” International Communications in Heat and Mass Transfer, vol. 31, no. 4, pp. 487–500, 2004.
11. T. Fang, “A note on the incompressible couette flow with porous walls,” International Communications in Heat and Mass Transfer, vol. 31, no. 1, pp. 31–41, 2004.
12. J. Wang, Z. Gao, G. Gan, and D. Wu, “Analytical solution of flow coefficients for a uniformly distributed porous channel,” Chemical Engineering Journal, vol. 84, no. 1, pp. 1–6, 2001.
13. L. Oxarango, P. Schmitz, and M. Quintard, “Laminar flow in channels with wall suction or injection: a new model to study multi-channel filtration systems,” Chemical Engineering Science, vol. 59, no. 5, pp. 1039–1051, 2004.
14. O. D. Makinde and A. Ogulu, “The effect of thermal radiation on the heat and mass transfer flow of a variable viscosity fluid past a vertical porous plate permeated by a transverse magnetic field,” Chemical Engineering Communications, vol. 195, no. 12, pp. 1575–1584, 2008.
15. O. D. Makinde and R. L. Maserumule, “Thermal criticality and entropy analysis for a variable viscosity Couette flow,” Physica Scripta, vol. 78, no. 1, Article ID 015402, 2008.
16. M. A. Chaudhary and J. H. Merkin, “Free-convection stagnation-point boundary layers driven by catalytic surface reactions: I the steady states,” Journal of Engineering Mathematics, vol. 28, no. 2, pp. 145–171, 1994.
17. J. H. Merkin and M. A. Chaudhary, “Free-convection boundary layers on vertical surfaces driven by an exothermic surface reaction,” Quarterly Journal of Mechanics and Applied Mathematics, vol. 47, no. 3, pp. 405–428, 1994.
18. M. A. Chaudhary, A. Liñan, and J. H. Merkin, “Free convection boundary layers driven by exothermic surface reactions: critical ambient temperatures,” Mathematical Engineering in Industry, vol. 5, no. 2, pp. 129–145, 1996.
19. R. O. Ayeni, “On the explosion of chain-thermal reaction,” Journal of the Australian Mathematical Society, vol. 24, pp. 194–202, 1982.
20. F. S. Dainton, Chain Reaction, An Introduction Wiley, New York, NY, USA, 1960.
21. B. K. Jha, A. K. Samaila, and A. O. Ajibade, “Transient free-convective flow of reactive viscous fluid in a vertical channel,” International Communications in Heat and Mass Transfer, vol. 38, no. 5, pp. 633–637, 2011.
22. A. Hazem Attia, “The effect of suction and injection on unsteady Couette flow,” Kragujevac Journal of Science, vol. 32, pp. 17–24, 2010.