buoyantfoam计算发散求助

[Miraitowa]

相同的设置在计算圆柱绕流时可以得出结果，但当把网格换成潜艇网格后迭代一次后就不收敛，潜艇网格通过了checkmesh，第一个时间步长库郎数为0.0003，但第二个时间步长库郎数激增到95，具体质量和计算结果具体如下，想求助网格质量需要多高呢，最好有什么办法可以在现有网格基础上改进计算结果，困扰很久的问题了，谢谢各位大佬!

[李东岳]

这种大网格算例都不太好debug，别人帮不上忙，跑一次要好久

[Miraitowa]

@李东岳 老师您好，我是在icem中画的网格导入的，在icem中检查Determinant 3x3x3质量时最小为0.4,但是of中使用checkMesh -allGeometry -allTopology -writeAllFields -writeSets vtk命令时，很多网格的Determinant 小于0.001，想请问一下是什么原因呢，网格质量一般最低要求是什么呢，谢谢老师

[李东岳]

不看Determinant，你的Quality是多少

[Miraitowa]

@李东岳 Quality是0.354，用pimplefoam时可以收敛，但是buoyantfoam就会发散，老师请问使用openfoam计算的话在icem中提高网格质量主要是通过参考Quality吗，我一直是在提高Orthogonal Quality，现在Orthogonal Quality是0.551

[李东岳]

Quality是0.354

网格质量比较差。一般是看综合的Quality

[Miraitowa]

@李东岳 老师我把Quality调到了0.602，但是其他质量有所下降，计算结果还是发散，想请问一下一般Quality需要调整到多少呢

[李东岳]

0.6可以了。还发散需要debug算例了。

[Miraitowa]

@李东岳 但是老师其他完全相同的设置我把网格换成圆柱绕流的0.7Quality可以正常跑出来结果，用这个网格的话计算三四步就会发散，这种情况也是需要debug算例嘛

[李东岳]

是的，越复杂的求解器越容易发散。

[Miraitowa]

@李东岳 好的谢谢老师，我之前一直以为是网格问题在修改网格

[Afr1yne]

@Miraitowa 检查下p_rgh的边界条件，如果p_rgh用的GAMG可以改成PCG试试

[Miraitowa]

@Afr1yne 您好，谢谢您，想请您看一下这样设置有问题吗，这个是p的设置：

internalField   uniform 101325;

boundaryField
{
    WLY_IN
    {
        type            calculated;
        value           uniform 101325;
    }
    WLY_OUT
    {
        type            calculated;
        value           uniform 101325;
    }
    WLY_CE
    {
        type            calculated;
        value           uniform 101325;
    }
    QIANTING
    {
        type            calculated;
        value           uniform 101325;
    }
}

这个是prgh的设置

internalField   uniform 0;

boundaryField
{
    WLY_IN
    {
        type            fixedFluxPressure;
        value           $internalField;
    }

    WLY_OUT
    {
        type            fixedValue;
        value           uniform 0;
    }

    QIANTING
    {
        type            fixedFluxPressure;
        value           $internalField;
    }
    WLY_CE
    {
        type            fixedFluxPressure;
        value           $internalField;
    }

}

[ShaneHEEE]

有米有试试p_rgh除了出口fixedvalue，其他都给zeroGradient

[Miraitowa]

@ShaneHEEE 试过的，一样也是发散的，就很怪

[Afr1yne]

@Miraitowa 我看你说的是能算几步才发散是吧，可以把图放上来看一下。我之前用buoyantFoam发散是因为尺度问题，你把整个计算域缩小一下试试能不能跑

[Miraitowa]

@Afr1yne 您好，这个是计算域的图：

大小是19.5m9.6m9m的，下面是计算两个时间步后的结果：

我仔细看过之后发现速度是从出口开始发散，明显大于理论值，这个是残差图：

[Afr1yne]

@Miraitowa 把出口也设成fixedFluxPressure; 然后设置个pref试下

[Miraitowa]

@Afr1yne 好的谢谢，我试一下

[Miraitowa]

@Afr1yne 您好，我试过了，还是发散的，但是变成入口处速度明显大了