- 1、本文档共16页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
偏微分方程组解法
偏微分方程组解法
某厚度为10cm平壁原温度为20,现其两侧面分别维持在20和120,试求经过8秒后平壁内温度分布,并分析温度分布随时间的变化直至温度分布稳定为止。
式中为导温系数,;。
解:
模型转化为标准形式:
初始条件为:
边界条件为:
,
函数:
pdefun.m
%偏微分方程(一维动态传热)
function [c,f,s]=pdefun(x,t,u,dudx)
c=1/2e-4;f=dudx;s=0;
icbun.m
%偏微分方程初始条件(一维动态传热)
function u0=icbun(x)
u0=20;
bcfun.m
%偏微分方程边界条件(一维动态传热)
function [pl,ql,pr,qr]=bcfun(xl,ul,xr,ur,t)
pl=ul-120;ql=0;pr=ur-20;qr=0;
命令:
x=linspace(0,10,20)*1e-2;
t=linspace(0,15,16);
sol=pdepe(0,@pdefun,@icfun,@bcfun,x,t);
mesh(x,t,sol(:,:,1)) %温度与时间和空间位置的关系图
%画1、2、4、6、8、15s时刻温度分布图
plot(x,sol(2,:,1)) 1s时刻,(因为本题sol第一行为0时刻)
hold on
plot(x,sol(3,:,1))
plot(x,sol(5,:,1))
plot(x,sol(7,:,1))
plot(x,sol(9,:,1))
plot(x,sol(16,:,1))
计算结果:
%第8秒时温度分布
x
sol(9,:,1)
经过8秒时的温度分布为:
/cm 0 0.5263 1.0526 1.5789 2.1053 2.6316 3.1579 / 120.0000 112.5520 105.1653 97.8994 90.8100 83.9477 77.3562 /cm 3.6842 4.2105 4.7368 5.2632 5.7895 6.3158 6.8421 / 71.0714 65.1202 59.5200 54.2784 49.3930 44.8518 40.6338 /cm 7.3684 7.8947 8.4211 8.9474 9.4737 10.0000 / 36.7095 33.0419 29.5877 26.2982 23.1207 20.0000 或者求第8秒时,x=0,2,4,,6,8,10cm处的温度
[uout,duoutdx]=pdeval(0,x,sol(9,:,:),[0,2,4,6,8,10]*1e-2)
120.0000 92.2279 67.5007 47.5765 32.3511 20.0000
不同时刻温度分布图
将上图的视角转至xt平面也得到本图,从本图可知当时间达到15s时平壁内的温度分布已近稳定。
某厚度为20cm钢板原温度为20,现将其置于1000的炉中加热,平壁导热系数为,对流传热系数,导温系数为;试分析温度分布随时间的变化及钢板表面温度达到500时所需的时间。
解:
模型转化为标准形式:
初始条件为:
边界条件为:
(平壁中心坐标为0,绝热),
函数:
pdefun1.m
%偏微分方程(一维动态平壁两侧对流)
function [c,f,s]=pdefun1(x,t,u,dudx)
c=1/0.555e-5;f=dudx;s=0;
icbun1.m
%偏微分方程初始条件(一维动态平壁两侧对流)
function u0=icbun1(x)
u0=20;
bcfun1.m
%偏微分方程边界条件(一维动态平壁两侧对流)
function [pl,ql,pr,qr]=bcfun1(xl,ul,xr,ur,t)
pl=0;ql=1;pr=174*(ur-1000);qr=34.8; %平壁两侧置于同一流体中具有对流传热,平壁中心为绝热
命令:
%600s内的温度分布变化
x=linspace(0,10,20)*1e-2;
t=[0:60:600];
sol=pdepe(0,@pdefun1,@icfun1,@bcfun1,x,t);
mesh(x,t,sol)
%2160s内的温度分布变化
t=[0:60:2160];
sol=pdepe(0,@pdefun1,@icfun1,@bcfun1,x,t);
mesh(x,t,sol)
%60、120、180、240、300、360、420s时刻温度分布图
plot(x,sol(2,:,1)) 60s(t网格为0:60:2160,其时间间隔为0,60,120,180,……2160,第二点为60s)
文档评论(0)