This work is licensed under the Creative Commons Attribution 4.0 Public License.
Introduction
Viscous flow through deformable porous material along with coupled fluid movements is studied widely by many researchers using analytical and numerical techniques. It has many significant applications such as flow through packed beds and ion-exchange beds, energy extraction from the geothermal regions, solid filtration from liquids. The fluid flows through deformable porous materials are driven by the several biological applications like the study of soft tissue, articular cartilage, skin and arterial walls etc. Most of the research works were available inundeformable porous media. Further the work done on deformable porous media is very limited. The coupled phenomenon of fluid flow and solid deformation in permeable materials is classical problems in the fields of geomechanics and biomechanics. The analysis of hemodynamic effect of the endothelial glycocalyx is an important application for interaction of free flow and deformable porous media. In sight of these applications, Terzaghi [11] initiated the study of flow through deformable porous materials. Biot [2] extended the work by considering soil consolidation and acoustic propagation. Atkin [3] and Bowen [4] presented some significant results to the theory of mixture. Later, Jayaraman [5] discussed water transport throughthe artery wall. Mow et al. [6] analyzed fluid flow with mechanical properties of articular cartilage. Barry et al.[7] investigated fluid flow over a thin deformable porous layer. Klubertanz et al. [8] discussed multiphase flow in deformable porous media. Ranganatha and Siddagamma [9] investigated the flow of Newtonian fluid trough a deformable porous channel walls. Wen et al. [10] investigated the characteristics of shear flows in deformable porous surface layer through cylindrical tube. Sreenadh et al. [11] analyzed the MHD free surface flow of a Jeffery fluid over a deformable porous layer. Later, Sreenadh et al. [1212, 13] have studied non-Newtonian fluid flows in deformable porous media including the effects of heat transfer and magnetic field.
Woods [14] studied thermodynamics of fluid systems. Second law analysis in heat transfer energy is explained by Bejan [15]1. Paoletti et al. [16] presented calculation of energetic losses in compact heat exchanger passages. Bejan [17] also studied entropy generation through heat and fluid flow. Ajibade et al. [18] analyzed the entropy generation with effect of suction/injection. Makinde et al.[19] analyzed of inherent irreversibility in a variable viscosity of Couette flow through permeable walls with MHD effects. Ramana Murthy et al. [20] made second law analysis for Poiseuille flow of immiscible micropolar fluids in a channel. Rashidi et al. [21] analyzed entropy generation in MHD and slip flow over a rotating porous disk. Das et al. [22] discussed entropy generation in a rotating Couette flow with suction/injection. Vikaskumar et al. [23] investigated entropy generation in Poiseuille flow through a channel partially filled with a porous material. Rashidi et al. [24] studied entropy generation analysis for stagnation point flow through a porous medium over a permeable stretching surface.
Motivated by the above studies, the present paper deals with the entropy generation in Couette flow with suction and injection through a deformable porous channel. The expressions for the fluid velocity, solid displacement and temperature distribution are found analytically. The effect of various physical parameters on flow characteristics are analyzed through graphs and tables.
Flow geometry and governing equations
Consider the Couette flow of a viscous fluid bounded by rigid walls through a deformable porous medium. The lower impermeable wall is stationary whereas the upper one is moving with constant velocity U0 . The width of the channel is taken as h . The fluid is injected into the porous channel perpendicular to the lower wall with a constant velocity V and is sucked out of the upper wall with same velocity V . The fluid velocity and the displacement in the porous region are assumed respectively as (v, 0, 0) and (u, 0, 0).
$$\begin{array}{}
\displaystyle
At \,\,\,\, y = 1:u = 0,\,v = {U_0},\,\theta = 1 + m
\end{array}
$$
Solution of the problem
The governing equations (5)-(8) are coupled differential equations that can be solved analytically by using the boundary conditions (9) and (10). The solid displacement, the fluid velocity and the temperature distribution are obtained as
The entropy generation equation (16) consists of two terms, the first term on the right hand side is the entropy generation due to the heat transfer across a finite temperature difference and the second term is the local entropy due to viscous dissipation.
The dimensionless entropy generation number may be determined by the following relationship:
where $\begin{array}{}
\displaystyle
Br = \frac{{2{\mu _a}{U^2}}}{{{K_0}({T_w} - {T_0})}}
\end{array}$ the Brinkman number and $\begin{array}{}
\displaystyle
\Omega = \frac{{{T_w} - {T_0}}}{{{T_0}}}
\end{array}$ the non-dimensional temperature.
In view of this Paoletti et al. [16] has defined an alternative irreversibility distribution parameter in terms of Bejan number as
where $\begin{array}{}
\displaystyle
{N_s} = {N_1} + {N_2} , {N_1} = {\left( {\frac{{d\theta }}{{dy}}} \right)^2} , {N_2} = \frac{{Br}}{\Omega }{\left( {\frac{{dv}}{{dy}}} \right)^2}
\end{array}$ and $\begin{array}{}
\displaystyle
\Phi = \frac{{{N_2}}}{{{N_1}}}
\end{array}$ is the irreversibility ratio.
Results and discussions
In this section, we analyzed the Couette flow through a horizontal deformable porous layer. The effects of various pertinent parameters such as the volume fraction of the fluid ϕ, drag δ, Brinkman number Br, suction velocity V and heat source parameter β on solid displacement, fluid velocity, temperature distribution, entropy generation and Bejan number profiles are discussed through graphs. The numerical computations are carried out by using Matlab software. In this study P = –8, ϕ = 0.6, δ = 1, β = 1, Br = 1, V = 0.5, m = 1, Re = 1, Fr = 1 and U0 = 1 are used for numerical computation. The variation of solid displacement u with y is calculated from Eq. (11) for different values of δ, U0, ϕ and V which are represented in figs. 2 - 5.From fig. 2, it is noticed that, an increase in the drag increases the solid displacement in the porous channel.
Figure 3 illustrates that that increasing upper plate velocity causes more solid displacement in the porous channel. It is noticed that the solid displacement reduces with the increasing ϕ and v. That is, the increase in the volume fraction of the fluid and suction/injection velocity decreases the solid displacement in the porous channel which is presented in figs. 4 and 5 respectively. The variation of fluid velocity v with y is calculated from Eq. (12) for various values of ϕ, U0, P, δ and V is shown in figs.6 - 10 respectively. It is found that the velocity increases with the increasing ϕ, U0, and P. This shows that, the increase in the movement of the upper plate results considerable increase is the fluid velocity in the porous layer along with the favorable pressure gradient. The same phenomenon is noticed even when the volume fraction of the fluid increases in the deformable porous channel. It is found that the velocity reduces with the growing values of δ and V. That is, the fluid velocity gets reduced due to an increase in the drag and suction/injection.
The variation of temperature θ with y is calculated from Eq. (13) for different values of U0, Br, m, β, ϕ, δ and v which are shown in figs. 11 - 17. From fig. 11, it is noticed that the temperature profiles increases with increasing upper plate velocity. The temperature increases with increasing Brinkman number Br in porous channel is represented in fig. 12. That is the larger values of Brinkman number Br are indicative of larger frictional heating system so increase in temperature. The effect of m on temperature is shown in fig.13. The increase in m enhances the temperature in the deformable porous channel i.e., m = 0 implies both the plates are maintained at constant temperature, m = 1 and m = 2 represents heating of the plates y = 0 . Hence, the temperature increases with increase in due to increase in convection. Figure 14 illustrates that the temperature increases with increasing heat source parameter β. It is observed from fig. 15 that the increase in volume fraction of the fluid increases the temperature in the deformable porous channel. The decrease in temperature with increasing drag and suction/injection velocity are presented in figs. 16 and 17 respectively.
The effects of the volume fraction of the fluid ϕ, viscous dissipation parameter or group parameter Br/Ω, Brinkmen numberBr, drag δ and suction velocity V on the entropy generation number Ns are shown in figs. 18-22 respectively. In this study P = –1, ϕ = 0.6, δ = 1, β = 1, Re = 1, Fr = 1, m = 1, U0 = 1, V = 0.5, (Br/Ω) = 1 and Br = 1 are used for numerical computation. From fig. 18, it is noticed that the entropy generation number Ns enhances near the fixed plate y = 0 while it declines near the moving plate y = 1 with increase in volume fraction of the fluid ϕ . Group parameter is an important dimensionless number for entropy generation rate analysis. Figure 19, shows that the entropy generation number increases with growing values of group parameter . The magnitude of entropy generation number takes higher values for larger values of group parameter, since the group parameter determines the relative importance of viscous effects. The increase in entropy generation number with increase in Brinkman number is depicted in fig. 20. The increase in Br enhances fluid frictional heating in the system which enhances the entropy generation. The influence of Br on entropy generation is more pronounced near the channel walls while the effect lessens towards the centerline from both porous plates. It is observed that there exists a fluid section within the channel where the influence of Br on entropy is absent. Figure 21 reveals that the entropy generation number decreases near the stationary plate y = 0 where it increases near the moving plate y = 1 with increasing drag parameter δ . It is seen from fig. 22 that the entropy generation number decreases near the stationary plate y = 0 and increases near the moving plate y = 1 with increase suction/injection velocity V.
Bejan number profiles for various values of Br/Ω and Br are presented in figs. 23 - 24 respectively. Figure 23 reveals that the increase in group parameter decreases Bejan number. This figure displays that for moderately small values of Br/Ω, heat transfer irreversibility dominates entropy generation near the stationary porous wall but the dominance is gradually being ceded to the fluid fraction irreversibility as Br/Ω increases. Near the moving porous wall, fluid fraction irreversibility dominates the entropy generation for all values of the group parameter. It is observed from fig. 24, that the Bejan number increases near the plate y = 0 and y = 1 decreases near the plate with increase in Brinkman number . That is for relatively small values of Br, heat transfer irreversibility dominates entropy generation near the stationary porous wall while fluid friction irreversibility dominates from the centerline to the moving porous wall. For large values of Br, heat transfer irreversibility dominates entropy generation near porous walls. The variation of Nusslet number along with volume fraction of the fluid for different values of drag parameter is shown in fig. 25. It is noticed that the enhancement of drag parameter reduces the Nusselt number.
From Table I, It is found that the skin friction at both the walls increases with increasing the upper plate velocity U0. We also observed that the skin friction at the stationary plate increases with an increase in ϕ whereas it decreases with increasing δ. Further an opposite behavior is noticed at moving plate (upper) of the porous channel.
Effects of ϕ, δ and U0 on skin friction τ0, τ1, for fixed values of P = –8, Re = 1, V = 0.5 and Fr = 1
ϕ
δ
U0
τ0
τ1
0.2
1
1
1.7616
0.2830
0.4
1
1
2.4418
-0.5161
0.6
1
1
3.1219
-1.3151
0.8
1
1
3.8021
-2.1141
0.6
1
1
3.1219
-1.3151
0.6
1.5
1
2.9890
-1.0667
0.6
2
1
2.8678
-0.8345
0.6
2.5
1
2.7568
-0.6167
0.6
1
1
3.1219
-1.3151
0.6
1
2
3.7782
0.2663
0.6
1
3
4.4345
1.8477
0.6
1
4
5.0907
3.4290
From Table II, It is noticed that the increase in the volume fraction of the fluid ϕ or upper plate velocity U0 increases the mass flux in the deformable porous channel. It is observed that the mass flux in the deformable porous channel decreases with increasing drag. Further the mass flux in the porous channel decreases with increasing suction/injection velocity V.
Effect of ϕ, δ, V, and U0 on Mass flux for fixed values of P = –8, Re = 1 and Fr = 1.
ϕ
δ
V
U0
M
0.2
1
0.5
1
0.6214
0.4
1
0.5
1
0.7422
0.6
1
0.5
1
0.8630
0.8
1
0.5
1
0.9837
0.6
1
0.5
1
0.8630
0.6
1.5
0.5
1
0.8295
0.6
2
0.5
1
0.7989
0.6
2.5
0.5
1
0.7706
0.6
1
1
1
0.8223
0.6
1
2
1
0.7388
0.6
1
3
1
0.6581
0.6
1
4
1
0.5849
0.6
1
0.5
1
0.8630
0.6
1
0.5
2
1.2881
0.6
1
0.5
3
1.7132
0.6
1
0.5
4
2.1383
Conclusion
In this article, entropy generation analysis for Couette flow through a deformable porous channel is examined. The expressions for the fluid velocity, solid displacement and temperature distribution are obtained analytically. The effects of volume fraction, drag, upper plate velocity and suction/injection velocity on mass flux are studied through tables. The influence of suction/injection velocity is significant on solid displacement, fluid velocity and temperature distribution which transitively affects the entropy generation with in the deformable porous channel. Some of the important observations are summarized as follows
The increase in upper plate velocity causes more solid displacement in the deformable porous channel. Solid displacement enhances with increasing drag where opposite behavior is observed in the case of increasing volume fraction of the fluid and suction/injection velocity.
It is clear that the fluid velocity increases with increasing ϕ, U0 and P where the fluid velocity reduced due to an increase in the drag and suction/injection velocity.
It is noticed that the temperature increases with increase in U0, Br, m, β and ϕ where as an opposite behavior is observed with increase in δ and V.
It is found that the entropy generation number decreases near the stationary plate y = 0 and increases near the moving plate y = 1 with increase V suction/injection velocity.
It is observed that the mass flux increases with increase in ϕ or U0 where as it decreases with increasing drag or the suction/injection velocity.