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，是国外的热力学教材，我横向对比下其他文章里的公式确认下

CFD中文网

Openfoam中的冷凝公式存在错误（珠状冷凝）