2-D DSMC algorithm based on Delaunay triangles

WANG Chuncai; CHENG Jia; JI Linhong; LU Yijia; SUN Yuchun; LIN Jia

doi:10.16511/j.cnki.qhdxxb.2015.22.010

PDF(1887 KB)

Journal of Tsinghua University(Science and Technology) ›› 2015, Vol. 55 ›› Issue (10) : 1079-1086,1097. DOI: 10.16511/j.cnki.qhdxxb.2015.22.010

MECHANICAL ENGINEERING

2-D DSMC algorithm based on Delaunay triangles

{{article.zuoZhe_EN}}

Author information +

History +

Abstract

Traditional direct simulation Monte Carlo (DSMC) models using rectangular elements inevitably have errors when dealing with complex computational domain boundaries. Unstructured grids are able to adapt to any complex geometric shape but are rarely applied because of the highly disordered topological structure, complex calculations and low execution efficiency. This paper describes a particle trajectory tracking algorithm based on an unstructured Delaunay triangle mesh. The computational domain is divided into various rectangular regions with only a small number of triangular elements with the particle trajectory tracking achieved by first searching the rectangular regions and then the triangular elements. The algorithm is used in a two-dimensional DSMC program with the results comparing well with results in the literature. The DSMC program is then used to predict the flow field in a showerhead aperture with argon as the neutral reference gas with an inlet pressure of 200 Pa and an inlet temperature of 300 K. The results show that the radial velocity and temperature distributions are less affected by the showerhead aperture than by the radial pressure distribution. Increasing the showerhead aperture improves the radial velocity and temperature distribution uniformity, while decreasing the showerhead aperture improves the radial pressure distribution uniformity.

Key words

vacuum chamber / showerhead / millipore / direct simulation Monte Carlo (DSMC) / numerical simulation / unstructured grid

Cite this article

EndNote

Ris (Procite)

Bibtex

Download Citations

WANG Chuncai, CHENG Jia, JI Linhong, LU Yijia, SUN Yuchun, LIN Jia. 2-D DSMC algorithm based on Delaunay triangles[J]. Journal of Tsinghua University(Science and Technology). 2015, 55(10): 1079-1086,1097 https://doi.org/10.16511/j.cnki.qhdxxb.2015.22.010

References

