This work is devoted to numerical approximation of radiative transfer models. On the one hand, we focus on the discrete ordinates model. In order to couple this phenomena with slower ones, accurate and efficient numerical methods for long times are recquired. From a double time-space approximation of the solution, a high order GRP type scheme is developped with unrestricted time steps for hyperbolic linear systems on unstructured meshes. This scheme is then extended to discrete ordinates model. On the other hand, we focus on moment models of radiative transfer. Actually, in many applications, they remain a lot of properties from the RTE and give a sufficient approximation of the solution. Once the Riemann problem of the grey M1 model is solved, the numerical approximation of the multigroupeM1 model is considered. A particular attention is paid on the calculation of opacity means and closure laws. A precalculation algorithm is developped. The last application is concerned with an extension of radiative transfer to estimate the dose in radiotherapy. Unlike the usual greyM1 model, the space dependence in the fluxes is not necessary smooth. Thanks to changings of variables, a backward HLL scheme is developped. Many examples illustrate the interest of the obtained schemes.