24小时热门版块排行榜

>论坛更新日志 (3617)
>虫友互识 (625)
>导师招生 (233)
>文献求助 (214)
>硕博家园 (119)
>考博 (101)
>公派出国 (98)
>论文投稿 (79)
>博后之家 (76)
>找工作 (73)
>基金申请 (72)
>休闲灌水 (71)
>考研 (64)
>论文道贺祈福 (62)
>教师之家 (54)
>招聘信息布告栏 (44)

返回列表

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

回复此楼

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

程序

» 猜你喜欢

求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有283人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复
求助NH4V4O10晶体的CIF文件已经有0人回复
英国全奖博士招聘-深度学习与量子物理已经有0人回复
间接带隙半导体有效质量求助已经有0人回复

不会的太多

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

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

FMStation

至尊木虫 (知名作家)

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