Dear kest,
you might be interested in implementing the same testcase with OpenLB as a reference.
There is also a double-distribution thermal model included in OpenLB.
Similar code snippets should appear within the dynamics which are for example showcased in some of the examples/thermal/ folders.
BR
stephan