@李东岳 感谢您的回复，我的计算的是和压力温度都相关的，但是主要是压力；在您的提醒下，这两天找到了一个相关的包含空化源代码的文章和一篇冷凝没有源代码的文章，目前正在学习，不知道是否是您所指的那篇链接文本链接文本

@hurricane007 这样做的确可以用。相当于使用linux的命令行，windows的UI。

@try-harder 谢谢，谢谢！:xiexie: :xiexie:

@李东岳 谢谢李老师的回复~ 教师节快乐~ 是一个标准的channel。但是雷诺数非常高，Re=80 million, 没有8000万网格，是NASA官网上的一个标准算例，独特的地方就是雷诺数非常大

Capture1.PNG
Capture2.PNG

如上2图，白色的点为测量数据点，每个点都有ux,uy,uz
先在inlet 上通过topoSet 专门设置出了一个小的区域，
图1是方形，图2是根据测量区域画出来的一个正方形加两个圆形的区域。

现在最大的问题是，你可以发现topoSet的边界特别明显。。
这种情况应该通过什么办法能解决呢？

为什么用topoSet的原因是如果直接在inlet 上设置为timeVaryingMappedfixedValue 这个边界条件的话，code 会直接插值到最边界。。结果看起来更有问题，topoSet这个解决方法一定程度上规避了interpolation到最边界的情况，但是会造成这个patch 和他原来关联的patch的边界问题。。

求大佬指点迷津

是不是模型的xyz搞翻了，可以先用paraview看看。

或者（16 2 1）试试

@李东岳 So cool!

对y+、u+感兴趣的，可以直接看看openfoam tutorials里面incompressible/boundaryFoam里面的自带算例，有处理y+、u+的图的脚本

@冠竹 自行尝试下载metis-5.1.0.tar.gz试试

@李东岳 我重新定义了k_, 将其强制变成const，但是还是报错，错误如下：

error: no matching function for call to ‘ddt(const alphaField&, const rhoField&, const volScalarField&)’ + fvc::ddt(alpha, rho, k_) ~~~~~~~~^~~~~~~~~~~~~~~~

可见k_前面已经有const，但是还是不行，可见应该还是要修改alpha与rho的类型？但是我重新赋值定义alpha， rho为volScalarField, 会报初始化错误，不能直接这样改，修改方式如下：
const volScalarField alphat = this->alpha_; 再在fvc中调用fvc::ddt(alphat, rhot, k_)

@tianyn 不客气。欢迎加入OpenFOAM开源大家庭 :duang:

@yingdong
问题已经解决
5df70411-c33b-43b4-a648-6c2eb39b319b-图片.png
（图示网格质量较差）
最开始我建立模型使用的一半区域（使用对称边界symmetryPlane），全区域就不存在这个问题了

请问老师有相关的dynamicMeshDict中dynamicFvMesh dynamicOversetFvMesh;相关的资料么？
感谢老师指点@队长别开枪

1630305503(1).png
我也尝试了其他的边界条件，只有这个运行出来的结果是最接近现实的，但是他的出口压力边界我设置的是固定值

@浪迹天大 在 运行时动态加载库 中说：

@Micro 我有个朋友写过相关的博文，感兴趣可以看看
https://zmeng.tech/2020/10/DynamicLibraries/

https://zmeng.tech/2020/10/DynamicLibraries/

@samuel-tu =等号类似的操作都需要correctBoundaryConditions()
@李东岳 请问岳哥，为啥 = 号等类似的操作都需要显式更新下边条呢？能不能稍微解释下，感谢！

自己顶一下，有没有哪位老师给点意见，应该是网格拉伸方向错误了，但是不知道如何下手

@tens 非常感谢，好好用:xinxin::146:

请问您解决了吗，想问一下怎么解决的，遇到了同样的问题～

接楼上的发言，Info在并行时只会输出master processor的信息，所以把Info改成Pout应该就可以了

@zp气固两相流 博主这个问题解决了吗？想知道后续