Turns out the viscous flux was never executed because the wrong parameter was set in the input file. After I turned it on the solution changed. I updated the reference. @fwplatzek Can you make sure that the result still looks like what you would expect?