Update solver_direct_mean.cpp

Heatflux computation fixed by using UnitNormal, rather than Normal