[求助]求助:如何用MATLAB的偏微分方程工具箱畫出不同區(qū)域內(nèi)的溫度分布圖?
初始懸賞金幣 5 個
求助:MATLAB的偏微分方程工具箱問題
問題:畫出下圖所示模型內(nèi)的溫度分布圖。
穩(wěn)定狀態(tài)下的溫度分布通過數(shù)值求解標準的熱傳導方程得到
-▽﹒(k▽T)=Q
其中各層的熱載Q及參數(shù)如下:
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 ] |