thank you once again for your help. Now I have created a simple python script to simulate such an inhomogeneous heat conductivity problem over time (see example_forum.py).
The example is about a brick element which is heated by a constant heat flow into the element.
However, I wonder about the solution which seems to be unphysical. If I choose the values
scal_x = 1
scal_y = 0.01 #1
scal_z = 0.1 #1
to account for direction dependent heat conductivity, I get temperatures (at intermediate points in time) which are below the initial temperature. This cannot be correct since the heat flow into the brick is always positive.
The more serious issue is the missing comparison principle of finite elements, even of of lowest order methods. Thus, also positive sources may lead to negativ temperature. This failure is much more serious for very anisotropic conductivities.
But it works for lowest order mixed or hybrid mixed methods (I think for anisotropic as well) !
Have a look into sections