Open AccessArticle

Effective Boundary Slip Induced by Surface Roughness and Their Coupled Effect on Convective Heat Transfer of Liquid Flow

Yunlu Pan

^1,*

Dalei Jing

^2,*,

He Zhang

¹ and

Xuezeng Zhao

Key laboratory of Micro-Systems and Micro-Structures Manufacturing, Ministry of Education and School of Mechatronics Engineering, Harbin Institute of Technology, Harbin 150001, China

School of Mechanical Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China

Authors to whom correspondence should be addressed.

Entropy 2018, 20(5), 334; https://doi.org/10.3390/e20050334

Submission received: 29 March 2018 / Revised: 22 April 2018 / Accepted: 25 April 2018 / Published: 2 May 2018

(This article belongs to the Special Issue Entropy Generation and Heat Transfer)

Download

Browse Figures

Figure 1
Schematic of the roughness texture in (a) the cone model and (b) the groove model, and the schematic of the flow on the (c) Cone model, (d) Groove-P model and (e) Groove-V model. "> Figure 2
Effective boundary slip (reversed) vs. Ra on the cone model with Rsm = 10 nm (Green); 20 nm (Black); 100 nm (Dark blue); 200 nm (Pink); 1000 nm (Red); 2000 nm (Light blue). The −beff and Ra are in log scale. "> Figure 3
Effective boundary slip (reversed) vs. Ra on three different surface models: Cone model; Groove-V model and Groove-P model. The Ra is in log scale. "> Figure 4
Error of corrected effective boundary slip (obtained experimentally from AFM) on cone model with different Ra and Rsm. "> Figure 5
Increasing rate of Nusselt number for a pressure-driven flow in a micro channel with one-side rough surface in Groove-V model. The Rsm is fixed at 100 nm. The effective boundary slip are also shown. ">

Versions Notes

Abstract

As a significant interfacial property for micro/nano fluidic system, the effective boundary slip can be induced by the surface roughness. However, the effect of surface roughness on the effective slip is still not clear, both increased and decreased effective boundary slip were found with increased roughness. The present work develops a simplified model to study the effect of surface roughness on the effective boundary slip. In the created rough models, the reference position of the rough surfaces to determinate effective boundary slip was set based on ISO/ASME standard and the surface roughness parameters including Ra (arithmetical mean deviation of the assessed profile), Rsm (mean width of the assessed profile elements) and shape of the texture varied to form different surface roughness. Then, the effective boundary slip of fluid flow through the rough surface was analyzed by using COMSOL 5.3. The results show that the effective boundary slip induced by surface roughness of fully wetted rough surface keeps negative and further decreases with increasing Ra or decreasing Rsm. Different shape of roughness texture also results in different effective slip. A simplified corrected method for the measured effective boundary slip was developed and proved to be efficient when the Rsm is no larger than 200 nm. Another important finding in the present work is that the convective heat transfer firstly increases followed by an unobvious change with increasing Ra, while the effective boundary slip keeps decreasing. It is believed that the increasing Ra enlarges the area of solid-liquid interface for convective heat transfer, however, when Ra is large enough, the decreasing roughness-induced effective boundary slip counteracts the enhancement effect of roughness itself on the convective heat transfer.

Keywords:

surface roughness; effective boundary slip; convective heat transfer

1. Introduction

The velocity boundary condition of fluid flow in a micro/nano fluidic system with gradually increasing surface to volume ratio can be a significant interfacial property. For the conventional hydrodynamics, it was believed that the velocity of the fluid close to the solid surface was zero, which means there is no relative motion between the solid and fluid at the boundary [1]. However, in the past decades, lots of studies [2,3,4,5,6,7,8] have found that there is a non-zero velocity condition between the solid surface and the adjacent fluid, presenting the so-called boundary slip condition. Although in macro scale, the exist of the boundary slip can be mostly neglected [9,10], there can be lots of potential applications in micro/nano fluidic systems [11,12,13]. The mechanisms of boundary slip are of interest, but there is still no clear understanding of the slip. Previous studies found that surface wettability is a possible reason to affect slip, however, the effect of wettability on slip can be complex and both hydrophobic and hydrophilic surface can induce either slip or no-slip boundary condition, which is dependent on the strength of solid-liquid interface [14,15,16,17,18]. Additionally, surface roughness including texture spacing and texture height can also significantly affect the degree of slip [18,19,20].

