自带的例子中，NS方程的源项 包括重力项和表面张力项都是放在压力泊松方程中迭代的
248e597c-c754-48af-980a-37713b5bb86d-7c072d73eab735a11c026b606429f8d.png
我在想如果把体积力源项直接放到速度方程那部分：
d07dffba-d5bb-406a-a8f0-115e2e67f7f4-7cff84172630dee75cf72c6111da8db.png

也就是说把源项的影响放入到UEqn.A和rUA中03723c61-3207-44b2-924b-dc418f5c021f-image.png 是不是也没问题？
这样就不需要把源项加到泊松方程中了会简单一点。
试了一下计算也能跑。

请问这种方法是否有问题？如果没问题的话，和放在泊松方程中处理相比，算法上有什么优劣的说法吗

@bestucan 多谢大佬这么详细的回答啊，非常非常感谢！我在用多相流模型的时候温度没有上升也没有产物，但是我在使用组分运输模型的有些情况温度有上升，我也不敢确定反应是否发生。不知道是不是这种比较新的化学反应不能使用fluent自带的化学机理呀，是不是要自己导入他的化学机理呢？

@liuzh 您好 想问一下您最后是怎么解决的？

@ivan 同学 问一下您后来解决counter-balance这一项的设置了么

@随风飘姚 您好 我正遇到这个问题，我设置了outflow出口，但是进出口质量就是达不到平衡。请问可以和您交流一下么？

@dongsc 您好 想问一下您的问题最后是怎么解决的？不知道能不能加您的QQ？

@cccrrryyy 我是Fluent软件学习使用者，我的理解是Fluent中的Velocity-inlet类似于OF中的inlet，outflow类似于OF中的outlet，实际上velocity-inlet设置为负值时也可以理解为出流。初学者的一点理解，还望大家批评指正！

注册omega如下，有地方不理解：

omega_ ( IOobject ( IOobject::groupName("omega", alphaRhoPhi.group()), this->runTime_.timeName(), this->mesh_, IOobject::MUST_READ, IOobject::AUTO_WRITE ), this->mesh_ )

IOobject::groupName("omega", alphaRhoPhi.group()) 知道是定义omage的name，但是 后面这个alphaRhoPhi.group()，是起到什么作用呢？

如果我定义成：

omega_ ( IOobject ( "omega", this->runTime_.timeName(), this->mesh_, IOobject::MUST_READ, IOobject::AUTO_WRITE ), this->mesh_ )

会有什么效果无法实现？:tishizi:

@李东岳
针对 ncompressible::turbulenceModel 模板特化后的 ShihQuadraticKE 模型 ，仔细研究了一下其怎么实现非线性项的实现，发现其是在 NonlinearEddyViscosity.H 中定义 nonlinearStress_ 项，然后用 ShihQuadraticKE 作为子类继承父类 NonlinearEddyViscosity 以便在 ShihQuadraticKE.C 中直接引用omega_。

但由于我修改的是SSTkomegabase.H( 继承于TurbulenceModel类 )与 linearViscousStrss.H ( 继承于BasicturbulenceModel原则类 ) 没有特例化，不能简单类似上面的ShihQuadraticKE 模型。需要从新仿造框架，是一项大工程！

我的策略是：在现有的文件上直接硬植入，简单很多，毕竟涉及到湍流模型的类错综复杂，针对找不到omega，按照我的分析应该是对的，运行求解器时linearViscousStrss.C先运行此时并未运行SSTkomegabase.C 所以没有omega场注册，找不到他，所以我在creatfield.H中IOobject一下omega场，目前可以算了，结果还未检查，不知道我的理解有没有问题。

计算一个心脏血流模型，在导出数据时报错：“ERROR
Surface '000001B4E04832A0' should be a region2d category.
”321ff4ea6a63fea6aa4ff4ed99c4e3e.png
请问有什么解决办法？

Data-Create Zone-Rectangular提取出想截取的位置
d2469623-64c0-45c3-952d-a6339a3b560b-image.png
但是提取出来的Zone没有数据，Data-Interpolate根据原数据对创建的Zone进行插值就可以
abd5cf18-caaf-4a49-a528-3b9afd696f35-image.png
这样就得到一个名为Rectangular zone的Zone
758ce629-1262-40bf-83b2-1477d15a474d-image.png

tecplot应该可以吧？
先提取出来流线，然后对流线积分

cfd151b6-cb6a-4521-a36a-a3fb27fe3b41-image.png

