@禾斗石开 fluent

自己写的代码终于可以出结果了，但是用cgns存坐标数据的时候发现只能存一半，很奇怪，代码举例如下
三维的面网格，节点数2x11，因此是1x10个网格, 坐标储存在x,y,z数组中

real, dimesion(2,11) :: x,y,z CALL CG_BASE_WRITE_F(INDEX_FILE, BASENAME, 2, 3, INDEX_BASE, IER) ISIZE(1,1) = 2 ISIZE(2,1) = 11 ISIZE(1,2) = 1 ISIZE(2,2) = 10 ISIZE(:,3) = 0 CALL CG_ZONE_WRITE_F(INDEX_FILE, INDEX_BASE, ZONENAME, ISIZE, Structured, INDEX_ZONE, IER) CALL CG_COORD_WRITE_F(INDEX_FILE,INDEX_BASE,INDEX_ZONE, RealSingle, & 'CoordinateX', x, INDEX_COORD, IER) YZ坐标同上

读取时代码为

real :: node(2,11) node = 0.0 call CG_COORD_READ_F(INDEX_FILE,INDEX_BASE,INDEX_ZONE, 'CoordinateX', RealSingle, (/1,1/),(/2,11/),node,IER ) write(*,*) node

node的读取结果为

x(1,1) x(2,1) ... x(1,6) 0.0 0.0 ... 0.0

请问有没有大佬能指点一下

@XIN 请问您解决了吗，我也遇到这个问题了

@李东岳 好嘞 ，谢谢李老师

应该和左下角的 tolerance 相关。
这tolerance都那么小了，还能12点覆盖

估计是把一个几何体修复成一个点了，所以崩了

@yingqing 不知道为啥，有一次成功了，在弄就不能成功了。节点的位置进行了设置，但网格的节点与我设置的节点不重合。就像这样
177924ad-d438-468f-8cf5-a612fe049465-image.png

周老师，这里我测试发现好像这种结构网格的加密方式对于三维模型可以，对于二维网格，在经过局部加密后再导入到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中

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

感谢各位大佬不吝赐教。

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

以上是我的两个文件

网格文件为ENSIGHT GOLD FORMAT

@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

有没有人用OpenFOAM里面的滑移网格AMI功能的，想了解一下AMI的并行做到什么程度了？对网格面的分区划分有限制吗？可以支持任意网格并行分区吗

@hustzzw 微信图片_20221101093137.png
就像这个图中里面圆柱区域和外面圆环柱区域

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

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

@李东岳 感谢李老师。我去学习下:140:

@bestucan 嗯嗯，感谢大佬。检查过，check和Quallity其实都没有问题，再次感谢~

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

@haining-luo
我现在用ansa处理极复杂三维模型，保存.mesh格式面网格，用fluent meshing或star ccm生成多面体网格或切割体网格，挺方便。
如果模型太复杂，搞科研可以用icem，工程上不太建议用ICEM生成非结构网格。

@coolhhh 好的，谢谢大佬~

@队长别开枪 您好，我在OF8中编译了ccm26ToFoam，想模拟波物相互作用，转化后的网格checkMesh也是ok的，转化之后的网格我放到OF v1912里面进行计算的(还没编译ccmToFoam)，但是相同的设置用OpenFOAM自带的sHM划分的网格可以计算，ccm26ToFoam转化的网格却在0.02s就突然发散了。想问一下有什么解决办法吗？:135:

有的时候icem就是会报错，只能没做一些操作之后多保存，然后再进行。。

ICEM画二维网格导出时报错：
WARNING: Mesh has uncovered edges. ANSYS Fluent needs a complete boundary (lines in 2D) or it will give a variety of errors and not read in the mesh! If this was 2D Hexa, perhaps your edges are not associated with perimeter curves

模型描述：二维管道流体的测试网格，使用ICEM导入模型后：
1、划分块
2、在preMesh 里 使用Edge Paras进行边界层的绘制，
3、preMesh
4、在preMesh里Convert to Unstruct Mesh
导出网格选二维，其他未改变，出现上述报错
新手从B站学了一点划分网格操作就尝试了这一个，不知道哪里出错了，还希望前辈们解答