Before analyzing the mechanism of boundary slip, two kinds of slip should be distinguished first. The first one is intrinsic boundary slip which represents the directly relative motion between solid and fluid molecular, and the other one is the effective boundary slip which represents the slippage at a “complex heterogeneous” surface evaluated by averaging of a flow over the length scale of the configuration. The understanding of effective slip can be the “calculated equivalent” slip. For example, the solid surfaces cannot be ideally smooth, then due to the roughness, the reference surface position to determinate the effective slip should be at somewhere between the peak and valley. As a result, the boundary condition will be complicated, then, the effective boundary slip should be calculated to describe the liquid flow on the rough surface. In macro scale, if the roughness is small enough, the effective boundary slip can be neglected, while in micro/nano scale, the effective boundary slip induced by the surface roughness may significantly affect the fluid transportation [21].

There were lots of studies on the effect of surface roughness on the effective boundary slip [22,23,24,25,26,27,28,29,30]. However, some of the results showed a larger slip with a larger roughness while some showed opposite. In most of the measurement techniques of boundary slip, such as the atomic force microscopy (AFM) measurement [31,32,33,34] which is believed to be the most accurate one, the surface roughness can affect both of the effective boundary slip and measurement process [30], thus the effective boundary slip induced by the roughness should be clarified.

Besides effective slip, the surface roughness can also affect the contact area of liquid and solid, that is the heat transfer area for convective heat transfer, thus it will finally affect the convective heat transfer of liquid flow, which is quite important in many applications in micro/nano fluidic system and heat exchangers [35,36,37]. Additionally, the roughness-induced boundary slip can also affect the convective heat transfer by changing the flow behavior. Though the effect of boundary slip on the convective heat transfer has been widely studied [38,39,40,41,42,43], the coupled effect of surface roughness and roughness-induced boundary slip on the convective heat transfer is still lacked. The change of the convective heat transfer with different surface roughness should be considered comprehensively.

To solve this problem, in this study, simulation method based on COMSOL 5.3 is first carried out to obtain the relationship between the roughness parameters and the effective boundary slip. Based on the result of surface roughness-induced effective slip, the convective heat transfer property in a rough microchannel with different surface roughness is studied. The underlying mechanisms are discussed.

2. Simulation Model

Most of the surface roughness is random texture, which is complicated with lots of parameters. These parameters can be mainly divided into three groups: amplitude parameters, spacing parameters and other parameters, among which, the amplitude and spacing parameters are mostly used. In this study, an amplitude parameter Ra, which is the arithmetical mean deviation of the assessed profile of a rough surface, and a spacing parameter Rsm, which is the mean width of the profile elements of a rough surface, are chosen to be the key parameters affecting the effective boundary slip. Thus, the simplified models of the surface roughness on flat surface shown in Figure 1, which only contain simple amplitude and pitch characteristic, are used. Figure 1a gives a cone model with two-dimensional roughness, and Figure 1b gives a groove model with one-dimensional roughness. To study the effect of surface roughness on the boundary slip, values of Ra and Rsm are pre-fixed. Then, the dimensions of roughness textures in the format of cone and groove can be reverse-calculated based on the values of Ra and Rsm to realize the fixed Ra and Rsm. It should be noticed that the reference surface position of the rough surface to define Ra and Rsm is not at the bottom of the profile. Based on ISO/ANSI standard, it should be at some distance above the bottom, the distance m_p of the reference surface position above the bottom can be calculated as [44,45]:

For the cone model:

m_{p} = \frac{1}{S} \int_{S} z d x d y = \frac{π}{12} h

(1)

For the groove model:

