Exploiting the Leaky-wave Properties of Transmission-line Metamaterials for Single-microphone Direction Finding

A transmission-line acoustic metamaterial is an engineered, periodic arrangement of relatively small unit-cells, the acoustic properties of which can be manipulated to achieve anomalous physical behaviours. These exotic properties open the door to practical applications, such as an acoustic leaky-wave antenna, through the implementation of radiating channels along the metamaterial. In the transmitting mode, such a leaky-wave antenna is capable of steering sound waves in frequency-dependent directions. Used in reverse, the antenna presents a well defined direction-frequency behaviour. In this paper, an acoustic leaky-wave structure is presented in the receiving mode. It is shown that it behaves as a sound source direction-finding device using only one sensor. After a general introduction of the acoustic leaky-wave antenna concept, its radiation pattern and radiation efficiency are expressed in closed form. Then, numerical simulations and experimental assessments of the proposed transmission-line based structure, implementing only one sensor at one termination, are presented. It is shown that such a structure is capable of finding the direction of an incoming sound wave, from backward to forward, based on received sound power spectra. This introduces the concept of sound source localization without resorting to beam-steering techniques based on multiple sensors.


I. INTRODUCTION
The advent of metamaterials has opened new doors to physical and engineering applications.Following the first investigations and developments in electromagnetics, applications in other fields of physics were envisaged for such artificial structures such as optics, acoustics, mechanics, and thermodynamics.In acoustics, many novel concepts are now proposed that may prove to be beyond what was previously thought possible.Cloaking, 1 super lensing, 2 super absorbers, 3 metasurfaces, 4 or non-reciprocal devices 5,6 are some such examples.
The study on locally resonant sonic materials by Liu et al. 7 may be considered as the first work on acoustic metamaterials with subwavelength unit-cells.They experimentally showed that lead balls coated with a layer of silicon rubber can be used to achieve effective negative elastic modulus.Then a similar structure was proved to exhibit doublenegative acoustic properties. 8However, the resonant behaviour of these materials, which are characterized by a short bandwidth and high losses, encouraged investigating a new class of transmission-line (TL) based acoustic metamaterials.Following the work of Fang et al., 9 a new unit-cell was proposed consisting of a tube with a side hole, achieving negative bulk modulus over a large bandwidth. 10Side holes along with clamped membranes which exhibit negative dynamic mass 11,12 provided necessary and sufficient means of obtaining composite acoustic media with simultaneously negative mass density and bulk modulus. 13Then, TL modelling was proposed to design acoustic metamaterials with no band gap, and leaky-wave antenna (LWA) applications were first envisaged. 14Later, the theory of acoustic LWA was experimentally validated 15 and new applications were reported, such as direction finding 16 and acoustic dispersive prism. 17ecently, Xie et al. have brought the idea of single sensor sound localisation into the limelight by introducing single-sensor cocktail party listener using acoustic metamaterials, capable of separating broadband signals produced by three separate sources. 18In the present article, we propose an alternative technique using leaky-wave antennas based on an acoustic TL metamaterial to locate the direction of sound sources with one single microphone.
The leaky-wave appellation refers to the principle of power leakage along a waveguide and it has been first developed in the 1940s (Ref.19) in the electromagnetic realm.Owing to the high directivity and frequency scanning capability of such an antenna, it became very attractive in electromagnetics.However, LWAs based on conventional righthand materials exhibit only positive wavenumbers and, as a consequence, have the drawback of scanning only the halfspace from broadside to endfire (0, 90 ).With the introduction of composite right/left-hand (CRLH) TL-metamaterials, which support both positive and negative wavenumbers, the backfire to endfire (-90 , 90 ) electromagnetic antenna has been made possible by design. 20n this article the acoustic counterpart of electromagnetic LWA, exhibiting backfire to endfire radiation/reception, is exploited for direction finding.The structure of this article is as follows: first, the acoustic LWA theory and mechanism are presented and mathematically formulated to get the desired radiation efficiency and directivity.Then, a new concept of single microphone source localisation is introduced, with theoretical and numerical simulations highlighting the phenomenon.Finally, a frequency-dependent directionality is derived from experimental assessments, followed by a discussion on the practical use of such a structure.

