高精度有限差分格式的色散优化及耗散控制

李妍慧; 陈琮巍; 任玉新; 孙振生; 王秋菊

doi:10.7638/kqdlxxb-2020.0153

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

English Version 期刊登录

您当前的位置：

首页 >

文章列表页 >

高精度有限差分格式的色散优化及耗散控制

专栏——计算流体力学高精度格式及其应用 | 更新时间：2026-05-28

- 高精度有限差分格式的色散优化及耗散控制
- The dispersion optimization and dissipation adjustment for high-order finite difference schemes
- 空气动力学学报 2021年39卷第1期页码：138-156
- 作者机构：
  
  1. 清华大学航天航空学院,北京,100084
  2. 火箭军工程大学,西安,710025
  3. 北京应用物理与计算数学研究所,北京,100094
- 作者简介：
  
  李妍慧(1996-), 女, 山西人, 博士研究生, 研究方向: 计算流体力学.E-mail: liyanhui_lyhn@163.com
  任玉新(1966-), 男, 河北人, 博士, 教授, 研究方向: 计算流体力学.E-mail: ryx@tsinghua.edu.cn
- 基金信息：
  
  国家重点研发计划（2016YFA0401200）;国家重大项目（GJXM92579）;国家自然科学基金（91752114，11672160）
- DOI：10.7638/kqdlxxb-2020.0153
  中图分类号： V211.3
- 纸质出版：2021
- 稿件说明：
移动端阅览
李妍慧, 陈琮巍, 任玉新, 孙振生, 王秋菊. 高精度有限差分格式的色散优化及耗散控制[J]. 空气动力学学报, 2021, 39(1): 138-156.

LI Yanhui, CHEN Congwei, REN Yuxin, et al. The dispersion optimization and dissipation adjustment for high-order finite difference schemes[J]. 2021, 39(1): 138-156.
李妍慧, 陈琮巍, 任玉新, 孙振生, 王秋菊. 高精度有限差分格式的色散优化及耗散控制[J]. 空气动力学学报, 2021, 39(1): 138-156. DOI： 10.7638/kqdlxxb-2020.0153.

LI Yanhui, CHEN Congwei, REN Yuxin, et al. The dispersion optimization and dissipation adjustment for high-order finite difference schemes[J]. 2021, 39(1): 138-156. DOI： 10.7638/kqdlxxb-2020.0153.

摘要

本文综述了我们在高精度有限差分格式的色散优化和耗散控制方面的研究进展。首先，我们提出了半离散有限差分格式色散和耗散相互独立的充分条件，实现了差分格式色散和耗散特性的独立调节。在此基础上，提出了色散最小、耗散可控的高精度差分格式，称为MDCD格式。MDCD格式已经得到了广泛应用，取得了很好的计算效果，但其主要缺点是耗散的调节依赖于经验。为了解决这一问题，我们进一步提出耗散的自适应调节方法。具有自适应耗散特性的高精度有限差分格式的基本特征是，差分格式的耗散能够随解的局部尺度自适应调节。为了构造这类格式，我们提出了一种新型的尺度识别器，它能够以等效无量纲波数的形式来定量衡量数值解的局部长度尺度。在此基础上，设计差分格式耗散参数与尺度识别器得到的等效无量纲波数之间的关系，从而构造了一类色散最小、耗散自适应的差分格式，称为MDAD格式。为了计算含有间断的问题，同时保持在光滑区的良好耗散特性，我们利用尺度识别器对一种经典的激波探测器进行改进，提出了一种新的激波探测器，并将自适应耗散格式与对应的WENO格式相混合，得到自适应耗散混合格式。近似色散关系显示该混合格式兼具高分辨率和鲁棒性。多个含间断流场的标准算例测试结果显示，自适应耗散混合格式具有良好的分辨率和激波捕捉能力。

Abstract

In this paper

our recent work on the dispersion optimization and dissipation adjustment for high-order finite difference schemes is reviewed. For compressible flows with a broad range of length scales

good dispersion and dissipation properties are crucial for numerical schemes to realize high-fidelity simulations. It has been recognized that the minimization of dispersion error is an effective measure to improve the resolution of schemes

while the dissipation should be adjusted according to the local scale of the solution. We found a sufficient condition for a semi-discretized finite difference scheme to have independent dispersion and dissipation properties. Based on this finding

the high-order finite difference scheme with minimized dispersion and controllable dissipation (MDCD) was proposed. The MDCD scheme has been widely used and achieved good performance. However

one drawback is that the dissipation is controlled empirically. To realize the automatic adjustment of the dissipation according to the local scale of the solution

