The convergence of finite volume method (FVM) or discrete ordinate method (DOM) is known to degrade for optical thickness greater than unity and large scattering albedo. The present article presents a convergence acceleration procedure for the FVM based on a full approximation storage (FAS) multigrid method. Among a variety of multigrid cycles, the V-cycle is used and the full multigrid algorithm (FMG) is applied to an analysis of radiation in irregular two-dimensional geometry. Solution convergence is discussed for the several cases of various optical thickness and scattering albedo. At small scattering albedo and optical thickness, there is no advantage to using the multigrid method for calculation CPU time. For large scattering albedo greater than 0.5 and optical thickness greater than unity, however, the multigrid method improves the convergence and the solution is rapidly obtained.