对于结构化网格上的3D流场,最精确的插值方法是什么?


9

我在3D结构化网格上求解多物种可压缩的Navier-Stokes方程。我已经在给定的网格上获得了一个解决方案(假设是一个相对粗糙的网格)。现在,我想在重新启动仿真之前优化网格并在新网格上插入以前的解决方案。当前,我们有一个插值工具,可以构建2个网格的kd树,然后可以使用2种不同的方法来计算新网格上的值:

  • 简单平均
  • 逆距离加权(IDW)
  • 移动最小二乘法(MLS)

我想专注于准确性,因为由于我处理的是大梯度,因此当我重新开始计算时,无法正确捕获它们会产生波动。我首先尝试了简单平均,但准确性不够好。

我认为带有2阶多项式的MLS方法会给我合理的结果,因为它应该是非振荡的。但是,当我查看插值字段时,我看到局部最小值/最大值超出了我的初始字段的值。这是否表示该程序中MLS的实现不正确?我应该小心模板的大小和多项式的顺序吗?您会推荐其他哪种方法?

提前致谢 !

Answers:


4

您可以使用单调三次样条:

http://en.wikipedia.org/wiki/Monotone_cubic_interpolation

有关在multi-D中执行此操作的说明,在这里:

http://dl.acm.org/citation.cfm?id=1285766

另一个选择是加权基本非振荡插值;舒志旺(Chi-Wang Shu)最近发表了一篇有关该主题的评论文章。


我检查了多维单调三次插值纸,然后找到了适用该方法的强大先决条件:>提供插值数据的节点等距分布,或者>从> [遵循严格单调,连续的一对一映射0,n]到插值间隔。显然,这对于我的一般3D流场不是正确的。我会挖掘其他参考,谢谢。
FrenchKheldar

2
这是我认为David所指的文章
Matt Knepley 2011年

是的,马特,就是那个。
大卫·凯奇森2011年
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.