CFD中文网

@搬运工

高斯迭代法处理20阶矩阵就已经非常慢了。通常和多重网格+高斯赛德尔求解器一起处理。并且OpenFOAM这种成熟的商软里面矩阵求解器算法非常精炼。通常由主求解器+预处理+光顺+重组构成。多重网格更是复杂。另外，在编程的时候使用动态、寻址的方法寄存系数。

所以如果你要写一个快速的求解器：

使用外挂包； 潜心研究稀疏线性系统求解理论；我建议花大量时间研读：YousefSaad, 2003。看看大型稀疏矩阵怎么求解的。而不是3/4阶矩阵的高斯塞德尔求解器。:cheeky:

@一介匹夫 小的线段没有识别成一个边界。

@txwwbz
还不确定

@murphy 我主要考虑的是过程时间非常短，传热可能并没有那么迅速，所以不施加两相传热模型反而更加符合实际情况

@wsy11 谢谢您~我现在已经可以成功构建网格啦^^

@Junren-Hou 在 ICEM划分网格 中说：

可以试试check block，有些时候其实没有问题，就它软件自己的毛病，然后block里最后有一个check block，可以让他自己检查一下。里面有好多选项，都试试，什么check/fixed block

回忆了一下，我的旋转区域用ICEM画的应该没有问题。出问题的应该是我画的外流场区域，一方面追求地面边界层，致使网格尺度过小，一方面为了平衡网格数量，过渡不均匀。猜测应该是外流区域的地面边界层设置导致了相关报错的产生。

@东岳 i其实是第一个网格，n是最后一个网格
也就是
先求每一个网格中的alpha * x 再对整个计算域求和
然后求和所有计算域中alpha的值
两者做比值

欧拉多相流的话不是应该是网格尺寸大于颗粒尺寸吗，和气泡也会有具体的关系吗？

topoSet有一些比如寻找点最近的面，寻找距离stl最近的面之类的，你试试

icem也有这个功能，只是我没用过，

https://cfd-china.com/topic/6069/搅拌反应器icem结构网格划分方式/7

Screenshot from 2024-09-28 14-40-10.png

最近还在写自己的solver，但是感觉并没有用到封装这类C++的高级应用，还是当年fortran风格的写公式。。。可能还是初期吧~

请教一下你是每一层加密都对三个方向进行加密吗，shm中能能不能只对z方向进行加密呢。@一朵苍凉

@pippowong
可以简单说说这种非常大尺寸的大气模拟和普通的流场模拟的区别么？以及相关的侧重点？

还没见过教材里面讨论low-Mach-number assumption，气动声学这块我还没研究过 :zoule:

Source[kth].ref() += (1 - kth)/(kth + 1)* pow(a[n], kth + 2)*g*w[n];

上面要改为

Source[kth].ref() += (1.0 - kth)/(kth + 1.0)* pow(a[n], kth + 2)*g*w[n];

否则 (1 - kth)/(kth + 1)在kth大于1的时候为0

东岳大哥您好，算几十万的网格是没有问题的，我扩展了虚拟机的内存到64g和虚拟内存到10g都不行。 。。也是奇怪

ｏｐｅｎｆｏａｍ　边界条件引入到离散控制方程求解中，实际是把边界节点的值引入到了源项来看待的吧？最想不明白的就是再靠近壁面的体积单元Ｐ，再其通过惯性力/压力差/剪切力形成动量方程的时候，是不是还要把壁面剪切应力考虑进来。

AFAIK, you are the only one working in this area. Any update?

多面体明显质量好过关 过基础计算关 相对正六面体 同样参数控制 数量变小了