流固耦合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单域算例测试结果如下：
使用零梯度边界条件
使用alphatJayatilleWallFunction壁面函数
最后有两个疑问。1：湍流模型计算方腔自然对流的alphat应该使用哪个边界条件，哪个结果更合理。
2：如果温度方程我不考虑湍流普朗特数Prt对热扩散系数的影响,直接使用固定值，然后使用零梯度计算（绕开这个alphatJayatilleWallFunction）。温度方程的浮升力影响动量方程，而湍流模型只用来求解动量方程，对温度方程不影响。是否仍然算作湍流自然对流。因为我目前看到的利用湍流模型来求解自然对流的文章中，都还极少提到Prt。是否意味着他们使用湍流模型计算自然对流也是不考虑Prt.可能问题较多，在这先谢谢各位。祝各位在身体健康，一夜暴富！
最后附上我的算例与求解器(链接: https://pan.baidu.com/s/1m8LzT6iJU84UEq2qOp4iuw 提取码: pj6w)
