各位老师好，

我用openfoam4.x - pimpleDyMFoam - sixDoFMotionSolver计算刚体运动，我的几何由body1和body2两部分组成，两个body都是刚体，相连并同时运动，因此也可以看做一个刚体。我目前将其看做一个刚体，因此在dynmaicMeshDict文件中写作：

dynamicFvMesh dynamicMotionSolverFvMesh; motionSolverLibs ( "libsixDoFRigidBodyMotion.so" ); solver sixDoFRigidBodyMotion; sixDoFRigidBodyMotionCoeffs { ... patches (BODY); ... }

计算后能得到刚体上的升阻力，但是我要如何根据现有结果分别得到两个body上的升阻力呢？（如果我理解正确的话，BODY1上的气动力 + BODY2上的气动力 = BODY上的气动力）

还是说只能将其看做两个刚体并重新计算，即dynamicMeshDict中改为patches (BODY1 BODY2);

谢谢！

@冠竹simplefoam、pimplefoam、rhosimplefoam和rhopimplefoam的层流计算是不是也可以称为准直接数值模拟呢？

在v2012版本中，TEqn.H中有rhoCpPhi的定义为：

const surfaceScalarField rhoCpPhi ( "rhoCpPhi", rhoPhi*(Cp1 - Cp2) + phi*rho2*Cp2 );

这个定义是如何得到的? 如果Cp1 = Cp2，那rhoCpPhi = phi*rho2*Cp2，这个和 phi*rho*Cp有差距。

如果直接把中心点的rhoCp插值到面上和phi相乘应该表示为 fvc::interpolate(rho1*alpha1*Cp1 + rho2*alpha2*Cp2)*phi。

@杨英狄 这个也是计算粒子与壁面的作用，所以就是在particleerosion函数的基础上改的

@李东岳 李老师，fluent和其他机器计算时长的对比有吗

@cccrrryyy 让运维换了个of2106版本，速度有所提升，相当于自己电脑的四分之一速度，离期望值还是低点

@李东岳 烦请东岳老师有空看下，期待您的回复。谢谢！

@bestucan 谢谢大佬，已经解决了

看下tutorials里的Allrun和Allmesh操作步骤

@wenke 你好，你这个问题解决了吗？我也遇到了这种问题，现在也没有解决。

@qing 是这样的，openfoam完全可以这样任意折腾。它是基于文件的，只要能读取、识别到相应的文件就可以进行计算。

我思考了下，

matrix relax 虽然具体实现方式不同，但是本质不变。

你说的是对的。一楼我说的应该有问题（我更正一下）。不仅仅field relax，matrix relax，也会影响计算结果。我当时把matrix relax理解成不会影响计算结果了。

matrix relax是使得计算结果更慢的趋向于真实值，但换来一个对角占优矩阵

field relax是使得计算结果更慢的趋向于真实值

两个都使得计算结果更慢的趋向于真实值。

当时看这个公式的时候：$\phi^{n+1}=\phi^n+\beta(\phi^{n+1}-\phi^n)$，以为$\beta$变小时候，$\phi^{n+1}$还是原来那个值。然而并不是:136: :136: 好尴尬

@李东岳 李老师，您好！现在使用foam-ext 4.0,已经编译完成，计算成功一次，后续计算时一直报错（可以划分网格，计算中读取信息时出现问题）。
bed30e22-ed0d-4079-b2df-86b3002c91bf-image.png
7502d3ba-1a3d-4601-bd2a-3de6d249d09f-image.png

@zhouxu 谢谢，我尝试下。

@李东岳 东岳老师，请问这个速度能不能从10-2下降到10-4或者-6？

这里有介绍：link text

类似的需求：

How do I add some other package to the container?

The easiest way is to modify your Dockerfile and rebuild. Particularly easy if the software you want to add has an Ubuntu deb package.

里面也有 docker file 的示例。看起来这个比较可行

没有经验，一直在 linux 上跑

@李东岳 懂了，谢谢李老师

应该研究下为什么transonic需要提升对角占优

@李东岳 回答一下自己的问题。transonic下的压力方程存在对流项，对流项不利于对角占优矩阵的产生，因此transonic的压力方程要增加方程松弛。

http://dyfluid.com/rhoSimpleFoam.html

@ruibin_tongji
Green-Gauss Cell Based
9f30b62e-ad51-4bd3-ad87-c12ce78f2194-image.png
Least Squares Cell Based
8e14d4c0-98f9-4866-b57a-05a79fe25109-image.png
Green-Gauss Cell Based是最粗暴最直接的，但也许恰好适用于结构网格？我看这个Theory Guide里面是这么说的
38584698-16af-4fa1-8327-194895b34791-image.png

有没有测试过对流项的格式？

在计算多相流动网格案例时，计算一段时间后库朗数突然增大并导致计算发散。起初以为是网格变形过大导致的发散，但查看日志文件和结果发现，计算并非是结构物受到波浪冲击后达到最大变形时终止的，而是在逐渐回到原位置时终止的，此时相分数alpha.water超过了1并在急剧增大。这时结构物附近有一定的越流（如下图所示，截图截得不太整齐请大家见谅），且计算停止时（44.2s）的气液交界面显然是有很大问题的。因此我认为应该是相方程计算不收敛进而导致的整个计算终止。
微信图片_202203281625522.png

微信图片_202203281625521.png

微信图片_20220328162552.png

日志文件中相分数在44.076s开始大于1：44.076.png
最后计算崩溃：44.202.png

因此，想请问大家是否在多相流、动网格的计算中是否遇到过类似的问题，有什么办法可以解决这种相分数越界的问题吗？