A. General description and properties
An acoustic LWA is a waveguiding structure that favours the sound power leakage through radiating openings as it propagates along its axis.Leaky-waves are usually employed in antenna applications, since the leakage phenomenon is generally associated with high directivity and therefore, their size is not constrained by the operating wavelength, but rather by the directivity specifications. 21A simple leaky-wave structure is schematically represented in Fig. 1.
In the transmission mode (for which a source is emitting sound inside the waveguide), the sound pressure field leaking out of the waveguide, in free-field conditions, follows the form 20 pðy; zÞ ¼ p 0 e Àcy e Àjk z z ; (1) where c ¼ a þ jb is the complex propagation constant of the wave along the y direction, and k z is the propagation constant perpendicular to this direction, related to b by (Fig. 1) where k 0 ¼ x=c is the wavenumber, c being the sound velocity in the host medium, and x the angular frequency of the wave.Equation (2) yields the following observations: if b > k 0 (the propagating wave inside the waveguide is slower than the sound velocity v p < c), then k z is imaginary and no sound is radiated outside the structure.Conversely, if b < k 0 (v p > c, fast-wave), k z is real and the wave leaks out.Thus, a slow-wave, characterized by b > k 0 , is a guided wave, whereas a fast-wave, characterized by b < k 0 , is a leakywave.The region of the dispersion diagram where the condition b < k 0 is satisfied is called the radiation region, and any waveguiding structure presenting a dispersion curve xðbÞ inside the radiation region is considered a leaky-wave structure.Using Eq. ( 2), the propagation constant b is linked to the angle h 0 of the main lobe by For a dispersive medium, the quantity b=k 0 depends on frequency and therefore the main beam angle can vary as a function of frequency, this property is called frequency scanning.

B. Radiation pattern
For a 1D leaky-wave with continuous leakage along the structure (Fig. 1 with w ¼ d), the radiated sound pressure towards the upper half space reads 22 p 0 e Àcy e Àjk 0 r 2pr dy; where r is the distance between the observation point and a position y along the LWA.In the far field, distance r can be approximated by r ' r 0 À y sin h (in the phase term e Àjk 0 r ), and by r ' r 0 for the geometric attenuation 1=2pr, where r 0 is the distance between the origin (y; z ¼ 0) and the observation point.After Eq. ( 4), the radiated sound pressure reads Using Eq. ( 5) the acoustic power radiated by an infinitely long leaky-wave structure (L ! 1) is proportional to Using the radiation power pattern of Eq. ( 6) the half-power beamwidth Dh of the structure can be derived from Let us now consider a periodic arrangement of open holes (Fig. 1) of width w situated between nd and nd þ w, where n ¼ 0; …; i is the index of the hole and d is the distance between two consecutive holes.The pressure radiated at an observation point ðr 0 ; hÞ is given by where N is the number of radiating holes.Comparing Eqs. ( 5) and ( 8), it appears that the arrangement of periodic holes behaves as an array of small radiating apertures.Moreover, the frequency at which the magnitude of the radiated pressure field of Eq. ( 5) and Eq. ( 8) is maximum for a given angle h verifies Eq. ( 3), provided the length of the antenna or the number (N) of open holes is sufficiently large.

C. Radiation efficiency
The radiation efficiency of the leaky-wave structure is another important parameter for the design.It can be calculated by considering that the sound power sent at the input of a lossless leaky-wave structure is conserved in addition to the transmitted, reflected and radiated power (W in ¼ W ref þW tran þ W rad ).Defining the radiation power efficiency of the leaky-wave antenna as the ratio of the radiated power to the input power, for a lossless structure with no reflection (W ref ¼ 0) the radiation efficiency is calculated as 21 where a is the real part of the complex propagation constant (c ¼ a þ jb) and accounts for the leaked power due to the radiation.The leakage factor (a) can be calculated by where W 1 and W 2 are the average power measured in two different cross sections of the leaky-wave that are separated by length l 12 .Using Eq. ( 9) and specifying a radiation efficiency of 90%, it yields Given a leakage factor (a), which depends on the design of the leaky-wave structure, Eq. ( 11) defines the necessary length to obtain a 90% radiation efficiency.Thus, the structure can be designed based on the framework introduced above to achieve a prescribed radiation efficiency.

