We present a parameter-free method for high Reynolds number turbulent flow. The method is based on adaptive mesh optimization using adjoint techniques. We connect the numerical approximation to dissipative weak solutions. Turbulent boundary layers are left unresolved, no boundary layer mesh is needed. We highlight benchmark problems for which the method is validated.