Key Comparison CCAUV.A-K4

Final Report on the Key Comparison CCAUV.A-K4 Abstract The Key Comparison CCAUV.A-K4 has been carried out under the auspices of the Consultative Committee for Acoustics, Ultrasound and Vibration (CCAUV) of the International Bureau of Weights and Measures (BIPM). The participating NMI s are the Instituto Nacional de Metrologia, Normalizacao e Qualidade Industrial (INMETRO, Brazil), the Centro Nacional de Metrología (CENAM, Mexico), the Laboratoire National de Métrologie et d'essais (LNE, France), the National Metrology Institute of Japan (NMIJ, Japan), the Korea Research Institute of Standards and Science (KRISS, Korea), and the Physikalisch-Technische Bundesanstalt (PTB, Germany). The role of the Pilot laboratory was undertaken by the DFM and CENAM. The measurements took place between January 007 and February 008. The time schedule was organized in a single star configuration. Two LSaP microphones were circulated. This report includes the measurement results from the participants, information about their calibration methods, and the analysis leading to the assignation of the degrees of equivalence. Prepared by Salvador Barrera-Figueroa, Lars Nielsen and Knud Rasmussen Danish Primary Laboratory of Acoustics Danish Fundamental Metrology Ltd. Matematiktorvet 307 DK-800 Kgs. Lyngby, Denmark Andrés E. Pérez Matzumoto and José Noé Razo Razo Centro Nacional de Metrología Carretera a Los Cués km 4.5 El Marqués, Querétaro, México. CP 7646

Content 1 Introduction 3 Comparison protocol 3 3 Travelling microphones 4 4.1. Calibration methods 5 4 Calibration methods and reported results 5 4.. Microphone parameters 13 4.3. Microphone sensitivities 14 5 Analysis of the results 17 5.1. Degrees of equivalence 17 5.. Equivalence and consistence of results 18 6 Results 0 6.1. Reference values 0 6. Degrees of equivalence per laboratory 0 6.3 Degrees of equivalence 5 7 Conclusions 7 8 References 7 Appendix A Determination of the correlation coefficients 8 Appendix B INMETRO corrections to reported results 9 Appendix C Uncertainty budgets 30

1 Introduction The CCAUV.A-K4 is a Key comparison realised under the auspices of the CCAUV-BIPM. The standards circulated among the laboratories were two LSaP microphones of the type Brüel & Kjær 4180. Both microphones were also standards used for the CCAUV.A-K3. The microphones had to be calibrated using the reciprocity technique under free-field conditions in the frequency range from 1 khz to 31.5 khz (optionally up to 40 khz), and the open-circuit free-field sensitivity level had to be reported at the nominal preferred 1/3 rd octave band. The participating NMIs are the Instituto Nacional de Metrologia, Normalizacao e Qualidade Industrial (INMETRO, Brazil), the Centro Nacional de Metrología (CENAM, Mexico), the Laboratoire National de Métrologie et d'essais (LNE, France), the National Metrology Institute of Japan (NMIJ, Japan), the Korea Research Institute of Standards and Science (KRISS, Korea), and the Physikalisch-Technische Bundesanstalt (PTB, Germany). The National Physical Laboratory (NPL, UK) and the National Institute of Standards and Technology (NIST, USA) were included in the original time schedule, however they resigned participating during the course of the comparison. The time schedule was organized in a single star configuration. The role of the Pilot laboratory was undertaken by the DFM and CENAM. The measurements took place between February 007 and February 008. This report includes the measurement results from the participants, information about their calibration methods, uncertainty budgets, and the analysis leading to the assignation of the degrees of equivalence. Comparison protocol The September 006 Draft technical protocol for the key comparison CCAUV.A-K4 was presented during the 5th meeting of CCAUV on September 006, and the subsequently prepared October 006 Draft was used for this key comparison. This protocol is based on the CCAUV.A-K3 protocol with some additional clarifications. It was emphasized that the calibrations had to be performed at the preferred nominal frequencies rather than the exact frequencies. It was also recommended that the laboratories should report the results using three decimals. In this report, results are presented in the same way they were declared. The microphones were transported in all cases by means of a courier company selected by each participant. The container used under the courier transportation was provided by DFM. The actual participants are listed in table 1: Table 1. Participants in the CCAUV.A-K4 key comparison Participant Acronym Country Country code Regional Metrology Organization Danish Primary Laboratory for Acoustics. DFM Denmark DK EURAMET Danish Fundamental Metrology LTD. Centro Nacional de Metrología CENAM Mexico MX SIM Korea Research Institute of Standards KRISS Korea KR APMP and Science National Metrology Institute of Japan NMIJ Japan JP APMP Laboratoire national de métrologie et LNE France FR EURAMET d'essais Physikalisch-Technische Bundesanstalt PTB Germany DE COOMET (EURAMET) Instituto Nacional de Metrologia, Normalizacão e Qualidade Industrial INMETRO Brazil BR SIM

3 Travelling microphones Two LS microphones, Brüel & Kjær 4180, serial numbers 1395455 and 1395456, were supplied by DFM and circulated among the participants. The DFM has monitored the microphones since the beginning of the comparison. The monitoring consisted of free-field reciprocity calibrations and also pressure reciprocity calibrations. The DFM also performed free-field reciprocity measurements from 1 khz as stated in the protocol. Figure 1 shows DFM s results of the control calibrations for some selected frequencies during the comparison span. Figure 1 also shows the difference with respect to a common average for the two microphones over the whole frequency range. No significant trend was observed in the control calibrations. Sensitivity difference (db) Normalised sensitivity (db) 0. 0.1 0-0.1 Microphone 4180.1395455-0. 10 0 10 1 Frequency (khz) Microphone 4180.1395455 0. 0.1 0-0.1-0. 0/07 06/07 09/07 1/07 04/08 date Sensitivity difference (db) Normalised sensitivity (db) 0. 0.1 0-0.1 Microphone 4180.1395456-0. 10 0 10 1 Frequency (khz) Microphone 4180.1395456 0. 0.1 0-0.1-0. 0/07 06/07 09/07 1/07 04/08 date Figure 1. Upper charts: Difference of the sensitivity level obtained at each calibration at DFM with respect to their common average for each microphone as a function of frequency (DFM s uncertainty bounds are shown). Lower charts: Changes in the sensitivity level of the microphones as a function of time. The sensitivity level has been normalized to the level of the first measurement. Blue line: 1 khz; red line khz; black line 5 khz; and magenta line 10 khz.

4 Calibration methods and reported results 4.1. Calibration methods A description of the calibration methods and the reporting of results of each laboratory is given below as submitted by each laboratory. No editing of the contents but a small amount of reformatting has been performed. Consequently, table, equation and figure numbers and references do not follow the numbering of the report. DFM: Free-field calibration of LS microphones. The free-field sensitivity of LS microphones are determined by conducting a reciprocity calibration in accordance with IEC 61094-3 [1]. The calibrations were performed in a small, specially designed anechoic room of about 1.4 m 1.50 m 1.7 m free dimensions. The reflection coefficient of the lining material is below -40 db above 1 khz. Measurement method The electrical transfer impedance is measured in linear frequency steps of 10 Hz from 3 khz (initially 960 Hz) to 51 khz. Due to the unavoidable poor signal-to-noise ratio and reflections from the environment a special technique has been developed based on narrow-band filtering in both frequency- and time domain [] Filtering in the frequency domain is applied by using a B&K Pulse Analyzer in the SSR-mode when measuring the transfer impedance. To enter the time domain by applying an inverse FFT-technique to the measurement result, the missing low frequency values down to 0 Hz are estimated as follows: The average pressure response (normalized at 50 Hz) has been determined for a large number of BK type 4180 microphones. The free-field to pressure sensitivity diffraction has been calculated as well using the Boundary Element Method [3]. The missing values of the transfer impedance from 10 Hz to 3 khz (resp. 1kHz) are then estimated from the above normalized responses, the absolute pressure sensitivity of the involved microphones at the actual environmental conditions and assuming an acoustic centre of 4.74 mm. To avoid sudden jumps in the transfer impedance the calculated values are blended into the measurement results over a frequency range of 100 Hz. Above 51 khz the measurements are extended to 70 khz by a linear extrapolation of the measurement results from 48 51 khz in order to avoid an abrupt high-frequency cut-off when applying the inverse FFT-technique. Next, in the time domain, a time window of 0.9 ms is applied to eliminate the unwanted part of the signal and the remaining signal is converted back to the frequency domain. The final free-field sensitivities are then calculated in accordance with IEC 61094-3 using the above manipulated transfer impedances.

Acoustic centres The acoustic centres used are derived from a linear regression analysis of the above manipulated transfer impedances at four distances (corrected for air attenuation). This has been done for a large number of microphones and the average value has been expressed as a polynomial slightly different from the values in reference [4]. Correction to reference environmental conditions In general, calibrations will be performed at the following environmental conditions: + a static pressure of 101.35 ± 3 kpa + a temperature of 3 ± 1 C and + a relative humidity of 50 ± 15 %RH The calibration results are corrected to reference environmental conditions by applying the individual static pressure- and temperature coefficients valid for the pressure sensitivity, i.e., the additional influence of the radiation impedance has been neglected. Correction to nominal frequencies The correction of the measured sensitivities in steps of 10 Hz to the required nominal frequencies has been performed by a simple linear interpolation between the measured sensitivities. References [1] IEC 61094-3, 1995. Measurement microphones Part 3: Primary method for free-field calibration of laboratory standard microphones by the reciprocity technique. [] S. Barrera Figueroa, New methods for transducer calibration: Free-field reciprocity calibration of condenser microphones, PhD thesis, Ørsted DTU, Technical University of Denmark, 003. [3] S. Barrera Figueroa, K. Rasmussen, and F. Jacobsen, On experimental determination of the free-field correction of laboratory standard microphones at normal incidence, Metrologia 44, pp 57-63 (007). [4] S. Barrera Figueroa, K. Rasmussen, and F. Jacobsen, The acoustic center of laboratory standard microphones, Journal of the Acoustical Society of America 10, pp 668-675 (006). CENAM: a) Procedure: 510-AC-P.105: Calibración primaria de micrófonos patrón mediante la técnica de reciprocidad en campo libre (Primary calibration of standard microphones by the reciprocity technique in free-field). b) Instrumentation: + Microphone B&K Type 4180, SN: 167796 + Microphone B&K Type 4180, SN: 1893458 + Reciprocity apparatus B&K Type 5998, SN: 04046. + Polarization voltage meter B&K Type 5981, SN: 1468333. + Signal Analyzer B&K Type 3560D/3110, SN: 415164. + Measuring Amplifier B&K Type 610, SN: 919975

+ Calibrator Fluke Type 5700A, SN: 630304. + Environmental Monitoring System CENAM SN 063. c) Measurement: c.1) Determination of the free-field sensitivity modulus The open circuit sensitivity is determined by means of the reciprocity technique according to IEC 61094-3, [1], using a set of three microphones of the same type, {1,, 3}. The microphone under test is sequentially acoustically coupled in pairs, {1-, 1-3 and -3}, at known distances within the two microphones in a field essentially free of acoustic reflections (free-field). The sensitivity of the microphone is obtained from the measurements of the electrical transfer impedance function, and the calculation of the acoustical transfer impedance between the three pairs of microphones. Electrical transfer impedance is measured in the frequency range from 1 15 Hz to 5 99 Hz with a 18 Hz resolution, using a dual channel digital signal analyzer. A cleaning procedure, based on time selective window, is applied to reduce acoustic reflections, microphone interference and cross-talk between the measurement channels. Sensitivity levels for 1/3 octave frequencies above 3 khz are calculated by linear interpolation from the set of measurement results in the frequency range given above, according to: M f, f M high M f, freq = M f, f + ) (1) low f f high f, flow ( freq flow low where: freq : frequency of interest in [Hz]; f low : frequency immediately below freq, in [Hz] in the set of measurement result; f high : frequency immediately above freq, in [Hz] in the set of measurement result; M f,freq : sensitivity level corresponding to the frequency freq, in [db re 1V/Pa]; M f,flow : sensitivity level corresponding to the frequency f low, in [db re 1V/Pa] M f,fhigh : sensitivity level corresponding to the frequency f high, in [db re 1V/Pa] Free-field sensitivity levels for frequencies below 3 khz are calculated through a free-field correction, determined from measured free-field and pressure sensitivity levels. A ninth degree polynomial fitting is calculated for each free-field correction, obtained from repeated measurements. For frequencies below 500 Hz it is assumed that both sensitivity levels, free field and pressure, are equal. Then the free-field correction calculated, through the polynomial fitting, for frequencies from 1 khz to 3 khz is used to determine the corresponding free-field sensitivity level. Four different separation distances, between each pair of microphones, are used in each measurement; these are given in table A: Table A: Separation distances used during the electrical transfer impedance measurements, for each pair of microphones. Number 1 3 4 Distance / mm 160.03 189.83 0.03 49.40 Corrections for the position of the acoustic centre are applied, and the corresponding correction factors for static pressure and temperature are used to obtain the sensitivity at reference conditions (101.35 kpa and 96.15 K (3 ºC)). Polarizing voltage of 00 V ± 0.04 V is applied to both, transmitter and receiving, measuring channels.

c.) Acoustic centre and environmental coefficients Positions for acoustic centre are obtained from reference []. Static pressure coefficients are determined for each microphone using the pressure sensitivity levels obtained by the reciprocity technique for five static pressures from 81 kpa to 101.4 kpa. Later a linear regression fitting is used to obtain the corresponding coefficient for frequencies from 15 Hz up to 31.5 khz. Temperature coefficient is calculated using the typical value given by the manufacturer and the polynomial approximation given in [3]. The Environmental conditions of measurement are shown in table B Table B: Environmental conditions of measurement for each standard Microphones 4180.1395455 4180.1395456 Static Pressure: 80.97 ± 0.35 kpa 81.00 ± 0.31 kpa Temperature:.47 ± 1.89 C.0 ± 1.44 C Relative Humidity: 3.50 ± 9.55 % 36.30 ± 7.09 % d) Results a e) Uncertainty: The expanded measurement uncertainty is obtained multiplying the standard combined uncertainty by a coverage factor, k =, which gives a confidence level of al least 95 %, under the assumption that the probability distribution function for the measurand is approximately normal. The measurement uncertainty is estimated according to the Mexican National Standard NMX-CH- 140-IMNC-00: Guía para la Expresión de la Incertidumbre en las Mediciones, which is equivalent to the Guide for the Expression of Uncertainties in Measurement: BIPM, IEC, IFCC, ISO, IU- PAC, IUPAP, OIML (1995). The uncertainty value given in this certificate, does not take into account the contributions due to long term stability, drift and transportation for the calibrated instrument. f) Traceability: Calibration results are traceable to the kilogram (kg), Kelvin (K), second (s), metre (m) and ampere (A), base units of the SI, through the Mexican National Standards for Barometric Pressure (CNM-PNM-4), Temperature (CNM-PNE-), Time (CNM-PNE-1), Length (CNM-PNM-) and Intensity of Alternating Current (CNM-PNE-10), maintained by CENAM. References: [1] IEC International Standard 61094/3. Measurement microphones Part 3: Primary method for freefield calibration of laboratory standard microphones by the reciprocity technique. a Note: CENAM results are presented in the next sub-section.

[] Salvador Barrera-Figueroa, Knud Rasmussen and Finn Jacobsen, The acoustic center of laboratory standard microphones, J. Acoust. Soc. Am. 10 (5), 668-675 (006). [3] Brüel & Kjær. Technical Review No. 1-001. Properties and calibration of laboratory standard microphones & Uncertainties in microphone frequency response. ISSN 0007-61. 001. KRISS: The dimensions of the anechoic chamber used for the free-field calibration are.3 m (depth),. m (width), and.4 m (height). The microphones are mounted on 1.7 mm diameter vertical rods. The distance between two microphones is adjusted by a computer controlled servo motor and encoder with the resolution of 1 mm. The voltage ratio functions are measured using Reciprocity Calibration Apparatus (B&K 5998) and PULSE (B&K 3560C). The frequency is swept from 0.8 khz to 51 khz in linear frequency steps of 10 Hz. The number of measurement at each frequency is increased until the statistical uncertainty of the measurement response (i.e., twice the standard deviation) is less than 0.01 db. A perfect free-field facility of no interference from walls and obstacles is not readily available in practice, so it has led to the use of signal processing techniques to eliminate the reflections from boundaries. In addition, an electrical interference noise due to imperfect insulation of measuring instruments, so-called cross-talk, should be eliminated with the same manner. The time selective windowing technique for effective removal of multi-path noises is widely utilized for reciprocity calibrations, but it is still short in overcoming the limitations imposed by overlapped signals, in particular at low frequencies. Based on FFT analysis, the leakage due to the limited frequency range makes the overlap problem worse, not to be perfectly separated. Instead of conventional low pass filter designed to have a flat response in the frequency range of interest, Chebyshev window shaped function is selected because of its low sidelobe levels. After multi-path noise removal, the voltage ratio functions are restored by compensation with inversing filtering. The free-field sensitivity level at the nominal preferred 1/3rd octave centre frequencies is calculated with the linear interpolation between those at two adjacent frequencies. The free-field sensitivity levels near 1 khz are estimated by the interpolation technique between the pressure sensitivity levels at the low frequencies and the free-field sensitivities at the high frequencies. The distances between two microphones are varied from 100 mm to 300 mm with 50 mm step. The final results are the average of these five measurements. The acoustic centres used for the sensitivity calculation are the nominal values given in IEC 61094-3:1995. The acoustic centres for the higher frequencies were obtained by extrapolation method. NMIJ: The free-field sensitivity was determined in compliance with IEC 61094-3:1995, using a reciprocity calibration system developed by NMIJ. In the system, both the signal generation and the signal processing were executed by a dual-channel FFT analyzer, model CF-50 of ON0 SOKKI Co. Signal to noise ratio was improved by the synchronous waveform averaging method. The insert voltage technique was used to cancel the effect of the gain and impedance of an electrical circuit. The calibration was performed by using software of our own making. Sound absorbing wedges made of glasswool are attached on the inside walls, ceiling and floor of an anechoic chamber. The inner space of the anechoic chamber is (W) 9.5 m by (L) 8.0 m by (H) 7. m. The cutoff frequency of the anechoic chamber is 40 Hz. The temperature and humidity inside the anechoic chamber are well controlled (3.0 ± 0.5 ºC and 50 ± 5 %, respectively).