III. CONCEPT OF SINGLE-MICROPHONE DIRECTION FINDING A. Proposed design
The proposed single-microphone direction finding concept is illustrated in Fig. 2, where a microphone is wallmounted close to one termination, the two terminations presenting anechoic conditions, with a layer of mineral wool.
The structure is composed of an acoustic TL metamaterial with subwavelength unit-cells made of small ducts of radius r waveguide and length d unit-cell , comprising a membrane (of same cross-section, and thickness t membrane ), and a transverse cylindrical open-channel connected on one side of the duct (also denoted "stubs") of radius r stub and length L stub .Each unit-cell can be modelled as series mass and parallel compliance (supporting the right-hand or positive index propagation), combined with series compliance and parallel mass (supporting left-hand or negative index propagation), as explained in Refs.14 and 17.Thus, it behaves as an acoustic bandpass filter.The transverse open ducts, which are modelled by radiation impedances, make the structure behave as a leaky-wave antenna, allowing the radiation of an acoustic wave to the ambient medium.Schematically, a leakage can be modelled as an acoustic radiation resistance, in series with the acoustic mass representing the acoustic impedance of the stub. 17For the sake of simplicity, we do not consider here mutual coupling between consecutive stubs, assuming they are distant enough.Although the radiation resistance does not affect the TL metamaterial behaviour, it physically provides the appropriate condition for power leakage along a direction of radiation that depends on the frequency of sound.Here, the leaky-wave mechanism favours a given radiation direction depending on the refractive indices assigned to each frequency of the input wave by the TL metamaterial, after Eq. ( 3).
This structure is similar to the one proposed in Ref. 14. Here, however, flanged cylindrical ducts are substituted for the axisymmetric open channels. 16,17This modification has been mainly proposed to prevent the monopolar behaviour of the axisymmetric structure, which is expected to present an omnidirectional radiation pattern.Indeed, the radiation pattern of such an axisymmetric leaky-wave antenna is "doughnut-shaped," with constant directivity around the FIG. 2. 2D schematic representation of the single-microphone direction finding concept, based on the periodic structure inspired by Refs.13 and 14.A microphone is inserted at the left termination of the waveguide, and absorbing material (mineral wool) fills the two terminations.symmetry axis of the structure.The proposed geometrical modification of the new LWA structure, with open channels facing only the positive z direction, is then intended to sense sound waves only from one half-space.Moreover, the LWA structure is flanged to on an acoustically hard surface, yielding two effects.First, the flanged hard surface suppresses the back lobes, thus preventing the structure from radiating to the back side.Second, it increases the radiated power towards the upper half-space.

Presentation of the numerical model
Based on the design procedure explained in Refs.14 and 17, the structure of Fig. 2 with 10 unit-cells and an approximate length of 2k is designed.Note that the input and output cross sections of the waveguide are designed smaller than the main duct to ensure perfect impedance matching with the leaky-wave structure.Moreover, the LWA presents a periodic structure, except for the first and last stubs.Indeed, in order to comply with the unit-cell definition, the first and the last transverse ducts should be geometrically half of the other stubs. 14evertheless, in order to improve impedance matching, the first and last stubs' dimensions have been modified, leading to the final dimensions presented in Table I (the material properties are also given in Table II).The finite-element simulation (FEM) of the structure is carried out using COMSOL MULTIPHYSICS where the acoustic-shell interaction physics module in frequency domain is used for simulating the full structure.The open end of the stubs are flanged to an acoustic hard surface and the full structure, as well as the acoustic domain represented by a halfsphere are bounded by a perfectly matched layer (PML) to realize non-reflecting boundaries.An acoustic background pressure field is applied to the acoustic domain, and the sound pressure computed at a given position near the input of the leaky-wave is the output of the model.The FEM results are then compared to the analytical model presented in Sec.II, where the value of c is processed from the TL model. 14,17

