Field experiment on ground-to-ground sound propagation from a directional source

Similar documents
An experimental evaluation of a new approach to aircraft noise modelling

Environmental Noise Propagation

6-channel recording/reproduction system for 3-dimensional auralization of sound fields

The spatial structure of an acoustic wave propagating through a layer with high sound speed gradient

Multi-spectral acoustical imaging

A Road Traffic Noise Evaluation System Considering A Stereoscopic Sound Field UsingVirtual Reality Technology

Please refer to the figure on the following page which shows the relationship between sound fields.

Field noise measurement in the huge industrial plants for accurate prediction

From concert halls to noise barriers : attenuation from interference gratings

MEASURING SOUND INSULATION OF BUILDING FAÇADES: INTERFERENCE EFFECTS, AND REPRODUCIBILITY

Pre-Construction Sound Study. Velco Jay Substation DRAFT. January 2011 D A T A AN AL Y S IS S OL U T I ON S

METHODOLOGY FOR VERIFICATION OF SOFTWARE FOR NOISE ATTENUATION CALCULATION ACCORDING TO ISO STANDARD

Problems with the INM: Part 2 Atmospheric Attenuation

INTRODUCTION. Reducing noise annoyance. Aircraft noise is a global problem. First, we have to know how sound is emitted and propagated

Noise attenuation directly under the flight path in varying atmospheric conditions

Validation and evolution of the road traffic noise prediction model NMPB-96 - Part 1: Comparison between calculation and measurement results

Validation of lateral fraction results in room acoustic measurements

ENGINEERING STAFF REPORT. The JBL Model L40 Loudspeaker System. Mark R. Gander, Design Engineer

Analysis on Acoustic Attenuation by Periodic Array Structure EH KWEE DOE 1, WIN PA PA MYO 2

HIGH FREQUENCY INTENSITY FLUCTUATIONS

JBL Professional Application Note. Loudspeaker Array Low-Frequency Pattern Control using Filtered Array Technology

Modeling Diffraction of an Edge Between Surfaces with Different Materials

inter.noise 2000 The 29th International Congress and Exhibition on Noise Control Engineering August 2000, Nice, FRANCE

Performance of Roadside Sound Barriers with Sound Absorbing Edges

Sound source localization accuracy of ambisonic microphone in anechoic conditions

Convention Paper 6274 Presented at the 117th Convention 2004 October San Francisco, CA, USA

Unidirectional Sound Signage for Speech Frequency Range Using Multiple-Loudspeaker Reproduction System

Influence of ground reflections and loudspeaker directivity on measurements of in-situ sound absorption

group D DSA250 Specifications 2-WAY FULL-RANGE DIGITALLY STEERABLE ARRAY See TABULAR DATA notes for details CONFIGURATION Subsystem Features

Rec. ITU-R P RECOMMENDATION ITU-R P PROPAGATION BY DIFFRACTION. (Question ITU-R 202/3)

19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, 2-7 SEPTEMBER 2007 ACOUSTICAL ASPECTS OF THE SAGRADA FAMILIA CHURCH

Is the 2.5 db(a) Façade Correction for Road Traffic Noise Correct?

Investigation of Noise Spectrum Characteristics for an Evaluation of Railway Noise Barriers

Audio Engineering Society. Convention Paper. Presented at the 119th Convention 2005 October 7 10 New York, New York USA

Localizing Noise Sources on a Rail Vehicle during Pass-by

FINAL REPORT. On Project Supplemental Guidance on the Application of FHWA s Traffic Noise Model (TNM) APPENDIX L Tunnel Openings

ENHANCED PRECISION IN SOURCE LOCALIZATION BY USING 3D-INTENSITY ARRAY MODULE

Black. LWECS Site Permit. Stearns County. Permit Section:

ROOM SHAPE AND SIZE ESTIMATION USING DIRECTIONAL IMPULSE RESPONSE MEASUREMENTS

DISTANCE CODING AND PERFORMANCE OF THE MARK 5 AND ST350 SOUNDFIELD MICROPHONES AND THEIR SUITABILITY FOR AMBISONIC REPRODUCTION

SOUND PROPAGATION OUTDOORS. T.F.W. Embleton, J.E. Piercy, N. Olson Division of Physics, National Research Council Ottawa, Ontario K1A 0S1

Active Control of Energy Density in a Mock Cabin

