 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 985 | 回復(fù): 2 | |||
fanghd2005鐵蟲 (小有名氣)
|
[求助]
matlab有限差分法解Couette流出錯(cuò) 已有2人參與
|
|
各位前輩，我用有限 差分法解一維couette流時(shí)計(jì)算結(jié)果與教材中的結(jié)果不一樣， 我用matlab寫的程序，我對(duì)比了一下中間計(jì)算結(jié)果，與教材中的結(jié)果相同，矩陣與線性方程右端系數(shù)都相同，只有求解之后的結(jié)果不同。 請(qǐng)各位前輩幫我看看我的程序有什么問(wèn)題，請(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ì)算過(guò)程 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ù)此樓» 收錄本帖的淘帖專輯推薦
 程序 |
» 猜你喜歡
計(jì)算與實(shí)驗(yàn)科學(xué)工程國(guó)際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有4人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有92人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來(lái)能源學(xué)院
已經(jīng)有29人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有4人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來(lái)能源學(xué)院
已經(jīng)有26人回復(fù)
-
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
UJN物理學(xué)專業(yè)調(diào)劑
已經(jīng)有10人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有0人回復(fù)
-
北京納米能源與系統(tǒng)研究所王中林院士/曹南穎研究員課題組招收2026碩士研究生1名
已經(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
專家顧問(wèn) (正式寫手)
-
專家經(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ā)表 | |
|---|---|---|---|---|
|
[考研] 285求調(diào)劑 +5 | AZMK 2026-04-02 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 環(huán)境工程297分求調(diào)劑一志愿杭高院 +12 | GENJIOW 2026-03-31 | 13/650 |
|
|
|
[考研] 總分328生物與醫(yī)藥考數(shù)學(xué)求調(diào)劑 +3 | aaadim 2026-04-02 | 3/150 |
|
|
|
[考研] 初試301，代碼085701環(huán)境工程，本碩一致，四六級(jí)已過(guò)，有二區(qū)一作，共發(fā)表5篇論文 +6 | axibli 2026-04-01 | 6/300 |
|
|
|
[考研] 266分，一志愿電氣工程，本科材料，求材料專業(yè)調(diào)劑 +10 | 哇呼哼呼哼 2026-04-01 | 11/550 |
|
|
|
[考研] 271求調(diào)劑 +15 | 勒布朗@ 2026-03-31 | 20/1000 |
|
|
|
[考研] 07生物學(xué)求調(diào)劑 一志愿同濟(jì)大學(xué)359分 +3 | LAMC. 2026-03-30 | 3/150 |
|
|
|
[考研] 一志愿安徽大學(xué)計(jì)算機(jī)科學(xué)與技術(shù)學(xué)碩，331分求調(diào)劑 +5 | 蔣昌鵬qtj 2026-04-01 | 5/250 |
|
|
|
[考研] 求生物學(xué)調(diào)劑 +8 | 15172915737 2026-04-01 | 8/400 |
|
|
|
[考研] 279求調(diào)劑 +7 | 莫xiao 2026-04-01 | 7/350 |
|
|
|
[考研] 332求調(diào)劑 +8 | Lyy930824@ 2026-03-29 | 8/400 |
|
|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +7 | 阿姨復(fù)古不過(guò) 2026-03-27 | 7/350 |
|
|
|
[考研] 材料與化工調(diào)劑一志愿大連海事085600，349 +9 | 吃的不少 2026-03-30 | 9/450 |
|
|
|
[考研] 生物學(xué)296求調(diào)劑 +10 | 湯圓包 2026-03-29 | 14/700 |
|
|
|
[考研] 267求調(diào)劑 +13 | uiybh 2026-03-31 | 13/650 |
|
|
|
[考研] 求0861交通運(yùn)輸專碩or材料專碩調(diào)劑 +4 | 勒布朗@ 2026-03-31 | 4/200 |
|
|
|
[考研] 一志愿西交大080500材料學(xué)碩349 +6 | jqx1258 2026-03-31 | 7/350 |
|
|
|
[考研] 375求調(diào)劑 +7 | 雨夏整夜 2026-03-29 | 7/350 |
|
|
|
[考研] 材料工程專碩求調(diào)劑 +10 | hyl3153942 2026-03-29 | 10/500 |
|
|
|
[考研] 08工科求調(diào)劑286 +5 | tgs_001 2026-03-28 | 5/250 |
|