CFD中文网

想要计算的物理问题是跨临界流体（二氧化碳，或者是氟利昂制冷剂）在大范围压降条件下的相变过程，说人话就是比如入口是超临界的CO2，经过Laval nozzle膨胀到亚临界；
对于这种情况必须考虑真实流体热力学方程了，目前我用的实现方法是通过调用NIST REFPROP生成的real gas property lookup table的办法来计算气液相流体的热物性，这个方法在计算亚临界范围内，大压降的工况也可以计算，收敛情况相对满意，但是一旦计算跨临界状态就发散了，在外国的CFD论坛上找过，类似的答复是NIST流体数据库使用的方程在临界点附近似乎存在不收敛的问题，如图所示:

1c5761fa-95be-44a0-a527-4e2f418b4098-9bacb2a3045cb6950be3e5620394144.png
现在很苦恼，一直卡在这个问题上面，请教一下各位大神，在fluent里面计算跨临界流体热物性还有什么别的实现方法吗

算不准大概率是网格问题，不仅仅是数量，还有质量，另外，系数值的获取是如何计算，自己手动计算还是直接通过软件监测？若通过软件监测，则需注意参考值的设定。

@李东岳
懂了，谢谢李老师！

@李东岳 好的，谢谢李老师耐心解答

周老师，这里我测试发现好像这种结构网格的加密方式对于三维模型可以，对于二维网格，在经过局部加密后再导入到OF里面是会报错的，错误是：--> FOAM FATAL ERROR:
face 11528151 does not belong to cell 210100. Face owner: -1 neighbour: 210099

From Foam::cellShape Foam::extrudedQuadCellShape(Foam::label, const labelList&, const faceList&, const labelList&, const labelList&, Foam::label, Foam::faceList&) in file extrudedQuadCellShape.C at line 105.

FOAM aborting
感觉好像是在z方向拉伸出来的面没有识别到具体的cell中

顺便说一下，liqiwen的上述理解是错误的。网格的空间尺寸和时间步长的要求，由解的波长和频率来决定。没有这样的结论：时间平均后，湍流的空间尺度一定变大。

但是，今天大多数人都结受它。提醒的是，它没有科学的依据。

本人小白，有个疑问：为什么在做流体仿真的时候设置完边界条件后，还要进行初始化后才能开始计算。而做固体有限元分析的时候设置完边界条件后就可以直接开始计算了？照理说微分方程有边界条件后就可以确定唯一的解了，想不明白初始化的目的是什么，各位老师可不可以深入浅出的给我讲解一下？谢谢

@刀尔東 在 向各位大佬请教一下网格划分的基础问题 中说：

感谢各位大佬不吝赐教。

链接: https://pan.baidu.com/s/11z9O_mzDLuqEp30mWRr5Ow?pwd=1234 提取码: 1234 复制这段内容后打开百度网盘手机App，操作更方便哦
--来自百度网盘超级会员v6的分享

以上是我的两个文件

网格文件为ENSIGHT GOLD FORMAT

用fluent Meshing划分的多面体网格转换成OpenFOAM格式出现错误
/---------------------------------------------------------------------------
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: v1912 |
| \ / A nd | Website: www.openfoam.com |
| \/ M anipulation | |
*---------------------------------------------------------------------------*/
Build : _f3950763fe-20191219 OPENFOAM=1912
Arch : "LSB;label=32;scalar=64"
Exec : fluent3DMeshToFoam 000.msh
Date : Jul 03 2020
Time : 09:25:27
Host : King
PID : 514
I/O : uncollated
Case : /mnt/d/OpenFOAM/run/test
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Reading header: "ANSYS(R) TGrid(TM) 3D, parallel 18.1.0"
Dimension of grid: 3
Number of points: 34807
--> FOAM Warning : Found unknown block of type: "11"
on line 16
Number of faces: 39160
Number of cells: 6053
PointGroup: 1 start: 0 end: 15. Reading points...done.
PointGroup: 8 start: 16 end: 23. Reading points...done.
PointGroup: 216 start: 24 end: 307. Reading points...done.
PointGroup: 218 start: 308 end: 3281. Reading points...done.
PointGroup: 59 start: 3282 end: 9500. Reading points...done.
PointGroup: 61 start: 9501 end: 34806. Reading points...done.
--> FOAM Warning : Found unknown block of type: "11"
on line 34838
--> FOAM Warning : Found unknown block of type: "11"
on line 34852
--> FOAM Warning : Found unknown block of type: "11"
on line 34866
--> FOAM Warning : Found unknown block of type: "11"
on line 35156
FaceGroup: 2 start: 39148 end: 39159. Reading uniform faces...done.
FaceGroup: 7 start: 0 end: 3253. Reading mixed faces...done.
FaceGroup: 62 start: 3254 end: 39147. Reading mixed faces...done.
CellGroup: 63 start: 0 end: 6052 type: 1
--> FOAM Warning : Found unknown block of type: "316"
on line 74613
--> FOAM Warning : Found unknown block of type: "71"
on line 74632
Zone: 63 name: zone000:000--- type: fluid. Reading zone data...done.
Zone: 62 name: interior-p-62 type: interior. Reading zone data...done.
Zone: 7 name: geom-zone000:000---:7 type: wall. Reading zone data...done.
Zone: 2 name: geom-zone000:000--- type: geometry. Reading zone data...done.
Zone: 64 name: edge-p-64 type: boundary-edge. Reading zone data...done.
Zone: 214 name: edge-p-214 type: boundary-edge. Reading zone data...done.
Zone: 9 name: _zone000:000---:feature:9 type: boundary-edge. Reading zone data...done.
Zone: 3 name: _zone000:000---:feature type: boundary-edge. Reading zone data...done.
--> FOAM Warning : Found unknown block of type: "73"
on line 78743

FINISHED LEXING

Creating patch 0 for zone: 2 name: geom-zone000:000--- type: geometry
Adding polyPatch for unknown Fluent type geometry
Creating patch 1 for zone: 7 name: geom-zone000:000---:7 type: wall
Creating cellZone 0 name: zone000:000--- type: fluid
Creating faceZone 0 name: interior-p-62 type: interior
faceZone from Fluent indices: 3254 to: 39147 type: interior
patch 0 from Fluent indices: 39148 to: 39159 type: geometry
patch 1 from Fluent indices: 0 to: 3253 type: wall

--> FOAM FATAL ERROR:
Illegal cell label -1 in neighbour addressing for face 39148

From function void Foam::polyMesh::initMesh() in file meshes/polyMesh/polyMeshInitMesh.C at line 67.

FOAM exiting


@昝磊 好像是做凝固，李喜老师那的:huahua:

@Wayne 对，电磁场里面会用到

@bestucan 感谢前辈的热心回复，目前问题已经解决，问题出现在类phaseChangeTwoPhaseMixture的make文件的files文件里面，

LIB = $(FOAM_USER_LIBBIN)/libphaseChangeTwoPhaseMixtures //会报错 LIB = (FOAM_USER_LIBBIN)/libtwoPhaseChange //随便叫一个名字均可，均可编译通过，但不能是libphaseChangeTwoPhaseMixtures

自己也是偶然发现的，但也比较困惑为什么这么神奇:mihu:

不过最终问题还是解决了，很感谢前辈耐心的解答:xinxin: ，自己后面也应该好好补一点C++知识

@jdzhang 你的颗粒会动么？不会动，为什么颗粒间还有距离？会动，那对着“多相流”查比对着“多孔介质”查更容易找到解决方案些。

这有一个画网格的：
https://www.bilibili.com/video/BV1eS4y1F7Wm/?vd_source=c01e8d4fff0db460c75b32e0d0191b24

这有些相关的：
https://www.researchgate.net/publication/308996206_Pore-scale_numerical_simulation_of_mud_erosion_in_the_subsea_sand-mud_alternate_layer_using_lattice_Boltzmann_method

https://ascelibrary.org/doi/10.1061/(ASCE)HY.1943-7900.0001489

@东岳 是的，还在朋友圈看见老师说以后邀请您来杭州:chouchou:

https://portal.nersc.gov/cfs/visit/third_party/libccmio-2.6.1.tar.gz

连接挂掉了

IMG_20231127_193223.jpg
在计算单相流作为初场，可以收敛，但是湍流粘度比超限，继续设置为多相流空化计算的时候，开始发散，不知道有哪些解决办法，已经重新画过网格。
边界条件是速度入口，压力出口，湍流sst k-w，稳态计算

@李东岳 谢谢:threaten: :qinqin:

@coolhhh 好的，谢谢大佬~

已经解决。
（在一个视频中提到这一点）网格统计问题有三种情况：
第一种，是显示问题的网格数目很少，就几个，占总体的0.0002%。对于这种情况直接更改限定值。
第二种，是显示问题的网格数目很多，占总体的超过10%，甚至更多，对于这种情况需要更改MeshData的参数设置，重新对网格进行绘制。
必须保证，不能出现负网格。

@小考拉 alpha方程如下：61c83717-3135-4a66-9cdd-c1af5ba33159-image.png