有限差分法的坐标变换诱导误差

doi:10.7527/S1000-6893.2020.24397

流体力学与飞行力学

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

前一篇 | 后一篇

有限差分法的坐标变换诱导误差

刘君, 魏雁昕, 韩芳

大连理工大学航空航天学院, 大连 116024

收稿日期:2020-06-11 修回日期:2020-09-18 出版日期:2021-06-15 发布日期:1900-01-01
通讯作者: 刘君 E-mail:liujun65@dlut.edu.cn
基金资助:
国家重点研发计划（2018YFB0204404）；国家自然科学基金（11872144）

Coordinate transformation induced errors of finite difference method

LIU Jun, WEI Yanxin, HAN Fang

School of Aeronautics and Astronautics, Dalian University of Technology, Dalian 116024, China

Received:2020-06-11 Revised:2020-09-18 Online:2021-06-15 Published:1900-01-01
Supported by:
National Key Research & Development Program of China (2018YFB0204404); National Natural Science Foundation of China (11872144)

摘要/Abstract

摘要： 有限差分法应用于具有复杂外形的网格时需要进行坐标变换，在此过程中经常会引入坐标变换诱导误差。在柱坐标系下使用均匀网格进行均匀流场计算，计算结果表明，即使物理坐标对计算坐标的变换函数连续可导、计算过程中坐标变换系数直接采用准确的解析式、网格完全正交并且充分光滑，也无法避免坐标变换诱导误差。理论分析表明，产生坐标变换诱导误差的机理是笛卡尔坐标系下的守恒型欧拉方程变换至贴体坐标系下后增加了源项。针对该问题，目前国内外学者通常采用几何守恒律，构建与差分格式相匹配的坐标变换系数计算方法来消除源项。本文介绍了从包含源项的离散等价方程基础上直接进行离散的新算法，在此基础上针对非等距网格条件下MUSCL类格式重构过程进行误差分析，理论推导表明重构中需要考虑非等距插值公式的影响系数，将变量转换至计算空间内进行MUSCL重构才能保证该过程具有均匀网格下的插值精度。通过理论分析及数值实验证明新算法对于均匀流场完全不会引入坐标变换误差。

关键词: 有限差分法, 几何守恒律, 离散等价方程, 几何诱导误差, 坐标变换

Abstract: Coordinate transformation is required when the finite difference method is applied to the mesh with complex geometries, and the errors induced by the coordinate transformation are often introduced in this process. These errors are proved to be inevitable in the uniform flow field calculation with uniform grids in cylindrical coordinate systems, even if the transformation function of the physical coordinates to the calculated coordinates is continuously derivable, or the coordinate transformation coefficients in the calculation process are calculated by the accurate analytical formula, or the grid is completely orthogonal and fully smooth. Theoretical analysis shows the mechanism of the coordinate transformation induced errors:when the conservative Euler equation is transformed from the Cartesian coordinate system to the body fitted coordinate system, a source term is added. Currently, scholars usually use the geometric conservation law to construct a method based on coordinate transformation coefficients, which are matched with the format of the finite difference, to eliminate the source term. In this work, we introduce a new algorithm that processes the direct discretization from the discrete equivalent functions including the source term. Based on the above new algorithm, error analysis is carried out for the reconstruction process of MUSCL format under non-equidistant grid conditions. Theoretical derivation shows that the influence coefficient of the non-equidistant interpolation formula needs to be considered in reconstruction, only when the variables are transformed into the computational space for MUSCL reconstruction can the interpolation accuracy be guaranteed under uniform grid. Our theoretical analysis and numerical experiments have proven that this algorithm will not introduce coordinate transformation errors to the uniform flow field calculations.

Key words: finite difference method, geometric conservation law, discrete equivalence equation, geometrically induced errors, coordinate transformation

中图分类号:

V211.3

刘君, 魏雁昕, 韩芳. 有限差分法的坐标变换诱导误差[J]. 航空学报, 2021, 42(6): 124397-124397.

LIU Jun, WEI Yanxin, HAN Fang. Coordinate transformation induced errors of finite difference method[J]. ACTA AERONAUTICAET ASTRONAUTICA SINICA, 2021, 42(6): 124397-124397.

