OpenFOAM的原生多孔介质模型是否可以在变形网格中的静态网格区域实现功能？

[liujm]

大家好！

如题，OpenFOAM的原生多孔介质模型是否可以在变形网格中的静态网格区域实现功能？我的目的是想让多孔介质模型像六自由度刚体一样运动。如果我有描述不清楚的地方，大佬们可以瞄一眼这个链接，我想做的是让这里的多孔介质区域动起来。

所以，想问问大佬们这个功能是否可以在不改代码的情况下实现？

现存问题：我在设置sixDoFRigidBodyMotionCoeffs的时候，需要指定patches对象（如下图），可是patch需要在blockMeshDict中boundary中创建，一旦创建了boundary，那么就是固壁边界了，而不是多孔介质模型了。希望大佬们能指出几条解决方案~

如果需要改代码，也请大佬们简单指出需要更改的库和注意点~~

[liujm]

手动艾特佬中佬 @李东岳

烦请李老师解答~

[liujm]

以我对代码的理解，多孔介质模型的区域是基于cellZone划定的，源项的计算也是在cellZone中进行的，与网格号是绑定的，因此想要实现我的目标，多孔介质模型代码应该是不需要做更改的。
问题应该出现在一楼图中的

diffusivity quadrtic inverseDistance (cylinder)

以及sixDoFRigidBodyMotionCoeffs中的

patches (cylinder)

部分。这里需要给变形网格区域指定一块交界面，所以这个交接面该如何给定是要解决的问题。

[李东岳]

动网格的时候，需要计算面的受力，因此需要patch边界。如果你用多孔介质，不需要patch。但是如何计算边界的受力？

[liujm]

@李东岳 李老师，那不可以按照多孔介质模型本身的计算方法计算吗？就像这张图，D赋值很大，效果等同于固壁边界。我的提问不够专业，还请老师耐心解答~

[liujm]

@李东岳 李老师，如果没有办法避免创建patch，那我用多孔介质模型把它包起来，这个思路可行吗？我这边尝试模拟一个二维圆柱在水中自由落体的算例，但貌似算了一下就报错了，从报错信息上看不出哪里有问题，我把报错log和算例关键设置也一并放在下边~
图1

图2

log文件
log.mpirun

[李东岳]

你可以试试创建计算域网格，然后多孔介质内部全部都是多孔介质，这样也存在了patch

[liujm]

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