24小时热门版块排行榜

返回列表

fanghd2005

铁虫 (小有名气)

应助: 1 (幼儿园)
金币: 722.3
散金: 20
红花: 1
帖子: 60
在线: 33小时
虫号: 3065504
注册: 2014-03-18
性别: GG
专业: 塑料

[求助] 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年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
物理学I论文润色/翻译怎么收费? 已经有202人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有23人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复

不会的太多

1楼 2016-08-26 16:30:36

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

FMStation

至尊木虫 (知名作家)

应助: 591 (博士)
贵宾: 0.03
金币: 18115.7
红花: 97
帖子: 8891
在线: 1485.3小时
虫号: 2400059
注册: 2013-04-04
专业: 计算机应用技术

【答案】应助回帖

感谢参与，应助指数 +1

u'

ans =

  Columns 1 through 13

   0    0    0    0    0    0    0    0    0    0    0    0    0

  Columns 14 through 21

   0    0    0    0    0    0    0    1