| 查看: 1014 | 回复: 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 |
回复此楼» 收录本帖的淘帖专辑推荐
 程序 |
» 猜你喜欢
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有210人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
-
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有1人回复
-
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生
已经有3人回复
-
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究
已经有12人回复
FMStation
至尊木虫 (知名作家)
- 应助: 591 (博士)
- 贵宾: 0.03
- 金币: 18701.7
- 红花: 97
- 帖子: 8891
- 在线: 1493.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
5