dayovan

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 3069.2
散金: 380
帖子: 228
在線: 125.2小時(shí)
蟲號(hào): 2321676
注冊(cè): 2013-03-05
性別: GG
專業(yè): 傳遞過(guò)程

[求助] LBM編程求解二維非穩(wěn)態(tài)對(duì)流擴(kuò)散方程已有2人參與

用其他方法也可以，如有限元，有限體積，有限差分等�。。。。。。。≈x謝

二維區(qū)域?qū)α鲾U(kuò)散機(jī)制研究界面?zhèn)髻|(zhì)動(dòng)力學(xué)無(wú)因次方程.jpg

回復(fù)此樓

» 猜你喜歡

2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有194人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有5人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問(wèn)還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有30人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有60人回復(fù)

Nothingisimpossible

1樓 2014-09-03 17:44:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

asiayazhou

金蟲 (正式寫手)

應(yīng)助: 3 (幼兒園)
金幣: 994.1
帖子: 890
在線: 87.4小時(shí)
蟲號(hào): 2561592
注冊(cè): 2013-07-23
性別: GG
專業(yè): 水文地質(zhì)學(xué)(含地?zé)岬刭|(zhì)學(xué))

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
dayovan: 金幣+40, ★★★很有幫助, 能用LBM編程更好 2014-09-04 09:30:55

我用有限差分算好了，MATLAB程序見：鏈接: http://pan.baidu.com/s/1GnHVK 密碼: yte6
請(qǐng)查收。謝謝！

贊一下(3人)

回復(fù)此樓

獻(xiàn)身流體，獻(xiàn)身航空發(fā)動(dòng)機(jī)

4樓2014-09-03 21:03:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhp-jxu

木蟲 (小有名氣)

應(yīng)助: 7 (幼兒園)
金幣: 2292.3
散金: 18
紅花: 1
帖子: 73
在線: 438.3小時(shí)
蟲號(hào): 1587212
注冊(cè): 2012-01-26
性別: GG
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
dayovan: 金幣+50, ★★★很有幫助 2015-03-29 10:29:41

引用回帖:

5樓: Originally posted by dayovan at 2014-09-04 15:08:15
有限差分計(jì)算濃度場(chǎng)出現(xiàn)負(fù)數(shù)，有辦法解決不...

一階導(dǎo)數(shù)應(yīng)該用迎風(fēng)離散格式，效果有明顯改善
clc
clear
delta_x=pi/100;
delta_y=pi/100;
delta_t=0.001;%定義時(shí)間步長(zhǎng)；
x=0:delta_x:pi;
y=0:delta_y:pi;
t=0;
[xx,yy]=meshgrid(x,y);
N=102;
c_new=zeros(N,N);
c_old=zeros(N,N);
for i = 1:N-1
for j = 1:N-1
      p(i,j)=0.5*pi^2*sin(x(j))*cos(y(i));
      q(i,j)=-0.5*pi^2*cos(x(j))*sin(y(i));
