Heat and mass transfer and thermodynamic analysis of power-law fluid flow in a porous microchannel

Transfer of heat and mass and thermodynamic irreversibilities are investigated in a porous, parallel-plate microreactor in which the working fluid is non-Newtonian. The investigated microreactor features thick flat walls with uneven thicknesses, which can be subject to different thermal loads. The dimensionless governing equations of the resultant asymmetric problem are first derived theoretically and then solved numerically by using a finite volume technique. This results in two-dimensional solutions for the velocity, temperature and concentration fields as well as the distributions of Nusselt number and local and total entropy generations. The results clearly demonstrate the significance of the numerical value of the power-law index and departure from Newtonian behavior of the fluid. In particular, it is shown that by increasing the value of power-law index the Nusselt number on the wall decreases. This leads to the intensification of the temperature gradients in the system and therefore magnifies the local and total entropy generations. Also, it is shown that the wall thickness and thermal asymmetry can majorly affect the heat transfer process and thermodynamic irreversibility of the microreactor. It is noted that the current work is the first comprehensive study of heat transfer and entropy generation in porous micro-chemical reactor with non-Newtonian, power-law fluid.


List of symbols a sf
Interfacial area per unit volume of porous media (m −1 ) Bi Biot number h 1 Half-thickness of the microchannel (m) h 2 Half-height of microchannel (m) h sf Interstitial heat transfer coefficient (W K −1 m −2 ) H w Wall heat transfer coefficient (W K −1 m −2 ) k Effective thermal conductivity ratio of the fluid and the porous solid k 1 Thermal conductivity of wall 1 (W K −1 m −2 ) k 2 Thermal conductivity of wall 2 (W K −1 m −2 ) k e1 Ratio of thermal conductivity of wall 1 and thermal conductivity of the porous solid k e2 Ratio of thermal conductivity of wall 2 and thermal conductivity of the porous solid k ef Effective thermal conductivity of the fluid phase (W K −1 m −2 ) k es Effective thermal conductivity of the solid phase of the porous medium (W K −1 m −2 ) k R Reaction rate constant on the walls (m s −1 ) L Length of the microchannel (m) m Power law consistency index n Power law index N DI Dimensionless diffusive irreversibility N FF Dimensionless fluid friction irreversibility 1 3 N

