Investigating simultaneous effects of temperature, surface heterogeneity and geometry on fluid mixing in electroosmotic flow considering temperature dependent properties by Nernst-Planck Poisson method.