Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2016 / Article

Research Article | Open Access

Volume 2016 |Article ID 2568785 |

Syed Tauseef Mohyud-Din, Naveed Ahmed, Umar Khan, Asif Waheed, Saqib Hussain, Maslina Darus, "On Combined Effects of Heat Transfer and Chemical Reaction for the Flow through an Asymmetric Channel with Orthogonally Deformable Porous Walls", Mathematical Problems in Engineering, vol. 2016, Article ID 2568785, 10 pages, 2016.

On Combined Effects of Heat Transfer and Chemical Reaction for the Flow through an Asymmetric Channel with Orthogonally Deformable Porous Walls

Academic Editor: Francesco Franco
Received19 Nov 2015
Revised26 Jan 2016
Accepted21 Feb 2016
Published17 Mar 2016


The combined effects of heat transfer and chemical reaction are studied for the flow through a semi-infinite asymmetric channel with orthogonally deformable porous walls. The similarity transforms have been used to reduce the conservation laws to a corresponding system of nonlinear ordinary differential equations. The resulting equations are solved, both analytically and numerically, by using Homotopy Analysis Method (HAM) and the fourth-order Runge-Kutta (RK-4) method, respectively. The convergence of the analytical solution is assured through the so-called total squared residual error analysis. The optimal values of auxiliary parameters are obtained by minimizing the total squared residual error.

1. Introduction

The flow between porous channels can rightly be considered as one of the most important problems in fluid mechanics. Its practical applications in diversified fields of science and technology have been the basic reason why scientists from all over the world have tried to understand these flows in a better way. These applications include uniformly distributed irrigation, interbody biological transport, ablation processes, sublimation mechanisms, propellant combustion, and industrial cleansing systems.

Formative work regarding these kinds of flows can be traced back to Berman [1], who initiated a mathematical model to describe the flows between parallel plates. He established that the normal component of the velocity must be independent of streamwise coordinate. This in fact enabled him to reduce Navier-Stokes equations to a single, fourth-order, nonlinear ordinary differential equation. Following his footsteps, numerous studies have been carried out later on. The detailed history about these investigations can be found in [2].

Heat transfer coupled with chemical reaction plays an important role in many fields of science and engineering. Its practical applications may include combustion systems, atomic reactor safety, dying of cloths, metallurgy, and chemical engineering. A chemical reaction is said to be of first order if the rate of reaction varies directly with respect to the concentration. Almost in all chemical industries, a relatively inexpensive raw material is put through some specifically designed chemical processes to obtain high value products. In most of the situations, these chemical processes are based on chemical reactions in the presence of heat transfer. Bridges accompanied by Rajagopal [3] investigated the pulsatile flow of a chemically reacting fluid, whose viscosity was dependent on the concentration of the constituents. The flow was governed by a convection-reaction-diffusion equation and the velocity gradient, which could thicken or thin the fluid. The effects of chemical reaction and the space porosity on mixed convective MHD flow in an asymmetric channel were investigated by Srinivas and Muthuraj [4]. The details of heat transfer analysis can be found in [511].

In a very recent article, Reddy et al. [12] discussed the combined effects of heat and chemical reaction on an asymmetric laminar flow between slowly expanding and contracting walls. To make their perturbation solution valid, they had to impose physical constraints such as lower wall deformation rate and weak permeability. Besides, the presence of other nonperturbed physical parameters makes their solution more vulnerable as the perturbation solution strictly relies on the existence of very small parameters in equation.

In this study, we have tried to remove those physical constraints, which have been imposed in [13], only to insure the convergence of their perturbation solution. A probable invalidity of their solution for concentration profile has also been pointed out even for smaller values of wall deformation rate and the permeation Reynolds number.

We have used a well-known analytical method called Homotopy Analysis Method (HAM) [1418] to guarantee the convergence of the analytical solution. To further ensure the validity of the solution, a numerical solution (fourth-order Runge-Kutta coupled with shooting method) has also been sought. The aforesaid numerical method has effectively been used in several studies [1921]. We have compared the analytical and numerical results and have found an excellent agreement between them. Finally, the effects of physical parameters on temperature and concentration profiles are displayed and analyzed with the help of graphs accompanied by comprehensive discussions.

2. Mathematical Analysis

The laminar flow of a viscous incompressible fluid is taken into account, in a rectangular domain of infinite length, and having two deformable porous walls. The fluid may enter or exit, during the successive contractions/expansions, through the permeable walls. The leading end of the container is closed with a solid but flexible membrane, which deforms incorporation with the movement of the walls. A narrow gap assumption has been imposed; that is, the height () of the channel is taken to be very small as compared to its width. It enables us to confine the whole problem in to a half domain and a planer cross section of the simulating domain is presented in Figure 1.

