24小时热门版块排行榜

>论坛更新日志 (4069)
>文献求助 (423)
>虫友互识 (350)
>导师招生 (339)
>考博 (176)
>博后之家 (168)
>硕博家园 (156)
>休闲灌水 (143)
>论文投稿 (121)
>基金申请 (91)
>招聘信息布告栏 (57)
>考研 (57)
>教师之家 (54)
>公派出国 (43)
>绿色求助(高悬赏) (38)
>找工作 (37)

返回列表

crystal_ling

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 1456.1
红花: 2
帖子: 80
在线: 51.3小时
虫号: 1977534
注册: 2012-09-05
专业: 光学

[求助] 如何用MATLAB的偏微分方程工具箱画出不同区域内的温度分布图？

[求助]求助：如何用MATLAB的偏微分方程工具箱画出不同区域内的温度分布图？
初始悬赏金币 5 个
求助：MATLAB的偏微分方程工具箱问题
问题：画出下图所示模型内的温度分布图。

稳定状态下的温度分布通过数值求解标准的热传导方程得到
-▽﹒(k▽T)=Q
其中各层的热载Q及参数如下：
W层：
Q1=(2*n1*P*a1/pi*w^2)exp(-2*r^2/w^2)exp(-a1(z01-z1))
n1=1, P=10W, w=50μm, a1=0.242μm^-1,  z01=0μm
kz1=0.000015W/(μm﹒K), kr1=0.000013W/(μm﹒K)
MQW层：
Q2=(2*n2*P*a2/pi*w2)exp(-2*r^2/w^2)exp(-a2(z02-z2))
n2=0.2231, P=10W, w=50μm, a2=0.925μm^-1,  z02=-0.627μm
kz2=0.000028W/(μm﹒K), kr2=0.000028W/(μm﹒K)
DBR层：
Q3=(2*n3*P*a3/pi*w^2)exp(-2*r^2/w^2)exp(-a3(z03-z3))
n3=1, P=10W, w=50μm, a3=0.21μm^-1,  z03=-2.931μm
kz3=0.000069W/(μm﹒K), kr3=0.00006W/(μm﹒K)

我的具体做法如下：
第一步：打开MATLAB偏微分工具箱，选模型：Heat Transfer
第二步：画三个矩形分别代表window层、MQW层、DBR层
第三步：设定边界条件:上、左、右绝缘，所以设定dT/dn=0，下层T=300K
第四步：分段输入k值和Q值。（如何输入具体的值？不太会）

但画出的图温度分布不对，可能是不太会使用偏微分工具箱，恳求哪位高人能指点，最好将每一步操作步骤列出来，小女子将不胜感激！

模型图.jpg

[ Last edited by 华丽的飘过 on 2012-12-7 at 00:47 ]

回复此楼