In this paper the entropy generation is studied due to micropolar fluid flow through an inclined channel of parallel plates with constant pressure gradient. The lower plate is maintained at constant temperature and upper plate at a constant heat flux. The governing equations are solved by applying the spectral quasilinearization method. The velocity, microrotation and temperature profiles are obtained numerically and are used to calculate the entropy generation number. The influence of pertinent parameters on velocity, microrotation, temperature, entropy generation and Bejan number is discussed with the help of graphs. The results reveal that the entropy generation number increases with the increase in Brinkman number and angle of inclination. Further, it is observed that the increase in coupling number, Prandtl number and Reynolds number reduces the entropy generation number.
Micropolar fluid; Spectral quasilinearization; Entropy; Bejan number
The production of thermal and engineering devices is concerned by irreversible losses that lead to increase the entropy and decrease system efficiency. Thus it is important to find the factors that minimize the entropy generation and maximize the flow system efficiency. To analyze the irreversibilities in the form of entropy generation, the second law of thermodynamics is applied. Factors that are responsible for the irreversibility are heat transfer across finite temperature gradients, characteristics of convective heat transfer and viscous dissipation. Most of the energy related applications such as cooling of modern electronic systems, solar power collectors, and geothermal energy systems depend on entropy generation. Several investigations ([1], [2] and [3]) were carried out on entropy generation under various flow configurations.
Fluid flow and heat transfer inside channels with simple geometry and different boundary conditions is one of the fundamental areas of research in engineering. It has wide range of applications such as thermal insulation engineering, electronics cooling and water movement in geothermal reservoirs. Recently, a wide literature on fluid flow, heat transfer and entropy generation in various channels has been developed. Havzali et al. [4] investigated the entropy generation of incompressible, viscous fluid flow through an inclined channel with isothermal boundary conditions. They observed that the entropy generated in a small section is dominant on the total entropy produced in the entire system. Komurgoz et al. [5] discussed the entropy generation of porousfluid layers contained in an inclined channel by using the differential transform method. Karamallah et al. [6] studied the consequence of differential heat isothermal walls on entropy generation of a vertical square channel saturated with porous media. Liu and Lo [7] performed a numerical analysis on the entropy generation within a mixed convection magneto hydrodynamic (MHD) flow in a parallelplate vertical channel. They observed that the minimum entropy generation number and the maximum Bejan number occur at centerline region of the channel under asymmetric heating conditions. Makinde and Eegunjobi [8] addressed the irreversibility analysis in the flow of couple stress fluid through a porous medium. Damseh et al. [9] studied the presence of transverse magnetic field on entropy generation due to laminar forced convection flow in a channel. Das and Jana [10] investigated the combined effects of Navier slips, suction/injection and magnetic field on entropy generation in a porous channel by considering the constant pressure gradient. Chen and Liu [11] numerically investigated the viscous dissipation effect on entropy generation due to mixed convection nanofluid flow within a vertical channel.
The study of heat transfer has much importance in high temperature processes such as nuclear plants, gas turbines, thermal energy storage. Different kinds of boundary conditions are applied in heat transfer processes. Thermal boundary conditions, heat flux boundary conditions and convective boundary conditions are commonly encountered in heat transfer processes. The effects of these boundary conditions on entropy generation for any fluid flow of different geometries have been studied by several authors. Mahmud and Fraser [12] applied second law of thermodynamics to analyze the heat transfer and fluid motion in rotating concentric cylinders using heat flux boundary conditions. Iman [13] investigated the importance of thermal boundary conditions of the heated/cooled walls in the development of flow and heat transfer, and observed the characteristics of entropy generation in a porous enclosure. Antar and ElShaarawi [14] investigated the effects of uniform heat flux boundary conditions on entropy generation over a rotating solid sphere with forced convection. Butt et al. [15] presented the effects of hydrodynamic slip on entropy generation in a viscous flow over a vertical plate with convective boundary conditions. Anandalakshmi and Basak [16] studied the effect of various heating patterns (different heating and Rayleigh–Benard convection) on entropy generation in a porous rhombic enclosure for different Pr and inclination angles. Anand [17] discussed the velocity slip effect on heat transfer and entropy generation of fully developed power law fluid flow in a micro channel. Mostafa and Ali [18] presented the effect of slip boundary condition on entropy generation for Newtonian and nonNewtonian fluid flows through a parallel plates channel. Ibanez [19] studied the combined effects of magnetic field, suction/injection Reynolds number and hydrodynamic slip on the entropy generation subjected to convective boundary conditions.
The majority of these studies deal with the traditional Newtonian fluids. The study of Newtonian fluids is not sufficient to characterize the flow properties of polymeric fluids, animal blood, coal slurries, mine tailings and mineral suspensions. Such properties are described in nonNewtonian fluid flow model. Many fluids in nature and industrial processes show a nonNewtonian fluid behavior. A number of mathematical models were proposed to explain the rheological behavior of nonNewtonian fluids. One among them is micropolar fluid introduced by Eringen [20] who exhibits certain microscopic effects arising from the local structure and microrotations of fluid elements. The Micropolar fluids accurately simulate the flow characteristics of geomorphological sediments, polymeric additives, colloidal suspensions, liquid crystals, lubricants, hematological suspensions, etc. The rotation of fluid particles in micropolar fluid model is governed by independent kinematic vector called the microrotation vector, which makes it different from other nonNewtonian fluids. Heat transfer in micropolar fluids is also important in the context of industrial manufacturing processes and aerospace engineering and also in chemical engineering.
Although the analysis of entropy generation in a micropolar fluid is important, very little work has been reported in the literature. Ramana Murthy and Srinivas [21] considered the problem of First and Second Law Analysis for the MHD Flow of Two Immiscible Couple Stress Fluids between Two Parallel Plates. Jangili and Murthy [22] analyzed the entropy generation on the steady Poiseuille flow of two immiscible incompressible micropolar fluids between two horizontal parallel plates of a channel with constant wall temperatures in terms of entropy generation. Though most of the work has been done on entropy generation, to the best of the authors’ knowledge, entropy generation analysis for micropolar fluid flow with heat flux boundary condition has not yet been addressed in the literature. Therefore, the objective of the present paper is to investigate the characteristics of micropolar fluid flow, heat transfer and entropy generation inside an inclined channel with constant pressure gradient and constant flux on the upper plate. The governing equations are solved using spectral quasilinearization method. The behavior of flow characteristics with pertinent flow parameters is discussed through graphs.
Consider a steady laminar incompressible micropolar fluid flow bounded by two infinite inclined parallel plates separated by a distance . The channel is inclined at an angle α. Choose the cartesian coordinate system with xaxis aligned at the center of the channel in the direction of the flow and the yaxis perpendicular to the plates. The plates are of infinite length in xdirection so that the physical quantities do not depend on x. Therefore the physical variables are depending on y only. Assume that the properties of the fluid are constant other than the variations of density in the buoyancy force term. Let be the fluid velocity vector. The transpiration cross flow velocity is constant, where is the injection velocity and is the suction velocity. In this study the upper plate of the channel is subjected to a constant uniform heat flux q (isoflux) while the lower plate of the channel is kept at constant temperature (isothermal).
Under these assumptions and Boussinesq approximations the governing equations for the steady flow of micropolar fluid are