Both the walls are assumed to have different (from each other) permeability and to deform uniformly at a time dependent rate (=). The origin is taken at the center of channel as shown in Figure 1. The temperature and concentration both possess higher values at the lower wall. Following the aforementioned assumptions, the governing equations take the following form [12]:In the above equations, the velocity components in and directions are denoted by and , respectively. Moreover, , and denote pressure, density, specific heat, kinematic viscosity, temperature, thermal conductivity, concentration, coefficient of mass diffusivity, and time, respectively. Further, represents the first-order chemical reaction rate taken to be positive for destructive reaction, negative for constructive reaction, and zero for no chemical reaction at all.

Suitable boundary conditions for the problem are as follows [12]:where and denote the measure of permeability of lower and upper wall, respectively. Temperature and concentration are taken to be greater on the lower wall as compared to the upper wall; that is, and .

We can abridge the governing system of equations by eliminating the pressure terms from (2) and (3) and utilizing (1). Introducing vorticity , we get whereConservation of mass enables us to develop a similarity solution with respect to as follows:where represents . Using (9), we obtain Using (10) in (7), we have while is the nondimensional wall expansion or contraction rate, taken to be positive for expansion.

We can set by taking to be a constant or a quasi constant in time [3].

The auxiliary conditions can also be transformed as where and are the permeation Reynolds numbers taken to be positive for injection.

Aforementioned formulation can be made nondimensional by introducing the following normalizing parameters: Consequently, we have with where . The primes in the above equations represent the differentiation with respect to .

The equations describing the temperature and concentration are [12] where and are the power law indices of temperature and concentration, respectively. Besides, and are the fluid constants.

Substituting (16) in (4) and (5) and by using and , we have dimensionless equations for the temperature and concentration distributions as follows:The boundary conditions become in the above equation denotes Prandtl number, is Schmidt number, is the chemical reaction parameter, and . It is also appropriate to mention that the similarity transforms given in (16) are only valid if . We still retain the terms involving and in (17) to make them comparable with the work done by Reddy et al. [12]. However, the forthcoming analysis is done after correcting the values of these parameters.

3. Homotopy Solution

Pursuing the technique suggested by Liao [13, 14], we can define the following initial guesses: The linear operators can be chosen asThese operators satisfy the following properties: where are the arbitrary constants.

3.1. Zero-Order Deformation Problem

Representing as the embedding parameter, zero-order deformation problem can be constructed as where , and are the nonzero auxiliary parameter.

The nonlinear operators are

3.2. th-Order Deformation Problem

The th-order problems satisfy whereAlso,For and , we have Using Taylor’s series in terms of , one can get Substituting in the above equations, we obtain where , and can be obtained by solving set of (24) using computer software, Mathematica. Substituting these back into (29) gives us the final solution.

3.3. Convergence of the Solution

The solution obtained in the preceding subsection contains the auxiliary parameters , and ; these are very essential to control and adjust the convergence of the series solution. To determine the optimal values of these parameters so that the so-called total residual error is minimized, we have used a Mathematica package named BVPh 2.0 [14]. This is in accordance with the suggestions given by one of the respected anonymous reviewers. The resulting optimal values along with the total residual error, corresponding to different sets of physical parameters, are displayed in Table 1.

αTotal error


The convergence of the solution, for a specific set of physical parameters, is shown in Figure 2. A decline of total residual error with increasing the order of HAM iterations is evident. As the residual approaches zero quickly, it ensures the convergence of the solution.

The optimal values of the auxiliary parameters, which are far from , reveal that the solutions obtained by other methods such as homotopy perturbation method and Adomian’s decomposition method are not suitable for this problem as they are suspected to diverge. The possible divergence is also an issue with the regular perturbation method. The solution presented in [12] is only valid in very limiting cases where the physical parameters are taken to be too small. Besides, they have only perturbed the wall deformation parameter () and permeation Reynolds number . As there are other several physical parameters involved, the solution portrayed by [12] becomes more vulnerable. The solutions we obtained in this study are free from existence of small or large parameters, and also the convergence is assured by the use of optimal values of , and . Therefore, it seems right to proclaim that ours are the improved results.

The results are also compared with the numerical results obtained by using well-known fourth-order Runge-Kutta method. An excellent agreement between both the results has been documented in the form of Table 2.

, , , , , , and

HAM Numerical HAM Numerical HAM Numerical