We construct and analyze robust strong stability preserving IMplicit-EXplicit Runge-Kutta (IMEX RK) methods for models of flow with diffusion as they appear in astrophysics and in many other fields where equations with ...