The sensitivity was corrected to the reference environmental conditions by using K. Rasmussen's method [I]. At 31.5 khz, only pressure dependency has been corrected because of no reliable temperature coefficients. At 40.0 khz, neither pressure nor temperature dependency has been corrected because of no reliable pressure and temperature coefficients. The calibrations were normally carried out using only B&K type 4180 microphones. At low frequencies below 5000 Hz, however, type WS microphones (B&K type 4165 and type 4190) with cavity rings (B&K UA085) were also used to achieve high signal to noise ratio because the sensitivity of these microphones was about 10 db higher than that of type 4180 microphones. B&K type 4165 microphones were just used as the transmitter microphones and B&K type 4 190 microphones were used as the reciprocal microphones. To minimize the influence of sound reflection, the virtual pulse technique was applied [, 3]. In this technique, firstly the voltage transfer function (VTF) between the input terminal of the transmitter microphone and the output terminal of the receiver microphone was measured by using swept-sine signal. The bandwidth of the swept-sine signal was 00 Hz (1 khz - 1.6 khz), 400 Hz (.0 khz - 5 khz), khz (6.3 khz - 1.5 khz) and 4 khz (16 khz - 40 khz), respectively. The centre frequency of swept-sine signal was equal to each measurement frequency. Secondly, the pulse response of VTF was computationally calculated by multiplying the virtual pulse signal by the VTF in the frequency domain. The influence of sound reflection could be minimized by the convolution of the pulse response and a window function [, 3]. The acoustic centres were calculated from the VTFs by using least square method. In the VTFs measurements, the actual distances were set to 10, 15 and 0 cm (1000 Hz - 5000 Hz) and 15, 0 and 5 cm (6300 Hz - 40000 Hz), respectively. The microphones used for determining the acoustic centre were the same as the calibration. The influence of reflection was similarly minimized by the virtual pulse method [, 3]. Values of the air attenuation coefficients were obtained according to IEC 61093-3 Annex B. Density of air was calculated in compliance with IEC 61094-, Annex F.1 (Density of humid air). Reference [1] K. Rasmussen, "The static pressure and temperature coefficients of laboratory standard microphones," Metrologia 36 (1999). [] T.Fujimori, "Free-field calibration of "1/-in" microphone," 3 rd. joint meeting of ASA and ASJ (1996). [3] T. Fujimori et. al, "Measurement of acoustic center position of LSP Laboratory standard microphones," J. Acoust. Soc. Jpn. ( J ) 58, 579-585 (00). LNE: According to the technical protocol of the key comparison CCAUV.A-K4 (October 006), this report completes the two standard calibration certificates issued from LNE after the key comparison CCAUV.A-K4 in order to provide additional information on the calibration method. All the pieces of information required in the technical protocol ( reporting results) are listed below in the following paragraphs.

The environmental conditions during the calibrations. Pressure : 1000.0 to 1004.5 hpa Temperature :.8 to 3. C Relative humidity : 44 to 50 % The actual distance(s) used between the microphones during the calibrations. The results given in the calibration certificate correspond to the mean values of three calibrations obtained using three distances. The used distances between microphones body (not including frontal cavity depth) are: 5 cm, 30 cm and 35 cm. Values of the temperature and pressure coefficients of the microphones used in the calculations. The temperature and pressure coefficients of the microphones are calculated using the polynomial constants given by the manufacturer (technical review N 1 001) and using a constant microphone resonant frequency of 3 khz. A summary of the uncertainty calculation and the overall uncertainties calculated for a coverage factor of k=. The budget was sent before the comparison starts. The effect on the reported results. Air density is calculated using the polynomial formula given in Report PL11b, 1997 DTU. The effect on sensitivity calculation seems to be lower than 0.0015 db. Air attenuation coefficient is calculated using the polynomial formula given in Report PL11b, 1997 (DTU). The effect on sensitivity calculation seems to be lower than 0.0005 db. The values of the acoustic centres used in the calculation. Acoustic centres used in calculation are those given in IEC 61094-3:1995. For frequencies greater than 0 khz, this value is set to 0 mm. Calibration results using three decimals (if not already given in the certificate). Calibration results are given with 3 decimals in the calibration certificate and with a full resolution in the additional results spreadsheet. The interpolation technique used if the actual measurements are not performed at the required frequencies. None interpolation technique is used: calibrations are carried out at the required frequencies. The technique used to estimate the free-field sensitivity at low frequencies if a direct calibration is not available. Free field calibrations are carried out at the all the nominal preferred 1/3-octave frequencies from 1000 to 31500 Hz. However in order to provide acoustical disturbances filtering, calibration requires additional frequency points which are measured for frequencies greater than 950 Hz and estimated above. That estimation involves the individual pressure sensitivity of each microphone and the diffraction correction given in IEC/TS 61094-7.

The effect of the diffraction correction is lower than 0.01 db on the free-field sensitivity at 1 khz and negligible for the other frequencies. PTB: Methodology for free field calibration CCAUV-A-K4 The calibration was performed as a reciprocity calibration according to IEC 61094-3 in an electrically shielded anechoic chamber with a net volume of about 15 m 3 and free dimensions of about.0 x.0 x 3.7 m 3. The microphone mounts were vertically placed half-inch hollow tubes holding the receiver preamplifier (B&K 645) and the transmitter housing. An essentially constant voltage of approximately 6 V (rms) was directly fed to the terminals of the transmitter housing and measured by one channel of a dual channel FFT analyzer type Ono Sokki 6400 which also provided the transmitter voltage and performed a phase-locked averaging of the input signals. The other channel of the analyzer was connected to the receiver chain, consisting of the preamplifier type B&K 645, a measuring amplifier type B&K 636 and a third-octave filter type B&K 1617. The gain of the receiver chain, including the loading of the receiver microphone by the preamplifier, was determined using the insert voltage technique. The analyzer measured the ratio of the FFT spectra of both channels. In order to calculate the transmitter current, the effective capacitance of the transmitter microphones was measured for all frequencies by means of a capacitance bridge type QuadTech 7600 using a specific microphone housing. The electrical transfer impedance was then calculated from the spectral ratio, the receiver channel gain and the effective transmitter microphone capacitance. The polarization voltage for the transmitter microphone was provided by a modified microphone front-end type Norsonic 336, for the receiver microphone it was generated within the B&K 636 measuring amplifier. The polarization voltages were checked by a differential voltmeter type Fluke 893A. For each measurement frequency several electrical transfer impedance measurements were performed in linear distance steps, distance ranges and number of steps are given in the certificates. The electrical transfer impedance per frequency was finally derived from a linear regression analysis of the electrical transfer impedance over the microphone distance. The open circuit free-field sensitivity was calculated using the formulae and data given in the current edition of IEC 61094-3 and transformed to the reference conditions, see certificate. Because all calibrations were performed at the required frequencies, no interpolation technique was applied. INMETRO: The basic procedure to calibrate microphones in free-field is described in the international standard IEC 61094-3. The calibration method used by INMETRO employs digital signal processing techniques to identify and to minimize the effect of the components that interfere in the determination of the sensitivity of the microphones. Among these components, crosstalk, multiple backscattering and sound reflections inside the anechoic chamber are considered. The crosstalk was suppressed using the subtraction technique. A correction file including the magnitude and phase response of a dummy microphone with similar capacitance to the B&K4180 is first

obtained. Then, all voltage measurements of the receiver microphone are corrected by the subtraction of this correction file. This procedure allows the suppression of this specific component, which interferes most strongly in the result of the reciprocity method in free field. The multiple backscattering between the microphones and the other sound reflections present in the impulse response of the voltage output of the receiver microphone were suppressed by using a time-selective signal processing technique. A hybrid non-symmetric time-window function formed by the combination hanning-rectangular-hanning with width of respectively ms, 11 ms, 0.85 ms was used for the time-selection of the direct sound. The measurements of the electrical impedance were carried out at the distances between the microphones of 170 mm, 00 mm, 40 mm and 300 mm. The analyzer was configured for a resolution of 1.47 Hz between frequency bins. The interpolation technique used to obtain the sensitivities at the preferred nominal frequencies was the simple arithmetic mean between the results measured at the adjacent lower and upper frequency bins. Table 1 presents the values of the used acoustic centres for the determination of the sensitivity of the B&K4180 microphones serial numbers 1395455 and 1395456. Table 1 Acoustic centres in metres. a The measured sensitivities at the environmental conditions of the laboratory were corrected to the reference environmental conditions. The correction method used is the one described in IEC 61094-3 (1995) for measurements up to 31.5 khz. This same correction method was also applied up to the frequency of 40 khz. INMETRO used the values [of the temperature and static pressure coefficients] described in Rasmussen, Metrologia 36 (1999). 4.. Microphone parameters Besides the open-circuit, free-field sensitivities of the microphones, the protocol requested that the laboratories should submit additional microphone parameters: Values of the temperature and pressure coefficients used in the calculations, and values of the acoustic centres used in the calculation, The values of the temperature and pressure coefficients reported by each laboratory are given in the following tables and 3. The values of the acoustic centre are given in the table 4 below. Table. Microphone parameters 4180.1395455 DFM CENAM KRISS NMIJ LNE PTB INMETRO Front cavity depth (mm) 0.488 0.50 0.516 N.R. N.R. N.R. 0.41 Temp. coeff. at 50 Hz (db/k) -0.0015-0.00-0.00546-0.00-0.001-0.001-0.001 Static press. coeff. at 50 Hz (db/kpa) -0.0057-0.00578-0.00116-0.006-0.005-0.005-0.005 Resonance frequency (Hz) 1800 1603 000 N.R. 3000 1400 N.R. a Note: This table is omitted here. The acoustic centres are presented in the following section.

