How to run axisymmetric heat conduction simulation

[richmondodufisan]

Hi, I am trying to run a simple heat conduction problem where I reduce the full 3D model to an axisymmetric 2D one. I am sure I'm missing something but I'm not sure what. I thought I'd ask here first before I get into the weeds of radial coordinates/working it out manually.

The problem setup is simple: a semi-infinite box with homogeneous material properties (thermal conductivity) has a gaussian heat flux applied at the top surface. It is steady state too, so no worry about the time terms for now.

Below is the temperature profile from the top of the box to the bottom at r = 0 (i.e x=0, y=0).

I implemented the exact same simulation in 2D with a mesh that should be a slice of the 3D version. But my result differs. Not just in magnitude but the temperature profile itself is different. edit, it's based on the way derivatives are calculated in the rz coordinate system.

Is there an existing model/kernel that can help solve a problem like this/is there something simple I'm missing?

[richmondodufisan]

Solved, MOOSE has built in functionality, no need to reinvent the wheel

[Problem]
coord_type = RZ
[]

Although I got a warning that this is deprecated and should be under the mesh block, but I got an error when I did that

[GiudGiud]

this should work

[Mesh]
  coord_type = RZ
  [other_mg]
      ...
  []
[]

[richmondodufisan]

That worked thanks!