请教大佬们关于后处理中的阻力计算

[Noob]

大佬们好，我想来请教一下关于OpenFOAM后处理计算阻力的问题。
我用simpleFoam计算了个稳态问题，并打算利用后处理计算阻力，使用的control如下：

但是有人回复说对于不可压缩求解器rhoinf要始终设置为1，具体请看此图片

但是我设置rhoinf是998.8的时候得到的阻力值确实是我想要的，我的理解是simplefoam求解的时候用的是归一化的rho，然后这个rhoinf只是后处理的一个缩放因子用来得到真实密度下的阻力的，不知道这样理解是否有误

迷糊好久了还请各位大佬们指点迷津，祝大家科研顺利，生活愉快，双节快乐！

[xpqiu]

不可压缩求解器计算 force 时rhoInf设置为实际密度是对的。

[Noob]

@xpqiu 非常感谢大佬的回复！也就是说我设置的没问题呗，那上面那个英文回复是啥意思呢，再次感谢啦！

[xpqiu]

@Noob 需要看一下上下文才能判断，可能是很久以前的版本跟现在不一样。

[Noob]

@xpqiu
上下文大概是这样，我用的是OpenFOAMv1812，我还是感觉rhoinf对总阻的计算就是后处理时的缩放因子，貌似不会影响求解过程

[xpqiu]

@Noob 对，不会影响计算过程。
不可压计算中的p，粘度（nu，nut）都是除以了密度的，所以直接用p和nu(以及nut)计算得到的压差力和剪切力也都是除以了密度的。但是使用 libforces 计算受力时，如果检测到是不可压求解器，则会乘以一个 rhoInf，来得到单位为 N 的力。

[Noob]

@xpqiu 感谢大佬回复！祝大佬中秋快乐，万事顺遂啊！也就是说rhoinf要填流体的真实密度而且单位是kg/m³吧，那我runscript脚本写的998（水）是合理的吧

[xpqiu]

@Noob 是的，合理

[Noob]

@xpqiu 感谢感谢