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



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

    与论文的速度对比图:9bbdc124-23a5-438c-adef-9c534082b8fe-image.png
    我尝试使用blockMesh自己画楔形结果也是在某一处发散。
    两个楔形网格如图所示:
    2a0cc4e0-d3fe-413f-907d-79c4ef192d34-image.png

    
    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
    

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



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


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2