end
end
for j=1:N
   c_old(1,

=1;
end
while t<=1
for i=2:N-1
for j=2:N-1
      if p(i-1,j-1)>=0
         div_x=(c_old(i,j)-c_old(i,j-1))/delta_x;
      else
         div_x=(c_old(i,j+1)-c_old(i,j))/delta_x;
      end
      if q(i-1,j-1)>=0
         div_y=(c_old(i,j)-c_old(i-1,j))/delta_y;
      else
         div_y=(c_old(i+1,j)-c_old(i,j))/delta_y;
      end
      diffu=(c_old(i,j+1)+c_old(i,j-1)-2*c_old(i,j))/power(delta_x,2)+...
         (c_old(i+1,j)+c_old(i-1,j)-2*c_old(i,j))/power(delta_y,2);
      temp=0.006*diffu-q(i-1,j-1)*div_y-...
         p(i-1,j-1)*div_x;
      c_new(i,j)= c_old(i,j)+delta_t*temp;
end
end
for j=1:N
c_new(1,j)=1;
c_new(N,j)=0;
end
for i=2:N-1
c_new(i,1)=1/3*(4*c_new(i,2)-c_new(i,3));
c_new(i,N)=1/3*(4*c_new(i,N-1)-c_new(i,N-2));
end
  for i=1:N
   for j=1:N
   c_old(i,j)=c_new(i,j);
   end
  end
t=t+delta_t;

for i=2:N-1
   for j=2:N-1
   c(i,j)=1/4*(c_new(i+1,j)+c_new(i-1,j)+c_new(i,j-1)+c_new(i,j+1));
   end
end
c(1,

=1;
c(2:N-1,1)=1/3*(4*c(2:N-1,2)-c(2:N-1,3));
end

figure('color','white');
surf(xx,yy,c);
view(2)
xlabel('x—橫坐標(biāo)','fontsize',14)
ylabel('y—縱坐標(biāo)','fontsize',14')
shading interp
colorbar
title('濃度場(chǎng)分布','fontsize',14)

贊一下(2人)

回復(fù)此樓

7樓2015-03-28 16:47:25

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

yishichao

銀蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 228.4
散金: 4
帖子: 53
在線: 43.9小時(shí)
蟲號(hào): 2745040
注冊(cè): 2013-10-22
性別: GG
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

是控制方程嗎？

[ 發(fā)自小木蟲客戶端 ]

回復(fù)此樓

2樓2014-09-03 19:13:50

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dayovan

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 3069.2
散金: 380
帖子: 228
在線: 125.2小時(shí)
蟲號(hào): 2321676
注冊(cè): 2013-03-05
性別: GG
專業(yè): 傳遞過(guò)程

引用回帖:

2樓: Originally posted by yishichao at 2014-09-03 19:13:50
是控制方程嗎？

就是上面無(wú)因次方程數(shù)值求解呀

贊一下

回復(fù)此樓

Nothingisimpossible

3樓2014-09-03 19:57:46

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dayovan

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 3069.2
散金: 380
帖子: 228
在線: 125.2小時(shí)
蟲號(hào): 2321676
注冊(cè): 2013-03-05
性別: GG
專業(yè): 傳遞過(guò)程

引用回帖:

4樓: Originally posted by asiayazhou at 2014-09-03 21:03:39
我用有限差分算好了，MATLAB程序見：鏈接: http://pan.baidu.com/s/1GnHVK 密碼: yte6
請(qǐng)查收。謝謝！

有限差分計(jì)算濃度場(chǎng)出現(xiàn)負(fù)數(shù)，有辦法解決不

贊一下

回復(fù)此樓

Nothingisimpossible

5樓2014-09-04 15:08:15

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yishichao

銀蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 228.4
散金: 4
帖子: 53
在線: 43.9小時(shí)
蟲號(hào): 2745040
注冊(cè): 2013-10-22
性別: GG
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

???????:

5?: Originally posted by dayovan at 2014-09-04 15:08:15
???????????????????????а??????...

是邊界上??????？

[ ????С??????? ]

贊一下

回復(fù)此樓

6樓2014-09-05 08:13:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dayovan

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 3069.2
散金: 380
帖子: 228
在線: 125.2小時(shí)
蟲號(hào): 2321676
注冊(cè): 2013-03-05
性別: GG
專業(yè): 傳遞過(guò)程

引用回帖:

7樓: Originally posted by zhp-jxu at 2015-03-28 16:47:25
一階導(dǎo)數(shù)應(yīng)該用迎風(fēng)離散格式，效果有明顯改善
clc
clear
delta_x=pi/100;
delta_y=pi/100;
delta_t=0.001;%定義時(shí)間步長(zhǎng)；
x=0:delta_x:pi;
y=0:delta_y:pi;
t=0;
=meshgrid(x,y);
N=102;
c_new=zeros ...

能幫忙設(shè)置一下二階迎風(fēng)格式嗎？？當(dāng)N取較大時(shí)，尤發(fā)散了

贊一下

回復(fù)此樓

Nothingisimpossible

8樓2015-04-10 22:32:43

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dayovan

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 3069.2
散金: 380
帖子: 228
在線: 125.2小時(shí)
蟲號(hào): 2321676
注冊(cè): 2013-03-05
性別: GG
專業(yè): 傳遞過(guò)程

引用回帖:

能解決一下N取較大時(shí)，求解發(fā)散的問(wèn)題嗎？

贊一下

回復(fù)此樓

Nothingisimpossible

9樓2015-04-10 22:34:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhp-jxu

木蟲 (小有名氣)

應(yīng)助: 7 (幼兒園)
金幣: 2292.3
散金: 18
紅花: 1
帖子: 73
在線: 438.3小時(shí)
蟲號(hào): 1587212
注冊(cè): 2012-01-26
性別: GG
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

引用回帖:

9樓: Originally posted by dayovan at 2015-04-10 22:34:49
能解決一下N取較大時(shí)，求解發(fā)散的問(wèn)題嗎？

...

目前程序中的N的取值為N=100+2,如果增加N的值，前面的所有出現(xiàn)100的地方也要改

贊一下

回復(fù)此樓

10樓2015-04-11 10:22:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 求調(diào)劑 +3	Mqqqqqq 2026-03-19	3/150	2026-03-19 14:11 by peike
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠(chéng)心求調(diào)劑 +4	吃吃吃才有意義 2026-03-19	4/200	2026-03-19 13:48 by houyaoxu
[考研] 281求調(diào)劑（0805） +9	煙汐憶海 2026-03-16	19/950	2026-03-19 11:42 by laoshidan
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +4	生物工程調(diào)劑 2026-03-16	12/600	2026-03-19 11:10 by 生物工程調(diào)劑
[考研] 304求調(diào)劑 +6	司空. 2026-03-18	6/300	2026-03-18 23:03 by 星空星月
[考研] 295求調(diào)劑 +3	一志愿京區(qū)211 2026-03-18	5/250	2026-03-18 17:03 by zhaoqian0518
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +7	手機(jī)用戶 2026-03-17	7/350	2026-03-18 14:34 by vgtyfty
[考研] 0703化學(xué)求調(diào)劑總分331 +3	ZY-05 2026-03-13	3/150	2026-03-18 10:58 by macy2011
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 275求調(diào)劑 +4	太陽(yáng)花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 285求調(diào)劑 +6	ytter 2026-03-12	6/300	2026-03-16 15:05 by njzyff
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +6	zlingli 2026-03-13	6/300	2026-03-15 20:00 by ryzcf
[考研] 070305求調(diào)劑 +3	mlpqaz03 2026-03-14	4/200	2026-03-15 11:04 by peike
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流 +4	自由的_飛翔 2026-03-13	5/250	2026-03-14 19:39 by hmn_wj
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 333求調(diào)劑 +3	152697 2026-03-12	4/200	2026-03-13 07:08 by Iveryant
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3	崔wj 2026-03-12	4/200	2026-03-12 19:33 by 求調(diào)劑zz