华中科技大学能源与力工程学院.ppt

------------- 第三章 热传导问题的计算 第一节 一维稳态传热问题 二、控制方程 2、方程离散化 c、线性方程组 四、求解方法 五、计算实例1 计算结果比较 六、进一步的思考 2、界面的导热系数 3、界面的温度梯度 4、相邻点系数之间的关系 七、计算实例2 八、最后的细节 2、非均匀网格 3、热流量边界条件 4、流体换热边界条件 总 结 第二节 非稳态一维导热问题 二、方程的离散化 对时间的插值函数 三、显式与隐式模式 不同模式的物理意义 关于计算稳定性的一个示例 当导热系数与温度有关时的迭代过程 四、算例3 计算结果 第三节 二维导热问题 二、边界条件 第一种网格边界条件的处理 3、流体换热边界条件 第二种网格边界条件的处理 周期与对称边界条件 周期边界条件的处理 对称边界条件 三、算例4 算例4（续） 还有什么问题没有解决？ 第四节 线性方程组的求解 二、方程组的求解 收敛及收敛条件 收敛过程（1） 收敛过程（2） 松弛因子 第二章 总结 ------------- 华中科技大学能源与动力工程学院 周期性边界条件 ------------- 华中科技大学能源与动力工程学院 PR PL PL PR ------------- 华中科技大学能源与动力工程学院 ------------- 华中科技大学能源与动力工程学院 边界条件如图 ?＝150 k＝50 绝热 100? 400? 对流换热 Tf=50? ------------- 华中科技大学能源与动力工程学院 ------------- 华中科技大学能源与动力工程学院 线性方程组的解法 复杂边界形状时网格的生成 ------------- 华中科技大学能源与动力工程学院 一、方程组的特点 a11 a12 a13 a14 a15 a16 …… a21 a22 a23 a24 a25 a26 …… a31 a32 a33 a34 a35 a36 …… a41 a42 a43 a44 a45 a46 …… a51 a52 a53 a54 a55 a56 …… a61 a62 a63 a64 a65 a66 …… …… …… …… …… T1 T2 T3 T4 T5 T6 … = b1 b2 b3 b4 b5 b6 … 系数矩阵 变量 常数 系数矩阵为大型稀疏矩阵 ------------- 华中科技大学能源与动力工程学院 T1=( aE1T2+b1)/aP1 T2=(aW2T1+aE2T3+b2)/aP2 T3=(aW3T2+aE3T4+b3)/aP3 T4=(aW4T3+aE4T5+b4)/aP4 ……… 高斯－赛德尔迭代法 解法： 消去法 迭代法 存储量＝k×n2 存储量＝k×n ------------- 华中科技大学能源与动力工程学院 收敛：所得的解逐渐接近精确解 对所有方程 收敛条件： 对至少一个方程 残差：方程两端数值的差 ------------- 华中科技大学能源与动力工程学院 * ------------- 华中科技大学能源与动力工程学院 一、物理模型 t q 无限大平板(一维） 两面温度恒定 材料导热系数恒定 ------------- 华中科技大学能源与动力工程学院 三、离散化方法 x P W E w e ?xw ?xe ?x 1、计算域离散化 控制体积的划分 有限体积法 单元（cell） ------------- 华中科技大学能源与动力工程学院 a、单元内积分 x P E w e ?x ?xw ?xe b、插值函数 T w e E W P T w e E W P W ------------- 华中科技大学能源与动力工程学院 ------------- 华中科技大学能源与动力工程学院 1、边界条件 迭代法： 假定T，得到k 计算T，重算k 再计算T…… 2、非线性问题 3、方程求解 ------------- 华中科技大学能源与动力工程学院 1、题目： 板厚 ?＝8 T1=300 T2=100 k= 0.15(1+?T) ?＝0.3 S=0 2、计算条件： 单元数8 温度边界条件 3、解析解 ------------- 华中科技大学能源与动力工程学院 ------------- 华中科技大学能源与动力工程学院 1、数学近似与物理真实 数学近似方法必须符合物理规律 所有符合物理规律的数学近似方法都是可行的 ------------- 华中科技大学能源与动力工程学院 算术平均法 P E ?xe- ?xe+ e 几何平均法—从热流量平衡出发