Article Combined Effects of Pipe Diameter, Reynolds Number and Wall Heat Flux and on Flow, Heat Transfer and Second-Law Characteristics of Laminar-Transitional Micro-Pipe Flows

Fluid flow, heat transfer and entropy generation characteristics of micro-pipes are investigated computationally by considering the simultaneous effects of pipe diameter, wall heat flux and Reynolds number in detail. Variable fluid property continuity, Navier-Stokes and energy equations are numerically handled for wide ranges of pipe diameter (d = 0.50–1.00 mm), wall heat flux (q''= 1000–2000 W/m2) and Reynolds number (Re = 1 – 2000), where the relative roughness is kept constant at e/d = 0.001 in the complete set of the scenarios considered. Computations indicated slight shifts in velocity profiles from the laminar character at Re = 500 with the corresponding shape factor (H) and intermittency values (γ) of H = 3.293→3.275 and γ = 0.041→0.051 (d = 1.00→0.50 mm). Moreover, the onset of transition was determined to move down to Retra = 1,656, 1,607, 1,491, 1,341 and 1,272 at d = 1.00, 0.90, 0.75, 0.60 and 0.50 mm, respectively. The impacts of pipe diameter on friction mechanism and heat transfer rates are evaluated to become more significant at high Reynolds numbers, resulting in the rise of energy loss data at the identical conditions as well. In cases with low pipe diameter and high Reynolds number, wall heat flux is determined to promote the magnitude of local thermal entropy generation rates. Local Bejan numbers are inspected to rise with wall heat flux at high Reynolds numbers, indicating that the elevating role of wall heat flux on local thermal entropy generation is dominant to the suppressing function of Reynolds number on local thermal entropy generation. Cross-sectional total entropy generation is computed to be most influenced by pipe diameter at high wall heat flux and low Reynolds numbers.


