Spatial Distribution Imaging of Steady-state Sound Field with the Ray-tracing Method
Jiang Zhong-Jin,Qiu Xiao-Jun
DOI: https://doi.org/10.3321/j.issn:0469-5097.2007.01.007
2007-01-01
Abstract:In the traditional ray-tracing prediction, only the sound field at a single point is predicted through one simulation. It is difficult to predict the whole sound field in a two-dimensional plane or a three-dimensional space. The sound field at the test point is predicted by receiving the sound rays passing nearly by the point with a receiving sphere. However, how to determine the receiving radius is still a debatable problem, because unsuitable receiving radius will result in prediction error and different sound fields need different receiving radius. In this paper a new ray-tracing method is proposed, which is able to predict and image the spatial distribution of a steady-state sound field through once simulation. The receiving sphere is not used because the prediction is not for a single point, for the whole sound field. First the sound intensity spatial distribution in the whole space is predicted and memorized in a three-dimensional energy array. Then the sound intensity level (SIL) and the sound pressure level (SPL) along a line or in a two-dimensional plane can be calculated directly from the energy array and then imaged. In the new ray-tracing method, the initial sound rays are uniformly emitted from the point source with prefixed elevation angle and azimuth angle. When a directional sound source is simulated, the directivity coefficient can be used to modify the initial energy of the sound rays. The propagation of the sound rays includes specular reflection and diffuse reflection as in the Lambert diffusion model. Every sound ray is terminated by the energy discontinuity percentage (EDP) instead of the number of reflections. When a sound space has been established in the new ray-tracing method, a three-dimensional array is also established for the sound space, which can be defined as energy array here. The numbers of the rows, columns and layers of the energy array are proportional to the width, length and height of the sound space respectively. The proportion coefficient can be determined by the needed spatial resolution in prediction. Every element in the energy array is used to memorize the sound intensity of a certain little cube in the sound space, and the position of the element in the energy array is equal to the position of the cube in the sound space. When a sound ray propagates in the sound space, its travel path is traced. The path will be discretized into a series of evenly and compactly spaced nodes, and the coordinates of every node will be calculated to find the element with the same position in the energy array. The sound intensity at this node will be recorded in this element. When a sound ray has been traced to termination, the sound intensity along the travel path has been all memorized into the energy array. It seems as if the sound ray had traveled in the energy array and the sound intensity has been recorded along its travel path. If an element is passed through by the sound ray for two times, the latter sound intensity should be added to the former sound intensity. Different sound rays should be traced in the same way, and the later sound intensity should be added to the previous sound intensity in the energy array. When all the sound rays have been traced over, the spatial distribution of the sound intensity has been predicted and restored in the energy array. Each element in the energy array represents the sound intensity of a little cube in the sound space. This is an approximation method based on the ray-tracing statistical idea, while the prediction accuracy can still be very high when the proportion coefficient of the energy array is very large. When the sound intensity in a certain plane, along a certain line or at a certain point is interesting, it can be extracted from the energy array directly and then imaged. In the new ray-tracing prediction, the energy array is only used to record the spatial distribution of the steady-state sound field, so it can not show the time behaviour of the sound field. In fact the new ray-tracing method can be expanded in different ways to simulate the non-steady sound field. For example, the three-dimensional energy array can be used to record the energy distribution in a two-dimensional plane during a time span. The row and column represent the X and Y spatial axis respectively, and the layer represents the travel time. When a sound ray is traced, the coordinate, travel time and energy are all calculated and recorded. So all the functions of the conventional ray-tracing method can be recovered and the impulse response of every point in the whole plane can be calculated through only once simulation. Besides, a four-dimensional array can also be used to simulate the spatial distribution in a sound space during a time span, while the computation load will be larger. The steady-state sound field spatial distribution in a long rectangular space was measured and predicted with the new ray-tracing method. The predicted results and the measured results were compared with each other. Two experiments were performed. In the former experiment the long space was an enclosure with reflecting ends, in which the sound field was nearly reverberating. In the latter experiment the long space was with absorbent ends to study the sound field of tunnel spaces. The measured results in the long space were consistent well with the predicted results, so the validity of the new method was verified.