Improving initialization and evolution accuracy of cosmological neutrino simulations