热传导方程的差分解法公开课一等奖优质课大赛微课获奖课件.pptx
《热传导方程的差分解法公开课一等奖优质课大赛微课获奖课件.pptx》由会员分享,可在线阅读,更多相关《热传导方程的差分解法公开课一等奖优质课大赛微课获奖课件.pptx(29页珍藏版)》请在咨信网上搜索。
第九章 热传导方程差分解法第1页第1页n n9.1 9.1 热传导方程概述热传导方程概述 考虑三维空间温度改变情况考虑三维空间温度改变情况,设设 t t 时刻点时刻点(x,y,z)(x,y,z)处温度为处温度为u(x,y,z,t),u(x,y,z,t),则则 t t 时间内通过横截面积时间内通过横截面积 S S 传导热量为传导热量为(沿沿 n n 方向方向):):其中其中:K(x,y,z,t):K(x,y,z,t)是介质热传导系数是介质热传导系数,为温度梯度法向量分量为温度梯度法向量分量.取空间中一个小区域取空间中一个小区域 V,V,其边界面其边界面 S S 为一封闭曲面为一封闭曲面.则则 t t1 1 到到 t t2 2 时刻通过包面时刻通过包面 S S 传入传入 V V 热量为热量为:第2页第2页由高斯公式由高斯公式:n n为哈密顿算子为哈密顿算子:设介质比热容为设介质比热容为 c,c,密度为密度为 ,则则 V V 内温度改变消耗热量内温度改变消耗热量:设设 V V 内部热源密度为内部热源密度为 F(x,y,z,t),F(x,y,z,t),则内部热源产生热量为则内部热源产生热量为:第3页第3页依据能量守恒原则依据能量守恒原则:Q:Q2 2=Q=Q1 1+Q+Q3 3即即:亦即亦即:若若 F(x,y,z,t)F(x,y,z,t)0,c,0,c,K,K,为常数为常数,则则:第4页第4页其中其中:为拉普拉斯算子为拉普拉斯算子:因此热传导方程为因此热传导方程为:其中其中:K K c c.第5页第5页n n 9.2 一维热传导方程差分解法一维热传导方程一维热传导方程:n n初值问题初值问题初值条件初值条件:n n初边值混合问题初边值混合问题初值条件初值条件:边值条件边值条件:(:(关于边界点关于边界点x=0 x=0和和x=l)x=l)第一类第一类.第6页第6页n n第二类:n n第三类:n n其中g1(t),g2(t),1(t),2(t)为给定函数,要求1(t),2(t),且不同时为零.第7页第7页设空间步长为设空间步长为 h,h,时间步长为时间步长为 .把空间和时间离散化把空间和时间离散化:近似微分近似微分:故可定义故可定义:对空间一阶向前插商对空间一阶向前插商:第8页第8页对空间一阶向后插商对空间一阶向后插商:对空间二阶中心差商对空间二阶中心差商:对时间一阶向前插商对时间一阶向前插商:第9页第9页代入热传导方程代入热传导方程:迭代公式迭代公式:第10页第10页 t t i-1 i i+1 x k+1 k第11页第11页第一类初边值条件第一类初边值条件:已知已知:第12页第12页第二类初边值条件第二类初边值条件:已知已知即即:第13页第13页计算过程计算过程:第14页第14页第三类初边值条件第三类初边值条件:已知已知:即即:第15页第15页例例1:1:差分方程差分方程:初边值条件初边值条件:第16页第16页function u=rcd(lamda,tao,h,H,T)function u=rcd(lamda,tao,h,H,T)x=0:h:H;x=0:h:H;t=0:tao:T;t=0:tao:T;a=tao*lamda/h2;a=tao*lamda/h2;N=length(x);N=length(x);M=length(t);M=length(t);u(:,1)=(4*x.*(1-x);u(:,1)=(4*x.*(1-x);u(1,2:M)=0;u(1,2:M)=0;u(N,2:M)=0;u(N,2:M)=0;for k=1:M-1for k=1:M-1 for i=2:N-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);u(i,k+1)=a*u(i+1,k)+(1-2*a)*u(i,k)+a*u(i-1,k);end endendend第17页第17页h1=line(Color,1 0 h1=line(Color,1 0 0,Marker,.,MarkerSize,20,EraseMode,xor);0,Marker,.,MarkerSize,20,EraseMode,xor);for i=1:length(t)for i=1:length(t)set(h1,Xdata,0:0.1:1,Ydata,u(:,i);set(h1,Xdata,0:0.1:1,Ydata,u(:,i);pause(tao);pause(tao);endend第18页第18页X,Y=meshgrid(x,0:0.01:0.2);X,Y=meshgrid(x,0:0.01:0.2);Z=repmat(u(:,1),size(X,1),1);Z=repmat(u(:,1),size(X,1),1);h2=surface(X,Y,Z);h2=surface(X,Y,Z);shading interp,axis equal;shading interp,axis equal;set(h2,EraseMode,xor);set(h2,EraseMode,xor);for i=1:length(t)for i=1:length(t)CD=repmat(u(:,i),size(X,1),1);CD=repmat(u(:,i),size(X,1),1);set(h2,Cdata,CD);set(h2,Cdata,CD);pause(tao);pause(tao);endend第19页第19页n n9.3 二维热传导方程差分解法内部无热源均匀介质中二维热传导方程内部无热源均匀介质中二维热传导方程:初值条件初值条件:边值条件视详细情况而定边值条件视详细情况而定.设空间步长为设空间步长为 h,h,时间步长为时间步长为 .设设Nh=l,Mh=s,Nh=l,Mh=s,把把时间和空间离散化时间和空间离散化:第20页第20页即即:微分近似微分近似:第21页第21页代入热传导方程得代入热传导方程得:k+1ki-1 i i+1j+1jj-1第22页第22页例例:初值条件初值条件:即即:恒温边界恒温边界:绝热边界绝热边界:即即:第23页第23页差分公式差分公式:第24页第24页function u=rcd2(lamda,tao,h,T,L,S)function u=rcd2(lamda,tao,h,T,L,S)%二维热传导方程二维热传导方程t=0:tao:T;t=0:tao:T;x=0:h:L;x=0:h:L;y=0:h:S;y=0:h:S;a=tao*lamda/h2;a=tao*lamda/h2;if a0.25if a0.25 error(lamda*tao/h20.25);error(lamda*tao/h20.25);endend第25页第25页D=length(t);D=length(t);N=length(x);N=length(x);M=length(y);M=length(y);M1=ceil(M/2)-3;M1=ceil(M/2)-3;M2=ceil(M/2)+3;M2=ceil(M/2)+3;u=zeros(N,M,D);u=zeros(N,M,D);u(:,:,1)=0;u(:,:,1)=0;u(:,1,:)=0;u(:,1,:)=0;u(:,M,:)=0;u(:,M,:)=0;u(1,M1:M2,2:D)=1;u(1,M1:M2,2:D)=1;第26页第26页for k=1:D-1for k=1:D-1 for i=2:N-1 for i=2:N-1 for j=2:M-1 for j=2:M-1 u(i,j,k+1)=(1-4*a)*u(i,j,k)+a*(u(i+1,j,k)+u(i-u(i,j,k+1)=(1-4*a)*u(i,j,k)+a*(u(i+1,j,k)+u(i-1,j,k)+u(i,j+1,k)+u(i,j-1,k);1,j,k)+u(i,j+1,k)+u(i,j-1,k);u(N,j,k+1)=u(N-1,j,k+1);u(N,j,k+1)=u(N-1,j,k+1);if(jM2)if(jM2)u(1,j,k+1)=u(2,j,k+1);u(1,j,k+1)=u(2,j,k+1);end end end end end endendend第27页第27页X,Y=meshgrid(x,y);X,Y=meshgrid(x,y);Z=u(:,:,1);Z=u(:,:,1);h2=surface(X,Y,Z);h2=surface(X,Y,Z);shading interp,axis equal;shading interp,axis equal;set(h2,EraseMode,xor);set(h2,EraseMode,xor);for i=1:length(t)for i=1:length(t)CD=u(:,:,i);CD=u(:,:,i);set(h2,Cdata,CD);set(h2,Cdata,CD);pause(tao);pause(tao);endend第28页第28页rcd2(1,0.001,0.1,0.2,2,1);rcd2(1,0.001,0.1,0.2,2,1);第29页第29页- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 热传导 方程 解法 公开 一等奖 优质课 大赛 获奖 课件
咨信网温馨提示:
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,个别因单元格分列造成显示页码不一将协商解决,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【w****g】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【w****g】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,个别因单元格分列造成显示页码不一将协商解决,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【w****g】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【w****g】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。
关于本文