OpenFOAM膨胀喷管边界条件设置

[xuluxi18]

一个二维膨胀偏流喷管的算例。我先画了一个粗网格，喷管内流动算的还可以，但外部流动发散了，并且和其他Fluent计算的偏差很大。求助各位大佬，出口的边界条件应该怎么设置呢？
使用的是rhoCentralFoam求解器，喷管入口为总压边界条件4MPa，总温2000K，外界为大气压。
这是我的边界条件：

算出来的外流演化过程很奇怪：

这是其他用Fluent算的，出口设置的是pressureoutlet。

[李东岳]

网格处理到3万以下，可公开的话，上传上来我给你debug

[xuluxi18]

谢谢李老师，我用的是of9版本。penguan_RANS_org.zip

[李东岳]

好 明天给你debug

[李东岳]

你这个我跑起来了，但是瞬态跑的太慢了，你那个outlet，有参考的远场值么

[xuluxi18]

文献里给了一个大气压

[Lubing]

这里应该是有一个明显的设置错误，喷管是轴对称结构，但是你这里算的时候使用的是二维平面结构，应该改为wedge结构

[xuluxi18]

哦哦 是的，感谢指正。我改一改

[Lubing]

是的，二维和三维的能量耗散是不同的，所以结果上会有很大的差异。另外，你可能要重新画个数模以满足OpenFOAM中的wedge类型的使用要求

[李东岳]

还有，你的网格我看了一下，长度1000米，你要调整尺度

[xuluxi18]

我在ICEM里尝试将已经画好的二维网格旋转成三维网格，很容易出现一些WARNING。我试了一下最简单的矩形二维网格，网格输出和导入of没有问题。可能是之前画的二维网格的问题T~T

[xuluxi18]

我重新创建了一个楔形网格

但是计算的时候显示wedge不是一个平面，这个问题之前有人遇到过吗？

[Lubing]

关联的问题，线没有关联好，放到了wedge类中了

[yingqing]

温度和压力的边界条件与内场压力差非常大。要不要试着从小的入口压力开始跑，稳定之后在增加压力。这种操作方法，我在别人的喷管案例里面看见过

[xuluxi18]

感谢大家，应该是跑通了，主要是没有使用wedge边界的原因。
用ICEM二维旋转小角度生成三维网格很容易出现报错，使用of7的makeAxialMesh工具实现楔形体网格很方便。

[xuluxi18]

大佬们，这个算例用rhoCentralFoam能跑通。因为后面想加入颗粒计算，现在在同样的边界条件下我用rhoPimpleFoam计算时，却一直卡在这里，入口的压力也不对，总压给定的是3e6。这是什么问题呢？