关于结构入水的问题有些疑问

kennyadapter

各位老师们好，我最近使用OpenFOAM做楔形自由落体入水，OpenFOAM版本是 1906。我通过修改overInterDyMFoam例子中floatingBody案例将一个方块自由落体进入水中，发现可以计算，然后我将方块通过snappyHexMesh换为楔形，发现在计算到一半就发散了。
发散前后的速度场如图所示：

与论文的速度对比图：
我尝试使用blockMesh自己画楔形结果也是在某一处发散。
两个楔形网格如图所示：


motionSolverLibs    ("libsixDoFRigidBodyMotion.so");

dynamicFvMesh       dynamicOversetFvMesh;

solver          sixDoFRigidBodyMotion;

sixDoFRigidBodyMotionCoeffs
{

    patches         (floatingObject);
    innerDistance   100.0;
    outerDistance   101.0;

    centreOfMass    (3 1.5 1.2);
	momentOfInertia (40 921 921);

    // Density of the solid
    rhoSolid        466.7;

    // Cuboid mass
    mass           94.08;

velocity (0 0 -6.5788);
    report          on;
    accelerationRelaxation 0.6;
    accelerationDamping    0.9;

    solver
    {
        type Newmark;
    }

    constraints
    {

        /*
        fixedPoint
        {
             sixDoFRigidBodyMotionConstraint point;
             centreOfRotation (0.5 0.45 0.1);
        }
       */

        fixedLine
        {
            sixDoFRigidBodyMotionConstraint line;
            //centreOfRotation (0.5 0.45 0.1);
            direction (0 0 1);
        }
         fixedorentation
         {
            sixDoFRigidBodyMotionConstraint orientation;
             centreOfRotation (0.5 0.45 0.1);
         }
/*
        fixedAxis
        {
            sixDoFRigidBodyMotionConstraint axis;
            axis (0 1 0);
        }
        */
    }
}

这是我动网格的代码，也尝试改过

accelerationRelaxation

数值，但是效果也不明显，希望各位老师给出一些建议，非常感谢！

李东岳

为何不固定物体让水面上移？