基于速度空间非结构网格和守恒修正的改进离散速度方法

doi:10.7527/S1000-6893.2022.27033

临近空间高速飞行器数值模拟技术专栏

本期目录 | 过刊浏览 | 高级检索

前一篇 | 后一篇

基于速度空间非结构网格和守恒修正的改进离散速度方法

杨鲤铭, 吴杰, 董昊, 杜银杰

南京航空航天大学航空学院, 南京 210016

收稿日期:2022-02-14 修回日期:2022-03-17 发布日期:2022-03-30
通讯作者: 杨鲤铭,E-mail:lmyang@nuaa.edu.cn E-mail:lmyang@nuaa.edu.cn
基金资助:
江苏省自然科学基金（BK20210273）；江苏高校优势学科建设工程资助项目（PAPD）

Improved discrete velocity method based on unstructured mesh in velocity space and conservative correction

YANG Liming, WU Jie, DONG Hao, DU Yinjie

College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China

Received:2022-02-14 Revised:2022-03-17 Published:2022-03-30
Supported by:
Natural Science Foundation of Jiangsu Province (BK20210273); Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD)

摘要/Abstract

摘要： 传统离散速度方法在求解跨流域流动问题时，通常只求解动理学模型方程，即Boltzmann-BGK方程。与传统方法不同，改进离散速度方法同步求解了动理学模型方程和相应的宏观伴随方程。通过这种方式，可以将分子碰撞影响考虑到宏观伴随方程的通量计算中，同时宏观方程预估得到的结果可以用于预估平衡态，从而实现Boltzmann-BGK方程的全隐式离散。这两点改进可以有效克服传统方法在连续和近连续流区域计算效率低、计算精度差的缺陷。为了进一步减少速度空间的网格量和避免数值求积时的Runge现象，采用了非结构网格结合矩形律来离散速度空间并引进守恒修正来强制满足相容性条件。算例测试表明，采用速度空间非结构网格和守恒修正可以有效减少改进离散速度方法的计算量和内存花销。

关键词: 离散速度方法, 全隐式格式, 非结构网格, 守恒修正, 计算效率

Abstract: Conventional Discrete Velocity Method (DVM) only resolves the kinetic equation (i.e., Boltzmann-BGK equation) in each time step for fluid flow problem simulation in all flow regimes. Different from the conventional method, the Improved Discrete Velocity Method (IDVM) solves both the kinetic equation and the corresponding macroscopic governing equations so that the collisional effect at the cell interface can be involved in the calculation of numerical flux and the fully implicit discretization of the kinetic equation realized by the predicted results of macroscopic governing equations. These two improvements can effectively overcome the defects of low efficiency and poor accuracy of the conventional DVM in the continuum flow regime. To further reduce the number of discrete points in the velocity space and avoid the Runge phenomenon for numerical quadrature, the unstructured mesh combined with the rectangle rule is utilized to discretize the velocity space, and the conservative correction introduced to enforce the compatibility condition. Numerical results show that these strategies can significantly reduce the computational cost and memory consumption of the present method.

Key words: discrete velocity method, fully implicit scheme, unstructured mesh, conservative correction, computational efficiency

中图分类号:

V211.3

杨鲤铭, 吴杰, 董昊, 杜银杰. 基于速度空间非结构网格和守恒修正的改进离散速度方法[J]. 航空学报, 2022, 43(12): 627033-627033.

YANG Liming, WU Jie, DONG Hao, DU Yinjie. Improved discrete velocity method based on unstructured mesh in velocity space and conservative correction[J]. ACTA AERONAUTICAET ASTRONAUTICA SINICA, 2022, 43(12): 627033-627033.

参考文献

