各位老师好，
我用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);
谢谢！
@李东岳 谢谢李老师的回复！抱歉我的上个残差图绘制有问题，正确的应该是这样
从上到下依次是速度残差，压力残差和湍流统计量残差。
上面那个code是我随意找的，我目前使用的code是：
PIMPLE
{
correctPhi yes;
nOuterCorrectors 2;
nCorrectors 2;
nNonOrthogonalCorrectors 1;
turbOnFinalIterOnly false;
}
relaxationFactors
{
fields
{
p 0.3;
}
equations
{
"(U|k|omega)" 0.7;
"(U|k|omega)Final" 1.0;
}
}
我主要做了以下修改：
（1）调整 nOuterCorrectors=2，结果几乎不变
（2）turbOnFinalIterOnly=false，即每次PIMPLE loop都计算湍流，结果几乎不变
（3）去掉全部的松弛因子，结果几乎不变。
在模拟中，CFL保持在0.5-1.5之间，但是和实验结果对的不好。李老师觉得还有可能是哪里出了问题呢？
各位老师好，我关于fvSolution中PIMPLE算法的参数设置有几点疑问。下面是我找到的一个例子，有以下问题：
PIMPLE
{
momentumPredictor yes;
correctPhi yes;
nOuterCorrectors 1000;
nCorrectors 1;
nNonOrthogonalCorrectors 1;
pRefCell 0;
pRefValue 0;
moveMeshOuterCorrectors yes;
consistent yes;
residualControl
{
U 1e-6;
p 1e-6;
}
}
relaxationFactors
{
fields
{
p 0.3;
}
equations
{
"(U|k|omega)" 0.7;
"(U|k|omega)Final" 1.0;
}
}
另外，我有一个2D结果湍流统计量k的残差逐渐增大，不太清楚是什么原因造成的，能请各位老师帮忙看下吗？我猜测大概率是PIMPLE算法设置的问题。。。谢谢！
@cresendo 我重新查看了source code，aRel的松弛方式确实如您所说，感谢指出错误！我对比了不同aRel的结果，发现Rel = 0.4-0.8都可以正常运行且结果几乎相同，但是aRel >0.8之后会发散。我不确定为什么wikidict上说太小的aRel会造成不准确的响应，至少我的结果并不受aRel取值的影响。相反，恰当的aRel可以使本来发散的case收敛。我之前没有指定过aRel，也就是取默认值1.0，因此比较好奇为什么您会设定aRel=0.4而不是其他值，是根据模拟的经验吗？谢谢！
@cresendo 我按照您的建议，调整了accelerationRelaxation 0.4，计算可以正常进行，CFL数小于1， 流场看上去也正常。
但有一个小疑问，通过查询dynamicMeshDict文件的官方描述，得知accelerationRelaxation直接减小刚体的加速度，比如在某时间步solver求出的加速度为10m/s^2，如果设置accelerationRelaxation=0.4，那么实际应用的加速为4m/s^2，这是一个很大的松弛。而官方描述中建议的accelerationRelaxation范围是0.9-1，并且提到 “Be careful with this accelerationRelaxation. Too low of a value will mean that the Body does not respond to the fluid forces correctly”。因此我担心设置较小的加速度松弛会不会对结果产生太大影响？因为没有实验结果作为参照，所以无法判断结果的准确性。不知道您一般是如何处理这种误差的？谢谢！
@cresendo 请问您说的dynamicMeshDict里的松弛因子是 accelerationRelaxation 0.9 和 accelerationDamping 0.95 吗？
@cresendo 谢谢您的建议！去掉p松弛之后，算例确实可以计算下去，但是CFL保持在20左右，且速度场和压力场也明显有问题。所以还是不太对，不过还是感谢您的帮助！
@李东岳 谢谢李老师的回复！我尝试使用高版本的OF计算，发现相同设置下算例可以正常计算，不确定这是不是低版本的OF涉及动网格时存在bug。不过有个case在所有版本OF下都发散。。。能麻烦李老师有空的话帮忙看一下吗？算例文件case.zip （大约在time=1.43发散）十分感谢！！