a4a297dc-0192-46eb-8194-342bf20663cd-image.png

@李东岳 在 集群上并行测试OpenFOAM，并行效率并没有比单节点提升 中说：

-genv FI_PROVIDER tcp你这个去掉不能跑的话。你如何确定走的是infiniband，而不是以太网模式。我们这面跑openfoam不需要这个参数。我们之前测试也出现过你这种情况。后来我们换交换机硬件了。但目前我还不确定现在我们这5节点能到什么样，得下周能出个测试结果。另外，openfoam离散设置差异（比如GAMG那个），我个人感觉不会引起特别大的差异。不过你可以实测看看，我也不100%确定

@xpqiu 这位大佬之前好像测试过2048个核心 :mianmo:

老师，想问一下，用openmpi在自己组里面的集群上用pbs跨节点并行（10g以太网交换机）命令用的mpirun --mca btl_tcp_if_include <ip地址> -np reactingTwoPhaseEulerFoam -parallel，运行的时候发现计算的节点上cpu的用户进程占比us只有50-60%，系统进程占比sy有40-50%，这个问题有没有什么好的解决办法？

Info << patchCentre.size() << endl; // 若注释掉上一条Info语句，在这里输出patchCentre的size得到的值变为了0

换成这个呢？

Info << patch.faceCentres();

https://www.cfd-china.com/topic/4790/smagorinsky模型不输出k_sgs场/5

最后的方法你试一下

没理解，按照你这个并行计算的话，后面reconstructPar为啥会有问题呢

可以在 OpenFOAM cavity 算例里执行 blockMesh 后，查看 constant/plymMesh 文件夹。就是这些东西。
不过记录形式向会有区别。但都是典型的“非结构网格”。

网格系统，自己开发的通用性不强，处理起来全都得自己手写一套。当然好处就是不用去学成熟的网格系统。便于上手，但越往后，想要的功能越多，越麻烦。

如果是 OpenFOAM 的用的网格系统，paraview 就可以看，也可以用 blender 结合 blockMesh 生成。snapHexMesh 也可以。对网格、流场的各种处理也可以调用 VTK 库直接在程序里进行。

看了一下你上传的数据，没有像样的文件头，估计是刚接手一套祖传代码:136:

@dxl 噢，那可以在相应的位置加一些测试语句试试。

@李东岳 感谢大佬，我再去了解一下这个URANS

@bestucan 好的，谢谢您。

@firejet C_STRAIN_RATE_MAG，你想问的是这个吗？

会不会是质量有问题？sum变成了0.5就不对了吧

69cdd51f-33c7-4640-9840-dea844d49ba6-image.png
请问，想把这个设置成以红点为转轴，左右摆动的摆片，该怎么设置？

@刀尔東 在 向各位大佬请教一下网格划分的基础问题 中说：

感谢各位大佬不吝赐教。

链接: https://pan.baidu.com/s/11z9O_mzDLuqEp30mWRr5Ow?pwd=1234 提取码: 1234 复制这段内容后打开百度网盘手机App，操作更方便哦
--来自百度网盘超级会员v6的分享

以上是我的两个文件

网格文件为ENSIGHT GOLD FORMAT

@李东岳

@李东岳 在 能不能设置两组颗粒（一组one-way，一组four-way）同时模拟？ 中说：

我目前有个想法是添加两个不同的拉格朗日cloud，看看能否分开处理，但我现在在出差手头没有openfoam还不确定这么做是否straightforward。
欧拉场这面可以算一下速度的物质导数就是加速度$D \bfU/D t$

你好李老师，我最近也面临这个问题，我现在是尝试用solidparticle做one-way，kinematicCloud用作four-way，但是在求解的时候报错说找不到U文件，两个粒子的cloud分开用的时候都没问题，一起调用的时候出错了

更新一下
发现了问题的源头，是网格的问题。
我使用ICEM画网格，或者OpenFOAM的blockMesh工具画网格，都是先画背景网格，然后对近壁面进行多次加密。加密是八叉树类型的加密，所以在网格尺寸变化的地方，大网格的尺寸直接是小网格的二倍。这样的网格确实不好，可能一些拓扑结构上的不好导致了singleGraph程序进行插值时找不到插值的基点。
解决办法：
举例说距离壁面0~1的网格尺寸是1,1~2的网格尺寸是2，那么就使用singleGraph在0~1这个范围内提取一条线，在1~2再提取一条线，这时的type可以是uniform，也可以是midPoint或者midPointAndFace。一句话，就是分段提取信息。
我是重新画了网格， 老师说之前画的这网格根本不能用:mihu:

