基于离散等价方程的非结构网格有限差分法

doi:10.7527/S1000-6893.2019.23248

流体力学与飞行力学

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

前一篇 | 后一篇

基于离散等价方程的非结构网格有限差分法

刘君, 陈洁, 韩芳

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

收稿日期:2019-06-26 修回日期:2019-08-06 出版日期:2020-01-15 发布日期:2019-09-02
通讯作者: 刘君 E-mail:liujun65@dlut.edu.cn
基金资助:
国家自然科学基金（11872144）

Finite difference method for unstructured grid based on discrete equivalent equation

LIU Jun, CHEN Jie, HAN Fang

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

Received:2019-06-26 Revised:2019-08-06 Online:2020-01-15 Published:2019-09-02
Supported by:
National Natural Science Foundation of China (11872144)

摘要/Abstract

摘要： 激波装配法把解析关系式嵌入流场避免了间断引起的理论问题，使全场一致高精度的实现成为可能，有望解决超声速流动转捩研究中的感受性模拟难题。但装配复杂激波结构以后，被分割的流场空间常出现不规则的几何形状，这给常规的基于结构网格的有限差分法应用带来困难。基于离散等价方程理论，提出了一种新的基于非结构网格的有限差分法，在空间二维离散点附近仅用3条网格线就可以构造出一阶迎风格式。数值算例表明收敛过程对网格质量不敏感，解决了激波装配法模拟激波相交出现小夹角后使用结构网格进行计算存在的难题。根据这种方法的特点展望了未来的应用前景。

关键词: 有限差分法, 非结构网格, 激波装配法, 离散等价方程, 流场

Abstract: The shock-fitting method embeds the analytical relationship into the flow field to avoid the theoretical problems caused by discontinuity, making it possible to achieve consistent high order scheme in the whole flow field. It is expected to solve the susceptibility simulation problem in the research of supersonic flow transition. However, after complex shock wave structures are fitted, the segmented flow field space often has irregular geometric shapes, bringing difficulties to the conventional finite difference method based on structural mesh. In this paper, based on the theory of discrete equivalence equations, a new finite difference method for unstructured mesh is proposed. In the case of two-dimensional space, a first-order upwind scheme can be constructed by only using three mesh lines at discrete points. Numerical examples show that the convergence process is not sensitive to the quality of the mesh. The proposed scheme is used in the shock-fitting simulation to solve the problem that the structural mesh appears after the shock wave intersects with a small angle. Finally, this paper forecasts the application prospects of the proposed scheme according to its characteristics.

Key words: finite difference method, unstructured mesh, shock-fitting method, discrete equivalence equation, flow field

中图分类号:

V211.3

刘君, 陈洁, 韩芳. 基于离散等价方程的非结构网格有限差分法[J]. 航空学报, 2020, 41(1): 123248-123248.

LIU Jun, CHEN Jie, HAN Fang. Finite difference method for unstructured grid based on discrete equivalent equation[J]. ACTA AERONAUTICAET ASTRONAUTICA SINICA, 2020, 41(1): 123248-123248.

参考文献 29

[1]	MORETTI G. Thirty-six years of shock fitting[J]. Computers & Fluids, 2002, 31:719-723.
[2]	HARTEN A. High resolution schemes for hyperbolic conservation laws[J]. Journal of Computational Physics, 1983, 49(3):357-393.
[3]	HARTEN A, ENGQUIST B, OSHER S, et al. Uniformly high order accurate essentially non-oscillatory schemes III[M]. Berlin Heidelber:Springer, 1987:231-303.
[4]	SHI J, HU C, SHU C W. A technique of treating negative weights in WENO schemes[J]. Journal of Computational Physics, 2002,175(1):108-127.
[5]	LEVY D, PUPPO G, RUSSO G. On the behavior of the total variation in CWENO methods for conservation laws[J]. Applied Numerical Mathematics, 2000, 33(1):407-414.
[6]	QIU J, SHU C W. On the construction, comparison, and local characteristic decomposition for high-order central WENO schemes[J]. Journal of Computational Physics, 2002, 183(1):187-209.
[7]	WANG Z J, CHEN R F. Optimized weighted essentially nonoscillatory schemes for linear waves with discontinuity[J]. Journal of Computational Physics, 2001, 174(1):381-404.
[8]	宗文刚, 邓小刚, 张涵信. 双重加权实质无波动激波捕捉格式[J]. 空气动力学学报, 2003, 21(2):218-225. ZONG W G, DENG X G, ZHANG H X. Double weighted essentially non-oscillatory shock-capturing schemes[J]. Acta Aerodynamica Sinica, 2003, 21(2):218-225(in Chinese).
[9]	沈孟育, 李海东, 刘秋生. 用解析离散法构造WENO-FCT格式[J]. 空气动力学学报, 1998, 16(1):56-63. SHEN M Y, LI H D, LIU Q S. Analytic discrete WENO-FCT scheme[J]. Acta Aerodynamica Sinica, 1998, 16(1):56-63(in Chinese).
[10]	DENG X G, MAEKAWA H. Compact high-order accurate nonlinear schemes[J]. Journal of Computational Physics, 1997, 130:77-91.
[11]	DENG X G, MAO M. Weighted compact high-order nonlinear schemes for the Euler equations[C]//Computational Fluid Dynamics Conference, 2006.
[12]	ZHONG X L. High-order finite-difference schemes for numerical simulation of hypersonic boundary-layer transition[J]. Journal of Computational Physics, 1998, 144(2):662-709.
[13]	PRADEEP S R, ZHONG X L. High-order shock-fitting and front-tracking schemes for numerical simulation of shock-disturbance interactions[J]. Journal of Computational Physics, 2010, 229(19):6744-6780.
[14]	ZHONG X L, WANG X W. Direct numerical simulation on the receptivity, instability, andtransition of hypersonic boundary layers[J]. Annual Review of Fluid Mechanics, 2012, 44(1):527-561.
[15]	周恒, 张涵信. 有关近空间高超声速飞行器边界层转捩和湍流的两个问题[J]. 空气动力学学报, 2017, 35(2):151-155. ZHOU H, ZHANG H X. Two problems in the transition and turbulence for near space hypersonic flying vehicles[J]. Acta Aerodynamica Sinica, 2017, 35(2):151-155(in Chinese).
[16]	MARCELLO O, RENATO P. Shock-fitting:Classical technique, recent developments, and memoirs of Gino Moretti[M]. Springer, 2017.
[17]	CHUNG K C. A generalized finite-difference method for heat transfer problems of irregular geometries[J]. Numerical Heat Transfer, 1981, 4:345-357.
[18]	BATINA J T. A gridless Euler/Navier-Stokes solution algorithm for complex-aircraft applications:AIAA-1993-0333[R]. Reston, VA:AIAA, 1993.
[19]	RAINALD L, CARLOS S, EUGENIO O, et al. A finite point method for compressible flow[J]. International Journal for Numerical Methods in Engineering, 2002, 53:1765-1779.
[20]	SRIDAR D, BALAKRISHNAN N. An upwind finite difference scheme for meshless solvers[J]. Journal of Computational Physics, 2003, 189:1-29.
[21]	DING H, SHU C, YEO K S. Development of least-square-based two-dimensional finite difference schemes and their application to simulate natural convection in a cavity[J]. Computers & Fluids, 2004, 33:137-154.
[22]	SHU C, DING H, CHEN H Q, et al. An upwind local RBF-DQ method for simulation of inviscid compressible flows[J]. Computer Methods in Applied Mechanics and Engineering, 2005, 194:2001-2017.
[23]	TOTA P, WANG Z J. Meshfree Euler solver using local radial basis functions for inviscid compressible flows:AIAA-2007-4581[R]. Reston, VA:AIAA, 2007.
[24]	KATZ A, JAMENSON A. A comparison of various meshless schemes within a unified algorithm:AIAA-2009-0596[R]. Reston, VA:AIAA, 2009.
[25]	SU X R, YAMAMOTO S, NAKAHASHI K. Analysis of a meshless solver for high Reynolds number flow[J]. Journal of Computational Physics, 2013, 72:505-527.
[26]	SUNDAR D S, YEO K S. A high order meshless method with compact support[J]. Journal of Computational Physics, 2014, 272:70-87.
[27]	XUE L L, YU X R, W L. Construction of the high order accurate generalized finite difference schemes for inviscid compressible flows[J]. Computers & Fluids, 2019, 25(2):481-507.
[28]	刘君, 韩芳, 夏冰. 有限差分法中几何守恒律的机理及算法[J]. 空气动力学学报, 2018, 36(6):918-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):918-926(in Chinese).
[29]	刘君, 韩芳. 有限差分法中的贴体坐标变换[J]. 气体物理, 2019(5):18-29. LIU J, HAN F. Body-fitted coordinate transformation for finite difference method[J]. Physics of Gases, 2019(5):18-29(in Chinese).

