使用dynamicMeshDict.sixDoF出现could not find rho问题

[jxzs_zy]

请问下，我现在使用dynamicMeshDict.sixDoF计算小球在气体中的振荡但出现了以下报错，
然而我在transportProperties中定义了rho：
也在fvSolution定义了rho的求解：

。

外文论坛有好多人也在问这个问题，有说版本问题或者bug的：https://www.cfd-online.com/Forums/openfoam-solving/178985-problem-calculating-forces-blades-2-finding-forces.html

请问各位这是什么问题呢？

[李东岳]

你用的什么求解器？

[jxzs_zy]

@李东岳 老师，您好。使用的是pimpleFoam。

[李东岳]

pimleFoam本身没有密度。算例要是不大的话可以发上来我给你看一下。

[jxzs_zy]

@李东岳 Thank you, Dongyue teacher, the calculation is relatively large. It seems to be this problem, can you help see, how should this be changed? I have already changed the first half of the sentence, I don't know how to change the second half?
谢谢您，东岳老师。文件有点大。
这应该是红线出的问题，我改了前半句，但是后半句不知道怎么改，您能帮忙看看吗？
https://cpp.openfoam.org/v11/forces_8C_source.html

这个rho应该是刚体的密度，我看V9\V10版本中只有dynamicMeshDic也使用sixDoFRigidBodyMotion来求解6DOF问题，这里的rho应该是刚体密度？
tutorials/incompressible/pimpleFoam/RAS/wingMotion/wingMotion2D_pimpleFoam/constant/dynamicMeshDict

老师：1.只用 dynamicmeshdict文件可以求解6DOF问题吗？
2. 刚体密度后半句应该怎么改？
万分感谢老师！！！

[jxzs_zy]

李老师能帮忙看看吗？

[李东岳]

@jxzs_zy 这个最好能搞个小算例，然后我给你debug一下

[jxzs_zy]

好的，老师。我用少一点的网格，把文件搞小点

[jxzs_zy]

老师，这个是算例文件，辛苦您看下。
Exm1.7z
感谢老师。

[李东岳]

@jxzs_zy 用哪个求解器跑？pimpleFoam么? pimpleFoam不需要重力，还是rhoPimpleFoam？用rhoPimpleFoam的话边界条件为什么是fluxPressure?

[jxzs_zy]

pimpleFoam,老师。抱歉，刚学of没多久。刚体会受到重力，在.sixDoF中定义g吗？

[李东岳]

@jxzs_zy 在 使用dynamicMeshDict.sixDoF出现could not find rho问题 中说：

我现在使用dynamicMeshDict.sixDoF计算小球在气体中的振荡

好的我知道了，我回归你最原始的问题给你调试一下

[jxzs_zy]

@李东岳 好的，辛苦老师了，我在.sixDoF中也定义了g

[李东岳]

你这个需要改的太多了，首先你这个边界条件，为什么全部都是速度固定值0？ 你这样吧，参考wingMotion2D_pimpleFoam这个算例，你自己尝试调一下，都整差不多了你再发给我我再弄。目前你发给我的这个算例，需要改的太多了。

[jxzs_zy]

老师， 是全封闭状态，球初始时刻静止的，流体运动都是球带动引起的

[李东岳]

@jxzs_zy

流体运动都是球带动引起的

那岂不是球是主动运动？

[jxzs_zy]

@李东岳 是的，老师。球在重力的作用下主动运动，并且被一根弹簧拉着

[jxzs_zy]

@jxzs_zy因为弹簧力体积很小，我做了简化，用函数代替了老师

[jxzs_zy]

@jxzs_zy 抱歉老师，是弹簧的体积很小，没有画出弹簧