Skip to content
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]
皮肤
  • 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

  • 默认(不使用皮肤)
  • 不使用皮肤
折叠
CFD中文网

CFD中文网

  1. CFD中文网
  2. OpenFOAM
  3. 搅拌桨叶附近压力异常

搅拌桨叶附近压力异常

已定时 已固定 已锁定 已移动 OpenFOAM
21 帖子 5 发布者 12.1k 浏览
  • 从旧到新
  • 从新到旧
  • 最多赞同
回复
  • 在新帖中回复
登录后回复
此主题已被删除。只有拥有主题管理权限的用户可以查看。
  • 青 离线
    青 离线
    青山
    在 中回复了 tidedrinker 最后由 编辑
    #4

    @tidedrinker 你好,搅拌是可以用simpleFoam模拟稳态的,of9自带的simleFoam案例中就有搅拌案例,另外我用multiphaseEulerFoam也算了瞬态情况,压力出现同样问题。

    tidedrinkerT 1 条回复 最后回复
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    写于 最后由 编辑
    #5

    是加的多重参考系么?

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

    青 1 条回复 最后回复
  • tidedrinkerT 离线
    tidedrinkerT 离线
    tidedrinker
    在 中回复了 青山 最后由 编辑
    #6

    @青山 我去看了一下那个二维的搅拌器,我手头没有装好的OF9,所以看的是OF10,应该是一样的。但是有几个问题:

    第一个,你的算例中出现的这种现象,在官方算例里同样出现了,只是没有你的那么明显。

    第二个,官方算例里是二维的,你这里是三维的,流动复杂程度天差地别,能不能简化成稳态计算很成问题。比如说,搅拌器运行的时候,流体会同时从叶片上下和尖端流过,而且可能形成复杂的非定常的脱落涡。

    第三个,官方算例里是用的MRF,多重参考系,而且划分了不同的计算域。但是在你给的网格图里,我不太能想象你是怎么划分网格区域的,或者说,你的MRF部分是怎么设置的?

    第四个,也是我认为可能影响最大的,你的网格太粗了。

    这四个原因都有可能,不太能确定那个影响最大。另外,如果要做三维搅拌器,

    更高的主频,更多的核心

    青 1 条回复 最后回复
  • 青 离线
    青 离线
    青山
    在 中回复了 李东岳 最后由 编辑
    #7

    @李东岳 是的东岳老师,用的MRF方法

    1 条回复 最后回复
  • 青 离线
    青 离线
    青山
    在 中回复了 tidedrinker 最后由 编辑
    #8

    @tidedrinker 感谢你的认真解答,首先网格粗细问题的话,我本来是在精细网格下计算多相流时发现的这个问题,然后用粗网格再次验证了下,最开始的网格如下:9a8f2d8b-93f8-427a-9ae0-b782917fa167-image.png
    压力也出现了相同的现象508c517b-9781-4e6a-9575-b7379b6e399e-image.png
    另外MRF区域使用topoSet进行设置,区域直径约为1.2倍的桨叶直径和高度。
    兄弟你最后一句话是啥,没写完。

    1 条回复 最后回复
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    写于 最后由 编辑
    #9

    把你的云图,不要显示插值场,显示间断场看一下。

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

    青 1 条回复 最后回复
  • 青 离线
    青 离线
    青山
    在 中回复了 李东岳 最后由 编辑
    #10

    @李东岳 东岳老师,这是用multiphaseEulerFoam算液固两相的压力场和速度场:09cac503-d901-4395-a2f0-c9540a28ac33-image.png
    82d126ce-9f96-4857-8894-9494e1c2b1f2-image.png
    这是simpleFoam算的压力场和速度场:f82d666b-2ba5-42d5-b91f-03a8aa24ba30-image.png
    14629460-0d31-4367-bed8-6521dde35d21-image.png

    1 条回复 最后回复
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    写于 最后由 编辑
    #11

    好像没转起来,MRF转的方向对么

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

    青 1 条回复 最后回复
  • 青 离线
    青 离线
    青山
    在 中回复了 李东岳 最后由 编辑
    #12

    @李东岳 老师,MRF区域应该没啥问题,沿Y轴旋转,650rmp。8e9a8d80-2f54-4d1b-b6f3-37f0d071d10d-image.png
    然后速度场看上去也是在转的:20cb61fd-47da-449b-bc33-cbcfbfbd63b0-image.png

    1 条回复 最后回复
  • 知 离线
    知 离线
    知行
    写于 最后由 编辑
    #13

    1.首先感觉你的网格正交性不是太好,请查验
    2.请告知你的网格物理尺寸
    3.是否做了网格无关性呢
    4.感觉大概率是网格数不足,及网格质量不好的原因
    5.请参考这篇文章 https://journals.sagepub.com/doi/pdf/10.1260/1757-482X.3.1.13
    6.条件许可的话网格建议上千万,湍流方法建议用LES
    7.所有壁面处建议都加密并注意查看y+

    收敛收敛

    青 1 条回复 最后回复
  • 青 离线
    青 离线
    青山
    在 中回复了 知行 最后由 编辑
    #14

    @知行 感谢耐心解答!首先网格的质量如下:

    Mesh stats
        points:           182724
        faces:            532288
        internal faces:   516536
        cells:            174804
        faces per cell:   6
        boundary patches: 3
        point zones:      0
        face zones:       1
        cell zones:       3
    
    Overall number of cells of each type:
        hexahedra:     174804
        prisms:        0
        wedges:        0
        pyramids:      0
        tet wedges:    0
        tetrahedra:    0
        polyhedra:     0
    
    Checking topology...
        Boundary definition OK.
        Cell to face addressing OK.
        Point usage OK.
        Upper triangular ordering OK.
        Face vertices OK.
        Number of regions: 1 (OK).
    
    Checking patch topology for multiply connected surfaces...
        Patch               Faces    Points   Surface topology
        walls               10840    10915    ok (non-closed singly connected)
        impeller1           1324     1326     ok (closed singly connected)
        outlet              3588     3663     ok (non-closed singly connected)
    
    Checking geometry...
        Overall domain bounding box (-0.109995 -1.11019e-16 -0.109997) (0.109995 0.22 0.109997)
        Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
        Mesh has 3 solution (non-empty) directions (1 1 1)
        Boundary openness (5.83486e-16 -2.30396e-16 4.46812e-17) OK.
        Max cell openness = 3.29007e-16 OK.
        Max aspect ratio = 5.02531 OK.
        Minimum face area = 1e-06. Maximum face area = 2.73916e-05.  Face area magnitudes OK.
        Min volume = 4.64288e-09. Max volume = 1.25165e-07.  Total volume = 0.0083523.  Cell volumes OK.
        Mesh non-orthogonality Max: 39.7405 average: 7.14766
        Non-orthogonality check OK.
        Face pyramids OK.
        Max skewness = 0.506949 OK.
        Coupled point location match (average 0) OK.
    
    Mesh OK.
    
    End
    

    最大非正交性为39,平均为7,正交性还不错。
    其次,几何是一个直径和高都是22cm的釜,桨叶直径8cm。
    关于网格数量的话,由于目前还在案例调试阶段,所以先用17w的网格算起来比较快,五六十万的网格我也算过,也会出现这种现象。LES的计算代价有点大,所以还是想考虑RANS。
    最后,在所有壁面处我使用了湍流模型kEpsilon的壁面函数,所以是不是对壁面处的网格厚度依赖性不是很大?我尝试把壁面加密一下试试。谢谢!

    知 2 条回复 最后回复
  • 知 离线
    知 离线
    知行
    在 中回复了 青山 最后由 编辑
    #15

    @青山 客气,都是几点非常微小的建议哈。
    这篇文章可能和你的算例接近,doi:10.1017/flo.2022.27
    你的网格量感觉不是太容易出较为理想的效果,不连续的地方正是网格正交性不好的地方。你这个网格,即使加密到1000万,正交性不好的地方仍然会出现奇观的点(因为我算过),建议你用CFX等其他工具画网格。浅见。

    收敛收敛

    知 1 条回复 最后回复
  • 知 离线
    知 离线
    知行
    在 中回复了 知行 最后由 编辑
    #16

    @知行 说错了,是ICEM

    收敛收敛

    1 条回复 最后回复
  • 知 离线
    知 离线
    知行
    在 中回复了 青山 最后由 编辑
    #17

    @青山 壁面和搅拌浆区域必须加密

    收敛收敛

    青 1 条回复 最后回复
  • 青 离线
    青 离线
    青山
    在 中回复了 知行 最后由 编辑
    #18

    @知行 ok,我的网格都是ICEM画的,我加密试试,谢谢

    1 条回复 最后回复
  • 1 离线
    1 离线
    15735420621
    写于 最后由 编辑
    #19

    您好,这个旋转的文件是在哪里设置的呢,目前正在做旋转喷射的案例,但是一直不知道怎么设置

    青 1 条回复 最后回复
  • 青 离线
    青 离线
    青山
    在 中回复了 15735420621 最后由 编辑
    #20

    @15735420621 可以参考自带案例mixerVessel2D,旋转域使用topoSet进行设置,旋转参数在MRFProperties中进行设置

    1 1 条回复 最后回复
  • 1 离线
    1 离线
    15735420621
    在 中回复了 青山 最后由 编辑
    #21

    @青山 好的谢谢!

    1 条回复 最后回复

  • 登录

  • 登录或注册以进行搜索。
  • 第一个帖子
    最后一个帖子
0
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]