CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    流固耦合alphat使用alphatJayatilleWallFunction报错跪求大佬指点

    OpenFOAM
    流固耦合 湍流热扩散系数 alphat边界条件 湍流自然对流
    2
    4
    926
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • 一
      一颗鸭蛋 最后由 编辑

      我是用湍流模型Komega求解流固耦合方腔自然对流(高宽比7.5),使用流固耦合求解器进行求解。温度方程中的alphat=turbulence.nut()/Prt. alphat使用zeroGradient或者Calculated边界条件,能够正常计算出来,但温度迭代极快完成,从高温壁面直接热传导到冷壁面。所以个人觉得应该使用标题所示的壁面函数,况且我在单域算例能够使用alphatJayatilleWallFunction壁面函数,结果也明显与零梯度边界条件不同,呈现壁面温度梯度高,中间温度均匀。现在我想在流固耦合中用alphatJayatilleWallFunction,就报如下错误。

      Region: fluid Courant Number mean: 0 max: 0
      Region: fluid Courant Number mean: 0 max: 0
      Time = 0.01
      
      Solving for fluid region fluid
      DILUPBiCG:  Solving for Ux, Initial residual = 0, Final    residual = 0, No Iterations 0
      DILUPBiCG:  Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0
      DILUPBiCG:  Solving for Uz, Initial residual = 0, Final residual = 0, No Iterations 0
      GAMGPCG:  Solving for p, Initial residual = 0, Final residual = 0, No Iterations 0
      time step continuity errors : sum local = 0, global = 0, cumulative = 1(0)
      DILUPBiCG:  Solving for omega, Initial residual = 0.3364115726, Final residual   = 9.041867288e-08, No Iterations 81
      DILUPBiCG:  Solving for k, Initial residual = 1, Final residual = 9.098095726e-08, No Iterations 158
      bounding k, min: -36.19837748 max: -14.87872514 average: -18.05165843
      [6] 
      [6] 
      
      [1] --> FOAM FATAL ERROR: 
      [1] 
      request for dictionary transportProperties from objectRegistry fluid failed
      available objects of type dictionary are
      
      5
      (
      RASProperties//
      fvSchemes
      fvSolution
      data
      turbulenceModel
      )
      [1] 
      [1] 
      [1]     From function objectRegistry::lookupObject<Type>   (const word&) const
      [1]     in file /home/kdd/OpenFOAM/OpenFOAM-2.3.1/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 198.
      [1] 
      FOAM parallel run aborting
      

      单域算例测试结果如下:
      使用零梯度边界条件

      d823d5e5-a257-44f2-ba92-dbed34af49b0-图片.png

      d45793eb-aefd-4b81-9c24-83b5f02aef2b-图片.png

      使用alphatJayatilleWallFunction壁面函数

      30760a26-a7fa-40b1-a817-70c9d3e60ce4-图片.png

      322e03f7-2512-477d-a906-3675b3d8841b-图片.png

      最后有两个疑问。1:湍流模型计算方腔自然对流的alphat应该使用哪个边界条件,哪个结果更合理。
      2:如果温度方程我不考虑湍流普朗特数Prt对热扩散系数的影响,直接使用固定值,然后使用零梯度计算(绕开这个alphatJayatilleWallFunction)。温度方程的浮升力影响动量方程,而湍流模型只用来求解动量方程,对温度方程不影响。是否仍然算作湍流自然对流。因为我目前看到的利用湍流模型来求解自然对流的文章中,都还极少提到Prt。是否意味着他们使用湍流模型计算自然对流也是不考虑Prt.

      可能问题较多,在这先谢谢各位。祝各位在身体健康,一夜暴富!
      最后附上我的算例与求解器(链接: https://pan.baidu.com/s/1m8LzT6iJU84UEq2qOp4iuw 提取码: pj6w)

      1 条回复 最后回复 回复 引用
      • 一
        一颗鸭蛋 最后由 编辑

        使用alphatJayatillWallFunction壁面函数的alphat结果。上面贴错了,不好意思。
        65665c8f-4fab-4812-9144-dbeb64bbb9bf-图片.png

        1 条回复 最后回复 回复 引用
        • D
          dzw05 教授 最后由 编辑

          你这是在耦合计算中没有找到transportProperties这个字典文件。可以仔细看一下,耦合计算和单场计算时,流体域是不是用的相同的求解器。

          自主匠心,普惠仿真。

          一 1 条回复 最后回复 回复 引用
          • 一
            一颗鸭蛋 @dzw05 最后由 编辑

            @dzw05 耦合计算和单域计算的流体域是使用上面提到的模型求解,求解器则是不同的,一个是壁面导电耦合求解器,一个壁面是绝缘求解器。我的耦合也只是壁面和流体的电势耦合。

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