Analysis of transport from cylindrical surfaces subject to catalytic reactions and non-uniform impinging flows in porous media

This paper investigates forced convection of heat and mass from the catalytic surface of a cylinder featuring non-uniform transpiration and impinging flows in porous media. The non-equilibrium thermodynamics including Soret and Dufour effects and local thermal non-equilibrium are considered. Through employing appropriate change of variables, the governing equations in cylindrical coordinate are reduced to nonlinear ordinary differential equations and solved using a finite difference scheme. This results in the calculation of the temperature and concentration fields as well as the local and surface-averaged Nusselt and Sherwood numbers. The conducted analyses further include evaluation of the rate of entropy generation within the porous medium. It is shown that internal heat exchanges inside the porous medium, represented by Biot number, dominate the temperature fields and Nusselt number. This indicates that consideration of local thermal non-equilibrium is of highly important. It is also demonstrated that Dufour and Soret effects can significantly influence the development of thermal and concentration boundary layers and hence modify the values of Nusselt and Sherwood numbers. In particular, it is shown that small variations in Soret and Dufour numbers can lead to noticeable changes in the average Nusselt and Sherwood numbers. Such modifications are strongly dependent upon the type of transpiration and characteristics of the impinging flow. The present work is the first analysis of non-equilibrium effects upon transport by stagnation flows around the curved surfaces embedded in porous media.

Velocity components along (r -uz)-axis U 0 (u) Transpiration z Axial coordinate Greek symbols a Thermal diffusivity c Modified conductivity ratio h(g, u) Non-dimensional temperature k Permeability parameter,k ¼ a 2 4k 1