Introduction
Catalytic porous microchannels, where catalytic channels of small sizes filled with porous media, have received increased interest in traditional and advanced chemical and energy technologies in the last decade. Featuring high surface-tovolume ratios, catalytic porous microchannels can readily achieve the desired performance on transport and chemical reactions [1][2][3]. Therefore, extensive studies have been reported in recent years on applications of catalytic porous microchannels [4][5][6][7]. However, due to the relatively complex natures in catalytic porous microchannels, our understanding of relevant fundamental mechanisms is still far from satisfactory and there are considerable voids for efficient utilization of catalytic porous microchannels. Thus, accurately quantifying the inherent characteristics in catalytic porous microchannels is the prerequisite for their subsequent efficient utilization and optimization. Among the massive aspects that require in-depth investigation for catalytic porous microchannels, transport processes and entropy generation are of vital importance as they are directly related to the overall performance of a catalytic porous microchannel.
Owing to the coupled complexities of multi-physics and geometrical topology in a typical catalytic porous microchannel, a comprehensive investigation of intrinsic characteristics and engineering principles through experiments may be too expensive [8][9][10]. Additionally, accurate measurements at the small spatiotemporal scales in catalytic porous microchannels pose great challenges to current measuring techniques [11,12]. Considering the dramatical increase in computational power and rapid development of theoretical models, numerical modeling has played greater roles in scientific research and engineering optimization for catalytic porous microchannels. Although the general transport [13][14][15] and entropy generation [16][17][18] characteristics of microreactors have been investigated to some extent, their behavior under non-Newtonian fluids is still largely unknown. As a matter of fact, there is still a significant knowledge gap in understating convective heat and mass transfer in porous media and microchannels with non-Newtonian fluid flow [19][20][21]. In the following, some of the recent studies in these areas are briefly discussed.
Neffah et al. [22] conducted a numerical study of heat and mass transfer in a non-Newtonian fluid in a parallel-plate channel partly filled with an anisotropic porous medium. They investigated the effects of various parameters pertinent to the porous medium, the non-Newtonian fluid and the chemical reaction on the velocity, temperature and concentration distributions, as well as mean Nusselt and Sherwood numbers. Neffah et al. [22] showed that compared with the isotropic case, the anisotropy of a porous medium can lead to significant improvements in heat and mass transfer. The shear-thickening fluids exhibit the highest values of mean Nusselt and Sherwood numbers at large Darcy number values. The work of Neffah et al. [22] clarified that an increment in the chemical reaction parameters mitigates heat and mass transfer rates.
In a recent numerical work, Wang et al. [23] investigated entropy generation in thermal natural convection with differentially discrete heat boundary conditions at various Rayleigh numbers. This study showed that the thermal, viscous and total entropy generation increase with the increase in Rayleigh number. It also confirmed that the system's availability improves in the presence of effective heat source at the bottom of the boundary.
Saeed et al. [24] analyzed the methods of biofuel processing in diesel engines to evaluate non-Newtonian biofuel flow in a magnetic microreactor. It was shown that variations in the magnetic field strength and the thermo-physical properties of the fluid could cause considerable temperature differences. It was also observed that the alterations in the properties of Casson fluid as well as the change in the magnitude and angle of the magnetic field could affect Nusselt number. It was further found that the extent of the variations is strongly dependent on the wall thickness.
Gholamalizadeh et al. [25] performed a two-dimensional numerical study of a convective flow in a nanofluid water/ FMWCNT-coated microchannel. The slip velocity boundary condition was applied for the solid walls. The results showed that as compared to Reynolds number, Darcy number has a more pronounced effect upon the velocity profile. It was also found that increasing porosity did not affect the velocity profile growth in any way.
It has been already reported that working fluids have significant impact on the overall performance of a catalytic porous microchannel [1,2,7]. In a numerical study, Maleki et al. [26] investigated flow and heat transfer in non-Newtonian nanofluids over porous surfaces. They investigated the effects of parameters, such as power-law index, volume fraction of nanoparticles, nanoparticles type and permeability parameter on the flow and heat transfer of the desired nanofluid. The study found that the hydro-thermal properties of non-Newtonian nanofluid flows exhibited different behaviors as compared with the common working fluids. For instant, the injections and permeability in the plate showed higher heat transfer performance for non-Newtonian nanofluids in comparison with Newtonian nanofluids.
In their numerical work, Animasaun and Pop [27] considered the geometry and position of the surface for the study of the desired fluid flow in order to evaluate the flow performance on the target parameters. The results showed that the non-Newtonian fluid temperature distribution was higher than that of the Newtonian fluid. It was also found that by applying the non-Newtonian fluid, the rate of local heat transfer decreases more rapidly. The apparent disagreement of this finding with that reported by Maleki et al. [26] is a clear indication that the field of heat convection in non-Newtonian fluids has not settled yet.
Al-Rashed et al. [28] conducted a study to evaluate the thermal and entropic properties of a non-Newtonian nanofluid containing CuO particles in the substrate (MCHS). This article discusses the effects of concentration of nanoparticles, the Reynolds number as well as the geometric size intended for MCHS based on the first and second laws of thermodynamics. The results showed that increasing Reynolds number improves the performance of MCHS by increasing the fluid convection heat transfer coefficient of the working fluid, which in turn results in uniform temperature of the substrate. They further revealed that the use of non-Newtonian nanofluid instead of the base ordinary fluid leads to an increase in heat transfer efficiency against increasing the pressure drop. Also, by using non-Newtonian nanofluid, the least amount of entropy generation in the system can be achieved.
So far, most studies and applications of catalytic porous microchannels have used Newtonian fluids, while effects of non-Newtonian fluids on thermal transport and entropy generation have been seldomly touched [29,30]. However, non-Newtonian fluids are more common in practice and our knowledge of their influences on performance of catalytic porous microchannels should be advanced. To partially close the gap of clarifying the effects of non-Newtonian fluids, a catalytic porous microchannel with power-law fluids has been studied in this study. The reasons to choose powerlaw fluids in this work are twofold. First, the mathematical description of power-law fluids is relatively simple, which can greatly benefit the formulation of conservation equations for modeling of thermal transport and entropy generation in porous microchannels. Second, as a widely encountered type of non-Newtonian fluids, power-law fluids are rather representative in the engineering.
In this work, numerical analysis of thermal transport and entropy generation of a typical catalytic porous microchannel filled with power-law fluids has been conducted. To focus on the effects of power-law fluids, catalytic reactions only happen at the sidewalls of channels and Soret effect is taken into consideration. In the following, governing equations of the studied system are described and discussed. Then, analytical solutions of relevant governing equations are derived and solved numerically. After validating the developed theoretical solution, the differences of thermal transport and entropy generation between Newtonian and power-law fluids are compared and discussed. Finally, a comprehensive parametric analysis was conducted to study the effects of physical properties of power-law fluids.

Assumptions
The assumptions considered in this study are presented as the following.
• The flow is laminar, steady and hydrodynamically fully developed. • The fluid is non-Newtonian power law. • The flow regime is assumed to be continuous with Knudsen number lower than 0.3. • All chemical reactions on the catalyst are assumed to be of zeroth order. • The effect of temperature on the surface reactions is neglected. • The porous medium is isotropic, homogeneous, saturated with fluid and under thermally non-equilibrium conditions (LTNE).

Governing equations and boundary conditions
The rheological equation for an isotropic, incompressible flow of a power-law fluid is [31] where m, u, y and n are the power-law consistency index, fluid velocity, dimensional transverse coordinate and powerlaw behavior index, respectively. The Darcy-Brinkman equation for the momentum equation [2,32,33] is considered as the follows: The following equations represent the energy transfer in the solid and fluid phases: The mass transfer involves the effect of thermal diffusion (Sorter effects) on the Fickian diffusion in the advective-diffusive system, as follows [1]: In Eqs. (2)-(4), p is thepressure, x is the dimensional axial coordinate, and eff are the viscosity of the fluid and the effective viscosity of the porous medium, respectively, is the permeability of the porous insert, h 1 and h 2, refer to the half-thickness and half-height of the microchannel. Further, k ef and k es are the effective thermal conductivity of the fluid phase and solid phase of the porous medium, k 1 and k 2 are the thermal conductivity of wall 1 and wall 2 (see Fig. 1), h sf and a sf are the interstitial heat transfer coefficient and interfacial area per unit volume of porous media. Furthermore, T s , T f , T 1 and T 2 denote the temperature of solid phase, fluid phase, lower and upper wall, respectively. Also, f and C p,f show the density and specific heat capacity of the fluid, C is the mass species concentration, D denotes the effective mass diffusion coefficient and D T represents the coefficient of thermal mass diffusion. Although most solutes feature positive coefficients of thermal mass diffusion representing that the mass can be moved down the thermal gradient, the other have the negative coefficients and are able to be diffused up the thermal gradient. Therefore, the negative or positive sign of the last term in Eq. (4) is acceptable.
For energy and momentum equations, the boundary conditions are The concentration boundary conditions are expressed as follows: Next, the dimensionless parameters for better physical analysis are stated in Table 1.
The dimensionless temperature follows the given chain, i = 1, 2,s, f.
In Table 1, q ′′ 1 and q ′′ 1 are the lower and upper wall heat flux, is the dimensionless temperature, k R is the reaction rate constant on the walls, X and Y are the dimensionless axial and transverse coordinates, Da, Bi , , Sr, Pr, Re and Pe are the Darcy, Biot, Damköhler, Soret, Prandtl, Reynolds and Peclet numbers, respectively. Also, Q represents the wall heat flux ratio, S is the shape factor of the porous medium, shows the aspect ratio of the microchannel, is the porosity of the porous medium, and U denote the dimensionless concentration and fluid velocity, respectively. Further, k e1 is the ratio of thermal conductivity of wall 1 and thermal conductivity of the porous solid, k e2 is the ratio of thermal conductivity of wall 2 and thermal conductivity of the porous solid and k refers to the effective thermal conductivity ratio of the fluid and the porous solid.

