版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

lkbiandou

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1230.5
帖子: 208
在线: 94.6小时
虫号: 945069
注册: 2010-01-19
专业: 化学反应工程

[交流] 有限差分法求解偏微分方程已有1人参与

最近碰到一个问题，在求解哦偏微分方程的时候，采用了有限差分法。但是结果总是不如人意。能不能麻烦各位帮我按照你们的思路写一下呢？不甚感激，当然如果有人是有偿代写，可以跟我联系，我们交流一下！

然后我自己写了点，
function crack_2D
clear all; clc
Ramda = 0.45;    % W/(m K)
%d33=2;d44=3;
%T(1:11,1:11,1:11)=0;
T(1:10,1:10)=483;
x0=T(1:10,1:10);
%X=fsolve(@temperature,x0);
%function f=temperature(T)
m=10; n=10;K=10;
Kw=1.2;Tw=513;R=0.08; Ramda = 0.45;    % W/(m K)
x0=0.3959; dx=0.06041;df=pi/6/10;dz=0.1;
P=5;w=34.572; R0=0.0125;L=10;Dr=1;G=500;Cp=1;M=8.2636;rho=1440;Cor=0.9;
a5=L*M*Ramda/R0^2/G/Cp; a6=rho*Cor*L*M/G/Cp*0.32; %(-r1(i,j,k)*HCO(i,j,k)-r2(i,j,k)*HCO2(i,j,k));
%*************************************************************************%
for h=1:100
for i=1:m;
T(i,1)=483; %这个是入口条件
%T(i,n)=(4*T(i,n-1)-T(i,n-2))/3;
T(i,n+1)=(4*T(i,n)-T(i,n-1))/3; %这些表达式都是参考别的文献得来
end
for j=2:n;
T(1,j)=1/(1+Ramda/(Kw*R*dx))*(Tw+Ramda/(Kw*R*dx)*(T(2,j))); %这个是紧挨加热管的温度表达式
T(m+1,j)=(4*T(m,j)-T(m-1,j))/3;  %边界条件
end
for j=1:n;
for i=2:m;
      d44=(T(i+1,j)-2*T(i,j)+T(i-1,j)+T(i+1,j+1)-2*T(i,j+1)+T(i-1,j+1))/(2*dx^2)+1/(x0+i*dx)*(T(i+1,j)-T(i-1,j)+T(i+1,j+1)-T(i-1,j+1))/(4*dx);
      T(i,j+1)=T(i,j)+dz*a5*d44+dz*a6;  %这个是除边界外，也就是微分方程表达式
      %T(i,n+1)=T(i,n-1);
end
end
end
f=T

希望高手能按照你们的思路帮我看看，不要直接看我的程序，我的很垃圾的。
谢谢各位了。

未命名.jpg

[ Last edited by lkbiandou on 2013-9-27 at 16:34 ]

回复此楼

» 猜你喜欢

澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
深圳大学应用技术学院招聘凝聚态物理博士后已经有124人回复
物理学I论文润色/翻译怎么收费? 已经有257人回复
求助中医院临床中医中标的青年标书学习，谢谢已经有5人回复
浙大物院表面物理课题组长期招收普博及直博学生已经有259人回复
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry 已经有47人回复
25届硕士申请，材料计算方向已经有2人回复
课题组招1-2名硕士已经有1人回复
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上已经有0人回复
斯坦福大学与Elsevier世界前2%学者新排名查询已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2013-09-27 10:32:10

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

SEEDEN

银虫 (正式写手)

应助: 35 (小学生)
金币: 997.9
散金: 200
红花: 8
帖子: 528
在线: 169.1小时
虫号: 2157719
注册: 2012-11-30
专业: 传递过程

★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
dingd: 金币+1, 谢谢交流。 2013-09-28 21:56:23

你的Z是时间的无量纲话、X是长度的无量纲话吧
初始条件好像不全，(0，X0)范围内的初始条件是什么呢？
这个是最常见的抛物偏微分方程了
跟热传导方程的求解差不多，你找一下传热学的课本，看看热传导数值求解那一章，就明白了，真心不难

赞一下(2人)