[1] Quirk M, Serda J. 半导体制造技术 [M]. 韩郑生, 译. 北京: 电子工业出版社, 2004.Quirk M, Serda J. Semiconductor Manufacturing Technology [M]. HAN Zhengsheng, trans. Beijing: Publishing House of Electronics Industry, 2004. (in Chinese)
[2] 谢欣峰. 高密度电浆用于十二吋晶圆0.22微米线宽铝内连线结构高选择比蚀刻研究 [D]. 台中: 逢甲大学, 2007.XIE Xinfeng. High Density Plasma Used in High Sensitivity Etching of Aluminum Inner Connection Structure with 0.22 Micrometer CD of 12 Inches Wafer [D]. Taichung: Feng Chia University, 2007. (in Chinese)
[3] Chen I. Mass transfer analyses of the plasma deposition process [J]. Thin Solid Films, 1983, 101(1): 41-53.
[4] Caquineau H, Despax B. Influence of the reactor design in the case of silicon nitride PECVD [J]. Chemical Engineering Science, 1997, 52(17): 2901-2914.
[5] Nienhuis G J, Goedheer W. Modelling of a large scale reactor for plasma deposition of silicon [J]. Plasma Sources & Technology, 1999, 8(2): 295-298.
[6] Sobbia R, Sansonnens L, Bondkowski J. Uniformity study in large-area showerhead reactors [J]. Journal of Vacuum Science & Technology, 2005, 23(4): 927-932.
[7] Sansonnens L, Howling A A, Hollenstein C. A gas flow uniformity study in large-area showerhead reactors for RF plasma deposition [J]. Plasma Sources & Technology, 2000, 9(2): 205-209.
[8] Howling A A, Legradic B, Chesaux M, et al. Plasma deposition in an ideal showerhead reactor: A two-dimensional analytical solution [J]. Plasma Sources & Technology, 2012, 21: 0150051-01-0150051-15.
[9] Bird G A. Approach to translational equilibrium in a rigid sphere gas [J]. Physics of Fluids, 1963, 6(10): 1518-1519.
[10] Celenligl M C, Moss J N. Application of the DSMC method to hypersonic flow about a delta wing [C]// International Symposium on Rarefied Gas Dynamics. Aachen, Germany, 1990.
[11] Wagner W. A convergence proof for Bird's direct simulation Monte Carto method for the Boltzmann equation [J]. Journal of Statistical Physics, 1992, 66(3/4): 1011-1044.
[12] Pulvircnti M, Wagner W, Zavelani M B. Convergence of particle schemes for the Boltzmann equation [J]. European Journal of Mechanics B: Fluids, 1994, 13(3): 339-351.
[13] Bird G A. Molecular Gas Dynamics [M]. Oxford: Clarendon Press, 1976.
[14] Bird G A. Molecular Gas Dynamics and the Direct Simulation of Gas Flows [M]. Oxford: Clarendon Press, 1994.
[15] Bird G A. Application of the DSMC Method to the Full Shuttle Geometry [R]. AIAA 90-1692. Washington DC, 1990.
[16] 樊菁, 彭世锍, 刘宏立, 等. DSMC 位置元方案中的表面元的程序标识及分子表面反射的确定论判据 [J]. 力学学报, 1999, 31(6): 671-676.FAN Jing, PENG Shiliu, LIU Hongli, et al. Program making of surface elements and the determination of molecular surface reflection in position element algorithm of DSMC method [J]. Chinese Journal of Theoretical and Applied Mechanics, 1999, 31(6): 671-676. (in Chinese)
[17] 樊菁, 刘宏立, 沈青, 等. 直接统计模拟位置元方法中的分子表面反射确定论判据 [J]. 空气动力学学报, 2000, 18(2): 180-187.FAN Jing, LIU Hongli, SHEN Qing, et al. A molecular reflection deterministic criterion used in the position element algorithm of direct statistical simulation [J]. Acta Aerodynamica Sinica, 2000, 18(2): 180-187. (in Chinese)
[18] 沈青, 樊菁, 胡振华, 等. 过渡领域三维绕流直接统计模拟中的一种新方案 [J]. 空气动力学学报, 1996, 14(3): 295-303.SHEN Qing, FAN Jing, HU Zhenghua, et al. A new version of position element space discretization in direct statistical simulation of three-dimensional flow in transitional regime [J]. Acta Aerodynamica Sinica, 1996, 14(3): 295-303. (in Chinese)
[19]Stanton T K. Simple approximates formulas for back scattering of sound by spherical and elongated objects [J]. Journal of the Acoustical Society in America, 1989, 86 (4): 1499-1510.
[20]Wilmoth R G, LeBeau G J, Carlson A B. DSMC Grid Methodologies for Computing Low Density, Hypersonic Flows about Reusable Launch Vehicles [R]. AIAA 96-1812. Washington DC, 1996.
[21]Laux M, Fasoulas S, Messerschmid E W. Development of a DSMC Code on Planar Unstructured Grids with Automatic Grid Adaptation [R]. AIAA 95-2053. Washington DC, 1996.
[22]王学德, 伍贻兆, 夏健. 二维非结构网格DSMC方法实现及其应用 [J]. 南京航空航天大学学报, 2004, 36(6): 704-707.WANG Xuede, WU Yizhao, XIA Jian. Implementation of 2-D unstructured DSMC method and its application [J]. Journal of Nanjing University of Aeronautics & Astronautics, 2004, 36(6): 704-707. (in Chinese)
[23]Berg M. Computational Geometry: Algorithms and Applications [M]. 北京: 世界图书出版公司, 2013.Berg M. Computational Geometry: Algorithms and Applications [M]. Beijing: World Publishing Corporation, 2013.
[24]刘汝佳, 黄亮. 算法艺术与信息学竞赛 [M]. 北京: 清华大学出版社, 2004.LIU Rujia, HUANG Liang. Algorithm Art and Informatics Competitions [M]. Beijing: Tsinghua University Press, 2004. (in Chinese)
[25]Ikegawa M, Kobayashi J. Development of a rarefied flow simulator using the direct simulation Monte Carlo method [J]. JSME International Journal, 1990, 30: 463- 467.
[26]Nance R P, Hash D B, Hassan H A. Role of Boundary Conditions in Monte Carlo Simulation of MEMS Devices [R]. AIAA 97-0375. Washington DC, 1997.[27]Liou W W, Fang Y. Implicit boundary conditions for direct simulation Monte Carlo method in MEMS flow predictions [J]. Computer Modeling in Engineering & Science, 2000, 4: 119-128.
[28]王沫然, 王金库, 李志信. DSMC方法的压力边界条件实现 [J]. 计算物理, 2004, 21(4): 316-320.WANG Moran, WANG Jinku, LI Zhixin. New implement of pressure boundary conditions for DSMC [J]. Chinese Journal of Computational Physics, 2004, 21(4): 316-320. (in Chinese)
[29]Alexeenko A A, Levin D A, Fedosov D A, et al. Performance analysis of microthrusters based on coupled thermal-fluid modeling and simulation [J]. Journal of Propulsion and Power, 2005, 21(1): 95-101.
[30]Hash D B, Mihopoulos T, Govindan T R, et al. Characterization of showerhead performance at low pressure [J]. Journal of Vacuum Science & Technology, 2000, 18(6): 2808-2813.
[31]麦克劳德. 光学薄膜技术 [M]. 周九林, 尹树百, 译. 北京: 国防工业出版社, 1974.Macleod H A. Optical Thin Film Technology [M]. ZHOU Jiulin, YIN Shubai, trans. Beijing: National Defence Industry Press, 1974.(in Chinese)