[1] RANA A S, MOHAMMADZADEH A, STRUCHTRUP H. A numerical study of the heat transfer through a rarefied gas confined in a microcavity[J]. Continuum Mechanics and Thermodynamics, 2015, 27(3):433-446.
[2] 王卫东. 面向MEMS设计的微流体流动特性研究[D]. 西安:西安电子科技大学, 2007. WANG W D. Study on flow characteristics of microfluidics for MEMS design[D]. Xi'an:Xidian University, 2007(in Chinese).
[3] 樊菁. 稀薄气体动力学:进展与应用[J]. 力学进展, 2013, 43(2):185-201. FAN J. Rarefied gas dynamics:Advances and applications[J]. Advances in Mechanics, 2013, 43(2):185-201(in Chinese).
[4] 梁杰, 李志辉, 李齐, 等. 返回舱再入跨流域气动及配平特性数值研究[J]. 空气动力学学报, 2018, 36(5):848-855. LIANG J, LI Z H, LI Q, et al. Numerical simulation of aerodynamic and trim characteristics across different flow regimes for reentry module[J]. Acta Aerodynamica Sinica, 2018, 36(5):848-855(in Chinese).
[5] 李锦, 耿湘人, 陈坚强, 等. DSMC量子动理学模型在火星再入流动中的应用[J]. 航空学报, 2020, 41(7):123240. LI J, GENG X R, CHEN J Q, et al. Application of DSMC quantum kinetic model in re-entry flow of Mars[J]. Acta Aeronautica et Astronautica Sinica, 2020, 41(7):123240(in Chinese).
[6] BROADWELL J E. Study of rarefied shear flow by the discrete velocity method[J]. Journal of Fluid Mechanics, 1964, 19(3):401-414.
[7] YANG L M, SHU C, WU J, et al. Comparative study of discrete velocity method and high-order lattice Boltzmann method for simulation of rarefied flows[J]. Computers & Fluids, 2017, 146:125-142.
[8] YANG J Y, HUANG J C. Rarefied flow computations using nonlinear model Boltzmann equations[J]. Journal of Computational Physics, 1995, 120(2):323-339.
[9] MIEUSSENS L. On the asymptotic preserving property of the unified gas kinetic scheme for the diffusion limit of linear kinetic models[J]. Journal of Computational Physics, 2013, 253:138-156.
[10] WANG P, HO M T, WU L, et al. A comparative study of discrete velocity methods for low-speed rarefied gas flows[J]. Computers & Fluids, 2018, 161:33-46.
[11] KUDRYAVTSEV A N, SHERSHNEV A A. A numerical method for simulation of microflows by solving directly kinetic equations with WENO schemes[J]. Journal of Scientific Computing, 2013, 57(1):42-73.
[12] DIAZ M A, CHEN M H, YANG J Y. High-order conservative asymptotic-preserving schemes for modeling rarefied gas dynamical flows with Boltzmann-BGK equation[J]. Communications in Computational Physics, 2015, 18(4):1012-1049.
[13] MIEUSSENS L. Discrete velocity model and implicit scheme for the BGK equation of rarefied gas dynamics[J]. Mathematical Models and Methods in Applied Sciences, 2000, 10(8):1121-1149.
[14] YANG L M, SHU C, YANG W M, et al. An improved discrete velocity method (DVM) for efficient simulation of flows in all flow regimes[J]. Physics of Fluids, 2018, 30(6):062005.
[15] 李志辉, 彭傲平, 方方, 等. 跨流域高超声速绕流环境Boltzmann模型方程统一算法研究[J]. 物理学报, 2015, 64(22):224703. LI Z H, PENG A P, FANG F, et al. Gas-kinetic unified algorithm for hypersonic aerothermodynamics covering various flow regimes solving Boltzmann model equation[J]. Acta Physica Sinica, 2015, 64(22):224703(in Chinese).
[16] 李志辉, 彭傲平, 马强, 等. 大型航天器离轨再入气动融合结构变形失效解体落区数值预报与应用[J]. 载人航天, 2020, 26(4):403-417. LI Z H, PENG A P, MA Q, et al. Study and application of deformation failure disintegration algorithm for large-scale spacecraft reentry aerodynamic coupling structure[J]. Manned Spaceflight, 2020, 26(4):403-417(in Chinese).
[17] LI Z H, HU W Q, PENG A P, et al. Gas-kinetic unified algorithm for plane external force-driven flows covering all flow regimes by modeling of Boltzmann equation[J]. International Journal for Numerical Methods in Fluids, 2020, 92(8):922-949.
[18] XU K, HUANG J C. A unified gas-kinetic scheme for continuum and rarefied flows[J]. Journal of Computational Physics, 2010, 229(20):7747-7764.
[19] 刘沙, 王勇, 袁瑞峰, 等. 统一气体动理学方法研究进展[J]. 气体物理, 2019, 4(4):1-13. LIU S, WANG Y, YUAN R F, et al. Advance in unified methods based on gas-kinetic theory[J]. Physics of Gases, 2019, 4(4):1-13(in Chinese).
[20] 张华波, 周瑞睿, 李思达, 等. 基于气体动理学方法的多尺度稳态热辐射计算[J]. 航空学报, 2021, 42(S1):726411. ZHANG H B, ZHOU R D, LI S D, et al. Multiscale steady-state thermal radiation calculation based on gas kinetics method[J]. Acta Aeronautica et Astronautica Sinica, 2021, 42(S1):726411(in Chinese).
[21] GUO Z L, XU K, WANG R J. Discrete unified gas kinetic scheme for all Knudsen number flows:Low-speed isothermal case[J]. Physical Review E, 2013, 88(3):033305.
[22] 姚博, 张创, 郭照立. 考虑分子转动自由度的离散统一气体动理学格式[J]. 航空学报, 2019, 40(7):122914. YAO B, ZHANG C, GUO Z L. Discrete unified gas kinetic scheme for diatomic gas with rotational degrees of freedom[J]. Acta Aeronautica et Astronautica Sinica, 2019, 40(7):122914(in Chinese).
[23] YANG L M, SHU C, YANG W M, et al. An improved three-dimensional implicit discrete velocity method on unstructured meshes for all Knudsen number flows[J]. Journal of Computational Physics, 2019, 396:738-760.
[24] YANG L M, ZHAO X, SHU C, et al. Parametric reduced order modeling-based discrete velocity method for simulation of steady rarefied flows[J]. Journal of Computational Physics, 2021, 430:110037.
[25] CHEN S Z, GUO Z L, XU K. Simplification of the unified gas kinetic scheme[J]. Physical Review E, 2016, 94(2-1):023313.
[26] 皮兴才, 朱炼华, 李志辉, 等. 基于宏观方程数值本构关系的气体动理论加速收敛方法[J]. 物理学报, 2020, 69(20):287-298. PI X C, ZHU L H, LI Z H, et al. Method of accelerating convergence for gas kinetic algorithm based on digital constitutive relation of macroscopic equations[J]. Acta Physica Sinica, 2020, 69(20):287-298(in Chinese).
[27] FORNBERG B, ZUEV J. The Runge phenomenon and spatially variable shape parameters in RBF interpolation[J]. Computers & Mathematics With Applications, 2007, 54(3):379-398.
[28] SHAKHOV E M. Generalization of the Krook kinetic relaxation equation[J]. Fluid Dynamics, 1972, 3(5):95-96.
[29] YANG L, WANG Y, CHEN Z, et al. Lattice Boltzmann and gas kinetic flux solvers:Theory and applications[M]. Hackensack:World Scientific Publishing Co. Pte. Ltd., 2020.
[30] YUAN R F, ZHONG C W. A conservative implicit scheme for steady state solutions of diatomic gas flow in all flow regimes[J]. Computer Physics Communications, 2020, 247:106972.
[31] CHEN J F, LIU S, WANG Y, et al. Conserved discrete unified gas-kinetic scheme with unstructured discrete velocity space[J]. Physical Review E, 2019, 100(4-1):043305.
[32] TITAREV V A. Implicit high-order method for calculating rarefied gas flow in a planar microchannel[J]. Journal of Computational Physics, 2012, 231(1):109-134.
[33] YANG L M, SHU C, YANG W M, et al. Numerical investigation on performance of three solution reconstructions at cell interface in DVM simulation of flows in all Knudsen number regimes[J]. International Journal for Numerical Methods in Fluids, 2019, 90(11):545-563.
[34] GHIA U, GHIA K N, SHIN C T. High-Re solutions for imcompressible flow using the Navier-Stokes equations and a multigrid method[J]. Journal of Computational Physics, 1982, 48(3):387-411.
[35] HUANG J C,XU K,YU P B.A unified gas-kinetic scheme for continuum and rarefield flows Ⅱ:Multi-dimensional cases[J]. Communications in Computational Physics, 2012,12(3):662-690.
[36] JIANG D W, MAO M L, LI J, et al. An implicit parallel UGKS solver for flows covering various regimes[J]. Advances in Aerodynamics, 2019, 1(1):8.
[37] LI Z H, ZHANG H X. Gas-kinetic numerical studies of three-dimensional complex flows on spacecraft re-entry[J]. Journal of Computational Physics, 2009, 228(4):1116-1138.

