Article AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES 14 (2021) 066–074 Contents lists available at http://qu.edu.iq Al-Qadisiyah Journal for Engineering Sciences Journal homepage: http://qu.edu.iq/journaleng/index.php/JQES * Corresponding author. E-mail address: donia.ali@utq.edu.iq (Dunya A. Mohammad) https://doi.org/10.30772/qjes.v14i1.736 2411-7773/© 2021 University of Al-Qadisiyah. All rights reserved. This work is licensed under a Creative Commons Attribution 4.0 International License. Numerical investigation of the electric double-layer effect on the performance of microchannel heat exchanger at combined electroosmotic and pressure-driven flow Dunya A. Mohammada*, Mushtaq I. Hasana, Ahmed J. Shkaraha aMechanical Engineering Department, College of Engineering, University of Thi-Qar –Iraq. A R T I C L E I N F O Article history: Received 4 April 2021 Received in revised form 8 May 2021 Accepted 23 May 2021 Keywords: microchannel heat exchanger (MCHE), electric double layer (EDL) thickness, electroosmotic flow, Poisson-Boltzmann equation. A B S T R A C T Numerically investigated the electric double layer (EDL) Effects on the performance of the square microchannel heat exchanger (MCHE) at combined electro-osmotic and pressure-driven flow with compared pure pressure-driven with a hydraulic diameter (10 – 50) μm. We defined at any size (Dh) of microchannel heat exchanger become the impact of EDL very slight with the studied effect of electric double layer thickness λ. The diluted water 1:1 potassium chloride (KCl) solution is used as a working fluid at an ionic concentration (10−4, 10−6) M, silicon microchannel at zeta potential of surface -0.2 volt. A three-dimensional (3D) Poisson-Boltzmann equations and Naiver-stoke equations with applied electric field solved by using the finite volume scheme in this work. The results show an increase in pressure drop of the microchannel heat exchanger at combined flow electroosmotic and pressure-driven flow with a percentage of 31.09 % at an ionic concentration 10−4 M and 42.71 % at 10−6 M, increase in pumping power, especially at low ionic concentration. Slight enhancement in average heat transfer rate and effectiveness due to an increase in average temperature difference. Decrease in overall performance at combined electroosmotic and pressure-driven flow compared with pure pressure driven. © 2021 University of Al-Qadisiyah. All rights reserved. 1. Introduction In recent decades, rapid progress is observed in the manufacturing of microsystems (have a length between 1 μm – 1 mm) especially in micro- electronic mechanic’s system (MEMS) which integrates mechanical and electrical systems. The development in MEMS caused an augmentation in understanding the structure and flow of fluid. Manufacture of microelectronic mechanic’s system (MEMS) by using different techniques such as bulk silicon micromachining, surface silicon micromachining, lithography, electro-discharge machining (EDM) and plastic modeling … etc. Mushtaq I. Hasan et al. [1]. Basic phenomena of electro-kinetic in the microchannel can be classified into electro-osmosis, Electrophoresis, electro-migration, and Streaming potential. The process of transporting fluid in a microchannel subjected to an electric field is called electro- osmotic flow or pressure gradient is called pressure-driven flow or transport of fluid by an applied electric field and pressure-driven together this called combined electroosmotic and pressure-driven flow Berli [2]. A number of studies investigated the effect of EDL at pressure-driven or combined electroosmotic and pressure-driven flow. For instance, Chun and Li (1997) [3] numerically studied the impact of the electrical double layer (EDL) on the flow of liquid, pressure distribution, streaming potential, and friction coefficient for rectangular microchannel. 2D Poisson-Boltzmann equation and the exact solution to the motion equation were solved by recruiting the green function formulation. They found that high zeta potential and law value of ionic concentration lead to a great effect of electric double layer field on the flow of liquid. The increase in http://qu.edu.iq/ mailto:donia.ali@utq.edu.iq https://doi.org/10.30772/qjes.v13i http://creativecommons.org/licenses/by/4.0/ DUNYA A. MOHAMMAD, MUSHTAQ I. HUSAN AND AHMED J. SHKARAH /AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES 14 (2021) 066–074 67 pressure drop caused increased induced electro-kinetic potential, while the increase in the ionic concentration of solution leads to a decrease in induced electro-kinetic. The EDL effect caused an increment in friction coefficient and viscosity apparent compared with no electric double layer effect. Prashanta D. and Ali B. (2001) [4] conducted an analytical solution of steady Newtonian flow through a 2D straight channel with combined pressure-driven flow PDF and electro-osmotic flow EOF at constant zeta potential of wall and constant buffer concentration. Pressure gradient, mass flow rate, vorticity in combined PDF/EOF were investigated. The shear stress of the wall depended on three parameters (effective electric double layer thickness, EDL vorticity thickness, and EDL displacement thickness. Their results indicate to (1) the distribution of electro-osmotic potential Ψ* is a relationship (function) of α (the ionic energy parameter) only. (2) They defined an effective electric double layer thickness as a function of α (ionic energy parameter) and clarify that the EDL effects are confined to a region x≤ EDL thickness. (3) Defect in mass flow rate because of the distribution of velocity within EDL, this defect depended on EDL displacement. (4) For combined EOF/PDF, they explained that using Helmholtz velocity as an appropriate condition between the electric double layer is completed and bulk flow. Ng and Tan, (2006) [5] Numerically investigated a single rectangular microchannel with electrical double layer (EDL) impact. Two models Poisson–Boltzmann model (PBM) and Nernst–Planck model (NPM) additional to Navier–Stock equations were used to analyze the EDL effect in the channel. In comparison between the two models, it was found that the Poisson- Boltzmann model is better than the Nernst-Planck model in terms of the cost (RAM and CPU). In addition, the P-B-M showed less contradiction of modeling electric double layer a far about surface of the microchannel. They concluded that reducing friction confidence at an increase the Schmidt number (Sc), the electric double layer (EDL) effects lead to an increase in the Schmidt number and hence the decrease in friction coefficient. Thus, the Poisson–Boltzmann model(PBM) is still an attractive model to compute the electrical double layer (EDL) effect in the microchannel. Dayong Y. (2011) [6] analytical solution of velocity distribution and potential for combined pressure-driven flow and electroosmotic flow in microchannel. The mathematical model to simulate fluid flow consisted of the Navier Stokes and the Poisson Boltzmann equations, which were solved by using the finite element method in Matlab software. Their results showed that at combined EOF/PDF the distribution of velocity is parabolic and compound of plug-like when the fluid is steady-state, the distribution of velocity is observed to be plug-like and similar to the EDL potential profile even in the case of pure EOF. Their results supply the guidelines to applications of combined flow (EOF/PDF) in microchannel chips. Nandy et al. (2013) [7] experimental study to improve the heat transfer performance by decreasing hydraulic diameter size or by using better thermal conductivity for working fluid in microchannel heat exchanger addition to studying the testing and design of microchannel heat exchanger (MCHE). Working fluids are distilled water, water nanofluid (Al2O3) at volume concentration 1%, 3%, and 5%, water nanofluids (SnO2) at volume concentration 1%. They concluded to enhance the heat transfer rate add Nanoparticle in the fluid, water 50% (Al2O3) and water-nano fluids 1% (SnO2) greater absorption the heat with 9% and 12% comparison with base fluid also higher overall heat transfer coefficient than the base fluid. Mushtaq et. al. (2017)[8] applied numerical investigation of two electrolyte solutions are water and PBS (Phosphate- Buffered-saline solution) through the square microchannel by solved Navier stokes equation, Laplace and Poisson equation with the pure electro-osmotic flow. It was found that the water solution gives lower flow and velocity compared with PBS, and it was noticed a very small increase in temperature because of the simple effect of Joule heating. Moreover, it was concluded that zeta potential and electric field have a great effect on flow rate. Mushtaq et. al. (2017) [9] studied the pure electro-osmotic flow through a microchannel with contractions to contribute to an increase in the biomolecules flow corresponding to the used different shapes (rectangular, trapezoidal, curved, and triangular) of contractions. COSMOL software to analyze the flow behavior was used. Their results showed that the distribution of velocity for triangular contraction has higher values of flow velocity, and highest compared to other geometries of contractions. In addition, it was concluded that the velocity of flow increased with the increasing of zeta potential and electric field. Ahmed A. et al. (2019) [10] numerically investigated of scaling effect (surface roughness effects) of parallel flow in the square microchannel heat exchanger on thermal and hydrodynamic performance at various hydraulic diameters. It was concluded that the impact of surface roughness leads to improve thermal performance and increase in pressure difference. The effect of roughness reduced with increment Re and hydraulic diameter. Ahmed A. et al. (2019) [11] numerically studied of scaling effect (slip flow effects) in counter flow MCHE on performance at different geometries (square, trapezoidal, triangle, and circular). In addition, changing the size of a channel on thermal and hydrodynamic performance were studied. Their results indicated that the trapezoidal channel has higher overall performance compared to other shapes. Moreover, triangle shape has the higher effect on slip flow and hence the square shape. It was found that greater slip flow is related to small Nomenclature A Cross-section area (𝑚2) Greek symbols Cp Specific heat (J/kg.K) p Density Dh Hydraulic diameter (µm) µ Viscosity H Height (µm) ζ Zeta potential K Thermal conductivity (J/k) ψ Electric potential L Length channel of heat exchanger (µm) κ Debye-Huckel parameter M Mass flow rate (Kg/s) λ EDL thickness P Pressure (Pa) ε Permittivity P.P Pumping power (watt) Subscripts Q Heat transfer rate (watt) EOF Electro-osmotic flow Re Reynold number PDF Pressure driven flow T Temperature (K) EDL Electric double layer t Separating wall thickness (µm) i inlet U Velocity in x-direction (m/s) o outlet V Velocity in y-direction (m/s) c Cold fluid w Velocity in z-direction (m/s) h Hot fluid 68 DUNYA A. MOHAMMAD, MUSHTAQ I. HUSAN AND AHMED J. SHKARAH /AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES 14 (2021) 066–074 velocity and hydraulic diameter. Thus, the slip flow causes an increase in pressure difference and pumping power. Qingkai et al. (2020) [12] theoretically investigated heat transfer and flow of nano fluid in a horizontal microchannel in combination with the effect of electric double layer and magnetic field. momentum equation with the presence of EDL effect were solved . In addition, nanoparticle volume fraction, temperature, velocity, the Brinkman number (Br), Hartmann number Ha, the impact of parameter k, and various physical quantities were discussed thoroughly. The finding indicates that the electric double layer and magnetic field may be used to control the heat transfer and flow in the microchannel. The heat improvement depends on the temperature applied on the wall and the Brinkman number. The magnetic field impact on friction coefficient and nanoparticle volume fraction can disregard the impact of magnetic field on Sherwood number and Nusselt number. Ahamed C. et. al. (2020) [13] numerically studied of electro-osmotic /pressure-driven flow with triangular obstacles (block) through rectangular microchannel by alternative technique CFD and modified IBM to solve Navier stokes equation with Poisson Boltzmann equation and Nernst Plank equation. They used hybrid technique to improve mixing efficiency by the Passive and active integration ways. Their findings refer to (a) The mixing efficiency depended on a change in several triangular blocks, the zeta potential of surface, and EDL thickness (b) Mixing efficiency reduced with an increase in an applied electric field, Reynold number, and Peclet number. (c) Mixing efficiency is higher from 28.2 – 50.2 % at increment the number of a block from (1- 5). 2. Problem description: In this paper, a three-dimension parallel flow square microchannel heat exchanger (MCHE) as shown in Fig. 1 is modeled. study one unit of microchannel heat exchanger consists of hot and cold fluid as shown in Fig. 2. The length of the heat exchanger is L=1mm and its hydraulic diameter (10-50µm). The thickness of the wall between the hot channel and cold channel equal t=3µm, while the ratio between the velocity of pressure-driven to the velocity of electro-osmotic flow is 0.96. the temperature of inlet hot and cold are 373 k, 293 k receptivity. The zeta potential of the channel surface (-0.2 volt), and the applied field electric Ez on working fluid ( 𝟏𝟎𝟔) volt/m. (d) The electric double layer overlaps effect on mixing efficiency. Their findings are useful in a lot of fields (Biomedical, cooling of microchips, deoxyribonucleic acid hybridization and biotechnological … etc. Most of the previous papers theoretical and practical investigated heat transfer and flow characteristics in the single microchannel. This study investigates the effect of electric double layer (EDL) at two values of ionic concentration through parallel flow microchannel heat exchanger (PMCHE) and determines at any size of microchannel heat exchanger stops EDL effect or is the slight effect. Compared the results with no electric double layer (EDL) effect (pure pressure-driven). Figure .1 Schematic model of microchannel heat exchanger Figure . 2 A unit of heat exchange consists of the hot and cold channel [14]. 2.1. Electric Double Layer (EDL) The electrical potential surface for a solid wall of the microchannel and the liquid has quantities of ions, this leads to attraction of the counter ions in dilute liquid to electrostatic charges of surface to create an electric field. The order of electrostatic charge and ions of liquid is called electric double layer (EDL) as shown in Fig. 3. There are two types of ions in term of motion, for the compact layer (the layer near of surface) the ions are immobile, and in the diffuse double layer (DDL) the electric field less influence on the ions (mobile). At the dilute liquid flow through the microchannel, the mobile ions of the electric double layer create an electric current (streaming current) to flow with liquid flow. In case the streaming potential of electro-osmotic is very small compared to the applied external electric field. The gathering of the ions in direction of flow sets up an electric potential and electric field together known streaming potential, this caused in creates a current is called conduction current. The conduction current direction is counter to the direction of liquid flow. The maximum value of the thickness of the electric double layer (EDL) 1µm depending on the properties of the liquid are (a) ionic concentration (b) temperature of the dilute liquid, and (c) zeta potential of the surface. To show characterize of the electric double layer (EDL) effects used the Debye-Huckel parameter k [13]. 1 𝑘 represent the thickness of the electric double layer (EDL) Kumar [15]. DUNYA A. MOHAMMAD, MUSHTAQ I. HUSAN AND AHMED J. SHKARAH /AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES 14 (2021) 066–074 69 Figure . 3 Scheme formed of electric of double layer on the surface wall [16] 3. Mathematical formulation Governing Equations 3D steady-state, incompressible and laminar flow, the following equations are solved to calculate the distributions of velocity, temperature, and EDL distribution for parallel flow microchannel heat exchanger. Poisson’s equation According to the theory of electrostatics, the relationship between Ψ and ρe is given by the Poisson’s equation Sheikhizad et al. [17], which for a rectangular channel 𝜕2𝛹 𝜕𝑥2 + 𝜕2𝛹 𝜕𝑦2 + 𝜕2𝛹 𝜕𝑧2 = - 𝜌𝑒 𝜀𝜀ₒ (1) ε is the dielectric constant of the medium εₒ is the electric permittivity of a vacuum. 𝜌𝑒 is the net volume charge density, we have? ni=nₒi 𝑒𝑥𝑝(− 𝑧𝑖𝑒𝛹 𝑘𝑏𝑇 ) (2) ρ𝑒 =(𝑛+ − 𝑛−) = −2𝑛ₒ𝑧𝑒𝑠𝑖𝑛ℎ( 𝑧𝑒𝛹 𝑘𝑏𝑇 ) (3) 𝑛+ 𝑎𝑛𝑑 𝑛− are a concentration of cations and anions, respectively Kb Boltzmann’s constant =1.3805 ∗ 10−23𝐽𝑚𝑜𝑙−1𝐾−1 e electron charge = 1.6021 ∗ 10−19𝐶 nₒi bulk concentration Poisson-Boltzmann eq. (1) become ∂2Ψ ∂x2 + ∂2Ψ ∂y2 + ∂2Ψ ∂z2 = 2nₒze εεₒ sinh zeΨ KbT Continuity equation: [5] ∂u ∂x + ∂v ∂y + ∂w ∂z = 0 (4) Momentum equation x-direction u ∂u ∂x + v ∂v ∂y + w ∂w ∂z = μ ρ ( ∂2u ∂x2 + ∂2u ∂y2 + ∂2u ∂z2 ) − μ ρ ∂p ∂x (5) y- direction u ∂v ∂x + ∂v ∂y + ∂v ∂z = μ ρ ( ∂2v ∂x2 + ∂2v ∂y2 + ∂2v ∂z2 ) − μ ρ ∂p ∂y (6) The applied electric field (𝐸𝑧 𝜌ₑ ) u ∂w ∂x + v ∂w ∂y + w ∂w ∂z = μ ρ ( ∂2w ∂x2 + ∂2w ∂y2 + ∂2w ∂ z2 ) − 1 ρ ∂p ∂z − Ezρₑ (7) Energy equation: u ∂T ∂u + v ∂T ∂v + w ∂T ∂w = k ρCp ( ∂2T ∂x2 + ∂2T ∂y2 + ∂2T ∂z2 ) (8) The heat transfer rate of a microchannel heat exchanger between hot and cold [1]: 𝑄 = 𝑚𝐶𝑝∆𝑇 (9) The effectiveness (ε)= Qact Qmax = Ch(Thi−Tho) Cmin (Thi−Tic) = Cc(Tco−Tci) Cmin(Thi−Tci) (10) The overall performance of MCHE can be calculated from eq. [1] η = ε ∆Pt (11) ∆Pt=∆Pc+∆Ph=(Pin,c-Po,c)+(Pin,h-Po,h) (12) The boundary conditions [17] At inlet of the channel Z = 0, then U = 0 V = 0 W= Uin & T = Tin. At the outlet of the channel Z = Lch and the flow is fully developed, then 𝜕𝑢 𝜕𝑧 = 0 𝜕𝑣 𝜕𝑍 = 0 𝜕𝑤 𝜕𝑧 = 0 & 𝜕𝑇 𝜕𝑧 = 0 For three cold walls adiabatic (sides walls and up the wall) and three hot walls adiabatic (sides walls and down wall) conditions are applied: 𝑢 = 0 𝑣 = 0 𝑤 = 0 𝜕𝑇 𝜕𝑧 = 0 At x=0, 𝛹 = 𝜁ₒ x=w, 𝛹 = 𝜁ₒ At y=0, 𝛹 = 𝜁ₒ At y=h, 𝛹 = 𝜁ₒ For the up wall of the hot channel and down a wall of the cold channel 𝛹 = 𝜁ₒ 4. Properties of material: Water is a working fluid with aqueous KCL solution at the ionic concentration ( 10−4 M, 10−6 M ), vacuum permittivity (F/m) 8.85 ∗ 10−12 , dielectric constant ε 80 and 6.39 ∗ 10−10 The permittivity of water (C/V*m). The data of the fluid (water) input by creating a User Defined Function (UDF). 70 DUNYA A. MOHAMMAD, MUSHTAQ I. HUSAN AND AHMED J. SHKARAH /AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES 14 (2021) 066–074 Table 1. shows the properties of dilute liquid and silicon. material Density (𝐾𝑔 𝑚3⁄ ) Cp (𝐽 𝐾𝑔. 𝐾)⁄ K (𝑊 𝑚. 𝐾)⁄ µ (𝐾𝑔 𝑚. 𝑠)⁄ Dilute liquid 997 4182 0.605 0.9e-03 Silicon 2329 700 148 - 5. Numerical solution: The system of governing equations and the Poisson-Boltzmann equation and boundary conditions are numerically solved using the finite volume method (FVM). SIMPLE algorithm is used to solve the problem of velocity-pressure distribution and UDF to solve Poisson-Boltzmann equation to study effect electric double layer (EDL). A Fluent19 software has been used to do the numerical solution. Table 2 shows independent grid size for hydraulic diameter (20μm) and effect on the solution results. Table 2. mesh independent Mesh size Outlet temperature (K) of hot channel Outlet temperature (K) of cold channel Mesh 1 (size element =5μm) 318.253 347.748 Mesh 2 (size element =3μm) 318.287 347.573 Mesh 3 (size element =1μμm) 318.204 347.302 Mesh 4 (size element = 0.8μm) 318.201 347.30 6. Results and discussions Numerical simulations of heat transfer and fluid flow problem were executed by using the commercial computational fluid dynamic (CFD) package. The governing equations for liquid flow were solved by finite volume method (FVM), and the SIMPLE algorithm has been used to solve the problem of velocity-pressure coupling. User-Defined Function (UDF) has been employed to solve the Poisson-Boltzmann equation to study the combination EOF/PDE flow through MCHE. The thickness of the electric double layer plays a significant role in combined electro- osmotic / pressure-driven flow due to its effect on the flow, temperature difference, and total pressure drop. The thickness of the electric double layer can be estimated from Debye Huckel parameter κ, such refer 1/κ to EDL thickness λ. The relation between EDL thickness and ionic concentration M is reversed which the increase in the ionic concentration of solution leads to a decrease in EDL thickness (λ). The parameters will be used in the investigation shown in Table 3. The parameter using in this paper are shown in Table 4. Table 3. the parameters using in this paper. The name Value the dielectric constant of the fluid 80 permittivity of vacuum 8.85 * 10−12 C/V*m The permittivity of fluid (water) 6.93*10−10 C/V*m Electric field 106 volt/m Zeta potential -200 mV Ionic concentration 10−4 & 10−6 M Charge of electron -1.602*10−19 C Valence (z) 1 Thickness of EDL 32 nm, 324 nm Table 4 shows the percentage change (pure PDF and combined PDF /EOF) for square MCHE of the variation parameters with the hydraulic diameter at two values of concentration. The percentage change of pressure drop calculated from the following equation: The percentage change= (∆P𝐸𝑂𝐹&𝑃𝐷𝐹−∆𝑃𝑃𝐷𝐹) ∆𝑃𝐸𝑂𝐹&𝑃𝐷𝐹 The other parameter calculated with the same method. Table 4. average percentage change of parameters parameters The percentage change % at c=10−4 The percentage change % at c=10−6 Total pressure drop 31.09 42.71 effectiveness 19.26 35.11 Overall performance 20.46 6.97 To validate the work of user-defined function (UDF) that was used to solve Poisson Boltzmann equation and applied electric field compared the present model with pure electro-osmotic in microchannel [8]. The model presented in [8] is a microchannel that has a hydraulic diameter of 100 μm, the width of channel 100 μm, the height of microchannel 100 μm, and the length 500 μm. They used a glass microchannel with zeta potential -0.1 volt and the water as a working fluid with ionic concentration 3.727 ∗ 10−6 M and applied electric field 100000 volt/m. The agreement between a present model and that for [8] is very good, and the percentage of the error 2 %. Fig. 4 shows the comparison between results of the data of [8] for the velocity of electro- osmotic with a width of microchannel and result of the present model. Figure .4 For the present model and that for [8]. To check the numerical model validity, verification was made by solving the model presented in [18] a compared the results. The model presented in [18] is a square microchannel at combined flow (electro- osmosis flow and pressure-driven flow) with a hydraulic diameter of 25 μm, channel height 25 μm, channel width 25 μm, and length 100 μm with ∆P/dx = 105 pa/m , an inlet temperature of 289 K, and molar concentration 10−4 M, 10−5 M, zeta potential ζ = 150, 200 mV. Fig. 5 shows the comparison between results of the data of [18] for velocity distribution with a width of microchannel at different ionic concentrations. It can be seen that from this Fig. 5 there is a good agreement between the 0 0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008 0.009 0 10 20 30 40 50 60 70 80 90 100 110 v e lo ci ty o f e le ct ro o sm o ti c width of microchannel (cm) mushtaq [8] present model DUNYA A. MOHAMMAD, MUSHTAQ I. HUSAN AND AHMED J. SHKARAH /AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES 14 (2021) 066–074 71 results of the present model and that for [18] and the error between the model presented and a result of [18] less than 5 %. Figure .5 For the present model and that for [18]. To show the distribution of velocity, pressure, and temperature using contours on created planes with length and width of the microchannel heat exchanger. Fig.6 explains the velocity distribution along the microchannel of PFMCHE (y-z) plane at the small size of channel Dh=10 μm due to the ratio of EDL thickness to channel height is larger especially when higher zeta potential (-200 volt) and low ionic concentration. It can be noticed from this figure that the flow velocity in a center is maximum and minimum toward the wall as a result of the friction and the EDL field is opposite of the flow direction, this causes a disabled to flow. Figure .6 contour explains the distribution of velocity (m/sec) at PFMCHE on (y-z) plane at Dh=10 μm and Re=30. Fig. 7 represents the distribution of flow velocity for many planes (x-y) along hot and cold channels of parallel flow MCHE when Dh=10 μm and Re=50. From Fig. 6 it is observed that the velocity at the center of the channel is higher and decreasing toward walls due to friction effects and the EDL field. Maximum velocity at the inlet flow and minimum value of the velocity at the outlet channel is observed due to losses. Figure .7 contour represents the variation of velocity (m/s) along the hot and cold channels at parallel flow MCHE, Dh=10 𝛍m, and Re=30. Fig. 8 shows (z-y) plane for the pressure drop distribution contour at the same value of hydraulic diameter and Re along cold and hot channels of PFMCHE. It can be concluded from the figure that, the pressure drop reduced with the flow as a result of losses due to friction and dynamic friction. Figure .8 contour pressure distribution along the cold channel and hot channel of parallel flow MCHE for Dh=10 𝛍m and Re=30. Fig. 9 displays the distribution of temperature contours along the length of the channel with a parallel flow microchannel heat exchanger. The inlet fluid of hot channel 373 K, and inlet fluid of cold channel 293 K at hydraulic diameter 10 μm. It can be seen from Fig. 9 that, there is an increment in heat transfer and temperature difference with the dilute liquid flow, and reach a high value in the channel end due to heat exchange from the hot liquid to cold liquid. 0 5 10 15 20 25 30 0 0.002 0.004 0.006 0.008 0.01 w id th o f ch a n n e l μ m velocity m/s result of present model at 10 E -05 result of [18] 10 E -05 result of present model at 10 E-04 result of [18] 10 E-04 72 DUNYA A. MOHAMMAD, MUSHTAQ I. HUSAN AND AHMED J. SHKARAH /AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES 14 (2021) 066–074 Figure .9 contour of temperature distribution (K) for many planes of parallel flow MCHE Dh=10 𝛍m and Re =50. Fig. 10 explains the variation of the total pressure difference ( ∆𝑃𝑡) with hydraulic diameter for unit parallel microchannel heat exchanger at two cases of driven flow (electro-osmotic/pure pressure-driven flow) with two value of ionic concentration (10−6, 10−4) M. From this figure, it can be observed that the total pressure drop declined with the increase of hydraulic diameter for two cases as a result of an increment in the cross- sectional area of a MCHE. In addition, it can be shown that higher pressure drop at combined EOF/PDF comparison with pure pressure- driven due to increase in velocity of liquid flow which caused larger losses in pressure. A low concentration leads to higher pressure drop due to an increase in EDL thickness at low concentration can be seen from Fig. 10 that a low concentration leads to higher pressure drop due to an increase in EDL thickness at low concentration, with notice that the impact of EDL lower with increasing in hydraulic diameter depended on the ratio between the thickness of EDL to the height of the channel, at Dh > 40 μm begin the effect of EDL stop due to the large height of microchannel compared with the thickness of EDL. Figure. 10 variation of total pressure drop (KPa) with a hydraulic diameter (μm). Fig. 11 shows the temperature difference variation ( ∆𝑇) for pure PDF and combined flow with a hydraulic diameter at two values of ionic concentrations (two values of EDL thickness). From the figure, it can be concluded that the temperature difference decrement with increasing hydraulic diameter for all cases. Can be noticed from Fig. 11 that, increase in average temperature difference at combined flow compared with PDF especially at higher EDL thickness as a result of the flow of electric double layer field is the opposite direction to the flow in the channel center, this lead to obstructing inflow at the surface of a channel at combined flow (EOF/PDF). Figure .11 variation of average temperature difference with the hydraulic diameter (μm). The relation between the effectiveness (ε) and hydraulic diameter is shown in Fig. 12. It can be observed that the effectiveness is decreased with the increase of hydraulic diameter due to a decrease in temperature difference with the increased size of the channel. Moreover, the effectiveness in the case of combined flow is higher than that for pure PDF and the effectiveness increased with decreasing the concentration. The increase in effectiveness in presence of EDL is due to an increase in temperature difference when mixed EOF/PDF and the slight effect on effectiveness Figure .12 Variation of effectiveness with a hydraulic diameter (μm). 0 100 200 300 400 500 600 700 800 900 1000 5 10 15 20 25 30 35 40 45 50 55 ∆ P ( K p a ) Dh (μm) PD EOF/PDF, C=10E-4 EOF/PDF, C=10E-6 0 5 10 15 20 25 30 35 40 5 10 15 20 25 30 35 40 45 50 55 ∆ T ( K ) Dh PD EOF/PDF, C=10^4 EOF/PDF, c=10^-6 0 0.1 0.2 0.3 0.4 0.5 0.6 5 10 15 20 25 30 35 40 45 50 55 ε Dh PD EOF/PDF, C=10E-4 EOF/PDF, C=10E-6 DUNYA A. MOHAMMAD, MUSHTAQ I. HUSAN AND AHMED J. SHKARAH /AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES 14 (2021) 066–074 73 Fig. 13 represents the variation of the performance index ( 𝜂) with a size of a microchannel. It can be seen the performance index increased with increasing hydraulic diameter for all cases due to the decrease in total pressure with increased channel size. Furthermore, the impact of electric double layer thickness has a significant impact on overall performance, which leads to a decrease in performance index due to the increase in total pressure difference at combined driven. The performance index is observed to be decreased in small size of microchannel compared with pure pressure-driven and improvement in performance at large size of microchannel heat exchanger because lower the ratio between a thickness of EDL to the height of microchannel. Figure. 13 Variation of overall performance with a hydraulic diameter(𝛍m). The relationship between an average heat transfer rate ( 𝑄𝑎𝑣𝑔. ) , hydraulic diameter for unit microchannel heat exchanger for two cases pure pressure-driven and mixed driven flow is shown in Fig. 14. It can be concluded that the heat transfer rate is increased with the increasing of hydraulic diameter for all cases as a result of an increment in the cross- sectional area of a MCHE. In this figure, the average heat transfer due to the lost heat is not equal to the heat gain a combined PDF/EOF due to the formation of the EDL, which means an increase in thickness of the separator wall. Fig. 14 shows Enhancement in the average heat transfer rate in combined PDF/EOF at higher EDL thickness (low ionic concentration) due to an increase in the lost heat of the hot channel. Figure .14 variation of average heat transfer (watt) with a hydraulic diameter Fig. 15 shows the variation of the pumping power ( P. P) with hydraulic diameter. In addition, it can be observed that the relation of pumping power with channel size is steady as a horizontal line for all cases. This means the change in hydraulic diameter has no impact on pumping power. Can be seen from Fig. 15 that an increase in the pumping power at combined flow due to an increase in total pressure drop and flow rate of liquid and higher pumping power at low ionic concentration. Figure .15 variation of pumping power (watt) with a hydraulic diameter. Fig. 16 expresses the variation of the performance factor ( η ∗) with hydraulic diameter. It can be seen from this figure that the performance factor increment with increasing hydraulic diameter due to an increase in the area of the microchannel. Moreover, it shows an increase in the heat exchanged with remained pumping power at a constant value with changed channel size. For mixed PDF/EOF flow, it can be seen that the relationship between performance factor and size of the microchannel is linear, also notice that reduced in performance factor especially when the small size of microchannel due to the higher effect of EDL with enhancement in performance factor at a large hydraulic diameter and low concentration. Figure .16 variation of overall performance with a hydraulic diameter 0 0.000001 0.000002 0.000003 0.000004 0.000005 0.000006 0.000007 5 10 15 20 25 30 35 40 45 50 55 η Dh PD EOF/PDF, C=10E-4 EOF/PDF, C=10E-6 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 5 10 15 20 25 30 35 40 45 50 55 Q ( w a tt ) Dh PD EOF/PDF, C=10E-4 EOF/PDF, C=10E-6 0 0.00001 0.00002 0.00003 0.00004 0.00005 0.00006 0.00007 0.00008 0.00009 0.0001 5 10 15 20 25 30 35 40 45 50 55 P .P Dh PD EOF/PDF, C=10E-4 EOF/PDF, C=10E-6 0 500 1000 1500 2000 2500 5 10 15 20 25 30 35 40 45 50 55 η * Dh PD EOF/PDF, C=10E-4 EOF/PDF, C=10E-6 74 DUNYA A. MOHAMMAD, MUSHTAQ I. HUSAN AND AHMED J. SHKARAH /AL-QADISIYAH JOURNAL FOR ENGINEERING SCIENCES 14 (2021) 066–074 Fig. 17 explains the variation of a temperature difference as a comparison between hot and cold channels with a hydraulic diameter at two values of ionic concentration. From the figure, it can be noticed that for pressure- driven flow temperature difference is similar for the two-channel (hot and cold). However, the temperature difference for a hot channel is higher than cold channel especially for low ionic concentration depending on EDL thickness. The reason is attributed to the electric double layer thickness of an increase with an increase in temperature inlet of liquid. The consequences, leads to make the loss of heat is not equal to the heat gain in the case of presence the electric double layer depended on ionic concentration and temperature inlet of liquid. Figure .17 variation of temperature different for hot and cold channel with a hydraulic diameter (μm). Conclusions: A numerical solution to study the effect of the electrical double layer (EDL) on the thermal and hydraulic performance of parallel flow microchannel heat exchanger by solving 3D Navier-stokes and 3D Poisson-Boltzmann equation was conducted. The results can be concluded that: 1. Increasing in total pressure drop at combined flow (electro- osmotic and pressure-driven flow) for parallel flow microchannel heat exchanger especially at a higher thickness of the electric double layer (low ionic concentration and higher zeta potential. 2. Notice the effects of EDL are very slight at large size of hydraulic diameter Dh > 40 μm depending on the ratio of a thickness of electric double layer to the height of microchannel. 3. Enhancement in average heat transfer rate and effectiveness. 4. Higher in pumping power at combined driven flow. 5. The thickness of the electric double layer has an important role in the effect on the flow rate and performance of the microchannel heat exchanger. 6. For pure pressure driven, the temperature difference for hot channel is equal to the temperature difference for cold channel (the heat gain is equal to heat loss). 7. For combined flow (PDF/EOF), the temperature difference for hot channel is higher than cold channel. REFERENCES [1] Mushtaq I. Hasan, A. Rageb, M. Yaghoubi, and Homayon Homayoni, Influence of channel geometry on the performance of a counter flow microchannel heat exchanger, International Journal of Thermal Sciences, 48 (2009) 1607–1618 [2] Claudio Berli, Theoretical modelling of electrokinetic flow in microchannel networks, Colloids and Surfaces A: Physicochemical. Eng. Aspects 301 (2007) 271– 280l. [3] Chun Yang and Dongqing Li, Electrokinetic effects on the liquid flow in rectangular microchannels, Journal of Colloid and Interface Science, 149 (1997) 95- 107. [4] Prashanta Dutta and Ali Beskok, Analytical Solution of Combined Electroosmotic/ Pressure Driven Flows in Two-Dimensional Straight Channels: Finite Debye Layer Effects, Analytical Chemistry J., 73 (2001) 1979-1986. [5] Ng and Tan, Numerical analysis of EDL effect on heat transfer characteristic of 3-D developing flow in a microchannel, Numerical Heat Transfer, Part A, 49 (2006) 991–1007. [6] Dayong Yang, Analytical solution of mixed electroosmotic and pressure-driven flow in rectangular microchannels, Key Engineering Materials , 483 (2011) 679-683. [7] Nandy Putra1, Wayan Nata Septiadi, Gerry Julian, Ary Maulana and Ridho Irwansyah, An experimental study on the thermal performance of Nano fluids in microchannel heat exchanger, International Journal of Technology, 2 (2013) 167‐ 177. [8] Mushtaq I. Hasan, Alaa M. Lafta, Ahmed J. Shkarah, Numerical investigation of electroosmotic flow in microchannel, International Review of Mechanical Engineering (I.R.M.E), 11(8) (2017) . [9] Mushtaq I. Hasan, Alaa M. Lafta, Ahmed J. Shkarah, Study the electroosmosis flow in microchannel with contractions of different geometries, proc. Of the 2nd Int. sci. conf., Southern Technical university, (march 2017). [10] Ahmed A. Ali, Mushtaq I. Hassan and Ghassan Adnan, Study surface roughness on the overall performance of parallel flow microchannel heat exchanger, the University of Thi-Qar Journal for Engineering Sciences, 10 (1), 2019. [11] Ahmed A. Ali, Mushtaq I. Hassan, and Ghassan Adnan, Numerical investigation of the slip flow effect in counter-flow microchannel heat exchanger with different channels geometries, First International scientific conference of Al- Ayen University (ISCAU), 2019. [12] Qingkai Zhao, Hang Xu, and Longbin Tao, Flow and heat transfer of nanofluid through a horizontal microchannel with magnetic field and interfacial electrokinetic effects, European Journal of Mechanics/B Fluids, 80 (2020) 72-79. [13] Ahamed Saleel, Asif Afzal, Mostafa Abdel mohimen, Manzoore Elahi, H. Fayaz, Numerical investigation on pressure-driven electroosmatic flow and mixing in a constricted microchannel by triangular obstacle, International Journal of Numerical Methods for Heat & Fluid Flow, (9 June 2020). [14] Mushtaq I. Hasan, Effect of variable fluid properties on the hydrodynamic and thermal characteristics of parallel flow microchannel heat exchanger, Journal of University of Thi-Qar, 10 (4). Dec. (2015). [15] S. Kumar, The EDL Effect in Microchannel Flow: A Critical Review, International Journal of Advanced Computer Research, 3 (4), (2013). [16] Ferrás, A. Afonso, M. Alves, J. Nóbrega and F. Pinho, Electro-osmotic and pressure-driven flow of viscoelastic fluids in microchannels: Analytical and semi- analytical solutions, PHYSICS OF FLUIDS 28 (2016) 093102 [17] Mehdi Sheikhizad and Mohammad Kalteh, Heat transfer investigation of combined electroosmotic/pressure driven nanofluid flow in a microchannel: Effect of heterogeneous surface potential and slip boundary condition, Journal Pre-proof, (2019). [18] Reza Monazami and Mehrdad T. Manzari, Analysis of combined pressure- driven electroosmotic flow through square microchannels, Microfluid Nanofluid, 3 (2007) 123–126.