0
Research Papers

Multiscale Study of Gas Slip Flows in Nanochannels

[+] Author and Article Information
Quy Dong To

Laboratoire Modélisation et
Simulation Multi Échelle,
UMR-CNRS 8208,
Université Paris-Est,
5 Boulevard Descartes,
Marne-la-Vallée Cedex 2 77454, France
e-mail: quy-dong.to@univ-paris-est.fr

Thanh Tung Pham, Vincent Brites, Céline Léonard, Guy Lauriat

Laboratoire Modélisation et
Simulation Multi Échelle,
UMR-CNRS 8208,
Université Paris-Est,
5 Boulevard Descartes,
Marne-la-Vallée Cedex 2 77454, France

1Present address: Université d'Evry Val d'Essonne, Laboratoire Analyse et Modélisation pour la Biologie et l'Environnement, LAMBE CNRS UMR 8587, Boulevard F. Mitterrand, Evry Cedex 91025, France.

Manuscript received January 6, 2014; final manuscript received March 17, 2014; published online May 14, 2015. Assoc. Editor: Yogesh Jaluria.

J. Heat Transfer 137(9), 091002 (Sep 01, 2015) (8 pages) Paper No: HT-14-1009; doi: 10.1115/1.4030205 History: Received January 06, 2014; Revised March 17, 2014; Online May 14, 2015

A multiscale modeling of the anisotropic slip phenomenon for gas flows is presented in a tree-step approach: determination of the gas–wall potential, simulation and modeling of the gas–wall collisions, simulation and modeling of the anisotropic slip effects. The density functional theory (DFT) is used to examine the interaction between the Pt–Ar gas–wall couple. This potential is then passed into molecular dynamics (MD) simulations of beam scattering experiments in order to calculate accommodation coefficients. These coefficients enter in an effective gas–wall interaction model, which is the base of efficient MD simulations of gas flows between anisotropic surfaces. The slip effects are quantified numerically and compared with simplified theoretical models derived in this paper. The paper demonstrates that the DFT potential is in good agreement with empirical potentials and that an extension of the Maxwell model can describe anisotropic slip effects due to surface roughness, provided that two tangential accommodation parameters are introduced. MD data show excellent agreement with the tensorial slip theory, except at large Kundsen numbers (for example, Kn 0.2) and with an analytical expression which predicts the ratio between transverse and longitudinal slip velocity components.

FIGURES IN THIS ARTICLE
<>
Copyright © 2015 by ASME
Your Session has timed out. Please sign back in to continue.

References