Analytical derivation of the dimensionless governing equations
Given the dimensionless equations defined in Table 1, Eq. (2) can be expressed as Similarly, no-slip boundary conditions and axial symmetry conditions at Y = 0 can be expressed as follows: Solution of Eq. (7) can be expressed with Using Eq. (8), the average dimensionless velocity across the channel equals By combining Eqs. (9) and (10), the following ratio of absolute and mean is defined in Eq. (11): Assuming the flow to be fully developed, as given the assumptions in Section "Assumptions", the conditions stated are [1,32] Subsequently, the temperature equation is [1] In Eq. (13), the function f i (y) is obtained by solving Eq. (3a) with the resultant equation from combination of Eq. (3b) and (3c). By placing the heat flux boundary conditions in Eq. (5b) and (5d), we have Substituting the velocity from Eq. (11) into Eq. (14) and using the dimensionless parameters defined in Table 1, the following equation is taken: where the average non-Newtonian fluid temperature is given as follows: Using the dimensionless parameters of Table 1 and rearranging Eqs. (11) to (15), the dimensionless form of Eq. (3a) to (3d) is obtained.
By applying Eq. (19a) to (19h), it is possible to solve analytically the dimensionless cross-sectional temperature profiles using Eqs. (17a), (17b) and (18a), as the following: where alpha equals By placing dimensionless cross-sectional temperature profiles of Eq. (20) inside Eq. (13), the final dimensionless temperature profiles are expressed by The following equation is developed by applying Leibniz's rule [34] on Eq. (4): The velocity profile along the microchannel is symmetrical, and therefore half of the channel is considered. The average velocity in half of the microchannel is given by The average flow velocity is defined as By setting C 0 as the initial concentration, the following equation is revealed.
Using the dimensionless parameters defined in Table 1, the concentration distribution yields Replacing Eq. (26) into Eq. (4), the distribution of dimensionless concentration function is given as Applying the dimensionless parameters of Table 1 conditions of Eq. (6a) and (6b) is required to solve Eq. (29). These are as follows: By applying parts (a) and (b) of Eq. (30), it is possible to determine analytically the second derivative of dimensionless concentration profiles on the cross section, which is Finally, substituting Eq. (28) into Eq. (26), the dimensionless concentration equation, containing two transverse and axial cross sections, can be expressed as follows:

Nusselt number
The heat convection coefficient on the upper and lower walls of the channel can be expressed by Considering dimensionless parameters presented in Table 1, the Nusselt number is taken by The non-Newtonian fluid average temperature is obtained by the dimensionless form of Eq. (16) along the channel. That is

Entropy generation
The entropy generation in the current system raises from various irreversible sources. The contribution of heat and mass transfer and friction between the fluid and solid surface to exergy destruction cause lower energy quality and weaker hydrothermal performance of the current system. The potential sources of entropy generation in the typical microchannel are presented as the following: The dimensionless results of Eq. (37a-37f) are in which Ṡ ′′′ w1 and Ṡ ′′′ w2 are the volumetric entropy generation rate from the lower and upper walls. Further, Ṡ ′′′ s , Ṡ ′′′ f and Ṡ ′′′ DI are the volumetric entropy generation in the porous solid, fluid and due to mass transfer, respectively. Further, R denotes the specific gas constant, N refers to the dimensionless irreversibility, is the dimensionless heat flux and is the irreversibility distribution ratio.
In order to provide a simple substrate for investigating the entropy generation and evaluating its distribution for different irreversible sources, the equations for the solid and fluid phases in the porous medium are described as follows: By adding Eqs. (41b) and (42c), the volumetric entropy generation expressions are defined as follows: The volumetric entropy generation in the porous media is equal to the sum of the equations applicable to the porous media.
The total entropy generation defines the entropy produced by the heat transfer, viscous dissipative and concentration gradient for each point in the X and Y directions. That is The total entropy generation on the whole microchannel is obtained by sweeping two-segmental integration on the physical domain, which equals

Numerical simulation and validation
The equations derived analytically in Section "Theoretical and numerical methods" provide the general solution of the dimensionless governing equations. Although it is possible to derive fully analytical solutions for the current problem, the non-Newtonian nature of the considered fluid makes such solution immensely complicated and therefore of little use. Hence, to evaluate the unknown coefficients and find the particular solutions, a numerical approach was taken. The details of numerical methodology along with the specification of the system configuration are given in this section. Figure 1 shows a schematic of the geometry considered in the current study. The microchannel is fully filled with porous materials, with two thick walls, subjected to constant but non-equal heat fluxes with a catalytic layer on the inner microchannel surface. A non-Newtonian, power-law fluid flow enters the left side of the microchannel and, after passing through the porous media, leaves from the right side of the system. The problem outlined in Section "Theoretical and numerical methods" was solved using various non-dimensional and dimensional parameters that are presented for quantitative analysis by the default values in Tables 2 and  3, respectively. It is noted that the values of parameters in Table 2 were chosen in a way that the results become comparable with the existing data in the literature for Newtonian fluids [1,3].
The present numerical data are compared with the closed-form analytical results of Hunt et al. [1] for Nusselt number and entropy generation (Figs. 3 and 4) in the limit case that power-law fluid approaches a Newtonian fluid at different values of wall thickness (Y 1 ). Evidently, the two data sets exhibit very good agreement and thus the analytical solutions are deemed to be valid.

