Research Papers

DSMC Simulation of Low Knudsen Micro/Nanoflows Using Small Number of Particles per Cells

[+] Author and Article Information
Ehsan Roohi

e-mail: e.roohi@ferdowsi.um.ac.ir

Hamid Niazmand

Department of Mechanical Engineering,
Faculty of Engineering,
Ferdowsi University of Mashhad,
P.O. Box 91775-1111,
Mashhad, Iran

Stefan Stefanov

Institute of Mechanics,
Bulgarian Academy of Science,
Acad. G. Bontchev Street,
Sofia 1113, Bulgaria

1Corresponding author.

Contributed by the Heat Transfer Division of ASME for publication in the Journal of Heat Transfer. Manuscript received July 13, 2012; final manuscript received February 1, 2013; published online September 11, 2013. Assoc. Editor: Sushanta K. Mitra.

J. Heat Transfer 135(10), 101008 (Sep 11, 2013) (8 pages) Paper No: HT-12-1367; doi: 10.1115/1.4024505 History: Received July 13, 2012; Revised February 01, 2013

Direct simulation Monte Carlo (DSMC) method in low Knudsen rarefied flows at micro/nanoscales remains a big challenge for researchers due to large computational requirements. In this article, the application of the simplified Bernoulli-trials (SBT)/dual grid collision scheme is extended for solving low Knudsen/low speed and low Knudsen/high gradient rarefied micro/nanoflows. The main advantage of the SBT algorithm is to provide accurate calculations using much smaller number of particles per cell, i.e., 〈N〉 ≈ 2, which is quite beneficial for near continuum DSMC simulations where the requirement of fine meshes faces the simulation with serious memory restrictions. Comparing the results of the SBT/dual grid scheme with the no time counter (NTC) scheme and majorant frequency scheme (MFS), it is shown that the SBT/dual grid scheme could successfully predict the thermal pattern and hydrodynamics field as well as surface parameters such as velocity slip, temperature jump and wall heat fluxes. Therefore, we present SBT/dual grid algorithm as a suitable alternative of the standard collision schemes in the DSMC method for typical micro/nanoflows solution. Nonlinear flux-corrected transport (FCT) algorithm is also employed as a filter to extract the smooth solution from the noisy DSMC calculation for low speed/low Knudsen number DSMC calculations.

Copyright © 2013 by ASME
Your Session has timed out. Please sign back in to continue.


