CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    Openfoam中的冷凝公式存在错误(珠状冷凝)

    OpenFOAM
    openfoam-冷凝蒸发
    2
    8
    673
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • 四
      四季之中的隐星 最后由 编辑

      在文件

      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去套用的冷凝模型,但是,通过查阅文献,可知:
      5888ae1e-e4d1-4676-b2f5-adb20cc1fd5e-image.png

      BERGMAN T L, LAVINE A S. Fundamentals of Heat and Mass Transfer, 8th Edition [M]. 2017.
      第632页

      文献中使用的是摄氏度,这里是否是源程序的错误呢?我查找资料,发现这个模型从Openfoam 的2.0版本好像就有了,然后一直延续的现在,一直没人发现?还是说我对这里方程的理解有误。跟各位老师确认下。

      1 条回复 最后回复 回复 引用
      • 四
        四季之中的隐星 最后由 编辑

        顺便,我想跟各位老师请教一个问题,就是他算例文件下

        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文件夹下的热工水力属性文件统一定义了呢,想顺便请教下点进来的各位老师(让标题引了下流)

        1 条回复 最后回复 回复 引用
        • 李东岳
          李东岳 管理员 最后由 编辑

          摄氏度与K转换一下不就可以了么?源代码没问题啊

          CFD高性能服务器 http://dyfluid.com/servers.html

          四 1 条回复 最后回复 回复 引用
          • 四
            四季之中的隐星 @李东岳 最后由 编辑

            @李东岳 换热系数的返回值有问题,比如27℃,300.15 K的时候,换热系数的返回值在教材中应该是

            51104+2044*27=106292
            

            但是根据代码计算规则,得到的换热系数是:

            51104+2044*300.15=664610.6
            

            这换热系数的差距是很大的

            李东岳 1 条回复 最后回复 回复 引用
            • 李东岳
              李东岳 管理员 @四季之中的隐星 最后由 编辑

              @四季之中的隐星 这个是bug,你自己改一下重新编译试试。

              CFD高性能服务器 http://dyfluid.com/servers.html

              四 1 条回复 最后回复 回复 引用
              • 四
                四季之中的隐星 @李东岳 最后由 编辑

                @李东岳 是的,我的意思就是这里存在这么个错误,一直没人发现还挺奇怪的:mihu:

                李东岳 1 条回复 最后回复 回复 引用
                • 李东岳
                  李东岳 管理员 @四季之中的隐星 最后由 编辑

                  @四季之中的隐星 openfoam里面有些模型非常非常小众,用的人特别少,再加上即使有人用了,也没看到源代码,同时即使有人看了源代码,也未必看到错误。只有非常仔细的人,对了原来的公式,才能找出来bug(比如你这种)。

                  另外,你确认下你看的SCI里面是不是对的,有时候SCI可能也有笔误:134:

                  CFD高性能服务器 http://dyfluid.com/servers.html

                  四 1 条回复 最后回复 回复 引用
                  • 四
                    四季之中的隐星 @李东岳 最后由 编辑

                    @李东岳 好的,明白,这个不是SCI,是国外的热力学教材,我横向对比下其他文章里的公式确认下

                    1 条回复 最后回复 回复 引用
                    • First post
                      Last post