m_{p} = \frac{1}{S} \int_{S} z d x d z = \frac{h}{2}

(2)

where S is the projection area of the texture, z is the coordinate along the direction of texture height, and h is the height. Then, the dimensions of cone model and the groove model can be respectively reverse-calculated based Ra and Rsm as follows.

For cone model,

r = \frac{1}{2} R s m, h = \frac{10368}{π {(12 - π)}^{3}} R a \approx 4.75 R a

where r is the is the radius of the bottom circle of the cone.

For groove model,

l = R s m, h = 4 R a

where l is the width of the groove.

Based on these rough surfaces, pressure driven flow through a straight channel confined by a lower side of the created rough surface and a upper side of a fully smooth surface was simulated. Furthermore, it is assumed that there is no boundary slip on both the lower and upper surfaces. The flow on the cone model is shown as Figure 1c, while for the groove model, two situations are set with different direction of flow on it, one’s flow is parallel with the groove which is labeled groove-p model shown as Figure 1d while the other one’s flow is vertical with the groove which is labeled groove-v model shown as Figure 1e.

Models with different Ra and Rsm are built to study the relationship between the roughness and boundary slip and their effect on heat transfer. The distance between the upper wall and the bottom profile of the rough lower wall H* is set as 20,000 nm, and the thickness of the two walls is 2000 nm. The external reference pressure is set as 1 bar (1.01325 × 10⁵ Pa), and room temperature is 293.15 K. Then the density of water ρ = 0.9982071 g/cm³, and dynamic viscosity μ = 1.0050 × 10⁻³ Pa⋅s [46]. Based on the setup, the flow field can be assumed to remains laminar. p_in is set as a constant at the inlet of the channel while p_out is set as 0 Pa at the outlet of the channel. Moreover, the lateral sides of the fluid flow field are defined by symmetric plane. Based on the simulation by COMSOL 5.3, the flow rate Q can be obtained for each model.

The surface roughness will avoidably affect the flow rate of the fluid flow, and this surface roughness-dependent flow rate can be characterized by the effective boundary slip induced by the surface roughness. Assuming the surface roughness-induced effective slip on the lower rough surface is b_eff, then the relationship between effective slip b_eff and flow rate Q can be described as [47]:

b_{e f f} = - \frac{12 μ H Q + (d p / d x) H^{4} W}{4 (d p / d x) H^{3} W + 12 μ Q}

(3)

where dp/dx is the pressure gradient. Therefore, the effective slip b_eff can be obtained with Equation (3). The height H of the channel need to be calculated for each model by H = H^* − m_p considering the surface roughness.

Then, for the cone model:

H = H^{*} - \frac{π}{12} h

(4)

for the groove model:

H = H^{*} - \frac{h}{2}

(5)

The heat transfer model is built by adding the heat transfer properties including convection and conduction to the liquid and solid of the groove-v model, shown as Figure 1e. It is a simulation of a typical water cooling system. The heat source with a constant temperature of 393.15 K is applied on the outside surface of the upper wall of the channel. The flow between the two plates is assumed to be non-isothermal flow, and it will cool the upper plate by convective heat transfer. The other parameter is as same as the previous simulation. The Rsm of the rough surface is fixed at 100 nm while Ra is variable. Based on the simulation by COMSOL 5.3, the bulk convective heat transfer power P can be obtained, which is related to the Nusselt number as shown in Equation (6).

N u = \frac{P H}{k A Δ T}

(6)

where k is the thermal conductivity of the fluid, A is the surface area,

Δ T

is the temperature difference.

3. Results and Discussion

Based on the created models, the effective boundary slip on surfaces with different Ra and Rsm are obtained. Then the convective heat transfer performance of the pressure-drive flow is obtained on groove model with different Ra. To keep the accuracy of the present analysis using COMSOL 5.3, the grid sensitivity is carefully checked. Through the grid sensitivity analysis, the grid in free tetrahedral shape is used and the number of gird is about 2 × 10⁵. The results are shown as following.

3.1. Effective Boundary Slip