参考文献

[1] TRULIO J G, TRIGGER K R. Numerical solution of the one-dimensional Lagrangian hydrodynamic equations:UCRL-6267[R]. Oak Ridge:Office of Scientific and Technical Information (OSTI), 1961.
[2] STEGER J. Implicit finite difference simulation of flow about arbitrary geometries with application to airfoils[C]//10th Fluid and Plasmadynamics Conference. Reston:AIAA, 1977:665.
[3] STEGER J L. Implicit finite-difference simulation of flow about arbitrary two-dimensional geometries[J]. AIAA Journal, 1978, 16(7):679-686.
[4] PULLIAM T, STEGER J. On implicit finite-difference simulations of three-dimensional flow[C]//16th Aerospace Sciences Meeting. Reston:AIAA, 1978:2514.
[5] HINDMAN R. Geometrically induced errors and their relationship to the form of the governing equations and the treatment of generalized mappings[C]//5th Computational Fluid Dynamics Conference. Reston:AIAA, 1981:81-1008.
[6] HINDMAN R G. Generalized coordinate forms of governing fluid equations and associated geometrically induced errors[J]. AIAA Journal, 1982, 20(10):1359-1367.
[7] VIVIAND H, GHAZZI W. Numerical solution of the compressible Navier-Stokes equations at high Reynolds numbers with applications to the blunt body problem[C]//Proceedings of the 5th International Conference on Numerical Methods in Fluid Dynamics. Enschede:Twente University, 1976:434-439.
[8] THOMAS P, LOMBARD C. The Geometric Conservation Law-A link between finite-difference and finite-volume methods of flow computation on moving grids[C]//11th Fluid and Plasma Dynamics Conference. Reston:AIAA, 1978:1208.
[9] THOMAS P D, LOMBARD C K. Geometric conservation law and its application to flow computations on moving grids[J]. AIAA Journal, 1979, 17(10):1030-1037.
[10] GAITONDE D, VISBAL M. Further development of a Navier-Stokes solution procedure based on higher-order formulas[C]//37th Aerospace Sciences Meeting and Exhibit. Reston:AIAA, 1999:557.
[11] DENG X G, MAO M L, TU G H, et al. Geometric conservation law and applications to high-order finite difference schemes with stationary grids[J]. Journal of Computational Physics, 2011, 230(4):1100-1115.
[12] DENG X G, MIN Y B, MAO M L, et al. Further studies on Geometric Conservation Law and applications to high-order finite difference schemes with stationary grids[J]. Journal of Computational Physics, 2013, 239:90-111.
[13] NONOMURA T, TERAKADO D, ABE Y, et al. A new technique for freestream preservation of finite-difference WENO on curvilinear grid[J]. Computers & Fluids, 2015, 107:242-255.
[14] ZHU Y J, SUN Z S, REN Y X, et al. A numerical strategy for freestream preservation of the high order weighted essentially non-oscillatory schemes on stationary curvilinear grids[J]. Journal of Scientific Computing, 2017, 72(3):1021-1048.
[15] 朱志斌, 杨武兵, 禹旻. 满足几何守恒律的WENO格式及其应用[J]. 计算力学学报, 2017, 34(6):779-784. ZHU Z B, YANG W B, YU M. A WENO scheme with geometric conservation law and its application[J]. Chinese Journal of Computational Mechanics, 2017, 34(6):779-784(in Chinese).
[16] 闵耀兵. 高阶精度有限差分方法几何守恒律研究[D]. 绵阳:中国空气动力研究与发展中心, 2015. MIN Y B. The studies on geometric conservation law for high order finite difference method[D]. Mianyang:China Aerodynamics Research and Development Center, 2015(in Chinese).
[17] 张来平. 计算流体力学网格生成技术[M]. 北京:科学出版社, 2017:4. ZHANG L P. Mesh generation techniques in computational fluid dynamics[M]. Beijing:Science Press, 2017:4(in Chinese).
[18] 刘君, 韩芳, 夏冰. 有限差分法中几何守恒律的机理及算法[J]. 空气动力学学报, 2018, 36(6):917-926. LIU J, HAN F, XIA B. Mechanism and algorithm for geometric conservation law in finite difference method[J]. Acta Aerodynamica Sinica, 2018, 36(6):917-926(in Chinese).
[19] 刘君, 韩芳. 有限差分法中的贴体坐标变换[J]. 气体物理, 2018, 3(5):18-29. LIU J, HAN F. Body-fitted coordinate transformation for finite difference method[J]. Physics of Gases, 2018, 3(5):18-29(in Chinese).
[20] 刘君, 韩芳. 有关有限差分高精度格式两个应用问题的讨论[J]. 空气动力学学报, 2020, 38(2):244-253. LIU J, HAN F. Discussions on two problems in applications of high-order finite difference schemes[J]. Acta Aerodynamica Sinica, 2020, 38(2):244-253(in Chinese).
[21] VAN LEER B. Towards the ultimate conservative difference scheme:V. A second-order sequel to Godunov's method[J]. Journal of Computational Physics, 1979, 32(1):101-136.

