感谢分享！

好的，谢谢东岳老师指点！

@东岳 感谢老师回复，耦合求解器后面放弃了😂

目前想做一个俯仰翼型的算例，使用dynamicMotionSolverFvMesh网格运动求解器和displacementLaplacian动网格求解器。看到一个圆柱振荡的算例，想跟着这个算例进行修改，振荡算例的pointDisplacement文件的边界条件如图：69ecd26a-81fe-4165-9daa-fcad778d8ba4-image.png 如何把type进行修改，变为转动的类型，对应需要的参数又有哪些，请大神指点一下

@东岳 :140:

@东岳 好的

我在计算圆柱涡激振动，想要设置随时间变化的弹簧刚度（在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的过程没有表现出来。

请高手指教！！

@东岳 上次问完问题，就没用并行计算，这两天要用了，哈哈哈

我要刷火箭

矩阵求解器发散了，原因很多，不好判断

@程迪 scalar t 没有用到？

可以尝试自己写代码试试啊

这种由发布者发布的视频资料，类似文章的通讯作者，应该具备解答问题的责任，所以你最好联系视频发布者，在这里问，我不好评价别人的工作。
小提示：你的看起来比较粘。

@东岳 瞬态问题的话，是不是就是内残差作为收敛标准，初始残差不看的？

问题背景：将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参数？请大神指教！！

库朗数如果放在80年代可以引一引，不过现在库朗数基本成为了一个必备内容，不引用应该也很安全

:144: Yes 我同意 @huoyanxiadebing 说的。每个湍流模型有自己适用的流型。这方面需要特定的研究领域的专家才能有很深的理解。

@kimy 在 Code production of k ? 中说：

@东岳 Many thanks. In fact, I understood this procedure. One confusion is that the formula implemented in openfoam is not consistent with that mentioned in the Spalding' paper.WeChat Screenshot_20200504110942.png

Nobody is interested in this issue?

@wwzhao 好东西！谢谢

小白刚接触openfoam，因为做振荡浮子有关的课题，所以在学习floatingObject算例。
有几个问题想向诸位前辈请教：
1、constant文件夹下存在两个动网格文件，是两个文件同时起作用吗？还是只有dynamicMeshDict文件控制动网格？
77152af7-c361-4b22-a3bd-a9778793fe8b-image.png
2、dynamicMeshDict中joint部分的参数含义是什么？我感觉与六自由度的constraint有关，不知道对不对？
f4c7d673-3d16-46df-b663-e0acf8965d52-image.png
3、如何通过pointDisplacement获取浮体的运动位移？是否与forces一样可以在controldict文件中设置functionObject来获取？