| 查看: 989 | 回复: 2 | |||
fanghd2005铁虫 (小有名气)
|
[求助]
matlab有限差分法解Couette流出错 已有2人参与
|
|
各位前辈，我用有限 差分法解一维couette流时计算结果与教材中的结果不一样， 我用matlab写的程序，我对比了一下中间计算结果，与教材中的结果相同，矩阵与线性方程右端系数都相同，只有求解之后的结果不同。 请各位前辈帮我看看我的程序有什么问题，请指正。 %定义模型 Ylength=1; %平板长度 detaY=0.05; My=Ylength/detaY; detat=12.5; %时间步长 t=15000; %结晶计算时间 Nt=t/detat; %定义初始条件 u=zeros(My+1,1); %各点速度 u(My+1)=1; %计算过程 for AC=1:Nt %流动场,Crank-Nicolson算法 ReD=5000; A=-detat/(2*ReD*detaY^2); %差分系数 B=1+detat/(ReD*detaY^2); for j=1:My-1 %循环矩阵 P(j,j)=B; if j>1 P(j-1,j)=A;P(j,j-1)=A; end end for i=2:My %右端常数项 f(i-1,AC)=(1-detat/(detaY^2*ReD))*u(i,AC)+detat/(2*detaY^2*ReD)*(u(i+1,AC)+u(i-1,AC)); end u(2:My,AC+1)=P\f(:,AC); %此处出错 u(1,AC+1)=0; u(My+1,AC+1)=1; end |
回复此楼» 收录本帖的淘帖专辑推荐
 程序 |
» 猜你喜欢
070200求调剂，一志愿某211，288分
已经有20人回复
-
法国博士后职位
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有59人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有8人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有9人回复
-
桂林理工大学物理学专业招收调剂
已经有4人回复
FMStation
至尊木虫 (知名作家)
- 应助: 591 (博士)
- 贵宾: 0.03
- 金币: 18663.7
- 红花: 97
- 帖子: 8891
- 在线: 1491.2小时
- 虫号: 2400059
- 注册: 2013-04-04
- 专业: 计算机应用技术
2楼2016-08-26 17:28:11
512095518
专家顾问 (正式写手)
-
专家经验: +52 - 应助: 118 (高中生)
- 金币: 7665.4
- 红花: 78
- 帖子: 712
- 在线: 186.1小时
- 虫号: 3195576
- 注册: 2014-05-10
- 性别: GG
- 专业: 大气环境与全球气候变化
- 管辖: 计算模拟
3楼2016-08-27 10:15:48
