应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 求助：如何用MATLAB的偏微分方程工具箱画出不同区域内的温度分布图？
61/1返回列表
查看: 3817  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-09-16 14:37:12
已阅   回复此楼   关注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的回帖

haoyali0311
xiegangmai: 专业版交流贴，请勿灌水！ 2012-09-25 07:06:29
6楼2012-09-24 11:07:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 crystal_ling 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 生命口面上会评结束了吗大家有知道结果吗？现在还没有结果正常不？ +3 碧水00 2026-07-16 5/250 2026-07-16 07:39 by jaasxulu
[基金申请] 生生命学部会评 +4 wuke100666 2026-07-15 4/200 2026-07-16 00:12 by 啦啦小妮子
[有机交流] 伯胺上甲基求助 100+4 002tmumu 2026-07-14 4/200 2026-07-15 22:01 by 学员oN9JQD
[基金申请] H口，面上，时间戳7月7日。 +3 plumage5 2026-07-15 3/150 2026-07-15 21:06 by 6543yes
[基金申请] 小木虫没落了，除了祈祷帖子，几乎看不到有价值的帖子 +13 botar 2026-07-14 14/700 2026-07-15 19:28 by 寒流32
[基金申请] 会评 +5 无名者登山 2026-07-15 5/250 2026-07-15 18:15 by 杠就是你对
[基金申请] E口会评 +11 布布和一二 2026-07-13 13/650 2026-07-15 16:34 by Kittylucky
[基金申请] E口青C会评结束了吗？ 40+3 小小书虫c 2026-07-14 9/450 2026-07-15 16:04 by 安静的知了
[基金申请] 生命口C13面上会评时间 +6 luckinging 2026-07-14 9/450 2026-07-15 14:39 by 雨冰共舞
[教师之家] 内心匮乏 +13 水冰月月野兔 2026-07-12 13/650 2026-07-15 14:21 by jiduquegai
[文学芳草园] 你可曾有过这种感受？ +3 Jeanya 2026-07-09 6/300 2026-07-15 11:07 by mbtwt
[论文投稿] 农业机械学报投稿周期 +3 G0DLIN 2026-07-14 7/350 2026-07-15 09:26 by xs74101122
[基金申请] 国自然面上D口祈祷 +4 monkeynana 2026-07-14 4/200 2026-07-15 09:14 by 白水煮
[论文投稿] 跨出版社商投稿 20+4 倪好520 2026-07-13 4/200 2026-07-15 08:35 by bobvan
[基金申请] 没收到消息，再次陪跑 +3 天空星辰fly 2026-07-14 3/150 2026-07-14 19:14 by 6543yes
[基金申请] 不要再数国自然申请书的 filecode 的分隔符个数了 +13 sadpencil 2026-07-12 20/1000 2026-07-14 13:12 by 杠就是你对
[基金申请] 明天E口面上会评 +8 布布和一二 2026-07-12 9/450 2026-07-13 22:59 by QTSorange
[基金申请] 祈祷青基必中 50+6 hadengry 2026-07-09 16/800 2026-07-12 20:31 by 登山虫虫
[考博] 27届辽宁大学应届毕业生申博 +3 可乐微风大海 2026-07-09 3/150 2026-07-10 21:23 by 贪玩的研究僧
[基金申请] 关于如何从代码看上不上会 +17 且听虎啸 2026-07-09 23/1150 2026-07-10 19:51 by sdfapple719
信息提示
关闭
请填处理意见
关闭