一种求解瞬态热传导方程的无条件稳定方法

摘要

瞬态热传导问题普遍存在于航空航天、土木和冶金等领域中,对这类问题精确、高效的数值求解方法一直备受关注. 为此,本文提出了一种无条件稳定的单步时间积分方法. 在所提出的方法中,拉格朗日插值函数被用来近似真实的温度场及其一次导数场,之后,加权残量法被用来建立二者之间的关系. 通过对算法参数的巧妙设计,本文提出的方法具有二阶精度和L型数值耗散,其中,L型耗散使得本文方法能够快速过滤掉虚假的高频振荡. 多数现有时间积分方法仅对线性系统具有无条件稳定性,对非线性系统则是条件稳定的. 为此,本文改进了Hughes针对一阶非线性热传导问题提出的时间积分方法稳定性评估理论,并将改进的理论用于方法的参数设计中. 理论分析的结果表明本文方法对线性和非线性热传导系统均是无条件稳定的. 即使对于著名的Crank-Nicolson方法失稳的非线性热传导问题,本文方法仍能给出稳定且精确的预测. 数值测试结果显示,所提出的方法相较于当前流行的方法具有明显的精度、耗散和稳定性优势.