CFD中文网

@李东岳
老师好，用的interFoam和twoLiquidMixingFoam，是VOF？

想设置一组边界条件，就是两相流压力驱动流，我设置的边界条件总是导致发散，想请教一下有没有稳定的边界条件？
这个边界条件算单向流的时候好像没有问题，怎么到两相流就不行了呢？

Inlet：
U:   type    zeroGradient;
p:   type    fixedValue;
     value   uniform 10;

Outlet：
U:   type    zeroGradient;
p:   type    fixedValue;
     value  uniform 0;

@OItoCFD 非常感谢，我们也准备加密网格，但是具体加密的程度还不知道怎么把握，您这边是使所有网格大小尽量均一吗？

@李东岳

非常感谢老师，虽然库朗数还是很小，但是其实已经发散了是吗？

我用的是层流求解，用的还是upwind格式，下一步是不是只能降低进口流速了？

@李东岳
老师好，是我没有表达清楚，是类似于这样
库朗数和deltaT都非常非常小，计算非常慢，到最后deltaT基本变成0了
我设置的最大库朗数是0.9，但是计算出来才0.0012，不应该自动缩小时间步长呀

请问老师，这个应该怎么解决，如何提高计算速度？

Courant Number mean: 3.58164e-05 max: 0.901511
Interface Courant Number mean: 2.94865e-08 max: 0.00119605
deltaT = 1.7962e-05
Time = 2.01626

MULES: Solving for alpha.Sarin
Phase-1 volume fraction = 6.16252e-05  Min(alpha.Sarin) = 0  Max(alpha.Sarin) = 0.109368
MULES: Solving for alpha.Sarin
Phase-1 volume fraction = 6.16252e-05  Min(alpha.Sarin) = 0  Max(alpha.Sarin) = 0.109368
smoothSolver:  Solving for alpha.Sarin, Initial residual = 5.46198e-10, Final residual = 5.46198e-10, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 0.00020754, Final residual = 9.85112e-08, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.000110383, Final residual = 4.3286e-08, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 0.000191943, Final residual = 1.79758e-06, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.00248148, Final residual = 1.41595e-05, No Iterations 6
time step continuity errors : sum local = 3.84475e-13, global = 1.89675e-14, cumulative = -3.04519e-09
GAMG:  Solving for p_rgh, Initial residual = 0.00152746, Final residual = 5.47895e-06, No Iterations 7
time step continuity errors : sum local = 1.48826e-13, global = -6.48062e-15, cumulative = -3.04519e-09
ExecutionTime = 680.02 s  ClockTime = 680 s

Courant Number mean: 3.57663e-05 max: 0.90151
Interface Courant Number mean: 2.94446e-08 max: 0.00119381
deltaT = 1.79319e-05
Time = 2.01628

不好意思老师，是我说的不清楚，发现一直在计算同一个时间，比如一直计算0.1111秒

各位老师好：

写了一个非常简单的pisoFoamTemp，测试了一下发现一切正常，用elbow.msh测试的

然后在实际网格中，发现根本就没有计算出来，一直在入口处，都没有进入计算域，但是速度压力求解都正常，请问这是什么原因？
不好意思图片不是很清晰，但是可以看到温度就一直在壁面处没有变化，但是速度计算看上去比较正常，这个是计算150s的结果

附上算例文件
链接: https://pan.baidu.com/s/1RaX15gCuychiMvbJBc9MMg 密码: 8gha

各位老师好：

用自动调节时间步的时候经常发现一直在循环计算某一个时间，这是什么原因，是发散了吗？
因为一开始测试为了简单和利于收敛，设置的层流条件，实际上有部分地方应该是湍流，这个会影响吗？

如果降低进口速度，deltaT，有助于改善吗？

各位老师好：

运行decomposePar的时候显示如下错误：

decomposePar: error while loading shared libraries: libscotch.so: cannot open shared object file: No such file or directory

此外，如果是单核运算的话，是没问题的。

请问各位老师，这是什么问题，怎样解决啊？

@bestucan 非常感谢老师的提醒，让网格变稀疏之后发现了两处可疑的质量较差的地方。尝试之后发现，有一些长方体网格划的太瘦了，ICEM的网格质量检查没有问题，修改之后OpenFOAM中就也没有问题了。

@bestucan 非常感谢老师的指点，非常明白，我会再检查一下模型转换的问题。

然后还是想了解一下，请问有没有什么能够快速 定位这些有问题的网格单元的位置 的方法呢？

