三维复杂界面非结构网格N-S/DSMC耦合方法

陈飞同; 王学德

doi:10.13224/j.cnki.jasp.20240329

三维复杂界面非结构网格N-S/DSMC耦合方法

doi: 10.13224/j.cnki.jasp.20240329

陈飞同,
王学德^*, ,

南京理工大学能源与动力工程学院，南京 210094

详细信息

作者简介:
陈飞同（1999-），男，硕士生，主要从事稀薄气体动力学研究

通讯作者:
王学德（1977-），男，副教授，博士，主要从事稀薄气体动力学研究。E-mail：wangxuede2000@njust.edu.cn

中图分类号: V411.4
计量
- 文章访问数: 423
- HTML浏览量: 313
- PDF量: 35
- 被引次数: 0
出版历程
- 收稿日期: 2024-05-21
- 网络出版日期: 2025-01-09

N-S/DSMC coupling method using three-dimensional unstructured mesh for complex interfaces

CHEN Feitong,
WANG Xuede^{*
, ,}

School of Energy and Power Engineering，Nanjing University of Science and Technology，Nanjing 210094，China

摘要

摘要:
针对近连续流过渡区多尺度绕流问题，学界发展了N-S（Navier-Stokes）/DSMC（direct simulation Monte Carlo）耦合方法，大多数此类求解器面临着耦合界面波动失稳的难题，因此对复杂外形和耦合界面的通用性值得重点研究。鉴于非结构网格面向复杂外形高度的贴体性、优良适应性以及工程领域对过渡流区高效通用型计算方法的需求，提出并实现了一套三维复杂界面四面体非结构网格N-S/DSMC耦合方法用于模拟高超声速过渡流。该方法使用局部克努森数作为连续失效参数划分连续/稀薄区域，并生成三维复杂N-S/DSMC耦合界面，沿分界面两侧分别推进一层或多层界面信息传递单元，基于边界状态法进行信息耦合。该耦合方法无需对复杂不规则分界面作光滑和修型处理，具备对复杂过渡流区工程问题数值模拟的通用性。分别对三维高超声速圆球和钝锥绕流进行模拟，数值结果显示：与参考文献中的DSMC方法相比，激波处数值和壁面特征值基本一致，最大误差不超过8%，但计算效率分别提高了1.74倍和2.28倍，验证了该耦合方法的正确性和高效性。
- N-S/DSMC耦合方法 /
- 四面体非结构网格 /
- 高超声速流 /
- 过渡流区 /
- 复杂界面耦合方法
Abstract:
In view of the highly conformable and adaptive nature of unstructured mesh to complex geometries and the demand in engineering for efficient and versatile computational methods, a new implementation of the N-S (Navier-Stokes)/DSMC (direct simulation Monte Carlo) coupling method using tetrahedral unstructured mesh for complex interfaces was presented with the aim of resolving the fluctuations in the positioning of the coupling interfaces in most N-S/DSMC coupling method for multi-scale transitional flows in the near-continuum regime. This implementation utilized local Knudsen number as a continuum breakdown parameter to partition the continuum/rarefied regions and generate three-dimensional complex N-S/DSMC coupling interfaces. Along each side of the interface, one or multiple layers of interface information transfer meshes were advanced, and information coupling was achieved by the state-based approach. According to this method, there was no need for smoothing or reshaping treatments applied to complex irregular interfaces, with the general applicability to numerical simulations of complex transitional flow regions. Simulations of three-dimensional hypersonic flow around a sphere and a blunt cone showed that, compared with the DSMC method, the shock wave and wall characteristics were in close agreement, with a maximum error of less than 8%. At the same time, computational efficiency was improved by 1.74 and 2.28 times, respectively, validating the method’s accuracy and efficiency.
- N-S/DSMC coupling method /
- tetrahedral unstructured mesh /
- hypersonic flow /
- transitional flow /
- complex interfaces coupling method

HTML

图 1 三维分子搜索示意图

Figure 1. Diagram of three-dimensional molecular search

下载: 全尺寸图片幻灯片

图 2 信息传递界面的生成示意图

Figure 2. Diagram of information transfer interface generation

下载: 全尺寸图片幻灯片

图 3 界面信息传递单元生成示意图

Figure 3. Diagram of interface information transfer meshes generation

下载: 全尺寸图片幻灯片

图 4 整体网格结构示意图

