@楼道男孩 这就不太好弄了 要是openfoam的还能传上来让大家看看

@楼道男孩

求解是采用的自行编译的程序

什么程序？不是openfoam？

@田畔的风 大佬，参考这种风格， 链接文本，如果你可以写一个比较容易让人理解的文章的话，可以发给我我在CFD界推一下

@田畔的风 大佬这个有点意思，感兴趣写个CFD界不

看起来是固定梯度边界条件？是这样么？

我见过有设置固定值为0的。他们做吸附的有做这个东西，之前oks课我也植入过这个东西，参考 https://www.jianguoyun.com/p/DbQV8p4Q9s3ZBhjhg9EEIAA 但是为什么会是负值？

这个网格多少？用的层流求解？

@bestucan 太厉害了，这你都能找到

感谢分享，1024！

这个是用哪个求解器计算的?

@星星星星晴

@myjimp 可以，没问题

plot over line应该检测不了粒子场的信息哦

是的。这也是限制云计算最大的一个问题。来回上传下载就够受的了。相当麻烦。

另外paraview处理上亿网格，这个硬盘读写就是很慢。感觉无解了。你的是什么盘，机械盘还是固态盘。现在有的盘能干到7000的读取 https://item.jd.com/100015879487.html#crumb-wrap

@haozhibin 稳态问题没有库朗数

嗯，我也想知道大佬都怎么弄的，这个属实比较慢。看结果都要好久。

有时候可能是卡在碰撞算法了。如果没碰撞的话可能是粒子算丢了。

嗯，正负号写错了

你用meld对比一下两个mesh的polyMesh文件夹一样？

@田畔的风 我去，那面都DDR5内存啦？

如果网格完全一样的，不应该有任何差异

算的什么流场？

@香柏树 在 sprayFoam计算等压环境下点火延迟时间 中说：

因此参考chemFoam的源代码，想把泊松方程改为完全气体状态方程P=rhoRT(见修改的pEqn.H)

你说的点火之类我不太熟悉。但针对你的这方面描述，chemFoam不涉及到整个场的连续性。sprayFoam要满足质量守恒。你这样做从算法上有点奇怪。

后续的截图看起来发散了，除以0了，当然原因肯定是非常本质的。

磨损的话属实拉格朗日更好些。不过你要不通过拉格朗日的话，那应该从物理的角度，怎样在NS方程中，构建一个杂质的影响。

这个杂质对流动的影响，是如何影响的？

什么问题截图看看

https://www.cfd-china.com/topic/4790 看文章最后

单节点4核47秒，双节点每个2核72秒？

4年前竟然翻译过

openfoam有个专门的tutorial研究这个卷吸回流。flowWithOpenBoundary就这个tutorial。不过你这个，你要不要试试inletOutlet边界条件。强制无回流。但好像实际上也是会存在。理论上，如果你这个太短的话，回流不可避免。如果你不想加太多网格。那就延长出口，做非均一网格。不关注的区域减少网格数量。

一开始的time step continuity errors就很大

稳态是可能得。瞬态要保证很小。

时间步长连续性误差会随步长增长，

这个看起来设置不太对。可以先用CPU跑没问题了之后，再跑GPU

22x这个太老了。另外GPU版本我这面也没编译过。没啥经验咯..

smoothed_aggregation预条件

这个能不能换一个别的

之前论坛里面有个人残差设置的比较松，结果不好。这个帖子我一直在找，一直没找到。也是个瞬态的自由射流。他还发了云图。你设置1e-8 -9试试。不过这个我也不太确定。按道理来说不应该差太多。你跑PISO就行。库朗数小于1

简单来说，这个可以通过边界条件组合来实现一定的控制。但是最好的办法就是延长出口。

这个密度差异看起来差了不少。LES每个时间步都收敛了么

@星星星星晴 org这面一直在整合。粒子那面自行will加入之后就变动很大。org那面湍流模型、热物理模型也是变动特别大。主要是架构方面的。

@yu_tian 这个是多计算域求解器里面的互相影响的两个边界条件

你那个是LES瞬态某个时间步的结果么

@strive_dut LES相比RANS的nut、alphat都要小，普适性来讲是这样的

@星星星星晴 esi那面更新更倾向于类似项目制度，就是有个什么项目就更新一下。所以那面很多都是一些特别新的求解器，彼此也没啥关系。org这面赞助方比较稳定，经常抓住一个算法挖好几年。org这面求解器倾向于整合。esi那面由于项目制度，不太能整合。

这是什么流动密度从450变到100，液体?

@thomasshi 主要看相分数多大，固体那面最大好像在0.63左右。如果小于这个，会直到趋向于这个数为止

还没这么玩过。简单百度了一下是。你试一下是不是这样？

是不是没到packing limit

最近大家说要这个求解器解析。最近正在备课。就回顾了一下。这个求解器非常简单。跟InterFoam非常非常类似。仅有以下区别：

1）interFoam求解alpha方程添加了人工压缩项来保证尖锐，twoLiquidMixingFoam就是普通的对流方程；
2）twoLiquidMixingFoam的alpha方程存在扩散项，interFoam没有；

目前我只能看出这两个区别。所以写这个求解器解析的动力不是很足。大部分内容都相同。在这里就简单说一下吧

interFoam解析： http://dyfluid.com/interFoam.html

2204阿里源

deb http://mirrors.aliyun.com/ubuntu/ jammy main restricted universe multiverse
deb-src http://mirrors.aliyun.com/ubuntu/ jammy main restricted universe multiverse
deb http://mirrors.aliyun.com/ubuntu/ jammy-security main restricted universe multiverse
deb-src http://mirrors.aliyun.com/ubuntu/ jammy-security main restricted universe multiverse
deb http://mirrors.aliyun.com/ubuntu/ jammy-updates main restricted universe multiverse
deb-src http://mirrors.aliyun.com/ubuntu/ jammy-updates main restricted universe multiverse
deb http://mirrors.aliyun.com/ubuntu/ jammy-proposed main restricted universe multiverse
deb-src http://mirrors.aliyun.com/ubuntu/ jammy-proposed main restricted universe multiverse
deb http://mirrors.aliyun.com/ubuntu/ jammy-backports main restricted universe multiverse
deb-src http://mirrors.aliyun.com/ubuntu/ jammy-backports main restricted universe multiverse

新版OpenFOAM-10里面driftFluxFoam更新挺多的，加了几个新的粘度模型，你可以看下

你应该在算例下面加functions函数就可以：

functions
{
    #includeFunc scalarTransport
}

可以，没问题，你看下你的沉降速度模型，这个需要模化

是的，看起来除以0了

更简单的方法就是不用动求解器了，在算例文件下面的controlDict下面把这个放进去：

cacheTemporaryObjects
(
    k
);

functions
{
    #includeFunc writeObjects("k")
}

然后直接跑就行