Introduction
With the progress of scientific investigations on energy management by developing new power supervision strategies, aiming to increase efficiency values and to suppress power loss amounts, the focus on this research branch is mainly directed towards the design and integration of micro-systems in actual industrial processes. The simultaneous advance of scientific efforts and the related manufacturing technologies recently made it possible for the micro-flow applications to take part in various engineering fields such as heat exchangers [1], fuel cells [2], pumps [3] and even in gas turbines [4]. Showing parallelism, both research outputs and industrial feedback signified the importance of surface roughness in micro-pipe flows by associating the roughness based viscous effects on the development of velocity and temperature profiles, flow-transition activity, entropy generation and power loss rates. Thus, due to these scientific fundamentals and industrial necessities, a research study, focusing on the roles of surface roughness in micro-pipe flows, must be structured on the three main issues: (i) frictional activity and flow characteristics, (ii) heat transfer mechanism and (iii) second-law analysis.
The roles of surface roughness on the friction characteristics, flow behavior and transitional activity in micro-pipe flows have been in the agenda of several researchers and there exist several published experimental and numerical reports. The role of frictional activity on the early transition from laminar to turbulent flow and on the grow of surface friction with surface roughness was reported by Guo and Li [5]. Celata et al. [6] determined, for the micro-pipe diameter range of d < 100 μm with the non-dimensional roughness of ε * < 0.01, that friction factor tends to deviate from the Poiseuille law when Re > 1300, which can be attributed to the acceleration associated with compressibility effect. As energy conversion of near-wall microfluidic transport was numerically investigated for slip-flow conditions, including different channel aspect ratios, pressure coefficients and slip flow conditions, by Ogedengbe et al. [7], Engin et al. [8] studied wall roughness effects in microtube flows and evaluated considerable shifts from the conventional laminar flow theory. Laminar and transitional flow experiments in dimpled tubes were performed by Vicente et al. [9], who reported a relatively low transition Reynolds number of 1,400 where the roughness induced friction factors were 10% higher than those of the smooth tubes. In a similar work, Kandlikar et al. [10] studied surface roughness effects on pressure drop in circular small hydraulic diameter tubes, where the transition onset appeared at Reynolds number values much below 2,300. Wen et al. [11] experimentally investigated the augmentation of pressure drop characteristics due to different strip-type inserts in small tubes. As Koo and Kleinstreuer [12] numerically and experimentally studied the significance of viscous dissipation on the friction factor, the roles of surface roughness on viscous dissipation, increased friction factor values and the corresponding earlier transitional activity was expressed by Celata et al. [13,14]. Morini [15] handled the concept of viscous dissipation for various duct cross-section types, and pointed out the minimum Reynolds number for which viscous dissipation effects can not be neglected. A comprehensive review, on the interaction of surface roughness with frictional activity and momentum rates in microchannels, was carried out by Obot [16]. He reported that a wide collection of available literature identified the definition of microchannels with a hydraulic diameter of d h ≤ 1,000 μm (1 mm), moreover, the limiting Reynolds number for the onset of transition in smooth microchannels was above 1,000.
The role of surface roughness in thermal systems was also taken into consideration in recent studies; some dealt with the roughness definition and others focused on the heat transfer mechanism. Wang et al. [17] investigated the influence of two-dimensional roughness on laminar flow in microchannels between two parallel plates by applying the regular perturbation method. Sheikh et al. [18] developed a model to eliminate the discrepancy in the fouling measurements by characterizing fouling as a correlated random process. In a complementary work, the role of fouling on entropy generation, as compared to that for clean surface tubes, was studied by Sahin et al. [19]. Cao et al. [20] employed non-equilibrium molecular dynamics simulation to investigate the effect of the surface roughness on slip flow of gaseous argon. More comprehensive research comprised both flow and heat transfer characteristics driven by roughness, such as Wu and Cheng [6], who reported on the Nusselt number augmentations with surface roughness, Kandlikar et al. [10] findings on heat transfer manipulation due to roughness, Wen et al. [11] work on the effects of the imposed wall heat flux, mass flux and strip inserts on heat transfer rates in small diameter tubes, Koo and Kleinstreuer's [12] evaluations on the temperature field with viscous dissipation, a review by Obot [16] indicating the variation of laminar Nusselt numbers with the square root of the Reynolds number and finally Ozalp's [21,22] numerical compressible flow investigations for aerospace propulsion applications.
Since the essential concern of thermo-fluid based research studies is to generate a complete overview on the schematic efficiency features, second-law analysis recently became prominent due to its potential to describe the overall system performance. Kotas et al. [23] indicated that for calculating loss of exergy, or process irreversibility, an exergy balance or the Gouy-Stodola theorem can be used. Zimparov [24] reviewed passive heat transfer augmentation techniques for single-phase flows. The effects of the spanwise fin tube spacing on the second-law performance of wavy plate fin-and-tube heat exchangers were investigated by Lin and Lee [25]. Ratts and Raut [26] obtained optimal Reynolds numbers for single-phase, convective, fully developed internal laminar and turbulent flows with uniform heat flux, by employing the entropy generation minimization method. Entropy generation for a fully developed laminar viscous flow in a duct subjected to constant wall temperature was analytically studied by Sahin [27]. His computations not only show that temperature dependency of viscosity becomes essentially important in accurately determining the entropy generation, but also indicate that for low heat transfer conditions the entropy generation due to viscous friction becomes dominant. Rakshit and Balaji [28] reported the results of a numerical investigation of conjugate convection from a finned channel, where the vertical rectangular fins were mounted outside the horizontal channel. The effects of streamwise variation of fluid temperature and rib height to diameter ratio on the entropy production of a tubular heat exchanger with enhanced heat transfer surfaces were investigated by Zimparov [29]. In singly connected microchannels with finite temperature differences and fluid friction, Richardson et al. [30] seeked the optimum laminar flow regime based on second-law analysis. Ko [31] computationally carried out the thermal design of plate heat exchanger double-sine ducts, from the point of entropy generation and exergy utilization. Laminar forced convection and entropy generation in a helical coil with constant wall heat flux was also numerically investigated by Ko [32]. He suggested that the optimal Reynolds number is to be chosen according to the flow operating condition, for least thermal irreversibility and best exergy utilization The past research categorized above clearly indicates the wide scientific prospects and the industrial importance and actuality of the micro-pipe flows with surface roughness. As the available literature designates, real-time operational and structural parameter ranges in micro-pipe flows are broad, which as a consequence significantly promotes the possible process and output scenarios due to the individual and combined parameter influences. Recently, Ozalp [33][34][35] computationally investigated the momentum transfer, heat transfer and second-law characteristics of air flow in a circular micro-pipe with a diameter of d = 1 mm. The analyses were carried out in the Reynolds number and non-dimensional surface roughness ranges of Re = 1-2,000, ε * = 0.001-0.05. The author continued his research by expanding the micro-pipe diameter and non-dimensional surface roughness limits to d = 0.10-1.00 mm and ε * = 0.001-0.1 for water flow; this content is recently supported by the Uludag University Research Fund. In the present article, the author reports the initial finding set of the ongoing project, comprising the first and second-law characteristics of roughness induced laminar-transitional micro-pipe flow, for the diameter range of d = 0.50-1.00 mm with the fixed non-dimensional surface roughness condition of ε * = 0.001. Radial velocity and temperature profiles, boundary layer parameters, friction coefficients, Nusselt numbers, frictional energy loss and temperature rise values are discussed to identify the flow and heat transfer characteristics. Additionally, second-law records are interpreted through radial profiles of thermal and frictional entropy values and cross-sectional total entropy generation rates and local and cross-sectional average Bejan numbers.

