Openfoam中的冷凝公式存在错误(珠状冷凝)
-
在文件
src\thermophysicalModels\thermophysicalPropertiesFvPatchFields\liquidProperties\humidityTemperatureCoupledMixed\humidityTemperatureCoupledMixedFvPatchScalarField.C
中,定义了一个可以用来计算潮湿空气冷凝的模型。该模型下,冷凝换热系数可以通过温度计算得到,源代码为:
Foam::scalar Foam::humidityTemperatureCoupledMixedFvPatchScalarField::htcCondensation ( const scalar Tsat, const scalar Re ) const { if (Tsat > 295 && Tsat < 373) { return 51104 + 2044*Tsat; } else { return 255510; } }
可以发现,这里程序应该是使用的温度为K去套用的冷凝模型,但是,通过查阅文献,可知:
BERGMAN T L, LAVINE A S. Fundamentals of Heat and Mass Transfer, 8th Edition [M]. 2017.
第632页文献中使用的是摄氏度,这里是否是源程序的错误呢?我查找资料,发现这个模型从Openfoam 的2.0版本好像就有了,然后一直延续的现在,一直没人发现?还是说我对这里方程的理解有误。跟各位老师确认下。
-
顺便,我想跟各位老师请教一个问题,就是他算例文件下
tutorials\heatTransfer\chtMultiRegionFoam\windshieldCondensation\0.orig\cabin\H2O
FoamFile { version 2.0; format ascii; class volScalarField; object H2O; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 0 0 0 0 0 0]; internalField uniform 0.01;//如果是体积或者质量分数,那为什么不是1?
这个H2O文件是定义的什么,mol量吗?还是说是质量分数或者体积分数。如果是分数,我在算例文件夹下无法找到定义的其他组分,还是说这些“其他组分”是被constant文件夹下的热工水力属性文件统一定义了呢,想顺便请教下点进来的各位老师(
让标题引了下流) -
摄氏度与K转换一下不就可以了么?源代码没问题啊
-
@李东岳 换热系数的返回值有问题,比如27℃,300.15 K的时候,换热系数的返回值在教材中应该是
51104+2044*27=106292
但是根据代码计算规则,得到的换热系数是:
51104+2044*300.15=664610.6
这换热系数的差距是很大的
-
@四季之中的隐星 这个是bug,你自己改一下重新编译试试。
-
@李东岳 是的,我的意思就是这里存在这么个错误,一直没人发现还挺奇怪的
-
@四季之中的隐星 openfoam里面有些模型非常非常小众,用的人特别少,再加上即使有人用了,也没看到源代码,同时即使有人看了源代码,也未必看到错误。只有非常仔细的人,对了原来的公式,才能找出来bug(比如你这种)。
另外,你确认下你看的SCI里面是不是对的,有时候SCI可能也有笔误
-
@李东岳 好的,明白,这个不是SCI,是国外的热力学教材,我横向对比下其他文章里的公式确认下