This work presents a sensitivity study on the use of an optical feedback interferometer to measure acoustic pressure from plane waves. The sensitivity is established by linearising the interferometer's governing equations. It is shown to be independent of the acoustic wave frequency but dependent on configuration parameters such as the optical feedback parameter or the length of the laser through which the acoustic wave passes. Experimental validation is carried out using three acoustic waveguides in the 0.5–18 kHz range. The sensitivity obtained enables broadband acoustic pressure measure with a low mean relative error in comparison with a reference condenser microphone.

Acousto-optic methods for the detection and the characterization of acoustic waves have been extensively developed and studied.1–4 They rely on the interferometric measurement of optical index variation caused by the acoustic waves: the acousto-optic effect. Recently, an optical feedback interferometer (OFI), also referred as a self-mixing interferometer, has been proposed for acoustic pressure field visualization.5,6 This interferometric technique relies on the use of a semiconductor laser that is sensitive to light feedback effects. Indeed, if a fraction of the light emitted by a laser diode returns in its cavity by reflection or scattering, its behavior changes in terms of emitted power and optical wavelength according to the optical path followed by the laser beam.7 In practice, it is implemented by using a laser diode which targets a retro-reflective surface and interferences are measured with an embedded photodiode. Thus, the OFI has the advantage of requiring little equipment and being self-aligned compared to other types of interferometers. As well as other optical methods, the OFI performs an integrated measurement along the laser path. Therefore, knowledge of the radiation from the acoustic source is required to solve this integral.2,3 Only a few studies have been conducted on the performance of the OFI for quantitative acoustic pressure measurement. In addition, the detection specifications in terms of amplitude, signal-to-noise ratio, and frequency content of the acoustic source are not clearly established for acoustic measurement with an OFI. Knudsen et al. studied the lower detection limit for 3 kHz plane sinusoidal waves in a waveguide and observed a proportionality between the OFI output signal and the sound pressure amplitude.8 However, in the laser self-mixing theory7 this proportionality is not highlighted by the equations and remains to be established. Also, the OFI sensitivity is supposed to be independent of the acoustic wave frequency, which is to be measured.

The purpose of this study is to characterize the sensitivity of the OFI to the acoustic pressure in a waveguide. A model of the OFI sensitivity for the measurement of acoustic plane waves is established. Measurements in acoustic waveguides validate this model. The paper is organized as follows: the equations that model the optical feedback and the response of an OFI to plane acoustic waves are detailed in Sec. 2. Then, the experimental apparatus is described in Sec. 3, it aims at measuring plane waves with an OFI and a microphone. The results are presented and discussed in Sec. 4. Finally, the conclusion and future works are described in Sec. 5.

In this section, a proportionality relationship is established between the acoustic pressure of a plane wave and the output signal of an OFI.