Table 3. Microphone parameters 4180.1395456 DFM CENAM KRISS NMIJ LNE PTB INMETRO Front cavity depth (mm) 0.514 0.508 0.539 N.R. N.R. N.R. 0.500 Temp. coeff. at 50 Hz (db/k) 0.0004-0.00-0.00546-0.00-0.001-0.001-0.001 Static press. coeff. at 50 Hz (db/kpa) -0.0059-0.00613-0.00116-0.006-0.005-0.005-0.005 Resonance frequency (Hz) 100 1191 000 N.R. 3000 1400 N.R. Table 4. Acoustic centre (in mm) as a function of frequency reported by each laboratory. Frequency (Hz) DFM CENAM KRISS NMIJ LNE PTB INMETRO 1000 4.56 4.4 5.00 4.50 5.00 5.00 150 4.50 4.37 4.90 4.40 4.90 4.90 1600 4.41 4.8 4.80 4.30 4.80 4.80 000 4.30 4.18 4.70 4.0 4.70 4.70 500 4.14 4.04 4.40 4.00 4.40 4.00 4.40 3150 3.9 3.84 4.0 3.80 4.0 4.80 4.0 4000 3.61 3.57 3.90 3.50 3.90.10 3.90 5000 3.3 3.3 3.60 3.0 3.60.70 3.60 6300.74.80 3.0.80 3.0 3.10 3.0 8000.15.8.80.40.80.60.80 10000 1.57 1.76.30 1.90.30 1.90.30 1500 1.03 1.7 1.80 1.30 1.80 1.40 1.80 16000 0.60 0.83 1.0 0.70 1.0 0.90 1.0 0000 0.37 0.51 0.50 0.40 0.50 0.60 0.50 5000 0.01 0.03-0.3-0.30 0.00 0.0-0.01 31500-0.9-0.88-1. -1.10 0.00-1.00-0.87 40000-1.4-1.30 -. -.40-1.0-1.37 4.3. Microphone sensitivities Each laboratory reported their results after their own measurement round. Once all data was gathered, a preliminary analysis was made, and it showed that there were some discrepancies among laboratories. A few laboratories were asked to review their submission, and to resubmit if it was necessary, with an explanation of the reason for the review. All but one laboratory (NMIJ) declined to resubmit. NMIJ sent the reviewed data with the following explanation. NMIJ found two errors and corrected the free-field sensitivity levels and acoustic centres accordingly. (1) As described in Sec. 4.1, NMIJ used the virtual pulse method to minimize the influence of sound reflection. At 16 khz and above, the delay time to separate the direct sound from the reflected sound was wrongly set to be 5 ms instead of correct 0.5 ms. Therefore, the influence of sound reflection among the transmitter and receiver microphones still remained in the VTF (voltage transfer function). Acoustic centres were also incorrect due to the same reason. () At 31.5 khz, unintended data were acquired from the FFT analyzer because the frequency was specified by mistake. The following tables show frequency dependence of free-field sensitivity levels for each microphone and values of the acoustic centres before and after correction.

The open-circuit, free-field sensitivity level reported by each laboratory for each microphone at each frequency is given in the tables 5 and 6 below. The sensitivity level is given in db re 1 V/Pa. The uncertainties in db declared by each laboratory are listed in table 7. Table 5. Free-field sensitivity level of the microphone 1395455 in db re 1 V/Pa Frequency (Hz) DFM CENAM KRISS NMIJ LNE PTB INMETRO 1000-38.364-38.386-38.390-38.351-38.356-38.375 150-38.317-38.347-38.350-38.319-38.99-38.60 1600-38.36-38.7-38.80-38.50-38.05-38.313 000-38.10-38.159-38.180-38.18-38.093-38.137 500-37.938-37.981-38.00-37.957-37.93-38.017-37.875 3150-37.650-37.689-37.750-37.676-37.644-37.815-37.599 4000-37.18-37.33-37.310-37.43-37.07-37.06-37.58 5000-36.578-36.575-36.680-36.619-36.594-36.659-36.53 6300-35.57-35.593-35.740-35.603-35.587-35.701-35.531 8000-34.100-34.148-34.40-34.145-34.140-34.39-34.191 10000-3.31-3.378-3.500-3.36-3.360-3.45-3.74 1500-30.364-30.48-30.570 * -30.411-30.48-30.513-30.04 * 16000-9.100-9.077-9.300 * -9.197-9.146-9.35-9.169 0000-30.111-30.010 * -30.0-30.13-30.101-30.174-30.131 5000-33.103-3.964-33.070-33.048-3.993-33.153-33.154 31500-37.101-36.880 * -36.960-37.043-37.018-37.139-37.183 40000-4.39-4.11-4.080-4.31-4.517-4.46 * Note: The values marked with * have been excluded from the calculation of the Key Comparison Reference value (see 5.) Table 6. Free-field sensitivity level of the microphone 1395456 in db re 1 V/Pa Frequency (Hz) DFM CENAM KRISS NMIJ LNE PTB INMETRO 1000-38.359-38.417-38.370-38.356-38.345-38.350 150-38.313-38.383-38.340-38.30-38.97-38.50 1600-38.31-38.315-38.70-38.34-38.06-38.307 000-38.114-38.06-38.170-38.133-38.09-38.11 500-37.939-38.06-38.000-37.963-37.97-37.965-37.857 3150-37.671-37.731-37.730-37.683-37.653-37.75-37.593 4000-37.35-37.55-37.90-37.36-37.00-37.136-37.9 5000-36.574-36.579-36.670-36.590-36.560-36.536-36.471 6300-35.560-35.58-35.70 * -35.569-35.548-35.556-35.457 8000-34.054-34.103-34.170-34.076-34.068-34.070-34.073 10000-3.0-3.68-3.370-3.30-3.7-3.18-3.079 1500-30.171-30.5-30.370 * -30.194-30.07-30.17-9.909 * 16000-8.853-8.819-9.040 * -8.908-8.881-8.840-8.856 0000-9.951-9.85 * -30.040-9.970-9.937-9.99-9.977 5000-33.09-3.987-33.060-33.056-3.994-33.114-33.30 31500-37.196-36.994 * -37.100-37.15-37.13-37.185-37.385 * 40000-4.468-4.31-4.0-4.403-4.487-4.606 * Note: The values marked with * have been excluded from the calculation of the Key Comparison Reference value (see 5.)

Table 7. Expanded uncertainties in db using a coverage factor of k=, as declared by the participant laboratories Frequency (Hz) DFM CENAM KRISS NMIJ LNE PTB INMETRO 1000 0.07 0.86 0.4 0.08 0.4 0.13 150 0.07 0.5 0.0 0.08 0.19 0.13 1600 0.07 0.35 0.19 0.08 0.16 0.13 000 0.07 0.0 0.18 0.07 0.13 0.13 500 0.07 0.13 0.17 0.07 0.11 0.55 0.13 3150 0.07 0.10 0.17 0.07 0.10 0.49 0.1 4000 0.07 0.10 0.16 0.07 0.09 0.3 0.11 5000 0.07 0.10 0.16 0.07 0.08 0.0 0.11 6300 0.07 0.10 0.16 0.06 0.08 0.1 0.13 8000 0.07 0.10 0.16 0.06 0.08 0.1 0.14 10000 0.07 0.10 0.16 0.06 0.10 0.1 0.15 1500 0.07 0.10 0.16 0.07 0.13 0. 0.19 16000 0.08 0.10 0.16 0.09 0.16 0.6 0.19 0000 0.09 0.10 0.16 0.10 0.1 0.18 0.0 5000 0.10 0.11 0.18 0.10 0.6 0.18 0.0 31500 0.1 0.15 0.5 0.10 0.31 0.19 0. 40000 0.0 0.30 0.40 0.10 0. 0.4

