CFD中文网

@Afr1yne 好的好的非常感谢 我已经找到了

问题已解决

@李东岳 感谢李老师的回复，这个我看过，关于DMD里面每个文件夹的意思，我有点不太理解，这个有相关的教程吗

@AppleKiller
李老师的书《无痛苦N-S方程笔记》第226页有推导。http://www.dyfluid.com/theory.pdf

谢谢您的解答，在openFoam v6版本里icoUncoupledKinematicParcelFoam这个求解器里，有pairCollision(软球模型),我现在尝试用软球模型先试试。
sprayFoam里面是有O’Rourke碰撞的，但我做的东西和化学反应没啥关联，在后面做的过程中如果pairCollision没啥进展，可以借鉴您的建议把O’Rourke移植过去。

@JamesCheng 大神，我现在也在想做颗粒的碰撞与聚并，您现在有眉目了吗

@李东岳 在 可凝气体的冷凝数值模拟 中说：

@杨英狄 在 可凝气体的冷凝数值模拟 中说：

高温水蒸气在液态水中的冷凝及流动

不太好理解。水蒸气在水里冷凝然后流动？

李老师您好，我想做高温空气或燃气从喷管流出后水蒸气遇冷凝结雾化的仿真，您有没有了解的求解器比较适用这个问题？谢谢

@李东岳 感谢老师提醒！忘了改残差文件了哈哈哈哈40dc84c4-ea49-45e1-bec3-6948e6a46f0b-image.png

这俩体明明是有接触的，但就是合并不了0f0e9e47-e759-425d-890a-5397db277833-image.png

@Jasper-0 建议更换org的openfoam11+版本即可，改动没有想象的那么大

哈哈，多谢老铁帮吹流弊！你们都是真流弊，我是吹流弊 :140:

@李东岳 好的谢谢李老师

楼主解决了吗

@李东岳 好的，:xinxin: 感谢李老师解答

求助各位大佬们，我计算了一个冷态的流场，中心氮气携带颗粒，周围为旋转流入的氧气。在计算一段时间后，我想用setFields进行点火。一开始尝试给到1800K，但很快发散了。目前尝试先加到600K。
但是在测试中，会出现奇怪的升温（关闭了combustion、chemistry）。之后会一直升温，直到温度超出物理限制后报错。测试了很多方面也不知道是为什么。我观测到的是，这个升温的区域和速度场有一些近似，但不知道有没有什么关系。求助各位老师同学们指点一下迷津，实在是想不明白了ToT。

冷态流场速度图：
0234e30a-07ff-4e2e-be53-71bf8ed6a5fb-image.png
全场升到600K后，计算一步的温度图（入口和壁面温度边界都给的300K；出口给的zeroGradient但出口在很下游，图中并未截取到）：
5cc79fe6-f5c3-429a-97cf-1b13101abad9-image.png

如图，报错提示找不到edge213 214
b0083f92-0d11-4f6c-a54c-5290161735a9-380e52f758ab5917db2f5c22ab1bf81a.png

但我选择的edge编号是213 214-1
a7290fc7-3e13-4f61-a2c0-91abfa977032-a1c5e1665748315b0eba74a5a64ffe9a.png

@李东岳 在 OpenFOAM中移动接触线问题（接触角）的边界条件设置 中说：

我去，竟然来了！！

别人提问的话，我这边能收到邮件呢。工作太忙，顾不上逛帖子了。

另外的话呢，对于像OpenFOAM和fluent这类使用CSF模型来计算气液两相流的话，低速问题时虚拟流动（spurious current）影响较大，做学术研究的话精度恐怕不够。建议用Basilisk去做，精度很高。但是工程问题的话，商软或者OpenFOAM是可以的。其他看到的吧友可以借鉴一下。
我目前已经不做单纯的界面流动问题了，研究方向稍微有了变化。

看了一篇论文，说是在openfoam上搞了个dem扩展，原文是这个:
OpenHFDIB-DEM.pdf

然后对应的代码库在github是这个：
https://github.com/techMathGroup/openHFDIB-DEM

表面上看起来好像挺有用的，但是毕竟没用过也不大知道这东西的弊端，有没有大佬尝试过这东西，评价一下给小弟开开眼界。

@dabiu 在 关于icoFOAM中通量phi的更新 中说：

icoFOAM中更新通量phi采用的如下语句（注：我理解这个语句）

if (piso.finalNonOrthogonalIter()) { phi = phiHbyA - pEqn.flux(); }

为什么不在下面这条语句的后面，直接用phi = fvc::flux(U)来更新通量，这样更新跟上面那条语句有什么区别吗？用phi=fvc::flux(U)会有什么问题吗？

U = HbyA - rAU*fvc::grad(p); U.correctBoundaryConditions();

初学icoFOAM，求各位高手解答，OpenFOAM用的版本为OpenFOAM 7.

我觉得这里的phi是经过Rhie-Chow插值修正过的面通量。

@李东岳 对哦，谢谢您