We present a multilevel method for solving multigroup neutron transport k-eigenvalue problems in two-dimensional Cartesian geometry. It is based on the nonlinear diffusion acceleration (NDA) method. The multigroup low-order NDA (LONDA) equations are formulated on a sequence of energy grids. Various multigrid cycles are applied to solve the hierarchy of multigrid LONDA equations. The algorithms developed accelerate transport iterations and are effective in solving the multigroup NDA low-order equations. We present numerical results for model reactor-physics problems with a large number of groups to demonstrate the performance of different iterative schemes.