应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » LBM编程求解二维非稳态对流扩散方程
111/2返回列表12下一页
查看: 3236  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

dayovan

金虫 (小有名气)

[求助] LBM编程求解二维非稳态对流扩散方程 已有2人参与

用其他方法也可以，如有限元，有限体积，有限差分等！！！！！！！！谢谢

LBM编程求解二维非稳态对流扩散方程
二维区域对流扩散机制研究界面传质动力学无因次方程.jpg
回复此楼

» 猜你喜欢
Nothingisimpossible
1楼 2014-09-03 17:44:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

asiayazhou

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
dayovan: 金币+40, ★★★很有帮助, 能用LBM编程更好 2014-09-04 09:30:55
我用有限差分算好了，MATLAB程序见：链接: http://pan.baidu.com/s/1GnHVK 密码: yte6
请查收。谢谢！
赞一下(3人) 回复此楼
献身流体，献身航空发动机
4楼2014-09-03 21:03:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhp-jxu

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
dayovan: 金币+50, ★★★很有帮助 2015-03-29 10:29:41
引用回帖:
5楼: Originally posted by dayovan at 2014-09-04 15:08:15
有限差分计算浓度场出现负数，有办法解决不...

一阶导数应该用迎风离散格式，效果有明显改善
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—横坐标','fontsize',14)
ylabel('y—纵坐标','fontsize',14')
shading interp
colorbar
title('浓度场分布','fontsize',14)
赞一下(2人) 回复此楼
7楼2015-03-28 16:47:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

yishichao

银虫 (小有名气)
是控制方程吗？

[ 发自小木虫客户端 ]
回复此楼
2楼2014-09-03 19:13:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dayovan

金虫 (小有名气)
引用回帖:
2楼: Originally posted by yishichao at 2014-09-03 19:13:50
是控制方程吗？

就是上面无因次方程数值求解呀
赞一下 回复此楼
Nothingisimpossible
3楼2014-09-03 19:57:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dayovan

金虫 (小有名气)
引用回帖:
4楼: Originally posted by asiayazhou at 2014-09-03 21:03:39
我用有限差分算好了，MATLAB程序见：链接: http://pan.baidu.com/s/1GnHVK 密码: yte6
请查收。谢谢！

有限差分计算浓度场出现负数，有办法解决不
赞一下 回复此楼
Nothingisimpossible
5楼2014-09-04 15:08:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yishichao

银虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
???????:
5?: Originally posted by dayovan at 2014-09-04 15:08:15
???????????????????????а??????...

是边界上??????？

[ ????С??????? ]
赞一下 回复此楼
6楼2014-09-05 08:13:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dayovan

金虫 (小有名气)
引用回帖:
7楼: Originally posted by zhp-jxu at 2015-03-28 16:47:25
一阶导数应该用迎风离散格式，效果有明显改善
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 ...

能帮忙设置一下二阶迎风格式吗？？当N取较大时，尤发散了
赞一下 回复此楼
Nothingisimpossible
8楼2015-04-10 22:32:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dayovan

金虫 (小有名气)
引用回帖:
7楼: Originally posted by zhp-jxu at 2015-03-28 16:47:25
一阶导数应该用迎风离散格式，效果有明显改善
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 ...

能解决一下N取较大时，求解发散的问题吗？
赞一下 回复此楼
Nothingisimpossible
9楼2015-04-10 22:34:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhp-jxu

木虫 (小有名气)
引用回帖:
9楼: Originally posted by dayovan at 2015-04-10 22:34:49
能解决一下N取较大时，求解发散的问题吗？...

目前程序中的N的取值为N=100+2,如果增加N的值，前面的所有出现100的地方也要改
赞一下 回复此楼
10楼2015-04-11 10:22:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dayovan 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 295材料求调剂，一志愿武汉理工085601专硕 +4 Charlieyq 2026-03-19 4/200 2026-03-20 14:26 by 无懈可击111
[考研] 一志愿中南化学337求调剂 +4 niko- 2026-03-19 5/250 2026-03-20 13:49 by 促天成
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
[考研] 材料080500调剂求收留 +6 一颗meteor 2026-03-13 6/300 2026-03-20 10:41 by EBSD
[考研] 304求调剂 +5 曼殊2266 2026-03-18 5/250 2026-03-20 09:00 by ZHANG0tao
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
[考研] 一志愿 西北大学 ，070300化学学硕，总分287，双非一本，求调剂。 +3 晨昏线与星海 2026-03-19 3/150 2026-03-19 13:36 by houyaoxu
[考研] 材料考研调剂 +3 xwt。 2026-03-19 3/150 2026-03-19 11:22 by w沐阳w
[考研] 0854可跨调剂，一作一项核心论文五项专利，省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 0854，计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 考研化学学硕调剂，一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考研] 一志愿，福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 283求调剂 +3 听风就是雨； 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
信息提示
关闭
请填处理意见
关闭