ECMA-108. Measurement of Highfrequency. emitted by Information Technology and Telecommunications Equipment. 5 th Edition / December 2010

THE USE OF VOLUME VELOCITY SOURCE IN TRANSFER MEASUREMENTS

Composite square and monomial power sweeps for SNR customization in acoustic measurements

APPENDIX A SOUND PROPAGATION THEORY & METHODOLOGIES

THE PATTERNS OF THE SOUND INTENSITY DISTRIBUTION OF MIDRANGE LOUDSPEAKER

Acoustic Based Angle-Of-Arrival Estimation in the Presence of Interference

WIND SPEED ESTIMATION AND WIND-INDUCED NOISE REDUCTION USING A 2-CHANNEL SMALL MICROPHONE ARRAY

ECMA-108. Measurement of Highfrequency. emitted by Information Technology and Telecommunications Equipment. 4 th Edition / December 2008

Development of a sonic boom measurement system at JAXA

Experimental Investigation on the Effect of Origami Geometry on the Acoustic Characteristics

SODAR- sonic detecting and ranging

MODEL STUDY OF DOUBLE TRAFFIC NOISE BARRIERS

Selecting the right directional loudspeaker with well defined acoustical coverage

Proceedings of Meetings on Acoustics

ISO INTERNATIONAL STANDARD

inter.noise 2000 The 29th International Congress and Exhibition on Noise Control Engineering August 2000, Nice, FRANCE

IS INTERNATIONAL STANDARD. Acoustics - Attenuation of sound during propagation outdoors - Part 2: General method of calculation

FINAL REPORT. On Project Supplemental Guidance on the Application of FHWA s Traffic Noise Model (TNM) APPENDIX K Parallel Barriers

Antennas and Propagation

Measuring Instrument Combinations

ECE 476/ECE 501C/CS Wireless Communication Systems Winter Lecture 6: Fading

The influence of non-audible plural high frequency electrical noise on the playback sound of audio equipment (2 nd report)

Reducing the influence of microphone errors on in- situ ground impedance measurements

Acoustic Calibration Service in Automobile Field at NIM, China

Optimization of an Acoustic Waveguide for Professional Audio Applications

RECOMMENDATION ITU-R P Prediction of sky-wave field strength at frequencies between about 150 and khz

Antennas & Propagation. CSG 250 Fall 2007 Rajmohan Rajaraman

Airborne Sound Insulation

Antennas and Propagation. Chapter 5

GNSS Ocean Reflected Signals

inter.noise 2000 The 29th International Congress and Exhibition on Noise Control Engineering August 2000, Nice, FRANCE

A comparing overview on ECAC Doc.29 3 rd Edition and the new German AzB

TBM - Tone Burst Measurement (CEA 2010)

A Comparison of the Convolutive Model and Real Recording for Using in Acoustic Echo Cancellation

Speech and Audio Processing Recognition and Audio Effects Part 3: Beamforming

EFFECT OF ARTIFICIAL MOUTH SIZE ON SPEECH TRANSMISSION INDEX. Ken Stewart and Densil Cabrera

not overpower the audience just below and in front of the array.

DESIGN OF VOICE ALARM SYSTEMS FOR TRAFFIC TUNNELS: OPTIMISATION OF SPEECH INTELLIGIBILITY

THE RELATIONSHIP BETWEEN FILL-DEPTHS BASED ON GIS ESTIMATION, EARTHQUAKE DAMAGE AND THE MICRO-TREMOR PROPERTY OF A DEVELOPED HILL RESIDENTIAL AREA

ECE 476/ECE 501C/CS Wireless Communication Systems Winter Lecture 6: Fading

Development of multichannel single-unit microphone using shotgun microphone array

University Tunku Abdul Rahman LABORATORY REPORT 1

Radiowave Propagation Prediction in a Wind Farm Environment and Wind Turbine Scattering Model

Multiple Sound Sources Localization Using Energetic Analysis Method

ECE 476/ECE 501C/CS Wireless Communication Systems Winter Lecture 6: Fading

A White Paper on Danley Sound Labs Tapped Horn and Synergy Horn Technologies

Directivity Controllable Parametric Loudspeaker using Array Control System with High Speed 1-bit Signal Processing

A cellular automaton for urban traffic noise

RECOMMENDATION ITU-R M.1652 *

Multiple Audio Spots Design Based on Separating Emission of Carrier and Sideband Waves

