An accurate conservative interpolation method for the mixed gridbased on the intersection of grid cells
-
摘要: 基于网格切割思想,发展了二维/三维混合网格条件下的单元相交算法,可精确计算任意两个多边形/多面体的交集。在此基础上,实现了基于单元相交(CIB/DC)的精确守恒插值算法。二维和三维验证算例表明,该方法能够保证插值过程中计算域内物理量的严格守恒,且具有比常规二阶插值更高的精度。Abstract: A method is presented for conservatively transferring cell-centered physical quantities from one mesh to another with second-order accuracy, which is well-suited for finite-volume methods that rely on cell-centered variables. The proposed methodology implements the cell-intersection algorithm of 2D/3D mixed grid based on the local supermesh idea, and is able to accurately calculate the intersection of any two polygons or polyhedrons, thus providing a basis for accurate conservative interpolation. The accuracy and the efficacy of the new method are demonstrated with multiple 2D and 3D numerical experiments. The test results show that this method ensures strict conservation of physical quantities in the interpolation process, and achieves a higher accuracy than that achieved by conventional second-order interpolation methods.
-
Key words:
- fluid mechanics /
- conservative interpolation /
- cell-intersection /
- mixed grid /
- local supermesh
-
表 1 计算网格参数
Table 1. Parameters of the computational grid
序号 结构网格 非结构网格 节点数 单元数 1 33×33 695 1 260 2 65×65 2 593 4 928 3 129×129 10 044 19 574 4 257×257 39 687 78 348 表 2 不同算例的守恒误差
Table 2. Conservation errors indifferent numerical cases
序号 二阶插值 守恒型二阶插值 1 3.724×10-3 0 2 -1.447×10-3 0 3 -3.985×10-4 1×10-14 4 5.636×10-5 -2×10-14 表 3 计算网格参数
Table 3. Parameters of the computational grid
序号 网格1 网格2 dx 节点数 单元数 节点数 单元数 1 914 4 033 676 2 711 0.067 2 2 222 10 632 1 424 6 075 0.050 3 6 527 32 913 3 733 16 681 0.035 4 26 099 138 542 12 600 59 347 0.022 -
[1] 郭正, 刘君, 瞿章华.非结构动网格在三维可动边界问题中的应用[J].力学学报, 2003, 35(2):140-146. http://d.old.wanfangdata.com.cn/Periodical/lxxb200302003Guo Zheng, Liu Jun, Qu Zhanghua. Dynamic unstructured grid method with applications to 3d unsteady flows involving moving boundaries[J]. Acta Mechanica Sinica, 2003, 35(2):140-146. http://d.old.wanfangdata.com.cn/Periodical/lxxb200302003 [2] Zeeuw D D, Powell K G. An adaptively refined Cartesian mesh solver for the Euler equations[J]. Journal of Computational Physics, 1993, 104(1), 104:56-68. http://cn.bing.com/academic/profile?id=0b27e7efb53cada94a6d8ccfd529a00b&encoded=0&v=paper_preview&mkt=zh-cn [3] 白晓征.包含运动界面的爆炸流场数值模拟方法及其应用[D].长沙: 国防科技大学, 2009. http://cdmd.cnki.com.cn/Article/CDMD-90002-2010147319.htm [4] 徐春光, 白晓征, 刘瑜, 等.爆炸近区流场的数值模拟方法研究[J].兵工学报, 2012, 33(5):565-573. http://d.old.wanfangdata.com.cn/Periodical/bgxb201205010Xu Chunguang, Bai Xiaozheng, Liu Yu, et al. Research on numerical simulation method of near-field flows in air blast[J]. Acta Armmamentarii, 2012, 33(5):565-573. http://d.old.wanfangdata.com.cn/Periodical/bgxb201205010 [5] Pärt-Enander E, Sj green B. Conservative and non-conservative interpolation between overlapping grids for finite volume solutions of hyperbolic problems[J]. Computers and Fluids, 1994, 23(3):551-574. http://cn.bing.com/academic/profile?id=1b04887df045469780080f87398e4aca&encoded=0&v=paper_preview&mkt=zh-cn [6] Margolin L G, Shashkov M. Second-order sign-preserving conservative interpolation (remapping) on general grids[J]. Journal of Computational Physics, 2002, 184(1):266-298. [7] Farrell P E, Piggott M D, Pain C C, et al. Conservative interpolation between unstructured meshes via supermesh construction[J]. Computer Methods in Applied Mechanics and Engineering, 2009, 198(33/34/35/36):2632-2642. http://d.old.wanfangdata.com.cn/NSTLQK/NSTL_QKJJ0211806853/ [8] Farrell P E, Maddison J R. Conservative interpolation between volume meshes by local Galerkin projection[J]. Computer Methods in Applied Mechanics and Engineering, 2011, 200(1/2/3/4):89-100. http://cn.bing.com/academic/profile?id=aac95ff70f0ddb440fb234e1c139f5d9&encoded=0&v=paper_preview&mkt=zh-cn [9] Menon S, Schmidt D P. Conservative interpolation on unstructured polyhedral meshes: An extension of the supermesh approach to cell-centered finite-volume variables[J]. Computer Methods in Applied Mechanics and Engineering, 2011, 200(41/44):2797-2804. [10] 郭正.包含运动边界的多体非定常流场数值模拟方法研究[D].长沙: 国防科技大学, 2002. http://cdmd.cnki.com.cn/Article/CDMD-90002-2003097667.htm [11] 刘君, 白晓征, 郭正.非结构动网格计算方法及其在包含运动界面的流场模拟中的应用[M].长沙:国防科技大学出版社, 2009. 期刊类型引用(18)
1. 高飞,邓树新,张国凯,纪玉国,刘晨康,王明洋. 缩比模型弹侵彻岩石靶尺寸效应试验研究与理论分析. 兵工学报. 2023(12): 3601-3612 . 百度学术
2. 朱少平,王志亮,熊峰. 卵形弹丸对混凝土侵彻动力响应数值研究. 合肥工业大学学报(自然科学版). 2022(02): 243-250 . 百度学术
3. 邵伟,范锦彪,耿宇飞,王玮. 基于微元法的侵彻体弹头摩擦升温计算方法. 探测与控制学报. 2022(02): 34-40 . 百度学术
4. 章毅,胡枫,吴昊,王安宝. 高强混凝土介质侵彻系数计算方法. 防护工程. 2021(02): 31-38 . 百度学术
5. 张丁山,谷鸿平,徐笑,张博,吕永柱. 截卵形头部平台直径对初始侵彻弹道偏转的影响. 高压物理学报. 2021(05): 138-144 . 百度学术
6. 吴飚,任辉启,陈力,杨建超,黄家蓉,高伟亮,金栋梁. 弹体侵彻混凝土尺度效应试验研究与理论分析. 防护工程. 2020(02): 1-10 . 百度学术
7. 彭永,卢芳云,方秦,吴昊,李翔宇. 弹体侵彻混凝土靶体的尺寸效应分析. 爆炸与冲击. 2019(11): 58-68 . 本站查看
8. 张学伦,汪衡,谭正军,王昭明. 混凝土靶边界效应与弹丸长径比关联性的研究. 兵器装备工程学报. 2018(04): 11-13+18 . 百度学术
9. 刘宗伟,武海军,张学伦,刘俞平,熊国松,谭正军,曾令清. 高超弹体侵蚀机理及抗侵蚀设计研究. 兵器装备工程学报. 2017(04): 46-49 . 百度学术
10. 李艳,范文,赵均海,翟越. 中低速长杆弹侵彻半无限岩石靶的动态响应研究. 工程力学. 2017(09): 139-149 . 百度学术
11. 薛建锋,沈培辉,王晓鸣. 基于层裂机理的弹体侵彻混凝土的工程模型. 国防科技大学学报. 2017(03): 194-200 . 百度学术
12. 曹扬悦也,蒋志刚,谭清华,蒙朝美. 基于Hoek-Brown准则的混凝土-岩石类靶侵彻模型. 振动与冲击. 2017(05): 48-53+60 . 百度学术
13. 薛建锋,沈培辉,王晓鸣. 不同头部形状弹体侵彻混凝土的试验研究. 兵工自动化. 2016(02): 75-78 . 百度学术
14. 邓佳杰,张先锋,乔治军,郭磊,何勇,陈东东. 卵形弹体侵彻预开孔靶理论分析. 爆炸与冲击. 2016(05): 625-632 . 本站查看
15. 张学伦,刘宗伟. 弹丸CRH值对侵彻混凝土深度影响研究. 兵器装备工程学报. 2016(10): 31-34 . 百度学术
16. 薛建锋,沈培辉,王晓鸣. 钻地弹斜侵彻混凝土靶的工程计算模型. 航空学报. 2016(06): 1899-1911 . 百度学术
17. 张丁山,吕永柱,周涛,谷鸿平,张立建. 侵彻战斗部引信前后置过载的影响因素. 探测与控制学报. 2016(06): 41-45+50 . 百度学术
18. 彭永,方秦,吴昊,孔祥振,肖云凯. 不同头部形状弹体侵彻混凝土靶体的终点弹道参数分析. 兵工学报. 2014(S2): 128-134 . 百度学术
其他类型引用(11)
-