5 Analysis of the results 5.1. Degrees of equivalence The nature of this comparison, with seven participants following a single loop with two travelling standards, should make it possible to use simple methods for estimating the comparison reference value such as those used in the comparisons CCAUV.A-K1 and CCAUV.A-K3. The methodology used in this report is similar to the one used in the Key Comparison CCAUV.A-K3 [1]. In general, a linear model described by Ey ( )= X ahas to be solved for each frequency. In the model a is a vector of parameters of the model, E(y) is the expectation of the measurements and X is the design matrix; these values are the values that should have been measured in absence of uncertainty. The vector y is conformed by the n measurement values provided by the participants on at least one of the circulated measurement objects. The elements of the design matrix are known a priori with zero uncertainty. The parameters a 1 a k, k 1, are unknown and have to be estimated from the n measurement results provided by the participants y, and the associated covariance matrix Σ. The covariance matrix Σ is the sum of two matrixes: Σ meas, which contains the square of the uncertainties claimed by the participants (diagonal elements) and the covariances between the provided measurement results (off-diagonal elements), and Σ object, which contains only diagonal elements that describe the estimated variance of the value of the measurand due to random instability of the circulated objects. Once these matrices have been built, the unknown parameters in the model and the reference values of the comparison can be estimated following the rest of the procedure described in reference []. In the analysis presented in this report it has been assumed that the correlation for two results y i and y j provided by different participants is equal to 0. In the case of two results y i and y j provided by the same participant, a common correlation coefficient for all participants has been assumed. The common correlation coefficient ρ was determined by the pilot laboratory using its own data, as described in Annex A. Additionally, the circulated measurement objects were considered stable, so that Σ object = 0 in the analysis. The covariance matrix and the matrix of the results of the participants are constructed using the data supplied by the participants. Because no further linking to any existing comparison is to be made, no additional matrixes with data from such a comparison are built. The model looks then as E M E E M E ( y ) 1,1 ( y1,7 ) ( y ),1 ( y ),7 1 M 1 = 0 M 0 0 M 0 a1, 1 a M 1 (1) where E(y i,j ) is the expected value of the i-th standard being measured in the j-th laboratory in the loop. The elements of the vector a i are the parameters to be estimated for the standard i. The covariance matrix is a 14 x 14 matrix that contains the declared uncertainties from each laboratory. The correlation value for measurements carried in the same laboratory is frequency dependent, and it was determined from the measurements made at the pilot laboratory. See Appendix A for a detailed explanation of the calculation of the correlations.

Once these matrices are complete, the estimates of the unknown parameters of the model can be calculated using: aˆ = CX C = T r 1 r T 1 1 ( X Σ X ), r Σ r y r r, () where y r and X r are obtained from y and X given in (1) by deleting the rows associated with the laboratories excluded from the calculation of the reference values, and Σ r is the covariance matrix associated with the reduced data set y r The reference values of the comparison and the associated covariance matrix should be calculated using: yˆ = Xaˆ, V (ˆ) y = XV (ˆ) a X T. (3) The degrees of equivalence are obtained using: V D = A ( y yˆ), T ( D ) = A V ( y yˆ ) A, ii ii T (4) where D ii is the degrees of equivalence per laboratory, A is an averaging matrix. The interlaboratory degrees of equivalence can be estimated using [1]: D = D D ij ii jj, ( ij ) = ii + jj ij ji, V D u u u u (5) where D ij is the inter-laboratory deviation. The u xx elements are obtained from equation (4). 5.. Equivalence and consistency of results As the participating laboratories are measuring the same item, it can be assumed that the reported results are drawn from a normal distribution. This hypothesis is tested by means of an observed χ - distributed estimator, as described in []. The degrees of freedom of the χ distribution are ν = n-k, where n is the number of measurements (14 for most frequencies) and k the number of standards ( LSaP microphones). To accept the hypothesis with a significance of 5 %, the probability P{ χ (ν)>χ obs} has to be larger than 5 %. During the course of the comparison, any anomalous results need to be identified by the pilot laboratory. Any laboratory whose results show deviations judged to be irregular must be contacted in order to have an opportunity to revise their results. Details of this process are given in section 4.3. Identifying discrepant measurements Once all data was collected, and the participants reported any problem with their measurement, the analysis method in [] was implemented. This method uses the normalized deviations that were calculated for every frequency and measurement. The method assumes that the normalized deviations are distributed as N(0,1). Therefore, if the modulus of a particular normalized deviation is

greater than, the corresponding measurement can be considered discrepant with a significance of 5 %. Handling of discrepant measurements If a measurement has been judged to be discrepant, it should be excluded from the least squares calculation of the reference values. Otherwise the calculated reference values cannot be considered to be proper estimates of the SI values of the quantities represented by the circulated objects. If there is a discrepant measurement y ij in a key comparison and the uncertainty assigned to this measurement is very small, the key comparison reference value will be attracted to this discrepant result. As a consequence the remaining results might appear to be discrepant as well even if they are mutually consistent. Since one discrepant measurement will always have a value d ij larger than the values d kj, i k of two or more mutually consistent values, the value y ij with the larger value d ij should be excluded first. A repeated least squares adjustment of the reference values will then show if there are further discrepant results. The procedure of excluding discrepant results one by one can be summarized as follows: 1. Identify the result y ij with the largest value d ij > in the (reduced) set of results.. Exclude the result y ij from the reduced set of results. 3. Repeat the least squares adjustment of reference values for the reduced set of results. 4. If the results in the reduced set are not mutually consistent, continue at point 1. It has to be emphasized, that the discrepant results are excluded only from the calculation of the reference values and not from the key comparison as such. It is therefore still necessary to calculate the deviations of the discrepant results from the reference values and the uncertainties of these deviations. Reference [] gives a detailed explanation of the procedure.

6 Results 6.1. Reference values As in previous comparisons, there is not a single key comparison reference value, but one per frequency associated with each one of the standards circulated. The reference values are listed below in table 8. These values were obtained using the procedure described in the previous section. Some discrepant results were found at some frequencies. In all cases where a discrepant measurement was found, this measurement result or results have been excluded from the determination of the reference values (see tables 5 and 6 in section 4.3). Table 8. Reference values for the comparison. Free-field sensitivity levels M f in db re 1V/Pa and expanded uncertainty U (k = ) in db. Freq. 4180.1395455 4180.1395456 (Hz) M f U k = M f U k = 1000-38.36 0.047-38.357 0.047 150-38.31 0.046-38.309 0.046 1600-38.50 0.045-38.43 0.045 000-38.17 0.041-38.16 0.041 500-37.945 0.039-37.949 0.039 3150-37.663 0.037-37.676 0.037 4000-37.34 0.036-37.33 0.036 5000-36.594 0.035-36.569 0.035 6300-35.596 0.034-35.561 0.035 8000-34.143 0.034-34.076 0.034 10000-3.354 0.036-3.4 0.036 1500-30.403 0.041-30.19 0.041 16000-9.133 0.047-8.863 0.047 0000-30.136 0.054-9.966 0.054 5000-33.059 0.051-33.065 0.051 31500-37.073 0.064-37.17 0.065 40000-4.343 0.074-4.44 0.074 6. Degrees of equivalence per laboratory The deviation of the measurements, that is, the reported measurements in table 5 and 6 minus their respective reference value in table 8 are plotted in figure.

Figure. Deviations of the free-field sensitivity levels of each microphone reported by the participant laboratories: (a) microphone 4180.1395455, and (b) microphone 4180.1395456. The degrees of equivalence per laboratory were determined as mentioned in the previous section (equation (4)). Figure 3 below shows average deviation in db per laboratory as a function of frequency. Tables 9 and 10 present the degrees of equivalence and their uncertainties in tabular form. Figure 3. Average deviation in db per laboratory and frequency.

Table 9. Degrees of equivalence per laboratory as function of frequency: deviations (db). Frequency (Hz) DFM CENAM KRISS NMIJ LNE PTB INMETRO 1000-0.00-0.04-0.01 0.006 0.009 - -0.003 150-0.005-0.055-0.035-0.009 0.01-0.055 1600 0.013-0.047-0.09 0.004 0.041 - -0.064 000 0.009-0.056-0.049-0.004 0.033 - -0.003 500 0.008-0.057-0.063-0.013 0.0-0.044 0.081 3150 0.009-0.040-0.070-0.010 0.01-0.114 0.074 4000 0.007-0.011-0.067-0.006 0.030 0.06-0.010 5000 0.005 0.004-0.094-0.03 0.005-0.016 0.080 6300 0.01-0.009-0.15-0.008 0.011-0.050 0.084 8000 0.033-0.016-0.095-0.001 0.006-0.045-0.0 10000 0.03-0.034-0.146-0.007-0.004-0.014 0.113 1500 0.030-0.09-0.173-0.005-0.00-0.045 0.41 16000 0.0 0.050-0.17-0.055-0.016-0.040-0.015 0000 0.00 0.10-0.079 0.000 0.03-0.001-0.003 5000-0.036 0.086-0.003 0.010 0.069-0.07-0.130 31500-0.05 0.187 0.094 0.06 0.054-0.038-0.160 40000-0.046 0.07 0.33 0.01 - -0.119-0.133 Table 10. Degrees of equivalence per laboratory as function of frequency: uncertainty k = (db) Frequency (Hz) DFM CENAM KRISS NMIJ LNE PTB INMETRO 1000 0.049 0.808 0.1 0.061 0.1-0.114 150 0.050 0.487 0.183 0.06 0.173-0.114 1600 0.050 0.36 0.174 0.06 0.144-0.115 000 0.053 0.184 0.165 0.053 0.116-0.116 500 0.054 0.117 0.156 0.054 0.097 0.516 0.117 3150 0.056 0.087 0.156 0.056 0.087 0.459 0.107 4000 0.056 0.088 0.147 0.056 0.077 0.99 0.098 5000 0.056 0.087 0.145 0.056 0.067 0.18 0.097 6300 0.057 0.087 0.146 0.046 0.067 0.19 0.116 8000 0.057 0.087 0.145 0.046 0.067 0.19 0.16 10000 0.051 0.079 0.13 0.041 0.079 0.175 0.13 1500 0.048 0.077 0.140 0.048 0.104 0.183 0.164 16000 0.054 0.073 0.137 0.063 0.16 0.11 0.15 0000 0.060 0.095 0.15 0.070 0.169 0.143 0.160 5000 0.076 0.086 0.153 0.076 0.6 0.153 0.171 31500 0.095 0.153 0.6 0.07 0.84 0.167 0.199 40000 0.153 0.39 0.34 0.055-0.171 0.188