The effective boundary slip on the cone model with different Ra and Rsm are shown in Figure 2. Form the results in Figure 2, it can be found that the effective boundary slip induced by the surface roughness is always negative, which means the surface roughness of the microchannel wall in the present work always reduces the flow rate and increases the fluid drag compared to a microchannel confined by two fully smooth walls without slip. Additionally, it is obvious that with the same Rsm, increasing Ra will lead to a decrease of the effective boundary slip and with same Ra, increasing Rsm will lead to an increase of the effective boundary slip.

In the two groove models, the trends of effective boundary slip with Ra and Rsm are similar with that in cone model. However three different models have different slips. The effective boundary slips with the same Rsm of 1000 nm on three different models are shown in Figure 3. It is shown that firstly the slip has the same trend with Ra on three models, and secondly the groove models have the larger slip than cone model.

As a summary, the effective boundary slip induced by the roughness of the rough surface with zero intrinsic boundary slip is always negative, and it will decrease on rougher surface and it also depends on the shape of the surface roughness. This means that the surface roughness in the present work always increases the fluid drag and the fluid drag increase with the increasing roughness. These results are in agreement with the previous studies [28,29,31,48,49,50], and all these theoretical and experimental studies found that friction factor increased with the growing surface relative roughness.

Based on the present work regarding effect of surface roughness on effective boundary slip, the effective boundary slip is related to the surface roughness and the reference surface of the rough surface, thus, most of experimental measurement of slip on rough surface is not the real effective boundary slip and should be corrected. Correction methods were developed by groups [29,30]. Based on definition from ISO/ANSI, the mean plane should be used as the reference position with no doubt. However, in AFM measurement of boundary slip, the reference position was always set at the contact position of the probe and surface. The contact position should be at the top point of the surface morphology at the contact area, since there is a roughness parameter called Rpm which is defined as the distance between the average peak and the mean plane, then the correction can be simply subtract the Rpm from the measured boundary slip.

To check the validity of the correction method, in this study, the reference surface position was set at the top point of surface morphology to simulate the slip obtained experimentally with AFM, labeled as b_AFM. Then the results were corrected by simply subtract the distance between the two different reference surface position which is referred to the R_pm of the surface. The corrected effective boundary slip is marked as b_AFM-C. Then the error of the corrected slip are obtained and shown in Figure 4. In Figure 3, it is shown that with smaller Rsm, the errors are always smaller than ±5% which is acceptable for slip measurement. However, when the Rsm reaches 1000 nm or even 2000 nm, the error can be as large as 15%. Error as 15% for boundary slip measurement seems to be too large, however, it should be noticed that the error of the uncorrected experimentally measured slip is in range of 60–600%, as a result, the correction is still necessary and the error is acceptable, especially when Rsm is smaller than 200 nm.

3.2. Convective Heat Transfer

The bulk convective heat transfer is obtained on the groove-v model. For easier understanding, the bulk convective heat transfer is characterized by the Nusselt number Nu. Nusselt number for the ideally smooth surface is written as Nu₀. The variation of Nusselt number with different Ra is shown in Figure 5. From Figure 5, it can be found that when the Ra increases from 0, the Nusselt number increases faster when Ra is lower than 25 nm, followed by a slow and unobvious increase. The increasing Nusselt number with the increasing Ra when Ra is small is in well agreement with the previous theoretical and experimental studies [48,49,50], this is because that the increase of Ra extends the real contact area of solid and liquid, which can enhance the convective heat transfer. In the meanwhile, the effective boundary slip decreases continually. However, when Ra is larger enough, the decreased effective boundary slip induced by surface roughness will lead to a significant reduction of the liquid flow velocity, which will weaken the convective heat transfer. The enhancement of Ra on the heat transfer by increasing the real contact area of solid and liquid and the reduction of Ra on the heat transfer by decreasing the effective slip counteract each other.

4. Conclusions

