Why the solution is higher than the initial conditions in a model without heatsource?
1 view (last 30 days)
I want to simulate heat conduction and convection without internal heat source and look at the temperature distribution.The initial temperature is 50.But the result of the temperature is higher than the initial condition, why is that?
Here is my code.
Any help would be greatly appreciated.
model = createpde('thermal','transient');
gm = importGeometry(model,'1.STL');
[gm,ID] = addFace(gm,1);
tlist = linspace(1,60,60);
R1 = solve(model,tlist);
figure('units','normalized','outerposition',[0 0 1 1])
Tmax = max(max(R1.Temperature))
Tmin = min(min(R1.Temperature))
Riccardo Scorretti on 1 Apr 2022
In my opinion which is observed in the figure is just a uniform temperature T = 50, plus a numerical noise, which may be due to the mesh used to discretize the problem. Basically, the obtained result is that the temperature is constant. By the way, if you look carefully at your picture you can see that "hot spots" are centered in nodes.
I suggest you to check the formulation of your problem. For instance, try with a much higher convection coefficient, or add a source, or simulate over longer time -- just to see what happens. Play with the parameters of your formulation, and I'm sure you will find quickly the source of the problem.