有限差分和matlabpde求解一维稳态传热问题甄选.docxVIP

精品 精品word. PAGE PAGE # / 10doc可编辑 精品 精品word. PAGE PAGE # / 10doc可编辑 有限差分和Matlabpde求解一维稳态传热问题?（优选） 有限差分和pde函数求解一维定态热传导方程 初始分别用有限差分方法和pde函数求解一维定态热传导方程, 条件和边界条件，热扩散系数 a=0.00001, 初始 2T 7 (1) 求解过程: 用Tylaor展开法推导出FTCS格式的差分方程 首先对T进行泰勒展开得到如下两式子: Tjn1 T￡2!2t Tjn1 T ￡ 2! 2t t2 3? 3t t3 Tjn1 Tn 2X2? 2T ~2 X 3 X 3! 上述两个方程变换得: t2工t t2 工 t3 Tjn1 Tjn t T；(2)2 T； (2) 2t X2 3t X3 2t7x2X3! 2t 7 x2 X 3! 3t 2 4! X2 4t X7 精品word. 精品word. 1 PAGE PAGE # / 10doc可编辑 精品word. 精品word. 1 PAGE PAGE # / 10doc可编辑 精品word. 精品word. (3) PAGE PAGE # / 10doc可编辑 将上述式子2 Tjn1 2Tjn 将上述式子 2 Tjn1 2Tjn Tjn1 2 x 2T 2 2 0 X (2)( 3)代入(1)得: H n Tjn1 Tjn x2 j t Tn1 2Tn Tn1 j 1 j j 1 2 x 0( t, x2) (4) 方程的相容性和稳定性讨论: 上述方程截项为: 0( t, x)由于I肌 0( t, x) 由于I肌0 n n t T t T / 2 t2 j 3! t3 j … t, x 0所以方程有相容性 x 3! 3T 7 2 x 4! :T n x4 j 其经过傅里叶变换后，只有当t, x满足下列条件时，方程具有 较好的稳定性: km x 其中km -LN N为节点个数，L为边界长度由于: sin2 sin2 J 2L 所以当s — 0.5时方程具有稳定性 x 说明该方程的类型和定解条件，如何在程序中实现这些定解条 件。 由式(4)得： TjTj Tj Tj T-n1 2Tn Tn 1 t j i j j i t 2 x sTjn1 (1 2s)Tjn sTjn1 要想知道下一个时刻tn1下的温度，则必须知道在tn时刻下该点 以及相邻两点的温度，如下图所示: 所以需要知道的是在各个时刻X轴两侧的温度，以及初始时刻 X轴上各点温度。下述各题求解过程中，均采用 0.00001, Tj0 0(j 1,2... jmax 1) Tn 100，作为边界条件。边界 条件是由外界给定，本身对过程求解，精度无影响 编写M文件求解上述方程 M文件的编写见附录所示。 在 matlab 中调用 heat_conduct(0.00001,15000,500,0.1),得到如 下结果： rms error is rms = 0.7959 精品 精品word. PAGE PAGE # / 10doc可编辑 精品 精品word. PAGE PAGE # / 10doc可编辑 3000b900Q,15000a1刻的温度分布 100 90 SO 70 60 板50 理 40 30 20 10 0 ）（方向说明：1?上述结果采用的x, t分别为0.1和500,时间最终给定为 ）（方向 15000 当x, t为（0.01,5）时的结果 精品 精品word. PAGE PAGE # / 10doc可编辑 PAGE PAGE # / 10doc 可编辑 整体分布0 05 1 整体分布 0 05 1 夏方向 3000,9000.1500001刻的温廈分布 100 90 eo 70 60 fe 50 40 30 20 10 0 附录1，传热方程有限差分的 M文件编写 fun cti on heat_c on duct(alpha,Timemax,dt,dx) %方程需要输入的参数有，最大时刻数，alpha, dt, dx,实 际 Tmax二Timemax/dt+1;%时刻数 JMAX=1/dx+1;%x轴上的差分点个数 J=JMAX-1;MAXEX=J; T=zeros(Tmax,JMAX);构造温度矩阵 TE二zeros(Tmax,JMAX);实际温度矩阵 %设定边界条件 T(:,1)=100;T(:,JMAX)=100; s=alpha*dt/dx/dx; 精品 word. %计算各个点的温度 for i=1:(Tmax-1) for j=2:J T(i+1,j)=s*T(i,j+1)+(