A real-time algorithm for rendering of an inhomogeneous scattering media such as smoke under dynamic low-frequency environment lighting is described. An input media animation is represented as a sequence of density fields, each of which is decomposed into a weighted sum of a set of radial basis functions (RBFs) and an optional residual field. Source radiances from single and optionally multiple scattering are directly computed at only the RBF centers and then approximated at other points in the volume using an RBF-based interpolation. Unique approximation techniques are introduced in the computational algorithms to simplify and speed up the computation of source radiance contributed by single and multiple scattering. Using the computed source radiances, a ray marching technique using slice-based integration of radiance along each viewing ray may be performed to render the final image.