Results and discussion
In this section, an emphasis is put on the influences of non-Newtonian fluid characteristics upon the thermal behavior of the micro-system. In the interest of brevity, other effects are not discussed. deduced from the fact that the fluid temperature takes larger values of dimensionless temperature at higher values of power index. Since the dimensionless temperature expresses the difference between the temperature of fluid and that of the wall, higher values of dimensionless temperature imply larger difference between fluid and wall temperature. For the current iso-flux problem, this means weakened convective heat transfer. This argument will be later supported by calculation of Nusselt number. Figures 6 and 7 show the contours related to the local entropy generation inside the fluid phase of the microchannel. It is important to note that entropy generation is only calculated for the fluid phase and porous media and the solid wall are excluded from the analyses of local entropy generation. Also, as depicted by Eq. (41d) the mass transfer irreversibilities are included in the current analysis. Figure 6 shows the total entropy variations for the two power-law indices with two different Damköhler numbers for the fluid phase. It is clarified that the value of Damköhler number influences the rate of mass transfer and thus can affect the irreversibility. As can be seen in this figure, for a value of Damköhler number less than one, an increment in the value of the power-law index intensifies the generation of entropy. This is due to the enhancement of friction irreversibilities at higher values of power-law index. More importantly, as shown earlier, the rate of heat transfer decreases at higher values of power-law index, which in turn strengthens the temperature gradients and thus increases the generation of entropy. Interestingly, this behavior is inversed at higher value of Damköhler number. This is because of the dominance of mass transfer in entropy generation in the current problem [1,35,36], which can overrule other effects. Figure 7 illustrates the influences of changes in Soret number upon the local entropy generation in the fluid phase. This has been done for two different values of the power-law index corresponding to Newtonian and non-Newtonian fluids. Figure 7 indicates that for a constant power-law index, in either of Newtonian and non-Newtonian fluids, Soret number growth causes a reduction in the local irreversibility. This can be explained by noting that in the current problem irreversibility is dominated by that of mass transfer [34][35][36]. Magnifying the value of Soret number enhances the process of mass transfer and hence reduces the irreversibility associated with the transport of mass. Thus, the local entropy generation drops. Further, Fig. 7 shows that increasing the value of power-law index at a fixed Soret number leads to a slight intensification of the local entropy generation. As already discussed, this is mostly due to the magnification of thermal irreversibilities following the fluid temperature gradient strengthening.

Detailed comparison of transport phenomena and entropy generation between Newtonian and power-law fluids
According to Fig. 8, increasing the power-law index results in the reduction in Nusselt number. This reduction can be quite significant such that by increasing the power-law index from 0.4 to 1.8, a drop of around 50% is observed in Nusselt number for all the investigated cases. This is an important result and shows that deviation from the state of Newtonian fluid can majorly affect the rate of heat transfer. Hence, close attention should be paid to  Figure 8 further shows the pronounced effects of the wall thickness upon Nusselt number. In keeping with the results reported in the literature [33,34], it is shown that thicker walls act as thermal resistance and thus tend to reduce the Nusselt number. This figure clearly shows that alteration of wall thickness by around 30% can change the value of Nusselt number by more than 40%. Thus, wall thickness is an essential parameter dominating the rate of heat transfer in microchannels. Figure 9 illustrates the effects of power-law index on the total generation of entropy. Clearly, variations in the power-law index strongly affect the total entropy generation. According to this figure, an increase in the power-law index from 0.4 to 1.8 can boost the total entropy by more than 50%. As already discussed, decreasing Y 1 (thickening the walls) also contributes to this increase. Once again, impediment of heat convection process by increases in the value of power-law index (as already discussed) is responsible for this trend.

