Modeling of sound propagation in the closed space and its interaction with obstacles/Garso sklidimo uzdaroje erdveje bei saveikos su kliutimis modeliavimas.

Mikalauskas, R. ; Volkovas, V.

1. Introduction

When the problems of noise reduction in the dwelling, industrial and other premises are being solved, the investigation on sound propagation in the closed space and its interaction with obstacles becomes indispensable. As essentially the sound propagation is the distribution of sound pressure, thus when the acoustic projecting of closed space (here the premises of various purposes are considered) is being implemented, it is often important to know and secure certain level of sound pressure in separate points. The values of this parameter in the analyzed object depend on a number of factors: geometry of the room, sound absorption of its walls, ceiling, floors, present things, nature of noise source, etc.

All this shows that the investigation on sound propagation in the closed space and its interaction with obstacles (for example, acoustic screens that limit closed space) is a complex task that needs theoretical modeling and experimental tests to be solved. Good choice or creation of theoretical model accelerates a lot the solution of the aforementioned task and guarantees the sufficiently accurate determination of the level of sound pressure in the closed space at any point. There are lots of bibliographical sources, which model the origin of acoustic field, while character of distribution in space and interaction with obstacles are described using the computational models of finite elements (FE), boundary elements (BE), finite differences (FD) and analytical models. When the acoustic tasks are being solved using the analytical method, the Kirchof-Helmholtz equation, the theory of geometrical diffraction, are used the most frequently. When the analytical acoustic models are being discussed, it is possible to state that they are more universal. However they describe the interaction of acoustic medium and structure approximately, more empirically [1]. The acoustic field modeled by them is diffusive, while the excitation source is point source. The method of finite differences is widely used to solve the problems of visualization of sound propagation in the rooms and reflection from the obstacles [2, 3]. When this method is used, the sound wave is expressed by the differential equations with partial derivatives. The main advantage of this method is relatively small need of computer resources, thus it is often used to model transient acoustic processes in the homogeneous field.

