I am trying to validate my simulation by comparing it to the analytical solution but the results do not match and I cannot find my mistake.
I set up a 1D-case of a plate with an initial temperature of 0 °C. On the left and right side of the plate I set the temperature to 20 °C and calculate the transient temperatures (see attached forum_question.py).
These results do not match the analytic solution (see comparison.png).
Thank you very much!
*(sorry there was a error with the attachments before)