In this paper, the computer modelling application based on the modal expansion method is developed to study the influence of a sound source location on a steady-state response of coupled rooms. In the research, an eigenvalue problem is solved numerically for a room system consisting of two rectangular spaces connected to one another. A numerical procedure enables the computation of shape and frequency of eigenmodes, and allows one to predict the potential and kinetic energy densities in a steady-state. In the first stage, a frequency room response for several source positions is investigated, demonstrating large deformations of this response for strong and weak modal excitations. Next, a particular attention is given to studying how the changes in a source position influence the room response when a source frequency is tuned to a resonant frequency of a strongly localized mode.

A simple analytical method is developed to estimate frequencies of longitudinal modes in closed hard-walled ducts with discontinuities in a cross-sectional area. The approach adopted is based on a general expression for the acoustic impedance for a plane wave motion in a duct and conditions of impedance continuity at duct discontinuities. Formulae for mode frequencies in a form of transcendental equations were found for one, two and three discontinuities in a duct cross-section. An accuracy of the method was checked by a comparison of analytic predictions with calculation data obtained by use of numerical implementation based on the forced oscillator method with a finite difference algorithm.

A theoretical method has been presented to describe sound decay in building enclosures and to simulate the room impulse response (RIR) employed for prediction of the indoor reverberation characteristics. The method was based on a solution of wave equation having the form of a series whose time-decaying components represent responses of acoustic modes to an impulse sound source. For small sound absorption on room walls this solution was found by means of the method of variation of parameters. A decay function was computed via the time-reverse integration of the squared RIR. Computer simulations carried out for a rectangular enclosure have proved that the RIR function reproduces the structure of a sound field in the initial stage of sound decay suffciently well. They have also shown that band-limitedness of the RIR has evident influence on the shape of the decay function and predicted decay times.

Elżbieta M. Walerian, Ph.D., D.Sc., a retired employee of the Institute of Fundamental Technological Research of the Polish Academy of Sciences (IPPT PAN), passed away after a serious illness, on the 26th December 2013. She was one of the scientific leaders in the Section of Environmental Acoustics of IPPT PAN and her career, educational and organizational activities were inseparably linked with the acoustics. Elżbieta Walerian was born on August 9th 1950 in Poznań. She graduated from the Faculty of Mathematics, Physics and Chemistry of the Adam Mickiewicz University in Poznań, receiving her Master of Science degree in the environmental acoustics in 1973. Five years later, under the supervision of Professor Ignacy Malecki, she obtained her PhD title, in the physical acoustics, in IPPT PAN in Warsaw. In 1979 she began working at the Section of Environmental Acoustics of IPPT PAN, where she dealt with the diffraction of acoustic waves and a description of the sound field produced by vehicles moving in an urban area.

A new method for determining optimum dimension ratios for small rectangular rooms has been presented. In a theoretical model, an exact description of the room impulse response was used. Based on the impulse response, a frequency response of a room was calculated to find changes in the sound pressure level over the frequency range 20–200 Hz. These changes depend on the source and receiver positions, thus, a new metric equivalent to an average frequency response was introduced to quantify the overall sound pressure variation within the room for a selected source position. A numerical procedure was employed to seek a minimum value of the deviation of the sound pressure level response from a smooth fitted response determined by the quadratic polynomial regression. The most smooth frequency responses were obtained when the source was located at one of the eight corners of a room. Thus, to find the best possible dimension ratios, in the numerical procedure the optimal source position was assumed. Calculation results have shown that optimum dimension ratios depend on the room volume and the sound damping inside a room, and for small and medium volumes these ratios are roughly 1 : 1.48 : 2.12, 1 : 1.4 : 1.89 and 1 : 1.2 : 1.45. When the room volume was suitably large, the ratio 1 : 1.2 : 1.44 was found to be the best one.

Reverberant responses are widely used to characterize acoustic properties of rooms, such as the early decay time (EDT) and the reverberation times T20 and T30. However, in real conditions a sound decay is often deformed by background noise, thus a precise evaluation of decay times from noisy room responses is the main problem. In this paper this issue is examined by means of numerical method where the decay times are estimated from the decay function that has been determined by nonlinear polynomial regression from a pressure envelope obtained via the discrete Hilbert transform. In numerical experiment the room responses were obtained from simulations of a sound decay for two-room coupled system. Calculation results have shown that background noise slightly affects the evaluation of reverberation times T20 and T30 as long as the signal-to-noise ratio (SNR) is not smaller than about 25 and 35 dB, respectively. However, when the SNR is close to about 20 and 30 dB, high overestimation of these times may occur as a result of bending up of the decay curve during the late decay.

Accurate definition of boundary conditions is of crucial importance for room acoustic predictions because the wall impedance phase angle can affect the sound field in rooms and acoustic parameters applied to assess a room reverberation. In this paper, the issue was investigated theoretically using the convolution integral and a modal representation of the room impulse response for complex-valued boundary conditions. Theoretical considerations have been accompanied with numerical simulations carried out for a rectangular room. The case of zero phase angle, which is often assumed in room acoustic simulations, was taken as a reference, and differences in the sound pressure level and decay times were determined in relation to this case. Calculation results have shown that a slight deviation of the phase angle with respect to the phase equal to zero can cause a perceptual difference in the sound pressure level. This effect was found to be due to a change in modal frequencies as a result of an increase or decrease in the phase angle. Simulations have demonstrated that surface distributions of decay times are highly irregular, while a much greater range of the early decay time compared to the reverberation time range indicates that a decay curve is nonlinear. It was also found that a difference between the decay times predicted for the complex impedance and real impedance is especially clearly audible for the largest impedance phase angles because it corresponds approximately to 4 just noticeable differences for the reverberation metrics.