(1) 

(2) 

(3) 

(4) 
where u is the velocity component in the x direction, σ is the microrotation component in z direction, T is the temperature, ρ and are the fluid density and gyration parameters, and γ are the material constants (viscosity coefficients), p is pressure, is the acceleration due to gravity, is the coefficient of thermal conductivity, β is the coefficient of thermal expansion.
The boundary conditions are

(5a) 

(5b) 
Introducing the following dimensionless quantities:

(6) 
in Eqs. (1), (2), (3) and (4), we get the following nonlinear system of differential equations:

(7) 

(8) 

(9) 
where primes denote differentiation with respect to is the characteristic velocity, is the Prandtl number, is the Reynolds number, is the suction/injuction Reynolds number, is the coupling number, is the Grashof number, is the constant pressure gradient, is the micropolar parameter, is the microinertial parameter, is the Brinkman number.
The corresponding boundary conditions are as follows:

(10a) 

(10b) 
The system of Eqs. (7), (8) and (9) along with the boundary conditions (10a) and (10b) is solved using the spectral quasilinearization method [23], [24] and [25]. Initially quasilinearization iteration scheme is applied to linearize the system of Eqs. (7), (8) and (9). This quasilinearization method is a generalization of the Newton–Raphson method and was proposed by Bellman and Kalaba [26] for solving nonlinear boundary value problems. In this method the iteration scheme is obtained by linearizing the nonlinear component of a differential equation using the Taylor series expansion. Chebyshev spectral collocation method is then used to solve the resulting scheme.
Let the and be an approximate current solution and and be an improved solution of the system of Eqs. (7), (8) and (9). By taking Taylors series expansion of nonlinear terms in (7), (8) and (9) around the current solution and neglecting the second and higher order derivative terms, we get the linearized equations as follows:

