TEMPLATE FOR ACADEMICA SCIENCE JOURNAL AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 11, No. 4 ISSN: 1998-4456 Page 426 Copyright  2018 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. NUMERICAL INVESTIGATION OF ROUGNESS EFFECTS ON HYDERDYNMIC AND THERMAL PERFORMANCE OF COUNTER FLOW MICROCHANNEL HEAT EXCHANGER Ahmed A. Ali Mechanical Engineering Department, College of Engineering, Thi-Qar University, Iraq E mail: eng.ahmedameen91@gmail.com Mushtaq I. Hasan Mechanical Engineering Department, College of Engineering, Thi-Qar University, Iraq E mail: mushtaq76h@gmail.com Ghassan Adnan Mechanical Engineering Department, College of Engineering, Thi-Qar University, Iraq E mail: Ghassanadnan37@gmail.com Abstract: In this paper the effect of surface roughness on the performance of counter flow microchannel heat exchanger has been numerically investigated. The studied Microchannel heat exchanger is a square shape and made of aluminum as substrate material with different values of hydraulic diameters (20, 50, 110, 150 ) μm. The working fluid used is water at constant properties. Roughness- viscosity model has been used to study the roughness effect with 0.14 ratio of roughness to hydraulic diameter. The results obtained indicate that pressure drop of (CFMCHX) increased with increasing surface roughness and decrease hydraulic diameter also the results showed that there is a slight increasing in thermal performance with increasing the surface roughness. Keywords: surface roughness, Microchannel heat exchanger, Laminar flow, Numerical investigation, Thermal performance mailto:eng.ahmedameen91@gmail.com mailto:mushtaq76h@gmail.com mailto:Ghassanadnan37@gmail.com AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 11, No. 4 ISSN: 1998-4456 Page 427 Copyright  2018 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. NOMENCLATURE: Symbol Description SI unit Cp specific heat J\kg.K Dh hydraulic diameter m F friction factor H height m K height of the roughness element m . m Mass flow rate kg/s P.P Pumping poewer Watt T Temperture k Uk velocity at the top of roughness m/s uavg average velocity m/s U fluid x-component velocity m/sec Re Reynold number Rek local Reynold number RVM roughness viscosity model W width m ΔP pressure drop across heat exchanger pa Ρ density kg/m3 μR roughness viscosity m2/s Μp apparent viscosity m2/s  * Performance factor ɛ Heat exchanger effectiveness AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 11, No. 4 ISSN: 1998-4456 Page 428 Copyright  2018 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. INTRODUCTION In recent years, advanced micro fabricated technology as heat exchanger received large interest and increase the using in many applications due to heat transfer dissipation over small size this invites to use it in biomedical process, microelectronics ,metrology and aerospace application so the study of the overall performance of microchannel heat exchanger with different parameters is very important[1]. Decreasing the hydraulic diameter has a considerable effect on all parameters also the roughness effect play important role on hydraulic and thermal performance of microchannel heat exchanger (MCHX). There are many researches in literature studied many parameters affecting the performance of microchannel heat exchanger including roughness effect. Below a review of most related papers. Mala and Li (1999) [2] presented experimental work to investigate the performance on microchannel with diameter ranging from (50 - 254) μm by using roughness - viscosity model. [RVM] has been used to study the roughness effect on microchannel. [RVM ] treat the roughness effect as apparent viscosity in momentum equation. The apparent viscosity μp is the sum of viscosity of fluid μf and roughness viscosity μR. The authors solved the momentum equation analytically for cylinder with adding the roughness viscosity. [RVM] appears large agreement with experimental data in varying hydraulic diameter. Their results showed that, the early transition from laminar to turbulent which cause it the effect of surface roughness. Weilin et al (2000)[3] studied experimentally micro aching techniques fabricated the trapezoidal channel with relative roughness value (2.4 % - 3.5%) due to silicon surface finish and hydraulic diameter ranging ( 51-169 ) μm .they concluded that the friction factor and pressure drop in microchannel higher than those given by conventional theory due to the surface roughness effect. Qu et al (2000) [4] Studied the effect of surface roughness on heat transfer characteristics of water as working fluid flow through trapezoidal silicon microchannel with the hydraulic diameter ranging from ( 62-169) μm. They found that, the velocity near the wall surface is decreased due to the effect of surface roughness. They proposed a roughness – viscosity model [RVM] to interpret their experimental data based on this model. Hassan. et al (2009) [5] evaluated the effect of channel geometry fluid flow and heat transfer behavers of counter flow microchannel heat exchanger by using numerical simulation. The effects of different channel shapes on performances of (CFMCHX) are studied for five shapes (square, rectangular, triangular, trapezoidal and circle). The results showed that, for the same volume of heat exchanger, increase number of channel leads to increase both pressure drop and effectiveness also they developed a new correlations to predicted the value of heat exchanger effectiveness and performance index as function of relative size channel with Re, overall volume of heat exchanger and thermal conductivity ratio.Dang.et al (2010) [6]presented experimental testes and numerical simulation to study the fluid flow and heat transfer characteristics of counter microchannel heat exchanger for rectangular shaped has hydraulic diameter of 375 μm and heat flux applied was 17.38 w/cm2 was achieved for water from cold side having inlet temperature of 21.5 and mass flow rate 0.2906 g/sec and for hot side having inlet temperature 70 and mass flow rate 0.2503 g/sec. from result obtained that, AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 11, No. 4 ISSN: 1998-4456 Page 429 Copyright  2018 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. pressure drop increase as water temperature rises where pressure drop increases from 880 to 4400 pa. also obtained a good agreement between numerical and experimental results. Byongjoo kim (2016) [7] Investigated microchannel to study the flow resistance and thermal behavior in ten rectangular microchannel with aspect ratio ranging from (0.25 -3.8). He used FC 770 and deionized water as working fluids. The relative roughness measured by digital device was 3.7 mm. He observed through experimental demonstrate large deviation in poiseuille number(fRe) because the surface roughness as the channel size decrease. Also the result showed that ,the existing roughness in microchannel cause early transition to turbulent flow. In this paper the effect of the surface roughness on the performance of counter flow microchannel heat exchanger will be numerically studied by using roughness – viscosity model on thermo- hydraulic characteristics .In addition to relation between channel and roughness effect. 1. MATHMATICAL MODEL: The schematic drawing in fig (1) shows a physical model of problem. Which represents counter flow microchannel heat exchanger (CFMCHX) consists of square channels with hot and cold fluids where the working fluid ( water) inters in upper channel as cold fluid while in the lower as hot fluid. There are many diffcults with soling complet heat excanger, where need much time and super computers there fore one unit is used in this study. The selected unit consists of hot and cold channel and separating wall between them. For each channel the length (L) equal to 500 μm, width of the channel( W) , (ts) is the thickness of separating wall between cold and hot fluid of (5 μm) and the aspect ratio equals to one as seen in fig.2 The governing equations used in this paper based on the following physical and geometry assumptions :  The flow is laminar and steady.  The value of Knudsen number is smaller thane 0.001 for continuance fluid (no slip).  Water used in this paper is incompressible, constant properties and Newtonian.  There is no heat transfer to\from the ambient  The energy dissipation is negligible.  The pressure drop in axial direction only.  Three dimensional flow and heat transfer . The governing equations and boundary conditions in Cartesian are write below [8] [9][10] Continuity equation: AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 11, No. 4 ISSN: 1998-4456 Page 430 Copyright  2018 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. (1) Momentum equation: (2) Energy equation: (3) Mala and Li [2] modified roughness viscosity model where surface roughness in laminar flow considered in terms roughness – viscosity function . the apparent viscosity in momentum equations become viscosity of fluid sum roughness viscosity as equ.9 . the ratio of roughness to hydraulic diameter is 0.14 and k is average height of surface roughness. ) (4 2)) k r -exp(-(1 k r A )(4) =  R (5) k r=k ) r u   = (kU (6) /v 2k r =R) r u   = (kRe (7) = r =R) r u   ( A=0.1306 ( k R )0.3693 exp {Re (6 10-5 k R - 0.0029) (8) The apparent viscosity of fluid is µp=µ+ µR (9) 0V TkTVC p 2 )(  ).().( VPVV   AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 11, No. 4 ISSN: 1998-4456 Page 431 Copyright  2018 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. uci Tci Hc C ts H Hh L H Thi W uhi y Fig.1.Model of CFMCHX[10] Fig (2) unit cell of CFMCHX[10] AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 11, No. 4 ISSN: 1998-4456 Page 432 Copyright  2018 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. 2.BOUNDERY CONDITIONS: The boundary conditions used to complete the model are: For lower channel (hot fluid) Location B.c comments At z =0 u=uhi , vh= wh=0 , Th=Thi hot fluid in flow At z= L = vh = wh =0 , =0 hot fluid out flow (fully developed flow , End of channel ) At x =0 uh = vh = wh =0 , no- slip, adiabatic wall At x = w uh = vh = wh =0 , no- slip, adiabatic wall At y = 0 uh = vh = wh =0 , no- slip, adiabatic wall At y =Hh uh = vh = wh =0 ,-kh = - ks , fluid – solid interface Ts =Th (no –slip conjgat heat transfer) For upper channel ( cold fluid) (Hh + t  y  Hh + t +Hc ) Location B.c comments At x = 0 = vc =wc =0 , = 0 cold fluid outlet (fully developed Fully developed flow end of channel At x= w uc= uci , vc =wc = 0 , Tc= Tci cold fluid inflow At z = 0 uc= vc =wc =0 , = 0 No –slip adiabatic wall AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 11, No. 4 ISSN: 1998-4456 Page 433 Copyright  2018 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. At z =L uc= vc =wc =0 , = 0 No –slip adiabatic wall At y = Hh + t uc= vc =wc =0 , - kc = - ks , Tc =Ts Fluid –solid interface ( no slip, conjugate heat transfer) At y = Hh + t +Hc uc= vc =wc =0 , (no slip adiabatic wall ) For solid wall separating two channels (Hh  y  Hh + t) Location B.C comments At x = 0 and x=W = 0 adiabatic wall At z=0 =0 adiabatic wall At z=L =0 adiabatic wall At y= Hh ,-kh = - ks , Ts =Th fluid –solid interface At y = Hh + t - kc = - ks , Tc = Ts Fluid –solid interface By CFD code the above equations with their boundary conditions have been solved, and values of velocity, pressure, effectiveness, pumping power, friction factor and overall performance are calculated The inlet velocity is calculated based on Reynolds number. Re = (10) Where Dh is the hydraulic diameter. From Darcy equation can be friction factor [11] AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 11, No. 4 ISSN: 1998-4456 Page 434 Copyright  2018 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. f = uavg p 2  L D  2 (11) And the total pressure drop in heat exchanger unit is the summation of the pressure drop in the hot fluid channel and the cold fluid channel. ΔPu = ΔPch.h + ΔPch.c = (Ph.i – Ph.o ) + (Pc.i –Pc.o ) (12) Heat exchanger effectiveness is the ratio of the actual heat transfer to the maximum possible heat that can be transferred[10]: ) = ) = (13) For any heat exchanger the effectiveness is a function of NTU and the heat capacity ratio. ε = f (NTU, Cr) Where NTU is the number of transfer units NTU = And Cr is the heat capacity ratio Cr = (14) volumetric flow rate multiplying total drop pressure. [12] P.P = . V P (15) To calculate the performance factor of heat exchanger the performance factor has been used to correlate the hydrodynamic and thermal performance as define below  * = PP Q . (16) AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 11, No. 4 ISSN: 1998-4456 Page 435 Copyright  2018 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. 3.NUMERICAL SOLUTION: The above system of governing equations and boundary conditions are numerically solved using finite volume method. The flow is developing; the 3D continuity and 3D Naiver stock equations are solved numerical by fluent ANSYS . A mesh has been chosen in accepted size and a mesh refinement has been made to find out the appropriate mesh size that gives highly accurate solution as indicate in table (3).Table 3 shows that the solution becomes independent of grid size and from third mesh further increase in the grids will not have a significant effect on the solution and results of such arrangement are acceptable. Therefore the fourth mesh size is used for all calculations. Table 1 Grid independent study Mesh size Outlet temperature (K) mesh1 (number of elements=1320743) 308 mesh2 (number of elements=1528356) 308.93 mesh3 (number of elements=1963380) 308.95 mesh4 (number of elements=2575633) 308.96 The convergence criteria used to control the numerical solution for both momentum and energy equations are 10-6. The properties of material are list in table (4) below [13 ] Table (2) properties of selected material Property Material Water Aluminum Thermal conductivity (k) (W/m2.K) 0.643 273 Specific heat Cp (J / Kg .K) 4200 871 density(Kg / m3 ) 1000 2791 Viscosity (m2/s) 0.001 ـــــ AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 11, No. 4 ISSN: 1998-4456 Page 436 Copyright  2018 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. 4. RESULTS The model was solved with pure water as a working fluid with constant properties selected according to the mean temperature across heat exchanger. Two cases have been solved for each size with considering surface roughness and with neglecting it to explore the effect of roughness on the performance of (CFMCHX) To check the validity of the built numerical model, verification was made by solving the experimental model presented in [7] and the results were compared. The experimental model presented in [7] is a rectangular microchannel with aspect ratio 1.92 μm. The width (265  m), height (510  m), length (76 mm) and hydraulic diameter (348.8  m). The author used FC770 and deionized water used as working fluid. The surface roughness of aluminum measured was 3.7  m due to surface finish of machining. Fig. 3 shows the comparison between results of present numerical model and the experimental data of [7] for variation of friction factor with Reynolds number. From this figure it can be seen that, the agreement between numerical and experimental results is acceptable with the average error of 5% which may be due to difference in material properties and some assumptions in this model. To check the behaver of thermal side, another validation has been made with reference [14] as shown in fig. 4 shows the validation between results of present model with data in Al- bakhit and Fakeri 2006[14]. The figure explains that, the variation of dimensionless axial distance x* = with dimensionless mean temperature Tm = for the cold and hot fluids flows in rectangular microchannel heat exchanger at Reynold number is 100 and =100. It can be noted in this figure that there is a good agreement between the results of present model and that in [14] where the maximum error was 1% due to the difference in some assumptions between them. AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 11, No. 4 ISSN: 1998-4456 Page 437 Copyright  2018 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. Fig (5) shows the velocity variation on planes drawing a long the channels of (CFMCHX). The channels have hydraulic diameter (Dh = 50  m) and Re =100. It can observed in this figure that the velocity of flow in entering of cold and hot channels is maximum value and decreasing flow velocity with toward flow. Fig (3) variation of friction factor with Re as comparison between results of present model and that of [7] Fig(4) variation of mean temperature with dimensionless distance as comparison between results of present model and that of [14] AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 11, No. 4 ISSN: 1998-4456 Page 438 Copyright  2018 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. Fig.6 explains the distribution of pressure drop contour in (y-z) plane at same value of Re and hydraulic diameter of previous figure. Pressure drop in hot channel is larger than cold channel due to affect temperature which leads to increase it. Also it can be seen that the surface roughness increased the pressure drop due the restriction between molecules of fluid. Fig (5) planes of contour represent the velocity(m/sec) in cold and hot cells of CFMCHX at Dh = 50 μm and Re=100 Fig(6) pressure(pa) distribution contour along cold and hot channel of CFMCHX has Dh = 50 μm and Re = 100 AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 11, No. 4 ISSN: 1998-4456 Page 439 Copyright  2018 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. Fig.7 shows variation of pressure drop at many plane distributions a long hot and cold channel of (CFMCHX) with hydraulic diameter of ( Dh =50  m) and Re =100. It can be noted that in fig.7 the pressure drop increased with increasing velocity of flow which due to increase mass flow rate. Fig.8 shows contour of variation temperature distribution a long (CFMCHX). The hot fluid entering at Th =100 C and cold fluid Tc =50 C. As seen in this figure the temperature difference in cold channel increase toward fluid while in hot channel decreased due to heat transfer. Fig(7) contour show many planes of pressure (pa) along cold and hot cell of CFMCHX at Dh =50 μm and Re =100 Fig (8) Temperature(k) distribution contour along hot and cold cells of CFMCHX as Re =100 and Dh = 50 μm AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 11, No. 4 ISSN: 1998-4456 Page 440 Copyright  2018 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. Fig.9 shows the variation of temperature on contours on many planes distributed along the channels of (CFMCHX). As can be seen from this figure that the temperature increased with increase residence time and decreased the velocity of flow. Fig.10 shows the variation of effectiveness with Reynolds number at different values of hydraulic diameter (20,50,110,150) μm with considering roughness effect and neglecting it. It can be seen from this figure that the effectiveness decreased with increasing Reynold number for all diameters due to decreasing temperature difference with increase the velocity of flow. For all selected range of hydraulic diameters it can been noted in fig.10 decreasing the effectiveness with increasing values of hydraulic diameter because the area of heat transfer increasing with size. Also it can be observed that, the effect of surface roughness is small on the effectiveness which is due to independent of temperature on roughness as a result of constant properties assumption. Fig(9) many planes of Temperature(k) explaining as contour at Dh = 50 μm and Re =100 Fig (10) show variation of effectiveness ( ɛ ) with Re number with and without considering roughness effect. AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 11, No. 4 ISSN: 1998-4456 Page 441 Copyright  2018 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. Fig.11 shows the variation of changing in friction factor ( the difference between friction factor with roughness effect and friction factor without it ) with Reynolds number for selected range of hydraulic diameters (20,50,110,150) μm. From this figure it can be noted that, The change in friction factor decreased with increasing Reynolds number due to increase the velocity of flow according to Darcy equatio. Fig. 12 explains the relation between pressure drop difference and Reynolds number at different hydraulic diameters. The pressure drop difference in this figure represents the difference between value of pressure drop with considering surface roughness and with neglecting it. As can be observed in the fig.12 The pressure drop difference is maximum at small diameters due to decreasing the area of flow. This figure shows that, pressure drop difference is increased with increasing Reynolds number due to increasing the velocity of flow since the losses is maximized at higher velocity. The surface roughness has very important roll on increasing the value of pressure drop because of the resistance which is show against the fluid flow. Fig (11) explain variation of friction factor with Re for different values of hydraulic diameter AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 11, No. 4 ISSN: 1998-4456 Page 442 Copyright  2018 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. Fig.13 shows the variation of decrement of pumping power (the difference between the value of pumping power with roughness effect and its value with neglecting it) with Reynolds number through many values of hydraulic diameter. As can be seen in fig.13 the difference of pumping power increased with increasing the Reynolds number because of increasing the pressure drop and velocity of flow. Also it can be seen in this figure that, the pumping power is increased with decreasing hydraulic diameter due to high pressure generated in small diameter. The roughness effects work to increase pumping power value and hydraulic performance generally because of the resistance which causes it. Fig(12) variation pressure drop difference with Re for different value of hydraulic diameter Fig (13) variation of pumping power difference with Re for different hydraulic diameter AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 11, No. 4 ISSN: 1998-4456 Page 443 Copyright  2018 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. Fig.14 represent cooperation between performance factor( it is ratio of heat transfer rate to pumping power) and Reynolds number for two cases , one depending the surface roughness and another neglecting it. The fig (14) divided to (a) and (b) to explain clear variation with different values of hydraulic diameters. The fig (14a) plotted for diameters range (20,50 ) μm while (14 b) plotted for diameter values (110 , 150 ) μm. From these figures it can be seen that, the performance factor decreased with increasing Reynolds for all values of hydraulic diameter due to the increasing of pumping power very high comparison with increasing in convective heat transfer. The roughness effect becomes very clear on overall performance with larger hydraulic diameter. Also these figures reveal that, the effect of roughness increased at lower values of Re. Fig(14-a) variation performance factor with Re for two cases with and without considering roughness effect for hydraulic diameter (20,50) μm Fig(14-b) variation performance factor with Re for two cases with and without considering roughness effect for hydraulic diameter(110, 150) μm AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 11, No. 4 ISSN: 1998-4456 Page 444 Copyright  2018 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. CONCLUATIONS In this paper numerical investigation is made to study the effect of roughness on the performance of CFMCHX with different sizes and used water at constant properties as working fluid. From the results it can be conclude that: 1- The surface roughness has considerable effect on hydraulic performance and its effect is increased when hydraulic diameter and velocity of flow is small. 2- Surface roughness enhances the thermal performance slightly. 3- The roughness effect is decreased with increasing Reynold number and hydraulic diameter. 4- The surface roughness causes an extra increasing in pressure drop and pumping power. 5- In general the surface roughness leads to decrease the overall performance of CFMCHX. REFERENCES 1. Hasan, M. I. (2011). Influence Of Wall Axial Heat Conduction On The Forced Convection Heat Transfer In Rectangular Channels. Basrah Journal for Engineering Science, 1, 31-43. 2. Mala, G. M., & Li, D. (1999). Flow characteristics of water in microtubes. International journal of heat and fluid flow, 20(2), 142-148. 3. Weilin, Q., Mala, G. M., & Dongqing, L. (2000). Pressure-driven water flows in trapezoidal silicon microchannels. International journal of heat and mass transfer, 43(3), 353-364. 4. Qu, W., Mala, G. M., & Li, D. (2000). Heat transfer for water flow in trapezoidal silicon microchannels. International Journal of Heat and Mass Transfer, 43(21), 3925-3936. 5. Hasan, M. I., Rageb, A. A., Yaghoubi, M., & Homayoni, H. (2009). Influence of channel geometry on the performance of a counter flow microchannel heat exchanger. International Journal of Thermal Sciences, 48(8), 1607-1618. 6. Dang, T., Teng, J. T., & Chu, J. C. (2010). A study on the simulation and experiment of a microchannel counter-flow heat exchanger. Applied Thermal Engineering, 30(14-15), 2163-2172. 7. Kim, B. (2016). An experimental study on fully developed laminar flow and heat transfer in rectangular microchannels. International Journal of Heat and Fluid Flow, 62, 224-232. 8. Hasan, M. I. (2011). Numerical investigation of counter flow microchannel heat exchanger with MEPCM suspension. Applied Thermal Engineering, 31(6-7), 1068-1075. 9. Peiyi, W., & Little, W. A. (1983). Measurement of friction factors for the flow of gases in very fine channels used for microminiature Joule-Thomson refrigerators. Cryogenics, 23(5), 273-277. AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES Vol. 11, No. 4 ISSN: 1998-4456 Page 445 Copyright  2018 Al-Qadisiyah Journal For Engineering Sciences. All rights reserved. 10. Hassan, M. I. Numerical investigation of counter flow microchannel heat exchanger with different geometries working fluid , PhD thesis, Mechanical engineering department, college of engineers, Basrah university. 11. Qi, S. L., Zhang, P., Wang, R. Z., & Xu, L. X. (2007). Single-phase pressure drop and heat transfer characteristics of turbulent liquid nitrogen flow in micro-tubes. International journal of heat and mass transfer, 50(9-10), 1993-2001. 12. Hasan, M. I., Rageb, A. M. A., & Yaghoubi, M. (2012). Investigation of a counter flow microchannel heat exchanger performance with using nanofluid as a coolant. Journal of Electronics Cooling and Thermal Control, 2(2), 35-43. 13. Shakir, A. M., Mohammed, A. K., & Hasan, M. I. (2011). Numerical investigation of counter flow microchannel heat exchanger with slip flow heat transfer. International Journal of Thermal Sciences, 50(11), 2132-2140. 14. Al-Bakhit, H., & Fakheri, A. (2006). Numerical simulation of heat transfer in simultaneously developing flows in parallel rectangular ducts. Applied thermal engineering, 26(5-6), 596-603.