Following a Gallavotti's conjecture, stationary states of Navier-Stokes fluids are proposed to be described equivalently by alternative equations besides the Navier-Stokes equation itself. We discuss a model system symmetric under time reversal based on the Navier-Stokes equations constrained to keep the enstrophy constant. It is demonstrated through highly resolved numerical experiments that the reversible model evolves to a stationary state which reproduces quite accurately all statistical observables relevant for the physics of turbulence extracted by direct numerical simulations (DNS) at different Reynolds numbers. The possibility of using reversible models to mimic turbulence dynamics is of practical importance for the coarse-grained version of Navier-Stokes equations, as used in large-eddy simulations. Furthermore, the reversible model appears mathematically simpler, since enstrophy is bounded to be constant for every Reynolds number. Finally, the theoretical interest in the context of statistical mechanics is briefly discussed.