@李东岳 好的李老师，里边的数据已经拷出来了。多谢。

@cccrrryyy 感谢感谢！给我指了一个好的修改方向！我做过RANS模型的对比，确实影响不大，但我看文献里有的是Smagorinsky算的很好，比如文献1，有的就是算的误差比较大，比如文献2，一头雾水:136: 还是得自己测试呀！

@sibo 我靠，前辈，我也是做结冰的

不仅是64以上不线性。10多核以上就不线性了。

@bestucan谢谢老师不吝赐教[抱拳]

@李东岳 好的，谢谢李老师，我再认真研究研究。

[General] Export full spreadsheet from Paraview -- CFD Online Discussion Forums (11_11_2022 8_57_13 PM).html

@bestucan 谢谢这位老师了，我的这个问题的确是李老师说的那种情况。

https://www.cfd-online.com/Forums/openfoam-solving/147135-how-switch-off-combustion-reaction-reactingfoam.html

https://github.com/vitst/dissolFoam

上来就自定义，完了还要跨版本号编译代码。把门远远甩在身后的人干起来也是横垄地里拉车，一步一个坎。:chigua2:

先照着视频教程之类的，实现自定义函数什么的。

提供的信息很难判断问题在哪

s.component(0) = 123 s.component(1) = 456

看看这样行不，这样行的话，那就是ux、vx的值有问题了

找不到 cyclicAMILduInterface.H 这个文件

搜着看

https://cpp.openfoam.org/v10/

https://www.openfoam.com/documentation/guides/v2112/doc/

@李东岳 好的，谢谢李老师。我发现面上梯度可通过fvc::reconstruct这个函数重构出的体心梯度，我还是得研究一下这个函数。比如刚刚的重力源项，定义在面上的话，我还没掌握怎么把它放回控制方程里面:143:

@李东岳
谢谢老师，大概理解this用法，这个是指向name()，还以为这个命名是起到什么作用，既然一般情况下没区别，那我就不深究了。

@小狗狗 在 有没有人用颗粒求解器计算过带AMI面的问题的？ 中说：

AMI和拉格朗日追踪

这个不清楚。之前处理的不是很好。openfoam-10把ami全部替代了。暂时没有测试。

@李东岳 在 如果用DPMFoam求解稀相流会怎么样？误差大么 中说：

@Alvin 不知道你的问题解决没有。我在我自己写的求解器中，通过

将p方程换为p_rgh方程 将g去掉

解决了。因为轻颗粒的浮力在方程\eqref{gra}中包含了，重颗粒也如此。就像 @dzw05 说的，连续相可以认为是单向不可压缩流，类似simpleFoam。 @Alvin 你试试将DPMFoam中的phig去掉，运算看看？

这个1,2两项是二选一即可吗，phig指的是不是phicForces里面的rAUcf*（g&mesh.Sf()）

@李东岳 交界面两侧网格任意形式的并行分区方式下的多进程任意coupling也可以？

请问大家有用过Cardiff的流固耦合代码solids4Foam吗？请问固体部分的阻尼要如何施加？solids4Foam中的固体力学部分是基于有限体积法来求解的。在传统的固体力学有限元法中，可以通过Rayleigh阻尼来给定结构阻尼，Rayleigh阻尼的阻尼矩阵是质量矩阵和刚度矩阵的线性叠加，solids4Foam中的一个固体模型unsNonLinGeoTotalLagSolid 中给定了一个质量比例阻尼，但没有刚度比例阻尼。
请问怎么在solids4Foam中对固体模型指定Rayleigh阻尼，或者对固体模型给定刚度比例阻尼？

@chen_hao

const volScalarField& omega_ = this->omega_;

这一种

@guohuiqun https://blog.csdn.net/dsfsdffgfd 你可以参考这个

@random_ran 您好，想咨询一下，这么多的网格您是如何进行后处理的呢？paraview 会崩溃的吧？

@四季之中的隐星
是的，根据你的需要进行修改。因为多数时候用不上，或者仅用作演示，所以都精简了。

@tidedrinker 在 functions里的forces在二相流里能否输出基于p_rgh的压力和转矩？ 中说：

@李东岳 这个forces是写在controlDict/functions里的，在求解过程中每个（或每若干个）时间步输出一次。这个压力和转矩是随时间变化的，不是根据结果文件用postProcessor -func那一套程序得到的。所以修改结果文件的文件名应该没什么用处？