@bestucan 老师您好，非常感谢老师您的回答。经过几天的检查发现，没封闭的原因在于把stl文件导入fluent的时候有的线条和面的格式出现了问题。

但是第一个网格质量较差的原因却没有查找出来，按照fluent对于skewness的定义，skewness应该不能超过1才对，而且我在icem中的网格质量最低是0.6，已经达到较高的标准了。

翻了翻资料说OpenFOAM检查的网格单元和 面元对应顶点的坐标、两个网格单元中心连线向量、面元中心与网格单元中心连线向量、面元的面法相向量 相关，但是没有找到具体的算法

请问老师这个问题怎么排查？或者说怎么才能快速的找到有质量问题的网格？

各位老师好，大家好：
使用ICEM画了网格之后checkMesh时显示了好几个问题，但是我不是很能理解是什么意思，也不太清楚往什么方向检查或者修改。但是，我之前在icem里面网格质量检查都显示所有网格都有较高的质量。。。

第一个网格问题：
说有1296个高扭曲的面，但是同时我发现，这一点并不影响计算，请问这个问题是说有一些网格质量太差吗？是OpenFOAM对于网格质量的判定标准和icem不一样？

Checking geometry...
    Overall domain bounding box (-72.0987 -6.78501 -2.15365) (59.2013 6.71664 3.24987)
    ......
    Face pyramids OK.
 ***Max skewness = 4.54714, 1269 highly skew faces detected which may impair the quality of the results
  <<Writing 1269 skew faces to set skewFaces
    Coupled point location match (average 0) OK.

Failed 1 mesh checks.

第二个问题：
这个就有些严重了，直接导致计算无法进行，请问这个是说边界没有封闭吗:

Checking geometry...
    Overall domain bounding box (-89.7486 -38.8002 1.14635) (-27.4486 24.9498 21.3019)
    Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
    Mesh has 3 solution (non-empty) directions (1 1 1)
 ***Boundary openness (0.0173317 5.152e-10 -7.40693e-05) possible hole in boundary description.
 ***Open cells found, max cell openness: 1, number of open cells 2112
  <<Writing 2112 non closed cells to set nonClosedCells
   Minimum face area = 0.00318294. Maximum face area = 0.230028.  Face area magnitudes OK.
    Min volume = 0.000594917. Max volume = 0.0453085.  Total volume = 7140.58.  Cell volumes OK.
     ......
   
Failed 2 mesh checks.

还请各位老师不吝赐教，非常感谢！

@fangyuanaza 非常感谢！

@东岳 老师好，非常感谢您！
请问是要在输出网格的时候存成ASCII码格式吗？如果只有.msh文件可以用ICEM转换成ascii吗？

东岳老师说过，fluent3DMeshToFoam对于转换3D网格更好一些，对于我自己使用ICEM画的网格，都很好，这还是我第一次遇到转换失败的问题
另外，这个网格文件打开都是乱码：

(3010 (e 1 9211 2 3)
(   @J窯?    ?@   鄒W?   棱G?   @ef@   纔

是因为比较大的原因吗？
我自己画的网格文件打开的话都是一行行数字啊

大家好，各位老师好：

最近网上下载了一个网格，在用ICEM观察时没有发现错误，
但是使用fluent3DMeshToFoam转换网格时报错：

--> FOAM FATAL ERROR: 
Do not understand characters: 
    on line 19

    From function virtual int yyFlexLexer::yylex()
    in file fluent3DMeshToFoam.L at line 748.

FOAM exiting

请问这会是什么原因？怎样才能解决呢？

请问会是用gambit生成和ICEM生成导致的吗？

另外，如果改成fluentMeshToFoam转换的话，还是会出现乱码：

Found end of section in unknown:)
Embedded blocks in comment or unknown:
(
@�@��@@��@@��������������
@��Found end of section in unknown:)
�@@��@@�����������@��@@�`Found end of section in unknown:#
(����@@�Embedded blocks in comment or unknown:�>
��@@�Found end of section in unknown:>
��@@���@@����@@@����@@����@@��@@�}�@@����@@�Embedded blocks in comment or unknown:�>
���@@�Embedded blocks in comment or unknown:��

哎呀，名额满了吗？
没满加一个？
18994065595@163.com

希望对大家有帮助。
实时更新sci-hub

发帖混个积分

@东岳 非常感谢老师 ！
就是这样迂回的话好像有些难，我还要再思考思考，再尝试尝试