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

回复此楼

» 收录本帖的淘帖专辑推荐

程序

» 猜你喜欢

津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名已经有0人回复
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有280人回复
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）已经有6人回复
【修正版】单环用磁约束低速电子实验方案（简化版）已经有0人回复
桂林理工大学物理学专业招收调剂，还有三个名额！！！已经有22人回复
考博自荐已经有4人回复
山东大学第二批博士研究生招生已经有0人回复
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学已经有0人回复
《电磁学》教材推荐已经有1人回复
【急招】合肥工大核聚变材料计算方向2026级工程博士生已经有4人回复

不会的太多

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

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

FMStation

至尊木虫 (知名作家)

应助: 591 (博士)
贵宾: 0.03
金币: 18701.7
红花: 97
帖子: 8891
在线: 1493.2小时
虫号: 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