@希望先生 硫酸老师说他上传的版本坏掉了，缺东西，找时间再上传一个新的:135:

@ct 我处理时也出现了负值，好奇怪啊

@adcpk 超算每个节点20核，是因为只有20核还是只用了20核？如果是只用了20核，其他核在跑别的作业吗？如果只有20核，建议换成15+15这样来测试。

一般情况下超算平台的CPU性能肯定没有你工作站的好，工作站肯定比超算新吧，快一些的话其实也可以理解。

@stratification 你好请问这个问题解决了吗

这种边界条件不是特别稳定。很容易发散。所以总是报错。

@李东岳 谢谢李老师，我先去找点相关的文章看一看吧

@李东岳 感谢李老师！现在对这个求解器理解了一些，对于$\mathbf{U}_1$的控制方程, $K_d \mathbf{U}_2$一项作为显式项被直接纳入pEqn.H的界面通量的修正。

非常感谢！已更新

嘿，我还没关注这个历史。还真是。太折磨人了。

@韬智tz 在 OpenFOAM模拟空化气泡一些参数的问题 中说：

那最大气泡半径怎么得知呢

这个应该是随着模拟自动发展出来的吧？

你可以直接采用湍流模型来算试一下。用湍流模型计算层流，只要湍流粘度足够小，影响就很小。

@李东岳

sudo apt-get install libcgal-dev

这样就可以编译了

openFoam可以无量纲的，应该在etc/controlDict这个文件里里吧dimensionset 1改为0

@nanxuan 请问你解决这个问题了吗？能否分享一下？谢谢

@星星星星晴 谢谢星星前辈，这就去试试:146:

各位前辈好，我采用open foam并行计算，想用tecplot后处理。导入controlDict后只有模型显示，没有计算时间信息。另外在open foam中重组的时候出现以下报错Create time

--> FOAM Warning :
From function int main(int, char**)
in file reconstructPar.C at line 220
No times selected
不知道什么原因。请问有没有什么办法可以让我用tecplot后处理的。

@winsway_zero 千万级别笔记本开paraview比较难，还是得找台本地的服务器。

数据传回本地你可以写个脚本啥的自动去跑，不过这个得看集群的网络。我们是大量千万网格级别的计算随时回传的，感觉速度也还行。不过我们的计算是稳态不可压，算是最简单的了。

有效的后处理看你需求了，如果同类型的案例，只要一开始能够看到三维的展示就好，这样能明确知道哪些地方需要关注。后续写成各种function配到controldict里面，计算完成之后就直接跑一遍，这样生成的结果文件数据量小一些。

是的，MPPIC把拉格朗日粒子看成parcel了，泡泡的话，还是粒子比较好

@rachel0096 我没看过他这个 但基本都大同小异 其实 写一个ibm在openfoam就用不了多少代码 我也就看了几篇文章 找准一篇你觉得好实现的对着做 你可以先用最简单的方法 比如在pimplefoam里 对于2d圆柱绕流 写一个函数检测小于半径r范围内的cell 这样就检测出一个圆柱范围内的流体cell 标记为1 其他流体的cell标记为0 然后在压力方程那里 根据固体cell和desired velocity的速度差算出体积力 加给phi那里 进入压力迭代 这就是一个最简单的ibm 这个弄对了再加上二阶插值分布函数那些来考虑更高精度的界面识别

@李东岳 明白了，谢谢李老师耐心细致的解答。

背景是模拟固体火箭的射流问题。
考虑到固体火箭产生的射流中存在$Al_2O_3$颗粒，需要采用多相流的方式去模拟流动。目前我采用Fluent中的欧拉-拉格朗日方法初步开展了流动仿真，对比有无颗粒，流场还是有变化的。但是接下来采用DO方法引入颗粒辐射后，计算得到壁面处的辐射热流值很小，似乎都不影响总热流结果。
有人有这方面的经验不？

@李东岳 感谢东岳老师提供的思路，我在代码中把ShihQuadraticKE改了个名，目前在源文件中用Allwmake编译成功了。
不过孩子还是很不明白怎么单独拿出来编译这类非线性湍流模型，看来暂时只能通过这种方式代替了。

@李东岳 在 DES模型 中说：

这个找找他们做湍流的原始文献应该会有详细的介绍，就是这些DES模型为什么做这些更改。

好的 我再深入学习下，谢谢李老师

@bestucan 好的，谢谢老师