Abnormal Compressor Noise Diagnosis Using Sound Quality Evaluation And Acoustic Array Method

UNIT Derive the fundamental equation for free space propagation?

FREQUENCY RESPONSE AND LATENCY OF MEMS MICROPHONES: THEORY AND PRACTICE

Abstract. Propagation tests for land-mobile radio service

Reliability of aircraft noise evaluation by measurement for comparison with prediction

IMPULSE RESPONSE MEASUREMENT WITH SINE SWEEPS AND AMPLITUDE MODULATION SCHEMES. Q. Meng, D. Sen, S. Wang and L. Hayes

Method of Determining Effect of Heat on Mortar by Using Aerial Ultrasonic Waves with Finite Amplitude

Proceedings of Meetings on Acoustics

Transcription:

Field experiment on ground-to-ground sound propagation from a directional source Toshikazu Takanashi 1 ; Shinichi Sakamoto ; Sakae Yokoyama 3 ; Hirokazu Ishii 4 1 INC Engineering Co., Ltd., Japan Institute of Industrial Science, The University of Tokyo, Japan 3 Kobayasi Institute of Physical Research, Japan 4 Japan Aerospace Exploration Agency, Japan ABSTRACT When predicting sound propagation in long distance, meteorological effects resulting from wind and temperature should be taken into consideration. Regarding this problem, many studies have been carried out by field measurements, numerical analyses and physical experiments. To predict and assess traffic noise by vehicles, trains and aircrafts, the sound sources are roughly modeled as omnidirectional. In reality, however, the sources have their inherent directivities according to their shapes, and the directivity can affects noise propagation characteristics as well as the meteorological effects. In this study, a field experiment on ground-to-ground long distance sound propagation using an omnidirectional loudspeaker and two types of directional loudspeakers were conducted at a flat field which approximately satisfied hemi-free field condition. To examine the relationship between the directivity effects and the meteorological effects, sound propagation characteristics by an omnidirectional loudspeaker and those by directional loudspeakers were compared. In addition, sound propagation by the omnidirectional point source was analyzed using Crank-Nicholson Parabolic Equation analyses and the experimental results were validated. Consequently, it was confirmed that sound propagation from the directional sources showed the same trend for excess attenuation characteristics as that from the omnidirectional source by considering its directional characteristics. Keywords: Meteorological effects, Directional source, long distance sound propagation, excess attenuation I-INCE Classification of Subjects Number(s): 4.6 1. INTRODUCTION The influence of meteorological effects resulting from wind and temperature is remarkable to long distance sound propagation. Regarding this problem, many studies have been carried out by field measurements, numerical analyses and physical experiments(1-4). To predict and assess traffic noise by vehicles, trains and aircrafts, the sound sources are roughly modeled as omnidirectional. In reality, however, the sources have their inherent directivities according to their shapes, and the directivity can affects noise propagation characteristics as well as the meteorological effects. Therefore it is necessary to grasp the influence of the sound source directivity to long distance sound propagation. In this study, we conducted the field experiment on ground-to-ground long distance sound propagation, in which an omnidirectional loudspeaker and two types of directional loudspeakers were used in the flat field which satisfied hemi-free field condition, and effects of meteorological condition and the source directivity were experimentally investigated. In addition, verification of an experimental result and parametric study on the meteorological effects were performed using a 1 t_takanashi@inc.ihi.co.jp sakamo@iis.u-tokyo.ac.jp 3 sakae@kobayasi-riken.or.jp 4 ishii.hirokazu@jaxa.jp Inter-noise 14 Page 1 of 9

Page of 9 Inter-noise 14 numerical analysis by the Crank-Nicholson Parabolic Equation (CN-PE) method (5).. Directional sound sources Figure 1 (a) and 1 (b) show directional loudspeakers used in this experiment. In advance of the outdoor field experiment, their directional characteristics were measured in an anechoic room. Horizontal cross section Vertical cross section (a) Directional loudspeaker 1 (b) Directional loudspeaker Figure 1 Directional sound source db Directional loudspeaker 1 Directional loudspeaker.1 1. 1. Frequency [khz] Figure Frequency characteristics of directional loudspeakers Figure shows frequency characteristics of sound pressure level measured at a point m away from in front of the respective loudspeakers. Frequency components ranging from 5 Hz to khz could be reproduced from the loudspeakers. Figure 3 shows directivity characteristics of sound pressure level of the two directional loudspeakers. It is shown that the directional loudspeaker 1 has sharper directivity characteristics than the directional loudspeaker. Since the directional loudspeaker to is not rotational symmetry, directivity characteristics in to plains perpendicular to each other were measured(see Fig.3 b and Fig.3 c), Comparing these results, the horizontal directivity is a bit sharper than the vertical one. 45 45 45 9 9 9 135 135 135 18 a). directional loudspeaker1 3. Field experiment 18 b). directional loudspeaker 18 c). directional loudspeaker (Horizontal cross section) (Vertical cross section) 5 Hz 5 Hz 1 khz khz 1grid : 5dB Figure 3 Directivity characteristics of the two directional loudspeakers 3.1 Outline Field experiment on outdoor sound propagation was carried out in early-mid July 13 at a Page of 9 Inter-noise 14

