用Matlab模拟温度分布

Question

我目前正在使用内部热水盘管加热的储罐上。我需要知道储罐中的温度分布如何。

我曾在Matlab中尝试过稳态模型，但我发现它不是很精确（我只能输入产品和储罐材料的热导率+边界条件），在我看来，这使分布过于笼统而不是针对我的情况。我也尝试了瞬态模型，但是在运行代码时出现了一些错误。

这是我用于瞬态模型的代码（我是Matlab的新手，所以请原谅任何代码错误）

model = createpde('thermal','transient');
geo = multicylinder([3.4,3.5],7.1)
model.Geometry= geo;
pdegplot(model,'FaceLabels',"on")
generateMesh(model);
pdemesh(model)


 thermalProperties(model,"Cell",2,"ThermalConductivity",0.1,"MassDensity",576,"SpecificHeat",0.5)
 thermalProperties(model,"Cell",1,"MassDensity",932,"SpecificHeat",2.28,"ThermalConductivity",0.12)
 thermalIC(model,0)
 tlist=0:.1:5

Thermalresults = solve(model,tlist);
u= Thermalresults.Temperature
pdeplot3D(model,'ColorMapData', u)