请问一下forces在计算完成后，如果想要得到其他面上的forces，有什么办法后处理提取吗？看见您说不是根据结果文件用postProcessor -func那一套程序得到的,请教一下，这有什么比较好的处理办法吗？

请问在使用openfoam中工具forces工具，工具只能伴随计算中使用，如何通过后处理的办法使用forces工具（计算结束后，想要提取其他面上的力）？或者有什么后处理的工具能够任意提取自己所需要的面上的力？

@疏影横斜水清浅 网格太大导致的精度误差，不应该跟软件有关系，只跟你想要研究的问题、你关注的现象的尺度以及数值格式这些有关系。说白了，只要案例设置的好，精度误差不会跟别的软件有明显区别。

你提到的计算大尺度的开源软件和程序是指哪些？我这边做风场模拟，对这些还比较感兴趣。

各位大佬，我想在OpenFOAMv2206中植入我在OpenFOAM9中修改的模型。但是在编译过程中出现了错误：

9cb8d1c3-cfce-407c-a09c-cdecbf702dd3-图片.png

我的代码在OF9中是可以运行的，但是在OpenFOAMv2206中就报错了。看其他的代码感觉OpenFOAM9和OpenFOAMv2206没有太大的区别，这个错误是为什么呢？

Foam::autoPtr<Foam::viscoelasticLogLaw> Foam::viscoelasticLogLaw::New ( const word& name, const volVectorField& U, const surfaceScalarField& phi, const dictionary& dict ) { const word typeName(dict.get<word>("type")); Info<< "Selecting viscoelasticLog model " << typeName << endl; dictionaryConstructorTable::iterator cstrIter = dictionaryConstructorTablePtr_->find(typeName); if (cstrIter == dictionaryConstructorTablePtr_->end()) { FatalErrorIn ( "viscoelasticLogLaw::New(const word& name, const volVectorField&, " "const surfaceScalarField&)" ) << "Unknown viscoelasticLogLaw type " << typeName << endl << endl << "Valid viscoelasticLogLaw types are :" << endl << dictionaryConstructorTablePtr_->sortedToc() << exit(FatalError); } return autoPtr<viscoelasticLogLaw>(cstrIter()(name, U, phi, dict)); }

@李东岳 辛苦李老师了，方便的话可以试试看。后面那种场景，主要是公司内集群是所有专业混用，很容易出现这种情况。甚至单单只针对OpenFOAM来说，暂时没限制节点的使用，所以会导致一个节点上能同时跑好几个计算或者计算的一部分，分别是隶属于不同作业的。

@李东岳 老师，我自己觉得是可以，只是我没想好具体的实际方法，这也是为什么我选择在均匀的网格下弄。 我对于我现在均匀网格的就这么搞得 我自定义了我的固体活动范围里的网格编号和关系 均匀网格的连接关系很简单 我甚至连数组都不用开个去记录 我只是为了并行 预处理了固体活动范围内每个网格所在核心的编号和所处分块区域内我自定义编号的cell在对应核心对应分解域内的真实编号 计算时候根据连接关系，o(1)的运算时间取出临近网格的自定义标号 然后根据预处理记录的它在所属分块内的真实编号来计算 。

只是对于非结构网格，不均匀，有变形这种，我具体不太知道openfoam里编号的规律，也不知道老师您这里具体的应用场景是什么。我个人觉得，也是在进入时间步循环前，搞个预处理，把非结构网格的cell中心坐标提前循环一遍预处理，映射给自己自定义的一套均匀网格，然后一个数组记录，数组下标是自定义映射均匀网格的编号，该数组下标里数组对应内容是落在这个自定义网格内的非结构网格的真实编号，均匀网格的关系很简单，可以简单查找临近关系。而且这个映射网格的作用纯粹就是，映射个大致的位置关系。但是这个映射的均匀网格的尺寸要小于实际非结构最小cell的尺寸，大致我是这么想的。但是这样我不知道提前记录的内存开销有多大，我觉得甚至还不如，自己知道每个非结构网格要查找的周围感兴趣网格范围，每个模拟开始进入时间步前，提前预处理输出给文件打个表 （因为是预处理，所以无所谓，各种findcell也花不太久，因为只需要执行一次并记录），然后运算时候，对于某个非结构网格，直接读取预处理打表的那个文件，再取出周围符合关系的网格编号。这里读取文件，和字符串处理花一点运算时间。