Figure 4. Diagram of the overall mesh structure

下载: 全尺寸图片幻灯片

图 5 N-S界面信息传递单元传递信息示意图

Figure 5. Diagram of information transfer among N-S interface information transfer meshes

下载: 全尺寸图片幻灯片

图 6 DSMC界面信息传递单元传递信息示意图

Figure 6. Diagram of information transfer among DSMC interface information transfer meshes

下载: 全尺寸图片幻灯片

图 7 圆球绕流耦合界面与流区网格分布

Figure 7. Distribution map of coupling interfaces and meshes for the flow field around a sphere

下载: 全尺寸图片幻灯片

图 8 计算结果比较

Figure 8. Comparison of calculation results

下载: 全尺寸图片幻灯片

图 9 圆球壁面特征值比较

Figure 9. Comparison of wall characteristic coefficients for a sphere surface

下载: 全尺寸图片幻灯片

图 10 钝锥绕流耦合界面与流区网格分布

Figure 10. Distribution map of coupling interfaces and meshes for the flow field around a blunt cone

下载: 全尺寸图片幻灯片

图 11 钝锥壁面特征量比较

Figure 11. Comparison of wall characteristics for a blunt cone surface

下载: 全尺寸图片幻灯片

表 1 两种方法总计算耗时

Table 1. Time consumption of the two methods h

方法圆球耗时钝锥耗时

耦合算法 41.23 27.3

DSMC 71.55 62.12

下载: 导出CSV

参考文献(21)