In this paper, the effective boundary slip induced by the surface roughness is studied with different Ra and Rsm based on a simplified model by a quantitative simulation. The convective heat transfer power of the pressure driven flow on rough surfaces are studied in the same model as well. The results show that for totally wetting rough surfaces with zero slip, the effective boundary slip induced by the roughness is always negative which means the surface roughness will always increase the drag of liquid flow at the interface, in the meanwhile, the negative effective boundary slip will further decrease with increasing Ra or decreasing Rsm. It is also found that the shape of the roughness texture will affect the effective boundary slip. For the surfaces with smaller Rsm, the simple corrected method for the measured slip by AFM, which subtracts the distance between two reference surface positions from the final results, can effectively decrease the error from up to 600% to 5%. Furthermore, the convective heat transfer can be enhanced with larger Ra of the surface, however, the increase of the convective heat transfer or Nusselt number will be unobvious when the Ra keep increasing due to the decreased effective boundary slip.

Author Contributions

Conceptualization, Y.P. and D.J.; Methodology, Y.P. and D.J.; Software, H.Z.; Validation, Y.P., H.Z. and X.Z.; Formal Analysis, Y.P. and X.Z.; Investigation, X.Z.; Resources, Y.P and D.J.; Data Curation, Y.P., H.Z. and D.J.; Writing-Original Draft Preparation, Y.P. and H.Z.; Writing-Review & Editing, Y.P. and D.J.; Visualization, Y.P. and H.Z.; Supervision, Y.P., J.D. and X.Z.; Project Administration, Y.P.; Funding Acquisition, Y.P.

Funding

The authors gratefully acknowledge the financial support of the National Natural Science Foundation of China (No. 51505108).

Conflicts of Interest

The authors declare no conflict of interest.

References

