OF计算动网格的问题

WYing

各位老师好，我在OF7中使用 dynamicMotionSolverFvMesh 和 sixDoFRigidBodyMotion 计算三维圆柱的涡激振动问题。圆柱直径D=1m，预计振幅0.6D左右，遇到问题如下：

圆柱振动到0.2D左右时，算例会突然崩掉，报错：
此时CFLmax=0.7，流场也比较正常，应该是还没来得及崩溃。
通过checkMesh发现计算域底面有几个网格（标红）的skewness非常大，这应该是导致失败的原因：

此时圆柱是向+y轴运动，不太明白为什么靠近圆柱（-y轴方向）的这几层网格会凸起。因为这个范围属于设置的inner Distance之内，按理说网格是不变形的，而是直接随圆柱运动。(目前设置是 inner Distance = 1.5m， outer Distance = 2.5m)

请问有人遇到过类似问题吗？第一次遇到这个情况。。。怎么调都没用

李东岳

@wying 你设置contrain以及restrain了么？这个可以不让他偏移

WYing

@李东岳东岳老师，我在dynamicMeshDict里设置了constrints和restraints两个条目，控制圆柱仅在y方向移动，不转动，如下：

constraints
{
    yLine
    {
        sixDoFRigidBodyMotionConstraint line;
        centreOfRotation    (0 0 0.5);     
        direction           (0 1 0);
    }

    noRotation
    {
        sixDoFRigidBodyMotionConstraint orientation;
    }
}

restraints
{
    verticalSpring
    {
        sixDoFRigidBodyMotionRestraint linearSpring;

        anchor          (0 0 0.5);
        refAttachmentPt (0 0 0.5);
        stiffness       2.1055;
        damping         0;
        restLength      0;
    }
}

不知道您说的是不是这个意思？

CFD中文网

OF计算动网格的问题