Discussions on two problems in applications of high-order finite difference schemes

LIU Jun; HAN Fang

doi:10.7638/kqdlxxb-2019.0073

English Version Journal Login

您当前的位置：

首页 >

文章列表页 >

Discussions on two problems in applications of high-order finite difference schemes

更新时间：2026-05-28

- Discussions on two problems in applications of high-order finite difference schemes
- Vol. 38, Issue 2, Pages: 244-253(2020)
- 作者机构：
  
  大连理工大学航空航天学院,大连,116024
- 作者简介：
- 基金信息：
- DOI：10.7638/kqdlxxb-2019.0073
  CLC： V211.3
- Published：2020
- 稿件说明：
移动端阅览
刘君, 韩芳. 有关有限差分高精度格式两个应用问题的讨论[J]. 空气动力学学报, 2020, 38(2): 244-253.

LIU Jun, HAN Fang. Discussions on two problems in applications of high-order finite difference schemes[J]. 2020, 38(2): 244-253.
刘君, 韩芳. 有关有限差分高精度格式两个应用问题的讨论[J]. 空气动力学学报, 2020, 38(2): 244-253. DOI： 10.7638/kqdlxxb-2019.0073.

LIU Jun, HAN Fang. Discussions on two problems in applications of high-order finite difference schemes[J]. 2020, 38(2): 244-253. DOI： 10.7638/kqdlxxb-2019.0073.

摘要

激波装配法结合有限差分高精度格式是发展全场一致高精度算法的一种途径。在对流场内的间断全部进行装配后，对高精度格式的应用研究成为发展本算法的主要研究问题。本文将有限差分的高精度格式应用于贴体坐标系时发现，对均匀流场，高精度格式因不满足几何守恒律而产生的数值误差比一阶迎风格式大，初步分析认为是由于高精度格式所用的模板比一阶格式更宽，涉及的网格点数更多，从而引入了更多的误差。而作者提出的基于离散等价方程的相容性算法可消除这一误差。此外，本文在利用激波捕捉法求解正方形均匀网格上的正激波运动问题时发现因通量分裂格式的使用，在激波处会产生随着特征线传播的非物理波动，这一波动在激波与网格不完全匹配时表现为多维波动相互干扰的虚假"数值湍流"现象，高精度格式的高分辨率特性使得这一现象更加明显。这是因为激波捕捉法假设激波为空间连续函数，用于包含激波的流场时必然得到数值解表示的过渡区，导致可信度评估困难，使用激波装配法可以避免这一问题。

Abstract

One way to achieve high degree fidelity in the whole flow field is by developing shock-fitting methods combined with high-order finite difference schemes. After all the discontinuities in shock-embedded compressible fluid flow are assembled

the application of the high-order finite difference schemes becomes the main problem which we need to consider. It has been found that the finite difference schemes always make numerical errors under curvilinear coordinate system due to the unfulfilled geometric conservation law. Moreover

we found that the errors of high-order schemes are greater than that of first-order upwind scheme for uniform flow. The reason given in this paper is that the stencils of high-order schemes involve more mesh points than that of first-order scheme. The source-coupled algorithm based on discrete equivalence equations proposed by the authors can eliminate the errors. Besides

the problem of normal shock's motion on uniform grid is also discussed in this paper. Because of the flux splitting

two waves can be created after the moving normal shock. These two waves travel at two different characteristic velocities

and can lead to numerical turbulence when the shock does not exactly match to grid points. Since the shock-capture methods assume that shock is spatial continuous

transition regions represented by the numerical solution leads to difficulties in the credibility evaluation. Shock-fitting methods can avoid this problem.

关键词

Keywords

references

CHANG S Y, ZOU D Y, LIU J. Simulating hypersonic projectile launching process in the ballistic range by adaptive discontinuity fitting solver technique[J]. Journal of Experiments in Fluid Mechanics, 2019, 33(2):23-29. (in Chinese)

LIU J, LIU Y, CHEN Z D. Unstructured deforming mesh and discrete geometric conservation law[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(8):2395-2407. (in Chinese)

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.

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)

ZOU D Y, LIU J, ZOU L. Applications of shock-fitting technique for compressible flow in cell-centered finite volume methods[J]. Acta Aeronautica et Astronautica Sinica, 2017, 38(11):139-149. (in Chinese)

MARCELLO O, RENATO P. Shock fitting: classical techniques, recent developments, and memoirs of gino moretti[M]. Springer, 2017: 131-149.

ZHANG L P, WANG Z J. A block LU-SGS implicit dual time-stepping algorithm for hybrid dynamic meshes[J]. Computers Fluids, 2004, 33(7):891-916.

LI W N, REN Y X. The multi-dimensional limiters for solving hyperbolic conservation laws on unstructured grids Ⅱ:Extension to high order finite volume schemes[J]. Journal of Computational Physics, 2012, 231(11):4053-4077.

THOMAS P D, LOMBARD C K. The geometric conservation law-a link between finite-difference and finite-volume methods of flow computation on moving grids[C]//11th Fluid and PlasmaDynamics Conference, Seattle, WA, USA. Reston, Virigina: AIAA, 1978.

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.

ABE Y, ⅡZUKA N, NONOMURA T, et al. Symmetric-conservative metric evaluations for higher-order finite difference scheme with the GCL identities on three-dimensional moving and deforming meshs[J]. Seventh International Conference on Computational Fluid Dynamics (ICCFD7), Big Island, Hawaii, 2012.

ABE Y, ⅡZUKA N, NONOMURA T, et al. Conservative metric evaluation for high-order finite difference schemes with the GCL identities on moving and deforming grids[J]. Journal of Computational Physics, 2013, 232(1):14-21.

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.

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.

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.

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.

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)

JIANG G S, SHU C W. Efficient implementation of weighted ENO schemes[J]. Journal of Computational Physics, 1996, 126(1):202-228.

HENRICK A K, ASLAM T D, POWERS J M. Mapped weighted essentially non-oscillatory schemes:Achieving optimal order near critical points[J]. Journal of Computational Physics, 2005, 207(2):542-567.

BORGES R, CARMONA M, COSTA B, et al. An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws[J]. Journal of Computational Physics, 2008, 227(6):3191-3211.

FU L, HU X Y, ADAMS N A. A family of high-order targeted ENO schemes for compressible-fluid simulations[J]. Journal of Computational Physics, 2016, 305:333-359.

CAI X D, LADEINDE F. Performance of WENO scheme in generalized curvilinear coordinate systems[C]//46th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada. Reston, Virigina: AIAA, 2008: 36.

LIU J, HAN F. Body-fitted coordinate transformation for finite difference method[J]. Physics of Gases, 2018, 3(5):18-29. (in Chinese)

Views

下载量

CSCD

Alert me when the article has been cited

提交

Tools

Publicity Resources

Initial discontinuity induced non-physical fluctuation and its elimination algorithm

A constrained slip-wall boundary condition based on discrete equivalent equations

Positivity-preserving strategies and their comparision based on staggered CPR methods

Multi-factor analysis and simulation of particle erosion characteristics of fan blades

Numerical simulation on gust responses of high aspect ratio aircraft

Related Author

LIU Jun

HAN Fang

SUN Yixuan

LIU Yu

WEI Yanxin

LIU Jun

Related Institution

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

北京宇航系统工程研究所

宁波大学机械工程与力学学院

AI问答

⁰