Radiating mode
First used in the transmitting mode, a small sound source is substituted to the microphone of Fig. 2. The farfield radiation pattern of the structure is shown in Fig. 3 for monochromatic sound sources with three different frequencies: (a-1), (b-i) 900 Hz; (a-ii), (b-ii) 1000 Hz; and (a-iii), (b-iii) 1100 Hz, and it can be seen that the direction of radiation actually varies with frequency.Figure 3(a) gives a polar representation in the yz-plane, which compares the radiation patterns derived using FEM, theoretical model of Eq. ( 8) (the value of c is retrieved using TL method) 17 and experimental data of Sec.IV.These figures confirm that all three methods qualitatively agree, although small discrepancies between FEM and analytical method can be seen.Indeed, the theoretical relation of Eq. ( 8) is based on 1D model, whereas the FEM considers a 3D geometry which might result in disagreement between actual distribution of pressure field inside the waveguide and in the analytical model of Eq. ( 1), and consequently, discrepancies in the far-field radiation pattern.Figure 3(b) represents the 3D radiation pattern of the structure.It shows that, by increasing the frequency of the sound source, the radiation pattern actually scans the upper hemisphere from backward to forward.The amplitude of the radiated sound pressure is maximal at a given angle, monotonically varying with frequency, according to Eq. ( 3).Using Eqs. ( 2) and ( 3) and replacing the denominator of Eq. ( 7) by k 0 cos h 0 ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi k 2 0 À b 2 q yields the relation between half-power beamwidth (a 1/directivity) and frequency for LWA of given length: For a structure with constant Leaky factor (a), the directivity depends on ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi k 2 0 À b 2 q behaviour which tends to zero when k 2 0 % b 2 (the boundaries of the leaky region) and becomes maximum when b % 0 (see Ref. 17).Then, starting from backward the directivity first increases with frequency, having the maximum directivity near broadside, and decreases in the forward direction.Moreover, substituting Eq. ( 11) in Eq. ( 7), the directivity of the acoustic leaky-wave antenna can be derived as a function of length L, Here, the directivity increases by increasing the length of the leaky-wave structure

Receiving mode
Let us now consider the leaky-wave structure when used in the receiving mode, as presented in Fig. 2. A plane wave sound field with varying directions h in the yz plane is considered as the input of the model, and the sound power q m 1420 kg/m 3 spectrum sensed at the microphone position is the output of the model for each angular direction h.The normalized sound power spectra are derived from the FEM output data and illustrated in Fig. 4(a) for three different angles of incidence.These results are also compared to the ones processed with the theoretical model of Eq. ( 8).Based on reciprocity Eq. ( 8) holds for both transmitting and receiving modes with w ¼ 2r stub ; d ¼ d unitÀcell , and the value of c retrieved using the TL method. 17As discussed before, the small discrepancies can be justified by the differences between the 1D configuration used in the analytical formulations and the actual 3D configuration of the FEM.Moreover, the value of c retrieved with the TL modeling might also result in small inaccuracies, which may be another source of errors.It can be seen that, for each incident angle, the power spectra sensed at the microphone position presents a lobe with maximal amplitude at a given frequency and this frequency varies with the angle h according to Eq. ( 3).Based on the received power spectra of the LWA, the curves of Fig. 4(b) are derived and depict the dependence of frequency of the maximum amplitude of received power [denoted max(PWR) in the following] with respect to incident angle for FEM and theoretical data.This confirms the functionality of the LWA as a single microphone sound source localizer.It is noticeable that the received power spectra present, at angular directions near broadside, a double-bump [see the curve corresponding to the 0 direction in Fig. 4(a)].This phenomenon is due to the occurrence of a bandgap which is the consequence of unbalanced series and parallel resonances in  8).(b) Frequency-direction curve is presented, every 10 , for FEM and analytical method.Two different strategies are used to derive these data: frequency of the maximum received power amplitude [max(PWR)] versus incident angle and SC À3dB of the received power versus incident angle.
the TL metamaterial design.Indeed, the geometrical dimensions and material properties have been chosen so as to be feasible in practice with readily available materials, which makes the unbalanced condition unavoidable.To overcome the potential ambiguity in the direction estimation, resulting from the secondary lobes observed on the power spectra, a new method for processing the power spectra is proposed.As the antenna behaves as a selective bandpass filter, the spectral centroid (SC À3 dB ) of the received power can be substituted for the estimation of the frequency at which the maximum power is received.For each angular position, the frequency boundaries (f À3 dB;min ; f À3 dB;max ) are first determined as the frequencies for which half of the maximum power is received at the sensing position.Then the 3 dB spectral centroid is computed as This technique may be very useful in presence of noisy experimental data as in Sec.IV.Therefore, we propose to validate this method by comparing the direction vs frequency functions obtained with the FEM and analytical models as illustrated in Fig. 4(b).
The presented numerical and analytical results confirm that the designed structure is capable of steering towards the direction of the incoming sound signal using only a single microphone.This interesting method may challenge existing techniques relying on multiple sensors beamforming for sound source localisation.It may also be compared to interference-type (shotgun) microphones, which use interferences throughout a conventional cylindrical waveguide to increase the directivity of the sensor. 23In Sec.IV, this result will serve as a reference and compared with the frequencydirection mapping achieved with an experimental prototype.

