Skip to content
  • 最新
  • Categories
  • 东岳流体
  • 随机看[请狂点我]
Skins
  • Light
  • Cerulean
  • Cosmo
  • Flatly
  • Journal
  • Litera
  • Lumen
  • Lux
  • Materia
  • Minty
  • Morph
  • Pulse
  • Sandstone
  • Simplex
  • Sketchy
  • Spacelab
  • United
  • Yeti
  • Zephyr
  • Dark
  • Cyborg
  • Darkly
  • Quartz
  • Slate
  • Solar
  • Superhero
  • Vapor

  • Default (No Skin)
  • No Skin
Collapse
CFD中文网

CFD中文网

  1. CFD中文网
  2. OpenFOAM
  3. 关于结构入水的问题有些疑问

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

Scheduled Pinned Locked Moved OpenFOAM
3 Posts 3 Posters 3.0k Views
  • 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 Offline
    K Offline
    kennyadapter
    wrote on last edited by
    #1

    各位老师们好,我最近使用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
  • 李东岳李 Online
    李东岳李 Online
    李东岳 管理员
    wrote on last edited by
    #2

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

    9月CFD算法编程课: http://dyfluid.com/class.html

    需要帮助debug算例的看这个 https://cfd-china.com/topic/8018

    Z 1 Reply Last reply
  • Z Offline
    Z Offline
    z597288
    replied to 李东岳 on last edited by
    #3

    @李东岳 老师我想请问一下,设置了重力加速度g后,我将dynamicDict文件下将 velcity设为(0 0 0),物体就停着不动了,请问需要怎么设置一下呢。

    1 Reply Last reply

  • Login

  • Login or register to search.
  • First post
    Last post
0
  • 最新
  • Categories
  • 东岳流体
  • 随机看[请狂点我]