Dovesi, R., Orlando, R., Civalleri, B., Roetti, C., Saunders, V. R., and Zicovich-Wilson, C. M., 2005, “Crystal: A Computational Tool for the ab Initio Study of the Electronic Properties of Crystals,” Z. Kristallogr., 220, pp. 571–573. [CrossRef]
Léonard, C., Brites, V., Pham, T. T., To, Q. D., and Lauriat, G., 2013, “Influence of the Pairwise Potential on the Tangential Momentum Accommodation Coefficient: A Multi-Scale Study Applied to the Argon on Pt(111) System,” Eur. Phys. J. B, 86. [CrossRef]
Doll, K., 2004, “CO Adsorption on the Pt (111) Surface: A Comparison of a Gradient Corrected Functional and a Hybrid Functional,” Surf. Sci., 573(3), pp. 464–473. [CrossRef]
Nada, R., Nicholas, J. B., McCarthy, M. I., and Hess, A. C., 1996, “Basis Sets for ab Initio Periodic Hartree-Fock Studies of Zeolite/Adsorbate Interactions: He, Ne, and Ar in Silica Sodalite,” Int. J. Quantum Chem., 60(4), pp. 809–820. [CrossRef]
Kulginov, D., Persson, M., Rettner, C. T., and Bethune, D. S., 1996, “An Empirical Interaction Potential for the Ar/Pt (111) System,” J. Chem. Phys., 100(19), pp. 7919–7927. [CrossRef]
Ramseyer, C., Hoang, P. N. M., and Girardet, C., 1994, “Interpretation of High-Order Commensurate Phases for an Argon Monolayer Adsorbed on Pt(111),” Phys. Rev. B, 49, pp. 2861–2868. [CrossRef]
Maruyama, S., and Kimura, T., 1999, “A Study on Thermal Resistance Over a Solid–Liquid Interface by the Molecular Dynamics Method,” Therm. Sci. Eng., 7, pp. 63–68.
Spijker, P., Markvoort, A. J., Nedea, S. V., and Hilbers, P. A. J., 2010, “Computation of Accommodation Coefficients and the Use of Velocity Correlation Profiles in Molecular Dynamics Simulations,” Phys. Rev. E, 81(1), p. 011203. [CrossRef]
Head-Gordon, M., Tully, J. C., Rettner, C. T., Mullins, C. B., and Auerbach, D. J., 1991, “On the Nature of Trapping and Desorption at High Surface Temperatures. Theory and Experiments for the Ar–Pt (111) System,” J. Chem. Phys., 94(2), pp. 1516–1527. [CrossRef]
Svanberg, M., and Pettersson, J. B. C., 1996, “Survival of Noble Gas Clusters Scattering From Hot Metal Surfaces,” Chem. Phys. Lett., 263(5), pp. 661–666. [CrossRef]
Yamamoto, K., 2002, “Slip Flow Over a Smooth Platinum Surface,” JSME Int. J. Ser. B, 45(4), pp. 788–795. [CrossRef]
Smith, R. J., Kara, A., and Holloway, S., 1993, “A Molecular Dynamics Study for the Trapping and Scattering of Ar/Pt(111),” Surf. Sci., 281(3), pp. 296–308. [CrossRef]
Lahaye, R. J. W. E., Stolte, S., Kleyn, A. W., Smith, R. J., and Holloway, S., 1994, “Site Dependent Energy Loss in Ar Scattering From Pt(111),” Surf. Sci., 307–309, pp. 187–192. [CrossRef]
Zeppenfeld, P., Becher, U., Kern, K., David, R., and Comsa, G., 1990, “Van Hove Anomaly in the Phonon Dispersion of Monolayer Ar/Pt (111),” Phys. Rev. B, 41(12), pp. 8549–8552. [CrossRef]
Sutton, A. P., and Chen, J., 1990, “Long-Range Finnis–Sinclair Potentials,” Philos. Mag. Lett., 61(3), pp. 139–146. [CrossRef]
Pham, T. T., To, Q. D., Lauriat, G., Léonard, C., and Hoang, V. V., 2012, “Effects of Surface Morphology and Anisotropy on the Tangential-Momentum Accommodation Coefficient Between Pt(100) and Ar,” Phys. Rev. E, 86, p. 051201. [CrossRef]
Finger, G. W., Kapat, J. S., and Bhattacharya, A., 2007, “Molecular Dynamics Simulation of Adsorbent Layer Effect on Tangential Momentum Accommodation Coefficient,” ASME J. Fluid Eng., 129(1), pp. 31–39. [CrossRef]
Rettner, C. T., 1998, “Thermal and Tangential-Momentum Accommodation Coefficients for N2 Colliding With Surfaces of Relevance to Disk-Drive Air Bearings Derived From Molecular Beam Scattering,” IEEE Trans. Magn., 34(4), pp. 2387–2395. [CrossRef]
Schlick, T., 2010, Molecular Modeling and Simulation: An Interdisciplinary Guide, Springer, New York.
Dadzie, S. K., and Méolans, J. G., 2004, “Anisotropic Scattering Kernel: Generalized and Modified Maxwell Boundary Conditions,” J. Math. Phys., 45(5), p. 1804. [CrossRef]
Pham, T. T., To, Q. D., Lauriat, G., and Léonard, C., 2013, “Tensorial Slip Theory for Gas Flows and Comparison With Molecular Dynamics Simulations Using an Anisotropic Gas–Wall Collision Mechanism,” Phys. Rev. E, 87, p. 053012. [CrossRef]
To, Q. D., Bercegeay, C., Lauriat, G., Léonard, C., and Bonnet, G., 2010, “A Slip Model for Micro/Nano Gas Flows Induced by Body Forces,” Microfluid. Nanofluid., 8(3), pp. 417–422. [CrossRef]
Karniadakis, G., Beskok, A., and Aluru, N., 2005, Microflows and Nanoflows: Fundamentals and Simulation, Springer, New York.
Struchtrup, H., 2005, Macroscopic Transport Equations for Rarefied Gas Flows: Approximation Methods in Kinetic Theory, Springer, Berlin.
Liou, W. W., and Fang, Y., 2005, Microfluid Mechanics: Principles and Modeling, McGraw-Hill, New York.
Bhattacharya, D. K., and Lie, G. C., 1989, “Molecular Dynamics Simulations of Nonequilibrium Heat and Momentum Transport in Very Dilute Gases,” Phys. Rev. Lett., 62(8), pp. 897–900. [CrossRef] [PubMed]
Morris, D. L., Hannon, L., and Garcia, A. L., 1992, “Slip Length in a Dilute Gas,” Phys. Rev. A, 46(8), pp. 5279–5281. [CrossRef] [PubMed]
To, Q. D., Pham, T. T., Lauriat, G., and Léonard, C., 2012, “Molecular Dynamics Simulations of Pressure-Driven Flows and Comparison With Acceleration-Driven Flows,” Adv. Mech. Eng., 4, p. 580763. [CrossRef]
Bazant, M. Z., and Vinogradova, O. I., 2008, “Tensorial Hydrodynamic Slip,” J. Fluid Mech., 613, pp. 125–134. [CrossRef]
Tehver, R., Toigo, F., Koplik, J., and Banavar, J., 1998, “Thermal Walls in Computer Simulations,” Phys. Rev. E, 57, pp. 17–20. [CrossRef]
Bhattacharya, D. K., and Lie, G. C., 1991, “Nonequilibrium Gas Flow in the Transition Regime: A Molecular-Dynamics Study,” Phys. Rev. A, 43(2), pp. 761–767. [CrossRef] [PubMed]
Bird, G. A., 1994, Molecular Gas Dynamics and the Direct Simulation of Gas Flows, Clarendon, Oxford, UK.
Belyaev, A. V., and Vinogradova, O. I., 2010, “Effective Slip in Pressure-Driven Flow Past Super-Hydrophobic Stripes,” J. Fluid Mech., 652, pp. 489–499. [CrossRef]
Priezjev, N. V., 2011, “Molecular Diffusion and Slip Boundary Conditions at Smooth Surfaces With Periodic and Random Nanoscale Textures,” J. Chem. Phys., 135(20), p. 204704. [CrossRef] [PubMed]
Albertoni, S., Cercignani, C., and Gotusso, L., 1963, “Numerical Evaluation of the Slip Coefficient,” Phys. Fluids, 6(7), pp. 993–996. [CrossRef]