@韬智tz ESI版本的OFv2012之后有一个ccmToFoam工具但默认没编译,在这个网站搜一下能找到一个编译的教程

@下里巴人 看到这个帖子想起我之前做那个流化床模拟，加了跟管子就死活不收敛了，估计和你这个一样的问题。
话说有尺度差异的时候出来的这个网格好呆……明明圆形的楞给画成多边形了。
我想起Gmsh他们组最近有个博士答辩的东西就是做curve mesh的，不知道能不能用在计算上了。。

@shrine 没关系，无关紧要这些文件

@李东岳 把圆弧去掉，全部改成六面体网格，还是发散。
从细管流入水池就发散了。
附了log文件
Screenshot from 2022-03-01 09-14-34.png

Screenshot from 2022-03-01 09-12-07.png

Screenshot from 2022-03-01 09-12-19.png

log

@yingqing 遇到了，我是根据错误提示去排错，然后搞定了。我遇到的大部分的报错似乎是缺少必要的文件。

@李东岳 ansys18版本 linux
链接：https://pan.baidu.com/s/1QpXN_GhW5tCv1rNJrlBUYA
提取码：1234

你的图就是节点对齐。

有文献里是中心六面体网格，边界层为四面体问题网格。

你可以试试hex-core，在做体网格的时候，可以成圣hex-core

@李东岳 好的~谢谢东岳老师:xinxin:

QQ浏览器截图20220113163626.png
问题：2021r2版本，在封闭形成流体区域的过程中，无法选择面。有哪位前辈指点一下吗

对简单的长方体进行网格划分，再检查网格质量时发现在长方体的周围有很多有问题的face产生，也不知原因，很疑惑，求大神指点，万分感谢！

0T5OC4QVS4OWJO998JSUAU7.png

@bestucan 在别人电脑上是正常的 ，很奇怪，我最后重装电脑解决了

缺这个包：binutils

或者路径没配置对，可以用echo $PATH看看有没有/usr/bin这一项

网格数量太多了。15510 * 16000 * 22个网格？

@bestucan
感谢您的回复！“不同方向的流动计算差异就大”这个解释很直接明了，我突然就有点明白这个网格线出现的原因了:xinxin:

snappyHexMesh边界层这个问题有的时候确实不太好处理。目前我也没有什么太好的建议。看看有没有其他大佬关注关注 :jingya:

@samuel-tu
前辈好，请问交界面网格畸形的问题您是否解决了？我也遇到了类似的问题:135:

@aeron 一会儿试试，感谢分享~:146:

这个主要看cfd之道公众号流沙的教程就可以学会了（哈哈哈哈，你是我拿到注册码之后回复的第一个人）

@tens 好的，谢谢你，我试一下

@shrine

我上面说过的，字典文件中有一个坐标，这个坐标落在几何里面，生成的网格就在里面，落在几何外面，生成的网格就在外面
66148774-6b5c-42db-a783-f365dc9e147d-image.png
2.在生成几何时就对每一个面进行分part并命名，然后就可以对part分别命名，我是用ICEM画的stl几何，里面可以进行这个操作，其他软件我不了解

@cccrrryyy 有所耳闻，有空玩玩 :duang:

@李东岳 尖角部位，加一个竖直的翼墙，这样平分后角度好点，但是附近网格不再是水平的，这样对interFoam计算有影响吗？

@李东岳 谢谢东岳老师。一个朋友给的建议，先生成开尔文波区的几何stl，用sHM对该区域进行体加密~

@李东岳明白啦，谢谢东岳老师。

@如果随风 pointswise可以画

@ycx1042341217 @huningdong 请问两位，壁面第一层网格高度可以做到什么程度？我发现最小只能设置到大约5e-5m:136:

:136: 好吧找到原因了，是这个版本的link edge factor默认是0

好的好的谢谢大家的意见~！！！！非常感谢大家，根据大家的意见已经调整好了，最近再准备毕业事宜所以没有及时感谢大家！

@李东岳 同一个计算域的网格，有单元的节点信息，单元相连的单元信息，面信息。但是两个不同计算域网格之间没有这种拓扑信息。这种情况下如何并行剖分网格同时保证负载平衡，减少通信量