最近碰到一個(gè)問(wèn)題,在求解哦偏微分方程的時(shí)候,采用了有限差分法。但是結(jié)果總是不如人意。能不能麻煩各位幫我按照你們的思路寫一下呢?不甚感激,當(dāng)然如果有人是有償代寫,可以跟我聯(lián)系,我們交流一下!
![遇到一個(gè)matlab的問(wèn)題]()
然后我自己寫了點(diǎn),
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; %這個(gè)是入口條件
%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; %這些表達(dá)式都是參考別的文獻(xiàn)得來(lái)
end
for j=2:n;
T(1,j)=1/(1+Ramda/(Kw*R*dx))*(Tw+Ramda/(Kw*R*dx)*(T(2,j))); %這個(gè)是緊挨加熱管的溫度表達(dá)式
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; %這個(gè)是除邊界外,也就是微分方程表達(dá)式
%T(i,n+1)=T(i,n-1);
end
end
end
f=T
希望高手能按照你們的思路幫我看看,不要直接看我的程序,我的很垃圾的。
謝謝各位了。
未命名.jpg |