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

Authors: Basant K. Jha; Ahmad K. Samaila; Abiodun O. Ajibade
Journal: ISRN Thermodynamics
 (2012)

---

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


---

## Body

## 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 [1].The literature on the topic of unsteady as well as steady-state free convection flow problems is well surveyed by Singh et al. [2], Florio and Harnoy [3], Jha and Ajibade [4], Paul et al. [5], and Langellotto et al. [6]. Muhuri [7] presented natural convection Couette flow between two vertical porous plates when one of the plates moves with uniform acceleration. Singh [8] 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 [9] 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. [12], Oxarango et al. [13], and Makinde and Ogulu [14]. Makinde and Maserumule [15] 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 [16], Merkin and Chaudhary [17], and Chaudhary et al. [18]. 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 [19] and Dainton [20] studied a realistic mathematical description of thermal explosion which include the effects of Arrhenius temperature dependence with variable preexponential factor. Recently Jha et al. [21] investigated transient natural convection flow of reactive viscous flow in a vertical channel. Most recently Hazem  Attia [22] 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 Figure1. The flow is assumed transient and fully developed. The x′-axis is taken along the direction of the flow and parallel to the infinite vertical porous plates and y′-direction perpendicular to the flow direction. At time t′≤0, the fluid and the porous plates are assumed to be same temperature T0, and there is no fluid motion. At t′>0, one of the porous plates (y′=H) 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.Figure 1
Schematic diagram of the problem.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(1)∂u′∂t′-v0∂u′∂y′=υ∂2u′∂y′2+gβ(T′-T0),∂T′∂t′-v0∂T′∂y′=α∂2T′∂y′2+QCoAρCpexp⁡(-ERT′),
where T0 the initial fluid and wall temperature, T′ dimensional temperature of fluid, Q the heat of reaction, A the rate constant, E the activation energy, R the universal gas constant,  Co the initial concentration of the reactant species, v kinematic viscosity, and H gap between the porous plates.The initial and boundary conditions to be satisfied are(2)t′≤0:u′=0,T′=T0for0≤y′≤h,t′>0{u′=0,T′=T0,aty′=0,u′=U,T′=T0,aty′=H.
Equations (1) and (2) can be nondimensionalized using the variables
(3)y=y′H,t=t′vH2,u=u′U,θ=ERT02(T′-T0),λ=QCoAEH2κRT02exp⁡(-ERT0),ε=RT0E,s=V0Hυ,δ=gβRT02h2EUυ,Pr=υα.
The physical quantities used in (3) are defined in the nomenclature.Using the expressions (3) and (1) the dimensionless momentum and energy equations are
(4)∂u∂t-s∂u∂y=∂2u∂y2+δθ,∂θ∂t-s∂θ∂y=1Pr∂2θ∂y2+λPrexp⁡(θ1+εθ),
while the initial and boundary conditions in dimensionless form are(5)t≤0:u=θ=0for0≤y≤1,t>0:{u=0,θ=0aty=0,u=1,θ=0aty=1.


## 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 ∂u/∂t=0 and ∂θ/∂t=0 in (4) to get
(6)d2udy2+sdudy+δθ=0,d2θdy2+Prsdθdy+λexp⁡(θ1+εθ)=0.
The boundary conditions are
(7)u=0,θ=0aty=0,u=1,θ=0aty=1.
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 λ:
(8)θ=∑i=0∞θiλi,u=∑i=0∞uiλi.
Substituting (8) into (6) and collecting the coefficients of like powers of λ, the solution of the governing equations is obtained as
(9)θ=λ[1sPr(1-exp⁡(-sPry))(1-exp⁡(-sPr))-ysPr]+λ2{12sPr(y2sPr-2y(sPr)2+2(sPr)3)-y(sPr)2((1+exp⁡(-sPry))(1-exp⁡(-sPr)))+1(sPr)31(1-exp⁡(-sPr))+k1sPr12sPr+k2exp⁡(-sPry)(y2sPr-2y(sPr)2+2(sPr)3)},u=[1-exp⁡(-sy)][1-exp⁡(s)]+λ[δ{12sPr(y2s-2s{ys-1s2})-1z1{ys-1s2}-1z1exp⁡(-sPry)s2Pr(1-Pr)}+c1s+c2exp⁡(-sy){ys-1s2}{12sPr(y2s-2s{ys-1s2})]+λ2[G1s-δ6(sPr)2{y3s-3y2s2+6ys3-6s4}+δ{12(sPr)3+12(sPr)2(1-exp⁡(-sPr))}×{y2s-2ys2+2s3}-δ(sPr)4(ys-1s2)-δ(sPr)3(1-exp⁡(-sPr))×{ys(1-Pr)-1s2(1-Pr)2}exp⁡(-sPry)-δ(sPr)41(1-exp⁡(-sPr))exp⁡(-sPry)s(1-Pr)-δ(sPr)3(1-exp⁡(-sPr)){ys-1s2}-k1sPr(ys-1s2)+k2exp⁡(-sPry)s2Pr(1-Pr)+G2exp⁡(-sy){y3s-3y2s2+6ys3-6s4}k2exp⁡(-sPry)s2Pr(1-Pr)].
Steady-state skin frictions on the boundary plates is
(10)τ0=dudy|y=0=s[1-exp⁡(-s)]+λ[δ(Prz1s(1-Pr)-1s3Pr)-c2s]-λ2[δs5(1Pr2+1Pr3+1Pr4)+δs4Pr2(1-exp⁡(-sPr))(1+1(1-Pr)2+1Pr)+k1s2Pr-k2s(1-Pr)-G2s1s3Pr],τ1=dudy|y=1=[sexp⁡(-s)][1-exp⁡(s)]+λ[k2sPrexp⁡(-sPr)s2Pr(1-Pr)δ{12sPr(2s-2s2)+sPrz1exp⁡(-sPr)s2Pr(1-Pr)}-c2sexp⁡(-s){12sPr(2s-2s2)+sPrz1exp⁡(-sPr)s2Pr(1-Pr)}12sPr][1-exp⁡(s)]2s2[sexp⁡(-s)][1-exp⁡(s)]-λ2[k2sPrexp⁡(-sPr)s2Pr(1-Pr)δ6(sPr)2{3s-6s2+6s3}+δ{12(sPr)3+12(sPr)2(1-exp⁡(-sPr))}×{2s-2s2}-δs5(Pr)4+δsPrexp⁡(-sPr)(sPr)3(1-exp⁡(-sPr))1s(1-Pr)+δ(sPr)41(1-exp⁡(-sPr))sPrexp⁡(-sPr)s(1-Pr)-δs2(Pr)3(1-exp⁡(-sPr))-k1s2Pr-k2sPrexp⁡(-sPr)s2Pr(1-Pr)-G2sexp⁡(-s)δ6(sPr)2]{12sPr(2s-2s2)+sPrz1exp⁡(-sPr)s2Pr(1-Pr)}sPrz1.
The steady-state rate of heat transfer on the boundary plates is
(11)Nu0=dθdy|y=0=λ[1(1-exp⁡(-sPr))-1sPr]-λ2[1(sPr)3+2(sPr)2(1-exp⁡⁡(-sPr))+sPrk2],Nu1=dθdy|y=1=λ[exp⁡(-sPr)(1-exp⁡(-sPr))-1sPr]+λ2[1(sPr)3(sPr-1)-1(sPr)2(1-exp⁡(-sPr))×{1+(1-sPr)exp⁡(-sPr)}1sPr-sPrk2exp⁡(-sPr)1(sPr)3].
The constants  c1,c2,k1,k2,A,B,A1,B1,G1,G2 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 (i,j) 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:(12)∑|Ai,j+1-Ai,j|M|A|max⁡<10-5
with respect to temperature and velocity fields. Here Ai,j stands for the velocity or temperature fields, M is the number of interior grid points, and |A|max is the maximum absolute value of Ai,j.


## 5. Results and Discussion

The basic parameters that governed this flow are the reactant consumption parameter(λ), the activation energy parameter (ε), the Grashof number (δ), the Prandtl number (Pr) which is inversely proportional to the thermal diffusivity of the working fluid, and suction/injection parameter (s), which were simultaneously applied each to opposite porous plates of the channel at the same rate. For the purpose of discussion, some numerical calculations are carried out for dimensionless velocity  (u), temperature (θ), skin friction (τ), and the rate of heat transfer (Nu). In the entire computation the value of Grashof number (δ) is taken to be unity. The effect of various parameter on the flow are presented in graphical form in Figures 2 to 17 so as to clearly reveal the influence of each governing parameter on the flow. Figure 2 describes the influence of s and t on the velocity profile within the channel. This figure reveals that as time increases, the velocity increases. This is the consequence of temperature increase with time. This figure also reflects that when suction takes place on the porous plate at  y=0(s>0) fluid velocity is higher compared to a situation when injection takes place on the porous plate at y=0(s<0). This is physically true since suction on the porous plate y=0 happens concurrently with injection on the moving porous plate (y=1) which acts in support of the motion of the porous plate (y=1), and hence an increase in fluid velocity is the consequence. On the other hand (s<0) means injection on the plate y=0 which translates to a corresponding suction on the plate (y=1) and this act against the influence of the impulsive motion of the porous plate. The effect of reactant consumption parameter and time on velocity profile is depicted in Figure 3. This figure shows that as λ increases the velocity increases. This physical fact can be explained, as reaction increases in the fluid, the temperature increases which causes velocity increase. This figure also reflects that the time required reaching steady-state velocity increases as λ increases. From the figure it is observed that the influence of λ on velocity is almost negligible when the time is small. However, at steady state, the velocity is observed to increase as reactant consumption parameter increases. Figure 4 reflects unsteady as well as steady-state velocity profiles for different value of nondimensional time and Prandtl number. It is evident from the figure that as time increases the velocity increases. Figure 4 also shows that the Prandtl number has insignificant effect on unsteady and steady-state velocity. Figure 5 shows the influence of s on the temperature profiles of water (Pr=7.0) within the channel. It is seen from the figure that temperature increases with time. In addition, in case of suction at y=0, velocity is maximum near the left plate (y=0), and in case of injection the maximum is shifted towards the right porous plate (y=1), and hence we conclude that there is no symmetry in the temperature profile in the presence of suction/injection. Figure 6 describes the behaviour of unsteady as well as steady-state temperature profiles as a function of time for different value of s  (s=-0.5,-1.5). In this figure, temperature is observed to be high near the plate with injection (y=1). Figure 7 reveals the influence of λ on the temperature profile. It is reflected from this figure that as time increases the temperature increases. It is evident from Figure 7 that the time required to reach steady state temperature increases as reactant consumption parameter increases. Variation of temperature as a function of reactant consumption parameter and time is illustrated in Figure 8. It is clear from the figure that as λ increases the temperature increases. It is also evident that the time required to reach steady state is strongly dependent on λ (i.e., time required to reach steady state is high for high value of λ). Variation of skin friction for water (Pr=7.0) and air (Pr=0.71) at the plate y=0 and y=1 is shown in Figures 9, 10, and 11. These figures reflect that as time increases the skin friction increases and finally attains its steady-state value. The effect of suction/injection parameter and the nondimensional time on the skin friction is shown in Figure 9. It is clear from the figure that the skin friction is directly proportional to nondimensional time and suction/injection parameter until steady-state value is attained. This figure also reveals monotonic increase in the skin friction as sincreases. This is due to the fact that as s  increases the velocity increases and consequently there is higher skin friction at the boundary. Figures 10, and 11 show the influence of time and reactant consumption parameter on the skin-friction at y=0 and y=1. From the two figures it is clear that the skin-friction increases as time increases until steady-state value is attained. Figures 9, 10, and 11 also reflect that the skin friction is always higher in the case of air (Pr=0.71) than water (Pr=7.0). The rate of heat transfer which is expressed as Nusselt number (Nu) at y=0 and y=1 is shown in Figures 12, 13, 14, 15, 16, and 17. Figures 12 and 13 show the influence of s on the rate of heat transfer. Figure 12 reveals that the rate of heat transfer increases as nondimensional time and suction/injection parameter increases until a steady-state value is attained. Figure 13 shows monotonic decrease in rate of heat transfer as s increases. It is reflected from this figure that both the transient and steady-state rate of heat transfer coincide at large values of s and t. Figures 14 and 15 show the variation of reactant consumption parameter on the rate of heat transfer. These two figures reveal that there is a monotonic increase in the rate of heat transfer at (y=0,y=1) corresponding to water (Pr=7.0) and air (Pr=0.71). This is due to the higher temperature gradient between fluid and porous plates. It is also observed from the figure that as time increases the rate of heat transfer increases until steady-state Nusselt number is reached. Figures 16 and 17 show the influence of reactant consumption parameter with time on the Nusselt number. It is reflected from these figures that the rate of heat transfer is directly proportional to the reactant consumption parameter. These two figures also show that as time increase the rate of heat transfer increases until steady-state value is reached.Figure 2
Velocity profile(ε=0.01,Pr=7.0,λ=1.0).Figure 3
Velocity profile(ε=0.01,Pr=7.0,s=0.5).Figure 4
Velocity profile(ε=0.01,s=0.5,λ=1.0).Figure 5
Temperature profile(ε=0.01,Pr=7.0,λ=1.0).Figure 6
Temperature profile(ε=0.01,Pr=7.0,λ=1.0).Figure 7
Temperature profile for(ε=0.01,s=0.5,Pr=7.0).Figure 8
Temperature profile(ε=0.01,Pr=0.71,s=0.5).Figure 9
Variation of skin friction(ε=0.01,λ=1.0) at y=0.Figure 10
Variation of skin friction(ε=0.01,s=0.5,λ=1.0) at y=0.Figure 11
Variation of skin friction(ε=0.01,s=0.5) at y=1.Figure 12
Variation of Nusselt number(ε=0.01,λ=1.0,Pr=7.0) at  y=0.Figure 13
Variation of Nusselt number(ε=0.01,λ=1.0,Pr=7.0) at y=1.Figure 14
Variation of Nusselt number(s=0.5,ε=0.01,Pr=7.0) at y=0.Figure 15
Variation of Nusselt number(s=0.5,ε=0.01,Pr=0.71) at y=0.Figure 16
Variation of Nusselt number(s=0.5,ε=0.01,Pr=7.0) at y=1.Figure 17
Variation of Nusselt number(s=0.5,ε=0.01,Pr=0.71) at y=1.


## 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(Pr=0.71) than water (Pr=7.0).(ii)
It is also seen that increase inλ increases heat transfer on the porous plates.(iii)
The heat transfer is higher at left porous platey=0 where injection takes place in comparison to right plate where suction takes place y=1.(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 ofPr on the fluid flow is to increase the velocity and temperature as Pr increases when injection is considered (s<0) at y=0 while fluid velocity and temperature decreases as Pr increases in the presence of suction (s>0) at y=0, which is not true in transient case (i.e., velocity and temperature decreases as Pr increases).



---

*Source: 794741-2012-07-16.xml*