如题，有没有人在MPPICFoam中添加过多孔介质源项的？怎样实现多孔介质阻挡颗粒啊？

@oitocfd 我记得是网格有些问题，不过后来我用waves2Foam算了，自带的还有olaFlow我都用过，感觉还是waves2Foam好用一些

@是嘿嘿牛啊 感觉是网格太粗 你的数据采点更密 一个间距里采到两个一样的值 不太清楚

@dingcy 自问自答了，在tecplot里输出非压缩的视频文件，然后放到ps里，增加分辨率，然后再渲染视频，最后的视频就清楚啦。

@李东岳 谢谢李老师

@李东岳 非常感谢东岳老师的回复，我会更改检查一下的。

@bestucan

如果是计算时要实时计算边界条件，这个可能符合你的要求：
http://openfoamwiki.net/index.php/Contrib/groovyBC#Usage_of_the_boundary_condition

感谢回复，需要的确实是在计算中实时计算的边界条件，我先了解一下:xinxin:

又仔细看了一下，两个版本应该是差不多的。
相对速度是因为CdRe不是曳力系数 而是曳力系数×雷诺数 所以可以用粘度、直径还有密度换算成相对速度。

@李东岳 老师，内部也在动的，截图可能看不出来，只是这个边界更明显。
75376c62-5086-4608-879f-b7823f7ebaa8-1639391378.png
7a6ca8ef-4d95-4e82-b6e1-77c1203adfab-1639391406(1).png

@bestucan 是的，原因应该是相似的，感谢您的解答。

@bestucan n 谢谢教授的解答

@shouan 这是我的pointMotionUx文件
/--------------------------------- C++ -----------------------------------
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: v2006 |
| \ / A nd | Website: www.openfoam.com |
| \/ M anipulation | |
*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class pointScalarField;
object pointMotionUx;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 -1 0 0 0 0];

internalField uniform 0;

boundaryField
{
inlet
{
type uniformFixedValue;
uniformValue table
(
( 0.000000 0.001250 )
( 0.020000 0.001300 )
( 0.040000 0.001500 )
( 0.060000 0.001550 )
( 0.080000 0.001700 )
( 0.100000 0.001900 )
( 0.120000 0.002000 )
( 0.140000 0.002200 )
。。。
。。。
);
}

bottom { type zeroGradient; } reef { type zeroGradient; } outlet { type fixedValue; value uniform 0; } atmosphere { type zeroGradient; } Back { type empty; } front { type empty; }

}

@zhe 好久没来了才看到回复，感谢！！

@袁宝强 我记得可压缩求解器里面有个可压缩系数，psi，不知道是否有用

@队长别开枪 非常感谢您的帮助！

@李东岳 谢谢李哥

不清楚 过几天准备写个弹簧试试看能不能准

@wwq 在 请教关于k和epsilon壁面函数的问题 中说：

https://www.openfoam.com/documentation/guides/latest/doc/guide-bcs-wall-turbulence-epsilonWallFunction.html

非常感谢:146:

@李东岳 谢谢李老师！A fifth-order shock capturing scheme.pdf 李老师，我是想把这个文献里的第二部分模块数值算法植入进去，我昨天上传的是第二部分前面一块，后面那一块主要是进行重构，去计算q_L和q_R，然后通过前面的公式去进行计算，前面那块我看和您发的rhocentralfoam解析中很像，但后面的公式不知道应该怎么植入，没有思路，想请教一下老师:chouchou:

@李东岳 因为换个地方继续算，重新分块的时候就根据最后时刻的网格来划分，所以用scotch 就是每个块基本同样网格数（上下几十浮动），而用simple划就会是那种XYZ方向的拉升