Inter-noise 14 Page 3 of 9 measurement field having a runway with a length of 1 km and a width of 6 m located at Taiki aerospace research field. Figure 4 shows the layout of the measurement field. The runway extends in the east and west direction. The sound sources were set at the center position of the runway, and 11 receiving points were arranged on the centerline of the runway with equally a 1 m intervals. In order to grasp meteorological condition, two-dimensional ultrasonic anemometers were set at 3 points of west-3 m point, east-3 m point and center point as shown Fig.4. Figure 5 shows setting configurations of sound sources, microphones and anemometers. The microphones and anemometers were set at their heights of 1. m as shown in Fig. 5. V N U Definition of wind components west-3 m center east-3 m 1 km :Runway :Grass field :Windbreak forest :Sound source point :Sound receiving point :Ultrasonic anemometer Figure 4 Layout of the measurement field Directional loudspeaker Omni-directional loudspeaker. m and 1. m 1.4 m 1. m Sound level meter Ultrasonic anemometer West 5 m West 3 m Source point East 3 m East 5 m Figure 5 Setting configurations of equipments Figure 6 Sound receiving point Figure 7 Experiment scenery for source point In this experiment, three loudspeakers (an omnidirectional loudspeaker and two directional loudspeakers) were used. The omnidirectional loudspeaker was set at a point of the height of 1.4 m and two directional loudspeakers were set at two points of their heights of m and 1 m using an extensible pole, as shown in Figs. 6 and 7. The directional loudspeakers were pointed to the east and west along the runway. As the source signal, swept-sine signals were used in order to measure Inter-noise 14 Page 3 of 9

Page 4 of 9 Inter-noise 14 impulse responses from the source points to the receiving points. In order to secure sufficient signal-to-noise (S/N) ratio, time duration of the swept-sine signal was set to 6 sec as sufficiently long duration to obtain enough sound energy. The frequency components included in the source signals ranged 5 octave bands from 15 Hz to khz for the omnidirectional loudspeaker and 4 octave bands from 5 Hz to khz for the directional loudspeakers. 4. Experimental results 4.1 Characteristics of distance attenuation Figures 8 (a) and 8 (b) show characteristics of distance attenuation of sound pressure level when the omnidirectional loudspeaker was set at the height of 1.4 m(fig. 8 (a)) and the directional loudspeaker 1 was set at the height of 1 m (Fig. 8 (b) ). 163 data for the omnidirectional loudspeaker and 54 data for the directional loudspeaker 1 measured for 1 days were overdrawn. In the case of the directional loudspeaker 1, the direction of the loudspeaker was set in the east. In these figures, all of the measurement data is overwritten regardless of the wind conditions. From all figures shown here, it is clear that the variation of the sound pressure level became larger as the distance became larger. When the distance was 5 m, the variation extended to db. For the omnidirectional loudspeaker, the distance attenuation in the western side tended slightly gentler than that in the eastern side. This is because wind direction on the measurement period tended to be easterly and as the result the measurement points tended to be at the downwind side. Regarding the directional loudspeaker 1, the sound pressure level at the east side was obviously larger than that at the west side because the loudspeaker pointed to the east. In order to eliminate the effect of the directivity characteristics, the direction angle of the loudspeaker at the respective measurement points were calculated and the effects of the directivity characteristics were corrected from the measurement data based on the directivity characteristics shown in Fig. 3 (a), 3 (b) and 3 (c). The calculation results are shown in Fig. 9. The tendency of the variation of the sound pressure level became similar as those for the omnidirectional loudspeaker. db West side East side -6 db/d.d. 5 Hz db 1 khz 1. 1. 1. 1. 1. 1. 1. 1. (a) Omni-directional loudspeaker db 5 Hz 1 khz db 1. 1. 1. 1. 1. 1. 1. 1. (b) Directional loudspeaker 1 (height 1 m and east direction) Figure 8 Distance attenuation of sound pressure level Page 4 of 9 Inter-noise 14

