高精度WCNS格式加权策略改进与数值验证

doi:10.16183/j.cnki.jsjtu.2022.014

高精度WCNS格式加权策略改进与数值验证

杨强, 李伟鹏^,

上海交通大学航空航天学院,上海 200240

Improvement and Numerical Verification of Weighting Strategy for High Precision WCNS Scheme

YANG Qiang, LI Weipeng^,

School of Aeronautics and Astronautics, Shanghai Jiao Tong University, Shanghai 200240, China

通讯作者: 李伟鹏,教授,博士生导师,电话(Tel.):021-34208076;E-mail:liweipeng@sjtu.edu.cn.

责任编辑: 王一凡

收稿日期: 2022-01-17 修回日期: 2022-06-14 接受日期: 2022-07-12

基金资助:

国家自然科学基金(92052101)

Received: 2022-01-17 Revised: 2022-06-14 Accepted: 2022-07-12

作者简介 About authors

杨强(1995-),硕士生,从事飞行器设计研究.

摘要

为了揭示复杂的流动机理,国内外提出了一系列高阶精度格式,其中加权非线性紧致格式(WCNS)具有良好的激波捕捉能力,已广泛用于复杂流动数值模拟,然而在模拟小尺度流动时分辨率不足、耗散偏大.在WCNS的框架下,借鉴基于目标本质无振荡(TENO)格式的加权策略,将间断检测和模板加权的新方法引入WCNS构造中,发展了一种7阶精度的WCNS7-T格式,通过一维激波管问题和二维黎曼问题等开展算例测试,并与传统的WCNS7-Z格式计算结果进行对比来验证新格式的改善性能.数值实验表明,WCNS7-T格式更好地抑制了间断附近的数值振荡,在提高分辨率和激波捕捉能力的同时也进一步降低了耗散.

关键词： 可压缩流动; 高精度格式; 激波间断; 非线性加权

Abstract

In order to reveal the complex flow mechanism, a series of high-order precision schemes have been proposed at home and abroad, of which, the weighted compact nonlinear scheme (WCNS) has a good shock capture ability and has been widely used in the numerical simulation of complex flows. However, it has insufficient resolution and large dissipation in the simulation of small-scale flows. In the framework of the WCNS, by using the weighting strategy of the targeted essentially non-oscillatory (TENO) scheme for reference, this paper introduces the new methods of discontinuity detection and template weighting into the construction of the WCNS scheme, and develops a WCNS7-T scheme with a 7-order accuracy. Example tests are conducted through one-dimensional shock tube problem and two-dimensional Riemann problem. By comparing with the traditional WCNS7-Z scheme, the improved performance of the new scheme is verified. The numerical experiments show that the WCNS7-T scheme can better suppress the numerical oscillation near the discontinuity, improve the resolution and shock capture ability, and further reduce the dissipation.

Keywords： compressible flow; high precision scheme; shock discontinuity; nonlinear weighting

PDF (3360KB) 元数据多维度评价相关文章导出 EndNote| Ris| Bibtex 收藏本文

本文引用格式

杨强, 李伟鹏. 高精度WCNS格式加权策略改进与数值验证[J]. 上海交通大学学报, 2023, 57(6): 719-727 doi:10.16183/j.cnki.jsjtu.2022.014

YANG Qiang, LI Weipeng. Improvement and Numerical Verification of Weighting Strategy for High Precision WCNS Scheme[J]. Journal of Shanghai Jiaotong University, 2023, 57(6): 719-727 doi:10.16183/j.cnki.jsjtu.2022.014

可压缩流动广泛存在于超声速/高超声速飞行器中,激波间断和湍流共存且相互耦合作用给高可靠性数值模拟带来了巨大挑战.近40年,人们发展了众多高精度激波捕获方法,例如总变差不增(TVD)方法^[1]、ASUM^[2]、基本无振荡/加权基本无振荡(ENO/WENO)方法^[3⇓-5]、加权非线性紧致格式(WCNS)方法^[6]、基于目标本质无振荡(TENO)方法^[7]和AFTENO^[8]等.这些格式在激波间断附近增加局部耗散来抑制非物理数值振荡,而在光滑湍流区域减少耗散以保持对多尺度湍流结构的分辨率.经典激波捕捉格式,如TVD方法和ENO方法,虽然可以解决激波间断无数值振荡,但对湍流光滑区域产生了过多的数值耗散.WENO方法在ENO方法的基础上,引入子模板光滑因子,通过非线性加权方法来动态调整局部耗散^[5],使得高精度数值格式得到进一步发展.Borges等^[9]在非线性权值中引入更高阶的光滑指标,提出计算成本更低、精度更高的WENO-Z格式.虽然WENO格式的精度可以设计为任意高阶^{[10⇓⇓-13]},但WENO格式对光谱区域的分辨率仍然不能令人满意^[14],它的加权机制和光滑度指标设定依旧对小尺度湍流结构中引入过多的数值耗散^[15].1999年,Deng等^[16]提出非线性紧致格式(CNS),后续引入加权技术,构造了 WCNS^[6].近年来,Fu等^[7]提出了另一类非线性加权框架,即TENO 方法,通过类ENO方法的模板选择来控制非线性耗散,在无间断区域(包括光滑临界点)保证最优的高阶空间重构,在间断区域通过移除跨越间断的模板来避免数值振荡,这种方法有利于候选模板的使用和高阶格式的改进.当子模板中存在多个相邻的不连续间断,该方法可退化为三阶重构.即使对于高阶版本,数值鲁棒性也得到了明显改善,且避免了阶跃退化问题.Peng等^[17]在此基础上进一步改进了TENO格式的C_T自适应函数策略,发展了TENO5-LAD格式,实验结果证明新方法能够更好地抑制不连续点附近的数值振荡,并能以较低的额外计算成本保持TENO的低耗散特性.而Ye等^[8]结合AFWENO格式和TENO格式的优点构建了可以实现任意高阶精度的AFTENO格式,在激波捕获和波数分辨率方面,相较于AFWENO-JS和AFWENO-Z有着明显的优势.之后,在加权紧致格式方面,Zhang等^[18]将WCNS与TENO格式组合发展了TCNS,实验结果显示TCNS能够在模拟中捕获更丰富的波结构,特别是在对瑞利-泰勒不稳定性问题的模拟中,TCNS在较粗网格上取得了与WCNS-JS^[5]相似甚至更好的结果.

将TENO格式的模板选择和加权策略引入到WCNS的构造过程,发展了WCNS-T格式,并在此基础上测验7阶精度下的数值表现情况.新格式在扩展模板宽度的同时改变了模板选择过程,当全局模板存在多个间断时,WCNS7-T格式将不会使用所有候选子模板,通过引入截断函数,将光滑和非光滑候选模板完全分离,并将非光滑候选模板的权重直接设为0,彻底剔除非光滑模板的影响.该模板策略能使格式在高阶精度时相较原来能保持更高的数值稳健性.同时,光滑度因子的新评测策略较之前能更准确判定子模板的光滑程度,降低了TENO格式中复杂光滑指标所带来的不必要的计算量,提高了格式捕捉激波的能力.一维和二维数值算例结果表明,WCNS7-T格式相比于传统的WCNS7-Z格式具有更高的精度和分辨率.