The degrees of equivalence calculated for the frequencies,.5 khz, 4.0 khz, 8.0 khz, and 0.0 khz are shown in figure 4.

Figure 4 Inter-laboratory degrees of equivalence at.5 khz, 4 khz, 8 khz, and 0.0 khz. Table 11 shows the results of the χ test applied to the comparison results for specified values of the correlation coefficient between results provided by the same laboratory. As demonstrated in [1], when high correlation coefficients are used, the test tends to fail. This is simply because a high correlation forces the differences to be the same for the two standards, and that may not be the case. Although at some frequencies the result of the χ test is low, in general this test has been considered satisfactory when the correlation between measurements of the same laboratory takes the values tabulated in Appendix A. Table 11. Results of the χ test as described in []. The correlation coefficients used in the calculations are shown in the table in Appendix A. Frequency (Hz) P{ χ (ν)>χ obs} 1000 100.0% 150 100.0% 1600 99.6% 000 99.7% 500 95.6% 3150 94.0% 4000 98.7% 5000 74.5% 6300 57.1% 8000 87.5% 10000 45.7% 1500 90.7% 16000 84.6% 0000 99.1% 5000 76.0% 31500 94.8% 40000 41.1%

6.3 Degrees of equivalence inter-laboratory Tables 1, 13, 14 and 15 contain the inter-laboratory degrees of equivalence of this comparison. Table 1. Inter-laboratory degrees of equivalence at.5 khz. The upper triangle contains the differences (db); the lower triangle contains the uncertainties k = (db).5 khz DFM CENAM KRISS NMIJ LNE PTB INMETRO DFM - 0.065 0.07 0.0-0.014 0.053-0.07 CENAM 0.139-0.006-0.044-0.079-0.01-0.138 KRISS 0.173 0.01 - -0.050-0.085-0.019-0.144 NMIJ 0.093 0.139 0.173 - -0.035 0.031-0.094 LNE 0.13 0.160 0.190 0.13-0.066-0.059 PTB 0.51 0.531 0.541 0.51 0.57 - -0.15 INMETRO 0.139 0.173 0.01 0.139 0.160 0.531 - Table 13. Inter-laboratory degrees of equivalence at 4.0 khz. The upper triangle contains the differences (db); the lower triangle contains the uncertainties k = (db) 4.0 khz DFM CENAM KRISS NMIJ LNE PTB INMETRO DFM - 0.017 0.073 0.013-0.03-0.056 0.017 CENAM 0.115-0.056-0.005-0.040-0.073 0.000 KRISS 0.164 0.177 - -0.061-0.096-0.19-0.056 NMIJ 0.093 0.115 0.164 - -0.036-0.068 0.004 LNE 0.107 0.17 0.173 0.107 - -0.033 0.040 PTB 0.308 0.315 0.336 0.308 0.313-0.07 INMETRO 0.13 0.140 0.183 0.13 0.134 0.318 - Table 14. Inter-laboratory degrees of equivalence at 8.0 khz. The upper triangle contains the differences (db); the lower triangle contains the uncertainties k = (db) 8.0 khz DFM CENAM KRISS NMIJ LNE PTB INMETRO DFM - 0.048 0.18 0.033 0.07 0.077 0.055 CENAM 0.113-0.079-0.015-0.01 0.09 0.006 KRISS 0.16 0.175 - -0.095-0.101-0.051-0.073 NMIJ 0.085 0.108 0.158 - -0.006 0.044 0.01 LNE 0.098 0.119 0.166 0.093-0.050 0.08 PTB 0.05 0.15 0.45 0.0 0.08 - -0.0 INMETRO 0.145 0.159 0.197 0.141 0.149 0.34 -

Table 15. Inter-laboratory degrees of equivalence at 0.0 khz. The upper triangle contains the differences (db); the lower triangle contains the uncertainties k = (db) 0.0 khz DFM CENAM KRISS NMIJ LNE PTB INMETRO DFM - -0.100 0.099 0.00-0.011 0.01 0.03 CENAM 0.11-0.199 0.10 0.088 0.10 0.13 KRISS 0.153 0.157 - -0.079-0.111-0.079-0.076 NMIJ 0.11 0.118 0.157 - -0.03 0.000 0.003 LNE 0.190 0.194 0.0 0.194-0.03 0.035 PTB 0.167 0.171 0.00 0.171 0.30-0.003 INMETRO 0.18 0.186 0.13 0.186 0.41 0.4 - Page 6 of 44

7 Conclusions The aim of the Key Comparison CCAUV.A-K4 is to establish reliable equivalence between National Metrology Institutes for free-field calibration of Laboratory Standard microphones of the type LSaP. The initial plans involved nine participants. However, once the comparison was underway, two laboratories decided to decline from participating. Two LSaP microphones were circulated among the participants in a single star configuration. All participants reported the free-field sensitivity and other microphone parameters such as acoustic centre, environmental coefficients, etc. The reported sensitivities have been analysed using a least-squares technique similar to the one used in the Key Comparison CCAUV.A-K3. All laboratories were included in the estimation of the reference values. There is no preceding comparison with which the obtained degrees of equivalence can be compared. 8 References [1] V. Cutanda-Henríquez, and K. Rasmussen, Final Report on the Key Comparison CCAUV.A- K3, Centro Nacional de Metrología, México, Danish Primary Laboratory for Acoustics, Denmark, May 006. [] L. Nielsen, Identification and handling of discrepant measurements in key comparisons, Danish Institute of Fundamental Metrology, DFM-01-R8, 00. Page 7 of 44

Appendix A Determination of the correlation coefficients The correlation coefficients used in the analysis were determined from a relatively large set of freefield sensitivity data gathered at DFM from all the control measurements carried out on the travelling standards. The DFM s own uncertainty budget was also used in the determination of the correlation coefficients. In principle, this calculation should only be valid for DFM s correlation coefficients, but it has been chosen to use them for all laboratories. In the following, all operations are referred to this set of measurements. The covariance was determined as the sum of the contributors to the uncertainty of the measurement that affect the sensitivities of both microphones in a similar manner, including the reproducibility, that is, the type B contributors. Some of these contributors are: the environmental variables, the polarising voltage, the reference impedance, the effective distance between microphones, and the mathematical manipulations. ( M f M f ) = ui Cov,, B 1 (A.1) where M f1 and M f are the free-field sensitivities of the microphones, and u i is the i th uncertainty type B component in the uncertainty budget. In order to have a more realistic estimate of the experimental reproducibility, the experimental correlation coefficient of the sensitivities of the two microphones, ρ xy, was calculated, and used for determining the experimental covariance by means of: (, ) ρ xy σ σ, Cov A M f 1 M f = 1 (A.) where σ 1, and σ are the standard deviations of the free-field sensitivity of each microphone. Finally, the correlation coefficient at each frequency was determined using the sum of all covariances divided by the product of the combined variance of the two microphones (which is the same for the two microphones): Cov = ( M f 1, M f ) + Cov B ( M f 1, M f ). A ρ (A.3) umf 1 umf The correlation coefficients obtained by such a procedure are presented in the table below. Frequency (Hz) Correlation Frequency (Hz) Correlation 1000 0.769 8000 0.716 150 0.769 10000 0.47 1600 0.769 1500 0.47 000 0.769 16000 0.356 500 0.769 0000 0.384 3150 0.769 5000 0.573 4000 0.769 31500 0.750 5000 0.716 40000 0.358 6300 0.716 Page 8 of 44

