6ce0cdba-2578-4c0f-ab51-4ceef0b045cb-image.png
配置了一台双CPU的服务器，且每个CPU是52核，Windows系统运行Flunet发现一个任务只能调用一个CPU且不能超过32核，如何解决使其能够调用两个CPU并超过32核的限制？

@yuno 致敬，帮顶！

勾选Data Sampling for Time Statistics

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

@tinghaik 个人觉得Gridstar是十分值得学习的，虽有不少缺点也已十分优秀，而且免费不用考虑盗版问题，使用起来与pointwise比较接近。而flowstar主要面对高速气动方向，这个就要看个人的专业了复合性了。我认为这些国产软件是国家队的主流产品，已经在逐步推广，不会轻易虎头蛇尾，有基础的学起来难度也不大，如果专业对的上，何乐而不为。

较粗的网格：
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

Fulent里做一个旋转域，类似于一个大圆柱套一个小圆柱，旋转域和背景网格两个计算域都涉及到流固耦合传热，且由于计算域具有周期性，只取30度的扇形区域。请问这样可以实现吗？

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

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

@kennyadapter 感觉你可以把Oblock范围扩大，同时减小一下oblock内的网格层数就好了

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

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

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

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

@Yongbo :high: 谢谢

@李东岳 哦哦，没问题了。:142:

@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 找原作者拿到了 ` 谢谢回复

原来这样。

最近在算好多LES算例，直接把服务器当主机用了。日常用没问题。但是发现最大的问题就是word文件的编辑。因为我这面很多都是别人发来的word文件我在上面继续编辑。但是linux下面用各种word编辑软件都会出现格式问题。有些文档就两页，但是导出pdf的时候格式就不对。有些是商务文件，格式对不齐，感觉不是很专业。

我先试试wps，看看会不会好一些。

@五好青年 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），这证明植入的模型还是不容易收敛是吧？

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

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

@Ding 你需要一句一句的看。
具体方法是：从上至下把无关的先注释掉，一句一句往里加，每一次都用fluent进行编译、初始化、模拟，看是加入那一句之后才出现的问题，这样才能看出是哪里出现了问题。
你这样单独看是很难发现问题的

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

pisoFoam不可以考虑重力

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

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

@Zhy2022 :high:

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

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

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

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

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

@李东岳 老师的解答很详细。我有一个小疑问：不调用高阶格式，网格也不够细，但是模拟结果还是能得到涡旋。这样的模拟结果能称为quasi-DNS么？

@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 ) ),

请问有用Fluent做锅炉模拟的同行嘛？可以一起交流一下~

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

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

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

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

ipmitool -I lan -U ADMIN -H 192.168.1.18 sensor thresh FANB lower 100 150 200

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

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

参考文献
Demirdzic, I.; Peric, M., Space conservation law in finite volume calculations of fluid flow.