| 查看: 3281 | 回复: 10 | |||
dayovan金虫 (小有名气)
|
[求助]
LBM编程求解二维非稳态对流扩散方程 已有2人参与
|
||
用其他方法也可以，如有限元，有限体积，有限差分等！！！！！！！！谢谢 二维区域对流扩散机制研究界面传质动力学无因次方程.jpg |
回复此楼» 猜你喜欢
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有199人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有11人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有12人回复
-
桂林理工大学物理学专业招收调剂
已经有12人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
美国顶级物理期刊《应用物理快报》（APL）的编辑欺骗和歧视作者及AIP的官僚主义傲慢
已经有1人回复
asiayazhou
金虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 994.1
- 帖子: 890
- 在线: 87.4小时
- 虫号: 2561592
- 注册: 2013-07-23
- 性别: GG
- 专业: 水文地质学(含地热地质学)
4楼2014-09-03 21:03:39
zhp-jxu
木虫 (小有名气)
- 应助: 7 (幼儿园)
- 金币: 2292.3
- 散金: 18
- 红花: 1
- 帖子: 73
- 在线: 438.3小时
- 虫号: 1587212
- 注册: 2012-01-26
- 性别: GG
- 专业: 计算数学与科学工程计算
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
dayovan: 金币+50, ★★★很有帮助 2015-03-29 10:29:41
dayovan: 金币+50, ★★★很有帮助 2015-03-29 10:29:41
|
一阶导数应该用迎风离散格式，效果有明显改善 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) |
7楼2015-03-28 16:47:25
yishichao
银虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 228.4
- 散金: 4
- 帖子: 53
- 在线: 43.9小时
- 虫号: 2745040
- 注册: 2013-10-22
- 性别: GG
- 专业: 计算数学与科学工程计算
2楼2014-09-03 19:13:50
dayovan
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3069.2
- 散金: 380
- 帖子: 228
- 在线: 125.2小时
- 虫号: 2321676
- 注册: 2013-03-05
- 性别: GG
- 专业: 传递过程
3楼2014-09-03 19:57:46
dayovan
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3069.2
- 散金: 380
- 帖子: 228
- 在线: 125.2小时
- 虫号: 2321676
- 注册: 2013-03-05
- 性别: GG
- 专业: 传递过程
5楼2014-09-04 15:08:15
yishichao
银虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 228.4
- 散金: 4
- 帖子: 53
- 在线: 43.9小时
- 虫号: 2745040
- 注册: 2013-10-22
- 性别: GG
- 专业: 计算数学与科学工程计算
6楼2014-09-05 08:13:40
dayovan
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3069.2
- 散金: 380
- 帖子: 228
- 在线: 125.2小时
- 虫号: 2321676
- 注册: 2013-03-05
- 性别: GG
- 专业: 传递过程
8楼2015-04-10 22:32:43
dayovan
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3069.2
- 散金: 380
- 帖子: 228
- 在线: 125.2小时
- 虫号: 2321676
- 注册: 2013-03-05
- 性别: GG
- 专业: 传递过程
9楼2015-04-10 22:34:49
zhp-jxu
木虫 (小有名气)
- 应助: 7 (幼儿园)
- 金币: 2292.3
- 散金: 18
- 红花: 1
- 帖子: 73
- 在线: 438.3小时
- 虫号: 1587212
- 注册: 2012-01-26
- 性别: GG
- 专业: 计算数学与科学工程计算
10楼2015-04-11 10:22:12
