请问大家有用过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工具（计算结束后，想要提取其他面上的力）？或者有什么后处理的工具能够任意提取自己所需要的面上的力？

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

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

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

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

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

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

@李东岳 谢谢！: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的边界面上速度修正。