| 查看: 3957 | 回复: 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 |
回复此楼» 猜你喜欢
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
梦想
已经有8人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教用matlab求解一个非线性偏微分方程组的数值解
已经有8人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】用MATLAB求解偏微分方程的数值解 悬赏20
已经有8人回复
- 【求助】求助用matlab解偏微分方程
已经有15人回复
shizijun
木虫 (著名写手)
自定义头衔
- 应助: 20 (小学生)
- 金币: 3662.6
- 散金: 2205
- 红花: 6
- 帖子: 2323
- 在线: 422.6小时
- 虫号: 443442
- 注册: 2007-10-27
- 性别: GG
- 专业: 有机高分子功能材料
4楼2012-09-17 11:36:25
shizijun
木虫 (著名写手)
自定义头衔
- 应助: 20 (小学生)
- 金币: 3662.6
- 散金: 2205
- 红花: 6
- 帖子: 2323
- 在线: 422.6小时
- 虫号: 443442
- 注册: 2007-10-27
- 性别: GG
- 专业: 有机高分子功能材料
2楼2012-09-17 10:27:05
crystal_ling
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1456.1
- 红花: 2
- 帖子: 80
- 在线: 51.3小时
- 虫号: 1977534
- 注册: 2012-09-05
- 专业: 光学
3楼2012-09-17 10:39:23
crystal_ling
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1456.1
- 红花: 2
- 帖子: 80
- 在线: 51.3小时
- 虫号: 1977534
- 注册: 2012-09-05
- 专业: 光学
5楼2012-09-19 17:22:01
30