When an acoustic wave propagates through the air, the variation of acoustic pressure induces a local variation of the refractive index. The pressure p ( r , t ), at a given time t and position r can be written as
(1)
where p0 is the mean pressure of the air and p ( r , t ) is the acoustic pressure. Similarly, the optical refractive index n ( r , t ) can be written as follows:
(2)
where n0 is the mean refractive index and n ( r , t ) is the fluctuation of refractive index caused by the acoustic waves which is called the acousto-optic effect. The latter can be calculated using Ciddor's model9 as proportional to the acoustic pressure,
(3)
In Eq. (3), β is a quantity depending on: λ0 the electromagnetic wavelength, p0 the static pressure, T0 the temperature of the air, ϕ h the humidity rate, and c CO 2 the molar concentration of CO 2. It is then possible to calculate the acoustic pressure by measuring the refractive index of the medium by using an optical interferometer and Eq. (3). For laboratory conditions of λ 0 = 1309 nm, p 0 = 1 bar, T 0 = 20 °C, ϕ h = 0.5, and c CO 2 = 440 ppm, the value of β is 2.6 × 10 9 Pa–1.
As a laser beam goes through the medium in which acoustic waves propagate and reflects back from a surface into the laser cavity, the optical path L along the beam is altered due to changes in the optical index of the medium. Thus, the optical path L is expressed as follows:
(4)
where x is the coordinate along the laser beam, n(x, t) is the optical index at position x, and L(t) is the distance between the laser cavity and the reflecting surface. The latter is defined as L ( t ) = L 0 + L V ( t ), where L0 is the laser beam length at rest and L V ( t ) is the geometric length variation of the laser beam due to possible displacements of the reflecting surface. It is assumed that the laser beam is not deflected. By combining Eqs. (2), (3), and (4), the optical path L is rewritten as
(5)
In Eq. (5), L AO ( t ) is the apparent length variation of the laser beam caused by the acousto-optic effect.8 For the sake of simplicity, the time variable t is dropped in the remaining equations.
An OFI uses the self-mixing properties of laser diodes: if a fraction of the photons emitted by the laser diode returns into the laser cavity, the laser diode power and wavelength vary.7 This phenomenon is known as optical feedback effect. If the amount of photons returning to the laser cavity is much smaller than the amount leaving it, which in practice can be obtained when the laser beam is backscattered from rough surfaces,10 this phenomenon can be modeled as7 
(6)
In Eq. (6), L is the optical path outside the laser's cavity, λ0 is the wavelength of the laser diode without optical feedback, λ is the actual wavelength of the laser with feedback, C is the feedback parameter that depends on the amount of light reflected back into the laser cavity, and α is the linewidth enhancement factor.11 
In the case where L V is negligible with respect to L AO in Eq. (5) and based on numerical simulations of Eq. (6),10 the wavelength λ can be described as
(7)
where λ L 0 is the value of λ where L = L0 and γ is a coefficient depending on C and α. Note that λ L 0 γ L AO for acoustic measurements which in practice means that the variation of the optical wavelength caused by L AO is negligible compared to the value of λ L 0.
The effect of the optical feedback on the laser power is modeled by the equation12 
(8)
where P is the power of the laser diode, P 0 is the power of the laser diode without optical feedback, and m is the modulation index which is proportional to C. In a practical setup, the power of the laser diode P is measured by using an embedded photodiode. Its current is converted into voltage by using a transimpedance amplifier, resulting in the OFI signal. One denotes u OFI the alternative component (AC) of this signal.

For sinusoidal plane waves, it has been observed experimentally that u OFI is proportional to L AO.8 However, there is no proportionality between P and L AO when replacing L and λ in Eq. (8) by their expressions in Eqs. (5) and (7), respectively. Nevertheless, as detailed in the supplementary material one can obtain a proportional relationship under the following assumptions:

  • cos ( 4 π n 0 L AO / λ ) 1, which can be satisfied if L AO λ.

  • cos ( 4 π n 0 L 0 γ L AO / ( λ L 0 ) 2 ) 1, which can be satisfied if L AO ( λ L 0 ) 2 / ( L 0 γ ).

Using a first order Taylor series expansion, the relationship between u OFI and L AO is expressed as
(9)
where U0 is the direct component (DC) component of the OFI signal without optical feedback.
When an acoustic plane wave propagates perpendicularly to the laser beam, the latter lies along a wavefront whose pressure is constant. L AO can then be described as
(10)
where Π is the rectangular function, p is the acoustic pressure of the plane wave along the laser beam, L1 is the part of the laser beam that is disturbed by the plane wave centered in x = x0. Thus, u OFI can be written as
(11)
where the OFI sensitivity κ OFI (in V/Pa) is
(12)
According to this formula, the sensitivity of the OFI depends on L1 but not on the frequency of the acoustic waves.

The aim here is to capture acoustic waves in a section of a waveguide using a microphone and an OFI. One compares these two measurements to verify a proportionality relationship by linear regression and deduce κ OFI. A scheme of the setup is shown in Fig. 1.

Fig. 1.

(a) Scheme of the experimental setup. The OFI laser beam emitted by the laser diode (LD) crosses a section of the acoustic waveguide and points at the reflective tape. AMP is the power amplifier of the signal emitted by the white noise generator. CON is the microphone conditioner. TIA is the transimpedance amplifier that converts the current from the photodiode (PD) into a voltage which is the OFI signal. (b) Scheme of the waveguide section crossed by the laser beam.

Fig. 1.

(a) Scheme of the experimental setup. The OFI laser beam emitted by the laser diode (LD) crosses a section of the acoustic waveguide and points at the reflective tape. AMP is the power amplifier of the signal emitted by the white noise generator. CON is the microphone conditioner. TIA is the transimpedance amplifier that converts the current from the photodiode (PD) into a voltage which is the OFI signal. (b) Scheme of the waveguide section crossed by the laser beam.