E-mail：hkxb@buaa.edu.cn

关于我们

期刊社服务

专业学科

封面文章

友情链接

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

有限差分法的坐标变换诱导误差

Coordinate transformation induced errors of finite difference method

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

本文评价

[1]	薄钧天, 王国宏, 于洪波, 彭志刚. 基于平行线坐标变换的检测前跟踪算法[J]. 航空学报, 2022, 43(8): 325644-325644.
[2]	刘君, 韩芳, 魏雁昕. 应用维数分裂方法推广MUSCL和WENO格式的若干问题[J]. 航空学报, 2022, 43(3): 125009-125009.
[3]	刘君, 韩芳, 魏雁昕. 特定条件下高阶WENO格式计算结果误差[J]. 航空学报, 2022, 43(2): 124940-124940.
[4]	周印佳, 张志贤, 付新卫, 阿嵘. 再入飞行器烧蚀热防护一体化计算方法[J]. 航空学报, 2021, 42(7): 124520-124520.
[5]	刘君, 魏雁昕, 陈洁. 基于非结构网格有限差分法的扎染算法[J]. 航空学报, 2021, 42(7): 124557-124557.
[6]	张云光, 李志强, 王耀奇, 李红, 李淑慧. 2060铝锂合金冷模热成形界面换热系数确定的实验与算法[J]. 航空学报, 2021, 42(2): 423805-423805.
[7]	李超群, 李易, 张晨曦, 唐硕. 沿流向微结构沟槽流场直接数值模拟[J]. 航空学报, 2020, 41(11): 123628-123628.
[8]	刘君, 陈洁, 韩芳. 基于离散等价方程的非结构网格有限差分法[J]. 航空学报, 2020, 41(1): 123248-123248.
[9]	陈俊霖, 徐浩军, 魏小龙, 陈增辉, 吕晗阳. 感性耦合夹层等离子体隐身天线罩电磁散射分析[J]. 航空学报, 2018, 39(3): 321472-321472.
[10]	刘君, 刘瑜, 陈泽栋. 非结构变形网格和离散几何守恒律[J]. 航空学报, 2016, 37(8): 2395-2407.
[11]	齐振超, 张开富, 李原, 程晖. 大型壁板自动钻铆定位误差分析与优化[J]. 航空学报, 2015, 36(10): 3439-3449.
[12]	刘君, 徐春光, 张帆. 电磁阀颤振的流固耦合模拟研究[J]. 航空学报, 2014, 35(7): 1922-1930.
[13]	袁坚锋, 尼早, 陈保兴. 面内弯曲载荷作用下两边简支两边固支复合材料层合板的屈曲[J]. 航空学报, 2014, 35(4): 1026-1033.
[14]	魏东;周正干;. 固体中脉冲超声波传播的有限差分模拟[J]. 航空学报, 2010, 31(2): 388-392.
[15]	任元;孙玉坤;朱纪洪. 四相永磁容错电机的SVPWM控制[J]. 航空学报, 2009, 30(8): 1490-1496.