IV. EXPERIMENTAL VALIDATION A. Experimental prototype
An experimental prototype following the design presented in Fig. 2 is fabricated by stacking 10 unit-cells made of aluminum, as illustrated on Fig. 5(a).To assemble the full structure, each membrane (DuPont Kapton FPC) is glued on a 3D printed support, and squeezed between two adjacent unit-cells using four screws as in Fig. 5(b).The 3D printed support in Fig. 5(c) is designed to decouple the stress applied on the membrane by the adjacent unit-cells.Instead, uniform gluing of the membrane over the support ensures uniform stress over the circumference of each membrane, and almost uniform clamping condition can be achieved for all membranes.This method is more practical and straightforward compared to the one relying on individual alignment of resonant frequency of each membrane proposed in Ref. 17.To complete the assembly of the prototype, the two extremities of the fabricated structure are filled with layers of mineral wool to simulate anechoic boundary conditions, and a PCB Piezotronics 130D20 ICP microphone is flush-mounted at one termination of the structure, and used as the single sensor of the LWA.Finally, it is flanged to a wooden panel of dimensions 50 cm Â 50 cm, which mimics the hard boundary conditions considered in the FEM model, while suppressing the back lobes.

B. Experimental setup
To measure the performance of the fabricated prototype as a single-microphone direction finding, the LWA mounted on a panel is installed on a turntable (Bruel & Kjaer type 5960, remote-controlled through a controller type 5949) in an anechoic chamber, and a loudspeaker (sound source) is installed at 4 m from the center of the leaky-wave, as illustrated in Fig. 6.The sound generation and the microphone signal acquisition are operated with a Bruel & Kjaer type 3160 Pulse multichannel analyzer.The FFT analysis is set on 3200 points within the [600-1400 Hz] frequency range, thus a frequency resolution of 0.25 Hz, and the signal generator delivers a sweep sine signal within the frequency range [600-1400 Hz] to the loudspeaker.

C. Measurement results
For each angular position of the turntable, the power spectrum P(f) of the received signal is processed.Figure 7(a) presents the normalized power spectra measured for different orientations of the turntable.As explained in Sec.III B, the main lobe observed on the power spectrum measurement presents a frequency, that can then be mapped to the orientation of the turntable.It is noticeable that, for each orientation, the antenna presents a relatively narrow beam width, except around 0 where two lobes are visible.This can be explained by a misalignment of the elements of the LWA (membranes and stubs), resulting in the occurrence of a band-gap around the transition frequency and emergence of second lobe in the received power spectra.

D. Frequency-source direction mapping
With a view to assessing the single-microphone direction finding functionality of the LWA, the spectral centroid and max(PWR) are derived from the power amplitude received by the ICP microphone for each orientation h of the turntable, and compared to the theoretical value obtained from Eq. (3). Figure 7(b) presents the frequency-direction curve of the LWA which confirms the good agreement between theoretical and measured data.Moreover, in order to overcome the ambiguity resulting from the emergence of a second lobe near the broadside direction, the SC À3 dB is derived from experimental data using Eq. ( 14) and is plotted against the angular position of the turntable.As expected, the SC À3 dB derived from the measured data shows smoother variations compared to the max(PWR), and is in closer agreement with the theoretical curve.

