CFD中文网

各位老师好，

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

谢谢！

@WYing 在 fvSolution中PIMPLE算法的参数设置问题 中说：

@李东岳 谢谢李老师的回复！抱歉我的上个残差图绘制有问题，正确的应该是这样

Residuals_Ur5.png

从上到下依次是速度残差，压力残差和湍流统计量残差。

上面那个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之间，但是和实验结果对的不好。李老师觉得还有可能是哪里出了问题呢？

请问同学，你的问题解决了吗

@cresendo 在 计算钝体涡激振动时发散（pimpleDyMFoam） 中说：

songtao.chen@sjtu.edu.cn

谢谢！非常感谢您的建议！

@wying 我仔细看了一下，不是一个时间步长内，结构运动响应和流体域都求解多次，应该是一个时间步内，结构运动响应和流体域都求解一次。目前openfoam大部分求解器都是一个时间步内求解一次，只有少部分分为自时间步求解。

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

@李东岳 东岳老师，我在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; } }

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

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

@BlookCFD 都看的意思是说，如果我的目标是y+<1，那么y+ min max average 都要小于1是这样子吗？

@cccrrryyy 老师，我这个结果跑完。后续贴上，麻烦！

@xpqiu 谢谢您的回复！关于2，我设置的phi的tolerance是10^-6，solver是GAMG。把nNonOrthogonalCorrectors设为10之后，最后几次计算的No. Iteration 就随计算次数减小了，最后为0，我想应该是残差已经降到最低。谢谢！

@wying 要非常高的精度，比如普方法、有限差分那种7 8 9阶，有限体积法大体都不行

参考文献
Demirdzic, I.; Peric, M., Space conservation law in finite volume calculations of fluid flow.

那就在做O形网格的时候，“O穿”整个区域。就是假设那个圆柱两端无限延伸，对这个无限长的圆柱做O形网格。

如果不知道空心实心怎么选。做的时候；part里有个part，是专门放已经被删除的block；把这个part打开。就当是个全实心的圆柱。完了再关了那个part。

@bestucan 谢谢您的回复！改成其他边界条件后，顶部网格就没有加密了，确实是cyclic的原因。

@李东岳 谢谢东岳老师，我还是在win下安装paraview吧~

我在计算圆柱涡激振动，想要设置随时间变化的弹簧刚度（在dynamicMeshDict中设置stiffness）。

我采用timeActivatedFileUpdate实现这一目的，即在不同时间段使用不同dynamicMeshDict文件，比如
time = 0 - 200，使用 dynamicMeshDict_1 (stiffness_1)
time = 200- 400，使用 dynamicMeshDict_2 (stiffness_2)

输出结果表明，程序确实按照预先的设置，在不同时间段读取了不同的dynamicMeshDict文件，但是结果似乎不太对。
以圆柱振幅为例，正确的结果应该是：

使用stiffness_1时，圆柱在0-200期间已经产生周期振动，且最大振幅稳定在 a1； 当在 time=200 改变刚度为stiffness_2时，圆柱振幅有一个突变，并逐渐趋于周期振动，最大振幅稳定在a2.

但目前使用timeActivatedFileUpdate的结果是，圆柱振幅逐渐增大并最终稳定在a2，0-200的过程没有表现出来。

请高手指教！！

问题背景：将sixDoFRigidBodyMotion中linearSpring的参数stiffness从常数改成随时间变化的函数 stiffness = f(t)
f(t) 可以是分段函数也可以是任意函数，t为目前的计算时间。

因此我对$FOAM_SRC/sixDoFRigidBodyMotion/sixDoFRigidBodyMotion/restraints/linearSpring目录中的.C和.H文件进行改动。但是在获得当前计算时间这一步遇到了问题。。尝试了多种方法都无法编译，以下是我尝试的一些语句和报错：

(1) scalar current_time = db().time().value()
Error: ‘db’ was not declared in this scope
(2) scalar current_time = runTime().value()
scalar current_time = runTime_.time().timeOutputValue();
Error: ‘runTime’ was not declared in this scope
(3) scalar current_time = timeDirs.last().value()
Error: ‘timeDirs’ was not declared in this scope
(4) scalar current_time = U.runTime().value();
scalar current_time = U_.thisDb().time().value();
Error: ‘U’ was not declared in this scope
(5) scalar current_time = mesh_.time().value();
Error: ‘mesh_’ was not declared in this scope
(6) scalar current_time = this->db().time().timeOutputValue();
Error: ‘const class Foam::sixDoFRigidBodyMotionRestraints::linearSpring’ has no member named ‘db’

挺简单一个问题，在网上找了很多方法，都失败了。。。不知道是不是因为没有包含需要的头文件？或者除了修改库文件，是否还有其他方法修改常数参数为time-dependent参数？请大神指教！！

终于注册成功了，感谢李老师！