Initial discontinuity induced non-physical fluctuation and its elimination algorithm

LIU Jun; HAN Fang; SUN Yixuan; LIU Yu

doi:10.7638/kqdlxxb-2023.0203

English Version Journal Login

您当前的位置：

首页 >

文章列表页 >

Initial discontinuity induced non-physical fluctuation and its elimination algorithm

更新时间：2026-07-27

- Initial discontinuity induced non-physical fluctuation and its elimination algorithm
- Vol. 42, Issue 11, Pages: 18-32(2024)
- 作者机构：
  
  1. 大连理工大学航天航空学院,大连,116024
  2. 北京宇航系统工程研究所,北京,100076
  3. 宁波大学机械工程与力学学院,宁波,315211
- 作者简介：
- 基金信息：
- DOI：10.7638/kqdlxxb-2023.0203
  CLC： V211.3;O354.5
- Published：2024
- 稿件说明：
移动端阅览
刘君, 韩芳, 孙逸轩, 等. 初始间断诱导非物理波动及其消减算法[J]. 空气动力学学报, 2024, 42(11): 18−32.

LIU Jun, HAN Fang, SUN Yixuan, et al. Initial discontinuity induced non-physical fluctuation and its elimination algorithm[J]. 2024, 42(11): 18-32.
刘君, 韩芳, 孙逸轩, 等. 初始间断诱导非物理波动及其消减算法[J]. 空气动力学学报, 2024, 42(11): 18−32. DOI： 10.7638/kqdlxxb-2023.0203

LIU Jun, HAN Fang, SUN Yixuan, et al. Initial discontinuity induced non-physical fluctuation and its elimination algorithm[J]. 2024, 42(11): 18-32. DOI： 10.7638/kqdlxxb-2023.0203.

摘要

使用激波捕捉法求解以数学间断作为初始条件的无黏可压缩流场时，间断结构会逐渐演变为包含多个网格节点的数值过渡区，在此过程中会产生两个平行于间断且规则分布的非物理波动。使用捕捉法计算激波流场时，流场参数应满足修正方程，但作为初始条件进行计算的初始激波满足Euler方程，两者之间的矛盾是产生初始激波诱导非物理波动的原因。鉴于激波等间断是由特征型Euler方程定义的，从该方程出发构造了一种基于特征值迎风特性的通量计算格式（upwind flux scheme based on characteristics

UFSC），并采用若干种常规守恒型通量分裂格式作为参考。数值计算结果表明，UFSC可以消除Steger-Warming、Van Leer等矢通量分裂格式的初始接触间断诱导非物理波动，还可以减小初始激波诱导扰动的幅值。在光滑流场区域，UFSC和矢通量分裂格式的计算结果相近，但是在强激波附近会出现较大的压力尖峰。为克服这一缺陷，进一步构造了在激波区域和其他区域分别采用Steger-Warming格式和UFSC格式计算的混合通量计算格式UFSC+S，可以有效抑制初始激波诱导误差，提高计算精度。

Abstract

When the shock capturing method is used to solve inviscid compressible flow with a mathematical discontinuity as the initial condition

the discontinuity will gradually evolve into a numerical transition region across several grid points

and non-physical fluctuations in two directions parallel to the discontinuity can be induced in this process. The initial shock

which is calculated as the initial condition

satisfies the Euler equation

while the flow field parameters should satisfy the modified equations

such conflict might be the reason for the initial shock-induced non-physical fluctuations. As the discontinuities such as shock waves are defined by the characteristic lines of the Euler equation

an upwind flux scheme based on characteristics (UFSC) is proposed in the present study. Using several conventional conservation flux schemes as a reference

the UFSC can eliminate the non-physical fluctuations induced by the Steger-Warming and Van Leer schemes at the initial contact discontinuities

and can reduce the amplitude of the initial shock-induced disturbances. For the flow field in smooth regions

the computational results from the UFSC and conserved flux scheme are similar

but larger pressure peaks appear near the strong shock. To overcome this defect

a hybrid scheme UFSC+S is constructed

which uses the Steger-Warming for calculation in the shock region and the UFSC in other regions.

关键词

Keywords

references

LIU J, BAI X Z, ZHANG H X, et al . Discussion about GCL for deforming grids[J ] . Aeronautical Computing Technique, 2009 , 39 ( 4 ): 1 − 5 (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, CHEN Z D, CHANG S Y, et al . A shock detection method based on the theory of characteristics and its application in shock-fitting methods[J ] . Physics of Gases, 2018 , 3 ( 1 ): 41 − 49 (in Chinese).

LIU J, CHEN Z D, CHEN J, et al . A surface extraction method of flow feature for shockGfitting scheme[J ] . Acta Aerodynamica Sinica, 2019 , 37 ( 3 ): 373 − 381 (in Chine se).

CHANG S Y, BAI X Z, LIU J . A two-dimensional shock wave pattern recognition algorithm based on cluster analysis[J ] . Acta Aeronautica et Astronautica Sinica, 2020 , 41 ( 8 ): 123626 (in Chinese).

LIU J, ZOU D Y, XU C G . An unsteady shock-fitting technique based on unstructured moving grids[J ] . Acta Aerodynamica Sinica, 2015 , 33 ( 1 ): 10 − 16 (in Chinese).

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).