1 数值方法

1.1 控制方程

针对黎曼问题,差分格式的推导仅考虑一维守恒型Euler方程:

(1)∂Q/∂t+∂E/∂x=0

守恒变量Q和E为

(2) $\begin{aligned}&\boldsymbol{Q}=\left.\begin{array}{c|cc}\rho&&\rho\\\rho u&,&\boldsymbol{E}=&\left|\begin{array}{c}\rho u\\\rho u^2+\rho\end{array}\right|\end{array}\right|\\&e=\frac p{\gamma-1}+\frac{\rho u^2}2\end{aligned}$

式中:ρ为密度;u为来流速度;p为压力;γ为绝热指数.

在每个网格节点x_i处,对式(1)进行有限差分格式离散得:

(3)

{(\partial Q / \partial t)}_{i}

=-E'_i

式中:E'_i为空间导数的近似值,采用显式的节点和半节点混合通量求解方法^[19],

(4)E'_i=

\frac{α}{h}

(

{\hat{E}}_{i + 1 / 2}

{\hat{E}}_{i - 1 / 2}

$\frac{192 - 175 α}{256 h}$ ( ${\hat{E}}_{i + 1}$ - ${\hat{E}}_{i - 1}$ )+

$\frac{35 α - 48}{320 h}$ ( ${\hat{E}}_{i + 2}$ - ${\hat{E}}_{i - 2}$ )+

\frac{64 - 45 α}{3840 h}

(

{\hat{E}}_{i + 3}

{\hat{E}}_{i - 3}

)

式中:α为常数;h为网格单元长度; ${\hat{E}}_{i}$ 为节点边界上的通量.残差主项为 $(\frac{1}{8960} - \frac{5}{65536} α)$ h⁶ $E_{i}^{(6)}$ ,当α=175/256时格式具有8阶精度.对于半节点边界上的通量 ${\hat{E}}_{i \pm 1 / 2}$ =E( ${\hat{Q}}_{i \pm 1 / 2}$ )可以使用多种通量差分或分裂方法求解,如Roe方法^[20]、AUSM^[21]、van leer方法^[22]等. ${\hat{Q}}_{i \pm 1 / 2}$ 是Q_i_±1_/₂的高阶近似,可由紧致插值格式或显式插值格式得到.

1.2 WCNS7-Z格式构造

WCNS是在CNS基础上构造的,在解决CNS三对角线反演效率低于显式格式的问题同时,使用类似WENO格式加权的思想,充分利用计算模板的信息,对光滑模板和包含间断的模板分配对应的权值系数,从而使格式更高效.WCNS的构造方法为3个步骤:① 特征变量节点到半节点的非线性加权插值;② 半节点的通量求解;③ 节点-半节点到节点的混合差分求导.重点介绍第1步骤中的非线性加权插值过程.

与5阶格式的模板不同,在7阶WCNS^[19]的构造中使用4个4点子模板:(x_i_-3, x_i_-2, x_i_-1, x_i),(x_i_-2, x_i_-1, x_i, x_i₊₁),(x_i, x_i₊₁, x_i₊₂, x_i₊₃),(x_i, x_i₊₁, x_i₊₂, x_i₊₃),如图1所示.图中:L为全局模板;K₁、K₂、K₃、K₄分别是4个子模板.每个子模板上的4阶紧致格式设定为特定偏置形式:

(5)

\begin{array}{l} 5 {\hat{Q}}_{i - 1 / 2} + {\hat{Q}}_{i + 1 / 2} = - \frac{1}{4} {\hat{Q}}_{i - 2} + \frac{5}{2} {\hat{Q}}_{i - 1} + \frac{15}{4} {\hat{Q}}_{i} \\ {\hat{Q}}_{i - 1 / 2} + {\hat{Q}}_{i + 1 / 2} = \frac{1}{4} {\hat{Q}}_{i - 1} + \frac{3}{2} {\hat{Q}}_{i} + \frac{1}{4} {\hat{Q}}_{i + 1} \\ {\hat{Q}}_{i + 1 / 2} + \frac{1}{5} {\hat{Q}}_{i + 3 / 2} = - \frac{1}{20} {\hat{Q}}_{i - 1} + \frac{1}{2} {\hat{Q}}_{i} + \frac{3}{4} {\hat{Q}}_{i + 1} \\ {\hat{Q}}_{i + 1 / 2} + {\hat{Q}}_{i + 3 / 2} = \frac{1}{4} {\hat{Q}}_{i} + \frac{3}{2} {\hat{Q}}_{i + 1} + \frac{1}{4} {\hat{Q}}_{i + 2} \end{array}|

式中: ${\hat{Q}}_{i}$ 为节点处的数值通量.

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1 WCNS7-Z格式的模板构造

Fig.1 Candidate stencils of WCNS7-Z scheme

设定子模板的线性权值为

(6)d₀=

\frac{1}{56}

, d₁=

\frac{23}{56}

, d₂=

\frac{25}{56}

, d₃=

\frac{7}{56}

组合可得到7阶紧致插值格式:

(7)

\frac{1}{2} {\hat{Q}}_{i - 1 / 2}

{\hat{Q}}_{i + 1 / 2}

\frac{3}{14} {\hat{Q}}_{i + 3 / 2}

\frac{1}{224}

Q_i_-2+

\frac{1}{8}

Q_i_-1+

\frac{15}{16}

Q_i+

\frac{5}{8}

Q_i₊₁+

\frac{1}{32}

Q_i₊₂

经过非线性权w₀、w₁、w₂、w₃再组合后可得:

(8)\begin{gathered}\left|\begin{array}{c|c}5w_0+w_1&\hat{Q}_{i-1/2}\\w_0+w_1+w_2+w_3&\hat{Q}_{i+1/2}\\w_2/5+w_3&\hat{Q}_{i+3/2}\end{array}\right|= \\\left|\begin{array}{r|l}-w_0/4&\\(50w_0+5w_1-w_2)/20&Q_{i-1}\\(15w_0+6w_1+2zw_2+w_3)/4&Q_i\\(w_1+3w_2+6w_3)/4&Q_{i+1}\\w_3/4&Q_{i+2}\end{array}\right| \end{gathered}

对于非线性权w_k的使用策略,Jiang等^[5]设计的经典权函数方案对去除跨越间断模板带来的振荡现象有很好的效果,但也引入了较大的非线性误差,且在极值点附近格式的精度不能得到保证,Borges等^[9]引入更高阶的光滑指标,并因此提出了计算成本更低、精度更高的Z型权函数,显著地降低了由非线性权带来的格式耗散.权重方案采用形式简洁的Z型权函数:

(9)w_k=

\frac{α_{k}}{\overset{3}{\sum_{l = 0}} α_{l}}

, α _k=d_k

(1 + \frac{τ_{7}}{β_{k} + ε})

式中:k=0,1,2,3;ε=10^-12为防止分母出现零除引入的小量;τ₇为全局模板光滑度指标,

(10)τ₇=

|β_{0} + 3 β_{1} - 3 β_{2} - β_{3}|

=O(h⁷)

式中:O(h⁷)是数学上表示“与h⁷同阶的项”的一种记法,并且可明显地表示出截断误差的量级.

非线性权重中局部平滑度指标β_k的大小决定了每一个子模板的光滑度,通式为

(11)β_k=

c_{1}^{2}

c_{2}^{2}

c_{3}^{2}

对于k=0的子模板:

(12)

\begin{array}{l} c_{1} = (- 2 Q_{i - 3} + 9 Q_{i - 2} - 18 Q_{i - 1} + 11 Q_{i}) / 6 \\ c_{2} = - Q_{i - 3} + 4 Q_{i - 2} - 5 Q_{i - 1} + 2 Q_{i} \\ c_{3} = - Q_{i - 3} + 3 Q_{i - 2} - 3 Q_{i - 1} + Q_{i} \end{array}|

对于k=1的子模板:

(13)

\begin{array}{l} c_{1} = (Q_{i - 2} - 6 Q_{i - 1} + 3 Q_{i} + 2 Q_{i + 1}) / 6 \\ c_{2} = Q_{i - 1} - 2 Q_{i} + Q_{i + 1} \\ c_{3} = - Q_{i - 2} + 3 Q_{i - 1} - 3 Q_{i} + Q_{i + 1} \end{array}|

对于k=2的子模板:

(14)

\begin{array}{l} c_{1} = (- 2 Q_{i - 1} - 3 Q_{i} + 6 Q_{i - 1} + 11 Q_{i}) / 6 \\ c_{2} = Q_{i - 1} - 2 Q_{i} + Q_{i + 1} \\ c_{3} = - Q_{i - 1} + 3 Q_{i} - 3 Q_{i + 1} + Q_{i + 2} \end{array}|

对于k=3的子模板:

(15)

\begin{array}{l} c_{1} = (- 11 Q_{i} + 18 Q_{i + 1} - 9 Q_{i + 2} + 2 Q_{i + 3}) / 6 \\ c_{2} = 2 Q_{i} - 5 Q_{i + 1} + 4 Q_{i + 2} - Q_{i + 3} \\ c_{3} = - Q_{i} + 3 Q_{i + 1} - 3 Q_{i + 2} + Q_{i + 3} \end{array}|

1.3 WCNS7-T格式构造

不同于传统的权函数方案策略,TENO格式的模板加权策略在优化WENO格式模板选择的同时,更为准确地评估了模板光滑度,起到了间断检测器的作用,并在重构中以完全去除振荡模板来取代WENO格式通过非线性加权策略减小振荡模板影响的方式,使格式具有更加健壮的间断捕捉能力.在高阶的TENO方案中,如Hiejima^[23]所采取的策略一样,使用了递增宽度的模板序列.图2展示了WCNS7-T格式的模板构造,为了获得7阶空间精度,使用了3个3点模板和2个4点模板,从而使模板的选择更加紧致.

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2 WCNS7-T格式的模板构造

Fig.2 Candidate stencils of WCNS7-T scheme

传统WCNS中插值过程计算量大且复杂,对光滑区域来说这种插值没有必要.利用评估策略发现间断点,对于格式计算效率和精度都有非常重要的意义.将TENO格式的模板加权和宽度递增思想引进WCNS构造,在构造间断检测器的同时,更快地实现线性格式和非线性格式之间的重建切换.因此,在WCNS7-T格式中将表征插值模板光滑度的权重分配方案更替为TENO的光滑度评测策略,有利于增强格式的激波捕捉能力.优化后WCNS7-T格式的权重方案为

(16)

\begin{array}{l} γ_{k} = {(\frac{1}{β_{k, r} + ε})}^{q}, χ_{k} = γ_{k} / \overset{K - 3}{\sum_{k = 0}} γ_{k} \\ δ_{k} = \{\begin{array}{l} 0, & χ_{k} < C_{T} \\ 1, & 其 他 \end{array} \\ w_{k} = d_{k} δ_{k} / \overset{K - 3}{\sum_{k = 0}} d_{k} δ_{k} \end{array}|

式中:k=0, 1, …, K-3;γ_k、χ_k、δ_k、q皆为光滑因子指示参数,经截断函数C_T过滤得到最终非线性权值w_k.

光滑因子计算式与WCNS的不同,具体形式为

β₀=

\frac{13}{12}

(Q_i_-2-2Q_i_-1+Q_i)²+

\frac{1}{4}

(Q_i_-2-4Q_i_-1+3Q_i)²

\begin{array}{l} β_{1} = & \frac{13}{12} (Q_{i - 1} - 2 Q_{i} + Q_{i + 1})^{2} + \frac{1}{4} (Q_{i - 1} - Q_{i + 1})^{2} \\ β_{2} = & \frac{13}{12} (Q_{i} - 2 Q_{i + 1} + Q_{i + 2})^{2} + \\ \frac{1}{4} (3 Q_{i} - 4 Q_{i + 1} + Q_{i + 2})^{2} \\ β_{3} = & \frac{1}{36} (- 11 Q_{i} + 18 Q_{i + 1} - 9 Q_{i + 2} + 2 Q_{i + 3})^{2} + \\ \frac{13}{12} (2 Q_{i} - 5 Q_{i + 1} + 4 Q_{i + 2} - Q_{i + 3})^{2} + \\ \frac{781}{720} (- Q_{i} + 3 Q_{i + 1} - 3 Q_{i + 2} + Q_{i + 3})^{2} \\ β_{4} = & \frac{1}{36} (- 2 Q_{i - 3} + 9 Q_{i - 2} - 18 Q_{i - 1} + 11 Q_{i})^{2} + \\ \frac{13}{12} (- Q_{i - 3} + 4 Q_{i - 2} - 5 Q_{i - 1} + 2 Q_{i})^{2} + \\ \frac{781}{720} (- Q_{i - 3} + 3 Q_{i - 2} - 3 Q_{i - 1} + Q_{i})^{2} \end{array}

对于非光滑的模板,根据类ENO模板的选择,δ_k=0,这样振荡模板将会被剔除;而对于光滑模板,δ_k=1,格式将恢复为线性方案.显然,截断参数C_T的取值至关重要,而控制耗散特性的C_T取值依赖于不同的问题,且在实际应用中,参数调整通常非常耗时.Fu等^[7]给出了C_T的自适应函数取值公式.对于本文中的算例,控制耗散特性的C_T值在几组实验下对比测算后,选定为10^-7,其中K阶精度 ${τ_{K}}^{[7]}$ 对应各个阶数的计算量过大,尤其在高阶时计算成本过高,本文将其值依旧限定为1,ε=10^-12,其余光滑度指标参数与TENO给定的值保持一致.

2 数值测试算例

为了验证WCNS7-T格式计算表现和性能,开展了一维和二维问题的算例验证,计算结果与传统的WCNS7-Z格式进行对比.求解过程中采用了HLLC通量分裂方法^[24],时间推进采用3阶TVD Runge-Kutta格式^[25],实验的参考精确解由WCNS5-JS通过网格加密获得.

2.1 数值精度测试

为验证WCNS7-T格式的计算性能,首先以文献[8]中的数值精度实验为算例,给定初始条件为

(17)

\begin{array}{l} ρ (x, y, 0) = 1 + 0.2 s i n (x + y) \\ u (x, y, 0) = 1, v (x, y, 0) = 1 \\ p (x, y, 0) = 1 \end{array}\}

计算域[x, y]同样设定为 $[0, π]$ × $[0, π]$ ,u、v分别为x、y方向上的速度分量.各个方向上使用周期边界条件,精确的密度解为ρ(x, y, t)=1+0.2sin(x+y-2t),计算时间直到t=2 s,表1中给出了WCNS7-Z和WCNS7-T格式L₁和L_∞的误差(E₁和E₂)和精度(P₁和P₂)比较,网格单元划分N_x从10²到320².可知,相较于WCNS7-Z格式,WCNS7-T格式体现出更好的网格收敛性,且格式的设计精度也得到了较好的数值验证.

表1 二维欧拉方程的WCNS7-Z 和 WCNS7-T 格式的误差和精度

Tab.1 Errors and accuracy with two dimensional Euler equations for WCNS7-Z and WCNS7-T schemes

N_x	WCNS7-Z				WCNS7-T
N_x	E₁	P₁	E₂	P₂	E₁	P₁	E₂	P₂
10²	4.34×10^-3	0.00	3.12×10^-3	0.00	2.17×10^-3	0.00	1.47×10^-3	0.00
20²	5.74×10^-5	6.24	3.91×10^-5	6.32	1.86×10^-5	6.87	1.02×10^-5	7.17
40²	5.41×10^-7	6.73	3.46×10^-7	6.82	1.33×10^-7	7.13	8.08×10^-8	6.98
80²	4.17×10^-9	7.02	2.51×10^-9	7.11	8.15×10^-10	7.35	5.46×10^-10	7.21
160²	3.34×10^-11	6.96	1.86×10^-11	7.08	5.98×10^-12	7.09	3.90×10^-12	7.13
320²	2.19×10^-13	7.25	1.43×10^-13	7.02	4.04×10^-14	7.21	3.01×10^-14	7.02

新窗口打开| 下载CSV

2.2 Sod问题

一维Sod激波管问题是检验差分格式性质的经典算例之一.该问题包含一道激波、一个接触间断和一个膨胀波,问题的具体描述参考文献[26].给定初始条件为

(ρ, u, p)=

\{\begin{array}{l} (1, 0, 1), & - 0.5 \leq x < 0 \\ (0.125, 0, 1), & 0 \leq x \leq 0.5 \end{array}

计算域为[-0.5, 0.5],网格均匀划分为101个网格点,其中Δx=0.01.在边界处设置零梯度边界条件,以收敛条件判断数CFL为0.4和定时间步长Δt=0.01 s,计算到t=0.2 s.计算结果如图3所示,WCNS7-Z格式和WCNS7-T格式都能很好地捕捉间断,不会产生虚假的数值振荡,表明在间断处两种格式都引入了足够且稳定的耗散.然而,WCNS7-T格式要比WCNS7-Z格式在激波和间断处斜率更陡峭,对激波间断位置的分辨率更准确.计算结果表明,WCNS7-T格式引入的数值耗散相对较小,对激波捕获的精度更高.

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3 Sod激波管问题的的密度分布结果

Fig.3 Density distribution results for the Sod problem

2.3 Lax问题

Lax问题是左稀疏波、右激波类型的黎曼问题,同样能够用于检验数值算法捕捉间断能力.设定初始条件为

(ρ, u, p)=

\{\begin{array}{l} (0.445, 0.698, 3.528), & 0 \leq x < 0.5 \\ (0.500, 0.000, 0.571), & 0.5 \leq x \leq 1 \end{array}

计算域为[0,1],均匀划分为101个网格点,其中Δx = 0.01,CFL为0.4,Δt=0.01 s,计算到t=0.14 s,密度分布计算结果如图4所示.由放大视图可见,数值解未观察到过冲或振荡问题,与精确解较为吻合.WCNS7-T通过加权组合多个低阶紧致格式达到了高阶紧致格式,在间断附近采用TENO格式完全剔除振荡模板的方法,获得了比WCNS7-Z更加陡峭的间断捕获,表明WCNS7-T格式具有更弱的局部耗散特性使其能够更高效地捕捉间断.

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4 Lax问题的密度分布结果

Fig.4 Density distribution results for the Lax problem

2.4 Shu-Osher问题

Shu-Osher问题涉及向右移动的超音速激波与正弦分布密度场相互作用, 构造一个既有平滑结构又有不连续结构的流场,用于测试格式捕获流场间断和解析平滑多尺度结构的能力.设定初始条件为

(ρ, u, p)=

\{\begin{array}{l} (\frac{27}{7}, \frac{4 \sqrt{35}}{9}, \frac{31}{3}), & - 5 \leq x < - 4 \\ (1 + \frac{1}{5} s i n 5 x, 0, 1), & - 4 \leq x \leq 5 \end{array}

接触间断由两个强激波的碰撞产生,以 4 001 网格点计算结果为参考,终止时间为t=1.8 s,两种格式计算网格取为400,CFL数设为0.5,测试结果如图5所示.从WCNS7-Z和WCNS7-T格式的计算对比结果可见,WCNS7-T格式的分辨率更高,对高频波动的求解有着更好的波形和幅度,对可分辨频域的波峰波谷区域抹平也更少,表明WCNS7-T格式比经典的WCNS7-Z格式具有更优的间断捕捉能力和脉动分辨率.

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5 Shu-Osher问题的密度分布结果

Fig.5 Density distribution results for the Shu-Osher problem

2.5 Titarev-Toro问题

Titarev-Toro问题是对Shu-Osher问题的强扩展,初始条件给定为

(ρ, u, p)=

\{\begin{array}{l} (1.515 695, 0.523 326, 1.805), \\ - 5 \leq x < - 4.5 \\ (1 + \frac{1}{10} s i n 20 π x, 0, 1), \\ - 4.5 \leq x \leq 5 \end{array}

计算域为[-5, 5],均分为1 001个网格点,其中Δx=0.01,边界设置为零梯度边界条件.给定时间步长Δt=0.5 s,CFL数设为0.5,终止时间为t=5 s.由于没有理论上的精确解,所以参考结果由WCNS5-JS格式在 5 001 个网格点计算获得.图6给出了两种格式计算的密度分布, 与Shu-Osher问题测试算例相似,WCNS7-T格式在激波后解析出了更好的密度波形,而WCNS7-Z格式结果中小尺度波结构较为模糊,数值耗散也较强,分辨率也较低,表明WCNS7-T格式比WCNS7-Z格式能更优、更好地解析局部大振幅密度波动问题.

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6 Titarev-Toro问题的密度分布结果

Fig.6 Density distribution results for the Titarev-Toro problem

2.6 二维黎曼问题

Lax等^[27]使用的19个二维黎曼问题中的第3个问题和第6个问题所产生的流场都具有丰富的小尺度特征,非常适合用于测验数值格式求解精细流场结构的能力,因此选择其中的这两个问题来评估目前WCNS7-T格式的性能.

2.6.1 算例 3

初始条件给定为

(ρ, u, v, p)=

\{\begin{array}{l} (1.5, 0, 0, 1.5), \\ x > 0.8, y > 0.8 \\ (33 / 62, 4 \sqrt{11}, 0, 0.3), \\ x \leq 0.8, y > 0.8 \\ (77 / 558, 4 \sqrt{11}, 4 \sqrt{11}, 9 / 310), \\ x \leq 0.8, y \leq 0.8 \\ (33 / 62, 0, 4 \sqrt{11}, 0.3), \\ x > 0.8, y \leq 0.8 \end{array}

设定计算域为[0,1]×[0,1],网格尺寸划分为 1 001×1 001,终止时间为t=0.3 s,边界条件都设定为无反射边界条件.Kelvin-Helmholtz不稳定性在滑移线上产生的小尺度复杂结构情况可用来评估数值格式的耗散表现.图7所示为黎曼问题3的密度轮廓结果.从图7可以看出,对比于WCNS7-Z格式,WCNS7-T格式在滑移线上捕获的涡量更多,能够更好地捕捉到激波间断导致的流动不稳定性.此外,在这种数值情况下,WCNS7-T的结果与文献[28 ⇓-30]中一样在引用TENO非线性加权后出现了对称性不足的现象,本文推测可能是控制格式耗散表现的C_T值的取值问题所致.

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7 黎曼问题3密度轮廓结果

Fig.7 Density contour results for Riemann problem 3

2.6.2 算例6

第2个黎曼的计算域、网格划分和边界条件不变,终止时间设定为t=0.3 s.初始条件为

(ρ, u, v, p)=

\{\begin{array}{l} (1, 0.75, - 0.5, 1), & (x, y) \in [\frac{1}{2}, 1] \times [\frac{1}{2}, 1] \\ (2, 0.75, 0.5, 1), & (x, y) \in [0, \frac{1}{2}] \times [\frac{1}{2}, 1] \\ (1, - 0.75, 0.5, 1), & (x, y) \in [0, \frac{1}{2}] \times [0, \frac{1}{2}] \\ (3, - 0.75, - 0.5, 1), & (x, y) \in [\frac{1}{2}, 1] \times [0, \frac{1}{2}] \end{array}

图8所示为黎曼问题6的密度轮廓结果.由图可见,两种格式都解析出了接触线上丰富的小尺度流动结构,而WCNS7-T格式的结果有更丰富的小尺度结构,流场分辨率更高,产生的耗散也小于传统的WCNS7-Z格式.表明WCNS7-T格式经过TENO非线性加权后数值耗散相比之前有所减小,且模板选择策略的改变也进一步降低了整体的数值耗散,使格式在流场计算时能更准确模拟流场细节结构.

图8

新窗口打开| 下载原图ZIP| 生成PPT

图8 黎曼问题6密度轮廓结果

Fig.8 Density contour results for Riemann problem 6

3 结论

基于TENO格式的模板选择和加权策略,将其间断检测和光滑度因子引入WCNS构造中,发展了一种7阶精度的WCNS7-T格式,在抑制数值耗散的同时增强了激波捕捉能力,通过一维和二维问题开展算例测试,验证了WCNS7-T格式的良好性能,重要结论总结如下:

(1) WCNS7-T格式体现出了流场光滑区域耗散小、激波间断区域分辨率高、计算结果更准确的特点,表明经权重优化后格式不仅保持了原WCNS的稳定,也强化了小扰动分辨率和激波捕获能力.

(2) WCNS7-T格式以完全剔除振荡模板来取代传统的加权分配,而在光滑区域恢复最优非线性插值,不仅降低了格式在间断处需要更复杂的平滑指标所带来的较高计算量,保证了格式的数值精度,也提高了格式的间断检测能力,增强了加权方法抑制数值振荡的能力.

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

HARTEN

High resolution schemes for hyperbolic conservation laws

[J]. Journal of Computational Physics, 1997, 135(2): 260-278.

DOI:10.1006/jcph.1997.5713 URL [本文引用: 1]

[2]

LIOU

M S

A sequel to AUSM: AUSM⁺

[J]. Journal of Computational Physics, 1996, 129(2): 364-382.

DOI:10.1006/jcph.1996.0256 URL [本文引用: 1]

[3]

HARTEN

, ENGQUIST

, OSHER

et al.

Uniformly high order accurate essentially non-oscillatory schemes, III

[J]. Journal of Computational Physics, 1987, 71(1): 231-303.

DOI:10.1016/0021-9991(87)90031-3 URL [本文引用: 1]

[4]

SHU

C W

, OSHER

Efficient implementation of essentially non-oscillatory shock-capturing schemes

[J]. Journal of Computational Physics, 1989, 77(2): 439-471.

DOI:10.1016/0021-9991(88)90177-5 URL [本文引用: 1]

[5]

JIANG

G S

, SHU

C W

Efficient implementation of weighted ENO schemes

[J]. Journal of Computational Physics, 1996, 126(1): 202-228.

DOI:10.1006/jcph.1996.0130 URL [本文引用: 4]

[6]

DENG

, ZHANG

Developing high-order weighted compact nonlinear schemes

[J]. Journal of Computational Physics, 2000, 165(1): 22-44.

DOI:10.1006/jcph.2000.6594 URL [本文引用: 2]

[7]

, HU

X Y

, ADAMS

N A

A new class of adaptive high-order targeted ENO schemes for hyperbolic conservation laws

[J]. Journal of Computational Physics, 2018, 374: 724-751.

DOI:10.1016/j.jcp.2018.07.043 URL [本文引用: 3]

[8]

C C

, WAN

Z H

, SUN

D J

An alternative formulation of targeted ENO scheme for hyperbolic conservation laws

[J]. Computers & Fluids, 2022, 238: 105368.

DOI:10.1016/j.compfluid.2022.105368 URL [本文引用: 3]

[9]

BORGES

, CARMONA

, COSTA

et al.

An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws

[J]. Journal of Computational Physics, 2008, 227(6): 3191-3211.

DOI:10.1016/j.jcp.2007.11.038 URL [本文引用: 2]

[10]

MARTIN

M P

, TAYLOR

E M

, WU

et al.

A bandwidth-optimized WENO scheme for the effective direct numerical simulation of compressible turbulence

[J]. Journal of Computational Physics, 2006, 220(1): 270-289.

[本文引用: 1]

[11]

X Y

, WANG

, ADAMS

N A

An adaptive central-upwind weighted essentially non-oscillatory scheme

[J]. Journal of Computational Physics, 2010, 229(23): 8952-8965.

DOI:10.1016/j.jcp.2010.08.019 URL [本文引用: 1]

[12]

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.

DOI:10.1006/jcph.2000.6443 URL [本文引用: 1]

[13]

GEROLYMOS G A

. SÉNÉCHAL

, VALLET

Very-high-order WENO schemes

[J]. Journal of Computational Physics, 2009, 228(23): 8481-8524.

DOI:10.1016/j.jcp.2009.07.039 URL [本文引用: 1]

[14]

PIROZZOLI

On the spectral properties of shock-capturing schemes

[J]. Journal of Computational Physics, 2006, 219(2): 489-497.

DOI:10.1016/j.jcp.2006.07.009 URL [本文引用: 1]

[15]

JOHNSEN

, LARSSON

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

DOI:10.1016/j.jcp.2009.10.028 URL [本文引用: 1]

[16]

DENG

, MAEKAWA

Compact high-order accurate nonlinear schemes

[J]. Journal of Computational Physics, 1997, 130(1): 77-91.

DOI:10.1006/jcph.1996.5553 URL [本文引用: 1]

[17]

PENG

, LIU

, LI

et al.

An efficient targeted ENO scheme with local adaptive dissipation for compressible flow simulation

[J]. Journal of Computational Physics, 2021, 425: 109902.

DOI:10.1016/j.jcp.2020.109902 URL [本文引用: 1]

[18]

ZHANG

, ZHANG

, XU

Towards optimal high-order compact schemes for simulating compressible flows

[J]. Applied Mathematics and Computation, 2019, 355: 221-237.

DOI:10.1016/j.amc.2019.03.001 URL [本文引用: 1]

[19]

马燕凯, 刘化勇, 燕振国,

等.

基于 HWCNS 格式的紧致插值方法研究

[J]. 计算力学学报, 2015, 32(3): 388-393.

[本文引用: 2]

Yankai

, LIU

Huayong

, YAN

Zhenguo

et al.

Research on compact interpolation method based on HWCNS scheme

[J]. Chinese Journal of Computational Mechanics, 2015, 32(3): 388-393.

[本文引用: 2]

[20]

ROE

P L

Approximate Riemann solvers, parameter vectors, and difference schemes

[J]. Journal of Computational Physics, 1997, 135(2): 250-258.

DOI:10.1006/jcph.1997.5705 URL [本文引用: 1]

[21]

LIOU

M S

, STEFFEN

C J

A new flux splitting scheme

[J]. Journal of Computational Physics, 1993, 107(1): 23-39.

DOI:10.1006/jcph.1993.1122 URL [本文引用: 1]

[22]

VAN

LEER B

Flux-vector splitting for the Euler equations

[J]. Lecture Notes in Physics, 1982, 170(1): 507-512.

[本文引用: 1]

[23]

HIEJIMA

A high-order weighted compact nonlinear scheme for compressible flows

[J]. Computers & Fluids, 2022, 232: 105199.

DOI:10.1016/j.compfluid.2021.105199 URL [本文引用: 1]

[24]

, SPRUCE

, SPEARES

Restoration of the contact surface in the HLL-Riemann solver

[J]. Shock Waves, 1994, 4(1): 25-34.

DOI:10.1007/BF01414629 URL [本文引用: 1]

[25]

GOTTLIEB

, SHU

C W

Total variation diminishing Runge-Kutta schemes

[J]. Mathematics of Computation, 1998, 67(221): 73-85.

DOI:10.1090/mcom/1998-67-221 URL [本文引用: 1]

[26]

WONGA

M L

, LELEA

S K

High-order localized dissipation weighted compact nonlinear scheme for shock-and interface-capturing in compressible flows

[J]. Journal of Computational Physics, 2017, 339: 179-209.

DOI:10.1016/j.jcp.2017.03.008 URL [本文引用: 1]

[27]

LAX

P D

, LIU

X D

Solution of two-dimensional Riemann problems of gas dynamics by positive schemes

[J]. Siam Journal on Scientific Computing, 1998, 19(2): 319-340.

DOI:10.1137/S1064827595291819 URL [本文引用: 1]

[28]

ZHANG

, ZHANG

, LIU

et al.

A simple extended compact nonlinear scheme with adaptive dissipation control

[J]. Communications in Nonlinear Science and Numerical Simulation, 2020, 84: 105191.

DOI:10.1016/j.cnsns.2020.105191 URL [本文引用: 1]

[29]

SHI

, ZHANG

Y T

, SHU

C W

Resolution of high order WENO schemes for complicated flow structures

[J]. Journal of Computational Physics, 2003, 186(2): 690-696.

DOI:10.1016/S0021-9991(03)00094-9 URL [本文引用: 1]

[30]

ZHAO

, SUN

, XIE

et al.

Numerical dissipation control in an adaptive WCNS with a new smoothness indicator

[J]. Applied Mathematics and Computation, 2018, 330: 239-253.

DOI:10.1016/j.amc.2018.01.019 URL [本文引用: 1]

High resolution schemes for hyperbolic conservation laws

1997

... 可压缩流动广泛存在于超声速/高超声速飞行器中,激波间断和湍流共存且相互耦合作用给高可靠性数值模拟带来了巨大挑战.近40年,人们发展了众多高精度激波捕获方法,例如总变差不增(TVD)方法^[1]、ASUM^[2]、基本无振荡/加权基本无振荡(ENO/WENO)方法^[3⇓-5]、加权非线性紧致格式(WCNS)方法^[6]、基于目标本质无振荡(TENO)方法^[7]和AFTENO^[8]等.这些格式在激波间断附近增加局部耗散来抑制非物理数值振荡,而在光滑湍流区域减少耗散以保持对多尺度湍流结构的分辨率.经典激波捕捉格式,如TVD方法和ENO方法,虽然可以解决激波间断无数值振荡,但对湍流光滑区域产生了过多的数值耗散.WENO方法在ENO方法的基础上,引入子模板光滑因子,通过非线性加权方法来动态调整局部耗散^[5],使得高精度数值格式得到进一步发展.Borges等^[9]在非线性权值中引入更高阶的光滑指标,提出计算成本更低、精度更高的WENO-Z格式.虽然WENO格式的精度可以设计为任意高阶^{[10⇓⇓-13]},但WENO格式对光谱区域的分辨率仍然不能令人满意^[14],它的加权机制和光滑度指标设定依旧对小尺度湍流结构中引入过多的数值耗散^[15].1999年,Deng等^[16]提出非线性紧致格式(CNS),后续引入加权技术,构造了 WCNS^[6].近年来,Fu等^[7]提出了另一类非线性加权框架,即TENO 方法,通过类ENO方法的模板选择来控制非线性耗散,在无间断区域(包括光滑临界点)保证最优的高阶空间重构,在间断区域通过移除跨越间断的模板来避免数值振荡,这种方法有利于候选模板的使用和高阶格式的改进.当子模板中存在多个相邻的不连续间断,该方法可退化为三阶重构.即使对于高阶版本,数值鲁棒性也得到了明显改善,且避免了阶跃退化问题.Peng等^[17]在此基础上进一步改进了TENO格式的C_T自适应函数策略,发展了TENO5-LAD格式,实验结果证明新方法能够更好地抑制不连续点附近的数值振荡,并能以较低的额外计算成本保持TENO的低耗散特性.而Ye等^[8]结合AFWENO格式和TENO格式的优点构建了可以实现任意高阶精度的AFTENO格式,在激波捕获和波数分辨率方面,相较于AFWENO-JS和AFWENO-Z有着明显的优势.之后,在加权紧致格式方面,Zhang等^[18]将WCNS与TENO格式组合发展了TCNS,实验结果显示TCNS能够在模拟中捕获更丰富的波结构,特别是在对瑞利-泰勒不稳定性问题的模拟中,TCNS在较粗网格上取得了与WCNS-JS^[5]相似甚至更好的结果. ...

A sequel to AUSM: AUSM⁺

1996

Uniformly high order accurate essentially non-oscillatory schemes, III

1987

Efficient implementation of essentially non-oscillatory shock-capturing schemes

1989

Efficient implementation of weighted ENO schemes

1996

... [5],使得高精度数值格式得到进一步发展.Borges等^[9]在非线性权值中引入更高阶的光滑指标,提出计算成本更低、精度更高的WENO-Z格式.虽然WENO格式的精度可以设计为任意高阶^{[10⇓⇓-13]},但WENO格式对光谱区域的分辨率仍然不能令人满意^[14],它的加权机制和光滑度指标设定依旧对小尺度湍流结构中引入过多的数值耗散^[15].1999年,Deng等^[16]提出非线性紧致格式(CNS),后续引入加权技术,构造了 WCNS^[6].近年来,Fu等^[7]提出了另一类非线性加权框架,即TENO 方法,通过类ENO方法的模板选择来控制非线性耗散,在无间断区域(包括光滑临界点)保证最优的高阶空间重构,在间断区域通过移除跨越间断的模板来避免数值振荡,这种方法有利于候选模板的使用和高阶格式的改进.当子模板中存在多个相邻的不连续间断,该方法可退化为三阶重构.即使对于高阶版本,数值鲁棒性也得到了明显改善,且避免了阶跃退化问题.Peng等^[17]在此基础上进一步改进了TENO格式的C_T自适应函数策略,发展了TENO5-LAD格式,实验结果证明新方法能够更好地抑制不连续点附近的数值振荡,并能以较低的额外计算成本保持TENO的低耗散特性.而Ye等^[8]结合AFWENO格式和TENO格式的优点构建了可以实现任意高阶精度的AFTENO格式,在激波捕获和波数分辨率方面,相较于AFWENO-JS和AFWENO-Z有着明显的优势.之后,在加权紧致格式方面,Zhang等^[18]将WCNS与TENO格式组合发展了TCNS,实验结果显示TCNS能够在模拟中捕获更丰富的波结构,特别是在对瑞利-泰勒不稳定性问题的模拟中,TCNS在较粗网格上取得了与WCNS-JS^[5]相似甚至更好的结果. ...

... [5]相似甚至更好的结果. ...

... 对于非线性权w_k的使用策略,Jiang等^[5]设计的经典权函数方案对去除跨越间断模板带来的振荡现象有很好的效果,但也引入了较大的非线性误差,且在极值点附近格式的精度不能得到保证,Borges等^[9]引入更高阶的光滑指标,并因此提出了计算成本更低、精度更高的Z型权函数,显著地降低了由非线性权带来的格式耗散.权重方案采用形式简洁的Z型权函数: ...

Developing high-order weighted compact nonlinear schemes

2000

... [6].近年来,Fu等^[7]提出了另一类非线性加权框架,即TENO 方法,通过类ENO方法的模板选择来控制非线性耗散,在无间断区域(包括光滑临界点)保证最优的高阶空间重构,在间断区域通过移除跨越间断的模板来避免数值振荡,这种方法有利于候选模板的使用和高阶格式的改进.当子模板中存在多个相邻的不连续间断,该方法可退化为三阶重构.即使对于高阶版本,数值鲁棒性也得到了明显改善,且避免了阶跃退化问题.Peng等^[17]在此基础上进一步改进了TENO格式的C_T自适应函数策略,发展了TENO5-LAD格式,实验结果证明新方法能够更好地抑制不连续点附近的数值振荡,并能以较低的额外计算成本保持TENO的低耗散特性.而Ye等^[8]结合AFWENO格式和TENO格式的优点构建了可以实现任意高阶精度的AFTENO格式,在激波捕获和波数分辨率方面,相较于AFWENO-JS和AFWENO-Z有着明显的优势.之后,在加权紧致格式方面,Zhang等^[18]将WCNS与TENO格式组合发展了TCNS,实验结果显示TCNS能够在模拟中捕获更丰富的波结构,特别是在对瑞利-泰勒不稳定性问题的模拟中,TCNS在较粗网格上取得了与WCNS-JS^[5]相似甚至更好的结果. ...

A new class of adaptive high-order targeted ENO schemes for hyperbolic conservation laws

2018

... [7]提出了另一类非线性加权框架,即TENO 方法,通过类ENO方法的模板选择来控制非线性耗散,在无间断区域(包括光滑临界点)保证最优的高阶空间重构,在间断区域通过移除跨越间断的模板来避免数值振荡,这种方法有利于候选模板的使用和高阶格式的改进.当子模板中存在多个相邻的不连续间断,该方法可退化为三阶重构.即使对于高阶版本,数值鲁棒性也得到了明显改善,且避免了阶跃退化问题.Peng等^[17]在此基础上进一步改进了TENO格式的C_T自适应函数策略,发展了TENO5-LAD格式,实验结果证明新方法能够更好地抑制不连续点附近的数值振荡,并能以较低的额外计算成本保持TENO的低耗散特性.而Ye等^[8]结合AFWENO格式和TENO格式的优点构建了可以实现任意高阶精度的AFTENO格式,在激波捕获和波数分辨率方面,相较于AFWENO-JS和AFWENO-Z有着明显的优势.之后,在加权紧致格式方面,Zhang等^[18]将WCNS与TENO格式组合发展了TCNS,实验结果显示TCNS能够在模拟中捕获更丰富的波结构,特别是在对瑞利-泰勒不稳定性问题的模拟中,TCNS在较粗网格上取得了与WCNS-JS^[5]相似甚至更好的结果. ...

... 对于非光滑的模板,根据类ENO模板的选择,δ_k=0,这样振荡模板将会被剔除;而对于光滑模板,δ_k=1,格式将恢复为线性方案.显然,截断参数C_T的取值至关重要,而控制耗散特性的C_T取值依赖于不同的问题,且在实际应用中,参数调整通常非常耗时.Fu等^[7]给出了C_T的自适应函数取值公式.对于本文中的算例,控制耗散特性的C_T值在几组实验下对比测算后,选定为10^-7,其中K阶精度

{τ_{K}}^{[7]}

对应各个阶数的计算量过大,尤其在高阶时计算成本过高,本文将其值依旧限定为1,ε=10^-12,其余光滑度指标参数与TENO给定的值保持一致. ...

An alternative formulation of targeted ENO scheme for hyperbolic conservation laws

2022

... [8]结合AFWENO格式和TENO格式的优点构建了可以实现任意高阶精度的AFTENO格式,在激波捕获和波数分辨率方面,相较于AFWENO-JS和AFWENO-Z有着明显的优势.之后,在加权紧致格式方面,Zhang等^[18]将WCNS与TENO格式组合发展了TCNS,实验结果显示TCNS能够在模拟中捕获更丰富的波结构,特别是在对瑞利-泰勒不稳定性问题的模拟中,TCNS在较粗网格上取得了与WCNS-JS^[5]相似甚至更好的结果. ...

... 为验证WCNS7-T格式的计算性能,首先以文献[8]中的数值精度实验为算例,给定初始条件为 ...

An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws

2008

A bandwidth-optimized WENO scheme for the effective direct numerical simulation of compressible turbulence

2006

An adaptive central-upwind weighted essentially non-oscillatory scheme

2010

Monotonicity preserving weighted essentially non-oscillatory schemes with increasingly high order of accuracy

2000

Very-high-order WENO schemes

2009

On the spectral properties of shock-capturing schemes

2006

Assessment of high-resolution methods for numerical simulations of compressible turbulence with shock waves

2010

Compact high-order accurate nonlinear schemes

1997

An efficient targeted ENO scheme with local adaptive dissipation for compressible flow simulation

2021

Towards optimal high-order compact schemes for simulating compressible flows

2019

基于 HWCNS 格式的紧致插值方法研究

2015

... 式中:E'_i为空间导数的近似值,采用显式的节点和半节点混合通量求解方法^[19], ...

... 与5阶格式的模板不同,在7阶WCNS^[19]的构造中使用4个4点子模板:(x_i_-3, x_i_-2, x_i_-1, x_i),(x_i_-2, x_i_-1, x_i, x_i₊₁),(x_i, x_i₊₁, x_i₊₂, x_i₊₃),(x_i, x_i₊₁, x_i₊₂, x_i₊₃),如图1所示.图中:L为全局模板;K₁、K₂、K₃、K₄分别是4个子模板.每个子模板上的4阶紧致格式设定为特定偏置形式: ...

基于 HWCNS 格式的紧致插值方法研究

2015

... 式中:E'_i为空间导数的近似值,采用显式的节点和半节点混合通量求解方法^[19], ...

Approximate Riemann solvers, parameter vectors, and difference schemes

1997

... 式中:α为常数;h为网格单元长度;

{\hat{E}}_{i}

为节点边界上的通量.残差主项为

(\frac{1}{8960} - \frac{5}{65536} α)

h⁶

E_{i}^{(6)}

,当α=175/256时格式具有8阶精度.对于半节点边界上的通量

{\hat{E}}_{i \pm 1 / 2}

=E(

{\hat{Q}}_{i \pm 1 / 2}

)可以使用多种通量差分或分裂方法求解,如Roe方法^[20]、AUSM^[21]、van leer方法^[22]等.

{\hat{Q}}_{i \pm 1 / 2}

是Q_i_±1_/₂的高阶近似,可由紧致插值格式或显式插值格式得到. ...

A new flux splitting scheme

1993

... 式中:α为常数;h为网格单元长度;

{\hat{E}}_{i}

为节点边界上的通量.残差主项为

(\frac{1}{8960} - \frac{5}{65536} α)

h⁶

E_{i}^{(6)}

,当α=175/256时格式具有8阶精度.对于半节点边界上的通量

{\hat{E}}_{i \pm 1 / 2}

=E(

{\hat{Q}}_{i \pm 1 / 2}

)可以使用多种通量差分或分裂方法求解,如Roe方法^[20]、AUSM^[21]、van leer方法^[22]等.

{\hat{Q}}_{i \pm 1 / 2}

是Q_i_±1_/₂的高阶近似,可由紧致插值格式或显式插值格式得到. ...

Flux-vector splitting for the Euler equations

1982

... 式中:α为常数;h为网格单元长度;

{\hat{E}}_{i}

为节点边界上的通量.残差主项为

(\frac{1}{8960} - \frac{5}{65536} α)

h⁶

E_{i}^{(6)}

,当α=175/256时格式具有8阶精度.对于半节点边界上的通量

{\hat{E}}_{i \pm 1 / 2}

=E(

{\hat{Q}}_{i \pm 1 / 2}

)可以使用多种通量差分或分裂方法求解,如Roe方法^[20]、AUSM^[21]、van leer方法^[22]等.

{\hat{Q}}_{i \pm 1 / 2}

是Q_i_±1_/₂的高阶近似,可由紧致插值格式或显式插值格式得到. ...

A high-order weighted compact nonlinear scheme for compressible flows

2022

... 不同于传统的权函数方案策略,TENO格式的模板加权策略在优化WENO格式模板选择的同时,更为准确地评估了模板光滑度,起到了间断检测器的作用,并在重构中以完全去除振荡模板来取代WENO格式通过非线性加权策略减小振荡模板影响的方式,使格式具有更加健壮的间断捕捉能力.在高阶的TENO方案中,如Hiejima^[23]所采取的策略一样,使用了递增宽度的模板序列.图2展示了WCNS7-T格式的模板构造,为了获得7阶空间精度,使用了3个3点模板和2个4点模板,从而使模板的选择更加紧致. ...

Restoration of the contact surface in the HLL-Riemann solver

1994

... 为了验证WCNS7-T格式计算表现和性能,开展了一维和二维问题的算例验证,计算结果与传统的WCNS7-Z格式进行对比.求解过程中采用了HLLC通量分裂方法^[24],时间推进采用3阶TVD Runge-Kutta格式^[25],实验的参考精确解由WCNS5-JS通过网格加密获得. ...

Total variation diminishing Runge-Kutta schemes

1998

High-order localized dissipation weighted compact nonlinear scheme for shock-and interface-capturing in compressible flows

2017

... 一维Sod激波管问题是检验差分格式性质的经典算例之一.该问题包含一道激波、一个接触间断和一个膨胀波,问题的具体描述参考文献[26].给定初始条件为 ...

Solution of two-dimensional Riemann problems of gas dynamics by positive schemes

1998

... Lax等^[27]使用的19个二维黎曼问题中的第3个问题和第6个问题所产生的流场都具有丰富的小尺度特征,非常适合用于测验数值格式求解精细流场结构的能力,因此选择其中的这两个问题来评估目前WCNS7-T格式的性能. ...

A simple extended compact nonlinear scheme with adaptive dissipation control

2020

... 设定计算域为[0,1]×[0,1],网格尺寸划分为 1 001×1 001,终止时间为t=0.3 s,边界条件都设定为无反射边界条件.Kelvin-Helmholtz不稳定性在滑移线上产生的小尺度复杂结构情况可用来评估数值格式的耗散表现.图7所示为黎曼问题3的密度轮廓结果.从图7可以看出,对比于WCNS7-Z格式,WCNS7-T格式在滑移线上捕获的涡量更多,能够更好地捕捉到激波间断导致的流动不稳定性.此外,在这种数值情况下,WCNS7-T的结果与文献[28⇓-30]中一样在引用TENO非线性加权后出现了对称性不足的现象,本文推测可能是控制格式耗散表现的C_T值的取值问题所致. ...

Resolution of high order WENO schemes for complicated flow structures

2003

Numerical dissipation control in an adaptive WCNS with a new smoothness indicator

2018

〈

〉