Mesh adaptation and higher order extrapolation of the Reynolds Averaged Navier-Stokes equations using tau-estimation