搅拌桨叶附近压力异常
-
fvoption和fvschemes附上：
/*--------------------------------*- C++ -*----------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: 9 \\/ M anipulation | \*---------------------------------------------------------------------------*/ FoamFile { format ascii; class dictionary; location "system"; object fvSolution; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // solvers { p { solver GAMG; tolerance 1e-08; relTol 0.05; smoother GaussSeidel; nCellsInCoarsestLevel 20; } U { solver smoothSolver; smoother GaussSeidel; nSweeps 2; tolerance 1e-07; relTol 0.1; } k { solver smoothSolver; smoother GaussSeidel; nSweeps 2; tolerance 1e-07; relTol 0.1; } epsilon { solver smoothSolver; smoother GaussSeidel; nSweeps 2; tolerance 1e-07; relTol 0.1; } } SIMPLE { nNonOrthogonalCorrectors 1; pRefCell 0; pRefValue 0; } relaxationFactors { fields { p 0.3; } equations { U 0.5; k 0.5; epsilon 0.5; } }/*--------------------------------*- C++ -*----------------------------------*\ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org \\ / A nd | Version: 9 \\/ M anipulation | \*---------------------------------------------------------------------------*/ FoamFile { format ascii; class dictionary; location "system"; object fvSchemes; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // ddtSchemes { default steadyState; } gradSchemes { default Gauss linear; } divSchemes { default none; div(phi,U) bounded Gauss limitedLinearV 1; div(phi,k) bounded Gauss limitedLinear 1; div(phi,epsilon) bounded Gauss limitedLinear 1; div((nuEff*dev2(T(grad(U))))) Gauss linear; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; } snGradSchemes { default corrected; } -
我有个问题哈，就是simpleFoam一般是算稳态的？fvSchemes里ddt也是steadyState，但是搅拌器可以用稳态算吗？
-
@tidedrinker 你好，搅拌是可以用simpleFoam模拟稳态的，of9自带的simleFoam案例中就有搅拌案例，另外我用multiphaseEulerFoam也算了瞬态情况，压力出现同样问题。
-
@青山 我去看了一下那个二维的搅拌器，我手头没有装好的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中进行设置
