Forum Message

 

 

We have moved the forum to https://forum.ngsolve.org . This is an archived version of the topics until 05/05/23. All the topics were moved to the new forum and conversations can be continued there. This forum is just kept as legacy to not invalidate old links. If you want to continue a conversation just look for the topic in the new forum.

Notice

The forum is in read only mode.

Comparison 1D transient heat transfer with analytic solution wrong

More
3 years 3 months ago - 3 years 3 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 3 months ago by MaxxiKing. Reason: attachment missing?
More
3 years 2 months ago - 3 years 2 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 2 months ago by MaxxiKing.
Time to create page: 0.122 seconds