Inter-noise 14 Page 5 of 9 West side East side -6 db/d.d. db 5 Hz 1 khz db 1. 1. 1. 1. 1. 1. 1. 1. Directional loudspeaker 1 (height 1 m and east direction) Figure 9 Corrected distance attenuation regarding the source directivity 4. Excess attenuation Excess attenuation over a standard condition, in which no wind or no strong temperature profile is assumed, was calculated for sound propagation from the three loudspeakers. As a reference sound pressure level, only distance attenuation on a rigid surface was taken into consideration. Under a geometrical condition shown in Fig. 1, sound pressure level at a receiving point P is calculated as a summation of contributions for direct path and reflecting path in energy base as, LDirect LReflect 1 1 L = Dist 1log1 1 + 1, (1) where, L Direct and L Reflect are sound pressure levels for direct and reflect paths, respectively, and they are calculated as follows. L Direct = L R ( θ ) d 1 log 1 + L, () dir d R Rr LReflect = L + Ldir ( θ r ) R 1 log 1, (3) where, L denotes the sound pressure level at a reference point R m distant from the sound source, L dir (θ d ) and L dir (θ r ) denote the corrected values of the correction [db] due to the source directivity at the angle θ d and θ r, respectively. For the omnidirectional loudspeaker, L dir (θ d ) and L dir (θ r ) are db for all angles, and for the directional loudspeakers 1 and, measured results of the directivities of the loudspeakers described in the chapter were used as the values of L dir (θ d ) and L dir (θ r ). The excess attenuations were of obtained as differences between measures values and the reference sound pressure levels. S R d θ d P h s h r h r R r S' θ r l Figure 1 Geometrical configuration of the source, receive and a flat, rigid surface Excess attenuation levels for all measurement data were calculated as the procedure mentioned above and the results were arranged in relationship with the vector component of the wind speed at 1. m high, U. U is positive in the downwind direction and negative in the upwind direction. The calculation results are shown in Figs. 11, 1 and 13. The figures show that the excess attenuation is Inter-noise 14 Page 5 of 9

Page 6 of 9 Inter-noise 14 related to the vector wind speed, distance from the source, frequency of sound and height of the source. The excess attenuation becomes larger as the absolute value of the vector wind speed becomes larger in negative, as the distance from the source and receiving point becomes larger, and as the frequency is higher. 3 1-1 -8-6 -4-4 6 8 3 1 3 1-1 -8-6 -4-4 6 8 3 1 3 1 m Point 3 m Point 5 m Point 1 m Point 3 3 m Point 5 m Point 1-1 -8-6 -4-4 6 8 Omni-directional (height=1.4 m) Directional 1 (height=.m) Directional (height=.m) -1-8 -6-4 - 4 6 8 Figure 1 Excess attenuation 1 khz ( directional speaker height=.m ) 1-1 -8-6 -4-4 6 8 Omni-directional (height=1.4 m) Directional 1 (height=.m) Directional (height=.m) -1-8 -6-4 - 4 6 8 Figure 11 Excess attenuation 5 Hz ( directional speaker height=.m ) Page 6 of 9 Inter-noise 14

3 1 m Point 1-1 Page 7 of 9-8 -6-4 - 4 6 8 Inter-noise 14 3 3 m Point 1-1 3-8 -6-4 - 4 6 8 Omni-directional (height=1.4 m) 5 m Point 1 Directional 1 (height=1.m) Directional (height=1.m) -1-8 -6-4 - 4 6 8 Figure 13 Excess attenuation 5 Hz ( directional loudspeaker height=1.m ) 5. Comparison between measurement and PE analysis Crank-Nicholson Parabolic Equation (CN-PE) analyses were conducted for the same geometrical configuration of the source and receivers as the experiment for the omnidirectional loudspeaker, of which height is 1.4 m. The calculation results were compared with the experimental results. In the calculation, meteorological conditions of wind speed U and temperature T were set as, z U (z ) = a log 1 +, z T ( z ) = T + bz, (4) (5) where, z is the height [m], z is the roughness length [m], T is the temperature in Celsius degree [ C] on the ground and b is a rise rate of the temperature. The parameters z, T and b were set to.3 m, 17 C and.3, respectively. Then, sound speed at the height z, c(z), is expressed as follows. 1/ T (z ) c( z ) = 331.5 1 + 73 + U (z ), (6) Under the above conditions, sound pressure distribution for single frequencies of 5 Hz, 5 Hz and 1 khz was calculated. Figure 14 shows an example of calculation result for 5 Hz, which shows an influence of wind speed on the sound pressure distribution. 1 db Height [m] Height [m] 1 (a) U at 1. m is m/s 5 (b) U at 1. m is 4 m/s 5-6 db Figure 14 Calculation result of CN-PE analysis at 5 Hz Inter-noise 14 Page 7 of 9

