@李东岳 , 李老师您好，我已经在chtMultiRegionFoam中植入了第一张图所示的拉普拉斯及泊松方程，

QQ截图20210619102650.jpg

QQ截图20210619102636.jpg

实现了对电势在流固区域的耦合求解，但电荷的求解似乎不对，然后电场力的求解结果也和参考文章里面差的太多，但我边界设置理论上是符合文章要求的，第二张图和第三张图是我和参考文章计算结果的对比。参考文章及我的算例和编译的求解器链接放在这里：https://pan.baidu.com/s/1U8NBBon8CFVvDoF9Fsi9UA 提取码：mppz
麻烦李老师能抽空看一下，检查一下问题，不胜感激。一个头俩大:136: :136:

QQ截图20210619095212.jpg QQ截图20210619095249.jpg

@conghongchuan 提高一下热度，我可不可以理解成：把气泡替换成旋转机械，就是旋转速度了。

关于网格尺度小了以后interpolation和mapping的问题

这个算例太酷炫了 小液滴都算出来了 还特别圆

@chengan-wang 在 关于codedFixedValue使用 中说：

我定义了ymax，之后如果我使用if ( Cf[faceI].y() > ymax ) 设置相分数，时间步长越来越小。但如果改成if ( Cf[faceI].y() > 0.1 ) 一个具体数值，情况就好很多。我的使用方法是否有问题呢？

理论上不应该有任何问题 :jingya:

@浪迹天大 这样嘛。。。好的，谢谢！

如题，想用Foam-extend里的浸没边界法（Immersed Boundary Method）求解器，解决水与弹性杆的流固耦合数值模拟，有人有这方面经验吗？

@static 你好，我最近也在做三维飞机的相关数值模拟，可否给个联系方式，请教一些问题。谢谢，这是我的qq：3310792254

@李东岳 很抱歉才看到老师的回复。感谢老师，我已经自己解决啦:xinxin:

图挂了

@李东岳 东岳老师，我在dynamicMeshDict里设置了constrints和restraints两个条目，控制圆柱仅在y方向移动，不转动，如下：

constraints { yLine { sixDoFRigidBodyMotionConstraint line; centreOfRotation (0 0 0.5); direction (0 1 0); } noRotation { sixDoFRigidBodyMotionConstraint orientation; } } restraints { verticalSpring { sixDoFRigidBodyMotionRestraint linearSpring; anchor (0 0 0.5); refAttachmentPt (0 0 0.5); stiffness 2.1055; damping 0; restLength 0; } }

不知道您说的是不是这个意思？

@李东岳 感谢李老师百忙之中指点！这么来看correction(fvm::ddt(p))就代表前后时间步压力差值的时间变化率啦

各位大佬好，想请教一下，我想显示一个密度、体积不断变化的物体的质量，请问大家有什么好办法吗？

屏幕截图 2021-06-05 161820.png openfoam自带案例interPhaseChangeFoam求解器下的cavitatingBullet，该图为我跑完0.05s后的水相分布，和0s时没啥变化,请问我哪里错了吗？

@李东岳 感谢李老师回复，我先试下

@李东岳

我用下面的命令解决了 ，具体为啥不清楚..

sudo ntfsfix /dev/sdXY

https://askubuntu.com/questions/145902/unable-to-mount-windows-ntfs-filesystem-due-to-hibernation

@李东岳 好的李老师，我明白了，特别感谢您！！

@李东岳 谢谢老师的回复！

@鲸落 其实也不是只显示一层。我想你跟我当时一样，感觉摄入粒子很多，其实由于体积太小，所占的相分数其实不大。那么最后就会只有一层。如果有很多粒子的话，最后会形成很多层，因为只要满足你预设的最大相分数，它就不会再允许粒子进入，那么就会在其上方的网格中开始堆积。

@tens 学到了，感谢老师的指点！:chitang: