版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

dayovan

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 3069.2
散金: 380
帖子: 228
在线: 125.2小时
虫号: 2321676
注册: 2013-03-05
性别: GG
专业: 传递过程

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

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

二维区域对流扩散机制研究界面传质动力学无因次方程.jpg

回复此楼

» 猜你喜欢

南开大学物理学院张书辉副教授招收凝聚态物理理论方向博士生、硕士生已经有22人回复
杂志要求提供official language-editing certificate 已经有4人回复
物理学I论文润色/翻译怎么收费? 已经有50人回复
深圳大学2026年秋博士招生-物理学-活性胶体方向-高永祥课题组已经有17人回复
论物质与能量的统一模型及物理现象解释已经有0人回复
基于基元I统一理论的数学相关应用推导已经有0人回复
基元I统一理论：宇宙本质、层级演化与修炼文明的本源规律已经有1人回复
基元I理论下三大核心空间现象精准推导与细节解析已经有0人回复
基于基元 I 统一理论的反重力理论推导已经有0人回复
基于基元I统一理论的量子力学本源推导已经有1人回复
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费已经有11人回复

Nothingisimpossible

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

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

asiayazhou

金虫 (正式写手)

应助: 3 (幼儿园)
金币: 994.1
帖子: 890
在线: 87.4小时
虫号: 2561592
注册: 2013-07-23
性别: GG
专业: 水文地质学(含地热地质学)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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

木虫 (小有名气)

应助: 7 (幼儿园)
金币: 2292.3
散金: 18
红花: 1
帖子: 73
在线: 438.3小时
虫号: 1587212
注册: 2012-01-26
性别: GG
专业: 计算数学与科学工程计算

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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 (幼儿园)
金币: 228.4
散金: 4
帖子: 53
在线: 43.9小时
虫号: 2745040
注册: 2013-10-22
性别: GG
专业: 计算数学与科学工程计算

是控制方程吗？

[ 发自小木虫客户端 ]

回复此楼

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dayovan

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 3069.2
散金: 380
帖子: 228
在线: 125.2小时
虫号: 2321676
注册: 2013-03-05
性别: GG
专业: 传递过程

引用回帖:

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

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 3069.2
散金: 380
帖子: 228
在线: 125.2小时
虫号: 2321676
注册: 2013-03-05
性别: GG
专业: 传递过程

引用回帖:

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

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 228.4
散金: 4
帖子: 53
在线: 43.9小时
虫号: 2745040
注册: 2013-10-22
性别: GG
专业: 计算数学与科学工程计算

【答案】应助回帖

感谢参与，应助指数 +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

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 3069.2
散金: 380
帖子: 228
在线: 125.2小时
虫号: 2321676
注册: 2013-03-05
性别: GG
专业: 传递过程

引用回帖:

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取较大时，尤发散了