可以。我记得湍流方程右边都是只考虑网格点的作用，不需要考虑边界。因此你的G声明了内部场，没有问题。

@zy_feng 哎哟，那我还是在10下面装的paraFoam，

我看下2206

@mecyhl

过 apt-get download 命令下载的包拷进去无法安装

为啥呢，我这可以

@郑学习 一直到毕业也没合作成功呢。不过，如果你有什么问题，可以在这里留言，我还是会尽力帮忙的。

@李东岳 谢谢！:140: :140:

@李东岳 李老师你好，那麻烦你推荐几篇非牛顿流体+颗粒流方面的论文啊，多谢了

@李东岳 您能简单推测大概是什么地方出现了问题吗？

具体研究体系是气液固催化反应，催化剂粉末与液体作为拟均相混合，气液不互溶，固体只存在于液相，流动过程中夹杂催化剂的液体会被切碎成小液滴。本人小白看了论坛上很多帖子，目前感觉这个MPPICInterFoam最符合我的体系，不知道有没有用过的前辈可以指点一二

@李东岳 好的，李老师，我去调试看看

最近在做风浪相互作用的模拟时遇到一个问题：在气液剪切过程中，湍流的研究是非常重要的，我在加入湍流模型之后，计算结果耗散特别严重，但是没有气体作用（nowind）时，SSTk-w湍流模型能够较准确的模拟，从下面的水位图中也可以看出来。并且湍动能k和湍流粘度nut在出口outlet逐渐向计算域内部扩散导致波高衰减（见下图），但是不知道具体原因。单纯的平底地形结果也没问题，但是遇到到这种斜坡地形就有衰减。
目前我已经尝试过of自带的几乎所有RANS类的湍流模型，计算结果均耗散严重（以SSTk-w为例），也尝试过对边界条件和壁面边界层的各种处理方式，还是效果不理想，想请教一下各位大佬，想解决这个问题，我还应该从哪些地方入手？
二维模型如图：
模型图.png
水位图：
水位图.png
不同时刻湍动能k分布：
k-12s.png
k-30s.png
k-60s.png
不同时刻湍流粘度nut分布：
nut-12s.png
nut-30s.png
nut-60s.png

@李东岳 是的李老师，我理解成了这个边界用在入口时相当于固定值，由inletValue指定，用在出口相当于零梯度，应该是我理解错了。。。

@veen 您好，关于“Primary job terminated normally...”的报错您后面还有研究过吗？我最近也遇到了相同的报错信息，不知道该怎么处理

你可能需要仿照MRFZone类写一个平动的类了，但是相对来说要更加容易些，主要工作是修正面通量，推荐用三角分解计算。还要注意下include的边界面上速度修正。

@我是河滩 请问您这个问题解决了吗？我最近也遇到了同样的问题

或者你应该在intU[1]=...之前加个if(cellID > -1).

想问一下，在绘制二维截面的流线图时，在paraview中以列表形式导出了点的数据，存在以下疑问：
①列表中各列数据代表什么含义，比如normals,integrationtime,point ID?
②各点为什么会有block number之分?block number是什么含义？
③能否通过列表数据做出流线图？
④保存数据时，为什么每个时刻的列表文件会被保存成2个csv文件？
以上问题求指教
3e5026fe-ea05-4101-86ee-b5ad8e43bfad-1666618809(1).png

刚接触OpenFOAM，安装的是OpenFOAM-v2206，了解到需要用pisofoam或者pimpleFoam求解器，貌似在tutorials里面没有现成的NACA翼型低速绕流大涡模拟不可压并行算例。

@cyberk 抱歉，这个问题我没有遇到过，但我看了你之前用ICEM画的网格，这个六叶Ruston桨的搅拌槽的网格不太复杂的呀，这个你可以试试只画1/6个区域，然后对称过去，画完整体之后再去openfoam中用topoSet划分搅拌的区域

@李东岳
谢谢老师，虽然只是一个小建议和思路，却给我节约了大量的时间，目前已经在OF7中调通了mutiphaseInterFoam求解器。:xiexie:

@hotairballon1997

在您回帖的链接中好像没有看到关于这个问题的解答，

对呀，要分类讨论，就像你说的，是用 PBM，还是拉格朗日。
我发的那个链接里 MPPIC 和 DPMFoam 中用到的粒子都是粒子包（particle parcel），这种用粒子包去代替粒子的方法我没了解过，但肯定已经不是物理粒子了。所以要看哪个层面，哪种模型的碰撞聚合了。