Stokes, V.K. Couple stresses in fluids. Phys. Fluids 1996, 9, 1709–1715. [Google Scholar] [CrossRef]
Vinogradova, O.I. Drainage of a Thin Liquid-Film Confined between Hydrophobic Surfaces. Langmuir 1995, 11, 2213–2220. [Google Scholar] [CrossRef]
Ou, J.; Perot, B.; Rothstein, J.P. Laminar drag reduction in microchannels using ultrahydrophobic surfaces. Phys. Fluids 2004, 16, 4635–4643. [Google Scholar] [CrossRef]
Choi, C.-H.; Kim, C.-J. Large Slip of Aqueous Liquid Flow over a Nanoengineered Superhydrophobic Surface. Phys. Rev. Lett. 2006, 96, 066001. [Google Scholar] [CrossRef] [PubMed]
Steinberger, A.; Cottin-Bizonne, C.; Kleimann, P.; Charlaix, E. High friction on a bubble mattress. Nat. Mater. 2007, 6, 665–668. [Google Scholar] [CrossRef] [PubMed]
Lee, C.; Kim, C.J. Maximizing the Giant Liquid Slip on Superhydrophobic Microstructures by Nanostructuring Their Sidewalls. Langmuir 2009, 25, 12812–12818. [Google Scholar] [CrossRef] [PubMed]
Wang, Y.; Bhushan, B. Boundary slip and nanobubble study in micro/nanofluidics using atomic force microscopy. Soft Matter 2010, 6, 29–66. [Google Scholar] [CrossRef]
Pan, Y.; Bhushan, B.; Zhao, X. The study of surface wetting, nanobubbles and boundary slip with an applied voltage: A review. Beilstein J. Nanotech. 2014, 5, 1042–1065. [Google Scholar] [CrossRef] [PubMed]
Di Federicao, V.; Longo, S.; King, S.E.; Chiapponi, L.; Petrolo, D.; Ciriello, V. Gravity-driven flow of Herschel–Bulkley fluid in a fracture and in a 2D porous medium. J. Fluid Mech. 2017, 821, 59–84. [Google Scholar] [CrossRef]
Longo, S.; Chiapponi, L.; Di Federicao, V. On the propagation of viscous gravity currents of non-Newtonian fluids in channels with varying cross section and inclination. J. Non-Newton Fluid 2016, 235, 95–108. [Google Scholar] [CrossRef]
Choi, W.; Tuteja, A.; Chhatre, S.; Mabry, J.M.; Cohen, R.E.; McKinley, G.H. Fabrics with tunable oleophobicity. Adv. Mater. 2009, 21, 2190–2195. [Google Scholar] [CrossRef]
Brown, P.S.; Bhushan, B. Bioinspired, roughness-induced, water and oil super-philic and super-phobic coatings prepared by adaptable layer-by-layer technique. Sci. Rep. 2015, 5, 14030. [Google Scholar] [CrossRef] [PubMed]
Li, F.; Wang, Z.; Huang, S.; Pan, Y.; Zhao, X. Flexible, Durable, and Unconditioned Superoleophobic/Superhydrophilic Surfaces for Controllable Transport and Oil–Water Separation. Adv. Func. Mater. 2018, 28. [Google Scholar] [CrossRef]
Bakli, C.; Chakraborty, S. Slippery to Sticky Transition of Hydrophobic Nanochannels. Nano Lett. 2015, 15, 7497–7502. [Google Scholar] [CrossRef] [PubMed]
Vo, T.Q.; Park, B.; Park, C.; Kim, B. Nano-scale liquid film sheared between strong wetting surfaces: Effects of interface region on the flow. J. Mech. Sci. Technol. 2015, 29, 1681–1688. [Google Scholar] [CrossRef]
Bao, L.; Priezjev, N.V.; Hu, H.; Luo, K. Effects of viscous heating and wall-fluid interaction energy on rate-dependent slip behavior of simple fluids. Phys. Rev. E 2017, 96. [Google Scholar] [CrossRef] [PubMed]
Ghorbanian, J.; Beskok, A. Scale effects in Nano-channel liquid flows. Microfluid Nanofluid 2016, 20, 121. [Google Scholar] [CrossRef]
Jeong, M.; Kim, Y.; Zhou, W.; Tao, W.; Ha, M.Y. Effects of surface wettability, roughness and moving wall velocity on the Couette flow in nano-channel using multi-scale hybrid method. Comput. Fluid. 2017, 147, 1–11. [Google Scholar] [CrossRef]
Noorian, H.; Toghraie, D.; Azimian, A.R. Molecular dynamics simulation of Poiseuille flow in a rough nano channel with checker surface roughnesses geometry. Heat Mass Transf. 2014, 50, 105–113. [Google Scholar] [CrossRef]
Cao, B.; Chen, M.; Guo, Z. Effect of surface roughness on gas flow in microchannels by molecular dynamics simulation. Int. J. Eng. Sci. 2006, 44, 927–937. [Google Scholar] [CrossRef]
Nosonovsky, M.; Bhushan, B. Multiscale Dissipative Mechanisms and Hierarchical Surfaces: Friction, Superhydrophobicity, and Biomimetics; Springer: Heidelberg, Germany, 2008. [Google Scholar]
Baldoni, F. On slippage induced by surface diffusion. J. Eng. Math. 1996, 30, 647–659. [Google Scholar] [CrossRef]
Jabbarzadeh, A.; Atkinson, J.D.; Tanner, R.I. Effect of the wall roughness on slip and rheological properties of hexadecane in molecular dynamics simulation of Couette shear flow between two sinusoidal walls. Phys. Rev. E 2000, 61, 690–699. [Google Scholar] [CrossRef]
Pit, R.; Hervet, H.; Leger, L. Direct experimental evidence of slip in hexadecane: Solid interfaces. Phys. Rev. Lett. 2000, 85, 980–983. [Google Scholar] [CrossRef] [PubMed]
Ponomarev, I.V.; Meyerovich, A.E. Surface roughness and effective stick-slip motion. Phys. Rev. E 2003, 67, 026302. [Google Scholar] [CrossRef] [PubMed]
Bonaccurso, E.; Butt, H.-J.; Craig, V.S.J. Surface Roughness and Hydrodynamic Boundary Slip of a Newtonian Fluid in a Completely Wetting System. Phys. Rev. E 2003, 90, 144501. [Google Scholar] [CrossRef] [PubMed]
Zhu, Y.; Granick, S. Limits of the Hydrodynamic No-Slip Boundary Condition. Phys. Rev. Lett. 2002, 88, 106102. [Google Scholar] [CrossRef] [PubMed]
Schmatko, T.; Hervet, H.; Léger, L. Effect of nanometric-scale roughness on slip at the wall of simple fluids. Langmuir 2006, 22, 6843–6850. [Google Scholar] [CrossRef] [PubMed]
Guriyanova, S.; Semin, B.; Rodrigues, T.S.; Butt, H.-J.; Bonaccurso, E. Hydrodynamic drainage force in a highly confined geometry: Role of surface roughness on different length scales. Microfluid Nanofluid 2010, 8, 653–663. [Google Scholar] [CrossRef]
Pan, Y.; Jing, D.; Zhao, X. Effect of Surface Roughness on the Measurement of Boundary Slip Based on Atomic Force Microscope. Sci. Adv. Mater. 2017, 9, 122–127. [Google Scholar] [CrossRef]
Vinogradova, O.I.; Yakubov, G.E. Surface roughness and hydrodynamic boundary conditions. Phys. Rev. E 2006, 73, 045302(R). [Google Scholar] [CrossRef] [PubMed]
Cottin-Bizonne, C.; Cross, B.; Steinberger, A.; Charlaix, E. Boundary slip on smooth hydrophobic surfaces: Intrinsic effects and possible artifacts. Phys. Rev. Lett. 2005, 94, 056102. [Google Scholar] [CrossRef] [PubMed]
Bhushan, B.; Wang, Y.; Maali, A. Boundary Slip Study on Hydrophilic, Hydrophobic, and Superhydrophobic Surfaces with Dynamic Atomic Force Microscopy. Langmuir 2009, 25, 8117–8121. [Google Scholar] [CrossRef] [PubMed]
Jing, D.; Bhushan, B. The coupling of surface charge and boundary slip at the solid–liquid interface and their combined effect on fluid drag: A review. J. Colloid Interface Sci. 2015, 454, 152–179. [Google Scholar] [CrossRef] [PubMed]
Kandlikar, S.G.; Garimella, S.; Li, D.Q.; Colin, S.; King, M.R. Heat Transfer and Fluid Flow in Minichannels and Microchannels; Elsevier: Oxford, UK, 2006. [Google Scholar]
Li, D. Encyclopedia of Microfluidics and Nanofluidics; Springer Science & Business Media: Berlin, Germany, 2008. [Google Scholar]
Lin, B. Microfluidics: Technologies and Applications; Springer: Berlin, Germany, 2011. [Google Scholar]
Ngoma, G.D.; Erchiqui, F. Heat flux and slip effects on liquid flow in a microchannel. Int. J. Therm. Sci. 2007, 46, 1076–1083. [Google Scholar] [CrossRef]
Yavari, H.; Sadeghi, A.; Saidi, M.H.; Chakraborty, S. Combined influences of viscous dissipation, non-uniform Joule heating and variable thermophysical properties on convective heat transfer in microtubes. Int. J. Heat Mass Transf. 2012, 55, 762–772. [Google Scholar] [CrossRef]
Tan, D.K.; Liu, Y. Combined effects of streaming potential and wall slip on flow and heat transfer in microchannels. Int. Commun. Heat Mass 2014, 53, 39–42. [Google Scholar] [CrossRef]
Keramati, H.; Sadeghi, A.; Saidi, M.H.; Chakraborty, S. Analytical solutions for thermo-fluidic transport in electroosmotic flow through rough microtubes. Int. J. Heat Mass Transf. 2016, 92, 244–251. [Google Scholar] [CrossRef]
Jing, D.; Pan, Y. Electroviscous effect and convective heat transfer of pressure-driven flow through microtubes with surface charge-dependent slip. Int. J. Heat Mass Transf. 2016, 101, 648–655. [Google Scholar] [CrossRef]
Jing, D.; Pan, Y.; Wang, X. Joule heating, viscous dissipation and convective heat transfer of pressure-driven flow in a microchannel with surface charge-dependent slip. Int. J. Heat Mass Transf. 2017, 108, 1305–1313. [Google Scholar] [CrossRef]
Anonymous. Geometrical Product Specifications—Surface Texture: Profile Method-Nominal Characteristics of Contact (Stylus) Instruments. ISO3274. International Standardization Organization: Geneva, Switzerland, 1 December 1996. [Google Scholar]
Anonymous. Surface Texture (Surface Roughness, Waviness, and Lay). ANSI/ASME B46.1-2009. ASME: New York, NY, USA, 20 August 2010. [Google Scholar]
Haynes, W.M. Handbook of Chemistry and Physics, 96th ed.; Academic Press: New York, NY, USA, 2015. [Google Scholar]
Huang, Y.; Zhao, X.; Pan, Y.; Ahmad, K. Simulation of effective slip and drag in pressure-driven flow on superhydrophobic surfaces. J. Nanomater. 2016, 2016. [Google Scholar] [CrossRef]
Li, Z.; He, Y.L.; Tang, G.H.; Tao, W.Q. Experimental and numerical studies of liquid flow and heat transfer in microtubes. Int. J. Heat. Mass Transf. 2007, 50, 3447–3460. [Google Scholar] [CrossRef]
Wu, H.Y.; Cheng, P. An experimental study of convective heat transfer in silicon microchannels with different surface conditions. Int. J. Heat. Mass Transf. 2003, 46, 2547–2556. [Google Scholar] [CrossRef]
Guo, L.; Xu, H.; Gong, L. Influence of wall roughness models on fluid flow and heat transfer in microchannels. Appl. Therm. Eng. 2015, 84, 399–408. [Google Scholar] [CrossRef]