we devised a scale sensor to quantify the length scale of the numerical solution which can be calibrated to give the effective scaled wavenumbers. Then the dispersion-dissipation condition is used to construct the relationship between the dissipation parameter and the effective scaled wavenumber. A class of finite difference schemes with adaptive dissipation is thus obtained. To achieve the shock-capturing capability while maintaining the superior spectral properties in the smooth regions

we make an improvement to a classical shock detector by incorporating the proposed scale sensor

and construct a new hybrid scheme which consists of the adaptive dissipation scheme and the corresponding weighted essentially non-oscillation (WENO) scheme. The approximate dispersion relation (ADR) shows that the new hybrid scheme is accurate and robust. Several benchmark test problems with broadband length scales as well as discontinuities are presented to validate the high-resolution and the good shock-capturing capability of the proposed scheme.

关键词

Keywords

references

LIU X D, OSHER S, CHAN T. Weighted essentially non-oscillatory schemes[J]. Journal of Computational Physics, 1994, 115: 200-212.

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

JOHNSEN E, LARSSON J, BHAGATWALA A V, et al. Assessment of high-resolution methods for numerical simulations of compressible turbulence with shock waves[J]. Journal of Computational Physics, 2010, 229(4): 1213-1237.

HENRICK A K, ASLAM T A, POWERS J M. Mapped weighted essentially non-oscillatory schemes: achieving optimal order near critical points[J]. Journal of Computational Physics, 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: 3191-3211.

LIU P X, LI Q, ZHANG H X. A kind of third order central scheme based on mapping functions[J]. Acta Aerodynamics Sinica, 2017, 35(1): 71-77.

TAYLOR E M, WU M, MARTIN M P. Optimization of nonlinear error for weighted essentially non-oscillatory methods in direct numerical simulations of compressible turbulence[J]. Journal of Computational Physics, 2007, 223: 384-397.

MARTIN M P, TAYLOR E M, WU M, et al. A bandwidth-optimized WENO scheme for effective direct numerical simulation of compressible turbulence[J]. Journal of Computational Physics, 2006, 220: 270-289.

HU X Y, WANG Q, ADAMS N A. An adaptive central-upwind weighted essentially non-oscillatory scheme[J]. Journal of Computational Physics, 2010, 229: 8952-8965.

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.

ADAMS N A, SHARIFF K. A high-resolution hybrid compact-ENO scheme for shock/turbulence interaction problems[J]. Journal of Computational Physics, 1996, 127: 27-51.

PIROZZOLI S. Conservative hybrid compact-WENO schemes for shock-turbulence interaction[J]. Journal of Computational Physics, 2002, 178: 81-117.

REN Y X, LIU M, ZHANG H X. A characteristic-wise compact WENO scheme for solving hyperbolic conservation laws[J]. Journal of Computational Physics, 2003, 192: 365-386.

WU C H, ZHAO N, TIAN L L. An improved hybrid compact-WENO scheme[J]. Acta Aerodynamics Sinica, 2013, 31(4): 477-481. (in Chinese)

HARTEN A. The artificial compression method for computation of shocks and contact discontinuities. Ⅲ. Self-adjusting hybrid schemes[J]. Mathematics of Computation, 1978, 32(142): 363-389.

JAMESON A, SCHMIDT W, TURKEL E. Numerical solution of the Euler equations by finite volume methods using Runge-Kutta time stepping schemes[C]//Proceeding of the 14th AIAA Fluid and Plasma Dynamics Conference, 1981.

SJOGREEN B, YEE H C. Multiresolution wavelet based adaptive numerical dissipation control for shock-turbulence computations[J]. Research Institute for Advanced Computer Science, 2001, 1: 1-43.

HILL D J, PULLIN D I. Hybrid tuned center-difference-WENO method for large eddy simulations in the presence of strong shocks[J]. Journal of Computational Physics, 2004, 194(2): 435-450.

VISBAL M, GAITONDE D. Shock capturing using compact-differencing-based methods[C]//Proceeding of the 43rd AIAA Aerospace Sciences Meeting and Exhibit, 2005.

MOVAHED P, JOHNSEN E. A solution-adaptive method for efficient compressible multifluid simulations with application to the Richtmyer-Meshkov instability[J]. Journal of Computational Physics, 2013, 239: 166-186.

ZHAO G Y, SUN M B, PIROZZOLI S. On shock sensors for hybrid compact/WENO schemes[J]. Computers and Fluids, 2020, 199: 104439.

