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.