高精度WENO有限体积格式.pptVIP

研究现状 Harten和Osher等人在1987年提出的ENO格式，重点论述了一维有限体积的构造，自此，人们对于高精度格式的研究进入了一个新的阶段。ENO格式已被成功的用来解决双曲型守恒律问题和其它的对流扩散问题。后来，Casper 和Harten把ENO有限体积格式的构造推广到二维守恒律方程。在间断解区域，ENO格式能有效的捕捉到激波，但此格式也有缺陷，就是在光滑解区域不很理想。Liu 等提出了WENO格式, 可以克服这种不足，这种方法不是选择最光滑的网格模板，而是取它们的一个凸组合，对组合中的每组网格模板上构造的插值多项式设置一个系数，以此决定对数值通量逼近函数的贡献，在光滑区域可以达到 阶精度（在ENO格式达到 阶精度的情况下）。 * * 上述进行的研究，大多局限于结构网格上。近年来，在非结构网格上开展这方面的工作也有了一些进展，如Abgrall，Harten对ENO格式进行了研究,Hu等人对WENO格式进行了研究。但总起来说开展这方面的研究工作不多，还不很成熟，没有形成全面系统的理论体系。 近年来又出现了用点云离散计算区域的无网格算法，适合处理复杂的气动外形。 区域达到 阶精度，在ENO格式达到 阶精度的情况下)，在间断解区域权设置系数近似于零，可达到基本无震荡的效果。 Jiang对WENO格式进行了改进和提高，提出了一种新的光滑措施，能使这种格式在光滑的 考虑一维守恒方程 写成半离散的形式 其中 其中 一维守恒方程 ENO和WENO格式 ⅰ ⅱ ⅲ 时，节点模板选择如下图 节点模板通常有固定模板和优化模板两种选取方法。 本文采用优化模板，在每次增加节点扩充模板时，要进行Newton差商的绝对值比较。取其极小值的扩充作为新的模板。设选择的节点模板为： 共有 组节点模板，每组有 个节点。例如对于 精度；权系数之和为1，即 光滑因子 是表示插值多项式函数光滑程度的度量，即表明插值多项式函数振荡程度如何。 越大，振荡越剧烈。 WENO格式的构造是建立在ENO的基础之上，它不是选择最光滑的节点模板构造插值多项式，而是对所有的多项式进行一次凸组合。振荡因子 的选择是关键。要求在间断 解附近，要无振荡（ENO）； 在光滑的区域，有尽可能高的 有限体积格式可以写成 本文选择Lax-Friedrichs通量 为Jacobian矩阵特征值的上确界。 为在控制体 上构造 插值多项式在 处的值， 为与相邻的控制体 上构造的插 处的值 。 ，由 上平均值 构造插值多项式函数。 插值多项式在 有限体积法计算的是平均值 ，在 网格 上 我们以三角形网格作为计算网格，三角形网格的重心作 为赋值点，即所谓的格子中心型。 控制方程的有限体积格式 构造二次多项式 网格模板的选取 本网格模板。 选取十个三角形网格作为基 构造WENO二次插值多项式选取的典型的网格模板如下图 * *