非线性优化的保单调差分格式.pdfVIP

LHD 2012 年度夏季学术研讨会 非线性优化的保单调差分格式 李新亮，冷岩 中国科学院力学研究所，高温气体动力学国家重点实验室（筹），北京海淀区 100190 摘要 采用保单调思想，结合非线性优化，构造了六阶精度的保单调差分格式(OMP6). 格式构造过程 中采用非线性Forier 分析对格式的色散及耗散误差进行了优化。通过 Shu-Osher,双 Mach 反射以及 RT 不稳定性等算例对该格式进行了验证。算例表明，与常用的 WENO 等格式相比，该格式(OMP6)具有 格式简洁、计算量小、分辨率高等优点。 关键词 保单调，优化，差分格式 表现同样有很大的影响。对于真实的物理流动， 引 言 在低波数区域（光滑区），保单调格式不调用限 高精度高分辨率数值方法在湍流的精细计 制器，整个格式的数值表现完全由线性部分决 算（DNS/LES ）中发挥着重要作用。对于可压 定。然而即使是在高波数区域（激波附近或者 缩湍流，流场中的激波将给计算带来很大困难。 剧烈的干扰区域），线性部分对整个格式的数值 为了捕捉激波，需引入一定的数值耗散以抑制 结果同样有很大的影响。因此，给出优化的线 数值振荡，而过大的数值耗散会耗散掉湍流的 性部分对于保单调格式同样重要。 小尺度脉动，从带来较大的误差。因而需要在 作者基于非线性谱分析方法对保单调格式 激波捕捉与分辨小尺度之间进行平衡。 进行了优化，构造了六阶保单调差分格式 近年来，以WENO 格式[1-2]为代表的高精度 (OMP6) 。 通过多个算例对格式进行了验证。 激波捕捉格式得到了广泛的应用。WENO 格式 1. 数值格式 具有很好的鲁棒性及数值分辨率，在可压缩湍 流的精细模拟(DNS 及 LES) 中应用尤为普遍。但 考虑标量双曲守恒律方程 WENO 格式仍有计算量大，耗散偏大等不足之 u f (u)  0 （1） 处，需要进行改进。目前针对 WENO 格式的优 t x 化与改进是高精度数值格式研究的热点 ，比较 其中：f (u ) au , a 0 有代表性的包括混合（hybrid ）WENO 格式以及 方程(3.1)的半离散形式为： 优化的 WENO 格式等。 u 1 与 WENO 格式的思想不同，Suresh 和  (f j 1/ 2 f j 1/ 2 ) 0 (2) t x   Huynh[3] 发 展 了 一 系 列 保 单 调 格 式 （ MP 根据保单调格式[ ]的思想，数值通量具有 schemes ），这些格式与常用的 WENO 格式相比， 如下表达式： 具有对小尺度波的高分辨率和计算效率高等优 f f Linear MP （3 ） 势。保单调格式可以看成由两部分构成，一是