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



  • 我是用湍流模型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)



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


  • OpenFOAM教授

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



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



CFD中文网 2016 - 2020 | 京ICP备15017992号-2