用snappyHexMeshDict生成周期性边界问题

莫晓柒

@东岳

--> FOAM FATAL ERROR: 
face 5120 area does not match neighbour by 2.02625% -- possible face ordering problem.
patch:cyclic1 my area:6.14979e-05 neighbour area:6.02643e-05 matching tolerance:0.0001
Mesh face:3995462 fc:(0.0020217 0.0037764 5.81648)
Neighbour fc:(-0.00217946 0.00377493 5.81651)
If you are certain your matching is correct you can increase the 'matchTolerance' setting in the patch dictionary in the boundary file.
Rerun with cyclic debug flag set for more information.

    From function void Foam::cyclicPolyPatch::calcTransforms(const primitivePatch&, const pointField&, const vectorField&, const pointField&, const vectorField&)
    in file meshes/polyMesh/polyPatches/constraint/cyclic/cyclicPolyPatch.C at line 212.

FOAM exiting

莫晓柒

就会提示上述错误

Calf.Z-DNS

@莫晓柒在用snappyHexMeshDict生成周期性边界问题中说：

does not match neighbour

我用fluent3DMeshToFoam转化成OpenFOAM的网格的时候也遇到了does not match neighbour这个问题；我是手动改成cyclic边界的，貌似是要用createPatch的方法？

李东岳

要不用cyclicAMI
要么用手动
要么createPatch，把matchTolerance放的特别大

Calf.Z-DNS

@东岳谢谢东岳老师
我在cfd-online上面找到一种方法：
下面是原文：
“I encountered the same problem and I could solve it by defining rotational periodicity and then by declaring the corresponding vertices in ICEM as periodic. Then Right Click on Faces in the Blocking tree -> Periodic Faces -> ensure that really every periodic face is colored accordingly. createPatch went fine after this procedure.” ---by Fabian E.
然后我勾选了ICEM里面blocking里的faces的Periodic，再导入openfoam，发现就不会出现上述错误。
但我不太清楚为何这样会消除错误？

李东岳

ICEM里面的periodic执行非常严格的标准，它使得block的vertices严格的一一对应，如果能用ICEM做并附加periodic操作，尽可能用ICEM做

jinlinna

@Calf-Z-DNS 在用snappyHexMeshDict生成周期性边界问题中说：

@东岳谢谢东岳老师
我在cfd-online上面找到一种方法：
下面是原文：
“I encountered the same problem and I could solve it by defining rotational periodicity and then by declaring the corresponding vertices in ICEM as periodic. Then Right Click on Faces in the Blocking tree -> Periodic Faces -> ensure that really every periodic face is colored accordingly. createPatch went fine after this procedure.” ---by Fabian E.
然后我勾选了ICEM里面blocking里的faces的Periodic，再导入openfoam，发现就不会出现上述错误。
但我不太清楚为何这样会消除错误？

请问同学可以告知是怎么设置的吗，icem里面有face pericdic吗？我这里只设置了点对称

kuangyonglin

将 matchTolerance 放的特别大能够让仿真跑起来，但是计算结果显示两边周期面上的压力分布不是一样,如图所示是不同时间步，两周期面上的平均压力，其中有一面的平均压力稳定高于另一面的压力

不知道是求解器的问题还是这里设置的问题

有朋友遇见过类似情况吗?

李东岳

对一个网格，如果 matchTolerance 特别小的时候，有这个问题么

kuangyonglin

@李东岳网格是用fluentMeshing画的，严格按照周期性边界设置的，导入OpenFOAM之后，当matchTolerance很小的时候，使用createPatch设置周期性边界的时候会报错，

以下是 createPatch 代码

patches
(
    {
        name cyclic1;
        patchInfo
        {
            type cyclic;
            neighbourPatch cyclic2;
            transform rotational;
            rotationAxis (0 0 1);
            rotationCentre (0 0 0);
            matchTolerance 1e10;
        }
        constructFrom patches;
        patches (cyclic_1);
        set f0;
    }
    {
        name cyclic2;
        patchInfo
        {
            type cyclic;
            neighbourPatch cyclic1;
            transform rotational;
            rotationAxis (0 0 1);
            rotationCentre (0 0 0);
            matchTolerance 1e10;
        }
        constructFrom patches;
        patches (cyclic_2);
        set f0;
    }
);

以下是报错信息

--> FOAM FATAL ERROR: 
More than one patch accessing the same transform but not of the same sign.
patch:cyclic2 transform:0 sign:1  current transforms:1(1)

    From function Foam::label Foam::globalIndexAndTransform::addToTransformIndex(Foam::label, Foam::label, bool, Foam::scalar) const
    in file lnInclude/globalIndexAndTransformI.H at line 210.
FOAM exiting

CFD中文网

用snappyHexMeshDict生成周期性边界问题