Introduction
Non-equilibrium thermodynamics are often of significance in transport of heat and mass in chemically reactive systems [1,2]. Most specifically, the thermal diffusion of mass (Soret effect) and the transport of energy through diffusion of chemical species (Dufour effect) become noticeable in the presence of strong thermal and concentration gradients [2]. Such situations are frequently encountered in chemical reactors. In porous media, a common manifestation of nonequilibrium thermodynamics is through local thermodynamic non-equilibria, which include local thermal nonequilibrium [3,4]. Thus, porous thermochemical systems can involve Soret and Dufour effects and also feature local thermal non-equilibrium. The former has been already the topic of a large number of investigations, e.g. [5][6][7], and the latter has been confirmed in a few recent studies [8][9][10]. Nevertheless, comprehensive non-equilibrium analyses of chemically reactive porous media that take into account all the preceding effects are rare. Resolving this issue is the primary objective of the present work. Surface reactions in general and catalytic reactions in particular are essential in a wide range of chemical reactors. In such reactors, the catalyst might be deposited on the surfaces embedded in a porous medium. This is, perhaps, the reason for the existence of a series of studies on heat and mass transfer from a flat surface covered with porous material. A concise survey of these studies is put forward in the followings. Postelnicu [11] conducted an analysis of natural convection of heat and mass from the surface of a vertical flat plate embedded in saturated porous media. Soret and Dufour effects were considered in this work, and two-dimensional governing equations were solved numerically through employing a similarity solution [11]. Constant wall temperature and concentration were implemented as boundary conditions, and a magnetic field was also applied [11]. It was argued that by intensifying the magnetic field, the local Nusselt and Sherwood numbers would increase [11]. This work was then extended to the cases with homogeneous chemical reactions with varying reaction order and Darcian flow in porous media, [12]. Later, Postelnicu [13] investigated the same problem under a stagnation-point flow and added transpiration of mass to the flat plate. The effects of variable viscosity on the natural convection of heat and mass from a vertical, flat plate covered by porous materials were investigated by Afify [5]. Similar to the work of Postelnicu [11,12], the influences of a magnetic field as well as the effects of Soret and Dufour numbers were taken into account [5]. However, Afify considered a non-Darcian flow with temperature-dependent viscosity [5]. In another work, the effects of viscoelastic fluid were also included in the problem by Hayat et al. [14].
Chamkha and Ben-Nakhi [6] took two major steps in the analysis of magnetohydrodynamic (MHD) convection from a flat plate embedded in porous media with Soret and Dufour effects. These authors [6] introduced the added influences of forced convection through examining mixed convection and also included thermal radiation by employing Rosseland approximation. They reported that the local Nusselt and Sherwood numbers decreased in the presence of a magnetic field for the free and mixed convection regimes. However, Nusselt and Sherwood numbers decrease and then increased forming minima as the mixed convection parameter was increased from free convection regime to the forced convection regime [6]. The existence of thermal radiation was reported to reduce the local Nusselt number and to increase the local Sherwood number particularly when there is fluid suction through the flat surface [6]. A similar configuration was analysed by Tsai and Huang [15], who extended the classical Hiemenz flow through porous media to heat and mass transferring cases with Soret and Dufour effect. These authors also considered the effects of heat of reaction and thermal radiation and non-uniform wall temperature and concentration [15]. It was concluded that for some mixtures (e.g. hydrogen-air) with the low and medium molecular weights, the Soret and Dufour effects are significant and hence should be taken into consideration [15]. Mohammad Hemmat et al. [16,17] considered the problem of mixed convection inside liddriven cavities filled with nanofluids. They reported a comprehensive dataset on fluid flow and heat transfer in different configurations including square, rectangular, triangular and trapezoidal with varying governing parameters such as Rayleigh number, Hartmann number, Richardson number and solid volume fraction [16,17]. Other related works can be found in the literature, see, for example, [18][19][20][21] and references therein.
In the work of Prasad et al. [22], Soret and Dufour effects were included in the magnetohydrodynamic free convection analysis of a vertical plate embedded in porous media. Nonlinear and viscous effects were considered in the momentum transfer, while the wall temperature and concentration were assumed constant [22]. It was concluded that increasing Soret number and decreasing Dufour number intensify the local rate of heat transfer (local Nusselt number) on the surface of the plate with the opposite effect upon mass transfer rate (i.e. local Sherwood number) [22]. More recently, Mabood et al. [23] considered the problem of MHD forced convection of heat and mass transfer in a chemically reactive flow impinging on a flat plate. This analysis did not include Soret and Dufour effects but considered transpiration of mass from the surface of flat plate [23]. Amongst other findings, Mabood et al. demonstrated the significance of mass transpiration upon the transport of heat and mass [23]. In a recent work of Reddy and Chamkha [24], Soret and Dufour effects as well as thermal radiation were considered in the problem of MHD forced convection on a flat plate covered with a porous medium [24]. This investigation also included unsteady and temperature-dependent heat generation and considered two types of nanofluids [24]. It is worth mentioning that the general literature on impinging flows upon stretching surfaces with chemical reactions and Soret and Dufour effect is rather wide (see for example [23,[25][26][27]). However, as the focus of the current work is on the curved surfaces embedded in porous media, nonporous configurations are not discussed here.
A few points are well reflected by the preceding survey of the literature. First, the problem of combined heat and mass transfer by natural convection and by the inclusion of Soret and Dufour effects is well studied [6,15,[22][23][24]. However, the same under forced convection has received less attention. Second, MHD effects have been investigated extensively, e.g. [22][23][24], and influences of thermal radiation and mass transpiration have been also considered [6]. Third, almost all investigated cases have been concerned with flows impinging on a flat plate and curved surfaces have remained unexplored. Importantly, few recent works on curved surfaces in porous media did not consider stagnation flows [28,29]. Fourth, all the existing works have assumed that local thermal equilibrium (LTE) holds within the porous medium and local thermal non-equilibrium (LTNE) effects have been totally ignored. The current work aims at releasing the restrictive assumptions stated in the last two points. In practice, the surfaces in contact with porous media can be curved. This introduces a new level of complexity to the analysis of transport phenomena, which is still largely unexplored. Further, recent works in other configurations showed that forced convection in porous media with chemically reactive flows is most likely to violate the assumption of LTE [9,30,31]. Consequently, analysis of these flows demands an LTNE approach. Given this, the proceeding part of this work develops an LTNE model of forced convection of heat and mass over a cylinder with non-uniform suction/injection of fluid. This is then complemented by an evaluation of entropy generation. It is noted that while there are few studies on entropy generation in stagnation flows on stretching surfaces, e.g. [32]. The equivalent problem in porous media that includes Soret and Dufour effects does not currently exist.

Theoretical and numerical methods
Problem configuration, assumptions and governing equations Figure 1 depicts schematically the problem under investigation. This includes a cylinder with radius a centred at r = 0 covered with a layer of catalyst and embedded in a porous medium. The surface of the cylinder can include uniform or non-uniform transpiration with prescribed circumferential distributions, while the temperature of the external surface of the cylinder is maintained constant. An external axisymmetric radial stagnation-point flow of strain rate of k impinges on the cylinder. Due to the non-uniformity of transpiration, the flow configuration around the cylinder can be highly un-axisymmetric. The following assumptions are made through this work.
• The flow is steady, incompressible and laminar.
• The cylinder is assumed to be infinitely long.
• A zeroth-order, temperature-independent, catalytic reaction [33][34][35] takes place on the external surface of the cylinder. • Thermal diffusion of chemical species (Soret effect) and transport of energy through mass diffusion (Dufour effect) are considered.
• The porous medium is homogenous, isotropic and under local thermal non-equilibrium (LTNE). • The radiation heat transfer, gravitational effects and viscous dissipation of the kinetic energy of the flow are ignored. • Physical properties such as porosity, specific heat, density and thermal conductivity are assumed to be constant, and hence, the thermal dispersion effects are negligible. • A moderate range of pore-scale Reynolds number is considered in the porous medium, and therefore, the nonlinear effects in momentum transfer are negligibly small.
It should be noted that a zeroth-order surface chemical reaction is a fair representation of many catalytic reactions [36] and is therefore of practical significance. The governing equations and boundary conditions, in the cylindrical coordinate system shown in Fig. 1, can be summarised as follows.
The continuity of mass reads, The transport of momentum in the radial direction is while in the axial direction it takes the form of

Porous media
Surface reaction 1 The transport of thermal energy in the fluid phase is expressed by: The last term on the left-hand side of this equation represents the Dufour effect [33,34].
Further, the transport of thermal energy in the solid phase of the porous medium can be written as Mass transfer of chemical species is governed by the following advective-diffusive model, which considers the thermal diffusion of mass (Soret effect) in addition to the classical Fickian diffusion of species [14,32]: In Eqs. (4-6), the subscripts ''f'' and ''s'' refer to the fluid and solid properties, respectively. The velocity boundary conditions of the momentum equations are as follows: Also, the two boundary conditions with respect to u (angular coordinate) are expressed by Equation (7) denotes the no-slip conditions on the external surface of the cylinder. Further, Eq. (8) indicates that the viscous flow solution approaches, in a manner analogous to the Hiemenz flow, the potential flow solution as r ! 1 [37][38][39]. This can be verified by starting from the continuity equation in the followings. À 1 kz and integrating in r and z directions with boundary conditions, w = 0 when z = 0 and u = -U 0 (u) when r = a.
The boundary condition for the transport of thermal energy is given by and the two boundary conditions with respect to the angular coordinate, u, are in which T w is the cylinder surface temperature and T ! is the freestream temperature. The boundary condition for the transport of mass is given by in which D is the molecular diffusion coefficient and k R is the kinetic constant for a zeroth-order chemical reaction [28,29]. Further, C ! is the freestream concentration. The two boundary conditions with respect to angular coordinate, u, are

Self-similar solutions
A reduction of the governing Eqs. (1-6) is developed through applying the following similarity transformations.
where g ¼ r a À Á 2 is the dimensionless radial variable. It is important to note that it has been already demonstrated that Analysis of transport from cylindrical surfaces subject to catalytic reactions and non-uniform… Eq.
(3) has a negligible contribution with the flow field [40]. Hence, this equation is discarded in the similarity solution. Transformations (14) satisfy the two-dimensional version of Eq. (1) automatically, and their substitution into Eqs. (2) and (3) leads to the following system of coupled differential equations.

Nusselt and Sherwood numbers
For the current problem with isothermal boundaries, the local heat convection coefficient and rate of heat transfer for fluid phase are defined as and Hence, Nusselt number on the surface of the cylinder can be written as Similarly, the local mass transfer coefficient and rate of mass transfer are defined as and Hence, Sherwood number can be expressed as Entropy generation Considering the assumption stated in Sect. 3.1, the volumetric rate of local entropy generation in the problem is given by [41,43]: Using the similarly variables given in Eqs. (14) and (37), the local entropy generation reduces to in which are the non-dimensional entropy generations due to heat transfer, fluid friction, mass transfer and the characteristic entropy generation rate, respectively. The dimensionless form of the volumetric rate of local entropy generation (N GT , N GF , N GD ) can be written as follows.
is the diffusive constant parameter, and Þ is the Brinkman number. The Bejan number, Analysis of transport from cylindrical surfaces subject to catalytic reactions and non-uniform… defined as the ratio of entropy generation due to heat transfer to the total entropy generation, is used to facilitate understanding of the mechanisms of entropy generation.
Bejan number for the current problem is expressed as

Grid independency and validation
To ensure about the grid independency of the numerical simulations, the surface-averaged values of Nusselt, Sherwood and Bejan numbers were calculated for mesh sizes of 51 9 18, 102 9 36, 204 9 72, 408 9 144 and 819 9 288. Table 1 shows that there are no considerable changes in these quantities for the (g, u) mesh sizes of (204 9 72), (408 9 144) and (8169 9 288). Thus, a (408 9 144) grid in gu directions was used for the computational domain reported in the present work. A non-uniform grid was applied in g-direction to capture the sharp gradients around the external surface of the cylinder, and a uniform mesh was implemented in u direction. The computational domain extends over u max = 360°and g max = 15. In this expression, g max corresponds to g ? !, which for all investigated cases is located outside the momentum, thermal and concentration boundary layers. As the convergence criterion in the numerical simulations, when the difference between the two consecutive iterations became less than 10 -7 , the iterative process was terminated. On the basis of the implemented numerical scheme, the numerical error estimated to be of O(Dg) 2 [41,42]. Tables 2 and 3 show that in the limit of very large porosity and permeability (no porous material) and in the absence of mass transfer, the numerical solutions developed in Sect. 2 reproduce the results of Wang [44] and Gorla [45] for stagnation flow over a cylinder. Further, although not shown here, it was confirmed that in the limit of large Biot numbers the current LTNE results reduce to the LTE results reported in Ref. [40].

Temperature, concentration and entropy generation fields
In the current problem, the flow field is the result of interactions between the stagnation flow and the transpiration from the surface of the cylinder. Since the fluid was assumed to have constant density, the heat and mass transfer do not influence the flow and hence the hydrodynamics discussed in Ref. [40] remain unaltered. Nonetheless, because of the importance of the flow field it is briefly discussed here. Figure 2 shows the dimensionless radial velocity field (f) when the non-uniform transpiration, shown in Fig. 1, is in place. It is clear from this figure that in the regions of 0°[ u [ 135°there exists a low velocity region, whereas in the rest of the cylinder circumference there is a strong flow towards the centre. Transpiration of mass in the form of blowing and it is interaction with the impinging flow has resulted in the radially stagnant flow in the first quarter of the cylinder cross section. In other parts, however, suction of mass is in the direction of the impinging flow, which results in a uniform flow towards the centre of the sphere. Figure 2 further indicates that the extent of low-speed region is dependent upon the Reynolds number of the external flow and in general increases as the Reynolds number becomes smaller. It will be later shown that the state of the flow around the cylinder greatly influences the transport.  Fig. 3. It is clear from this figure that for low values of Biot number there are significant differences between the temperature distributions in the fluid and solid phases. Nonetheless, as the numerical value of Biot number increases, these differences diminish and for high Biot number (i.e. Bi = 100) the two temperature distributions are very similar. The influences of non-uniform transpiration are completely noticeable in Fig. 3. At low values of Biot number, there is a nearly uniform circumferential temperate distribution in the solid phase. Yet, this is clearly not the case in the fluid phase, in which there is a strong heat transfer in the region of the circumference with fluid blowing (see Fig. 1). In other regions, however, the thickness of thermal boundary layer is generally low. As Biot number increases, the temperature distribution within the solid phase becomes increasingly more non-uniform and variations become limited to small area near the surface of the cylinder. This is to be expected as the strong heat exchanges between the solid and fluid phases at high values of Biot number hinder diffusion of heat away from the hot surface of the cylinder.
In the limit of high Biot number, the extracted heat from the solid phase causes an increase in the thickness of the thermal boundary layer. This can be readily verified in Fig. 3a through comparing the high-temperature region at the back of the cylinder (90°\ u \ 360°) for different values of Biot number. Figure 4 shows the effects of Dufour number on the temperature field within the porous medium. Increasing Dufour number strengthens the thermal energy received by the fluid phase due to the diffusion of mass. It is therefore not surprising that at high Dufour numbers the region of hot fluid within 0 \ u \ 90°has been extended. Interestingly, however, this has little effects upon the temperature distribution in the solid phase. Reynolds number appears to have a pronounced effect on the temperature field (see Fig. 5). At low values of Reynolds number, representing a gentle impinging flow, Fig. 5a shows that the fluid temperature, in almost the entire investigated domain, has been influenced. However, as the flow impingement becomes stronger, two noticeable events happen. First, thickness of the thermal boundary layer in the region of 90°\ u \ 360°decreases. This is due to the well-known effect of flow velocity on the boundary layer thickness, which in this case is further intensified by the action of mass suction. Second, in the region of 0°\ u \ 90°the heat transfer is significantly enhanced, and thus, a larger volume of the fluid phase in this part of the domain   Analysis of transport from cylindrical surfaces subject to catalytic reactions and non-uniform… becomes hot. This is caused by the interactions between the high momentum impinging flow and the blowing of the fluid from the surface of the cylinder. These two opposing flows tend to generate thick hydrodynamic and thermal boundary layers. Thus, a large volume of fluid can almost reach thermal equilibrium with the surface of the cylinder. As indicated in Fig. 5b, the effects of thermal equilibrium are much less significant on the solid phase. This results in distinctively different solid and fluid temperature distributions at high Reynolds numbers. Figure 6 illustrates the influences of Damköhler and Schmidt numbers upon the concentration field. Figure 6a clearly shows that increasing Damköhler number results in significant intensification of the concentration field. This is to be expected as the chemical kinetics of the surface are stronger at higher Damköhler numbers. At low values of Schmidt number, which can be viewed as large mass diffusivity, the concentration filed is fairly uniform. This is due to the fact that at this limit the mass transfer is diffusion dominated, and therefore, the hydrodynamic non- uniformities have rather insignificant effects. However, as the Schmidt number grows in value and thus the mass diffusivity becomes relatively smaller than the kinematic viscosity, the concentration field becomes increasingly more affected by the flow field. Figure 6a and b demonstrates the effects of non-uniform transpiration upon the concentration field. Most noticeably, thickening of the concentration boundary layer for the part of the circumference with blowing transpiration is evident in these figures.
The spatial distributions of Bejan number with varying Brinkman number and permeability parameter are shown in Fig. 7. Part a of this figure shows that at low Brinkman numbers the numerical value of Bejan number is relatively high for most of the domain. Brinkman number is proportional to the square of the flow strain rate, and thus, low Brinkman number indicates a low momentum impinging flow. This allows for the development of thick thermal and concentration boundary layers and renders higher values of Bejan number. Figure 7a shows that as the value of Brinkman increases, the region with finite value of Bejan number becomes increasingly small and ultimately it becomes limited to a narrow region close to the surface of the cylinder. As already discussed, the existence of blowing around the first quarter of the cylinder circumference thickens the boundary layers and enlargers the region of finite Bejan number. Figure 7b demonstrates the pronounced effect of the permeability parameter on Bejan number. With high permeability of the porous medium (low values of permeability parameter), the transfer processes take place more conveniently and therefore a large fraction of the investigated area experiences temperature and concentration gradients. This results in finite values of Bejan number for most of the domain. Increasing the permeability parameter, or decreasing the permeability of the porous medium, affects Bejan number in two different experiencing blowing of mass is clear. This is due to the slow motion of the fluid in this region (see Fig. 2), which reduces the frictional entropy generation. The influences of Biot number upon the thermal entropy generation are shown in Fig. 8b. This figure indicates that regardless of the value of Biot number, there is little thermal entropy generation in the region of the domain with blowing. This can be explained by noting that, as shown in Figs. 2-5, the slow flow in this region is almost isothermal and thus does not include strong temperature gradients and thermal entropy generation. Figure 8b further shows that increasing the Biot number and the resultant increase in the heat exchanges between the fluid and solid phases widen the region with finite thermal entropy generation. Tables 2-5 provide the numerical values of the average Bejan number for different values of the pertinent parameters. In particular, Table 4 shows that by increasing Reynolds number the value of Bejan number increases substantially. This is consistent with the behaviour observed in Fig. 5, in which increasing Reynolds number results in thinner thermal boundary layers and thus more intense thermal gradients. Also Table 4 shows that enhancing the permeability parameter reduces Bejan number. Significant increases in the flow friction at lower values of permeability of the porous medium and the subsequent increase in the flow irreversibility are the reason of this behaviour.

Nusselt and Sherwood numbers
The circumferential distributions of Sherwood and Nusselt numbers are depicted in Figs. 9-12. Figure 9 shows the angular variations of Sherwood number for different values of Reynolds number and permeability parameter. The strong effect of Reynolds number on the distribution of Sherwood number is evident in Fig. 9a. For low values of Reynolds number, the distribution is more and less uniform with an increase at around 0°/360°. This is because in the non-uniform transpiration case shown in Fig. 1, 0°is a singular point in which transpiration changes from suction to blowing. Thus, this point effectively provides a  Table 4, on average, increasing the Reynolds number boosts the value of Sherwood number on the cylinder. This observation is in keeping with the physical anticipation, which correlates Sherwood number with Reynolds number in forced convection problems. Figure 9b reflects the effects of permeability parameter on the distribution of Sherwood number. In general, permeability appears to have a modest effect on the Sherwood number (also see Table 4). This is particularly the case at around 0°and 360°f or which there is a sharp increase in Sherwood number. Further, it seems that Sherwood number varies only within a specific range of permeability parameter. For instance, increasing the permeability parameter from 1000 to 5000 in Fig. 9b results in almost no increase in Sherwood number. Figure 10a reveals a strong positive correlation between Damköhler number and Sherwood number. This correlation appears to be nearly linear, in which increasing Damköhler number from 1 to 4 increases the Sherwood number by almost four times. It should be noted that due to the implementation of a zeroth-order surface reaction in the current problem, increases in Damköhler number can be viewed as increasing the surface mass flux. This intensifies the radial gradient of the concentration and through Eq. (34) enhances the mass convection coefficient and Sherwood number. Figure 10b shows that the effect of Soret number on Sherwood number varies with the sign of Soret number. In general, Soret number can be either positive or negative [46][47][48]. It is clear in Fig. 11b  small thickness of the concentration boundary layer and hence the high strength of the primary mass transfer process in this region, which overrules the secondary influences of Soret effect. Figure 11 illustrates the effects of permeability parameter and Biot number on the numerical value of Nusselt number. Similar to that discussed with respect to Sherwood number in Fig. 9, Nusselt number appears to be rather insensitive to the changes in the permeability of the porous medium (see Fig. 11a).Yet, Biot number does have noticeable effects upon the Nusselt number. For u [8 0°, increasing Biot number results in an increase in Nusselt number. However, in the rest of the circumference the numerical value of Nusselt number decreases through increasing Biot number. This behaviour is consistent with that observed in Fig. 3 on the temperature field of the fluid phase. Figure 3a shows that for 0°[ u [ 80°, which coincides with the blowing transpiration, the thickness of thermal boundary layer decreases by increasing the Biot number. This causes an increase in the radial temperature Analysis of transport from cylindrical surfaces subject to catalytic reactions and non-uniform… gradient and thus enhances the Nusselt number. In the rest of the circumference, increasing Biot number thickens the thermal boundary layer and therefore the value of Nusselt number decreases. Table 5 shows that the surface-averaged value of Nusselt number decreases by increasing Biot number. The extent of this reduction is quite considerable, rendering Biot number an important factor influencing the overall rate of heat transfer. Figure 12 demonstrates the very significant effects of Reynolds number on Nusselt number. As expected, there exists a strong positive correlation between Reynolds number and Nusselt number, which becomes particularly noticeable for Re [ 10 (see also Table 4). Figure 12b shows that increasing Dufour number results in the reduction of Nusselt number. This finding can be confirmed by noting that in Fig. 4 Table 7 shows that increasing Prandtl number leads to a    Table 7. The preceding analyses have consistently reflected the strong effects of transpiration on the different characteristics of the system under investigation. Table 8 reports the results of a more systematic study on these effects in which the average Nusselt, Sherwood and Bejan numbers have been calculated for two different types of transpiration and for different values of Reynolds number. It is clear from Table 8 that the numerical values of all these dimensionless numbers are smaller in the case of non-uniform transpiration in comparison with those under uniform transpiration. The difference increases with increasing Reynolds number and becomes quite significant for Re C 1. This is particularly true for the average Nusselt number, for which changes in the type of transpiration can alter the numerical values by a few folds. The trend can be explained by noting the strong influences of transpiration upon the thickness of the thermal and concentration boundary layers as shown in Figs. 4-7. It is well known that convection coefficients and therefore Nusselt and Sherwood number are directly affected by changes in the boundary layer thicknesses. Transpiration in the form and blowing and suction of fluid change these thicknesses and either help or hinder the transport process. In the nonuniform transportation case, the overall effect is to thicken the boundary layers and therefore to decrease Nusselt, Sherwood and Bejan numbers.

Conclusions
A set of semi-similar solutions were developed for the problem of forced convection of heat and mass from the surface of a cylinder embedded in porous media and subject to a stagnation flow. The problem includes zerothorder chemical reactions on the surface of the cylinder as well as a non-uniform transpiration and an impinging flow. The conducted analyses considered the non-equilibrium thermodynamics including the local thermal non-equilibrium in the porous medium and coupled heat and mass transfer through Soret and Dufour effects. To the best of authors' knowledge, none of these have been investigated in the curved surfaces embedded in porous media. The major findings of this study can be summarised as follows: • It was shown that the numerical value of Biot number can considerably influence the temperature fields of the fluid and solid phases in the porous medium and thus affect the average Nusselt number. This clearly reflects the importance of considering local thermal nonequilibrium in the problem under investigation. • Variation of Dufour number appears to have considerable effects on the fluid temperature but does not noticeably change the temperature of the solid phase of the porous medium. • Changes in Reynolds number strongly affect the fluid temperature. However, they have much limited influences upon the temperature of the solid phase. Also, increases in Brinkman number enhance Bejan number in a highly non-uniform way. • It was shown that small variations in Soret and Dufour numbers can lead to noticeable changes in Nusselt and Sherwood numbers. This reflects the significance of these secondary mechanisms of transport and also highlights the importance of predicting the temperature and concentration fields accurately. • Non-uniform transpiration has a strong effect upon the temperature and concentration fields and hence majorly influences the average Nusselt, Sherwood and Bejan numbers. Further, the circumferential distribution of Sherwood and Nusselt numbers is shown to be highly sensitive to transpiration on the surface of the cylinder.
It follows that consideration of non-equilibrium thermodynamics and particularly local thermal non-equilibrium is an important necessity in the analysis of the porous thermochemical systems such as that investigated in this study.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creative commons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.