@李东岳 好的我简要描述下问题，就是我按照完整的形式推导了不可压缩流动的能量守恒公式。
4422b726-0b92-4704-82cd-c145d75dcc1c-image.png 但是案例结果并不满足。
我对这一项有疑问4326e809-66fe-472b-a04c-686c03c3cb2b-image.png
我在程序里输出是有值的，但是之前的回复中提到不可压缩流动中应当为零。
这一项是按照下面这个式子，考虑了速度梯度为零，去掉右边第一项积分来的
e470d379-bbaf-47b1-ba05-c1afb5d0c8d9-image.png
我发现这个式子中的三项，积分后的形式都是一致的，感到很困惑。

在openfoam 5.x下编译你试试 我那个虚拟机里面有
如果编译你都不知道怎么编译的话 你需要先学一下openfoam 然后再弄ibScourFoam

谢谢李老师，问题解决了，他们的程序调数据都要走一下虚拟硬盘，不是直连的，中间文件乱码了。

@李东岳 好的，谢谢东岳老师！

第一次在这个网站上求助各位大佬，报错图片不知道怎么看不了了，不好意思![IWE)}@6(NH98E9LRH@PDEK.png

@wangfei9088 不好意思哈，我这边没有描述清楚。是这样的，我用of2012当中reactingTwoPhaseEulerFoam里面的 fluidisedbed，用reactingTwoPhaseEulerFoam求解器计算双流体，这个文件里面第22行到后面都没有报错，以及我用这个算例，增加一个颗粒相的文件，用multiPhaseEulerFoam求解器，计算多流体，也没有报错，但是我用reactmultiPhaseEulerfoam求解器计算就报错了，说是不识别RASModel kineticTheory.不知道什么原因，如何修改。大神有遇到这种情况嘛？屏幕截图 2024-07-23 113729.jpg

resample to image这个功能在对三维流场重新采样时，在网格加密的交接处，容易出现速度为零的情况，有没有大佬遇见过指点一下7dbacbe5-e936-47cf-a860-a96d063cfe09-image.png
原加密区域：946c0aa2-5862-4850-84bb-d16ef158fc4d-image.png

@李东岳 谢谢 已经找到了. 忘了取消帖子了.

已解决。

你这个应该是处理器数量1 内核可以给到10

@李东岳 好的好的，感谢老师！！！

@coolhhh 明白了。谢谢大佬解答

那太感谢了，李博士，我的邮箱：821954980@qq.com。

@hehaohao2013 多谢老铁支持。 有帮助就好。

@zhe 在 在interFoam中通过setFields设置水体速度，但结果没有 中说：

@李东岳 不好意思啊，东岳老师，因为mesh还挺多，并且是我改动过的求解器，所以没发上来一起算算。因为我这个求解器是对比MPPICInterFoam加入了粒子模型。我一直没尝试过给案例里通过设置setFields里的速度，这次一试，发现不对。

我返回去用interFoam算了一下，水流没有问题了。我又用OpenFOAM自带的MPPICInterFoam算了一下，一样的设置条件。0时刻，setFields起到了作用，速度场没有问题。在第一个时间布长，速度全部归零了。看了压强，跟interFoam的结果完全不同。这两个求解器没有特别大的不同，只是在MPPICInterFoam里考虑了粒子的体积占比，但即使不注入粒子，也一样速度归零。不知道有没有谁看过MPPICInterFoam，这个在算水流时，有什么大的不同吗？谢谢

你好，请问问题解决了么？我也遇到了相同的情况。
用修改的interFoam设置不同相速度，0时刻没问题，计算第一个结果时速度场也归0😅

计算域是350m乘以350m乘以60m 目前把计算域缩小后可以喷射出来 如果要维持之前的计算域 是否能够通过更改一些扩散系数之类的 让他喷射出来？

@李东岳 受教了, 吃了英语不好的亏了。

请问你现在算了大概几秒啊 你可以考虑用比较大的水流强度试一下吧 接近动床的强度（U/Uc接近1）来试试看吧

我翻了一下tutorial中的案例，只找到了如下的一段controlDic中的function，是否是对其进行一定修改来达成目的呢？求各位大佬指教
1721204958853.png

@疏影横斜水清浅 您好，最近我也在学习这个帖子的问题，我想在已经o切的圆柱侧面通过patchinject注射颗粒，所以我要定义这个注射的边界面，您现在知道如何去用topoSet和createPatch去实现在圆柱侧面上画边界了吗