Appendix B INMETRO corrections to reported results During the process of reviewing the draft A, INMETRO requested to introduce some corrections to the data originally reported. Table B1 contains the original and revised sensitivities of the two microphones. The reason given by INMETRO for this correction is the following: The sensitivity results reported in the cancelled documents (i.e. calibration certificates DIMCI 009/008 and 010/008) where erroneously referred to the nominal centre frequencies of 3rdoctave band, as requested by the measurement protocol of the comparison CCAUV.A-K4. It was detected that these sensitivity results refer to the exact centre frequencies of 3rd-octave bands of base, which are obtained by the equation 1000*N/3, where N is an integer. This situation implies that some frequencies match but others are different. A coincidence between the nominal and exact frequencies within the frequency range of the comparison is obtained at 1 khz, khz, 4 khz, 8 khz and 16 khz. Unfortunately, the INMETRO only perceived this fact after the confirmation sent to you. The mismatch between some of the frequencies reported and the ones requested in the protocol implies in significant differences that required this communication. Table B1 Original and revised sensitivity level data for the two microphones. Sensitivity level (db re 1 V/Pa) Freq. 4180.1395455 4180.1395456 (Hz) Original Revised Original Revised 1000-38.375-38.375-38.350-38.350 150-38.60-38.38-38.50-38.319 1600-38.313-38.5-38.307-38.34 000-38.137-38.137-38.11-38.11 500-37.875-37.951-37.857-37.933 3150-37.599-37.678-37.593-37.671 4000-37.58-37.58-37.9-37.9 5000-36.53-36.68-36.471-36.565 6300-35.531-35.639-35.457-35.565 8000-34.191-34.191-34.073-34.073 10000-3.74-3.411-3.079-3.17 1500-30.04-30.469-9.909-30.178 16000-9.169-9.169-8.856-8.856 0000-30.131-30.14-9.977-9.960 5000-33.154-33.037-33.30-33.097 31500-37.183-37.01-37.385-37.15 40000-4.46-4.99-4.606-4.476 Page 9 of 44

Appendix C Uncertainty budgets Uncertainty budgets were requested of all participants. The budgets are reproduced here as they were received from the participants. No editing of the contents but a small amount of reformatting has been performed. Table, equation and figure numbers do not follow the numbering in the report. Page 30 of 44

DFM The condensed uncertainty budget for a free-field reciprocity calibration of B&K type 4180 microphones are given in the table below. The free-field sensitivity is determined as the average of the results derived using four measuring distances between 170 30 mm for type B&K 4180. The background for the budget is given in the following remarks: 1. The figures represent the combined effects of the uncertainty on the distance between the microphone diaphragms (0.05 mm), front cavity depth (0.05 mm) and acoustic centres (0.5 1.0 mm) referred to the average distance of 40 mm.. The figures represent the combined effects of the measurement uncertainties on static pressure (50 Pa), temperature (0.5 K) and relative humidity (10 %). 3. The figures represent the combined effects of the measurement uncertainties on static pressure (50 Pa), temperature (0.5 K) and relative humidity (10 %). 4. The figures represent the effect of the basic uncertainty on the equations for calculating the absorption and density of humid air. 5. The figures represent the combined effects of the uncertainty on the measurement impedance and voltage ratios (3 ratios each derived from 4 voltage measurements). 6. The figures represent the uncertainty on the polarizing voltage (40 mv) and the non-linear relation between polarizing voltage and microphone sensitivity. 7. The figures represent the possible errors introduced by the signal manipulations, such as digital filtering, inverse FFT-applications and subsequent removal of reflections in the time domain etc. 10. The figures represent the uncertainty on applying a correction for dependence of static pressure and temperature on the microphone sensitivity. (Uncertainty on resonance frequency, on low-frequency value of static pressure coefficient and on temperature coefficient, allowed deviation from reference conditions 3 kpa respectively 1 K). Page 31 of 44

Uncertainty budget for free-field reciprocity calibration of B&K 4180 microphones Components of type B uncertainties in db*1000 Frequency khz 3 5 10 15 0 5 30 40 Source 1 Effective distance 5 9 39 Density of air 19 19 19 19 19 19 19 19 3 Air absorption 1 8 15 18 17 1 0 4 Environmental equations 0 1 4 6 9 11 0 5 Electrical transfer impedance 40 40 0 0 0 0 30 40 6 Polarizing voltage........ 7 Mathematical manipulations 18 7 8 Rounding of results 5 5 5 5 5 5 5 5 σ ( ) 1/ B = u /3 30.5 9.0 1..5 3.4 4.4 8.3 37.9 Components of type A uncertainties in db*1000 9 Allowed reproducibility 15 0 0 0 30 35 50 85 Overall uncertainty in db*1000 at measurement conditions ( ) 1/ σ = σ + σ A B 34.0 35.3 9.1 30.1 38.0 4.7 57.4 93.1 Uncertainty on corrections to reference environmental conditions in db*1000 10 Correction to reference conditions 3 5 9 13 18 0 30 Overall uncertainty in db*1000 at reference conditions (k=) 34.1 35.4 9.6 31.4 40. 46.3 60.8 97.8 Expanded uncertainty in db 0.068 0.071 0.059 0.063 0.080 0.093 0.1 0.0 Declared CMC values 0.07 0.07 0.07 0.08 0.09 0.1 0.1 0.0 * * Not included in the CMC Page 3 of 44

CENAM

KRISS 1000 150 1600 000 500 3150 4000 5000 6300 8000 10000 1500 16000 0000 5000 31500 40000 Electrical Transfer Impedance Symbol Uncertainty Unit 0.009 0.0035 0.0035 0.0036 0.0050 0.0050 0.0050 0.0050 0.0050 0.0050 0.0069 0.010 0.010 0.010 0.010 0.0149 0.0149 Series Capacitor C Table nf 0.0017 0.006 0.006 0.006 0.0044 0.0044 0.0044 0.0044 0.0044 0.0044 0.0044 0.0087 0.0087 0.0087 0.0087 0.0087 0.0087 Voltage Ratio VR Table 0.000 0.000 0.000 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.005 0.005 0.005 0.005 0.005 0.011 0.011 Inherent Noise Measured 0.0009 0.0009 0.0009 0.0009 0.0009 0.0009 0.0009 0.0009 0.0009 0.0009 0.0009 0.0009 0.0009 0.0009 0.0009 0.0009 0.0009 Distortion Measured 0.0004 0.0004 0.0004 0.0004 0.0004 0.0004 0.0004 0.0004 0.0004 0.0004 0.0004 0.0004 0.0004 0.0004 0.0004 0.0004 0.0004 Frequency f Table Hz 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0001 0.0005 0.0000 Reveiver Ground Shield B&K 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 Transmitter Ground Shield B&K 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 0.0003 Microphone Parameters 0.0307 0.0306 0.0306 0.0305 0.0304 0.0304 0.030 0.0301 0.099 0.097 0.095 0.093 0.090 0.086 0.08 0.073 0.056 Acoustic Center 1 mm 0.0307 0.0306 0.0306 0.0305 0.0304 0.0304 0.030 0.0301 0.099 0.097 0.095 0.093 0.090 0.086 0.08 0.073 0.056 Polarizing Voltage Uo 0.0064 V 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 0.0005 Measurement Conditions 0.0049 0.0049 0.0048 0.0048 0.0048 0.0048 0.0048 0.0047 0.0047 0.0047 0.0046 0.0047 0.0048 0.0053 0.0064 0.0300 0.006 Distance between Microphones l 0.1 mm 0.0048 0.0048 0.0048 0.0048 0.0048 0.0047 0.0047 0.0047 0.0046 0.0046 0.0046 0.0045 0.0045 0.0044 Static Pressure Ps 0.0137 kpa 0.0006 0.0006 0.0006 0.0006 0.0006 0.0006 0.0006 0.0006 0.0006 0.0006 0.0006 0.0006 0.0006 0.0006 0.0006 0.0005 0.0006 Temperature T 0.1733 K 0.0000 0.0000 0.0000 0.0000 0.0000 0.0001 0.0001 0.0001 0.000 0.0003 0.0005 0.0009 0.0016 0.008 0.0061 0.087 0.004 Relative Humidity RH 1.9370 % 0.0001 0.0001 0.0001 0.0001 0.0001 0.0001 0.0001 0.0001 0.000 0.000 0.0003 0.0004 0.0006 0.0009 0.0019 0.0086 0.0006 Processing of Results 0.1119 0.0944 0.0861 0.081 0.078 0.0744 0.0708 0.0708 0.0709 0.0691 0.0667 0.0677 0.0654 0.0640 0.083 0.1153 0.1966 Rounding Error 0.0050 db 0.009 0.009 0.009 0.009 0.009 0.009 0.009 0.009 0.009 0.009 0.009 0.009 0.009 0.009 0.009 0.009 0.009 Repeatability of Measurements Table 0.1000 0.0800 0.0700 0.0650 0.0600 0.0550 0.0500 0.0500 0.0500 0.0473 0.043 0.0437 0.0396 0.0356 0.0596 0.0984 0.194 Static Pressure Corrections Table 0.001 0.001 0.0011 0.0011 0.0011 0.0010 0.0009 0.0007 0.0003 0.0006 0.000 0.003 0.0004 0.0107 0.061 0.0308 0.09 Temperature Corrections Table 0.0009 0.0010 0.0011 0.001 0.0014 0.0017 0.001 0.006 0.0035 0.0051 0.0083 0.013 0.0143 0.0143 0.0059 0.016 0.008 Sjgnal Processing Error 0.0500 0.0500 0.0500 0.0500 0.0500 0.0500 0.0500 0.0500 0.0500 0.0500 0.0500 0.0500 0.0500 0.0500 0.0500 0.0500 0.,05 Sum 0.1161 0.0994 0.0916 0.0878 0.084 0.0807 0.0773 0.0773 0.077 0.0755 0.0734 0.0746 0.075 0.0710 0.0878 0.13 0.1988 Expanded Uncertainty (k=) 0.33 0.1988 0.1831 0.1755 0.1683 0.1614 0.1546 0.1545 0.1544 0.1510 0.1468 0.149 0.1449 0.141 0.1757 0.463 0.3976 Stated Uncertainty 0.4 0.0 0.19 0.18 0.17 0.17 0.16 0.16 0.16 0.16 0.16 0.16 0.16 0.16 0.18 0.5 0.40

