CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

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

    OpenFOAM
    2
    2
    701
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • K
      kennyadapter 最后由 编辑

      各位老师们好,我最近使用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
      

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

      1 条回复 最后回复 回复 引用
      • 李东岳
        李东岳 管理员 最后由 编辑

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

        CFD高性能服务器 http://dyfluid.com/servers.html
        2023年,线下CFD课,预热一下 http://dyfluid.com/class.html

        1 条回复 最后回复 回复 引用
        • First post
          Last post