格林函数法在计算带有非齐次边界条件或者内热源项的瞬态导热问题时具有非常好的适用性,但在使用时需要应对解析解中的反常积分项。对于图1所示的瞬态导热模型,研究人员首先通过格林函数法求出了温度函数的精确表达式(公式1),再将其中的积分项替换为求和项,优化了精确解的计算速度并保证了计算精度。

其中:



图1. 三维柱状瞬态导热模型切面示意图
为了验证格林函数法所得精确解的可行性,使用此精确解分别计算了壁面温度在轴向和周向上随时间变化的情况,并与此前所得准确结果进行了对比,结果如图2和图3所示。可以看出温度在瞬态响应初期与理论结果基本完全吻合。此模型中轴向热流区域选取 ,在轴向的温度计算中,当 时温度处于峰值状态,并逐步向未加热区域扩散,而随着时间的增加,在加热区与非加热区的温度也逐步升高。以上结论证明了所得解析解的准确性。

图2. 无量纲温度在轴向上随时间的变化

图3. 无量纲温度在周向上随时间的变化
为了验证格林函数法所得结果可以替换CFD数值计算方法,图4中计算并对比了两种方法所得无量纲温度在周向上的分布情况。其中CFD计算网格分别选取了coarse、medium和fine三种不同程度的细化方式。而解析解精度分别为 =10-2,10-3,10-4。当两种方法的精度都较低时,温度的变化趋势受到计算中吉布斯现象的影响呈现上下浮动的状态,而随着精度的不断提升,计算结果的收敛性也随之改善。结果证实了格林函数法与数值计算方法在高精度条件下具有很好的吻合性。虽然格林函数法不能处理比较复杂的几何结构,但是对于简单规则的几何模型进行瞬态导热分析时,可以使用这种方法替代CFD数值方法,避免重复性的网格划分及边界条件设置,在保证精度的前提下提高计算效率。同时,在没有实验结果支撑时,格林函数法可以作为一种有效的方法对数值计算结果进行验证。相关结论可以广泛应用于航空、材料焊接、地热等多个领域。本项研究受到国家重点研发项目 (No. 2016YFB0901402) 和国家自然科学基金 (No.51776198) 的支持,相关成果已发表在期刊Journal of Heat Transfer ASME上。

图4. 精确解与不同网格数CFD计算温度所得结果在无量纲时间 沿周向分布情况对比