E-mail：hkxb@buaa.edu.cn

关于我们

期刊社服务

专业学科

封面文章

友情链接

主管单位：中国科学技术协会主办单位：中国航空学会北京航空航天大学

基于速度空间非结构网格和守恒修正的改进离散速度方法

Improved discrete velocity method based on unstructured mesh in velocity space and conservative correction

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

本文评价

[1]	张宝振, 王汉平, 徐峰, 吴志青. VSV调节机构的仿真提速和精度补偿措施[J]. 航空学报, 2022, 43(9): 226034-226034.
[2]	陈坚强, 吴晓军, 张健, 李彬, 贾洪印, 周乃春. FlowStar:国家数值风洞(NNW)工程非结构通用CFD软件[J]. 航空学报, 2021, 42(9): 625739-625739.
[3]	任玉新, 王乾, 潘建华, 章雨思, 黄乾旻. 构建非结构网格高精度有限体积方法的新途径[J]. 航空学报, 2021, 42(9): 625783-625783.
[4]	刘君, 魏雁昕, 陈洁. 基于非结构网格有限差分法的扎染算法[J]. 航空学报, 2021, 42(7): 124557-124557.
[5]	刘君, 陈洁, 韩芳. 基于离散等价方程的非结构网格有限差分法[J]. 航空学报, 2020, 41(1): 123248-123248.
[6]	谭伟伟, 颜洪, 聂智军, 马涂亮, 梁益华. 大型客机涡扇发动机动力特性模拟[J]. 航空学报, 2019, 40(1): 522428-522428.
[7]	郭中州, 何志强, 赵文文, 陈伟芳. 高效非结构网格变形与流场插值方法[J]. 航空学报, 2018, 39(12): 122411-122411.
[8]	贡伊明, 刘战合, 刘溢浪, 张伟伟. 时间谱方法中的高效GMRES算法[J]. 航空学报, 2017, 38(7): 120894-120894.
[9]	刘君, 刘瑜, 陈泽栋. 非结构变形网格和离散几何守恒律[J]. 航空学报, 2016, 37(8): 2395-2407.
[10]	张扬, 张来平, 赫新, 邓小刚. 基于自适应混合网格的脱体涡模拟[J]. 航空学报, 2016, 37(12): 3605-3614.
[11]	崔鹏程, 邓有奇, 唐静, 李彬. 基于伴随方程的网格自适应及误差修正[J]. 航空学报, 2016, 37(10): 2992-3002.
[12]	张扬, 张来平, 赫新, 邓小刚. 基于非结构/混合网格的脱体涡模拟算法[J]. 航空学报, 2015, 36(9): 2900-2910.
[13]	燕振国, 刘化勇, 毛枚良, 马燕凯, 朱华君. 三阶HWCNS的构造及其在高超声速流动中的应用[J]. 航空学报, 2015, 36(5): 1460-1470.
[14]	李彬, 唐静, 邓有奇, 张耀冰. 并行的多重网格方法在离散伴随优化中的应用[J]. 航空学报, 2014, 35(8): 2091-2101.
[15]	黄飞, 沈清, 程晓丽, 赵波. 一种DSMC分子仿真下的权因子预定义方法[J]. 航空学报, 2014, 35(8): 2174-2181.