可能gcc版本不对

大佬问题解决了吗，看你贴的报错信息貌似是代码里面这段直接”=“出错了？
/home/dong/foam/foam-extend-5.0/src/foam/lnInclude/PtrList.C:54:18: error: invalid conversion from ‘Foam::Field<double>’ to ‘Foam::CoeffField<double>’ [-fpermissive]
54 | ptrs_[i] = (a[i]).clone().ptr();

各位老师好，我在使用interdymfoam计算船舶破舱入水时把六个自由度都打开了，结果发现三个方向的速度位移都一样，这到底是什么原因呢？fae8d152-b728-40b4-aec0-47f779d0d9ee-Screenshot from 2024-12-09 08-46-56.png

是的老师，打算评估一下各种非牛顿曳力模型作为小论文一部分内容

本人在使用openfoam自带的重叠网格进行计算。
按理来说在设置初场的时候两套网格重叠的区域是不参与计算的，但是本人在计算过程中，重叠网格区域好像也参与了计算，本人想就这个问题问一下各位大佬。
这个是被重叠的网格：
6384ab7c-e651-4dc9-9bd2-16298be04a95-432b47d32d1e3106b90bbfeb518ec39.png
这个是背景网格：
e2221afd-ff60-4bdc-b424-758b0087efb0-8d2e870096a5f44c86bbcf28b4318c3.png
![b6f02d62-3983-441a-9196-ca13443e58db-fb193eac98c7188e6dfc2a67c5ec6cf.png]
这个是重叠后的网格
4e960bd4-99a3-45d3-8c01-8f2c8888dabf-fb193eac98c7188e6dfc2a67c5ec6cf.png
这个是本人设置的初始条件：
b57ddbe1-1225-4cc4-847d-243257530ba2-913b4b3992b0d17eedda2894abcc303.png
由于初始条件只能给两个网格，并且重叠网格重叠区域不参与计算，所以就给了一个圆柱形的高温高压的初始条件。

计算过程中总是因为局部速度过大导致计算发散，这是计算结果。
这是没有进行Threshold的计算结果
777b8258-58a1-41ae-b84d-c6775eddfaca-4a3190545fd101c6ad980a5cbbf4455.png
这个是分离过后的结果
5fd95244-6427-43a8-a9d9-c8c1a0bb9b6f-bf547f2742f3a21d4d4fd8e035395ec.png
可以看出未进行后处理分离的结果出现了特别大的速度导致了计算的不稳定性，但是按理来说，重叠区域的网格是不参与计算的，上面的结果显示他好像参与计算，有没有大佬知道是为什么。
这个是这个算例的设置条件controlDict
dynamicMeshDict fvSchemes fvSolution
计算的求解器是自己仿写的，但是用他们自带的overcompressibleinterDyMFoam也是这样的结果，有没有大佬知道是为什么。

@katnisswen OpenFOAM有自带的主动消波技术，你自己看一下就知道了，

@范准 已解决

首先中间的这个面似乎不应该这样加，因为boundary目录里定义的应该都是边界，但多出来的这个面是在计算域内的，并不是边界，我觉得可以使用Sample或者facezone
其次这个报错应该是因为定义的这个面实际上跨越了两个面，实际上应该像这个“top”或者“bottom”一样把两个面分开来声明，应该就不会有问题

但是interFoam似乎只能读取p_rgh而不能读取p文件

interFoam求解的是p_rgh，所以读取p_rgh，但是p也读取，p都是calculated。interFoam的p都是calculated

@李东岳 我因为可能要做cfdem所以用的5.x跑的

5.x改了好多次网格尺寸和dt都在浮点溢出，后来看网上教程给piso改成pimple模式 倒是不溢出了，但是这个结果就很离谱了

感谢老师，这下在主程序能找到创建的场变量了:xinxin:

@学流体的小明 在 展向平均 中说：

@东方白杨

东方白杨 在 展向平均 中说：
这一段意思就是对流场同时X、Z方向平均，得到沿Y方向的分布吗？

是的

东方白杨 在 展向平均 中说：
如果我只想得到流场沿X方向平均，得到一个YZ面，就是流向平均的一个截面，这该怎么操作呀？

这两个工具都是同时对两个方向进行平均的，如果想得到一个截面，我觉得你可以研究一下代码然后自己改改程序，或者再找找，可能别人写过相关的工具。

好的，感谢！我研究一下

@李东岳 计算域网格是？老师可以给一条链接吗，我学习一下！

http://dyfluid.com/code.html

里面有个// 输出形变率，你可以参考一下

具体是在of2212里面的DPMFoam，如果不支持的的话能不能采用paraview去计算，主要是得到U.air的du_dx,du_dy......等这些量，由于已经计算完成就不考虑在求解器里面植入这个功能了。谢谢各位大佬给出建议！！如果在描述中有问题恳请大佬指正！

@李东岳 好的，李老师

@李东岳 老师这个题目就是我问的:135:

应该是不满足颗粒堆积的条件，太稀疏了

你不需要改求解器。你在算例层面，注入两种粒子就可以了

这是我历史以来看到过最简单的招聘需求。