Close modal

The laser diode (LD) of the OFI is a Thorlabs© L1310P5DFB with a maximum power of 5 mW and a wavelength λ 0 = 1309 nm. It is embedded with a photodiode (PD), which delivers a current that is proportional to the power P ( t ) of LD. The beam is collimated using a Thorlabs© C110TMD-C lens with a focal length of 6.24 mm. The LD is powered by a current driver Thorlabs© LDC205C and maintained at a constant temperature of 12 °C using a Thorlabs© TED200C temperature controller. The current from the PD is converted into a voltage u OFI ( t ) with a transimpedance amplifier Femto© DLPCA-200 whose gain is set to 104 V/A.

The laser light spot of the OFI is positioned on a rough reflective tape placed at L 0 40 cm from the LD. It allows the return of a certain quantity of photons in its cavity by backscattering so that Eq. (6) is valid.10 The reflective tape is glued on an accelerometer PCBpiezotronics© 352C65 linked to a PCBpiezotronics© 482C05 conditioner which delivers a voltage u acc ( t ) that is proportional to the sensor acceleration. It allows the measurement of mechanical vibrations L V ( t ) in Eq. (5).

The OFI feedback parameter C is adjusted by modifying the LD current. This parameter is related to the amount of light backscattered into the laser.7 It strongly depends on the alignment of the laser beam and L0. Adjusting the current, and thus the light intensity emitted by the laser diode, allows a finer and more controllable adjustment of the C parameter without changing the optical setup of the OFI.

Three rectangular section acoustic waveguides were used in the frequency range below their theoretical cut-off frequency (TCOF) to satisfy the plane wave hypothesis that was assumed in Eq. (10).13 The dimensions and acoustic properties of the waveguides are given in Table 1 with TCOF = c 0 / ( 2 L 1 ) (c0 is the speed of sound). The laser beam passes through the guide by two side holes. A reference microphone is flush-mounted above the beam in order to measure the acoustic pressure in the same section [see Fig. 1(b)]. The reference microphone is a 1/4″ GRAS© 40BH (waveguides Nos. 1 and 2) or a 1/8″ GRAS© 40DP (waveguide No. 3) connected to a Brüel & Kjaer© NEXUS 2690-A-0F2 conditioner. The latter delivers a voltage u mic ( t ) proportional to the acoustic pressure p ( t ) of the section. Each waveguide is excited by a different loudspeaker at one of its extremities. They are powered by a Visaton© AMP 2.2 LN amplifier. Hereinafter, the “acoustic source” refers to the amplifier, loudspeaker, and waveguide assembly.

Table 1.

Properties of the waveguides. See Fig. 1 for the nomenclature of the waveguide dimensions.

Waveguide No. Dimensions L 1 × L 2 × L 3 (mm) Waveguide end TCOF (Hz)
45 × 430 × 25  Closed  3 500 
35 × 350 × 35  Open  4 000 
10 × 200 × 10  Open  18 000 
Waveguide No. Dimensions L 1 × L 2 × L 3 (mm) Waveguide end TCOF (Hz)
45 × 430 × 25  Closed  3 500 
35 × 350 × 35  Open  4 000 
10 × 200 × 10  Open  18 000 

The excitation signal is a white noise generated in the range [0,20] kHz by a Siglab DSP 20–42 Dynamic Signal Analyzer controlled by a matlab program. The root mean square (RMS) value of the excitation signal is denoted u RMS. For each waveguide, the acoustic waves are measured sequentially by the microphone and the OFI. For the latter, an additional background noise measurement, denotes u noise ( t ), is done with the loudspeaker turned off. Then, the amplitude spectral densities (ASD) are calculated by analyzing 1000 signal frames of 200 ms duration. One denotes u ̃ mic ( f ) , u ̃ OFI ( f ) , u ̃ acc ( f ), and u ̃ noise ( f ) the ASD of u mic ( t ) , u OFI ( t ) , u acc ( t ), and u noise ( t ), respectively.

The first step is to identify the frequency range in which the contribution of vibrations L V can be neglected (Sec. 4.1). Then, the experimental method to estimate κ OFI is detailed and the results are discussed (Sec. 4.2).

