fluent 里的porous 边界如何在openfoam中实现

大家好，之前我发了一个帖子，大意是基于interFoam，在连续性方程中添加源项实现质量的输入和输出，进而造波：
http://www.cfd-china.com/topic/943/添加源项后相方程求解越界-请问这种情况有人遇到过吗

经过多次调整，质量源造波的程序已经基本上完成了。我按照上贴子的公式，把源项加入相方程中并开启了MULES修正，同时在压力方程中也添加了源项，这样求解的时候，质量源区域内的alpha值（默认为1，源区位于水下）大概在1.0*（质量输入）-0.9*（质量输出）之间波动（因为我的源项是类似于正弦波动的）。

然而，我采用相同的方法修改压力方程，使用动网格求解器interDyMFoam的时候，上述的方法出问题了，形成了类似于空化的现象（不懂空化，希望没说错），即相体积分数变得很小，最直观的描述是在水下质量源区域会形成一个气泡，后面相体积分数的求解救完全错误了。

我仔细看了一下求解器的文件，发现interFoam.C和interDyMFoam.C求解时主要的区别在于动网格求解程序中多了一些通量计算以及压力方程中的部分不同。但是由于之前没看过动网格的求解器，因此还是看的不太懂，请问有做过这方面的大牛吗？可否分享一下经验，为什么网格不变形的时候完全正常，但是用了动网格求解器就不行了。:big_mouth:

感激不尽:sunglasses:

是啊，相信CFD在国内也会有不错的发展！

因OpenFOAM现有格式不满足需求，需要实现新的插值格式，有没有熟悉的人说明一下思路，此外我不清楚MeshObjects文件夹下的 upwindCFCCellToFaceStencilObject和upwindCPCCellToFaceStencilObject的作用和区别是什么？谢谢

surfaceScalarField cPhi("cPhi", c & mesh.sf()); fvScalarMatrix c ( fvm::ddt(c) + fvc::div(cPhi) == source. );

@wwzhao 谢谢

@小龙 filter 搜一下就可以找到

继Mcgraw于1997年提出QMOM之后，Marchisio et al. 于2005年提出DQMOM，Yuan and Fox于2011年提出CQMOM，Yuan et al. 于2012年提出EQMOM，Marchisio and Fox于2013年提出ECQMOM，Fox et al. 于17天前又提出HyQMOM和CHyQMOM并发布预印本。

https://hal.archives-ouvertes.fr/hal-01632813

@李东岳
级测试，如您所说改为NO_READ后不再出现警告，谢谢

@random_ran 非常感谢！问题已解决

@李东岳 那我该怎么提取速度梯度中的量呢？

gradU.component(tensor::ZX)

这样对吗？

我记得这个好像自带的mixer2DAMI也是这样的？

是的你说的很正确。

@mohui
考虑过，但是刚开始压力发生这种令人不解的变化，之后压力又能够恢复正常，气泡脉动在流场产生的压力又在合理的范围内。但是我对程序进行定点调试，就发现经过我图3的代码之后，压力就突然变了。

5.0版本的可以正常安装，但是低版本如2.4.0安装就会报错：

-bash: /public/openfoam/OpenFOAM-2.4.0/etc/bashrc: line 140: syntax error near u nexpected token `;&' -bash: /public/openfoam/OpenFOAM-2.4.0/etc/bashrc: line 140: ` [ "$FOAM_V ERBOSE" -a "$PS1" ] && echo "Sourcing: $1" 1>&2'

四面体分解确定插值权重啊，用于相间耦合或者场平均（跟用多少耦合有关）。跟重心关联，还可以表示成局部坐标的形式。

能不能帮看一下这个问题，谢了！http://www.cfd-china.com/topic/1150/openfoam中cyclic周期性边界的问题

请教各位老师，有没有研究涉及OpenFOAM与PFC3D的耦合的。老板让用PFC3D研究DEM-CFD coupling，感觉很少有研究这个的，目前看的资料也就是github上一个itasca公司的老师编写的一系列代码和一篇今年6月份发表的文献，资料有点少，也不知道这是否是一个正确的方向，还望各位老师指点指点，谢谢各位老师！

@mohui 是的，interface处是有间断的，但是one-field形式的控制方程体现不出来间断
压力分布不对，看一下理论解就明白了～

不是少部分，是一点点负值，比如说-0.00006。如果你的方程有一些要对浓度做开方，那一点点负值都不行。