LIU J, ZOU D Y, DONG H B . Principle of new discontinuity fitting technique based on unstructured moving grid[J ] . Physics of Gases, 2017 , 2 ( 1 ): 13 − 20 (in Chinese).

ZOU D Y, XU C G, DONG H B, et al . A shock-fitting technique for cell-centered finite volume meth ods on unstructured dynamic meshes[J ] . Journal of Computational Physics, 2017 , 345 ( C ): 866 − 882 .

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 ): 121363 (in Chinese).

CHANG S, BAI X, ZOU D, et al . An adaptive discontinuity fitting technique on unstructured dynamic grids[J ] . Shock Waves, 2019 , 29 ( 8 ): 1103 − 1115 .

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).

JOHNSEN E . On the treatment of contact discontinuities using WENO schemes[J ] . Journal of Computational Physics, 2011 , 230 ( 24 ): 8665 − 8668 .

CLERC S . Accurate computation of contact discontinuities in flows with general equations of state[J ] . Computer Methods in Applied Mechanics and Engineering, 1999 , 178 ( 3-4 ): 245 − 255 .

NONOMURA T, MORIZAWA S, TERASHIMA H, et al . Numerical (error) issues on compressible multicomponent flows using a high-order differencing scheme: Weighted compact nonlinear scheme[J ] . Journal of Computational Physics, 2012 , 231 ( 8 ): 3181 − 3210 .

HE Z W, ZHANG Y S, LI X L, et al . Preventing numerical oscillations in the flux-split based finite difference method for compressible flows with discontinuities[J ] . Journal of Computational Physics, 2015 , 300 ( C ): 269 − 287 .

HE Z W, LI L, ZHANG Y S, et al . Consistent implementation of characteristic flux-split based finite difference method for compressible multi-material gas flows[J ] . Computers Fluids, 2018 , 168 : 190 − 200 .

XIAO F, HONMA Y, KONO T . A simple algebraic interface capturing scheme using hyperbolic tangent function[J ] . International Journal for N umerical Methods in Fluids, 2005 , 48 ( 9 ): 1023 − 1040 .

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).

HAN F, WEI Y X, LIU J . Numerical dissipation of upwind schemes in contact discontinuity and their induced error[J ] . Journal of National University of Defense Technology, 2024 , 46 ( 1 ): 51 − 62 (in Chinese).

STEGER J L, WARMING R F . Flux vector splitting of the inviscid gasdynamic equations with application to finite-difference methods[J ] . Journal of Computational Physics, 1981 , 40 ( 2 ): 263 − 293 .

WITHERDEN F D, JAMESON A . On the spectrum of the Steger-Warming flux-vector splitting scheme[J ] . International Journal for Numerical Methods in Fluids, 2018 , 87 ( 12 ): 601 − 606 .

LEVEQUE R J. Finite Volume methods for hyperbolic problems[M]. Cambridge, UK: Cambridge University Press, 2002.

HOU T Y, LE FLOCH P G . Why nonconservative schemes converge to wrong solutions: error analysis[J ] . Mathematics of Computation, 1994 , 62 ( 206 ): 497 − 530 .

GOTTLIEB S, SHU C W, TADMOR E . Strong stability-preserving high-order time discretization methods[J ] . SIAM Review, 2001 , 43 ( 1 ): 89 − 112 .

QUIRK J J . A contribution to the great Riemann solver debate[J ] . International Journal for Numerical Methods in Fluids, 1994 , 18 ( 6 ): 555 − 574 .

WANG Z J, FIDKOWSKI K, ABGRALL R, et al . High-order CFD methods: Current status and perspective[J ] . International Journal for Numerical Methods in Fluids, 2013 , 72 ( 8 ): 811 − 845 .

Views

下载量

CSCD

Alert me when the article has been cited

提交

Tools

Publicity Resources

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

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

Investigation on the aerodynamic characteristics of vertical axis airfoils and the influence mechanism of pitch moment

A study on fluctuating wind characteristics during strong wind processes in the surface layer over the Qinghai-Tibet Plateau, Xizang region

Wake model of a floating offshore wind turbine under surge motion

Related Author

WEI Yanxin

LIU Jun

HAN Fang

Related Institution

School of Aeronautics and Astronautics, Dalian University of Technology

AI问答

⁰