dayovan

金蟲 (小有名氣)

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

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

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

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

回復(fù)此樓

» 猜你喜歡

2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經(jīng)有0人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有6人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有80人回復(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ù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有30人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有62人回復(fù)
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣已經(jīng)有1人回復(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小時
蟲號: 2561592
注冊: 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
請查收。謝謝！

贊一下(3人)

回復(fù)此樓

獻(xiàn)身流體，獻(xiàn)身航空發(fā)動機(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小時
蟲號: 1587212
注冊: 2012-01-26
性別: GG
專業(yè): 計算數(shù)學(xué)與科學(xué)工程計算

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

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

引用回帖:

5樓: Originally posted by dayovan at 2014-09-04 15:08:15
有限差分計算濃度場出現(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;%定義時間步長；
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('濃度場分布','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小時
蟲號: 2745040
注冊: 2013-10-22
性別: GG
專業(yè): 計算數(shù)學(xué)與科學(xué)工程計算

是控制方程嗎？

[ 發(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小時
蟲號: 2321676
注冊: 2013-03-05
性別: GG
專業(yè): 傳遞過程

引用回帖:

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

就是上面無因次方程數(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小時
蟲號: 2321676
注冊: 2013-03-05
性別: GG
專業(yè): 傳遞過程

引用回帖:

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

有限差分計算濃度場出現(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小時
蟲號: 2745040
注冊: 2013-10-22
性別: GG
專業(yè): 計算數(shù)學(xué)與科學(xué)工程計算

【答案】應(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小時
蟲號: 2321676
注冊: 2013-03-05
性別: GG
專業(yè): 傳遞過程

引用回帖:

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;%定義時間步長；
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取較大時，尤發(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小時
蟲號: 2321676
注冊: 2013-03-05
性別: GG
專業(yè): 傳遞過程

引用回帖:

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

贊一下

回復(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小時
蟲號: 1587212
注冊: 2012-01-26
性別: GG
專業(yè): 計算數(shù)學(xué)與科學(xué)工程計算

引用回帖:

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

...

目前程序中的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ā)表

[考研] 265求調(diào)劑 +4	梁梁校校 2026-03-19	4/200	2026-03-20 00:03 by 23Postgrad
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考！ +3	lemonzzn 2026-03-17	3/150	2026-03-19 23:58 by 23Postgrad
[考研] 288求調(diào)劑 +15	于海海海海 2026-03-19	15/750	2026-03-19 22:41 by 學(xué)員8dgXkO
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +4	生物工程調(diào)劑 2026-03-16	12/600	2026-03-19 11:10 by 生物工程調(diào)劑
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx
[考研] 材料專業(yè)求調(diào)劑 +5	hanamiko 2026-03-18	5/250	2026-03-18 20:19 by 星空星月
[考研] 085600材料與化工 +5	安全上岸！ 2026-03-16	5/250	2026-03-18 15:33 by cmz0325
[考研] 311求調(diào)劑 +6	26研0 2026-03-15	6/300	2026-03-18 14:43 by haxia
[考研] 一志愿西南交大，求調(diào)劑 +4	材化逐夢人 2026-03-18	4/200	2026-03-18 14:22 by 007_lilei
[考研] 299求調(diào)劑 +5	△小透明* 2026-03-17	5/250	2026-03-18 11:49 by 盡舜堯1
[考研] 0703化學(xué)求調(diào)劑總分331 +3	ZY-05 2026-03-13	3/150	2026-03-18 10:58 by macy2011
[考研] 環(huán)境工程調(diào)劑 +8	大可digkids 2026-03-16	8/400	2026-03-18 09:36 by zhukairuo
[考研] 301求調(diào)劑 +9	yy要上岸呀 2026-03-17	9/450	2026-03-18 08:58 by 無際的草原
[考研] 278求調(diào)劑 +5	煙火先于春 2026-03-17	5/250	2026-03-18 08:43 by 星空星月
[考研] 293求調(diào)劑 +11	zjl的號 2026-03-16	16/800	2026-03-18 08:10 by zhukairuo
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 304求調(diào)劑 +4	ahbd 2026-03-14	4/200	2026-03-16 16:48 by 我的船我的海
[考研] 0703 物理化學(xué)調(diào)劑 +3	我可以上岸的對?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對?/a>