CFD中文网

0_1497945059409_2017-06-20 15-50-42屏幕截图.png

震荡已经非常非常小了，

都柏林大学发来贺电！

我想在Fluent中模拟单相超音速横流射流，已知横流入口马赫数为1.6，总压为241kPa,总温为295k，射流入口的马赫数为1，总压为476kPa,总温为295k,根据我查阅的资料，对于这种超音速可压缩流动应该压力入口和压力出口边界条件。
2bd2af9c-46b0-4bfe-af9e-0fa137ddbbd5-image.png
c4ee670e-4062-481a-b52d-db5cad6d6037-image.png
但是现在对于出口压力和温度设置我不清楚应该如何设置？除此之外，对于超音速可压缩流动模拟，还需要注意什么设置细节吗？

@bestucan 你是行家呀，厉害厉害，终于知道原因了

@bestucan hhhh谢谢您，我发现我哪里错了，我画网格的时候是xoz面上的，正确的应该是在xoy画个2D的网格，然后openfoam会在z轴增加一个面。

@赵一铭 不是用软件，是编程计算，现在的问题是一个不收敛，和出现负值

@东岳 是的VOF多相流，计算域三层，最上面是空气，红色是油，下面是水2020-10-22 21-10-44屏幕截图.png

搜你那个版本的fluent的《fluent udf manual》

里面有各个函数的使用方法和最小实现。

@yfclark 谢谢！

那就在做O形网格的时候，“O穿”整个区域。就是假设那个圆柱两端无限延伸，对这个无限长的圆柱做O形网格。

如果不知道空心实心怎么选。做的时候；part里有个part，是专门放已经被删除的block；把这个part打开。就当是个全实心的圆柱。完了再关了那个part。

@李东岳明白啦，谢谢东岳老师。

感谢感谢！

如果建模的时候分别建了左边和右边，中级的面就要设置为interior

天天净整这乱七八糟的

sad....

@ibelief 谢谢，我再看一看。

@赵一铭 整个程序我用的SIMPLE算法，所以连续性方程改成了压力修正方程来求解压力，不能直接求解出密度，怎么办？

@李东岳 尖角部位，加一个竖直的翼墙，这样平分后角度好点，但是附近网格不再是水平的，这样对interFoam计算有影响吗？

谢谢李老师的回复！因为雷诺应力不仅出现在动量方程里面，也出现在turbulent kinetic energy equation 与 turbulent frequency equation中，我的想法是如果我分别修改湍流模型的文件（kOmegaSSTBase）与动量方程的文件（linearEddyViscosity）就会要定义两次非线性项，所以我选择就在kOmegaSSTBase中更新turbulence->divDevRhoReff(U),在其中加入非线性项，参考的资料是OpenFoam里面的nonlineareEddyViscosity。由于这样修改需要在kOmegaSSTBase.C文件的constructors中加入eddyViscosity 也就是上图中蓝色后面的部分，我直接在原来的代码Turbulence后面加这一项会报错，想请问这怎么解决？非常感谢李老师耐心看我的问题，不知道我表述清楚没有？

请问有前辈知道gmesh如何能划分如图所示得网格吗，图1是图2的block

block-tuya.jpg model-tuya.jpg