Computational methodology for solving hydrodynamical equations for nematic liquid crystals