SUN Z S, REN Y X, LARRICQ C, et al. A class of finite difference scheme with low dispersion and controllable dissipation for DNS of compressible turbulence[J]. Journal of Computational Physics, 2011, 230: 4616-4635.

SUN Z S, LUO L, REN Y X, et al. A sixth order hybrid finite difference scheme based on the minimized dispersion and controllable dissipation technique[J]. Journal of Computational Physics, 2014, 270: 238-254.

LI Y H, CHEN C W, REN Y X. A class of high-order finite difference schemes with minimized dispersion and adaptive dissipation for solving compressible flows[J]. Submitted to Journal of Computational Physics.

PIROZZOLI S. On the spectral properties of shock-capturing schemes[J]. Journal of Computational Physics, 2006, 219: 489-497.

HU F Q, HUSSAINI M Y, MANTHEY J L. Low-dissipation and low-dispersion Runge-Kutta schemes for computational acoustics[J]. Journal of Computational Physics, 1996, 124: 177-191.

LI L, YU C P, CHEN Z, et al. Resolution-optimised nonlinear scheme for secondary derivatives[J]. International Journal of Computational Fluid Dynamics, 2016, 30(2): 107-119.

SHU C W. High order ENO and WENO schemes for computational fluid dynamics[M]. Berlin: Springer, 1999, 439-582.

SUN Z S, REN Y X. The development of the characteristic-wise hybrid compact-WENO scheme for solving the Euler and Navier-Stokes equations[J]. Computational Fluid Dynamics Review 2010, 2010, 1: 279-297.

SHEN Y Q, YANG B Y, ZHA G C. Implicit WENO scheme and high order viscous formulas for compressible flows[C]//Proceeding of the 25th AIAA Applied Aerodynamics Conference, 2007.

SHU C W. Advanced numerical approximation of nonlinear hyperbolic equations[M]. Berlin: Springer, 2006, 325-432.

WANG Q J, REN Y X, SUN Z S, et al. Low dispersion finite volume scheme based on reconstruction with minimized dispersion and controllable dissipation[J]. Science China Physics, Mechanics Astronomy, 2013, 56: 423-431.

WANG Q J, DEITERDING R, PAN J H, et al. Consistent high resolution interface-capturing finite volume method for compressible multi-material flows[J]. Computers and Fluids, 2020, 202: 104518.

ZENG W G, PAN J H, SUN Y T, et al. Turbulent mixing and energy transfer of reshocked heavy gas curtain[J]. Physics of Fluids, 2018, 30(6): 064106.

LI Z, ZHANG Y, CHEN H. A low dissipation numerical scheme for implicit large eddy simulation[J]. Computers Fluids, 2015, 117: 233-246.

CHEN H, LI Z, ZHANG Y. U or V shape: dissipation effects on cylinder flow implicit large-eddy simulation[J]. AIAA Journal, 2017, 55(2): 459-473.

JIANG S, FU S. Modifications to the SIMPLE algorithm with the MDCD approach for incompressible flow simulation[J]. International Journal of Numerical Methods for Heat Fluid Flow, 2018.

DUAN Z W, XIAO Z X, FU S. Direct numerical simulation of hypersonic transition induced by an isolated cylindrical roughness element[J]. SCIENCE CHINA Physics, Mechanics Astronomy, 2014, 57(12): 2330-2345.

DUAN Z W, XIAO Z X. Roughness element induced hypersonic boundary layer transition[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(8): 2454-2463. (in Chinese)

JIANG S Y, FU S. Numerical investigation of separation-induced transition in a low-pressure turbine cascade in a low-disturbance environment[J]. SCIENCE CHINA Physics, Mechanics Astronomy, 2020, 63(6): 1-17.

YANG J, ZHANG Y, CHEN H, et al. Unsteady flow control of a plane diffuser based on a Karman-vortex generator[J]. AIP Advances, 2020, 10(5): 055314.

WOODWARD P, COLELLA P. The numerical simulation of two-dimensional fluid flow with strong shocks[J]. Journal of Computational Physics, 1984, 54: 115-173.

TENAUD C, GARNIER E, SAGAUT P. Evaluation of some high-order shock capturing schemes for direct numerical simulation of unsteady two-dimensional free flows[J]. International Journal for Numerical Methods in Fluids, 2000, 33: 249-278.

浏览量

下载量

CSCD

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

提交

工具集

关联资源

非结构网格上的TV-van Leer混合迎风格式

基于LBM方法的LAGOON起落架噪声预测与降噪方法研究

某型电动跑车尾翼气动设计的仿真与试验研究

水滴撞击运动超疏水表面特性及其调控

风扇叶片颗粒侵蚀特性的多因素分析与模拟