面向三维激波问题的装配方法

doi:10.7527/S1000-6893.2020.24141

流体力学与飞行力学

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

前一篇 | 后一篇

面向三维激波问题的装配方法

邹东阳¹, 林敬周¹, 黄洁¹, 刘君²

1. 中国空气动力研究与发展中心超高速空气动力研究所, 绵阳 621000;
2. 大连理工大学航空航天学院, 大连 116024

收稿日期:2020-04-27 修回日期:2020-07-13 发布日期:2020-08-03
通讯作者: 刘君 E-mail:liujun65@dlut.edu.cn

Fitting algorithms for three dimensional flows with shock waves

ZOU Dongyang¹, LIN Jingzhou¹, HUANG Jie¹, LIU Jun²

1. Hypervelocity Aerodynamics Institute, China Aerodynamics Research and Development Center, Mianyang 621000, China;
2. School of Aeronautics and Astronautics, Dalian University of Technology, Dalian 116024, China

Received:2020-04-27 Revised:2020-07-13 Published:2020-08-03

摘要/Abstract

摘要： 给出了一种基于非结构动网格技术的三维激波装配方法。在该方法中，三维激波面由被标记为激波属性的网格点连接构成，标记为激波属性的网格点称为激波点。激波点具有两组参数分别代表激波的上下游，利用激波点上下游参数求解R-H关系式获得激波点运动速度。非结构动网格技术的使用允许激波大幅度运动，降低了对初始激波位置的要求。通过引入网格属性定义避免了对计算网格进行分区，增加了装配激波的灵活性。通过球柱体绕流问题验证了该三维装配方法的合理性，针对三维激波装配中比较困难的交点装配问题，通过对三维激波反射以及三维激波相交等算例进行研究找到了可用的三维激波交点运动速度的确定方法，保证了激波运动过程中交点运动与流场求解之间的相容性，获得了相应的装配结果。

关键词: 三维激波装配, 三维激波相交, 非结构动网格, 有限体积法, 超声速流动

Abstract: A three dimensional shock-fitting technique based on unstructured dynamic grids is proposed in this work. In this algorithm, the shock front comprises a series of grid nodes labeled as shock points which have two states, with one representing the upstream of the shock wave and the other the downstream. The R-H relations are solved using these two states on a shock point to obtain the velocity of the shock point. The use of unstructured dynamic grids enables the shock front to move in a large range, decreasing the requirement of the initial shock position. The shocks are labeled by the definition of grid nodes rather than using internal boundaries between different subdomains to improve the flexibility of fitting shock waves. The reliability of the proposed algorithm is proven by a test case of hypersonic flow past a hemisphere-cylinder, followed by the study of three-dimensional shock reflection and shock-shock interactions to solve the relatively complex shock interaction problems in three dimensional shock-fitting. To guarantee the compatibility of the shock point motion and the flow field, usable methods to determine the velocity of three dimensional shock interaction points are obtained, therefore ensuring the achievement of convergence shock-fitting results by this shock-fitting technique.

Key words: three dimensional shock-fitting, three dimensional shock interactions, unstructured dynamic grids, finite volume methods, supersonic flows

中图分类号:

邹东阳, 林敬周, 黄洁, 刘君. 面向三维激波问题的装配方法[J]. 航空学报, 2021, 42(3): 124141-124141.

ZOU Dongyang, LIN Jingzhou, HUANG Jie, LIU Jun. Fitting algorithms for three dimensional flows with shock waves[J]. ACTA AERONAUTICAET ASTRONAUTICA SINICA, 2021, 42(3): 124141-124141.

参考文献

