DNS 计算报错

冠竹

各位大佬，我现在正在用DNS求解器求解槽道湍流。
速度和压力设置如下：

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

internalField   uniform (0 0 0);

boundaryField
{
    INLET
    {
        type            fixedValue;
        value           uniform (0.1616 0 0);
    }
    OUTLET
    {
        type            zeroGradient;
    }

    wall_r
    {
        type            cyclic;
    }

   GRID
    {
        type            zeroGradient;
    }

    wall_l
    {
        type            cyclic;
    }

    wall_u
    {
        type            cyclic;
    }

    wall_d
    {
        type            cyclic;
    }

压力：

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    INLET
    {
        type            zeroGradient;
    }

    OUTLET
    {
        type            zeroGradient;
    }

    wall_r
    {
        type            cyclic;
    }

   GRID
    {
        type            zeroGradient;
    }

    wall_l
    {
        type            cyclic;
    }

    wall_u
    {
        type            cyclic;
    }

    wall_d
    {
        type            cyclic;
    }
}

由于出口想设置成outflow，故选用zeroGradient条件。
但是当我去求解的时候发现报错提示不满足连续性方程

--> FOAM FATAL ERROR: 
Continuity error cannot be removed by adjusting the outflow.
Please check the velocity boundary conditions and/or run potentialFoam to initialise the outflow.
Total flux              : 0.0390532
Specified mass inflow   : 5.18519e-225
Specified mass outflow  : 1616
Adjustable mass outflow : 8.9386e-209


    From function adjustPhi(surfaceScalarField&, const volVectorField&,volScalarField&)
    in file cfdTools/general/adjustPhi/adjustPhi.C at line 114.

FOAM exiting

但当我使用potentialFoam试图建立一个初始场的时候发现会出现同样的错误
各位大佬，请问一下我这个错误出现的原因是什么呢，应该如何修复

李东岳

@冠竹在 DNS 计算报错中说：

DNS求解器

哪个求解器？

冠竹

@李东岳这个是自己编译的一个求解器，老师，不是OF自带的

李东岳

@冠竹 openfoam里面的DNS求解器？openfoam这面很难做DNS啊，莫非是2阶DNS？

自己写的求解器别人很难debug问题在哪，一般只有自己知道

OItoCFD

@冠竹那问题就不好说了哎没有标准参考自己搞的话一般问题出在自己搞得地方不是因为用的零梯度边界的原因

冠竹

@李东岳这个是二阶DNS，是我师兄当时编译的，他做的是流动加传热。我把师兄的能量方程剔除掉之后计算的。当时他的这个求解器应该是没啥大问题吧压力求解就是用的简单的piso循环，速度就是N-S方程...
不知道为啥会出现这个报错...

冠竹

@oitocfd 大佬，也就是说我的边界这么设置符合求解要求嘛

OItoCFD

@冠竹出流用零梯度肯定没啥问题的呀肯定是你自己搞得这个求解器有地方是错的

冠竹

@oitocfd 感谢大佬，那我再检查监测

冠竹

@oitocfd 大佬，我这两天测试了一下，好像是边界条件的问题...我用这套边界条件，用OpenFOAM自带的DNS求解器、icoFoam、simpleFOAM求解，都是这个报错

OItoCFD

@冠竹那真不清楚了哎很奇怪

李东岳

Continuity error cannot be removed by adjusting the outflow.

这个肯定是边界条件的问题。你这个自己写的求解器？也有可能你的求解器边界条件没处理好。

这个是二阶DNS，是我师兄当时编译的，他做的是流动加传热。我把师兄的能量方程剔除掉之后计算的。当时他的这个求解器应该是没啥大问题吧压力求解就是用的简单的piso循环，速度就是N-S方程...

你为什么不直接用pisoFoam/pimpleFoam呢？这个就能算2阶DNS。这是我之前用pisoFoam算的quasi-DNS

替代文字

冠竹

@李东岳感谢东岳大佬，我去设置一下看看

CFD中文网

DNS 计算报错