As the sensitivity of the OFI is not measured yet, it is not possible to separate the contributions of L V and L AO in the OFI signal [see Eqs. (5) and (8)]. However, they can be estimated from the ASDs of the accelerometer and the microphone voltage
(13)
where L ̃ V ( f ) is the ASD of L V ( t ) and κ acc is the sensitivity of the accelerometer [in V/(m/s2)] and
(14)
where L ̃ AO ( f ) is the ASD of L AO ( t ) , p ̃ mic ( f ) = κ mic u ̃ mic ( f ) is the estimated ASD of p ( t ) by the microphone, and κ mic is the proportionality coefficient between u ̃ mic ( f ) and p ̃ mic ( f ) (in V/Pa). Note that Eq. (14) is only valid for plane waves. Measurements of L ̃ V ( f ) and L ̃ AO ( f ) for each waveguide are plotted in Fig. 2, showing their resonance frequencies.
Fig. 2.

Measurements of L ̃ V ( f ) (in solid blue) and L ̃ AO ( f ) (in dashed red) for (a) waveguide No. 1, (b) waveguide No. 2, and (c) waveguide No. 3 below their TCOF. u RMS = 0.5 V.

Fig. 2.

Measurements of L ̃ V ( f ) (in solid blue) and L ̃ AO ( f ) (in dashed red) for (a) waveguide No. 1, (b) waveguide No. 2, and (c) waveguide No. 3 below their TCOF. u RMS = 0.5 V.

Close modal

Beyond 500 Hz for the waveguides Nos. 1 and 2, and 1 kHz for the waveguide No. 3, L ̃ V ( f ) is for most frequencies at least 10 times lower than L ̃ AO ( f ). It allows one to consider that above these threshold frequencies, L ̃ V ( f ) is negligible compared to L ̃ AO ( f ) and that the variations in the OFI signal are mainly caused by the acoustic waves in the waveguide.

Measurements are made for several excitation levels u RMS, in order to have different values of p ̃ mic ( f ). The OFI feedback parameter C is set to its maximum value for the experimental setup. In Fig. 3, u ̃ OFI ( f ) is plotted as a function of p ̃ mic ( f ) for each frequency which respects
(15)
Fig. 3.

ASD of the OFI signal u ̃ OFI ( f ) as a function of p ̃ mic ( f ) in (a) waveguide No. 1, (b) waveguide No. 2, and (c) waveguide No. 3 for frequencies satisfying Eq. (15). Green circles are measurements made with u RMS = 0.5 V, red triangles with u RMS = 0.3 V, and blue crosses with u RMS = 0.1 V. The acoustic level (in dBRMS) of each measurement is given in the legend. The black dashed lines are linear regressions of slope κ OFI.

Fig. 3.

ASD of the OFI signal u ̃ OFI ( f ) as a function of p ̃ mic ( f ) in (a) waveguide No. 1, (b) waveguide No. 2, and (c) waveguide No. 3 for frequencies satisfying Eq. (15). Green circles are measurements made with u RMS = 0.5 V, red triangles with u RMS = 0.3 V, and blue crosses with u RMS = 0.1 V. The acoustic level (in dBRMS) of each measurement is given in the legend. The black dashed lines are linear regressions of slope κ OFI.

Close modal

The later condition, where the factor 2 is set arbitrarily, ensures that the points used for the linear regressions are sufficiently above the background noise u ̃ noise ( f ).8 The OFI sensitivity is then estimated by performing a linear regression on the remaining data. According to Eq. (11), the slope of the linear regression is κ OFI whose values for each waveguide are shown in Table 2. The linearization assumptions allowing the validity of Eq. (11) reach at least 0.94 for assumption (i) and 0.89 for assumption (ii).

Table 2.

Measured sensibility and MRE for each waveguide with the maximum value of C for this setup.

