foam-extend 4.0 浸没边界法 interIbFoam 边界条件报错

[ThomasShi]

我用interIbFoam求解有结构物的明渠流问题，在应用了k-epsilon湍流模型后，nut文件的浸没边界条件出现如下报错：

--> FOAM FATAL ERROR:

    evaluate() cannot be called for a genericFvPatchField (actual type immersedBoundaryWallFunction)
    on patch vegetation of field nut in file "/mnt/g/foam/thomaschi-4.0/run/cylBumpInterIbFoam/0/nut"
    You are probably trying to solve for a field with a generic boundary condition.

    From function genericFvPatchField<Type>::evaluate(const Pstream::commsTypes)
    in file fields/fvPatchFields/basic/generic/genericFvPatchField.C at line 760.

FOAM exiting

有人知道是什么原因吗？该怎么解决？

[bestucan]

@ThomasShi 在 foam-extend 4.0 浸没边界法 interIbFoam 边界条件报错 中说：

evaluate() cannot be called for a genericFvPatchField

这个函数是你自己定义的么？

[ThomasShi]

@bestucan 您好，感谢回帖。那个并不是我自己定义的。我只是参照simpleIbFoam的湍流算例里的边界条件文件，给我的浸没边界加了同样的边界条件。因为interIbFoam只有层流的算例没有湍流的，所以我只好照着别的类似的来做的，结果出现了边界条件的报错。nut文件浸没边界的边界条件设置如下

immersedBoundary
{
        type            immersedBoundaryWallFunction;
        patchType       immersedBoundary;
        refValue        uniform 1e-10;
        refGradient     uniform 0;
        fixesValue      false;
        setDeadCellValue yes;
        deadCellValue   1e-10;
        value           nonuniform 0();
}

[bestucan]

@ThomasShi 这个里有interIbFoam使用湍流的讲解，2.3.3。在里面搜immersedBoundaryWall，也可以搜到相关设置。
http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2016/MohsenIrannezhad/Final_Report.pdf

这个错误和你的很像，只是函数不一样，兴许他用的不是k~e
https://www.cfd-online.com/Forums/openfoam-solving/231984-immersed-boundary-method-error-turbulence-model.html
边界条件没对上，这个讲边界条件很清楚：
https://technodocbox.com/3D_Graphics/67910317-Immersed-boundary-method-in-foam.html

[ThomasShi]

@bestucan 万分感谢！