This thesis presents the development, coding and testing of the comprehensive beam tracing method (CBTM). The present beam tracing methods for simulation of sound waves propagation do not compute the refraction, but only the reflection of sound waves. CBTM is designed for simulation of propagation of sound in non‐homogenous environments ‐ environments where the wave propagates through several media, or through one media with slightly different physical characteristics. This method has different scene structure than traditional simulations ‐ in order to enable definition of scenes with several non‐convex media, scene structure is based on boundary surfaces, shells and entities. Geometry of scene is defined using triangular irregular networks. The beams used for tracing the propagation of sound in scene have triangular section. During propagation they are adapted to surrounding geometry ‐ beam is divided when it encounters a discontinuum. This algorithm is developed specially for CBTM, and is thoroughly tested. It is based on process of intersection and subtraction of two planar triangles. The algorithm for intersection and subtraction of two planar triangles in CBTM differs from usual Southerland‐Hodgson algorithm. Construction of refracted beams in CBTM has shown slight errors in geometry of beams, which occur because refraction is a non‐linear transformation. These errors were analyzed ‐ the analysis showed that errors can be tolerated for refraction occurring between two media whose physical characteristics do not differ much, and for beams that are not wide. The final result of the simulation is the plane scan of simulated scene, which presents the level of sound intensity in raster of points in a plane. Level of sound intensity is calculated in a different fashion from traditional beam tracing simulations. This new way of calculation gives an improved performance, and possibility for extension of simulated wave phenomena, but at cost of a slight error. These errors were analyzed theoretically, and measured, with the conclusion that they can be neglected, because they are smaller than errors already incorporated in a simulation model. The CBTM was tested on three scenes, in order to check the results of simulations, and also to measure time and space performance of this method. CBTM was compared with two commercial simulations based on ray‐tracing model. The comparison showed that CBTM is accurate, and that performances of CBTM are comparable to the faster of two ray‐tracing simulations. The only drawback of CBTM that comparison discovered was that for high order of reflected sound the performance decreased exponentially. The biggest advantage of method is the very good performance in case of big resolution of plane scan, or in case of large simulated scenes.