Waveguide No. 1 2 3
κ OFI (V/Pa)  6.69 × 10 6  5.71 × 10 6  1.71 × 10 6 
κ OFI / L 1 [V/(Pa m)]  1.48 × 10 4  1.63 × 10 4  1.71 × 10 4 
MRE (dB)  −6.5  −9.0  −4.6 
Waveguide No. 1 2 3
κ OFI (V/Pa)  6.69 × 10 6  5.71 × 10 6  1.71 × 10 6 
κ OFI / L 1 [V/(Pa m)]  1.48 × 10 4  1.63 × 10 4  1.71 × 10 4 
MRE (dB)  −6.5  −9.0  −4.6 
Using these κ OFI estimations, the acoustic pressure estimated by the OFI p ̃ OFI ( f ) is calculated by dividing u ̃ OFI ( f ) by κ OFI. Figure 4 compares the values of p ̃ mic ( f ) and p ̃ OFI ( f ) for each waveguide. Acoustic pressure estimations made with the microphone and the OFI are quite similar. For each waveguide, the mean relative error (MRE) over the N frequencies fi which respect Eq. (15) is shown in Table 2. Is is computed as follows:
(16)
Fig. 4.

Acoustic pressure estimated from the microphone p ̃ mic ( f ) (in dashed red) and the OFI p ̃ OFI ( f ) (in solid blue) as a function of the frequency of the acoustic waves for (a) waveguide No. 1, (b) waveguide No. 2, and (c) waveguide No. 3. p ̃ OFI ( f ) values which are not in the gray area are used for κ OFI estimation. u RMS = 0.5 V.

Fig. 4.

Acoustic pressure estimated from the microphone p ̃ mic ( f ) (in dashed red) and the OFI p ̃ OFI ( f ) (in solid blue) as a function of the frequency of the acoustic waves for (a) waveguide No. 1, (b) waveguide No. 2, and (c) waveguide No. 3. p ̃ OFI ( f ) values which are not in the gray area are used for κ OFI estimation. u RMS = 0.5 V.

Close modal

The low MRE values obtained suggest that κ OFI is independent of the frequency. In addition, for the frequencies where p ̃ OFI ( f ) is getting closer to 2 u ̃ noise ( f ) / κ OFI (the gray curves in Fig. 4), the similarity between p ̃ OFI ( f ) and p ̃ mic ( f ) decreases. Indeed, u OFI may go below the lower detection limit.8 It is therefore recommended to exclude these frequencies when estimating κ OFI.

As κ OFI is proportional to L1 [see Eq. (12)], values of κ OFI / L 1 are also calculated to remove the influence of the waveguide geometry on the OFI sensitivity and are presented in Table 2. It is observed that values of κ OFI / L 1 are of the same order of magnitude for each waveguide, which confirms a satisfactory κ OFI estimation with this measurement method. Table 3 shows different κ OFI values obtained for different values of C, set as explained in Sec. 3.1, with waveguide No. 1. These results show a reduction in the OFI sensitivity with the value of C as expected by Eq. (12). Differences between κ OFI / L 1 values of each waveguide may be due to the experimental method, or to the limit of validity of some hypotheses. In particular, the setting of the OFI-reflector path can be slightly modified when switching the waveguides, which can slightly change the values of C and κ OFI. The validity of assumptions such as neglecting acoustic wave radiation from the waveguides side holes through which the optical beam passes would also require more attention in future studies to improve sensitivity estimation.

Table 3.

Measured sensibility for different values of C such as C 1 < C 2 < C 3 on waveguide No. 1.

C C1 C2 C3
κ OFI (V/Pa)  2.60 × 10 6  4.95 × 10 6  6.69 × 10 6 
C C1 C2 C3
κ OFI (V/Pa)  2.60 × 10 6  4.95 × 10 6  6.69 × 10 6 

A study of the optical feedback theory has demonstrated that the output signal of the OFI can be considered proportional to the acoustic pressure with two assumptions. The calculation of the sensitivity has shown that it does not depend on the frequency of the acoustic waves to be measured. The measurement of acoustic plane waves with a reference microphone and an OFI has allowed to calculate the sensitivity of the latter. The similarity between the OFI and the microphone acoustic pressure estimations in the range 0.5–18 kHz suggests an independence of the OFI sensitivity to the frequency. To verify this property on higher frequencies, it will be necessary to change the acoustic source. Future work includes the measurements of other acoustic sources to study the ability of the OFI to measure ultrasonic waves. To measure acoustic waves of larger amplitudes, the linearization assumptions can be defeated. Future work will present how to estimate p from u OFI beyond linearization conditions.

See the supplementary material for the linearization of Eq. (8).