The method of boundary elements is also used to solve acoustic tasks [4, 5]. The environment is divided into direct (based on classical Helmholtz integral equation, where acoustic pressure and velocity of particles are the primary variables) and indirect (where primary variables are the differences of pressure and its direction's gradient beyond the boundary) [6, 7]. When this method is modeled, the discretization problem is encountered; in the area of high frequencies the number of elements has to be increased in order to receive accurate solution, and this increases significantly the usage of computer resources. When the interaction of partition with acoustic field is calculated, the plane of structure's center is used instead of its surface. In such a case when the tasks are solved by BE method, using the aforementioned indirect formulation, the thickness of the acoustic partition is not taken into account.

Another widely used method for formation of acoustic models is the method of finite elements. As well as in the case of finite differences, when this method is used, the wave's equation is being solved (taking into account boundary conditions) by dividing the space (in some cases, also the time) into elements. Then the wave's equation is expressed in the discreet set of linear equations for these elements. The method of finite elements could be also used to model the transfer of energy between separate surfaces, or energy exchange. The advantage of this method [8-10] is that it could relate directly the structural and acoustic mediums and under changing modeled environmental conditions evaluate their interaction, which is very important for formation of acoustic partition systems. This method is used to solve the three-dimensional tasks of acoustic medium, and the received results show completely the character of the acoustic field in the analyzed space. To note the disadvantages of this method, it shall be said that when the conditions of the modeled environment and excitation change, the model has to be made anew, and this usually needs a lot of time [11], while the solved tasks of free space are formed in the area of low frequencies. After the literature on the modeling on the basis of FEM had been reviewed, it became evident that this method is used to solve various acoustic tasks: investigation of acoustic properties of various materials, modeling of acoustic partition systems, investigation on sound propagation in various cavities, investigation on interaction between structural and acoustic mediums, etc. The main advantage of this method if compared to other ones, is that it could be used to model heterogeneous acoustic medium, to assess several excitation sources and to receive complete character of acoustic field in the analyzed space. The works [12] modeled two sources of different frequencies and sound pressures, which work separately and at the same time, as well as impulse (impact) excitation. While modeling the tasks of harmonic analysis and transient process were solved. The steady state acoustic fields received from separate excitation sources during harmonic analysis were summarized, according to the superposition principle, and in such a way the whole acoustic field was calculated that appears when several sources of different frequencies are acting at the same time. According to the received results, when the FE method formed by 2D model is used, it is possible to model the acoustic excitation that appears in the real conditions. However in order to model the sound propagation as adequately as possible in the closed space and its interaction with obstacles, besides the varying excitation conditions, the sound transmission through the obstacle should be taken into account, which is very important when the direct field of sound pressure is predominant if compared to the reflected field of sound pressure. Thus with regard to the aforementioned, the purpose of this work is

* to create the model of acoustic field that takes into account the sound propagation in the closed space and its interaction with obstacles in the real conditions, following FEM;

* to analyze the adequacy of this theoretical model for real fields and its application possibilities for designing of mobile and controlled systems of noise reduction.

2. Model of the sound propagation in the closed space on the basis of FEM

The interaction of structure and acoustic mediums in the formula of finite elements is described as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (1)

where [[M.sup.P.sub.e]], [[M.sub.e]] are matrixes of the mass of acoustic medium and structure accordingly; [[C.sup.P.sub.e]], [[C.sub.e]] are damping matrixes of acoustic medium and structure; [[K.sup.P.sub.e]], [[K.sub.e]] are stiffness matrixes of acoustic medium and structure; [[rho].sub.0][[[R.sub.e]].sup.T] is relation matrix of acoustic medium and structure; {[P.sub.e]} is vector of pressure in the nodes and its derivatives {[[??].sub.e]} {[[??].sub.e]} with regard to time; {[[??].sub.e]} is vector of nodal displacement and its derivatives {[[??].sub.e]} {[[??].sub.e]} with regard to time; {[F.sub.e]} is load vector; [[rho].sub.0] is density of air medium.

When the theoretical model was formed, the FEM software ANSYS 10 was used. The analyzed 2D model consisted of acoustic and structural mediums. In order to model them, the elements FLUID29 and PLANE42 were used. As the acoustic package FEM of the ANSYS 10 software does not take into account the loss of sound energy when the sound is transmitted through the obstacle, the methodology used for this model is specified in [13]. According to this methodology, when pressure of the incident sound wave is known, the loss of sound pressure is calculated when the wave passes from one medium to another, and the value of sound pressure is determined on the boundary of mediums--according to the scheme of the process shown in the Fig. 1, that would be the sound pressure on the junction of the second and third mediums. In such a way when the sound propagation in the closed space is modeled, firstly the system is excited by the sound source of certain size and frequency, and the field of sound pressure is determined in the closed space, as well as on the boundary between the incident wave and structural medium (boundary between the first and the second medium in the Fig. 1).

According to the present data, when the aforementioned methodology is used, the loss of sound pressure is calculated when the sound wave passes through the structural wave, as well as the values of sound pressure on the boundary of the second and third mediums (Fig. 1). These values are later used to excite and calculate sound pressure in the acoustic medium. Eventually the values of sound pressure calculated in the first stages are summarized using the principle of superposition, and acoustic field in the closed space with obstacle is calculated taking into account the loss of sound energy when the sound wave passes the obstacle. In order to computerize the calculations using the methodology [14], the ANSYS macro file was created in the programmable medium. This allowed making the calculations much faster.

[FIGURE 1 OMITTED]

The noise sources in this work are modeled in the area of moderate frequencies ([lambda]!1) i.e. in the area, to which the majority of industrial machines and equipment belongs and that is very topical in the machine acoustics. As the room of simple form (Testing laboratory of machine vibrations and acoustic noises of Technological Systems Diagnostics Institute) was used as the prototype of theoretical model, the method of point sources was used to model the noise sources. This method is simple and easily applicable. Two excitation sources were used, which frequencies were 1000 and 2000 Hz. The values of acoustic characteristics of walls, ceiling, floors and partition used in the model were taken from certain documentations of manufacturers: air density [rho] = 1.2 kg/[m.sup.3]; velocity of sound wave propagation c = 335 m/s; suppression coefficient of air sound [mu] = 0; density of sound suppression's partition [rho] = 750 kg/[m.sup.3]; elasticity module of sound suppression's partition E = 3.4 x [10.sup.9] Pa; velocity of sound propagation in partition's material [c.sub.P] = 675 m/s; coefficient of sound suppression's partition [mu] = 0.9. The received results of the theoretical experiment are presented below.

According to the received results of theoretical experiment, when the interaction of propagated sound with obstacles is being modeled taking into account the loss of sound pressure while passing through the obstacle, the distribution character of sound pressure in the area behind the partition changes (Figs. 2 and 3). These figures show the excitation when two sources of different frequencies are acting at the same time.

[FIGURE 2 OMITTED]

[FIGURE 3 OMITTED]

[FIGURE 4 OMITTED]

[FIGURE 5 OMITTED]

When the cases of excitation sources of different frequencies are compared (Figs. 4 and 5), it is seen that the loss of sound pressure generated by the source of 2000 Hz frequency while passing through the obstacle is bigger than in case of 1000 Hz source. Fig. 6 shows the values of sound pressure in different measurement points (the position of these points is shown in Fig. 7).

[FIGURE 6 OMITTED]

3. Experimental test

In order to check the adequacy of the theoretical model, the experiment was done, and its results were compared with the theoretical ones. In order to imitate the experimental test, the initial data were selected accordingly: value of sound pressure generated by sound source, its frequency, absorption coefficients of room's planes and partition, etc. The acoustic partition (acoustic screen) was used in the experimental test. During the experiment the values of sound pressure were measured in front of and behind the acoustic partition in certain points. The principal scheme of experimental measurement of sound pressure is shown in the Fig. 7.

Marking of points shown in this scheme means the coordinates of certain point, for example, the point marked as 3_025, has x coordinate of 3 m, and y coordinate -0.25 m, point 1_3 has x coordinate of 1 m, and y coordinate of 3 m, etc. The x coordinates of the first sound source shown in Fig. 7 are 5 m, y - 0.25 m; the x coordinate of the second sound source is 5 m, y - 1.2 m.

[FIGURE 7 OMITTED]

[FIGURE 8 OMITTED]

[FIGURE 9 OMITTED]

[FIGURE 10 OMITTED]

When the analysis of sound propagation in testing laboratory and its interaction with obstacles was done, the aforementioned method on the basis of FEM was used together with harmonic analysis, during which the harmonic excitation at certain values of sound pressure in the analyzed frequencies was performed. The acoustic field in the testing laboratory was created using two loudspeakers. The values of sound pressure in different points of testing laboratory in front of and behind the partition were measured using the device Investigator 2260 and analyzer of vibrations and noise PULSE 3560 [15]. The received results of theoretical and experimental test are presented below.

According to the Fig. 8, the results of the theoretical model that takes into account losses of sound pressure during sound transmission through the obstacle coincide significantly better with the experimental results, if compared to the previous model [12], where it was considered that a part of sound energy radiated by the source could only be absorbed or reflected. When the values of sound-pressure level are compared separately in the height of one and two meters (Figs. 9 and 10), the area of so called partition's "shadow" (place behind the partition) is clearly seen, where the values of sound pressure are significantly lower than these in front of the partition. The dependencies received during the experiment confirm this. Thus it is possible to state that the acoustic field created in the analyzed closed space is direct. When the values of sound pressure's level are compared at the measurement points 3_1 and 3_2, located immediately behind the partition at the height of one and two meters, it is seen that the level of sound pressure at point 3_1 is higher by several decibel than at point 3_2. This could be explained by the direction of sound emitted by the sources. That is, the point at the height of one meter, gets into the sound front emitted by the sources, and the waves that have passed through the obstacle and that has passed by it are summarized. The values of sound pressure's level of the points, which are furthest from the partition, which x coordinate is equal to one meter, at the heights of 0.25 m and 1 m are higher than of the points, which x coordinate is equal to 2 and 3 meters (Figs. 11 and 12).

[FIGURE 11 OMITTED]

[FIGURE 12 OMITTED]

The increase of the values of sound pressure level at the aforementioned points could be explained by the fact that the closer they are to the surfaces reflecting the sound well (in this case these are the floors of testing laboratory and one of its walls), the reflected sound wave is summarized with the waves that have passed through the obstacle and that have passed by it. The data of experiment (Fig. 11) confirm these results received with the help of theoretical model (Fig. 12). The theoretical test has analyzed additionally the interaction of the sound propagated in the closed space with the obstacle, depending on its thickness. The received results are presented below.

According to Fig. 13, different thickness of partition affects the character of sound pressure in the area behind the partition. When thickness of the partition is 0.1 m (Fig. 13, a), the gradient of sound pressure acquires bigger values than in case when the thickness is 0.3 m (Fig. 13, b). It is noticed that in the case of thinner partition the gradient of sound pressure is the biggest in the directions of the wave that has passed through the obstacle and that has passed by it, while in the case of thicker partition sound pressure varies the most in the direction of the wave that has passed by the obstacle. The values of sound pressure level calculated with the help of the model at certain points (Fig. 14) show that they are significantly smaller in the case of thicker partition.

This shows that changes of geometrical parameters of the partition can change noticeably the character of acoustic fields and the values of sound pressure at certain points. These results confirm the statements specified in the bibliographical source [15].

[FIGURE 13 OMITTED]

[FIGURE 14 OMITTED]

4. Conclusions

The model of acoustic field that evaluates sound propagation in closed space and its interaction with obstacles in real conditions was created on the basis of FEM. This model differs from the previous one [12], because it additionally evaluates the sound transmission through the obstacle when several excitation sources of different frequency are used at the same time. The test of adequacy of the theoretical model showed that it reproduces real fields adequately and that it can be used to analyze the sound propagated in closed space and its interaction with obstacles, as well as to apply it for designing of mobile and controlled systems of noise reduction.

Acknowledgement

This work was supported by Lithuanian State Scientific and Study fund, project T-87/09.

Received September 21, 2009 Accepted December 03, 2009

References

[1.] Gensei Matsumoto, Kyoji Fujiwara, Akira Omoto. Directivity of the sound radiated from a factory building.--Acoust. Sci. & Tech., 2001, 22. 6, p.434-436.

[2.] Takatoshi Yokota, Shinichi Sakamotoy, Hideki Tachibanaz. Sound field simulation method by combining finite difference time domain calculation and multi-channel reproduction technique.--Acoust. Sci. & Tech., 2004, 25. 1, p.15-23.

[3.] Takatoshi Yokota, Shinichi Sakamoto, Hideki Tachibana. Visualization of sound propagation and scattering in rooms.--Acoust. Sci. & Tech., 2002, 23. 1, p.40-46.

[4.] Ciskowski R. D. and Brebbia C. A. Boundary Element Methods in Acoustics.--New York: Elsevier Applied Science, 1991.-290p.

[5.] Filippi, D. Habault, J. P. Lefevre, and Bergassoli, A. Acoustics, Basic Physics, Theory and Methods.--New York: Academic, 1999.-317p.

[6.] Kludszuweit. Time iterative boundary element method ~TIBEM--a new numerical method of four-dimensional system analysis for the calculation of the spatial impulse response.--Acustica, 1991, 75, p.17-27 (in German).

[7.] Kopuz, S. and Lalor, N. Analysis of interior acoustic fields using the finite element method and the boundary element method.--Appl. Acoust., 1995, 45, p.193-210.

[8.] Everstine, G.C. Finite element formulation of structural acoustics problems.--Comp. & Structures, 1997, 65, p.307-321.

[9.] Ihlenburg, F. Finite Element Analysis of Acoustic Scattering.--New York: Springer, 1998.-132p.

[10.] Morand, H. J.-P. and Ohayon, R. Fluid Structure Interaction: Applied Numerical Methods.--Wiley, Chichester, UK, 1995.-224p.

[11.] Tsingos, N. and Gascuel, J.-D. Soundtracks for computer animation: Sound rendering in dynamic environments with occlusions.--GraphicsInterface'97, May 1997, p.9-16.

[12.] Mikalauskas, R., Volkovas, V. Investigation of adequacy of the acoustical field model.--Mechanika.--Kaunas: Technologija, 2009, Nr.2(76), p.46-49.

[13.] Randall, F. Barron. Industrial Noise Control and Acoustics.--Marcel Dekker, Inc.--New York, 2001, (electronic source).

[14.] Tumonis, L., Schneider, M., Kacianauskas, R., Kaceniauskas, A. Comparison of dynamic behaviour of EMA-3 railgun under differently induced loadings.--Mechanika.--Kaunas: Technologija, 2009, Nr.4(78), p.31-37.

[15.] Volkovas, V., Slavickas, E.S., Gulbinas, R.J. The investigation of effectiveness of acoustical isolation of noise sources.--Mechanika 2009: Proceedings of 14th International Conference, April 2-3, 2009, Kaunas, 2009, p.445-448.

R. Mikalauskas *, V. Volkovas **

* Kaunas University of Technology, Technological System Diagnostic Institute, Kestucio str. 27, 44312 Kaunas, Lithuania, E-mail: [email protected]

** Kaunas University of Technology, Technological System Diagnostic Institute, Kestucio str. 27, 44312 Kaunas, Lithuania, E-mail: [email protected]