We investigate dissipative full discretizations for the equations of motion associated with models of flow and radiative transport inside stars. We derive dissipative space discretizations for the parabolic part and demonstrate that together with specially adapted total-variation-diminishing implicit Runge-Kutta time integrators with adaptive step-size control this yields reliable and efficient integrators for the underlying high-dimensional nonlinear evolution equations.