两个壁面之间夹一个圆柱形杆的温度分布

两个壁面的温度已知且均高于流体温度，中间夹一个圆柱形杆

MATLAB代码为

%均为国际单位
t1=500;
t2=400;
tf=273;%流体温度
lam=1;
h=17;
l=0.3;
r=0.0625;
m=sqrt(2*h/lam/r);
theta1=t1-tf;
theta2=t2-tf;
f=@(x) theta1/ch(m*x)-theta2/ch(m*(l-x));
x1=fzero(f,[0,l]);%两个肋片的绝热面交点
x_1=0:0.01:x1;
x_2=x1+0.01:0.01:l;
y1=theta1.*ch(m.*(x_1-x1))./ch(m*x1)+tf;
y2=theta2.*ch(m.*(x1-x_2))./ch(m*x1)+tf;
plot([x_1,x_2],[y1,y2])
xlabel("x")
ylabel("温度")
function y=ch(x)
y=(exp(x)+exp(-x))./2;
end

效果如图