We thank Emmanuel Jondeau and Jean-Charles Vingiano for their help in the realization of the experimental setup and Edouard Salze for his readings. The present work is part of the program MAMBO “Méthodes Avancées pour la Modélisation du Bruit moteur et aviOn” (“Advanced methods for engine and aicraft noise modelling”) coordinated by Airbus SAS. It was supported by the Direction Générale de l'Aviation Civile (DGAC) under the Grant No. 2021–50. This work was supported by the Labex CeLyA of Université de Lyon, operated by the French National Research Agency (Grant No. ANR10-LABX-0060/ANR-11-IDEX-0007).

There are no conflicts of interest related to this work.

The data that support the findings of this study are available from the corresponding author upon reasonable request.

1.
L.
Zipser
and
H. H.
Franke
, “
Refracto-vibrometry—A novel method for visualizing sound waves in transparent media
,”
J. Acoust. Soc. Am.
123
(
5
),
3314
(
2008
).
2.
P.
Yuldashev
,
M.
Karzova
,
V.
Khokhlova
,
S.
Ollivier
, and
P.
Blanc-Benon
, “
Mach-Zehnder interferometry method for acoustic shock wave measurements in air and broadband calibration of microphones
,”
J. Acoust. Soc. Am.
137
(
6
),
3314
3324
(
2015
).
3.
P.
Lecomte
,
Q.
Leclère
, and
S.
Ollivier
, “
Equivalent source model from acousto-optic measurements and application to an acoustic pulse characterization
,”
J. Sound Vib.
450
,
141
155
(
2019
).
4.
K.
Ishikawa
,
Y.
Shiraki
,
T.
Moriya
,
A.
Ishizawa
,
K.
Hitachi
, and
K.
Oguri
, “
Low-noise optical measurement of sound using midfringe locked interferometer with differential detection
,”
J. Acoust. Soc. Am.
150
(
2
),
1514
1523
(
2021
).
5.
K.
Bertling
,
J.
Perchoux
,
T.
Taimre
,
R.
Malkin
,
D.
Robert
,
A. D.
Rakić
, and
T.
Bosch
, “
Imaging of acoustic fields using optical feedback interferometry
,”
Opt. Express
22
(
24
),
30346
(
2014
).
6.
P. F.
Urgiles Ortiz
,
J.
Perchoux
,
A. L.
Arriaga
,
F.
Jayat
, and
T.
Bosch
, “
Visualization of an acoustic stationary wave by optical feedback interferometry
,”
Opt. Eng.
57
(
05
),
1
(
2018
).
7.
K.
Petermann
,
Laser Diode Modulation and Noise
(
Springer Netherlands
,
Dordrecht
,
1988
), pp.
250
290
.
8.
E.
Knudsen
,
J.
Perchoux
,
T.
Mazoyer
,
F.
Jayat
,
C.
Tronche
, and
T.
Bosch
, “
Lower detection limit of the acousto-optic effect using optical feedback interferometry
,” in
2020 IEEE International Instrumentation and Measurement Technology Conference (I2MTC)
, Dubrovnik, Croatia (
IEEE
,
Piscataway, NJ
,
2020
), pp.
1
4
.
9.
P. E.
Ciddor
, “
Refractive index of air: New equations for the visible and near infrared
,”
Appl. Opt.
35
(
9
),
1566
1573
(
1996
).
10.
R.
Kliese
,
T.
Taimre
,
A. A. A.
Bakar
,
Y. L.
Lim
,
K.
Bertling
,
M.
Nikolić
,
J.
Perchoux
,
T.
Bosch
, and
A. D.
Rakić
, “
Solving self-mixing equations for arbitrary feedback levels: A concise algorithm
,”
Appl. Opt.
53
(
17
),
3723
3736
(
2014
).
11.
C.
Henry
, “
Theory of the linewidth of semiconductor lasers
,”
IEEE J. Quantum Electron.
18
(
2
),
259
264
(
1982
).
12.
S.
Donati
,
G.
Giuliani
, and
S.
Merlo
, “
Laser diode feedback interferometer for measurement of displacements without ambiguity
,”
IEEE J. Quantum Electron.
31
(
1
),
113
119
(
1995
).
13.
D.
Blackstock
,
Fundamentals of Physical Acoustics
(
Wiley
,
New York
,
2000
), pp.
421
424
.

Supplementary Material