Micro-Pipe and Roughness
The diameter and length of the micro-pipe (Figure 1a) analyzed in the present paper, are given as d and L, respectively. The present roughness model is based on the triangular structure of Cao et al. [20] ( Figure 1b), where the roughness amplitude is characterized by ε and period by ω. In all computations, the roughness periodicity parameter (ω' = ω/ε) is kept fixed to ω' = 2.31, which corresponds to the equilateral triangle structure (Cao et al. [20]). Equation 1 numerically characterizes the model of Cao et al. [20] with the implementation of the amplitude and period. The model function (f ε (z)) is repeated in the streamwise direction throughout the pipe length, where the Kronecker unit tensor (δ i ) attains the values of δ I = +1 and -1 for 0 ≤ z ≤ ε 2 2.31 and ε 2 2.31 ≤ z ≤ 2.31ε respectively.

Governing Equations
The density, dynamic and kinematic viscosity of water are defined with ρ, μ and ν (=μ/ρ); the average velocity of the flow at any cross-section of the duct is denoted by U o . As the surface and mean flow temperatures are denoted by T s and T o , thermal conductivity and convective heat transfer coefficient are characterized by κ f and h. Using these definitions, Reynolds number and Nusselt number are given by Equations 2a-b. Water has been selected as the working fluid in the present study, and it is well known that water properties like specific heat (C p ), kinematic viscosity (ν) and thermal conductivity (κ f ), are substantially dependent on temperature (Incropera and DeWitt [36]). To sensitively implement the property (ξ) variations with temperature into the calculations, the necessary water data of Incropera and De Witt [36] are fitted into 6 th order polynomials, which can be presented in closed form by Equation 2c. The uncertainty of the fitted water data is less than 0.03% and the temperature dependency is indicated by the superscript T throughout the formulation.
The problem considered here is steady ( ) 0 , fully developed and the flow direction is coaxial with pipe centerline ( ) , thus the velocity vector simplifies to . These justifications are common in several recent numerical studies on roughness induced flow and heat transfer investigations, like those of Engin et al. [8], Koo and Kleinstreuer [12] and Cao et al. [20]. The flow boundary conditions are based on the facts that on the pipe wall no-slip condition and constant heat flux exist, and flow and thermal values are maximum at the centerline. As given in Figure 1a, , pressure and temperature values are known at the pipe inlet and the exit pressure is atmospheric. Denoting U z =U z (r) and T=T(r,z), the boundary conditions can be summarized as follows: Since the viscous stress τ zz is characterized by the variation of flow velocity in the axial direction Application of wall heat flux ( ' ' q ) on the lateral surface of the micro-pipe results in the total surface heat transfer of Q s = ' ' q πdL. With these problem definitions and above implementations for a laminar incompressible flow with variable fluid properties, continuity, momentum and energy equations are given as: The cross-sectional average fluid velocity and temperature are obtained from and the shear stress and mass flow rate are defined as Gouy-Stodola theorem [23] considers the lost available work to be directly proportional to entropy generation, which results from the non-equilibrium phenomenon of exchange of energy and momentum within the fluid and at the solid boundaries. Due to the existence of the velocity and temperature gradients in the flow volume, the volumetric entropy generation rate is positive and finite. Computation of the temperature and the velocity fields through Equations 5-7 on the problem domain will produce the input data for Equation 10a. Equation 10a defines the local rate of entropy generation per unit volume ( ) As given in Equation 10b, the entropy generation due to finite temperature differences ( ) in axial z and in radial r directions is defined by the first term on the right side of Equation 10a, the second term defines the frictional entropy generation ( ) . Equation 11a stands for the cross-sectional total entropy generation rate; local and cross-sectional average values of Bejan numbers can be obtained by Equations 11b-c.

Computational Method
The author employed the present computational method in his former work [22,[33][34][35] and detailed the essentials and the sub-steps of the procedure extensively. Thus, here only the fundamental issues of the methodological structure are defined. To simultaneous handle the continuity, momentum and energy mechanisms of the flow, Equations 5-7 are assembled into the three-dimensional "Transfer Matrix", consisting of the converted explicit forms of the principle equations. The flow domain of Figure 1a is divided into m axial and n radial cells (m × n), where the fineness of the computational grids is examined to ensure that the obtained solutions are independent of the grid employed. After a series of successive runs the optimum axial and radial cell numbers are determined as m = 500→850 and n = 100→225 respectively, for d = 1.00→0.50 mm. Forward difference discretization is applied in the axial and radial directions in the two-dimensional marching procedure. To sensitively compute the velocity and temperature gradients on the pipe walls, the 20% of the radial region, neighboring solid wall, is employed an adaptive meshing with radial-mesh width aspect ratio of 1.1→1.05 (d = 1.00→0.50 mm). The influences of surface roughness and surface heat flux conditions, over the meshing intervals of the flow domain, are coupled by Direct Simulation Monte Carlo (DSMC) method. The author previously applied DMSC to compressible nozzle flow problems [21,22] and to micro-pipe flow scenarios with surface roughness [33][34][35]. The concept of triple transport conservation is as well incorporated into the DSMC, which makes it possible to sensitively evaluate the balance of heat swept from the micro-pipe walls, the energy transferred in the flow direction and also to perform accurate simulation for inlet/exit pressure boundaries. The "Transfer Matrix" scheme and the DSMC algorithm are supported by cell-by-cell transport tracing technique to guarantee mass conservation, boundary pressure matching and thermal equilibrium within the complete mesh. The resulting nonlinear system of equations is solved by using the Newton-Raphson method and in the case that the Reynolds number does not fit the required value, inlet pressure and temperature are both modified by DSMC, to regulate the Reynolds number of the former iteration step.

Results and Discussion
To generate a scientifically detailed and academically enlightening scheme for the combined effects of pipe diameter, Reynolds number and wall heat flux on the flow, heat transfer and second-law characteristics of laminar-transitional micro-pipe flows, computations are carried out in the ranges of d = 0.50-1.00 mm, Re = 1-2,000 and ' ' q = 1,000-2,000 W/m 2 . As the considered micro-pipe diameter range is consistent with the microchannel definition of Obot [16] (d ≤ 1.00 mm), the imposed wall heat flux values are decided in conjunction with the Reynolds number and the accompanying mass flow rate ranges, to bring about applicable and rational heating. The fixed parameters of the analyses are the length of the micro-pipe (L = 0.5 m), inlet temperature (T in = 278 K) and exit pressure (P ex = 0 Pa) of water and the non-dimensional surface roughness (ε * = 0.001), which is in harmony with those of Engin et al. [8] (ε * ≤ 0.08), Sahin et al. [19] (ε * ≤ 0.25), Ozalp [33][34][35] (ε * = 0.001-0.05) and Vijayalakshmi et al. [37] (ε * = 0.0004-0.0018). Flow and heat transfer characteristics are presented in terms of radial velocity and temperature profiles, boundary layer parameters, friction coefficients, Nusselt numbers, frictional energy loss and temperature rise values. Radial profiles of thermal and frictional entropy values, cross-sectional total entropy generation rates and Bejan numbers are discussed and identified to interpret the second-law records. Figure 2 displays the velocity profile (VP) formation in the Reynolds number range of Re = 100-2,000, for the pipe diameters d = 1.00, 0.75 and 0.50 mm with the fixed non-dimensional surface roughness of ε * = 0.001. To clearly display the deviations of the VPs from the characteristic styles of laminar and turbulent regimes and also to put forward the roles of d and Re on the transition mechanism, the characteristic laminar profile (Equation 12a) and the modified turbulent logarithm law for roughness (Equation 12b) (White [38]) are also displayed in Figure 2, where U * is the friction velocity (Equation 12c).

Radial schemes
To clarify the roles of micro-pipe diameter and Reynolds number on the frictional activity, the classical and normalized friction coefficient values are evaluated by Equations 13a-c [38]. To strengthen the discussions on transition, boundary layer parameters like shape factor (H) and intermittency (γ) [38] are estimated by Equations 14a and b, respectively.
As the laminar (H lam = 3.36) and turbulent (H turb = 1.70) shape factor values are computed with Equation 14a, by integrating the laminar (Equation 12a) and turbulent (Equation 12b) profiles, the shape factor data of the transitional flows were also calculated with Equation 14a, however with the computationally evaluated corresponding velocity profiles. The computed boundary layer and transition parameters are presented in Table 1. It can be inspected from Figure 2 that in the all pipe diameter cases investigated, the VPs for the Re = 100 flow fit with the laminar profile, where the shape factor values are evaluated as H ≥ 3.345 (Table 1) being very close to the characteristic laminar value of H = 3.36. In the Re = 500 case, the VPs are still in harmony with the laminar profile; however slight shifts from the laminar character can be detected, such that the center velocity values and the shape factor and intermittency data build up the ranges of U c /U o = 1.960→1.950, H = 3.293→3.275 and γ = 0.041→0.051 (d = 1.00→0.50 mm), respectively.
Moreover, in the scenarios with higher Re, recognizable deviations in the VPs are determined. Thus to identify the modifications, the near wall (0.998 ≤ r/R ≤ 1) and pipe-centerline (0 ≤ r/R ≤ 0.4) regions are demonstrated through zoomed plots. Computations put forward that in the systems with higher Re and/or lower d as C f elevates, the VPs depart significantly from the laminar profile causing the center velocity figures (U c /U o ) to decrease, where these evaluations can also be visualized from the zoomed plots and Table 1 Table 1). The near wall (0.998 ≤ r/R ≤ 1) zoomed plots of Figure 2 clarify that for flows with lower micro-pipe diameter the impact of Reynolds number on VPs is more noticeable. This evaluation is significant to interpret the combined effects of Re and d on the momentum transfer and transition mechanism in micro-pipe flows. In the light of several numerical and experimental research studies that examine the roles of micro-pipe diameter and surface roughness on the transition mechanism in internal flow problems, many of the involved scientists [5,9,10,16,39] accept the 10% rise of C f (C f * = 1. Re tra = 500 for d = 0.01 mm. The basis of this divergence may be due to the undefined/unavailable surface roughness information in many of the works. However, the present findings on transition onset still not only conform with the associated reports but also enriches the related research frame with scientific outputs on well defined scenarios. It can additionally be extracted from the present Re tra that, as the gap of (Re tra ) d=1.00 mm -(Re tra ) d=0.90 mm is 49, the corresponding one for (Re tra ) d=0.60 mm -(Re tra ) d=0.50 mm becomes 69. These values clearly pronounce that, the lower the micro-pipe diameter is reduced the more the sensed influence on Re tra is promoted. The tabulated outputs in Table 1 [36]. Moreover, the findings on the thermal behavior of micro-pipe flows are supported with the tabulated records of Ψ loss and flow temperature rise (ΔT) in Table 2.  It can be inspected from Figure 3 that more noticeable deviations in TPs from the laminar CHF profile (Equation 15) arise in the flow scenarios with lower micro-pipe diameter and heat flux and higher Reynolds number. This determination can be explained by the growth of Ψ loss values under these flow conditions ( Table 2). The higher energy loss outputs, with lower d and higher Re, count on the enhanced wall velocity gradients ( Figure 2) and the corresponding stronger friction activity (Table 1). On the other hand, the source of the raised Ψ loss with lower ' ' q can be enlightened with the decay of ΔT data ( Table 2) Analyses show that as ΔΨ loss gets the values of 7.6×10 -6 W and 1.9×10 -4 W for Re = 100 and 2,000 in the micro-pipe flow of d = 1.00 mm, these data rise to 1.41×10 -5 W and 3.4×10 -4 W at d = 0.75 mm and to 3.07×10 -5 W and 7.0×10 -4 W at d = 0.50 mm. These findings indicate that the potential of ' ' q on Ψ loss develops with lower d and higher Re. This determination is quite significant in identifying the scenarios where the TPs are more probable to deviate from the laminar constant heat flux profile. Since ΔT and Ψ loss respond inversely to the imposed parameters, they develop completely in a contrary mechanism under the influence of these operational factors. Furthermore, simultaneous inspection of the depended character of Ψ loss on the VPs and the affiliation of ΔT with TPs evidently put forward the strong interaction of the momentum and thermal characteristics in micro-pipe flows.  deviation of (C f * ) d=1.00 mm and (C f * ) d=0.50 mm at Re = 100 is only 0.1%, it becomes 0.97%, 1.13% and 2.57% at Re = 500, 1,000 and 1,500, respectively. These deviations put forward that the impact of micro-pipe diameter on momentum characteristics becomes more significant at high Reynolds numbers. The complete set of presentations and analysis on the characterization of micro-pipe flows with Reynolds number and micro-pipe diameter not only ascertained the decrease of H and U c /U o and increase of C f * and γ values, but also identified the shift process of VPs from fully laminar to the transitional structure. As tendency to transitional activity and promoted friction coefficients account to the augmentation of pressure loss amounts of the flow, the present findings on the individual and combined roles of Reynolds number and micro-pipe diameter are quite remarkable for micro-level industrial applications such as heat exchangers, pumps, electronics cooling, piping arrangements and lubrication systems. On the other hand, Figure 4 as well shows that although the C f * values of Yu et al. [40] were even below the laminar theory for 100 < Re < 2,000, the experimental friction coefficient data of Vijayalakshmi et al. [37] (ε * = 0.002), Choi et al. [41] (ε * = 0.0003), Kohl et al. [42] (ε * = 0.003-0.013) and Wu and Little [39] exposed strengthening trends for Re > 500, showing harmony with the present numerical evaluations.  To provide a more sensitive designation on the individual role of micro-pipe diameter on Nusselt number, the ratio of Nu d=0.50 mm /Nu d=1.00 mm is taken into consideration, which proposed the variation rates of 1.065, 1.126, 1.181 and 1.231 at the Reynolds number levels of Re = 500, 1,000, 1,500 and 2,000 respectively. These proportions clearly identify the growing influence of the micro structure on the heat transfer behavior at high Reynolds numbers. Besides, comparing the Nusselt numbers evaluated at the transition Reynolds numbers (Nu tra ) with the laminar typical value of Nu lam = 4.25 points out the augmentation rates of Nu tra /Nu lam as 1.30 and 1.41 at d = 1.00 mm and d = 0.50 mm, respectively. These values put forward that at the onset of transition, heat transfer activity is more potentially provoked at lower micro-pipe diameter cases. It can additionally be extracted from these proportions that, for the micro-pipe diameter range of d = 0.50-1.00 mm with the fairly low relative roughness of ε * = 0.001, the Nusselt numbers rise by a factor range of 1.30-1.40 at the transition onset.    Table 1), frictional entropy generation rates are as well computed to enhance with higher Re and lower d. This outcome can be enlightened through the elevated fluid flow-solid domain interaction in cases with narrower flow domains (lower micro-pipe diameter) or increased mass flow rate (higher Reynolds number). Present computational determinations are completely in harmony with the results of Ko [31] for double-sine ducts, on the variations of the thermal and frictional entropy generation rates with Reynolds number. The individual and combined roles of micro-pipe diameter and Reynolds number can be distinguished more specifically by comparing the wall (r/R = 1) frictional entropy generation values of different micro-pipe diameters and Reynolds numbers (Equations 19a-b). Numerical analyses show that α 1 has the values of 1.00, 5.62 and 16.90 for Re = 100, 1,000 and 2,000, respectively; besides α 2 is computed as 30.99 and 523.42 for the micro-pipe diameters of d = 1.00 mm and 0.50 mm. The grow of α 1 with higher Re clearly states that the influence of micro character on frictional entropy generation becomes more impressive at higher Reynolds numbers. Besides, the provoked view of α 2 at lower d  The above findings completely enlighten the chief configurative stance of heat flux among the characterizing parameters of Bejan number in micro-pipe flows. It can additionally be extracted from ζ 1 that, the overall scheme of the above defined characteristic actions on ' ' ' Be is more significant in the Reynolds number range of Re = 100-1,000; for Re > 1,000 the momentum of the augmenting role decreases. This evaluation puts forward that the exhibited functions is more decisive in laminar flow.

Analyses on Second-Law Characteristics
On the other hand, ζ 2 figures out the potential of micro-pipe diameter on Bejan number with the descriptive numbers of 1.00, 17.36 and 29.55 for Re = 100, 1,000 and 2,000 at ' ' q = 1,000 W/m 2 and1.00, 13.56 and 20.60 at ' ' q = 2,000 W/m 2 . ζ 2 puts forward complementary indicators, not only identifying the more noticeable role of micro-pipe diameter on Bejan number in scenarios with high Reynolds number and low heat flux but also pointing out the comprehensive motivation of micro-pipe diameter on Bejan number in the Reynolds number range of Re = 100-1,000.   Be with high Re (Figure 8). Since entropy generation is a fundamental scientific indicator of the amount of thermodynamic irreversibilities, the present foresights on the simultaneous actions and roles of Reynolds number, pipe diameter and surface heat flux on the second-law characteristics of micro-pipe flows can be considered as guiding tools in handling the second-law efficiency aspects of micro-level power generation systems or turbomachines, such as gas turbines, compressors, turbines and nozzles.  Variation of cross-sectional average Bejan numbers with surface heat flux, Reynolds number and micro-pipe diameter is given in Figure 10. As expected from the findings and discussions through Figure 8, Bejan number decreases with high Reynolds number and low micro-pipe diameter and heat flux. To provide a deeper inspection into the second-law mechanism of Be, based on flow (Re), thermal  Be variation for r/R = 0→1 is mostly matching. On the other hand, π 2 data differ from those of ζ 2 , especially at the high Reynolds number range of Re = 1,000-2,000.

Characteristic parameters
At ' ' q = 1,000 W/m 2 , π 2 are evaluated as 1.00, 16.77 and 28.05 for Re = 100, 1,000 and 2,000 and become 1.00, 12.86 and 19.35 at ' ' q = 2,000 W/m 2 . The slightly lower (<6%) view of π 2 than ζ 2 puts forward that the effectiveness of the extreme centerline Be , and also the cross-sectional average Bejan number of Be as well.

Conclusions
Temperature dependent property continuity, momentum and energy equations are simultaneously solved to computationally investigate the flow, heat transfer and second-law characteristics of laminar-transitional micro-pipe flows. The combined effects of pipe diameter, Reynolds number and wall heat flux are presented and discussed in detail.
Analyses on flow characteristics resulted in the following main findings: Friction coefficients increase with higher Reynolds number and/or lower micro-pipe diameter, with comprehensive velocity profile deviations from the laminar profile. As the center velocity values decrease, due to the rise of wall velocity gradients, non-dimensional friction coefficients transform towards the transition with the interpreting figures of C f * = 1.006→1.007 (d = 1.00→0.50 mm) at Re = 100, C f * = 1.028→1.038 at Re = 500, C f * = 1.063→1.075 at Re = 1,000 and C f * = 1.088→1.116 at Re = 1,500. Towards the transition regime, shape factor data drop below the laminar regime indicator of H = 3.36 and the intermittency elevates. As the combined influence of micro-structure and surface roughness cause earlier transition at the descriptive Reynolds numbers of Re tra = 1,656, 1,607, 1,491, 1,341 and 1,272 at d = 1.00, 0.90, 0.75, 0.60 and 0.50 mm respectively, the boundary layer parameters appear in the ranges of H = 3.133→3.159 and γ = 0.137→0.121 (d = 1.00→0.50 mm).
Analyses on heat transfer characteristics identified the following primary issues: The wall temperature gradients are evaluated to rise with lower micro-pipe diameter and higher Reynolds number; besides energy loss come out to be notably excited in these scenarios . Energy loss values additionally rise at lower heat flux applications; moreover the impact of surface heat flux on energy loss is determined to grow at higher Reynolds number and lower micro-pipe diameter. Due to the direct relation of Reynolds number and mass flow rate, Reynolds number inversely motivates the temperature rise of water in the streamwise direction. Heat transfer rates are determined to increase in flows with lower micro-pipe diameters and higher Reynolds numbers. The growing impact of the micro-structure on the heat transfer activity at high Reynolds numbers is revealed by the Nu d=0.50 mm /Nu d=1.00 mm ratios of 1.065, 1.126, 1.181 and 1.231 at Re = 500, 1,000, 1,500 and 2,000 respectively. Moreover, at the transition onset, Nusselt numbers are computed to augment above the typical laminar value (Nu lam = 4.25) by 30→40% (d = 1.00→0.50 mm), pointing out the rising strength of micro-flow affects on the heat transfer mechanism at lower micro-pipe diameters.
Analyses on second-law characteristics pointed out the following essential items: It is determined that, at low Reynolds numbers (Re = 100), the augmenting capacity of heat flux on local thermal entropy generation is almost fixed in the complete micro-pipe diameter range of d = 0.50-1.00 mm; however the influence grows with both high Reynolds number and low micro-pipe diameter. Due to the intensity of the velocity gradients, local frictional entropy generation rises towards the pipe walls.
Computations exhibited a synergy imposed on local frictional entropy generation with the combined action of higher Reynolds number and lower micro-pipe diameter. The promoted cross-sectional total entropy generation values in lower micro-pipe diameter flows can be attributed not only to the forced local frictional entropy generation with lower diameters but also to the overwhelming increase of local frictional entropy generation to the decrease rate of local thermal entropy generation in lower micro-pipe diameters. In scenarios with high Reynolds numbers and low micro-pipe diameters and heat flux levels, Bejan numbers are monitored to noticeably decrease. In flows with low micro-pipe diameter, local Bejan numbers, around the centerline vicinity, are comprehensively decisive on the cross-sectional average Bejan numbers, which is relatively important from the point of thermodynamic interactivity of the radial variations of local thermal and frictional entropy generation rates and local and cross-sectional average Bejan numbers.