CFD中文网

压差是分析很多工程设备的关键参数，在之前，主要都是针对多孔介质的Ergun模式来处理这一块。
也跟试验对比很好。但是，最近新开发的一个模型，虽然与试验的结果对比比较好，但感觉总是有问题。
于是不知道CFD中静压值到底是怎么得出来的？
对于一块多孔介质好像比较好处理，但是如果这个时候多孔介质表面沉积了灰，这个时候这块的整个压降会升高，然后灰越堆越多整个表面的压降都会升高。然而，这个局部压降和表面灰的分布十分相关，造成整个表面的这个灰分的压降相差一千多Pa。而最后一般的检测压降点，都是压差变松器的两个接口，一般都是设置在容器外壁的。但是如果我模拟检测点和试验一样的时候，这个时候的压差比我的表面计算压差平均值要小得多。（压差是1200多Pa，但是我的多孔介质有800多Pa，加上平均粉尘层的压降有640多Pa，合起来应该有1400多Pa，而不是1200多Pa.）
这就产生了一个问题，这个通过N-S方程求解的压力项到底是怎么处理的?
当他通过一层多孔介质的时候，多孔介质内部有压力源项，而在外部没有，那么当他气体从多孔介质外向内流通的时候，他的这层压力和速度是怎么计算的？如果这个时候我再添加一个源项在多孔介质层，而这个源项在多孔介质外是没有的，那么这层额外的源项是不是也会通过增加流体的静压来表达呢？（通过伯努利方程似乎很好理解，但是N-S是怎么实现的，他的这个\Delta P 是怎么获得的？）

@ejector 我试了一下，这样是没问题的。谢谢了。不过周期性时却少了，为啥啊？

现在解决了吗 我现在也是用toposet建立了cellzone 但是也是报错没有对应的cellzone

这个是二维的，plot-vector-variables，设置好两个方向的物理量比如x方向速度，y方向速度，就可以画出来通过任一点的流线了。因为流线本身实际上是y关于x的函数，在速度坐标系中等于y方向的速度关于x方向速度的函数。

@gemini 谢谢。刚刚试了您的方法，解决了！

以前简单看过，好像是研究表面活性剂的方法。大概是在三维曲面上用类似FVM的方法离散控制方程并求解。

各位专家好：
在多孔介质模拟的过程中，我发现了一个问题，虽然在发表文章中省略了。但是确实是一个直观的问题，我的模型，RANS方法中，不管采用哪一种，在多孔介质的近壁面处的拟合的速度与实测值总是具备一定的差距。这种差距我分析认为是在较高阻力降的多孔介质环境下，气流在流经多孔介质表面时，一部分流体会因为高阻降而产生流向变化，形成比较复杂的近壁面流场。而在这个近壁面流场中，因为Fluent默认是采用的多孔介质域或者是多孔跳跃面，应该是对这层近壁面没有处理的，依然按照充分发展层处理，造成这层近壁面的流速（甚至是流向）与实测值都有较大区别。
虽然能够较好的拟合实时压降（Endo based equations），但是如果在近壁面流场不符合实际的话，那么很显然对于最后的表面覆尘的预算难以符合实际情况。不知道大家对于这种多孔介质的近壁面流是否有过研究，或者能够有比较好的解决思路？欢迎大家提出宝贵的想法和建议，感谢，再次感谢！
下图为提出新模型对于滤管外壁的模拟和实测结果对比：
0_1544421734755_compare.png

@shrine 没关系，无关紧要这些文件

@bestucan 你是行家呀，厉害厉害，终于知道原因了

@下里巴人 看到这个帖子想起我之前做那个流化床模拟，加了跟管子就死活不收敛了，估计和你这个一样的问题。
话说有尺度差异的时候出来的这个网格好呆……明明圆形的楞给画成多边形了。
我想起Gmsh他们组最近有个博士答辩的东西就是做curve mesh的，不知道能不能用在计算上了。。

@李东岳 好的！终于知道我为啥模拟结果不对了！谢谢李老师！

@bestucan 先说声谢谢，我去仔细看看

用pisoFoam计算了两百雷诺数的层流圆柱绕流问题。但是计算出来的升阻力系数出现了拍现象，有大神可以指点一下是什么问题吗。以下是升阻力系数图和算例设置求解器.png 离散.png 层流.png U.png P.png Cl.png Cd.png

@东岳 刚性翼型用的pisoFoam，柔性翼型用的是fsiFoam。但是pisoFoam算出来的系数波动很大，而fsiFoam算出来系数波动就很小。不知道是为什么

都柏林大学发来贺电！

晚上又测试了一下，在linux系统下是可以顺利的进行计算，谢谢东岳老师的帮助:xinxin:

DEM加速比很高
CFD好像一般

moment method对于nucleation和condensation的模拟简单吗？

即为简单。但是溶解/蒸发数值上要复杂些。2017年我在HZDR的时候同事用fixed pivot method模拟增长，我用moment method，相比来看moment method要更简单，并且不需要处理相空间的drift：https://www.sciencedirect.com/science/article/pii/S1004954118301174/pdfft?md5=84216a8842ce3fd6d1a99d4f327c1ee1&pid=1-s2.0-S1004954118301174-main.pdf

@李东岳 谢谢:threaten: :qinqin:

已经解决。
（在一个视频中提到这一点）网格统计问题有三种情况：
第一种，是显示问题的网格数目很少，就几个，占总体的0.0002%。对于这种情况直接更改限定值。
第二种，是显示问题的网格数目很多，占总体的超过10%，甚至更多，对于这种情况需要更改MeshData的参数设置，重新对网格进行绘制。
必须保证，不能出现负网格。