E-mail：hkxb@buaa.edu.cn

关于我们

期刊社服务

专业学科

封面文章

友情链接

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

基于离散等价方程的非结构网格有限差分法

Finite difference method for unstructured grid based on discrete equivalent equation

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献 29

相关文章 15

编辑推荐

Metrics

本文评价

[1]	吕治国, 龚红明, 常雨, 廖振洋, 钟涌. 跨越第二宇宙速度的膨胀管风洞研制[J]. 航空学报, 2022, 43(S2): 138-151.
[2]	程剑锐, 施崇广, 瞿丽霞, 徐悦, 尤延铖, 朱呈祥. 二维弯曲激波/湍流边界层干扰流动理论建模[J]. 航空学报, 2022, 43(9): 125993-125993.
[3]	伍强, 徐浩军, 魏扬, 裴彬彬, 薛源. 结冰条件下飞机气动/运动耦合特性[J]. 航空学报, 2022, 43(8): 125566-125566.
[4]	齐龙舟, 赵鲲, 冯和英, 张俊龙, 覃晨. 不同冲击距离下斜板开槽对超声速射流冲击噪声的影响[J]. 航空学报, 2022, 43(8): 125712-125712.
[5]	陈浩, 毕林, 华如豪, 周清清, 唐志共, 袁先旭. 基于全线程树数据结构的笛卡尔网格高效生成技术[J]. 航空学报, 2022, 43(5): 125170-125170.
[6]	王国乾, 李寒松, 韩国峰, 李若琪, 张彦. 基于周期性换向冲液的掩膜电解加工技术[J]. 航空学报, 2022, 43(4): 525575-525575.
[7]	刘君, 韩芳, 魏雁昕. 特定条件下高阶WENO格式计算结果误差[J]. 航空学报, 2022, 43(2): 124940-124940.
[8]	杨鲤铭, 吴杰, 董昊, 杜银杰. 基于速度空间非结构网格和守恒修正的改进离散速度方法[J]. 航空学报, 2022, 43(12): 627033-627033.
[9]	刘佳豪, 李高华, 王福新. 倾转过渡状态旋翼-机翼气动干扰特性[J]. 航空学报, 2022, 43(12): 126097-126097.
[10]	战庆亮, 白春锦, 张宁, 葛耀君. 基于时程卷积自编码的机翼绕流特征识别方法[J]. 航空学报, 2022, 43(11): 526531-526531.
[11]	陈呈, 赵丹, 王岳青, 邓亮, 杨超, 苏铖宇, 王昉. NNW-TopViz流场可视分析系统[J]. 航空学报, 2021, 42(9): 625747-625747.
[12]	杨昌和, 李彦达, 张江, 王昉, 袁晓如. 高性能流场并行粒子追踪数据管理系统[J]. 航空学报, 2021, 42(9): 625748-625748.
[13]	陈坚强, 吴晓军, 张健, 李彬, 贾洪印, 周乃春. FlowStar:国家数值风洞(NNW)工程非结构通用CFD软件[J]. 航空学报, 2021, 42(9): 625739-625739.
[14]	任玉新, 王乾, 潘建华, 章雨思, 黄乾旻. 构建非结构网格高精度有限体积方法的新途径[J]. 航空学报, 2021, 42(9): 625783-625783.
[15]	周印佳, 张志贤, 付新卫, 阿嵘. 再入飞行器烧蚀热防护一体化计算方法[J]. 航空学报, 2021, 42(7): 124520-124520.