偏微分方程組解法-資料下載頁
2025-06-19 21:46本頁面
【正文】 外壁溫分別維持在為80、25，試確定該通道溫度分布隨時間的變化及溫度分布穩(wěn)定所需的時間。 解：pdetool工具各菜單中的參數(shù)設(shè)置略，其中solve菜單parameters選項的時間設(shè)置為0:100:8000。 首先在穩(wěn)態(tài)情況下（橢圓形）求解（三角形網(wǎng)格經(jīng)過兩次細化），導(dǎo)出解u；而后解以上非穩(wěn)態(tài)（拋物形）并導(dǎo)出解u1（三角形網(wǎng)格細化與穩(wěn)態(tài)相同）和三角形網(wǎng)格參數(shù)p，e，t；命令：pdemesh(p,e,t,u1(:,1)) 0spdemesh(p,e,t,u1(:,2)) 100spdemesh(p,e,t,u1(:,4)) 300spdemesh(p,e,t,u1(:,7)) 600spdemesh(p,e,t,u1(:,11)) 1000spdemesh(p,e,t,u1(:,21)) 2000spdemesh(p,e,t,u1(:,31)) 3000spdemesh(p,e,t,u1(:,51)) 5000spdemesh(p,e,t,u1(:,81)) 8000s穩(wěn)定時間（以某時刻的溫度分布與穩(wěn)態(tài)溫度分布的相對誤差來判斷）norm((uu1(:,21))./u) 結(jié)果：norm((uu1(:,31))./u) norm((uu1(:,51))./u) norm((uu1(:,81))./u) 可見在8000s以后的溫度分布已基本沒有變化。h=1 r=0 c=1 a=0 f=1誤差plot Height3D中選擇user entry : u(1x.^2y.^2)/4xyh=1 r=0 c=1 a=0 f=1xy需采用solve parameters中use nonlinear solver方可收斂
畢業(yè)設(shè)計相關(guān)推薦
鄂ICP備17016276號-1