一个有毒的算例……

[durm17]

下面两张图是在模拟：一定流量的空气注入到一个封闭空间中，看这个空间里的压力分布如何随时间发生变化。


这两个算例只有网格不同，其他条件均完全相同。
前一张图为：snappyHexMesh做出来的十万个六面体网格
后一张图为：gambit做出来的六十万个四面体网格
依我们对物理情景的理解，压力的分布应该比较均匀才对，因此似乎前一张图比较正确
可是后一张图网格更多，网格质量也比较高，为什么会出现如此不均匀的压力分布呢？
附一些细节：求解器为rhoPimpleFoam，湍流模型为kOmegaSST，时间步长为1e-4（足够小），除入口处采用速度边界条件之外，其他边界上都是墙壁（U-noSlip；P-fixedFluxPressure；k和Omega壁面函数）

[李东岳]

一定流量的空气注入到一个封闭空间中，看这个空间里的压力分布如何随时间发生变化。

应该会有压力波的传递吧？类似超级赛亚人的龟波气功

[durm17]

@东岳 没错，就是您说的那样。
snappyHexMesh那十万个网格可以很明显地模拟出龟派气功的效果，红色系的高压区在整个空间里不断震荡
gambit那六十万个网格也能大概地看出压力波的传递，但是由于压力分布特别不均匀，“高压区”不是很明显……
所以为什么网格多了，模拟结果反倒不对了呢

[李东岳]

为什么网格多了，模拟结果反倒不对了呢

理论上，网格越多，结果越好。但理论不一定代表实际。这就跟理论上，欧拉全隐无条件稳定的，但实际上并不是一个道理。

对于你这个，只能看出这些了。

[durm17]

@东岳 谢谢东岳老师

[李东岳]

这是模拟的什么东西？内部构件看起来很复杂

[wwzhao]

OpenFOAM对四面体网格的支持不是太好，需要对离散格式做一些特殊设置才行。

参考：http://www.tfd.chalmers.se/~hani/kurser/OS_CFD/HrvojeJasak/DiscretisationBestPractice.pdf

[张RP]

链接挂掉了，我找了一下
http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2015/HrvojeJasak/DiscretisationBestPractice.pdf

[李东岳]

@durm17 后来模拟的怎么样了？

[durm17]

@东岳 后来就只用snappyHexMesh画的网格做了一些简单的模拟……