The application of lenses combined with array antennas (also known as dome arrays or dome antennas) to the next generation of terrestrial and satellite communication systems brings a wide range of advantages in terms of improved radiation performance, reconfigurability in the use case, and reduction in power consumption. To facilitate the industrial implementation of dome antennas, highly efficient simulation tools are required. In this paper, we present a streamlined implementation of ray tracing for fast and efficient numerical analysis of the far-field radiation performance of 2D multilayer dielectric lenses combined with phased arrays. Unlike commercial physical-optical methods, our proposed ray-tracing method is capable of computing the effects of internal reflections in the dome in a multilayer configuration. In addition, the method estimates the absorption losses as a result of the Joule effect. To demonstrate the effectiveness of the proposed approach, we provide comparisons of the simulated radiation patterns using our proposed ray tracing with the results obtained from commercial full-wave simulation tools.