Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

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

    OpenFOAM
    2
    2
    547
    Loading More Posts
    • Oldest to Newest
    • Newest to Oldest
    • Most Votes
    Reply
    • Reply as topic
    Log in to reply
    This topic has been deleted. Only users with topic management privileges can see it.
    • K
      kennyadapter last edited by

      各位老师们好,我最近使用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 Reply Last reply Reply Quote
      • 李东岳
        李东岳 管理员 last edited by

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

        CFD高性能服务器 http://dyfluid.com/servers.html

        1 Reply Last reply Reply Quote
        • First post
          Last post

        CFD中文网 | 东岳流体 | 京ICP备15017992号-2
        论坛登录问题反馈可联系 li.dy@dyfluid.com