Parametric studies on the effects of physical properties of power-law fluids
In Fig. 10, temperature contours for non-Newtonian fluids for two different values of Y 1 are depicted. The righthand side contours correspond to the solid phase, and the left-hand side contours correspond to the fluid phase. It is observed that by increasingY 1 , that is by making the wall thinner, the deflection of temperature contours increases. Thus, at any cross section in the microchannel, the transversal temperature variation increases. As shown previously [1,34,35] and also further supported in the later sections, this causes an increase in the overall heat transfer rate of the system. Figure 11 shows the distribution of dimensionless fluid temperature inside the porous microchannel for different values of Reynolds number and the power-law index. This figure essentially compares the case of Newtonian fluid to that of non-Newtonian fluid. Evidently, regardless of the value of Reynolds number, the axial temperature increase along the microchannel for the non-Newtonian fluid is larger than that of the Newtonian fluid. For the current constant heat flux problem, this implies smaller convection coefficient for the non-Newtonian fluid flow. This finding will be later verified by calculating the Nusselt number for the two types of fluid. Further, the increase in Nusselt number has resulted in a general drop of the dimensionless temperature at any given point inside the microchannel. This is to be expected, as in general, increases in Reynolds number (on the basis of the microchannel height) lead to improved heat convection coefficient and thus result in lower temperatures. Figure 12a shows the effect of variation in the wall heat flux on the Nusselt number. As can be seen by increasing the value of wall heat flux ratio, Q, performed by strengthening the thermal asymmetry of the problem, the Nusselt number increases. This increment in Nusselt number is a consequence of changes in the shape of temperature contours and resultant modification of temperature gradient on the surface of the wall [33][34][35][36]. The effects of Darcy number on the Nusselt number are demonstrated in Fig. 12b.
As can be seen, the smaller the Darcy number, meaning the lower the permeability, the rate of heat transfer is increased. This is a well-known behavior of the Nusselt number [1,5,11,15] and shows the consistency of the current Non-Newtonian simulations with those conducted earlier on Newtonian fluids. In Fig. 12c, the effect of variations in the conductivity coefficient of the solid part of the porous media ( k s ) is presented. This graph shows that with augmentation of k s the value of Nusselt number increases. This can be attributed to the increase in heat transfer between the fluid and the solid part of the porous medium, which is enhanced by boosting the conductivity of the porous solid [37,38]. Figure 13 shows the influences of power-law index upon the values of Nusselt number and total entropy generation against the Reynolds number. As a general trend, both Nusselt number and total entropy generation increase monotonically with Reynolds number. It is to be expected that Nusselt number in forced convection is often strongly correlated with flow Reynolds number. However, the growth of total entropy generation with Reynolds number may sound counterintuitive. This is because Nusselt number rising relaxes the internal temperature gradients inside the microchannel, and this results in the thermal irreversibility reduction. These apparently contradicting observations can be justified by noting that the entropy generation in the current problem is dominated by flow friction and not heat transfer. Switching to non-Newtonian fluid intensifies the flow friction and thus magnifies the total generation of entropy [39]. Evidently, for all Reynolds numbers, the application of non-Newtonian fluid (n = 1.3) leads to a smaller value of Nusselt number (see Fig. 13a). The reduction in heat transfer and the resultant increase in the internal temperature gradient, as well as increases in the frictional losses of the flow due to non-Newtonian effects, are the reasons for augmentation of entropy generation in comparison with that of Newtonian fluid [40,41].

Conclusions
Transport of heat and mass and thermodynamic irreversibilities in a parallel-plate microreactor filled by a homogenous porous medium and with a non-Newtonian working fluid were investigated. The considered configuration was assumed to be geometrically and thermally asymmetric. The power-law behavior was assumed for the fluid, and the porous medium was set under local thermal non-equilibrium. A set of governing equations were first derived theoretically. The resultant nonlinear system of partial differential equations was then solved numerically through using a finite volume solver. This led to the development of two-dimensional solutions for the temperature concentration and local entropy generation fields and evaluation of Nusselt number and total entropy generation. The numerical results were validated against the existing fully analytical solutions for the Newtonian fluids. It was shown that increases in the value of power-law index reduce the rate of heat transfer and thus drop the value of Nusselt number. Hence, the local and total generations of entropy are both strongly affected by the value of power-law index in which increases in this parameter highly enhance the local and total irreversibilities. The other qualitative characteristics of the investigated system were found to remain consistent with those reported for the Newtonian fluids.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.