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

第九章 热传导方程的差分解法 第1页，共28页。 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 的热量为: 第2页，共28页。 由高斯公式: 为哈密顿算子: 设介质的比热容为 c, 密度为 ?, 则 V 内温度变化消耗的热量: 设 V 内部热源密度为 F(x,y,z,t), 则内部热源产生的热量为: 第3页，共28页。 根据能量守恒原则: Q2 = Q1 + Q3 即: 亦即: 若 F(x,y,z,t)?0, c, ?, K, 为常数,则: 第4页，共28页。 其中: ? 为拉普拉斯算子: 所以热传导方程为: 其中: ? ? K?c?. 第5页，共28页。 9.2 一维热传导方程的差分解法 一维热传导方程: 初值问题 初值条件: 初边值混合问题 初值条件: 边值条件:(关于边界点x=0和x=l) 第一类. 第6页，共28页。 第二类: 第三类: 其中g1(t), g2(t), ?1(t), ?2(t) 为给定函数, 要求?1(t)??, ?2(t) ??, 且不同时为零. 第7页，共28页。 设空间的步长为 h, 时间的步长为 ?. 把空间和时间离散化: 近似微分: 故可定义: 对空间一阶向前插商: 第8页，共28页。 对空间一阶向后插商: 对空间二阶中心差商: 对时间一阶向前插商: 第9页，共28页。 代入热传导方程: 迭代公式: 第10页，共28页。 t i-1 i i+1 x k+1 k 第11页，共28页。 第一类初边值条件: 已知: 第12页，共28页。 第二类初边值条件: 已知 即: 第13页，共28页。 计算过程: 第14页，共28页。 第三类初边值条件: 已知: 即: 第15页，共28页。 例1: 差分方程: 初边值条件: 第16页，共28页。 function u = rcd(lamda,tao,h,H,T) x = 0:h:H; t = 0:tao:T; a = tao*lamda/h^2; N = length(x); M = length(t); u(:,1) = (4*x.*(1-x)); u(1,2:M) = 0; u(N,2:M) = 0; for k=1:M-1 for i=2:N-1 u(i,k+1)=a*u(i+1,k)+(1-2*a)*u(i,k)+a*u(i-1,k); end end 第17页，共28页。 h1=line(Color,[1 0 0],Marker,.,MarkerSize,20,EraseMode,xor); for i=1:length(t) set(h1,Xdata,[0:0.1:1],Ydata,u(:,i)); pause(tao); end 第18页，共28页。 [X,Y]=meshgrid(x,[0:0.01:0.2]); Z = repmat(u(:,1),size(X,1),1); h2 = surface(X,Y,Z); shading interp,axis equal; set(h2,EraseMode,xor); for i=1:length(t) CD = repmat(u(:,i),size(X,1),1); set(h2,Cdata,CD); pause(tao); end 第19页，共28页。 9.3 二维热传导方程的差分解法 内部无热源均匀介质中二维热传导方程: 初值条件: 边值条件视具体情况而定. 设空间的步长为 h, 时间的步长为 ?. 设Nh=l, Mh=s, 把时间和空间离散化: 第20页，共28页。