一种数值稳定的低耗散激波捕捉格式

胡立军; 李龙

doi:10.7638/kqdlxxb-2023.0062

欢迎来到能源电力期刊网！

English Version 期刊登录

您当前的位置：

首页 >

文章列表页 >

一种数值稳定的低耗散激波捕捉格式

研究论文 | 更新时间：2026-07-27

- 一种数值稳定的低耗散激波捕捉格式
- A numerically stable and low dissipation shock-capturing scheme
- 空气动力学学报 2024年42卷第7期页码：35-47
- 作者机构：
  
  衡阳师范学院数学与统计学院,衡阳,421002
- 作者简介：
  
  胡立军*（1985—），男，湖南常德人，副教授，研究方向：计算流体力学. E-mail：hulijun@lsec.cc.ac.cn
- 基金信息：
  
  湖南省自然科学基金（2021JJ40009, 2021JJ30056 ）;湖南省教育厅资助科研项目（21B0626）
- DOI：10.7638/kqdlxxb-2023.0062
  中图分类号： O354.4
- 纸质出版：2024
- 稿件说明：
移动端阅览
胡立军, 李龙. 一种数值稳定的低耗散激波捕捉格式[J]. 空气动力学学报, 2024, 42(7): 35−47.

HU Lijun, LI Long. A numerically stable and low dissipation shock-capturing scheme[J]. 2024, 42(7): 35-47.
胡立军, 李龙. 一种数值稳定的低耗散激波捕捉格式[J]. 空气动力学学报, 2024, 42(7): 35−47. DOI： 10.7638/kqdlxxb-2023.0062

HU Lijun, LI Long. A numerically stable and low dissipation shock-capturing scheme[J]. 2024, 42(7): 35-47. DOI： 10.7638/kqdlxxb-2023.0062.

摘要

HLLC格式应用广泛且可以准确捕捉激波间断，其满足熵条件和保正性，易于推广到其他双曲系统。但是，在计算多维强激波时出现的数值不稳定现象大大限制了它在高马赫数流动问题中的应用。本文基于三维无黏可压缩欧拉方程组构造一种激波数值稳定的HLLC格式（SS-HLLC）：将传统的HLLC通量显式地改写成耗散的HLL通量和反扩散项的形式，利用充分光滑的控制函数调节激波附近横向通量中的反扩散项来改善格式的激波数值稳定性，分析证明了反扩散控制技术的有效性。此外，采用简单的THINC重构方法和BVD算法来减少数值耗散项中的密度差，从而构造了一种低耗散的激波数值稳定格式（LD-SS-HLLC），可进一步提高格式对接触间断的分辨率。经过一系列标准算例测试，结果证明了LD-SS-HLLC格式的低耗散性和激波数值稳定性。

Abstract

It is well known that the Harten-Lax-van Leer with contact (HLLC) scheme suffers from unstable performance in capturing strong multi-dimensional shock waves

which impedes its application in high-Mach-number flows. This paper presents a simple strategy to develop a stable HLLC-type scheme for three-dimensional inviscid Euler equations. The conventional HLLC flux is split into a dissipative HLL flux and an anti-diffusion term. The improvement of the numerical stability

justified by stability analyses and numerical experiments

is achieved by using a sufficiently smooth control function to adjust the numerical dissipation term of the transverse flux adjacent to the shock wave. In addition

the simple THINC (tangent of hyperbola for interface capturing) reconstruction method and the BVD (boundary variation diminishing) algorithm are employed to reduce the density difference in the numerical dissipation term

thus further improving the resolution for contact discontinuities. The numerical results of one- to three-dimensional benchmark examples demonstrate the proposed scheme’s high resolution for contact discontinuities and stability for multi-dimensional strong shock waves. The numerical results of a shock/boundary-layer interaction problem promise the future application of the proposed local anti-diffusion control method in simulating viscous flows. The modularization of the proposed scheme and its integration into Computational Fluid Dynamics software to simulate

for instance

hypersonic flows

multi-component chemical reaction flows

and turbulent combustion are worthy of further research in the future.

关键词

Keywords

references

GRESSIER J, MOSCHETTA J M . Robustness versus accuracy in shock-wave computations[J ] . International Journal for Numerical Methods in Fluids, 2000 , 33 ( 3 ): 313 − 332 .

HU L J, FENG Z S . Development of a carbuncle-free and low-dissipation Roe-type scheme: Applications to multidimensional Euler flows[J ] . Communications in Nonlinear Science and Numerical Simulation, 2023 , 116 : 106798 .

DUMBSER M, MOSCHETTA J M, GRESSIER J . A matrix stability analysis of the carbuncle phenomenon[J ] . Journal of Computational Physics, 2004 , 197 ( 2 ): 647 − 670 .

KITAMURA K, SHIMA E, ROE P L . Carbuncle phenomena and other shock anomalies in three dimensions[J ] . AIAA Journal, 2012 , 50 ( 12 ): 2655 − 2669 .

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

KIM S D, LEE B J, LEE H J, et al . Robust HLLC Riemann solver with weighted average flux scheme for strong shock[J ] . Journal of Computational Physics, 2009 , 228 ( 20 ): 7634 − 7642 .

SHEN Z J, YAN W, YUAN G W . A robust HLLC-type Riemann solver for strong shock[J ] . Journal of Computational Physics, 2016 , 309 : 185 − 206 .

FU L, GAO Z H, ZHANG X D . Construction and application research of HLL-HLLC scheme[J ] . Acta Aerodynamica Sinica, 2014 , 32 ( 1 ): 116 − 122 (in Chinese).

REN Y X . A robust shock-capturing scheme based on rotated Riemann solvers[J ] . Computers Fluids, 2003 , 32 ( 10 ): 1379 − 1403 .

NISHIKAWA H, KITAMURA K . Very simple, carbuncle-free, boundary-layer-resolving, rotated-hybrid Riemann solvers[J ] . Journal of Computational Physics, 2008 , 227 ( 4 ): 2560 − 2581 .

CHEN Z Q, HUANG X D, REN Y X, et al . General procedure for Riemann solver to eliminate carbuncle and shock instability[J ] . AIAA Journal, 2017 , 55 ( 6 ): 2002 − 2015 .

SIMON S, MANDAL J C . A simple cure for numerical shock instability in the HLLC Riemann solver[J ] . Journal of Computational Physics, 2019 , 378 : 477 − 496 .

CHEN S S, LIN B X, LI Y S, et al . HLLC+: Low-Mach shock-stable HLLC-type Riemann solver for all-speed flows[J ] . SIAM Journal on Scientific Computing, 2020 , 42 ( 4 ): B921 − B950 .

CHEN S S, YAN C, LIN B X, et al . Affordable shock-stable item for Godunov-type schemes against carbuncle phenomenon[J ] . Journal of Computational Physics, 2018 , 373 : 662 − 672 .

RODIONOV A V . Simplified artificial viscosity approach for curing the shock instability[J ] . Computers Fluids, 2021 , 219 : 104873 .

TAN S D, HU L J, YUAN H Z . Development of a shock-stable and contact-preserving scheme for multidimensional Euler equations[J ] . AIAA Journal, 2022 , 60 ( 9 ): 5232 − 5248 .

BALSARA D S . A two-dimensional HLLC Riemann solver for conservation laws: Application to Euler and magnetohydrodynamic flows[J ] . Journal of Computational Physics, 2012 , 231 ( 22 ): 7476 − 7503 .

HU L J, YUAN L, ZHAO K L . Development of accurate and robust genuinely two-dimensional HLL-type Riemann solver for compressible flows[J ] . Computers Fluids, 2020 , 213 : 104719 .

FLEISCHMANN N, ADAMI S, ADAMS N A . A shock-stable modification of the HLLC Riemann solver with reduced numerical dissipation[J ] . Journal of Computational Physics, 2020 , 423 : 109762 .

KEMM F . Heuristical and numerical considerations for the carbuncle phenomenon[J ] . Applied Mathematics and Computation, 2018 , 320 : 596 − 613 .

HU L J, YUAN H Z, ZHAO K L . A shock-stable HLLEM scheme with improved contact resolving capability for compressible Euler flows[J ] . Journal of Computational Physics, 2022 , 453 : 110947 .

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

SUN Z Y, INABA S, XIAO F . Boundary variation diminishing (BVD) reconstruction: A new approach to improve Godunov schemes[J ] . Journal of Computational Physics, 2016 , 322 : 309 − 325 .

TORO E F . The HLLC Riemann solver[J ] . Shock Waves, 2019 , 29 ( 8 ): 1065 − 1082 .

TORO E F, SPRUCE M, SPEARES W . Restoration of the contact surface in the HLL-Riemann solver[J ] . Shock Waves, 1994 , 4 ( 1 ): 25 − 34 .

COLELLA P . A direct eulerian MUSCL scheme for gas dynamics[J ] . SIAM Journal on Scientific and Statistical Computing, 1985 , 6 ( 1 ): 104 − 117 .

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

GARG N K, MARUTHI N H, RAGHURAMA RAO S V, et al . Use of Jordan forms for convection-pressure split Euler solvers[J ] . Journal of Computational Physics, 2020 , 407 : 109258 .

SAN O, STAPLES A E . Stationary two-dimensional turbulence statistics using a Markovian forcing scheme[J ] . Computers Fluids, 2013 , 71 : 1 − 18 .

HUANG K B, WU H, YU H, et al . Cures for numerical shock instability in HLLC solver[J ] . International Journal for Numerical Methods in Fluids, 2011 , 65 ( 9 ): 1026 − 1038 .

浏览量

下载量

CSCD

文章被引用时，请邮件提醒。

提交

工具集

关联资源

基于CFD仿真的漂浮式风力机尾流模型研究

考虑来流风速特征的风电机组超短期功率预测方法

考虑海气耦合的海洋风环境中-微尺度数值模拟

差分格式应用于非均匀网格的精度保持算法

纵荡运动下漂浮式风力机尾流模型