Figures

Grahic Jump Location
Fig. 2

Nanotextured surface with strips and gas beaming direction θ, ϕ in cartesian coordinate system. The geometric parameters are l1 = 11.76 Å, l2 = 19.6 Å, l3 = 7.84 Å, h = 0-5.88 Å.

Grahic Jump Location
Fig. 1

Pt(111)–Ar interaction potentials

Grahic Jump Location
Fig. 3

σtdir computed for striped walls versus azimuth angle ϕ for different roughnesses (Tw = 300 K, θ=45 deg). The solid, dashed, dashed-dotted lines are the analytical expressions (7) used to fit the present numerical results.

Grahic Jump Location
Fig. 4

Collision between gas atoms and a solid wall. The number of gas atoms going downward (velocity v') and upward (velocity v) within one time unit are denoted N and N+, respectively. If there is no gas accumulation at the wall, N+ = N-.

Grahic Jump Location
Fig. 6

Dimensionless slip length Ls/H as a function of ϕ at different Kn. Points are MD data which are fitted with solid lines corresponding to the analytical expression (30). The dashed lines represent quantitative estimation (31) based on αx,αy.

Grahic Jump Location
Fig. 7

Ratio of transverse and longitudinal components of fitted slip velocity as a function of ϕ for various Kn. Points are MD data, the solid and dashed lines are analytical expressions (32) and (33).

Grahic Jump Location
Fig. 8

Ratio of transverse and longitudinal components of real slip velocity as a function of ϕ for different Kn. Points are MD data and the dashed line is for the analytical expression (33).

Grahic Jump Location
Fig. 5

Longitudinal and transverse velocity profiles un, um for different values of ϕ and Kn = 0.104. The velocities are normalized with umax—the velocity at z = 0 for case ϕ = 90 deg.

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In