升阻力系数计算问题
-
在二维圆柱绕流计算中,在system/controlDict里添加如下代码:
functions { forces { type forceCoeffs; libs ("libforces.so"); writeControl adjustableRunTime; writeInterval 0.5; log yes; patches (CYLWALL); rho rhoInf; // Indicates incompressible rhoInf 1; // Redundant for incompressible origin (0 0 0); coordinateRotation { type axesRotation; e1 (1 0 0); e2 (0 1 0); } liftDir (0 1 0); dragDir (1 0 0); pitchAxis (0 0 1); magUInf 1; lRef 1; Aref 1; } }
得到的无量纲升阻力系数似乎不对,圆柱直径D=1,由于OpenFOAM将二维网格拉伸为三维,就又在paraview里测量了网格z向长度为2.68328,所以这里应该把lRef和Aref都改成2.68328吗?
我暂时没有改这两个参数,也没有重新运行,只是把刚开始计算得到的Cd和Cl都除以了2.68328,但是得到的结果还是不对,趋势应该是对的,就是比实际值放大了一个比例,下图是升力系数的时程和RMS值,雷诺数2000,正常结果应该是RMS(Cl)不大于0.1的。请问正确的步骤应该是除以哪个数呢?
-
补充:来流速度U=1,粘度设为1/2000=0.0005
-
lref设置为直径D=1,Aref=D×径向长度,改成2.68328就ok了,你使用什么软件导入的网格,fluent.mesh导入的时候会告诉你径向拉伸了多少,你再尝试一下吧https://blog.csdn.net/weixin_39124457/article/details/120187874
你可以看看这个网站设置的东西
2024年4月21日 02:41
4/4
2024年5月5日 08:47