关于多孔介质仿真的压力问题(Fluent)

LGGGD-F

本人根据FLuent里的教程，按照以下速度和压降的关系，算出粘性阻力系数和惯性阻力系数分别为2.37e+7和17.04，然后按照26.2m/s的速度进行验算时，得到的压降只有几百，在运算时没有加入能量及化学反应，压力云图如图片示，这是什么原因？请懂的朋友指点下，诚心请教，谢谢。
速度           压降
m/s                         Pa
20.4                        10500
23.3                        14500
26.2                        18000
29.1                        19000

suwei225

不知道你说的结果不正确是什么情况，有几个我可以想到的希望能帮到你：
1.速度是催化器载体入口面平均速度，不是气体在催化器里面的速度
2.delta n是载体长度，但是注意单位是m
3.催化器载体内部阻力系数各向异性，你在阻力系数公式里求解出来的是流通方向上的阻力系数，其他两个方向无穷大
3.催化器内部温度和压力耦合很紧密，你在阻力系统公式里面肯定用的是定常空气参数，但是这些参数（包括气体粘度、密度等）在fluent计算域内不是常数，所以需要在初步计算结果输出后反向对粘性系数进行修正
4.fluent输出压力数值有很多选择，你到底对比的是静压还是总压，这个需要认真确定
5.一个小细节，压力输出时记得把衬垫排除
我能想到的就这些了

cholocate5

粘性阻力系数与速度一次方关系，还与介质粘度有关，惯性阻力系数是二次方关系，与介质密度有关。看看fluent关于这一节的帮助就清楚了。有什么问题，感觉你没表述清楚

LGGGD-F

谢谢你的回复，我就是按fluent教程说的做，得到粘性阻力系数和惯性阻力系数后，用29.1m/s做仿真，结果得到的压差很小。

suwei225

因为不知道你的阻力区域尺寸，以下说法仅仅为大致估计：
速度20多，压降在10千帕以上，个人感觉粘性阻力系数不会是e7的数量级
另提醒楼主注意，在计算阻力系数的时候，应变量是“单位压降”，而不是“压降”

LGGGD-F

谢谢你的回复，我确实把压降直接用来计算了，请问“单位压降”是什么意思？

suwei225

就是“压降”/“沿程长度”

LGGGD-F

这样算出来的结果是大了些，但还是不是正确值。我有两个疑问：一是，我计算所用的速度是催化器入口的速度，不是用气体在催化器里的速度，应该没错吧？二是，看到教程里有个delta n 来表示多孔介质的厚度，教程里用1m，对于催化器来说，应用多少？我觉得我算不准的原因在这里。

LGGGD-F

谢谢你这么详细的解答，除以delta n之后，算出来的值就比较接近了。有一个问题：根据软件的教程，按照速度，压差的实验数据，算出来的1/a是粘性阻力系数还是 （1/a）*密度  为粘性阻力系数？