Page 8 of 9 Inter-noise 14 The excess attenuation for an omnidirectional source was calculated as a level difference between a case where the vertical distribution of wind speed, temperature and sound speed was set parametrically and a case where the vertical distribution of the sound speed was uniform. Figure 15 shows calculated excess attenuation overdrawn on the measurement results. Here, it should be noted that the measurement results were analyzed as 1/1 octave band values, whereas the calculation results were for single frequency. Calculation for single frequency may emphasize an influence of sound interference as sharp peaks or dips in the graph. Comparing the calculation and measurement for the omnidirectional source, similar tendency is seen in the results for 1 m point, whereas the difference between them is considerable in the results for 3 m. In the figure, calculated excess attenuation for an omnidirectional source located at 1 m high is shown as references. The calculated excess attenuation characteristics in relation to vector wind speed are similar with the measured ones for the directional loudspeakers 1 and. Omni-directional (height: 1.4 m) 3 1 3 1 m Point 3 m Point -1-1 -8-6 -4-4 6 8-8 -6-4 - 4 6 8 a). 1/1 octave band center frequency 5 Hz 3 1 Directional 1 (height: 1. m) Numerical analyses (height: 1.4 m) Numerical analyses (height: 1. m) -1-1 -8-6 -4-4 6 8-8 -6-4 - 4 6 8 b). 1/1 octave band center frequency 1 khz Figure 15 Comparison of the excess attenuation between the field experiment and numerical analyses 1 1 m Point 3 3 m Point 1 Directional (height: 1. m) 6. CONCLUSIONS In this study, long distance sound propagation from several types of sound sources, an omnidirectional loudspeaker and directional ones, were investigated experimentally to examine the relationship between the directivity effects and the meteorological effects. In addition, sound propagation by the omnidirectional point source was analyzed using the Crank-Nicholson Parabolic Equation analysis and the measurement results were validated. Consequently, it was confirmed that sound propagation from the directional sources showed the same trend for excess attenuation characteristics as that from an omnidirectional source. Page 8 of 9 Inter-noise 14

Inter-noise 14 Page 9 of 9 REFERENCES 1. T Yokota, K Makino, K Yamamoto, Y Okada and K Yoshihisa. Field measurements in Multi-purpose Aerospace Park at Taiki, Hokkaido : Comparison between numerical analysis by PE method and measurement on outdoor sound propagation. Proceedings of the spring meeting, the institute of Noise Control Engineering of Japan 6, 17-13, April 6.. K Yoshihisa, Y Okada and Y Itou. Field measurements in Multi-purpose Aerospace Park at Taiki, Hokkaido : Study on the effects of wind turbulence on complex sound propagation. Proceedings of the Autumn meeting, the institute of Noise Control Engineering of Japan 6, 13-16, September 6. 3. T Yokoi and T Oshima. Investigations of wind effects on road traffic noise propagation from a road embankment : Examination of the improvement proposal of ASJ RTN-Model. Summaries of Technical Papers of Annual Meeting Architectural Institute of Japan, 49-5, August 13. 4. T Yokota, K Makino and K Yamamoto. Finite difference time domain simulation of outdoor sound propagation under the influences of wind speed gradient. Proceedings of the autumn meeting, the Institute of Noise Control Engineering of Japan 1, 59-6, September 1. 5. M. West, K. Gilbert, R.A. Sack. A tutorial on the parabolic equation (PE) model used for long range sound propagation in the atmosphere. Applied Acoustics, Vol.37, pp.31-49, 199. Inter-noise 14 Page 9 of 9