the old version of the solution of 2.12 was obtained when the boundary conditions where switched to dirichlet for the energy equation.
Additionally old implementation of molar density is necessary to obtain solution of 2.12