- Thank you received: 0
Comparison 1D transient heat transfer with analytic solution wrong
3 years 10 months ago - 3 years 10 months ago #3498
by MaxxiKing
Dear ngsolve-Team,
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)
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)
Attachments:
Last edit: 3 years 10 months ago by MaxxiKing. Reason: attachment missing?
3 years 9 months ago - 3 years 9 months ago #3502
by MaxxiKing
Replied by MaxxiKing on topic Comparison 1D transient heat transfer with analytic solution wrong
Hello,
found the error: I mixed up 'region_hcs' and 'region_conductivities' in line 122.
Correct is:
Now the solution fits also the analytic solution.
found the error: I mixed up 'region_hcs' and 'region_conductivities' in line 122.
Correct is:
Code:
m = BilinearForm(fes, symmetric=False)
m += SymbolicBFI(CoefficientFunction(region_densities * region_hcs) * u * v)
m.Assemble()
Now the solution fits also the analytic solution.
Last edit: 3 years 9 months ago by MaxxiKing.
Time to create page: 0.098 seconds