请教关于浸入边界法（IBM方法）的问题

Rachel0096

请教各位，关于浸入边界法(IBM)（直接力法）的问题
我最近在用rhocentralFoam求解器做圆柱绕流，结果后处理速度云图时，会出现如图所示的波动，我采用的是下面的边界条件，想问一下是什么原因？如果我进一步缩小时间步长或者调整网格大小会有帮助嘛？

0/U

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0.3 0 0);

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform (0.3 0 0);
    }

    outlet
    {
        type            zeroGradient;
    }

    bottom
    {
        type            zeroGradient;
    }

    top
    {
        type            zeroGradient;
    }

    defaultFaces
    {
        type            empty;
    }
}

李东岳

你是调整物性把音速改了么？

Rachel0096

@李东岳是的，李老师。音速是1m/s。
下面是我的thermophysicalProperties

thermoType
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

// Note: these are the properties for a "normalised" inviscid gas
//       for which the speed of sound is 1 m/s at a temperature of 1K
//       and gamma = 7/5
mixture
{
    specie
    {
        nMoles          1;
        molWeight       11640.3;
    }
    thermodynamics
    {
        Cp              2.5;
        Hf              0;
    }
    transport
    {
        mu              0.00933;
        Pr              1;
    }
}

李东岳

你这个看起来就是压力波的传递导致的。近音速可能会这样。openfoam有无反射边界条件，不过用起来效果一般。你可以试试调整出口的边界条件，也可以试试把计算域延长一些

CFD中文网

请教关于浸入边界法（IBM方法）的问题