[1] 童秉纲, 孔祥言, 邓国华. 气体动力学[M]. 北京:高等教育出版社, 2011:86-169. TONG B G, KONG X Y, DENG G H. Gas dynamics[M]. Beijing:Higher Education Press, 2011:86-169(in Chinese).
[2] 吴子牛, 白晨媛, 李娟, 等. 高超声速飞行器流动特征分析[J]. 航空学报, 2014, 36(1):58-85. WU Z N, BAI C J, et al. Analysis of flow characteristics for hypersonic vehicle[J]. Acta Aeronautica et Astonautica Sinica, 2014, 36(1):58-85(in Chinese).
[3] RICHTMYER R D, MORTON K W. Difference methods for initial-value problems[M]. New York:Interscience, 1958:90-230.
[4] HARTEN A. High resolution schemes for hyperbolic conservation law[J]. Journal of Computational Physics, 1983, 49:357-393.
[5] YEE H C, KLOPFER G H, MONTAGNE J L. High resolution shock-capturing schemes for inviscid and viscous hypersonic flows[J]. Journal of Computational Physics, 1990, 88:31-61.
[6] HARTEN A, ENGQUIST B, OSHER S, et al. Uniformly high order essentially non-oscillatory schemes[J]. Journal of Computational Physics, 1987, 71:231-303.
[7] LIU X D, OHSER S, CHAN T. Weighted essentially non-oscillatory schemes[J]. Journal of Computational Physics, 1994, 115:217-237.
[8] 张涵信. 无震荡无自由参数耗散差分格式[J]. 空气动力学学报, 1988, 2:143-165. ZHANG H X. Non-oscillatory and non-free-parameter dissipation difference scheme[J]. Acta Aerodynamics Sinica, 1988, 2:143-165(in Chinese).
[9] SALAS M D. A shock-fitting primer[M]. New York:CRC Press, 2009:33-34.
[10] LEE T K, ZHONG X L. Spurious numerical oscillations in simulation of supersonic flow using shock-capturing schemes[J]. AIAA Journal, 1999, 37:383-394.
[11] MORETTI G. Thirty-six years of shock-fitting[J]. Computer & Fluids, 2002, 31:719-723.
[12] PARPIA I H, PARIKH P. A solution-adaptive mesh generation method with cell-face orientation control[C]//Aerospace Sciences Meeting and Exhibit,1994-416, 1994:1-13.
[13] ROSENDALE J V. Floating shock fitting via Lagrangian adaptive meshes:Report No. 94-89[R]. Hampton, VA:NASA Langley Research Center,1994.
[14] TREPANIER J Y,REGGIO M,CAMARERO R, et al. A conservative shock fitting method on unstructured grids[J]. Journal of Computational Physics, 1996, 126(2):421-433.
[15] ZAHR M J, PERSSON P O. An optimization based discontinuous Galerkin approach for high-order accurate shock tracking[C]//AIAA SciTech Forum. Reston:AIAA, 2018:1-11.
[16] ZAHR M J, PERSSON P O. An optimization-based approach for high-order accurate discretization of conservation laws with discontinuous solutions[J]. Journal of Computational Physics, 2018, 365:105-134.
[17] MIRHOSEINI M, PERSSON P O, ZAHR M J. A full space solver for optimization-based, high-order shock tracking using a discontinuous Galerkin discretization[C]//AIAA SciTech Forum. Reston:AIAA, 2019:1-9.
[18] CORRIGAN A, KERCHER A D, KESSLER D A. A moving discontinuous Galerkin finite element method for flows with interfaces[J]. International Journal for Numerical Methods in Fluids, 2019, 89:362-406.
[19] CORRIGAN A T, KERCHER A D, KESSLER D A, et al. Application of the moving discontinuous Galerkin method with interface condition enforcement to shocked compressible flows[C]//AIAA AVIATION Forum. Reston:AIAA, 2018:1-10.
[20] CORRIGAN A T, KERCHER A D, KESSLER D A. The moving discontinuous Galerkin method with interface condition enforcement for unsteady three-dimensional flows[C]//AIAA SciTech Forum. Reston:AIAA, 2019:1-13.
[21] PACIORRI R, BONFIGLIOLI A. A shock-fitting technique for 2D unstructured grids[J]. Computers & Fluids, 2009, 38:715-726.
[22] BONFIGLIOLI A, PACIORRI R. Convergence analysis of shock-capturing and shock-fitting solutions on unstructured grids[J]. AIAA Journal, 2014, 52(7):1404-1416.
[23] PACIORRI R, BONFIGLIOLI A. Shock interaction computations on unstructured, two-dimensional grids using a shock-fitting technique[J]. Journal of Computational Physics, 2011, 230:3155-3177.
[24] 邹东阳, 刘君, 邹丽. 可压缩流动激波装配在格心型有限体积中的应用[J]. 航空学报, 2017, 38(12):121363. ZOU D Y, LIU J, ZOU L. Applications of the shock-fitting technique for compressible flow in cell-centered finite volume methods[J]. Acta Aeronautica et Astronautica Sinica, 2017, 38(12):121363(in Chinese).
[25] ZOU D Y, XU C G, DONG H B, et al. A shock-fitting technique for cell-centered finite volume methods on unstructured dynamic meshes[J]. Journal of Computational Physics, 2017, 345:866-882.
[26] 刘君, 邹东阳, 徐春光. 基于非结构动网格的非定常激波装配方法[J]. 空气动力学学报, 2015, 1:10-16. LIU J, ZOU D Y, XU C G. An unsteady shock-fitting technique based on unstructured moving grids[J]. Acta Aerodynamics Sinica, 2015, 1:10-16(in Chinese).
[27] 刘君, 邹东阳, 董海波. 动态间断装配法模拟激波壁面反射[J]. 航空学报, 2016, 37(3):836-846. LIU J, ZOU D Y, DONG H B. A moving discontinuity fitting technique to simulate shock waves impinged on a straight wall[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(3):836-846(in Chinese).
[28] 常思源, 白晓征, 崔小强, 等. 一种改进的非定常激波装配方法[J]. 航空学报, 2020, 41(2):123498. CHANG S Y, BAI X Z, CUI X Q, et al. An improved unsteady shock-fitting algorithm[J]. Acta Aeronautica et Astronautica Sinica, 2020, 41(2) 123498(in Chinese).
[29] 刘君, 陈洁, 韩芳. 基于离散等价方程的非结构网格有限差分法[J]. 航空学报, 2020, 41(1):123248. LIU J, CHEN J, HAN F. Finite difference method for unstructured grid based on discrete equivalent equation[J]. Acta Aeronautica et Astronautica Sinica, 2020, 41(1):123248(in Chinese).
[30] BONFIGLIOLI A, GROTTADAUREA M, PACIORRI R, et al. An unstructured, three-dimensional, shock-fitting solver for hypersonic flows[J]. Computers & Fluids, 2013,73:162-174.
[31] 刘君, 徐春光, 白晓征. 有限体积法和非结构动网格[M]. 北京:科学出版社, 2016:42-163. LIU J, XU C G, BAI X Z. Finite volume methods and unstructured dynamic grids technique[M]. Beijing:Science Press, 2016:42-163(in Chinese).
[32] HOUMTMAN E M, BANNINK W J, TIMMERMAN B H. Experimental and numerical investigation of the high-supersonic flow around an axi-symmetric blunt cylinder-flare model[C]//Proceedings of the 2nd European Symposium on Aerothermodynamics for Space Vehicles, 1995:517-522.

E-mail：hkxb@buaa.edu.cn

关于我们

期刊社服务

专业学科

封面文章

友情链接

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

面向三维激波问题的装配方法

Fitting algorithms for three dimensional flows with shock waves

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

本文评价

[1]	黄依峰, 曾舒华, 江中正, 陈伟芳. 非线性耦合本构在高空侧向喷流中的数值研究[J]. 航空学报, 2022, 43(S2): 8-22.
[2]	刘君, 韩芳, 魏雁昕. 应用维数分裂方法推广MUSCL和WENO格式的若干问题[J]. 航空学报, 2022, 43(3): 125009-125009.
[3]	刘朋欣, 袁先旭, 孙东, 傅亚陆, 李辰. 高温化学非平衡湍流边界层直接数值模拟[J]. 航空学报, 2022, 43(1): 124877-124877.
[4]	任伟杰, 谢文佳, 田正雨, 张烨, 于航. 高超声速数值激波失稳的网格依赖性[J]. 航空学报, 2021, 42(S1): 726376-726376.
[5]	肖光明, 张超, 桂业伟, 杜雁霞, 刘磊, 魏东. 基于TLBM-FVM耦合的飞行器舱内热环境跨尺度预测方法[J]. 航空学报, 2021, 42(9): 625710-625710.
[6]	向星皓, 张毅锋, 袁先旭, 涂国华, 万兵兵, 陈坚强. C-γ-Re_θ高超声速三维边界层转捩预测模型[J]. 航空学报, 2021, 42(9): 625711-625711.
[7]	卢昱锦, 肖天航, 邓双厚, 支豪林, 朱震浩, 陆召严. 着水初始条件对水陆两栖飞机着水性能的影响[J]. 航空学报, 2021, 42(7): 124483-124483.
[8]	吴瀚, 王建宏, 黄伟, 杜兆波, 颜力. 激波/边界层干扰及微型涡流发生器控制研究进展[J]. 航空学报, 2021, 42(6): 25371-025371.
[9]	禹旻, 杨武兵, 沈清. 超声速尾迹-剪切流的混合增强[J]. 航空学报, 2021, 42(12): 625876-625876.
[10]	常思源, 白晓征, 崔小强, 刘君. 一种改进的非定常激波装配算法[J]. 航空学报, 2020, 41(2): 123498-123498.
[11]	李锦, 耿湘人, 陈坚强, 江定武, 李红喆. DSMC量子动理学化学反应模型的数值模拟[J]. 航空学报, 2019, 40(4): 122176-122176.
[12]	李佳伟, 王江峰, 杨天鹏, 李龙飞, 王丁. “Ⅳ型”激波干扰中流-热-固耦合问题一体化计算分析[J]. 航空学报, 2019, 40(12): 123190-123190.
[13]	金禹彤, 陈吉昌, 卢昱锦, 肖天航, 童明波. 楔形体入波浪水面数值模拟[J]. 航空学报, 2019, 40(10): 122854-122854.
[14]	江中正, 赵文文, 袁震宇, 陈伟芳. 基于非线性耦合本构关系的改进边界条件[J]. 航空学报, 2018, 39(10): 122057-122057.
[15]	刘君, 董海波, 刘瑜. 化学非平衡流动解耦算法的回顾与新进展[J]. 航空学报, 2018, 39(1): 21090-021090.