搅拌桨叶附近压力异常
-
@青山 我去看了一下那个二维的搅拌器,我手头没有装好的OF9,所以看的是OF10,应该是一样的。但是有几个问题:
第一个,你的算例中出现的这种现象,在官方算例里同样出现了,只是没有你的那么明显。
第二个,官方算例里是二维的,你这里是三维的,流动复杂程度天差地别,能不能简化成稳态计算很成问题。比如说,搅拌器运行的时候,流体会同时从叶片上下和尖端流过,而且可能形成复杂的非定常的脱落涡。
第三个,官方算例里是用的MRF,多重参考系,而且划分了不同的计算域。但是在你给的网格图里,我不太能想象你是怎么划分网格区域的,或者说,你的MRF部分是怎么设置的?
第四个,也是我认为可能影响最大的,你的网格太粗了。
这四个原因都有可能,不太能确定那个影响最大。另外,如果要做三维搅拌器,
-
@tidedrinker 感谢你的认真解答,首先网格粗细问题的话,我本来是在精细网格下计算多相流时发现的这个问题,然后用粗网格再次验证了下,最开始的网格如下:
压力也出现了相同的现象
另外MRF区域使用topoSet进行设置,区域直径约为1.2倍的桨叶直径和高度。
兄弟你最后一句话是啥,没写完。 -
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的壁面函数,所以是不是对壁面处的网格厚度依赖性不是很大?我尝试把壁面加密一下试试。谢谢! -
您好,这个旋转的文件是在哪里设置的呢,目前正在做旋转喷射的案例,但是一直不知道怎么设置
-
@15735420621 可以参考自带案例mixerVessel2D,旋转域使用topoSet进行设置,旋转参数在MRFProperties中进行设置