NMIJ Page 35 of 44

LNE Free field calibration of microphones by the reciprocity technique Uncertainty evaluation J-N Durocher, D. Rodrigues, LNE july 007 Microphone sensitivity product The module of the microphone sensitivity product is given by (microphone 1 as receiver and as emitter): M αd' U1 D' 0 1. M = 10 I ρ1 f U 1 = Z is the measured electrical transfert impedance e With 1 I - f is the frequency - D is the distance between the acoustic center of the microphones D =D-δ where δ is the acoustic center position (relative to the diaphragm) and D the distance between the microphone diaphragms (δ is assumed independent of the microphone). - ρ1 is the air density during the measurement of Z e1 - α is the air absorbtion (assumed constant during the 3 Z e measurements) Expressed in db, the above microphone sensitivity product become : ( M 1M ) = L1 + CD CR + C f + D K sys L + L = lg ' + 1 0 1 α with C = lg( D') C = 0lg( ρ ) C 0lg( / ) D 0 R1 1 f = f And where K sys corresponds to the systematic deviation from the ideal conditions (K sys is a null correction only used for uncertainties calculation) K + Sys = K pole + K polr + Kϑ E + KϑR + K Echos + K Elec K Air with : - Polarization voltage of the two microphones K pol - Angular deviation of the two microphones from the main axis K θ - Acoustical disturbances K echo - Electrical disturbances K elec - unknown deviation K unk Electrical transfert impedance measurement The elementary measurement is a voltage ratio measurement : U U R1 = ω E C U1g I g R1 E And thus Z e1 U = I 1 U R1g E = CωU g E R1 Page 36 of 44

Where : - g R1 corresponds to the gain of the receiver microphone amplifier (insert voltage with the microphone n 1) and G R1 = 0.lg(g R1 ) - g E corresponds to the gain of the emmiter microphone amplifier and G E = 0.lg(g E ) - Z e1 corresponds to the electrical transfert impedance with the microphone microphone n 1 and n ) and L 1 = 0.lg(Z e1 ) - C corresponds to the electrical capacitance used for the current measurement and CC is the corresponding correction (db) = 0lg( Cω) C C The electrical transfert impedance expressed may be written in db : L ( Ze 1 ) = G1 + GE GR Att CC 1 = 0 lg 1 Where Att corresponds to the attenuation (db) of the attenuator used for the insert voltage method. For the following microphone configurations 13 and 3 : L 13 = G13 + GE GR 1 Att CC L3 = G3 + GE GR Att CC and thus, the electrical transfert impedance combination : L + L L G = + G G G + G GR + Att C 1 13 3 1 13 3 E C R1 Individual microphone sensitivity The sensitivity of the microphone n 1 can be written as : L + L L C + C C + C + αd' + K 1 13 3 R1 R13 R3 D f sys L 1 = + + C C proc where C proc is the equivalent correction to process the results. ( G1 + G13 G3 ) + GE GR 1 + GR Att CC ( CR1 + CR 13 CR ) + CD + C f + D' + C proc K sys 1 3 α L = + Uncertainty budget : 4u ( L1 ) = 7u ( G) + 3u noise + u ( Att) + u ( CC ) + 3u ( CR ) + u ( CD ) + u ( C f ) + u ( αd') + u ( K sys ) + 4u proc Uncertainties evaluation Page 37 of 44

Basic parameters measurement Quantity Mean value Unit U k u Distance : D (between diaphragms) 0,3 m 1.10-3 3 0,58.10-3 Temperature : T 3 C 1 3 0,58 Pressure : P 1013 hpa 1 3 0,58 Relative humidity : HR 60 % 5 3,89 Voltage ratio measurement : G - db 0,01 3 0,0058 Noise component The effect of the acoustic noise depends on the frequency and on the distance between the microphones. For LS microphones, the standard uncertainty u noise is estimated as : u D ( db) = 0,133 0,3 noise F khz where F khz = f / 1000 Attenuator component The maximum uncertainty due to the attenuator is estimated as 0,01 db. u att ( db) = 0,01 3 = 0,0058 Capacitor component The calibration of the capacitor gives C = 4744,0 pf ± 0,01 % ± 0,1 pf (k=) It results : u C = 0,001 / = 0,00061 db Air density component Using the formulation of the air density given in IEC 61094- Annex F and the usual air conditions, we obtain : u( C Rmn ) 3 4 ( 8,57.10 ). u ( P ) + (,93.10 ). u ( T ) + ( 9,15.10 ). u ( Hr) = db hpa Page 38 of 44

and using the uncertainties given in 1 : u ( ) = 0,018 db C Rmn Distance between microphones component Using equation D' = D δ and D >> δ we find, u( D') D' u( D) D u( δ ) + D Using the approximation 0lg(1 + x) 8. 7x u D D' + u( D') u( D') u( D) u( δ ) = 0lg 8.7 8.7 + D' D' D D db Assuming that the maximum uncertainty of the acoustic center position is 1 mm, we have 0, 011 u D D = db Air absorbtion component u ( αd') = [ D' u( α )] + [ αu( D')] u( α ) D Assuming that the air absorption coefficient α(f) is defined with a maximum relative uncertainty of 10 %, while α 4 0,004 + 8,63.10 FkHz 0,1 D u + 3 0,3 [ 7 1,5 F ] 5 ( ) D = αd = 10 khz α db Result processing component 0,001 Rounding the results component : u = 0, 0056 db proc1 = 3 Repeatability is estimated as : u 001 proc = 0,005 + 0, F khz db Transposition to the reference conditions. The correction is C Cond = CPΔ P + CT ΔT with Δ P = P Pref ΔT = T Tref u Cond = C P u ( P) + u ( C ) u ( ΔP) + C u ( T ) + u ( C ) u ( ΔT ) P T T Page 39 of 44

where u(δp), u(δt) are the standard deviations corresponding to the ranges of static pressure or temperature allowed for the calibration. Using the raw approximations : C P = 0,0006 + 0,00001 F khz db/hpa and C T = 0,005 + 0,00067 F khz db/k And assuming that u(c P )/C P = u(c T )/C T = 10 % Using the uncertainties listed in 1, we obtain : u u 5 ( P) = 0,0017 + 3,5. FkHz proc3 = CPu 10 db 4 ( T ) = 0,006 + 7,7. FkHz proc4 = CTu 10 db And considering the extreme variations allowed : 00 3 u( ΔP) = = 115 hpa ( ΔT ) = = 1, 7 3 3 u C we obtain : u u 5 ( ΔP) = 0,004 + 8. FkHz proc5 = u( CP ). u 10 db 5 ( ΔT ) = 0,001+ 6,7. FkHz proc6 = u( CT ). u 10 db Polarization component Assuming that u( V 0, Pol ) = 3 V u VPol + u( VPol ) 0, = 0 lg 8,7 = 0, V Pol 00 3 Pol 005 db Angular deviation component Using a directivity law based on the cylindrical piston in a plane baffle, transposed on the size of the LS microphones, we obtain the approximation : u u( ϑ).36.10 6 F khz db ϑ with ϑ expressed in degree Using 6 u ( ϑ) = deg, we obtain uϑ 4.10 FkHz db 3 Acoustical disturbances component Page 40 of 44

Acoustical disturbances are not completely removed by the treatment and the treatment will cause a systematic error in the measured sensitivity product. Simulation with LS microphones give u echo = 0,01 db Electrical disturbances component The estimation for the total harmonic distortion is u = 0,01 3 = 0, 0058 db THD Crosstalk is partially removed by the acoustical treatment. The estimation for the residual component is : u 0,005 + 0, 005 = db cross F khz Unknown component Free field calibration comparisons carried out in the past showed systematic deviations in the high frequency range (f 10 khz). Causes are not well known but may be one of the following : Air composition in anechoic rooms Temperature increase near to the receiver microphone Crosstalk in the emitter microphone circuit In the absence of further explanations, deviations observed are take into account. The deviation follow approximately : ( F /8) u unk = 0,3. Lg khz db for frequencies > 8 khz Page 41 of 44

Uncertainties calculation Frequency Microphone sensitivity uncertainties (db) [k=] khz D = 5 cm D = 30 cm D = 35 cm 1,00 0,17 0,4 0,3 1,5 0,14 0,19 0,5 1,6 0,1 0,16 0,0,0 0,11 0,13 0,17,5 0,09 0,11 0,14 3,15 0,09 0,10 0,11 4,0 0,08 0,09 0,10 5,0 0,08 0,08 0,09 6,3 0,08 0,08 0,08 8,0 0,09 0,08 0,08 10 0,10 0,10 0,10 1,5 0,13 0,13 0,13 16 0,17 0,16 0,16 0 0,1 0,1 0,1 5 0,6 0,6 0,6 31,5 0,31 0,31 0,31 40 0,38 0,38 0,38 Page 4 of 44

PTB