| 查看: 983 | 回复: 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 |
回复此楼» 收录本帖的淘帖专辑推荐
 程序 |
» 猜你喜欢
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有294人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
FMStation
至尊木虫 (知名作家)
- 应助: 591 (博士)
- 贵宾: 0.03
- 金币: 18366.2
- 红花: 97
- 帖子: 8891
- 在线: 1490.6小时
- 虫号: 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
