第九章 热传导方程差分解法.pptVIP

第九章 热传导方程的差分解法 9.1 热传导方程概述 考虑三维空间的温度变化情况, 设 t 时刻点(x,y,z)处的温度为u(x,y,z,t), 则 ?t 时间内通过横截面积 ?S 传导的热量为(沿 n 方向): 其中: K(x,y,z,t) 是介质的热传导系数, 为温度梯度的法向量分量. 取空间中的一个小区域 V, 其边界面 S 为一封闭曲面. 则 t1 到 t2 时刻通过包面 S 传入 V 的热量为: 由高斯公式: 为哈密顿算子: 设介质的比热容为 c, 密度为 ?, 则 V 内温度变化消耗的热量: 设 V 内部热源密度为 F(x,y,z,t), 则内部热源产生的热量为: 根据能量守恒原则: Q2 = Q1 + Q3 即: 亦即: 若 F(x,y,z,t)?0, c, ?, K, 为常数,则: 其中: ? 为拉普拉斯算子: 所以热传导方程为: 其中: ? ? K?c?. 9.2 一维热传导方程的差分解法 一维热传导方程: 初值问题 初值条件: 初边值混合问题 初值条件: 边值条件:(关于边界点x=0和x=l) 第一类. 第二类: 第三类: 其中g1(t), g2(t), ?1(t), ?2(t) 为给定函数, 要求?1(t)??, ?2(t) ??, 且不同时为零. 设空间的