我真是没招了，有没有高人拯救我一下。这个模型本身就很简单，一个普通的圆柱腔，顶部低温壁面底部高温壁面，周围绝热，k-ωSST模型，Ra在1×10⁸～1×10⁹。这残差不管迭代多少次都下不去，调了半个月了一直这样，实在搞不懂了。右图是能量守恒后的竖直截面上温度分布结果，也就是上下热流密度相同我就当它收敛了，毕竟残差一直这样我实在没办法。这个结果和论文对比了一下整体分布趋势是对的，但中间的温度略高。但这些都不是最主要的，最头疼的是当我改变重力，比如从9.8改到1.6，虽然热流密度变了，但这个温度分布几乎没有变化。理论上来说重力减小应该更接近纯导热，温度梯度应该更明显。这个结果还没办法验证是不是对的，毕竟没有微重力环境下的实验。所以我现在有两个最主要的问题：
1.自然对流残差一直是这样合理吗
2.为什么重力减小温度分布几乎不变，是因为结果本身就误差太大还是说可能有其他原因

1.不合理，这是在错误解上残差收敛，建议优化壁面网格。2.没开boussinesq加热浮力同时没外部施加的动量源项下，不可压缩流调整重力本身就不会影响流场整体趋势，你整个多相流模型调整重力估计会有说法。先确认模型是否考虑了热浮力，随后确认热膨胀系数选取以及工作条件的boussinesq参考温度是否合理（正常Tref取材料液相线温度，没”改的话默认是不到300k），如果热体积膨胀系数极小（低于e-6量级）同时没外部动量源项作用，重力10跟1变化确实不大，因为本身就是热导率和htc在发力，而不是对流在发力

这玩意儿得用Fluent模拟，调好边界条件，多看看教程。