OpenFOAM重叠网格

Gentle

@stergopilot 您好，有一些重叠网格的问题想请教一下，能不能加下我qq或者微信呀541658583，谢谢！

lv1995

@wwzhao 你好，请问mergeMesh具体怎么用呢？需要添加那些信息或者有没有算例啊？

wwzhao

@lv1995 OpenFOAM v1706之后的版本提供了重叠网格的示例算例，里面有mergeMesh的用法。

kennyadapter

@wwzhao 您好，我在用overset网格做仿真，在背景网格通过边界条件生成波浪，但是当波浪到子网格中的时候就会出现锯齿形状甚至出现很乱的流动，明明波浪还未碰到物体。请问是不是overset本身的问题，还是我网格没配置好？背景网格和子网格的网格密度要什么关系呢？

谢谢您的解答。

wwzhao

@kennyadapter 应该是算例设置问题。能否看下你的网格、速度和压力场？

kennyadapter

@wwzhao 非常感谢您的回复。
我的两套网格是都是通过blockMesh建立的方形区域，配置文件如下：
背景网格：

vertices
(
    (0 0.09 -1.5)
    (2.5 0.09 -1.5)
    (2.5 0.11 -1.5)
    (0 0.11 -1.5)
    (0 0.09 1)
    (2.5 0.09 1)
    (2.5 0.11 1)
    (0 0.11 1)
);

blocks
(
    hex (0 1 2 3 4 5 6 7) (500 1 500) simpleGrading (1 1 1)
);

edges
(
);

boundary
(
    // Dummy patch to trigger overset interpolation before any other
    // bcs
    oversetPatch
    {
        type overset;
        faces ();
    }
    stationaryWalls
    {
        type wall;
        faces
        (
            (0 3 2 1)
        );
    }

    inlet
    {
        type patch;
        faces
        (
            (0 4 7 3)
        );
    }
    outlet
    {
        type patch;
        faces
        (
     (2 6 5 1)
        );
    }

    frontAndBack
    {
        type empty;
        faces
        (
            (1 5 4 0)
            (3 7 6 2)
         
        );
    }


    atmosphere
    {
        type patch;
        faces
        (
            (4 5 6 7)
        );
    }
);

mergePatchPairs
(
);

子网格：

convertToMeters 1;

vertices
(
    (0.275 0.09 -0.1)
    (0.875 0.09 -0.1)
    (0.875 0.11 -0.1)
    (0.275 0.11 -0.1)
    (0.275 0.09 0.4)
    (0.875 0.09 0.4)
    (0.875 0.11 0.4)
    (0.275 0.11 0.4)
);

blocks
(
    hex (0 1 2 3 4 5 6 7) inletChannel (230 1 100) simpleGrading (1 1 1)
);

edges
(
);

boundary
(
        sides
    {
        type overset;
        faces
        (
            (0 3 2 1)
            (2 6 5 1)         
            (0 4 7 3)
            (4 5 6 7)
        );
    }

    frontAndBack
    {
        type empty;
        faces
        (
            (1 5 4 0)
            (3 7 6 2)
        );
    }

    floatingObject
    {
        type wall;
        faces ();
    }
);

在子网格用snappyHexMesh生成楔形：

将两套网格结合起来：

速度边界条件：

boundaryField
{
    #includeEtc "caseDicts/setConstraintTypes"

    inlet
    {
        type            waveVelocity;
        value           uniform (0 0 0);
    }

    outlet
    {
        type            waveVelocity;
        value           uniform (0 0 0);
    }
    stationaryWalls
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    atmosphere
    {
        type            pressureInletOutletVelocity;
        value           uniform (0 0 0);
    }
    floatingObject
    {
        type            movingWallVelocity;
        value           uniform (0 0 0);
    }
}

压力：

boundaryField
{
    inlet
    {
        type            fixedFluxPressure;
        value           uniform 0;
    }

    outlet
    {
        type            fixedFluxPressure;
        value           uniform 0;
    }

    oversetPatch
    {
        type            overset;
    }
    stationaryWalls
    {
        type            fixedFluxPressure;
    }
    atmosphere
    {
        type            totalPressure;
        p0              uniform 0;
        U               U;
        phi             phi;
        rho             rho;
        psi             none;
        gamma           1;
        value           uniform 0;
    }
    floatingObject
    {
        type            fixedFluxPressure;
    }
    overset
    {
        patchType       overset;
        type            fixedFluxPressure;
    }
}

速度场：

压力场：

kennyadapter

@wwzhao 打扰您了！
alphawater：

kennyadapter

@wwzhao 老师方便的话能加一下QQ讨论么，我的QQ号码：1219854212，非常感谢您

wwzhao

@kennyadapter

我不用QQ。建议你做后处理的时候把洞单元过滤掉。

你可以试下用同样的算例、边界条件，不用重叠网格，看看波浪会不会出现这种情况。

kennyadapter

@wwzhao 谢谢您的回复！

VamosZhi

@kennyadapter 您好，您用的是2D overset网格吗？
是的话，SNAPPY之后不是要extrudemesh一下，这个有什么要注意的问题吗？
因为我merge之后，两套网格不在一个平面上。。

kennyadapter

@VamosZhi 如果是2D的话是要extrudemesh一下，不过要注意下方向和起始面。我一开始也遇到过这个问题，试着改一下

flipNormals true;

这个参数，这个跟方向有关系，如果不行

sourcePatches (front);
exposedPatchName Back;

尝试改一下这两个，，换一下顺序

netbuger

@VamosZhi 请问怎么去掉洞单元和后面的背景单元？paraview做了slice后背景单元的初始自由面总是存在，想把它去掉。

VamosZhi

@netbuger 哥，目前我还没有到后处理的那一步，到了跟你继续交流！

kennyadapter

@netbuger 在clip中clip type选择scalar然后在下面选择CELLtype，，调整为1就行

hongjiewang

@wwzhao 老师您好~我正在做包含多孔介质的传热仿真，不规则图形代表我的多孔介质，红色是网格部分。按照你的理解，我的整体网格是不是可以这样操作，如果我分两个case划分不同区域，是不是snappyHexMesh里的参数设置一样就可以保证他们的界面是完全可以耦合的~感谢老师

DY大世界

@kennyadapter 您好，打扰您了。请问您的问题解决了吗，我现在也需要模拟结构入水的案例，请问您的dynamicMeshDict是怎么写的呢，能否看一下您的dynamicMeshDict字典？

DY大世界

@wwzhao 老师您好，目前我也在用重叠网格模拟固体入水，但是我的水面是静止的没有波浪。我也出现了类似的问题，就是固体还没有接触到水面的时候，水面就已经发生了巨大的扰动，不知道为什么固体还没有接触到水面的时候就会对水体产生这么大的非正常影响。

Jingsheng.Yang

@wwzhao 请问在重叠网格中圆柱周围的贴体网格可以超出背景网格吗？我想做一个近壁圆柱的涡激振动模拟，但是我不太清楚圆柱周围的贴体网格越过背景网格是否合理。能麻烦您帮忙解答一下吗

CFD中文网

OpenFOAM重叠网格