边界层内WENO格式特征重构精度的理论分析

doi:10.7527/S1000-6893.2021.26362

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

前一篇 | 后一篇

边界层内WENO格式特征重构精度的理论分析

许珂¹, 高振勋^1,2, 蒋崇文^1,2, 李椿萱^1,2

1. 北京航空航天大学航空科学与工程学院, 北京 100083;
2. 中国航天科技集团航天飞行器气动热防护实验室, 北京 100191

收稿日期:2021-09-01 修回日期:2021-09-14 发布日期:2021-09-18
通讯作者: 高振勋 E-mail:gaozhenxun@buaa.edu.cn
基金资助:
国家数值风洞（NNW2018-ZT3A05）

Theoretical analysis of reconstruction accuracy of characteristic-wise WENO scheme within boundary layer

XU Ke¹, GAO Zhenxun^1,2, JIANG Chongwen^1,2, LEE Chunhian^1,2

1. School of Aeronautic Science and Engineering, Beihang University, Beijing 100083, China;
2. Aerospace Vehicle Aerodynamic Thermal Protection Laboratory, China Aerospace Science and Technology Corporation, Beijing 100191, China

Received:2021-09-01 Revised:2021-09-14 Published:2021-09-18
Supported by:
National Numerical Windtunnel Project (NNW2018-ZT3A05)

摘要/Abstract

摘要： WENO格式因其高精度以及良好的激波捕捉能力而被广泛地用于可压缩湍流的数值模拟过程中。为了在有强间断的流场中仍然能保证数值稳定，通常将流动物理量投影到特征空间后再进行WENO重构。然而这种特征重构方法在边界层内的精度问题仍需要进一步分析与研究。通过可压缩边界层自相似解构造平板边界层内法向的流动参数，并在该方向上进行WENO重构，得到求解半点数值通量的WENO非线性权重。通过分析实际得到的权重与理想权重的偏差大小发现在边界层流动中特征重构的精度明显低于分量重构。为了更深入地分析特征重构精度降低的原因，通过对曲线坐标系下的流动方程进行理论分析，推导得到特征重构时的特征变量形式。根据特征变量表达式以及自相似解中权重误差分布发现引起特征重构精度降低的主要原因是特征投影过程在半点处出现了额外的极值点，同时这种因特征投影过程而产生极值的特性将不会随着左右特征矩阵的不同选取而发生改变。进一步地，基于理论分析的普适性，对于任何光滑流场，只要使用特征重构的投影过程，都会在半点处产生极值从而使得WENO系列格式产生精度降低的现象。

关键词: WENO格式, 特征重构, 分量重构, 格式精度, 可压缩自相似解

Abstract: Weighted Essentially Non-Oscillatory (WENO) scheme has been widely used in the numerical simulation of compressible turbulent flows due to its high accuracy and excellent shock-capture capability. In order to ensure numerical stability in flow fields with strong discontinuities, the physical quantities of the flow are usually projected into the characteristic space before the WENO reconstruction. In this paper, we construct the flow parameters normal to the plate boundary layer by compressible boundary layer self-similar solutions and perform WENO reconstruction in this direction to obtain the nonlinear weights for solving the half-point numerical fluxes. By analyzing the magnitude of deviation of the actual obtained weights from the ideal weights, it is found that the accuracy of the characteristic-wise reconstruction in the boundary layer flow is significantly lower than that of the component-wise reconstruction. In order to analyze more deeply the reasons for the reduced accuracy of the characteristic reconstruction, the form of the characteristic variables in the characteristic reconstruction process has been derived by theoretical analysis of the flow equations in the curved coordinate system. Based on the expressions of the characteristic variables and the weight error distribution in the norm direction, it is found that the main cause for the degradation of the accuracy is the appearance of additional extreme points at the half-point of the projection process, and this characteristic of extreme values due to the projection process will not change with the different selection of the left and right eigenmatrices. Further, based on the generality of the theoretical analysis, for any smooth flow field, as long as the projection process of characteristic reconstruction is used, the extreme points will cause the accuracy of WENO scheme.

Key words: WENO scheme, characteristic-wise reconstruction, component-wise reconstruction, accuracy of method, compressible self-similar solution

中图分类号:

V211.3

许珂, 高振勋, 蒋崇文, 李椿萱. 边界层内WENO格式特征重构精度的理论分析[J]. 航空学报, 2021, 42(S1): 726362-726362.

XU Ke, GAO Zhenxun, JIANG Chongwen, LEE Chunhian. Theoretical analysis of reconstruction accuracy of characteristic-wise WENO scheme within boundary layer[J]. ACTA AERONAUTICAET ASTRONAUTICA SINICA, 2021, 42(S1): 726362-726362.

