CFD中文网

感谢分享

@cfd-china 没有必要手动进行分块，我只要保证有一个processor的INLET边界是完整的就好了，其他的怎么分块比较随意，实际手动分块的话也可以，我为了简单，就用了simple了。

记得之前好像有关于手动分块的方法的帖子，我找找看:happy:

interPlicFoam求解器已经并入扩展包geometricVofExt，仓库地址：https://github.com/daidezhi/geometricVofExt

@散漫守望2016 好贴，支持

哎呀我勒个去 当时自己可真能写 现在看了都想删了 不嫌害臊

PRESSURE GRADIENT TERM这个图中的$F=\frac{m}{V} \nabla p$，其中的$V$应该为$\rho$

@李东岳

常见的问题是求$\mathbf A \cdot \mathbf x=\mathbf b$中的$\mathbf x$，可以有直接解法（如LU），也可以用迭代法（如Gauss-Seidal），求逆矩阵肯定是下下之选，原因有几点： 操作太多，$~O(N^3)$; 稀疏矩阵的逆矩阵不是稀疏的； 你要求得$\mathbf x$还要再乘一次，$~O(N^2)$ 逆矩阵求了之后没有用，因为CFD方程组是非线性的，每次求的矩阵$\mathbf A$都不一样，所以你求得的逆矩阵只用一次，下次迭代又得重新求解。如果是针对线性问题，每次针对不同的$\mathbf b$求对应的$\mathbf x$，用逆矩阵还勉强说得过去。 现代的CFD程序似乎迭代法用得更多，JFNK, IDR(s)，最不济还有CG和BiCGSTAB.有的连矩阵都不需要组装，何苦去求逆...

@队长别开枪  好的，谢谢你的回复，我再看看

对了，这个问题我在2015年曾经汇报过bug，Henry直接表示建议采用高版本的OpenFOAM:sunglasses:

@arctan 非常感谢！

回复: 两个相切的圆如何进行关联
@kscandal
同学，我说的话，你确定，你懂？

9个月过去了，仔细看了一下一楼的结果，依然没看出什么区别..

@zxyf1 是的，我的模型是个圆柱体，只需要让高度方向尺寸越来越小就可以，我已经实现了网格移动，就是把它用于interFoam时算的alpha场不对，目前 还在找原因，有兴趣可以一块探讨一下。

本人石油相关专业，就读于石油大学（北京），硕士
现在油价萎靡不振，企业裁员信息不断，一线员工也是连连降薪。
虽然就读于石油专业牛逼的学府，但是深知自己行业知识不牢固，根本不懂什么是行业需要。
真是内忧外患，愁！

最近在家呆的闲，又没有实验室的资源，，想看看磁流体的流固耦合的代码和书最好中文的，，，求推荐啊～

@西湖冷月
nice兄弟，d=====(￣▽￣*)b

感谢鸭蛋分享

dyfluid@dyfluid:~/twoWayHyQMOMFoam/tutorials/PTC_2D/testCase2$ oneWayHyQMOMFoam oneWayHyQMOMFoam: symbol lookup error: oneWayHyQMOMFoam: undefined symbol: _ZN4Foam6HyQMOMC1ERKNS_7PtrListINS1_INS_14GeometricFieldIdNS_12fvPatchFieldENS_7volMeshEEEEEEE

原因：当时调整了src下的的构造函数，紧跟着solver那面的构造函数也加了一个参数进去。然后重新编译了src，忘记了编译solver，出错。
方案：重新编译solver后解决。

#0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigSegv::sigHandler(int) at ??:? #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"

原因：在class里面增加了一个私有成员，编译库之后，求解器可以对类进行初始化，随即报错。报错只有3行，没有任何其他信息
方案：编译class之后，尤其这种增加私有成员的时候，记得重新编译solver一下就可以解决

这实在是太详细了
虽然对我来说要能全部看懂还有些困难:zoule: