This project focuses on the propagation of waves in a medium where the properties are heterogeneous, randomly-fluctuating for instance. In general, the solution of such an equation in a given position displays two phases: (i) a coherent wave, somewhat similar to a wave propagating in a homogeneous medium, and (ii) an incoherent wave, seemingly more random. In some very heterogeneous media, like concrete, granular media or the Earth [1], that second phase of the solution is the most prominent. The amplitude of the solution can then be represented by a radiative transfer equation (RTE), obtained through an asymptotic expansion assuming random fluctuations of the properties in the appropriate regime [2]. The RTE only models the amplitude of the solution, not its phase, but has the advantage of homogenizing all the rapid fluctuations of the properties and solution. This paper reports a series of numerical experiments comparing solutions of the 3D acoustic wave equation in a heterogeneous medium and of the radiative transfer equations. Parameters of the two equations are chosen such that the radiative transfer solution is expected to provide an accurate approximation of the energy of the wave in the weak scattering regime. The comparisons indicate that the radiative transfer provides accurate approximations even quite far from that regime. A particular attention is devoted to analyzing the results close to boundaries [3,4], where the accuracy of the radiative transfer equation has not been evaluated before, while this is exactly where measurements are accessible.