Thermal PDE in cylindrical coordinates at the origin

9 次查看（过去 30 天）
Tobias
Tobias 2019-1-11
回答： Tobias 2019-7-29
Hi,
I have a problem with a thermal PDE in cylindrical coordinates. I followed the example in Heat Distribution in Circular Cylindrical Rod .
The script runs without errors, however the solution has peculiar spot at z = 0, r = 0 where the Temperature is constantly zero. The heat source at that spot is well larger than zero and all around the spot a temperature which is as expected is reached.
I tried changing several of the pde SolverOptions or choose the mesh differently, however nothing of what I tried was successful.
Does anybody have an idea what I can do to fix this issue?
Thanks in advance

请先登录，再回答此问题。

采纳的回答

Tobias
Tobias 2019-7-29
My error was in not transforming the boundary condition correctly. I had assumed that the BC in the example just had a curios form, but the *region.y originated from the transfer to cylindrical coordinates.
With the correct BC the solution makes sense on the whole domain.

更多回答（2 个）

Svetlana Pease
Svetlana Pease 2019-1-11
Hi Tobias,
Are you talking about the steady-state or transient solution? The only zero temperature I get is for the transient solution at t = 0 (time step is 1 for t = 0).
Regards,
Svetlana Pease
MathWorks Documentation Group
  5 个评论
显示 3更早的评论
Tobias
Tobias 2019-1-18
Is there a possibility to create the mesh with two different Hmax? I naively tried to assembly the mesh in two steps but I was not allowed to change the read-only FEMesh.
If I calculate my problem in full 3d the problem does not appear (however I cannot choose my grid fine enough to get rid of other oscillations)
Torsten
Torsten 2019-1-18
I'd simply start with the geometry at r = eps instead of r = 0 for eps being a small value.

请先登录，再进行评论。

Svetlana Pease
Svetlana Pease 2019-1-18
Hi Tobias,
No, there is no option to have two or more Hmax settings for the same geometry/mesh.
Regards,
Svetlana

类别

在 Help Center 和 File Exchange 中查找有关 Heat Transfer 的更多信息

标签

产品


版本

R2018b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by