参考文献

[1] JIANG G S, SHU C W. Efficient implementation of weighted ENO schemes[J]. Journal of Computational Physics, 1996, 126(1):202-228.
[2] REN Y X, LIU M E, ZHANG H X. A characteristic-wise hybrid compact-WENO scheme for solving hyperbolic conservation laws[J]. Journal of Computational Physics, 2003, 192(2):365-386.
[3] PENG J, ZHAI C L, NI G X, et al. An adaptive characteristic-wise reconstruction WENO-Z scheme for gas dynamic Euler equations[J]. Computers & Fluids, 2019, 179:34-51.
[4] 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:269-287.
[5] 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.
[6] BROWNING G L, KREISS H O. Comparison of numerical methods for the calculation of two-dimensional turbulence[J]. Mathematics of Computation, 1989, 52(186): 369.
[7] WEIRS V, CANDLER G, WEIRS V, et al. Optimization of weighted ENO schemes for DNS of compressible turbulence[C]//13th Computational Fluid Dynamics Conference. Reston: AIAA, 1997.
[8] DUMBSER M, BOSCHERI W, SEMPLICE M, et al. Central weighted ENO schemes for hyperbolic conservation laws on fixed and moving unstructured meshes[J]. SIAM Journal on Scientific Computing, 2017, 39(6): A2564-A2591.
[9] BALSARA D S, SHU C W. Monotonicity preserving weighted essentially non-oscillatory schemes with increasingly high order of accuracy[J]. Journal of Computational Physics, 2000, 160(2):405-452.
[10] QIU J X, 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.
[11] PIROZZOLI S, GRASSO F, GATSKI T B. Direct numerical simulation and analysis of a spatially evolving supersonic turbulent boundary layer at M=2.25[J]. Physics of Fluids, 2004, 16(3):530-545.
[12] SHU C W. Essentially non-oscillatory and weighted essentially non-oscillatory schemes[J]. Acta Numerica, 2020, 29:701-762.
[13] LIU X D, OSHER S, CHAN T. Weighted essentially non-oscillatory schemes[J]. Journal of Computational Physics, 1994, 115(1):200-212.
[14] 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.
[15] 刘巍,张理论,王勇献,等. 计算空气动力学并行编程基础[M].北京:国防工业出版社, 2013:13-19. LIU W, ZHANG L L, WANG Y X, et al. Foundations of computational aerodynamics parallel programming[M].Beijing: National Defense Industry Press, 2013: 13-19 (in Chinese).
[16] PUPPO G. Adaptive application of characteristic projection for central schemes[M]//Hyperbolic Problems: Theory, Numerics, Applications. Berlin: Springer Berlin Heidelberg, 2003: 819-829.
[17] 傅德薰,马延文,李新亮,等. 可压缩湍流直接数值模拟[M].北京:科学出版社,2010:64-67. FU D X, MA Y W, LI X L, et al. Direct numerical simulation of compressible turbulence[M].Beijing: Science Press, 2010: 64-67(in Chinese).
[18] 水鸿寿. 一维流体力学差分方法[M].北京:国防工业出版社, 1998: 86-103. SHUI H S. One-dimensional hydrodynamic difference method[M].Beijing: National Defense Industry Press, 1998: 86-103 (in Chinese).
[19] ANDERSON J D Jr. Hypersonic and high-temperature gas dynamics, third edition[M].Reston: AIAA, 2019.
[20] DON W S, BORGES R. Accuracy of the weighted essentially non-oscillatory conservative finite difference schemes[J]. Journal of Computational Physics, 2013, 250:347-372.

E-mail：hkxb@buaa.edu.cn

关于我们

期刊社服务

专业学科

封面文章

友情链接

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

边界层内WENO格式特征重构精度的理论分析

Theoretical analysis of reconstruction accuracy of characteristic-wise WENO scheme within boundary layer

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 4

编辑推荐

Metrics

本文评价

[1]	刘君, 韩芳, 魏雁昕. 应用维数分裂方法推广MUSCL和WENO格式的若干问题[J]. 航空学报, 2022, 43(3): 125009-125009.
[2]	刘君, 韩芳, 魏雁昕. 特定条件下高阶WENO格式计算结果误差[J]. 航空学报, 2022, 43(2): 124940-124940.
[3]	李超群, 李易, 张晨曦, 唐硕. 沿流向微结构沟槽流场直接数值模拟[J]. 航空学报, 2020, 41(11): 123628-123628.
[4]	叶舟, 徐国华, 史勇杰. 旋翼桨尖涡生成及演化机理的高精度数值研究[J]. 航空学报, 2017, 38(7): 520846-520846.