请教各位前辈一个问题：
268a789a-7dc8-4835-a14b-9a555865d688-image.png
phi是一个surfaceScalarField，和volScalarField不矛盾吗？

此外，现在需要对速度散度判断，以确定n的取值：
ab4b0384-55c3-4a21-b5aa-677b1b44c21f-image.png
n是一个volScalarField或者scalar也可以，要计算：
4da41b5d-9c03-426c-be01-c69940e0402a-image.png
写循环的话，应该用div（U），还是div（phi）？（ps和pc都是volScalarField）

@samuel-tu 主要这文章看看也很有帮助，它的角度偏代码整体架构而不是力学，对理解openfoam的设计还挺好的。

较粗的网格：
Screenshot from 2021-05-10 19-30-28.png
0e914623-9079-48c7-866e-45d35efac3d0-image.png
较细的网格：
Screenshot from 2021-05-10 19-19-03.png
05cabc1c-9a22-4303-bae7-6a4a24d7b145-image.png

之所以压力泊松方程离散出来之后得到的矩阵是对称的且不是严格对角占优，是因为通常泊松方程里面的 laplacian 项是用中心差分来离散的。楼主看一下 src/finiteVolume/finiteVolume/laplacianSchemes/gaussLaplacianScheme里面 fvmLaplacian 函数的定义就明白了。
动量方程为什么对角占优呢？主要来自于

对流项，特定的对流项格式可以保证离散之后对角占优 时间项，时间项的离散得到的系数只会进入系数矩阵的主对角。 隐式松弛，从构建系数矩阵的角度来看，隐式松弛和时间项达到的效果一样，都会增加对角占优 特殊的源项，比如 fvm::Sp 得到的隐式源项，也是只贡献于系数矩阵的对角元素，可以增加对角占优（不正确使用则会减弱对角占优）。

@小狗狗 谢谢您，我学习一下

请问大佬们，这个是什么东西？应该怎么加载？感谢感谢:140: e0dbfb1e-bf91-480c-8eba-929d4cf423d6-image.png

你好，我之前使用过interFoam求解器，在看这个多区域求解案例的时候看不懂了，这里的d58ea828-33c7-438e-8dc6-dd820a1625cd-image.png 是什么意思？是指的哪个边界条件？请大佬帮忙解释一下，感谢

@李东岳 确实~ 现在改了改代码 能跑完算例 但是结果还是有点问题 等我完全搞定了 把代码传上来

@Yongbo :high: 谢谢

@bestucan 后来我也是这么解决的。。相当于直接安了新的Paraview

@merlinizlf 我就是加密一点后 然后尽可能让网格均匀 网格别太细长 就好了 当然也有时候是其他原因导致的 但是毕竟改网格是最直观容易的第一个尝试选择

可以参考邱小平大佬的博客
链接文本

读完东岳老师关于icoFoam的文章后，有两个问题不太清楚：

icoFoam可以计算laminar和quasi-DNS，那么在设置计算参数的时候这两种情况是如何区分的呢？比如是离散格式的设置不一样吗？

使用icoFoam 和 使用pisoFoam + laminar 是等同的吗？

希望各位老师告知~

我想把Clausius-Clapeyron 方程引入，里面有一项 ln(p/pSat_)，编写为log(max((p/pSat_),1E-8))之后，编译成功，但是跑算例的时候就会显示下面的错误。同理，在对(p/pSat_)进行pow函数，甚至直接多次相乘时，都会报类似的错误。如果定义一个常数的volscalarField变量one，进行log（one）计算就是可以正常运行的。
请问是这些函数对volscalarField变量有特殊规定吗。

#0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in /lib/x86_64-linux-gnu/libpthread.so.0 #3 ? in /lib/x86_64-linux-gnu/libm.so.6 #4 Foam::log(Foam::Field<double>&, Foam::UList<double> const&) at ??:? #5 void Foam::log<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?

@starDust boundaryFieldRef()是非const的可以赋值

@bestucan 找原作者拿到了 ` 谢谢回复

@五好青年 emmm……这就不知道了，我没研究过。这只是一个线索。现成的应该是没有的。但是照着pdf里的线索查比从头做要快多了。

@Micro 关注过，但是没研究过。好像是8刚支持外挂矩阵求解器。还没机会学习:shangxue:

@shiyu 一个有点冗余的办法，定义一个vector场，每个网格都等于你那个gradP。这样就可以在别处lookup到了。

meanDiameter { type coded; libs ("libutilityFunctionObjects.so"); name error; codeExecute #{ const volScalarField& d = mesh().lookupObject<volScalarField>("d.alpha.oil"); scalar d32 = d.weightedAverage(mesh().V()).value(); if (Pstream::master()) { std::ofstream file; file.open ("d32", std::ofstream::out | std::ofstream::app); file << mesh().time().timeName() << " " << d32 << "\n"; file.close(); } #}; }

我最近也凑巧要用这个，我用的上面这种方式。

@李东岳 李老师，请问发散后继续计算又不发散了，比如如下报错：
a0c00201-ea8c-4d8d-bdad-e4e033c1ad9e-image.png
或者报错与solver有关（gauss-seidel），这证明植入的模型还是不容易收敛是吧？

单个叶片受力我之前听别人说过用犀牛建模可以对模型上的面命名，这样的话可以求单个叶片的受力

@李东岳 感谢东岳老师的回复，这模拟的是双流体，主要是关注离散相沉积

机翼某一截面（比如10%位置处），这个位置可以通过平面与模型相交，做出交线。
但是交线是一条曲线。
请问通过openfoam或者paraivew，怎么导出这条曲线上的压力值呢？
感谢谢

pisoFoam不可以考虑重力

interFoam可以考虑中立，你这么操作就是对的，没有别的操作了。

检查下网格边界条件吧。感觉有点问题。。

@Zhy2022 :high:

最后就想请教下，各位大佬

在进行OpenFOAM并行计算时，不管是虚拟机（还是真Ubuntu系统），CPU的核心可以100%满载吗？

比如电脑16核，设置16核。有人尝试过么？实际会乱跳么？

我这并行计算，实际只占用5%，这也太坑了叭！:136:

@Samuel-Tu: 好， 但仍谢谢Tu兄的回复，Tu兄，辛苦了

@zhangxc0223
Uc都是液相速度~
这一部分是生成一个cacheField，在后面calcCoupled中会根据parcel的坐标和cell中心做插值，得到这个DUcDt这个值，然后在计算parcel的受力

@浪迹天大 谢谢，我研究一下，这个file 我就是读了一下文件名字，然后放到IOobject里，原来在of4就完全没问题

positionsFile_(this->coeffDict().lookup("positionsFile")), positions_ ( IOobject ( positionsFile_, owner.db().time().constant(), owner.mesh(), IOobject::MUST_READ, IOobject::NO_WRITE ) ), PDFFileName_(this->coeffDict().lookup("PDFFile")), PDFsamples_ ( IOobject ( PDFFileName_, owner.db().time().constant(), owner.db(), IOobject::MUST_READ, IOobject::NO_WRITE ) ),

@cresendo 谢谢提醒，我就是做了一点点修改，有的时候会产生挺大的问题:mihu:

各位简直就是大佬中的大佬 厉害厉害 :chitang:

最近在看divDevSigma这个函数，请问老师标题指的“不可能完成的任务”是什么意思:jingya:

东岳老师，我最近在做一个类似的模型，颗粒粘度与其体积分数有关的，用的of8的multiphaseEulerFoam，现在完全没有思路，不知道要将其放在哪里:136: 望您指点一二

查到了一个类似的问题，是fluent用tecplot处理的时候也发现壁面速度不为0：
出现这样问题的原因为：fluent的计算数据是存储在每个单元的中心位置，而tecplot的云图显示是按照节点来显示的。如果是直接导入fluent的.cas文件和.dat文件的话，那么在tecplot中，节点上的值会通过单元中心点的值插值得到，从而造成壁面速度不为0的情况。
不知道OF这面有无特殊处理方式

是我没有解释清楚么？这里的genotype[i]就是输入的字符，因此minus 1.0 2.0 对应三个字符串，从i=0开始，对应minus, i++变为1， 返回第一个递归函数evalGeno对应的i=1,也就是进入else, i++ 为2，numValue为1.0，继续减号，然后到第二个递归evalGeno，此时i为2，所以按程序流程就是1-2，但是程序测试结果是2-1的结果，为什么会这样呢？请各位老师指点一下

这里稍微自挖自填下吧，我感觉OpenFOAM中ldu的矩阵存储格式只能支持 对称矩阵的存储和寻址（对称指的是非零元的位置对称，大小可以不相等）。当然，还没有经过验证。对于（非零元）不对称的情况，我实在想不出来怎么用LDU去存储，难道说，OF就不只是一些高阶迎风格式吗？可是我看OF中好像也支持二阶格式啊，咋回事啊。。。。要是有大佬清楚麻烦稍微告知下！谢谢！

@bestucan 感谢回复，这么一说确实是，就没法达到稳定状态。

不过如果按照您说的，用拉格朗日粒子，如果我设置粒子为passive scalar，用icoUncoupledKinematicParcelFoam对一个不稳定的场来求解粒子，这样粒子速度也会不断增加，最终导致发散吗？

https://dcwan.sjtu.edu.cn/userfiles/22.pdf

问题已解决~因为变量不能从字典中正确读取，修改读取语句后可以正常运行

有大佬在虚拟机上跑CFDEM吗？不太懂虚拟机，请问可以下载东岳老师分享的安装好OF-5x的虚拟机然后再安装LIGGGHTS耦合吗？

@bestucan 好的哈！谢谢版主！我仔细看看您发的网址

@bestucan 好的哈，谢谢教授回复！不过您说“但是经过相邻面的速度也反了”，速度的方向是（绝对）唯一的吧，就是说该是多少就是多少。感觉还是从高斯定理出发，法向量永远取外侧为正

@OItoCFD 在 新写的程序单核计算无问题，但无法并行计算？ 中说：

mpi并行那些需要考虑吗？

一般不需要考虑

除非你的surface场进行一些特殊的插值

第31页（页码13）最上面那段是不是你要的

https://www.jianguoyun.com/p/DYNwnmQQsLv2BRiS3fED

c02ffe4b-6329-4d30-ab7a-39442943d1ee-image.png 难道勾选这个不是显示圆柱实体吗？

@bestucan 感谢版主回复！我把farField设置为无反射，结果还不错；我也觉得几何模型画大一点，自由射流更好。