(11) 

(12) 

(13) 
where the coefficients are known functions calculated from previous iterations and are defined as

The above linearized Eqs. (11), (12) and (13) are solved using the Chebyshev spectral collocation method [27]. The approximations are done by using the Chebyshev interpolating polynomials for the unknown functions and further they are collocated at the GaussLobatto points represented as

(14) 
where J is the number of collocation points used.
The functions and are approximated at the collocation points by

(15) 
where is the th Chebyshev polynomial defined by
The derivatives of the variables at the collocation points are represented as

(16) 
where ‘’ is the order of differentiation and being the Chebyshev spectral differentiation matrix. Substituting Eqs. (15) and (16) into Eqs. (11), (12) and (13) leads to the matrix equation

(17) 
In Eq. (17), is a square matrix and and are column vectors defined by

(18) 
where

(19) 
Here represents identity matrix, zero matrix respectively.
The corresponding boundary conditions

(20a) 

(20b) 
The boundary conditions (20a) and (20b) are incorporated in the matrix system (17), and thus the solution is obtained as

(21) 
The initial approximations and are chosen to be functions that satisfy the boundary conditions (10a) and (10b) i.e.

(22) 
The volumetric rate of entropy generation for incompressible micropolar fluid is given as

According to [3], the dimensionless entropy generation number is the ratio of the volumetric entropy generation rate to the characteristic entropy generation rate. Thus the entropy generation number is given by

(23) 
where the characteristic entropy generation rate is and is the dimensionless value and it is a controlling parameter that depends on the heat flux, temperature of the fluid, length of the channel and the thermal conductivity. Eq. (23) can be expressed alternatively as follows:

(24) 
The first term on the right hand side of this equation denotes the entropy generation due to heat transfer irreversibility and the second term represents the entropy generation due to viscous dissipation.
To evaluate the irreversibility distribution, the parameter Be (Bejan number), which is the ratio of entropy generation due to heat transfer to the overall entropy generation (24) is defined as follows

(25) 
The Bejan number varies from 0 to 1. Subsequently, reveals that the irreversibility due to viscous dissipation dominates, whereas indicates that the irreversibility due to heat transfer is dominant. It is obvious that the is the case in which the irreversibility due to heat transfer is equal to viscous dissipation in the entropy production.
In this section we present the results that were found using the spectral quasilinearization method (SQLM). The number of collocation points used is . To check the accuracy of our numerical scheme, we have compared our results with analytical solution of Ariman and Cakmak [28] in the absence of with the fixed values of and . The comparison in the above case is found to be in good agreement as shown in Table 1.
η  Velocity  Microrotation  

Analytical solution [28]  Present  Analytical solution [28]  Present  
−1  0  0  0  0 
−0.8090  0.1557  0.1557  −0.0204  −0.0204 
−0.6129  0.2817  0.2817  −0.0275  −0.0275 
−0.4258  0.3696  0.3696  −0.0248  −0.0248 
−0.2181  0.4302  0.4302  −0.0147  −0.0147 
0  0.4518  0.4518  0  0 
0.2181  0.4302  0.4302  0.0147  0.0147 
0.4258  0.3696  0.3696  0.0248  0.0248 
0.6129  0.2817  0.2817  0.0275  0.0275 
0.8090  0.1557  0.1557  0.0204  0.0204 
1  0  0  0  0 
The solutions for the dimensionless velocity, microrotation, temperature, entropy generation and Bejan number are computed and are shown graphically in Figure 1, Figure 2, Figure 3, Figure 4 and Figure 5. The effects of coupling number N, angle of inclination α, Brinkman number Br, Reynolds number Re and Prandtl number Pr on the nondimensional governing parameters are discussed. To study the effects of and Pr computations were carried out by taking and .

Figure 1. Effect of coupling number on velocity, microrotation, temperature, entropy generation and Bejan number for .


Figure 2. Effect of angle of inclination on velocity, microrotation, temperature, entropy generation and Bejan number for .


Figure 3. Effect of Reynolds number on velocity, microrotation, temperature, entropy generation and Bejan number for .


Figure 4. Effect of Prandtl number on velocity, microrotation, temperature, entropy generation and Bejan number for .


Figure 5. Effect of Brinkman number on velocity, microrotation, temperature, entropy generation and Bejan number for .

The effect of the coupling number on velocity, microrotation, temperature, entropy generation and Bejan number of the fluid flow through an inclined channel is plotted in Fig. 1(a)–(e). The coupling of linear and rotational motion arising from the micromotion of the fluid molecules is characterized by coupling number. Hence, the coupling between the Newtonian and rotational viscosities is represented by N . The microstructure effect is significant as , and for a smaller value of N the substructure individuality is limited. The fluid is nonpolar as its micropolarity is lost at i.e. . Thus, for viscous fluid . Generally a parabolic velocity profile is observed, with maximum velocity along the channel centerline and minimum velocity at the plates.
Fig. 1(a) shows the decreasing nature of velocity with increasing coupling number. The velocity for viscous fluid is more compared to micropolar fluid. It is observed from Fig. 1(b) that the component of microrotation increases near the lower plate and decreases near the upper plate with increase in the coupling number N. It is observed from Fig. 1(c) and (d) that the temperature and entropy generation decrease with increase in the value of coupling number. It is clear from Fig. 1(e) that the Bejan number increase with increase in the value of N. We observe that heat transfer irreversibility dominates around the centerline region of the channel, and the fluid friction dominates at the lower plate.
Fig. 2(a)–(e) present the effect of angle of inclination α on velocity, microrotation, temperature, entropy generation and Bejan number. It is noticed from Fig. 2(a) that the velocity increases with the increase in angle of inclination. It is shown from Fig. 2(b) that the microrotation component decreases near the lower plate and increases near the upper plate with an increase in the value of angle of inclination, showing a reverse rotation near the two boundaries. It is observed from Fig. 2(c) and (d) that the temperature and entropy generation increases with increase in the value of angle of inclination. We observe that the entropygeneration rate is less at the lower plate of the channel and increases quite rapidly to its maximum values at the upper plate of the channel for all the parameter variations. It is observed from Fig. 2(e) that the Bejan number is not much effected with increase in the value of α at the center of the channel but Be decreases at the upper plate and increases at the lower plate.
From Fig. 3(a)–(c), we observe that as the Reynolds number increases the velocity, microrotation(numerically) and temperature decrease. Fig. 3(d) displays the entropy generation as a function of Reynolds number Re. The entropy generation decreases with increase in Re. This may be attributed to the decrease in the velocity which results in the decrease in heat transfer and hence decrease in entropy generation. Fig. 3(e) shows that the Bejan number decreases as Re increases at the lower plate.
It is observed from Fig. 4(a) that the velocity decreases with increase in Prandtl number. It is shown from Fig. 4(b) that the microrotation component increases near the lower plate and decreases near the upper plate with an increase in the value of Prandtl number. Fig. 4(c) represents the Prandtl number effect on temperature. It is observed that the increase in Prandtl number causes a decrease in temperature. Fig. 4(d) shows that the increase in Prandtl number causes a decrease in entropy generation. Since the increase in Prandtl number reduces the fluid temperature as highlighted in Fig. 4(c), the entropy generation decreases. Fig. 4(e) shows that the Bejan number decreases as Prandtl number increases.
The effect of Brinkman number on velocity, microrotation, temperature, entropy generation and Bejan number is displayed in Fig. 5. Fig. 5(a) depicts that the nondimensional velocity increases with an increase in the Brinkman number. It is seen from Fig. 5(b) that the microrotation component decreases near the lower plate and increases near the upper plate with an increase in the value of Brinkman number. It is observed from Fig. 5(c) and (d) that the temperature and entropy generation increase with increase in the value of Brinkman number. Brinkman number is related to heat conduction from a channel wall to a flowing fluid. The effect of viscous forces on entropy generation is significant in the region close to the channel walls. Ns profiles are similar in shape and are almost parallel to one another, for variation of all the parameters but they vary in magnitude. From Fig. 5(e) it is interesting to note that increasing values of Br result in an increase in dominant effect of fluidfriction irreversibility near the upper plate and decrease in dominant effect of fluid friction irreversibility near the lower plate.
In this study, the flow of micropolar fluid between two inclined parallel plates with the constant pressure gradient is studied. The solutions are obtained by applying the spectral quasilinearization method. The entropy generation number Ns at every point η between the channels is found. The effect of various parameters on entropy generation number (Ns) and Bejan number (Be) is studied through the graphs. The following are the observations from the present study:
Published on 12/04/17
Licence: Other
Are you one of the authors of this document?