Figure 1. Schematic of the roughness texture in (a) the cone model and (b) the groove model, and the schematic of the flow on the (c) Cone model, (d) Groove-P model and (e) Groove-V model.

Figure 2. Effective boundary slip (reversed) vs. Ra on the cone model with Rsm = 10 nm (Green); 20 nm (Black); 100 nm (Dark blue); 200 nm (Pink); 1000 nm (Red); 2000 nm (Light blue). The −b_eff and Ra are in log scale.

Figure 3. Effective boundary slip (reversed) vs. Ra on three different surface models: Cone model; Groove-V model and Groove-P model. The Ra is in log scale.

Figure 4. Error of corrected effective boundary slip (obtained experimentally from AFM) on cone model with different R_a and R_sm.

Figure 5. Increasing rate of Nusselt number for a pressure-driven flow in a micro channel with one-side rough surface in Groove-V model. The Rsm is fixed at 100 nm. The effective boundary slip are also shown.

© 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Pan, Y.; Jing, D.; Zhang, H.; Zhao, X. Effective Boundary Slip Induced by Surface Roughness and Their Coupled Effect on Convective Heat Transfer of Liquid Flow. Entropy 2018, 20, 334. https://doi.org/10.3390/e20050334

AMA Style

Pan Y, Jing D, Zhang H, Zhao X. Effective Boundary Slip Induced by Surface Roughness and Their Coupled Effect on Convective Heat Transfer of Liquid Flow. Entropy. 2018; 20(5):334. https://doi.org/10.3390/e20050334

Chicago/Turabian Style

Pan, Yunlu, Dalei Jing, He Zhang, and Xuezeng Zhao. 2018. "Effective Boundary Slip Induced by Surface Roughness and Their Coupled Effect on Convective Heat Transfer of Liquid Flow" Entropy 20, no. 5: 334. https://doi.org/10.3390/e20050334

APA Style

Pan, Y., Jing, D., Zhang, H., & Zhao, X. (2018). Effective Boundary Slip Induced by Surface Roughness and Their Coupled Effect on Convective Heat Transfer of Liquid Flow. Entropy, 20(5), 334. https://doi.org/10.3390/e20050334

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu