CFD中文网

感谢分享

@李东岳 谢谢李老师，附件是我的Case，麻烦老师您帮忙看看哪边出问题了。
代码的一些说明：
1）几何模型中Heatsink是通过读取0文件夹下gamma场进行生成的；
2）UEqn.H代码中增加了源项，其中si是阻流系数，用来区分Heatsink（si=alphamax）和Air（si=0）；
3）TEqn.H代码中的DT通过插值方式分别在Heatsink和Air区域赋予不同的值。

Case.zip
附件：自然散热Case

5eed8fd4-9397-41d2-b885-8959ce6871cf-image.png 图1
各位前辈，我将fluent官方案例的设置（欧拉临界热流沸腾模型）改成VOF+Lee模型后就模拟不出上图的温度趋势了。f0dd19fc-adb7-42f2-9245-496794390d68-image.png 图2
是我设置的不对还是VOF+Lee模型本身就模拟不出图1的这种效果呢？

@liujm 在 编译overWaveDyMFoam失败请求援助 中说：

@chenboyao 在 编译overWaveDyMFoam失败请求援助 中说：

@liujm 你对照当前的wave2foam中waveFoam求解器和interFoam求解器，然后你看看修改了哪儿，你就把修改的地方复制到overInterDyMFoam中，在把修改好的文件复制到wave2Foam存放求解器的地方，重新编译一下wave2foam就行

@chenboyao

Hello, 我说下我的编译方法，避免出错，您看看哪里有问题：

把multiphase/interFoam/overInterDyMFoam文件夹复制到solvers1912_PLUS目录下（与原有waveFoam同级），并将overInterDyMFoam改名为overWaveDyMFoam。

在overWaveDyMFoam.C的基础上对照waveFoam.C，填补缺少的代码。

做到第二步时我有两点疑问：

您提供的文件中没有完全将waveFoam.C中的代码补充到overWaveDyMFoam.C中，例如部分头文件以及如下等代码。请问略去的代码是不需要吗？ #include "readDyMControls.H" phi = mesh.Sf() & Uf(); MRF.update(); 两个文件夹中的部分同名文件也有较大区别，例如correctPhi.H等，这些该如何处理？对于waveFoam中有的文件，而overWaveDyMFoam没有的文件，是否应该复制到后者处？

我这个代码对应的OpenFOAM版本是2206，不知道1912版本适不适配

记录一下，由于fe40太老了，在ubuntu2204下面编译以下问题，给我卡住了1个小时..

sigFpe.C文件要改malloc scotch要用6.0.4，不能用6.0.0 scotch要自己编译 fe41默认lib创建so连接，要不然还是找不到

已更新：http://dyfluid.com/download.html

@李东岳 在 各位老师，请问我是算出棋盘型压强场了吗？ 中说：

知识盲区了老铁。运动大体上没问题就问题不大。不过看起来有一些小bug没处理好，也不能说肯定是压力期盼分布

运动大体上差不多，但是肯定不完全对:xinlei:

@李东岳 谢谢李老师回复，问题已经解决了:high:

@五好青年 大佬，这个问题解决了吗:135:

@lizhisongsjtu 我不太关心涡这个东西，有几个疑问，邀请李老师解惑，给我开开眼：

不同的涡判据，比如liutex、lambda、Q之类的，都是用来判别涡的，实际使用中，是哪个看起来更流弊就用哪个么？ :134:

如何判定那个涡识别方法更流弊呢？:134:

https://arxiv.org/pdf/1802.04099

刘老师Rotex的第一篇文章，就是另外一种写法的$\nabla\bfU$

@队长别开枪 太流弊了大佬，真羡慕啊

@李东岳 太好了！学习中！:142: :142: :142: :ok3:

我做的用部分泡沫金属强化石蜡相变材料的传热。模型为长60mm、宽80mm的矩形，上面为纯PCM，下面为泡沫金属石蜡复合材料，两者面积相等，设置初始条件297K；边界条件为左壁面为恒温360K，其余三个面绝热，泡沫金属孔隙率0.8。现在迭代发散，不知道什么情况，求各位解惑。附图为相关参数设置

1.png 2.png 3.png 4.png 微信截图_20240516104051.png 5.png 6.png 7.png 8.png

背景介绍的不完整，没太看明白是什么问题。

好的，谢谢李老师

@yuky 文章里看到过，https://doi.org/10.1016/j.csite.2022.102560 ，记忆中以前的老版本好像有人开发过，可以搜一下

@zhaoym 老铁，后来这个做的怎么样了，有进展么

interPlicFoam求解器已经并入扩展包geometricVofExt，仓库地址：https://github.com/daidezhi/geometricVofExt

@秋山结月 在 200万网格并行算力测试（OpenFOAM版本） 中说：

@tidedrinker 我之前也觉得可能是因为这个，但是事实上多插了两根内存条并没有改善现状。插满内存条的经费我现在也没有，只有先等等看以后能不能重新测个数据了。

4根内存条也不至于这么慢，插4根应该是每个三缓1根，我记得这个型号一共有4块三缓，8个核心共用一块。插4根内存，8核也应该在135s左右。没改善估计是位置插错了，内存通道插满，32核成绩应该在70s左右。

并行计算时候，每个进程掌握的网格就固定那些。
举例b点只存在于processor x1，不在processor x2当中，processor x2自然无法获取b点的信息。
要想processor x2获取b点信息，可以用reduce函数，也就是进程之间通讯一下。
首先初始化intU为零向量，然后赋值，再reduce。

point b (1,1,1); label cellID = mesh.findCell(b); vector intU(0,0,0); intU = U[cellID]; reduce(U,sumOp<vector>());

具体代码你再看一下，我这个是思路，代码可能哪里不对。