simpleFoam质量流出口出现奇点

[李子橙]

大家好，我在使用 simpleFoam 计算外流场时，利用质量流出口作为驱动，计算结果在出口处出现了奇点，以下为最新设置的 U 和 P，供各位参考，奇点的出现原因是什么？是否和边界类型有关？希望能收到大家的回复。

U文件如下

FoamFile
{
    version                                 2.0;
    format                                  ascii;
    class                                   volVectorField;
    location                                0;
    object                                  U;
}
dimensions                              [0 1 -1 0 0 0 0];
internalField                           uniform (0 0 0);
boundaryField
{
    outFlow_car
    {
        type                                    noSlip;
    }
    outFlow_inlet
    {
        type                                    pressureInletUniformVelocity;
        value                                   $internalField;
    }
    outFlow_wall
    {
        type                                    noSlip;
    }
    outFlow_bottom
    {
        type                                    noSlip;
    }
    outFlow_outlet
    {
        type                                    flowRateInletVelocity;
        massFlowRate                            constant -10;
        rho                                     rho;
        rhoInlet                                1;
    }
}

P 文件如下

FoamFile
{
    version                                 2.0;
    format                                  ascii;
    class                                   volScalarField;
    location                                0;
    object                                  p;
}
dimensions                              [0 2 -2 0 0 0 0];
internalField                           uniform 0;
boundaryField
{
    outFlow_car
    {
        type                                    zeroGradient;
    }
    outFlow_inlet
    {
        type                                    fixedValue;
        value                                   $internalField;
    }
    outFlow_wall
    {
        type                                    zeroGradient;
    }
    outFlow_bottom
    {
        type                                    zeroGradient;
    }
    outFlow_outlet
    {
        type                                    zeroGradient;
    }
}

谢谢！

[李东岳]

看起来是网格的问题。可以调节下格式，减少skewness带来的影响

[李子橙]

@李东岳 这个网格可简单了，应该不是网格的问题。

    Max skewness = 1.794448 OK.
    Coupled point location match (average 0) OK.
Mesh OK.

[李子橙]

问题解决了，之前是没发现 flowRateOutletVelocity ，把 flowRateInletVelocity 换成 flowRateOutletVelocity 即可