搅拌桨叶附近压力异常

青山

@李东岳东岳老师，这是用multiphaseEulerFoam算液固两相的压力场和速度场：

这是simpleFoam算的压力场和速度场：

李东岳

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

青山

@李东岳老师，MRF区域应该没啥问题，沿Y轴旋转，650rmp。
然后速度场看上去也是在转的：

知行

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

青山

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

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的壁面函数，所以是不是对壁面处的网格厚度依赖性不是很大？我尝试把壁面加密一下试试。谢谢！

知行

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

知行

@知行说错了，是ICEM

知行

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

青山

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

15735420621

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

青山

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

15735420621

@青山好的谢谢！

CFD中文网

搅拌桨叶附近压力异常