计算热物理大作业.doc

计算热物理大作业报告4.3 1.对原来的方程进行无量纲化处理 用无量纲方程化为： 初始条件化为： ， 2，对方程进行离散 （a）控制容积法 全隐式格式为： 其中： ， 显式格式为：， 其中： （b）ADI格式离散 3.Fourier分析稳定性 控制容积全隐式格式 可以求出放大因子为：，可见隐式格式稳定 全显式格式的放大因子为：，稳定 ADI格式的放大因子为， ，所以ADI格式稳定 4.离散格式分析，各种离散格式源程序附后 各种迭代格式的迭代次数：从表中可以出看出，就收敛速度而言： SOR点迭代G-S点迭代简单点迭代 ADI SOR线迭代G-S线迭代简单线迭代 收敛最快地当属ADI计算方法。松驰迭代的收敛收敛快慢与松驰因子有很大关系，选择最佳松弛因子可以大大加快收敛速度。 5.选取了G-S点迭代的不同时刻的图线 图中分别选择温度变化比较剧烈的5个时刻做等值线图 1， 图 SEQ 图 \* ARABIC 1 初始状态 图 SEQ 图 \* ARABIC 2时刻2 图 3 时刻3 图 4时刻4 图 5 时刻5 温度在开始的时候收敛比较快，接近收敛区时收敛比较慢，可以看出，等温线逐渐向高温度区偏移，证