一维热传导方程的差分格式.docxVIP

一维热传导方程的差分格式 《微分方程数值解》 学名1：徐惠清学号名2：项羽学号名3：邱文林学号名4：高军学号名5：赵玉恒学号号6：刘志刚学号院：理学院专业：级别14信息与计算科学讲师：陈洪斌 2021年6月25日 一维热传导方程的差分格式 一、《微分方程数值解》课程论文的格式1）引言：介绍研究问题的意义和现状2）格式：给出数值格式 3） 截断误差：给出数值格式的截断误差4）数值示例：根据给定的数值格式给出数值示例5）参考文献：论文涉及的文献和教材 二、《微分方程数值解》课程论文的评分标准1）文献综述：10分； 2） 研究方案可行性：10分；3） 数字格式：20分； 4）数值格式的算法、流程图：10分；5）数值格式的程序：10分； 6） 论文写作的组织性和完整性：10分；7） 论文工作量和题目难度：10分；8） 课程设计态度：10分；9） 自主创新：10分。 考虑如下一维非齐次热传导方程dirichlet初边值问题 =a2+f（x，t），c u(x,0)=?(x),c≤x≤d,(1.2) u（c，t）=α（t），u（d，t）=β（t），0 的有限差分方法,其中a为正常数,f(x,t),?(x),α(t), β（t）是一个已知常数，φ（c）=α（0）， ?(d)=β(0).称(1.2)为初值条件,(1.3)为边值条件. 本文给出了（1.1）（1.3）的前向Euler格式、后向Euler格式和crank-Nicolson格式，并通过比较给出了它们的截断误差和数值例子，发现crank-Nicolson格式误差最小，前向 euler格式次之,向后euler格式误差最大. 2差分方案的建立 2.1向前euler格式 将区间[C，D]分为m，[0，t]分为N，并写出H=（D-C）/m，τ=t/N， xj=c+jh,0≤j≤m,tk=kτ,0≤k≤n.分别称h和τ为空间步长和时间步长.用 x=xj，0≤J≤m、 t=tk，0≤K≤N 将ω分割成矩形网格.记ωh=xj|0≤j≤m,ωτ={tk|0≤k≤n},ωhτ=ωh?ωτ.称xj,tk为结点. ω=Uj | 0上ωhτ网格函数的定义≤ J≤ m、 0≤ K≤ n、 其中UJ=uxj，TK 在结点xj,tk处考虑方程(1.1),有 ?u（xj，tk）=a ?2u(xj,tk)+f(xj,tk),1≤j≤m-1,1≤k≤n-1.(2.1) 以u（XJ，TK+1）作为节点（XJ，TK） 为中心关于t运用泰勒级数展开,有 u（xu（xj，tk）τ2 j,tk+1)=u(xj,tk)+u(xj,tk)τ+ +o（τ2）。 u(x2j,tk+1)-u(xj,tk) =?u（xj，tk）τ?u（xj，tk）τ ?t+2?t +o（τ）。（2.2）分别将（XJ-1，TK）、（XJ+1，TK）作为节点（XJ，TK） 为中心关于x运用泰勒级数展开,有 u（xj，tk（）-h） j-1,tk)=u(xj,tk)+u(xj,tk)(-h) j、 tk（-h）+ u(x2! u（x3 u(4)(xj,tk)(-h) +o（h4），（xu（xu）（xj，tk）h2 u(xj,tk)h3 uj+1，tk）=j，tk）+u（xj，tk）h+ u(4)(xj,tk)h4 +o（h4）。 由上述两式可得 u（xj-1，tk）-2u（xj，tk）+u（xj+1，tk）?2u（xj，tkh2=）?x2+h2?4（xj，tk）12?十、 +o(h2 )（2.3）将方程（2.2）和（2.3）代入（2.1）以获得 u(xj,tk+1)-u(xj,tk) k） +u（xj+1，tk） u(xj-1,tk)-2u(xj,th +f（xj，tk）+rkj。(2.4) 2其中rkah?4u(xj,tk ) τ 2?2（xj，tk）j =-12?x4+2?t2 +o（τ+H2）是方程式（2.1）的截断误差，四舍五入截断误差并使用UK j代替u(xj,tk) ，得到以下差分方程 uk+1-ukjj ukkk