[1]	李志辉, 梁杰, 李中华, 等. 跨流域空气动力学模拟方法与返回舱再入气动研究[J]. 空气动力学学报, 2018, 36(5): 826-847. LI Zhihui, LIANG Jie, LI Zhonghua, et al. Simulation methods of aerodynamics covering various flow regimes and applications to aerodynamic characteristics of re-entry spacecraft module[J]. Acta Aerodynamica Sinica, 2018, 36(5): 826-847. (in Chinese doi: 10.7638/kqdlxxb-2018.0121 LI Zhihui, LIANG Jie, LI Zhonghua, et al. Simulation methods of aerodynamics covering various flow regimes and applications to aerodynamic characteristics of re-entry spacecraft module[J]. Acta Aerodynamica Sinica, 2018, 36(5): 826-847. (in Chinese) doi: 10.7638/kqdlxxb-2018.0121
[2]	HASH D B, HASSAN H A. Assessment of schemes for coupling Monte Carlo and Navier-Stokes solution methods[J]. Journal of Thermophysics and Heat Transfer, 1996, 10(2): 242-249. doi: 10.2514/3.781
[3]	HASH D, HASSAN H, HASH D, et al. Two-dimensional coupling issues of hybrid DSMC/Navier-Stokes solvers[R]. AIAA 1997-2507, 1997.
[4]	SCHWARTZENTRUBER T E, BOYD I D. A hybrid particle-continuum method applied to shock waves[J]. Journal of Computational Physics, 2006, 215(2): 402-416. doi: 10.1016/j.jcp.2005.10.023
[5]	SCHWARTZENTRUBER T E, SCALABRIN L C, BOYD I D. A modular particle-continuum numerical method for hypersonic non-equilibrium gas flows[J]. Journal of Computational Physics, 2007, 225(1): 1159-1174. doi: 10.1016/j.jcp.2007.01.022
[6]	ROVEDA R, GOLDSTEIN D B, VARGHESE P L. Hybrid Euler/direct simulation Monte Carlo calculation of unsteady slit flow[J]. Journal of Spacecraft and Rockets, 2000, 37(6): 753-760. doi: 10.2514/2.3647
[7]	SCHWARTZENTRUBER T, SCALABRIN L, BOYD I. Hybrid particle-continuum simulations of low Knudsen number hypersonic flows[R]. AIAA 2007-3892, 2007.
[8]	CAI Guobiao, TANG Zhenyu, WEI Xiaoqing. Investigation on the interface smoothing of coupled N-S/DSMC method using image processing filters[J]. Advances in Astronautics Science and Technology, 2018, 1(2): 221-227. doi: 10.1007/s42423-018-0023-3
[9]	STEGER J L, DOUGHERTY F C, BENEK J A. A chimera grid scheme[C]// ASME Mini-Symposium on Advances in Grid Generation. Houston, TX, US: ASME, 1982: 59-69.
[10]	WU Ziniu, XU Shanshu, GAO Bo, et al. Review of numerical computation of compressible flows with artificial interfaces[J]. Computers & Fluids, 2007, 36(10): 1657-1679.
[11]	XU Xiao, WANG Xuede, ZHANG Mu, et al. A parallelized hybrid N-S/DSMC-IP approach based on adaptive structured/unstructured overlapping grids for hypersonic transitional flows[J]. Journal of Computational Physics, 2018, 371: 409-433. doi: 10.1016/j.jcp.2018.05.021
[12]	许啸. 连续/稀薄混合流场动态重叠网格CFD/DSMC耦合方法研究[D]. 南京: 南京理工大学, 2015. XU Xiao. The research on dynamic overset grid coupled CFD/DSMC method of the hybrid continuum/rarefied flow field[D]. Nanjing: Nanjing University of Science and Technology, 2015. (in Chinese XU Xiao. The research on dynamic overset grid coupled CFD/DSMC method of the hybrid continuum/rarefied flow field[D]. Nanjing: Nanjing University of Science and Technology, 2015. (in Chinese)
[13]	WU J S, LIAN Y Y, CHENG G, et al. Development and verification of a coupled DSMC-NS scheme using unstructured mesh[J]. Journal of Computational Physics, 2006, 219(2): 579-607. doi: 10.1016/j.jcp.2006.04.013
[14]	BLAZEK J. Computational fluid dynamics: principles and applications[M]. 3rd. ed. Amsterdam: Butterworth Heinemann, 2015.
[15]	BIRD G A. Molecular gas dynamics and the direct simulation of gas flows[M]. Oxford: Clarendon Press, 1994.
[16]	王学德. 高超声速稀薄气流非结构网格DSMC及并行算法研究[D]. 南京: 南京航空航天大学, 2006. WANG Xuede. DSMC method on unstructured grids for hypersonic rarefied gas flow and its parallelization[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2006. (in Chinese WANG Xuede. DSMC method on unstructured grids for hypersonic rarefied gas flow and its parallelization[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2006. (in Chinese)
[17]	BOYD I D, CHEN Gang, CANDLER G V. Predicting failure of the continuum fluid equations in transitional hypersonic flows[J]. Physics of Fluids, 1995, 7(1): 210-219. doi: 10.1063/1.868720
[18]	WANG Wenlan, BOYD I D. Predicting continuum breakdown in hypersonic viscous flows[J]. Physics of Fluids, 2003, 15(1): 91-100. doi: 10.1063/1.1524183
[19]	SUN Quanhua, BOYD I D. Evaluation of macroscopic properties in the direct simulation Monte Carlo method[J]. Journal of Thermophysics and Heat Transfer, 2005, 19(3): 329-335. doi: 10.2514/1.12542
[20]	许啸, 王学德, 谭俊杰. 三维高超声速稀薄气体流场通量分裂式DSMC-IP方法研究[J]. 工程力学, 2016, 33(8): 238-247. XU Xiao, WANG Xuede, TAN Junjie. Research on a flux-splitting method in three dimensional hypersonic rarefied flow field[J]. Engineering Mechanics, 2016, 33(8): 238-247. (in Chinese doi: 10.6052/j.issn.1000-4750.2014.12.1091 XU Xiao, WANG Xuede, TAN Junjie. Research on a flux-splitting method in three dimensional hypersonic rarefied flow field[J]. Engineering Mechanics, 2016, 33(8): 238-247. (in Chinese) doi: 10.6052/j.issn.1000-4750.2014.12.1091
[21]	黄飞, 陈智, 程晓丽, 等. 一种基于自适应碰撞距离的DSMC虚拟子网格方法[J]. 空气动力学学报, 2014, 32(4): 506-510. HUANG Fei, CHEN Zhi, CHENG Xiaoli, et al. A virtual sub-cells technique with transient adaptive collision distance for the DSMC method[J]. Acta Aerodynamica Sinica, 2014, 32(4): 506-510. (in Chinese doi: 10.7638/kqdlxxb-2012.0161 HUANG Fei, CHEN Zhi, CHENG Xiaoli, et al. A virtual sub-cells technique with transient adaptive collision distance for the DSMC method[J]. Acta Aerodynamica Sinica, 2014, 32(4): 506-510. (in Chinese) doi: 10.7638/kqdlxxb-2012.0161