Thermal LBM implementation using current palabos libraries

This refers to the Thermal2D example presented in show cases
I would be approaching a similar problem for my research, however it will be a pure heat conduction problem in a porous media.
In the example the temperature field was modelled using an advective-diffusive behaviour, however I want to approach the heat conduction problem using a double distribution approach(http://www.sciencedirect.com/science/article/pii/S0021999198960570), which essentially approaches the thermal problem by solving an additional internal energy density distribution evolution equation along with the standard density distribution one. This involves evaluation of an additional IE distribution function.
I was wondering if palabos can handle that, or if I will need to make it extensible by adding/changing libraries, will that be possible or not. Plus I would like to know if I do approach the problem as in Thermal2D example, how much will that be valid for pure heat conduction.

Regards,

An LBM enthusiast