OpenFOAM 里计算Cp，Cf等参数

[莫晓柒]

我是用rhoSimpleFoam计算螺旋桨，但后处理时我需要得到螺旋桨的推力，扭矩，以及螺旋桨叶片表面（或者某一截面）的叶片表面压力系数Cp，阻力系数Cf等值。我是要在controlDict里面指定什么函数吗？在网上也没找到合适的教程，求各位大神帮忙解答一下。

[李东岳]

functions
{
    forceCoeffs1
{
    type            forceCoeffs;

    libs            ("libforces.so");

    writeControl    timeStep;
    timeInterval    1;

    log             yes;

    patches         (cylinder);
    rho             rhoInf;      // Indicates incompressible
    rhoInf          1;           // Redundant for incompressible
    liftDir         (0 1 0);
    dragDir         (1 0 0);
    CofR            (0 0 0);  // Axle midpoint on ground
    pitchAxis       (0 0 1);
    magUInf         1;
    lRef            2;        // Wheelbase length
    Aref            6;        // Estimated
    /*
    binData
    {
        nBin        20;          // output data into 20 bins
        direction   (1 0 0);     // bin direction
        cumulative  yes;
    }
    */
}

}

[莫晓柒]

李老师，这个function可以输出Cp吗？

[李东岳]

Cd、Cl是可以的，Cp我记不太清了，现在手头没有openfoam，你试试

[hoversoar]

@莫晓柒 不能输出Cp。你可以在后处理的时候对压力p进行处理直接获得Cp。
如果要获得Cpmean的话，我是运行过程输出压力平均场Pmean，然后同样后处理得到Cpmean。
希望有所帮助

[RebelYoung]

@东岳 请问，这个是后处理实现吗？postProcess执行后说找不到这个函数。。

[Zhoujq]

@hoversoar 你好，请问下这个Cp后处理是如何实现的呢，非常感谢

[zhouxu]

@rebelyoung 放在controlDict文件的functions中

[calrel]

@zhoujq 我之前算过Cp，你网上搜一下应该可以找到类似教程，大概就是导出pmean数据，在matlab中写几行代码可以把坐标转换成角度。很方便

[Zhoujq]

@calrel 非常感谢