 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 986 | 回復(fù): 2 | |||
fanghd2005鐵蟲 (小有名氣)
|
[求助]
matlab有限差分法解Couette流出錯(cuò) 已有2人參與
|
|
各位前輩，我用有限 差分法解一維couette流時(shí)計(jì)算結(jié)果與教材中的結(jié)果不一樣， 我用matlab寫的程序，我對(duì)比了一下中間計(jì)算結(jié)果，與教材中的結(jié)果相同，矩陣與線性方程右端系數(shù)都相同，只有求解之后的結(jié)果不同。 請(qǐng)各位前輩幫我看看我的程序有什么問題，請(qǐng)指正。 %定義模型 Ylength=1; %平板長(zhǎng)度 detaY=0.05; My=Ylength/detaY; detat=12.5; %時(shí)間步長(zhǎng) t=15000; %結(jié)晶計(jì)算時(shí)間 Nt=t/detat; %定義初始條件 u=zeros(My+1,1); %各點(diǎn)速度 u(My+1)=1; %計(jì)算過程 for AC=1:Nt %流動(dòng)場(chǎng),Crank-Nicolson算法 ReD=5000; A=-detat/(2*ReD*detaY^2); %差分系數(shù) B=1+detat/(ReD*detaY^2); for j=1:My-1 %循環(huán)矩陣 P(j,j)=B; if j>1 P(j-1,j)=A;P(j,j-1)=A; end end for i=2:My %右端常數(shù)項(xiàng) f(i-1,AC)=(1-detat/(detaY^2*ReD))*u(i,AC)+detat/(2*detaY^2*ReD)*(u(i+1,AC)+u(i-1,AC)); end u(2:My,AC+1)=P\f(:,AC); %此處出錯(cuò) u(1,AC+1)=0; u(My+1,AC+1)=1; end |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 程序 |
» 猜你喜歡
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來能源學(xué)院
已經(jīng)有27人回復(fù)
-
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有52人回復(fù)
-
UJN物理學(xué)專業(yè)調(diào)劑
已經(jīng)有12人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有0人回復(fù)
-
北京納米能源與系統(tǒng)研究所王中林院士/曹南穎研究員課題組招收2026碩士研究生1名
已經(jīng)有0人回復(fù)
-
山西大同大學(xué)物理學(xué)碩士研究生招收調(diào)劑生
已經(jīng)有0人回復(fù)
-
[調(diào)劑信息]211智能人工感知方向國(guó)家青年特聘專家課題組招收調(diào)劑研究生
已經(jīng)有0人回復(fù)
-
［調(diào)劑信息］211智能人工感知方向國(guó)家青年特聘專家課題組招收調(diào)劑研究生
已經(jīng)有0人回復(fù)
-
[調(diào)劑信息]211智能人工感知方向國(guó)家青年特聘專家課題組招收調(diào)劑研究生
已經(jīng)有0人回復(fù)
FMStation
至尊木蟲 (知名作家)
- 應(yīng)助: 591 (博士)
- 貴賓: 0.03
- 金幣: 18384.7
- 紅花: 97
- 帖子: 8891
- 在線: 1490.6小時(shí)
- 蟲號(hào): 2400059
- 注冊(cè): 2013-04-04
- 專業(yè): 計(jì)算機(jī)應(yīng)用技術(shù)
2樓2016-08-26 17:28:11
512095518
專家顧問 (正式寫手)
-
專家經(jīng)驗(yàn): +52 - 應(yīng)助: 118 (高中生)
- 金幣: 7665.4
- 紅花: 78
- 帖子: 712
- 在線: 186.1小時(shí)
- 蟲號(hào): 3195576
- 注冊(cè): 2014-05-10
- 性別: GG
- 專業(yè): 大氣環(huán)境與全球氣候變化
- 管轄: 計(jì)算模擬
3樓2016-08-27 10:15:48
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 fanghd2005 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-04-01 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 321求調(diào)劑 一志愿 浙江工業(yè)大學(xué)生物醫(yī)藥 +5 | 嘿嘿HC 2026-04-01 | 6/300 |
|
|
|
[考研] 298求調(diào)劑 +3 | zzz，，r 2026-04-02 | 6/300 |
|
|
|
[考研] 學(xué)碩化學(xué)工程與技術(shù)，一志愿中國(guó)海洋大學(xué)320+求調(diào)劑 +8 | 披星河 2026-04-02 | 8/400 |
|
|
|
[考研] 301求調(diào)劑 +14 | 駱駝男人 2026-04-02 | 14/700 |
|
|
|
[考研] 085601一志愿中山大學(xué)深圳材料工程330求調(diào)劑 +8 | pipiver 2026-03-30 | 8/400 |
|
|
|
[考研] 一志愿北交大材料工程總分358 +3 | cs0106 2026-04-02 | 5/250 |
|
|
|
[考研] 計(jì)算機(jī)265跨調(diào)環(huán)境 +5 | Yumeng_6 2026-03-27 | 5/250 |
|
|
|
[考研] 337求調(diào)劑 +11 | 《樹》 2026-03-29 | 11/550 |
|
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +12 | hanamiko 2026-03-27 | 12/600 |
|
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +7 | 哇呼哼呼哼 2026-03-26 | 7/350 |
|
|
|
[考研] 286求調(diào)劑 +5 | lim0922 2026-03-26 | 5/250 |
|
|
|
[考研] 材料專業(yè)調(diào)劑 +5 | 啦啦啦哭 2026-03-31 | 6/300 |
|
|
|
[考研] 322求調(diào)劑 +8 | 三水sss 2026-04-01 | 8/400 |
|
|
|
[考研] 289求調(diào)劑 +7 | BrightLL 2026-03-29 | 7/350 |
|
|
|
[考研] 一志愿西交大080500材料學(xué)碩349 +6 | jqx1258 2026-03-31 | 7/350 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-31 | 4/200 |
|
|
|
[考研] 322求調(diào)劑 +10 | 宋明欣 2026-03-27 | 10/500 |
|
|
|
[考研] 085600，專業(yè)課化工原理，321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
|
[考研] 一志愿吉大071010，316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|