雾化模拟（reactingFoam）在LES条件下细化网格后遇到的的报错

[mingyang]

我的计算域是50×50×75mm，使用的是OpenFOAM-10版本下的reactingFoam求解器。之前是使用的是OpenFOAM-8版本下的sprayFoam求解器，好像在OpenFOAM-10版本下，sprayFoam求解器合并到reactingFoam里面了，所以目前使用reactingFoam求解器。一开始使用blockMesh {hex (0 1 2 3 4 5 6 7) (100 100 150) simpleGrading (1 1 1)}，这样最小网格是0.5mm，此时还可以正常计算。因为网格优点粗糙，我看有的文献网格在0.15mm以下。接着我尝试细化网格，下面是我的网格情况：


可是这时候刚开始计算就会报以下错误：

中间区域的网格尺寸应该是0.2mm。没有细化网格之前的的时间步是1e-06，细化网格后我改成1e-07.这是基于单核计算的结果。喷雾速度大小为250m/s。
湍流模型的设置如下：

LES
{
    model          kEqn;

    delta          cubeRootVol;

    turbulence      on;

    printCoeffs     on;


    kEqnCoeffs
    {
        ck      0.07;
    }

    cubeRootVolCoeffs
    {
        deltaCoeff      1;
    }

    PrandtlCoeffs
    {
        delta           cubeRootVol;
        cubeRootVolCoeffs
        {
            deltaCoeff      1;
        }
        smoothCoeffs
        {
            delta           cubeRootVol;
            cubeRootVolCoeffs
            {
                deltaCoeff      1;
            }
            maxDeltaRatio   1.1;
        }
        Cdelta           0.158;
    }

    vanDriestCoeffs
    {
        delta           cubeRootVol;
        cubeRootVolCoeffs
        {
            deltaCoeff      1;
        }
        smoothCoeffs
        {
            delta           cubeRootVol;
            cubeRootVolCoeffs
            {
                deltaCoeff      1;
            }
            maxDeltaRatio   1.1;
        }
        Aplus            26;
        Cdelta           0.158;
    }

    smoothCoeffs
    {
        delta           cubeRootVol;
        cubeRootVolCoeffs
        {
            deltaCoeff      1;
        }
        maxDeltaRatio   1.1;
    }
}

麻烦有了解LES方面的大佬帮忙看看，感谢！感谢！感谢！（重要的事情说三遍！）

[李东岳]

从log来看，你这个就是发散了，判断不出其他原因

[mingyang]

@李东岳 谢谢李老师的回答，我再找找原因！