Ho, C., and Tai, Y., 1998, “Micro-Electro-Mechanical Systems (MEMS) and Fluid Flows,” Annu. Rev. Fluid Mech., 30, pp. 579–612. [CrossRef]
Beskok, A., and Karniadakis, G. E., 1999, “A Model for Flows in Channels, Pipes and Ducts at Micro and Nano Scale,” J. Microscale Thermophys. Eng., 3, pp. 43–77. [CrossRef]
Schaaf, S. A., and Chambre, P. L., 1961, Flow of Rarefied Gases, Princeton University Press, Princeton, NJ.
Karniadakis, G., Beskok, A., and Aluru, N., 2005, Microflows and Nanoflows: Fundamentals and Simulation, Springer-Verlag, New York.
Barber, R. W., and Emerson, D. R., 2006, “Challenges in Modeling Gas-Phase Flow in Microchannels: From Slip to Transition,” Heat Transfer Eng., 27 (4), pp. 3–12. [CrossRef]
Bird, G. A., 1976, Molecular Gas Dynamics, Clarendon Press, Oxford, UK.
Wagner, W., 1992, “A Convergence Proof for Bird's Direct Simulation Monte Carlo Method for the Boltzmann Equation,” J. Stat. Phys., 66, pp. 1011–1044. [CrossRef]
Fan, J., and Shen, C., 2002, “Micro-Scale Gas Flows,” Adv. Mech., 32, pp. 321–336.
Kaplan, C. R., and Oran, E. S., 2002, “Nonlinear Filtering for Low-Velocity Gaseous Microflows,” AIAA J., 40, pp. 82–90. [CrossRef]
Alexander, F. J., Garcia, A. L., and Alder, B. J., 1998, “Cell Size Dependence of Transport Coefficients in Stochastic Particle Algorithms,” Phys. Fluids, 10(6), pp. 1540–1542. [CrossRef]
Hadjiconstantinou, N. G., 2000, “Analysis of Discretization in the Direct Simulation Monte Carlo,” Phys. Fluids, 12(10), pp. 2634–2638. [CrossRef]
Bird, G. A., 1994, Molecular Gas Dynamics and the Direct Simulation of Gas Flows, Oxford University Press, Oxford, UK.
Ivanov, M. S., and Rogasinsky, S. V., 1988, “Analysis of the Numerical Techniques of the Direct Simulation Monte Carlo Method in the Rarefied Gas Dynamics,” Sov. J. Numer. Anal. Math. Model., 3(6), pp. 453–465.
Ivanov, M. S., Markelov, G. N., and Gimelshein, S. F., 1998, AIAA Paper No. 98-2669.
Stefanov, S. K., 2011, “On DSMC Calculation of Rarefied Gas Flows With Small Number of Particles in Cells,” SIAM J. Sci. Comput., 33(2), pp. 677–702. [CrossRef]
Amiri-Jaghargh, A., Roohi, E., Niazmand, H. and Stefanov, S., 2012, “Low Speed/Low Rarefaction Flow Simulation in Micro/Nano Cavity Using DSMC Method With Small Number of Particles per Cell,” J. Phys.: Conf. Ser., 362, p. 012007.
Roohi, E., Darbandi, M., and Mirjalili, V., 2009, “DSMC Solution of Subsonic Flow Through Micro-Nano Scale Channels,” ASME J. Heat Transfer, 131(9), p. 092402. [CrossRef]
Darbandi, M., and Roohi, E., 2011, “DSMC Simulation of Subsonic Flow Through Nanochannels and Micro/Nano Steps,” Int. Commun. Heat Mass Transfer, 38(10), pp. 1444–1449. [CrossRef]
Roohi, E., and Darbandi, M., 2012, “Recommendations on Performance of Parallel DSMC Algorithm in Solving Subsonic Nanoflows,” Appl. Math. Modell., 36(5), pp. 2314–2321. [CrossRef]
Ejtehadi, O., Roohi, E., and Abolfazli, J., 2012, “Investigation of Basic Molecular Gas Structure Effects on Hydrodynamics and Thermal Behaviors of Rarefied Shear Driven Flow Using DSMC,” Int. Commun. Heat Mass Transfer, 39, pp. 439–448. [CrossRef]
Mohammadzadeh, A., Roohi, E., Niazmand, H., Stefanov, S., and Myong, R. S., 2012, “Detailed Investigation of Thermal and Hydrodynamic Behavior in Micro/Nano Cavity Using DSMC,” Phys. Rev. E, 85, p. 056305. [CrossRef]
Mohammadzadeh, A., Roohi, E., and Niazmand, H., 2013, “A Parallel DSMC Investigation of Monoatomic/Diatomic Gas Flows in a Micro/Nano Cavity,” Numer. Heat Transfer, Part A, 63, pp. 305–325. [CrossRef]
Mizzi, S., Emerson, D. R., Stefanov, S. K., Barber, R. W., and Reese, J. M., 2007, “Effects of Rarefaction on Cavity Flow in the Slip Regime,” J. Comput. Theor. Nanosci., 4(4), pp. 817–822.
John, B., Gu, X. J., and Emerson, D. R., 2011, “Effects of Incomplete Surface Accommodation on Non-Equilibrium Heat Transfer in Cavity Flow: A Parallel DSMC Study,” Comput. Fluids, 45, pp. 197–201. [CrossRef]
John, B., Gu, X. J., and Emerson, D. R., 2010, “Investigation of Heat and Mass Transfer in a Lid-Driven Cavity Under Nonequilibrium Flow Conditions,” Numer. Heat Transfer, Part B, 58, pp. 287–303. [CrossRef]
Gu, X. J., and Emerson, D. R., 2011, “Application of the Moment Method in the Slip and Transition Regime for Microfluidic Flows,” RTO-EN-AVT-194 RTO AVT/VKI Lecture Series.
Jiang, J. Z., Fan, J., and Shen, C., 2003, “Statistical Simulation of Micro-Cavity Flows,” AIP Conference Proceedings Rarefied Gas Dynamics: 23rd International Symposium, 663, p. 784.
Lofthouse, A. J., 2008, “Nonequilibrium Hypersonic Aerothermodynamics Using the Direct Simulation Monte Carlo and Navier–StokesModels,” Ph.D. thesis, University of Michigan, Ann Arbor, MI.
Lengrand, J. C., Allegre, J., Chpoun, A., and Raffin, M., 1994, “Rarefied Hypersonic Flow Over a Sharp Flat Plate: Numerical and Experimental Results,” Rarefied Gas Dynamics: Space Science and Engineering, B. D.Shizgal and D. P.Weaver, eds., AIAA, Washington DC, p. 276.
Woronowicz, M. S., and Baganoff, D., 1994, “ Application of a Vectorized Particle Simulation to the Study of Flat Plates in Rarefied Flow,” Rarefied Gas Dynamics: Theory and Simulations, B. D.Shizgal and D. P.Weaver, eds., AIAA, Washington DC, p. 419.
Eggers, J., and Beylich, A. E., 1995, “Development of a Hybrid Scheme and Its Application to a Flat Plate Flow,” Rarefied Gas Dynamics, J.Harvey and G.Lord, eds., Oxford University Press, Oxford, UK, p. 1216.


Grahic Jump Location
Fig. 1

The simplified Bernoulli-trials collision procedure

Grahic Jump Location
Fig. 2

DSMC flowchart for staggered grid

Grahic Jump Location
Fig. 3

The majorant frequency scheme collision procedure

Grahic Jump Location
Fig. 4

Geometrical configuration of microcavity

Grahic Jump Location
Fig. 5

Grid independency test

Grahic Jump Location
Fig. 6

Effect of number of particle per cell in accuracy of the NTC collision scheme

Grahic Jump Location
Fig. 7

Effect of applying two sampling in each iteration

Grahic Jump Location
Fig. 8

Comparison of the SBT and NTC schemes in prediction of flow field

Grahic Jump Location
Fig. 9

Comparison of the SBT and NTC schemes in prediction of thermal pattern

Grahic Jump Location
Fig. 10

Comparison of SBT and NTC schemes in prediction of centerline velocity in low speed cavity flows (case 2)

Grahic Jump Location
Fig. 11

Velocity slip and Temperature jump along the driven lid computed by SBT and NTC schemes: (a) Case 1 and (b) case 2

Grahic Jump Location
Fig. 12

Accuracy of SBT scheme in prediction of wall heat flux (case 1)

Grahic Jump Location
Fig. 13

Accuracy of SBT scheme in prediction of wall shear stress (case 1)

Grahic Jump Location
Fig. 14

The geometry and flow configuration of nanoflat plate

Grahic Jump Location
Fig. 15

Effect of number of simulators on the accuracy of the results. The profiles are shifted to right for MFS and SBT schemes for more clarity.

Grahic Jump Location
Fig. 16

Accuracy of three schemes in evaluation of temperature profile with small number of particle per cell

Grahic Jump Location
Fig. 17

Comparison of SBT and NTC scheme in evaluation of thermal pattern




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