24小时热门版块排行榜    

查看: 3579  |  回复: 5

crystal_ling

铁虫 (小有名气)

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

求助: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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shizijun

木虫 (著名写手)

自定义头衔

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
xiegangmai: 金币+2, 谢谢提示! 2012-09-17 18:03:33
提供一个建议:
先化简一下Q,替代r 和z为x 和y;
然后输入就可以了
个人签名
2楼2012-09-17 10:27:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

crystal_ling

铁虫 (小有名气)

引用回帖:
2楼: Originally posted by shizijun at 2012-09-17 10:27:05
提供一个建议:
先化简一下Q,替代r 和z为x 和y;
然后输入就可以了

画图时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)  后Q还是一个表达式,r和z是用x和y代替,但画出的温度分布图还是不对,也检查不出哪儿错了,能否详细告诉我在MATLAB偏微分工具箱中如何输入具体的值。
3楼2012-09-17 10:39:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shizijun

木虫 (著名写手)

自定义头衔

【答案】应助回帖

xiegangmai: 谢谢参与!非应助请勿选择“应助回帖” 2012-09-19 23:39:26
引用回帖:
3楼: Originally posted by crystal_ling at 2012-09-17 10:39:23
画图时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)  后Q还是一个表达式,r和z是用x ...

麻烦你把生成的M文件上传上来。大家研讨研讨
个人签名
4楼2012-09-17 11:36:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

crystal_ling

铁虫 (小有名气)

有没有人可以帮帮我,非常感谢
5楼2012-09-19 17:22:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
xiegangmai: 专业版交流贴,请勿灌水! 2012-09-25 07:06:29
6楼2012-09-24 11:07:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 crystal_ling 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学调剂0703 +7 啊我我的 2026-03-11 7/350 2026-03-15 23:03 by 凌千颂111
[文学芳草园] 伙伴们,祝我生日快乐吧 +15 myrtle 2026-03-10 24/1200 2026-03-15 21:16 by 苏州_逗号
[考研] 267一志愿南京工业大学0817化工求调剂 +5 SUICHILD 2026-03-12 5/250 2026-03-14 14:53 by jean5056
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 求调剂 +5 鹤遨予卿 2026-03-09 5/250 2026-03-14 02:44 by JourneyLucky
[考研] 材料与化工 一志愿山大 321分 求调剂 +7 每天散步 2026-03-09 8/400 2026-03-14 02:18 by JourneyLucky
[考研] 一志愿天津大学,英一数二305分求调剂,四六级已过 +8 小小番的茄 2026-03-09 8/400 2026-03-14 01:53 by JourneyLucky
[考研] 341求调剂 +4 番茄头--- 2026-03-10 4/200 2026-03-13 23:12 by JourneyLucky
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 工科,求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[考研] 310求调剂 +3 【上上签】 2026-03-11 3/150 2026-03-13 16:16 by JourneyLucky
[考研] 考研调剂 +4 芬达46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[考研] 285求调剂 +4 ytter 2026-03-12 4/200 2026-03-13 14:48 by jxchenghu
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生,2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
[考研] 341求调剂 +4 捣蛋猪猪 2026-03-11 4/200 2026-03-12 14:47 by ruiyingmiao
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
[考研] 一志愿:武汉理工,材料工程,英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
信息提示
请填处理意见