我自己暂时的理解是，程序中fjfj.phase()/fj.x()和fkfk.phase()/fk.x()指代的是j和k的颗粒数目，也就是右边公式中的Nj和Nk。但是其中的fj.phase和fk.phase是什么含义就有点不太明白了。还有就是聚并速率的量纲是m3/s的，程序里面为什么还要再乘以i颗粒体积fi.x()？这样子的话量纲就不是每秒了。

各位老师好，

我用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

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