CFD中文网

@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

@bestucan 和论文原作者确认了，说我这样是等价的，调整了一下其他地方，目前算出来对上了。。

该书11.8节（中文版305页，英文版411页）中的例2，推导源项的表达式时，高阶格式减去一阶迎风格式，一阶格式的值取上上游的值φu，为什么这里不取φc呢？先谢谢了。
7a735b6e-01c9-4088-971a-359144da8189-image.png

管内气体流动，管壁定温加热，我在计算的时候设置了report：热流密度，对流换热系数，Nu等参数。发现用热流密度反推的对流换热系数和Nu和报告的对不上，相差很大，请教各位老师这是什么原因？谢谢

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

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

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

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

各位大佬，最近想用OF计算固定的多体装置，由于在大域内生成的结构过于粗糙 需在小域里生成 然后用mergeMesh合并到大域里去，但是小域的边界没办法删除
请问各位大佬，有没有什么设置建议？删除边界或指定特殊的边界？

您好，请问您这个问题解决了吗？能分享一下解决办法吗？

thinkpad t480s i7 8550u 16G内存
平时写论文、做ppt用，也可以开虚拟机改求解器算算小规模问题，真要数值计算还得远程服务器:140:

6b6b5aa0-11d9-4944-accc-20bb6eb70b36-image.png 表中的custom rate exponent是怎么确定的？感谢各位。73ea2e62-c36f-4287-9b1d-b62745bea525-image.png

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

9818ab58-ba59-4644-af8a-2b91924553cc-image.png
模拟环形截面流道内往复振荡流动，边界条件设置为正弦速度入口,余弦压力出口。理论及文献中曲线在180度时分界，而我的模拟结果发现曲线随周期呈整体右移趋势，在15周期分界点移动到192度。请教各位老师这是什么原因？

@袁宝强 啊我知道了 谢谢 :xiexie:

钝体绕流来说，升阻力的主频和流场监测点的速度主频应该是相同的。

感谢各位大佬！:chitang:

自己写的代码终于可以出结果了，但是用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

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

@veen 从源码编译of是个挺痛苦的事情，建议不把of安装在root，可以考虑试一试easybuild，它是基于module的环境管理，应该比较成熟了。我从0开始，最开始在PC上面编of-3.0，到后来在实验室的机群上home编译of-4.0, foam-extend等等，多版本的管理用module还是方便些。印象中，gcc（如果你用gcc的话）和mpi的版本很重要。

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