This paper calculates light scattering by a spherical water particle containing densely packed inclusions at a visible wavelength 0.55 μm by a combination of ray-tracing and Monte Carlo techniques. While the individual reflection and refraction events at the outer boundary of a sphere particle are considered by a ray-tracing program, the Monte Carlo routine simulates inteal scattering processes. The main advantage of this method is that the shape of the particle can be arbitrary, and multiple scattering can be considered in the inteal scattering processes. A dense-medium light-scattering theory based on the introduction of the static structure factor is used to calculate the phase function and asymmetry parameters for densely packed inclusions. Numerical results of the single scattering characteristics for a sphere containing multiple densely packed inclusions are given.