E. Discussion
The results of Figs.4(b) and 7(b) show that the theoretical, FEM and experimental data, present monotonicallyincreasing directions h as a function of frequency f, which can then be used for single-microphone direction finding applications.Furthermore, the frequency-direction mapping obtained with experimental and simulation data follows the theoretical trend, which also validates the theoretical formulas.The small discrepancies observed between measured data and the theoretical trend can be explained by two reasons: the first reason might be the differences in the fabricated prototype compared to the optimal design parameters considered in the theoretical formulation, especially the difficulty to ensure perfectly uniform clamped boundary conditions on the membranes in the experimental prototype.The second reason might be the differences between the simulation and experimental conditions.Indeed, the imperfect anechoic conditions in the waveguide extremities, as well as the likely imperfect frequency response of the loudspeaker may significantly differ from the simulation condition.Due to the listed problems and occurrence of second lobe near the broadside [particularly visible in the measured results in Fig. 7(a)], a statistical method of spectral centroid of the data is used to get smoother frequency-direction curves.

V. CONCLUSION
This study introduces a single-microphone direction finder as an interesting novel application of acoustic TL metamaterials.It consists in a periodic structure composed of vibrating membranes along a waveguide, and side open channels allowing sound leakage along the structure.In the receiving mode, it is shown that such device can be used as a single-microphone direction finder, without resorting to multiple-sensors beam-steering techniques.The concept is validated with an experimental prototype, the sound pressure  sensed by the single microphone showing an angledependent power spectrum, highlighting a monotonically increasing relationship between the angle of observation and the frequency.The frequency-direction mapping, also derived using spectral centroid (SC À3 dB ) of the received power, shows a smoother monotonic behaviour.It should also be emphasized that those results have been confirmed with several measurements, after mounting/dismantling of the prototypes, highlighting simplicity of the implementation and robustness of the LWA structure.Such device can be seen as an interesting alternative to the state-of-the-art techniques for sound source localization, and may find direct applications to sonar.Moreover, since acoustic and electromagnetic leaky-wave antenna are based on similar physics, they share the same waveguiding structure as a host medium.This suggests that a conventional (RH) electromagnetic leaky-wave radiating in the forward direction can turn into an acoustic leaky-wave antenna.However, achieving electromagnetic/acoustic leaky-wave antenna with both RH/LH behaviour radiating from backward to forward is a remaining challenge.This article may encourage future development with a view to merge acoustic/electromagnetic RH/LH leaky-wave antennas within the same structure, and proposing a new generation of mixed radar-sonar concepts based on leaky-wave concept.

FIG. 1 .
FIG. 1. Schematic representation of a leaky-wave structure.A 1D waveguide of length L along the y axis with openings of width w and a spatial spacing d.The solid and dashed lines above the structure represent the sound rays and wavefronts, respectively.

FIG. 3 .
FIG. 3. (Color online) (a) The polar representation of the far-field radiation pattern of the structure in the yz-plane for three distinct frequencies: (a-i) 900 Hz, (a-ii) 1000 Hz, and (a-iii) 1100 Hz.Each radiation pattern has been normalized to its maximum and then levelled to 8 dB.(b) Maximum-normalized 3D farfield radiation patterns of the structure, extracted from COMSOL simulations, provide information about the directivity of the structure as an acoustic antenna: (b-i) 900 Hz, (b-ii) 1000 Hz, and (b-iii) 1100 Hz.

FIG. 5 .
FIG. 5. (Color online) (a) Picture of the fabricated prototype: the holes on the top side correspond to the open channels (stubs), the ICP microphone used as the LWA single sensor being flush-mounted at one termination, on the opposite side (not visible here).(b) Picture of a disassembled LWA prototype, presenting a membrane glued to a support which is sandwiched between two unit-cells.(c) Picture of the membrane support as well as its schematic representation and geometrical dimensions.

FIG. 6 .
FIG. 6. (Color online) (a) Measurement setup in the anechoic chamber: the leaky-wave antenna is wall-mounted to a wooden panel, suppressing the back lobe and acting as a support, the whole being fixed on a turntable in front of a loudspeaker used as a sound source.(b) Schematic description of the experimental setup.

FIG. 7 .
FIG. 7. (Color online) (a) Measured power spectra sensed by the microphone in the LWA for different angular positions of the turntable h.(The experimental data are shown every five samples).Experimental results are compared to theoretical values of Eq. (8).(b) Frequency-direction curve of measured data is presented every 10 derived by two different methods: frequency of the maximum received power amplitude versus incident angle and received power SC À3dB versus incident angle.The analytical expected curve is also presented for comparison.

TABLE I .
Dimensions of the designed structure.