Comparison 1D transient heat transfer with analytic solution wrong

More
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)
Last edit: 3 years 10 months ago by MaxxiKing. Reason: attachment missing?
More
3 years 9 months ago - 3 years 9 months ago #3502 by MaxxiKing
Hello,

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