Palabos论坛笔记：升力系数的计算

地址：https://palabos-forum.unige.ch/t/lift-coefficient-of-a-wing/2445

提问者表示自己通过IBM方法计算一个翼的升力系数，并且使用如下代码：

forceConversion = param.rho * (param.dxFinest * param.dxFinest * param.dxFinest * param.dxFinest) / (param.dtFinest * param.dtFinest);
force = -reduceImmersedForce<T>(*container, iSurface) * forceConversion

他疑问这个force是不是他想要的，还是说是无量纲的力。同时表示自己用到了网格加密，maxOctreeLevel提升1，force几乎会翻倍。

用户回答：force是你想要的，应当使用如下代码：

forceConversion = 2* param.rho * (param.dxFinest * param.dxFinest * param.dxFinest * param.dxFinest) / (param.dtFinest * param.dtFinest);

提问者表示：如果他是3d的，翼的面积为1^2 * AspectRatio，是否代码应该如下所示。

forceConversion = 2 * param.rho * (param.dxFinest * param.dxFinest * param.dxFinest * param.dxFinest) / (param.dtFinest * param.dtFinest) / AspectRatio;

用户回答：如果想得到升力系数，就用" / AspectRatio"，不用就是force。

至于为什么网格加密会导致计算结果翻倍的问题暂未解决。