升阻力系数计算问题

sea

在二维圆柱绕流计算中，在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的。请问正确的步骤应该是除以哪个数呢？

sea

补充：来流速度U=1，粘度设为1/2000=0.0005

Joker

lref设置为直径D=1，Aref=D×径向长度，改成2.68328就ok了，你使用什么软件导入的网格，fluent.mesh导入的时候会告诉你径向拉伸了多少，你再尝试一下吧https://blog.csdn.net/weixin_39124457/article/details/120187874
你可以看看这个网站设置的东西

sea

@Joker 谢谢，我改一下试试

CFD中文网

升阻力系数计算问题