| 查看: 1075 | 回复: 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 |
回复此楼» 收录本帖的淘帖专辑推荐
 程序 |
» 猜你喜欢
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有195人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
FMStation
至尊木虫 (知名作家)
- 应助: 591 (博士)
- 贵宾: 0.03
- 金币: 18710.7
- 红花: 97
- 帖子: 8891
- 在线: 1495小时
- 虫号: 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
30
