# Negative reflection of elastic guided waves in chaotic and random scattering media

Benoît Gérardin, Jérôme Laurent, François Legrand, Claire Prada, Alexandre Aubry

► To cite this version:

Benoît Gérardin, Jérôme Laurent, François Legrand, Claire Prada, Alexandre Aubry. Negative reflection of elastic guided waves in chaotic and random scattering media. Scientific Reports, 2019, 10.1038/s41598-019-38480-3. hal-02335187

HAL Id: hal-02335187

https://hal.science/hal-02335187

Submitted on 28 Oct 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

# Negative reflection of elastic guided waves in chaotic and random scattering media

Benoît Gérardin, Jérôme Laurent, François Legrand, Claire Prada & Alexandre Aubry ID

Received: 22 May 2018

Accepted: 27 December 2018

Published online: 14 February 2019

The propagation of waves in complex media can be harnessed either by taming the incident wave-field impinging on the medium or by forcing waves along desired paths through its careful design. These two alternative strategies have given rise to fascinating concepts such as time reversal or negative refraction. Here, we show how these two processes are intimately linked through the negative reflection phenomenon. A negative reflecting mirror converts a wave of positive phase velocity into its negative counterpart and vice versa. In this article, we experimentally demonstrate this phenomenon with elastic waves in a 2D billiard and in a disordered plate by means of laser interferometry. Despite the complexity of such configurations, the negatively reflected wave field focuses back towards the initial source location, thereby mimicking a phase conjugation operation while being a fully passive process. The super-focusing capability of negative reflection is also highlighted in a monochromatic regime. The negative reflection phenomenon is not restricted to guided elastic waves since it can occur in zero-gap systems such as photonic crystals, chiral metamaterials or graphene. Negative reflection can thus become a tool of choice for the control of waves in all fields of wave physics.

Controlling the propagation of acoustic or electromagnetic waves, is of fundamental interest for many applications ranging from imaging the living and detecting hazardous components, to information processing and structural health monitoring. In the past decades, there has been many proposals in this regards, which can be separated within two approaches. On the one hand, the wave fields can be tamed in order to take advantage of the complexity of propagation media, for instance, to focus waves or image various objects. This is realized in the temporal domain using time reversal (TR) mirrors¹⁻³ or in the spatial domain using phase conjugation (PC)⁴⁻⁶ and wave-front shaping techniques⁷ developed in optics. On the other hand, one can force waves along desired paths through a careful design of man-made materials. This can be achieved using metamaterials, an arrangement of tailored sub-wavelength building blocks from which the material gains its unusual properties⁸⁻¹⁰. The advent of such structures has given rise to fascinating concepts such as negative refraction¹¹⁻¹³, transformation optics¹⁴,¹⁵ or metasurfaces¹⁶. Although the concepts of TR and negative refraction have been developed in an independent fashion, they are intimately linked processes¹⁷,¹⁸.

Here, we want to push forward this analogy by investigating the negative reflection (NR) phenomenon. A NR mirror is an interface at which light or sound is retro-reflected. There is a strong similarity with a PC mirror. In a PC experiment, if the incident wave is divergent, the PC wave is converging [Fig. 1(a)]. It follows the same path as the incident wave but in an opposite way, thus back-focusing exactly on the source location. Let us now consider a NR mirror and an incident forward wave with a Poynting vector, **P**ᵢ, and a wave vector, **k**ᵢ, pointing in the same direction (i.e. with a positive phase velocity, v<sub>φ</sub> > 0) [see Fig. 1(b,c)]. The NR mirror gives rise to a reflected backward wave (v<sub>φ</sub> < 0), i.e. with a Poynting vector, **P**<sub>r</sub>, and a wave vector, **k**<sub>r</sub>, of opposite direction [see Fig. 1(b,c)]. By virtue of the Snell-Descartes law, the incident and reflected wave vectors are strictly identical (k<sub>r</sub> = kᵢ) but their Poynting vectors shall be in opposite directions (P<sub>r</sub> = −Pᵢ) [Fig. 1(b)]. Negative reflection is a perfectly reciprocal phenomenon: A backward incident wave with anti-parallel **k**<sub>r</sub> and **P**<sub>r</sub> will be negatively reflected into a forward mode with parallel **k**<sub>r</sub> and **P**<sub>r</sub>.

If the incident wave is divergent, the NR wave back-converges towards the initial source position. We thus recover a similar effect as in PC, although the reflected wave vectors in each case are in opposite directions [Fig. 1(a,c)]. The PC concept and its temporal equivalent, TR, have shown to be particularly powerful in

ESPCI Paris, PSL Research University, CNRS, Univ Paris Diderot, Sorbonne Paris Cité, Institut Langevin, UMR 7587, 1 rue Jussieu, F-75005, Paris, France. Correspondence and requests for materials should be addressed to A.A. (email: alexandre.aubry@espci.fr)
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Analogy between phase conjugation and negative reflection. (a) A PC mirror has the property to reverse both the propagation direction and the phase of an incident wave field. (b) A wave packet exhibiting a positive/negative phase velocity, $v_φ$, is associated with parallel/antiparallel wave and Poynting vectors, respectively. (c) A NR mirror converts an incident wave of phase velocity $v_φ$ into a wave of opposite phase velocity $-v_φ$.</figcaption></figure>

reverberating or through inhomogeneous media¹,¹⁹,²⁰. If the phase accumulated by the waves following each scattering path is reversed, the PC waves re-accumulate the same phase during their trip back due to spatial reciprocity. They finally sum up coherently at the initial source location. This has led to spectacular focusing experiments in acoustics and optics⁷. Any scattering or reverberating medium can become a lens if the incident wave is properly designed by a PC process. In this paper, we show how NR can provide similar focusing abilities in complex environments.

NR may occur in a medium that can simultaneously support forward and backward waves. Such a coexistence is possible in, at least, two types of physical systems. The first situation corresponds to elastic guided waves propagating, for instance, in isotropic²¹⁻²⁴ or anisotropic²⁵ plates, pipes²², rods²⁶, strips²⁷, bi-layer structures²⁸, or even in the Earth's crust²⁹. Such systems can actually support an ensemble of modes, the so-called Lamb waves, which exhibit complex dispersion properties³⁰. In particular, some of these Lamb waves display a zero-group velocity (ZGV) point which results from the interference between forward and backward modes²¹⁻²⁵,³¹. The second situation corresponds to zero-gap systems, which actually support two bands that cross symmetrically, without forming a band gap. One of the most notable zero-gap material is graphene³², or its photonic/phononic crystal analogues³³. Others are quarter wave stack photonic crystals³⁴, chiral metamaterials³⁵,³⁶ or transmission-line devices³⁷. NR is therefore a very general phenomenon and this work can be, in principle, transposed to all the aforementioned systems. In this article, we have chosen to consider the NR of elastic waves in plates because they provide an excellent platform to investigate negative refraction³⁸⁻⁴⁰ and reflection⁴¹⁻⁴³. Indeed, laser interferometry allows to probe the wave-field over the whole plate surface with an excellent spatial, temporal and spectral resolution. The NR phenomenon and its analogy with PC will be investigated in two complex geometries: a chaotic billiard and a disordered plate. Whereas the reflection of waves at the boundaries of a cavity or by a random distribution of scatterers usually induces a fully random speckle wave field, NR will be shown to give rise to a purely deterministic wave field: the PC replica of the incident wave. At last, we will discuss about the manifestation of NR in the monochromatic regime. In particular, we will describe how it can lead to the super-focusing of elastic waves in a cavity or to wave trapping in a scattering medium.

## Results

**Forward and backward Lamb modes.** Whereas an infinite medium supports two longitudinal and transverse modes traveling at unique velocities, plates support two infinite sets of Lamb wave modes. Lamb waves are elastic waves whose particle motion lies in the sagittal plane that contains the direction of wave propagation and the plate normal. The deformation induced by these modes can be either symmetric (referred to as $S_1$) or antisymmetric (referred to as $A_1$) with respect to the median plane. All symmetric and antisymmetric Lamb modes (except the zero-order $S_0$, $A_0$, and first-order $A_1$ modes) can exhibit a ZGV point³¹. For the sake of simplicity, but without loss of generality, we will specifically focus on the first symmetric modes.

In this article, we study the propagation of elastic waves across a Duralumin plate (aluminium alloy) of density $d = 2790 \text{ kg/m}^3$, longitudinal wave velocity $c_L = 6400 \text{ m.s}^{-1}$, and transverse wave velocity $c_T = 3120 \text{ m.s}^{-1}$. The dispersion curves of Lamb modes in this material are displayed in Fig. 2(c). The symmetric zero-order mode $S_0$ is the extensional mode of the plate. It exhibits free propagation to zero frequency, whereas the higher order modes admit a cut-off frequency. The dispersion branch $S_1$ displays several crucial features. As illustrated by Fig. 2(c), its repulsion with the $S_2$-mode gives rise to a ZGV point ($f_{ZGV}, k_{ZGV}$)⁴⁴,⁴⁵. In the case of Duralumin, $f_{ZGV} \times 2h = 2.86 \text{ MHz.mm}$ and $k_{ZGV} \times 2h/(2\pi) = 0.26$, with $2h$ the plate thickness. From the ZGV resonance
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Experimental configuration. (a) Lamb modes are generated in a Duralumin plate of thickness $2h = 1.5$ mm by a 7 mm-diameter transducer. The Duralumin plate is here a 80 mm-diameter truncated chaotic billiard. (b) Lamb modes are excited in a disordered Duralumin plate of thickness $2h = 2$ mm by a piezoelectric strip. In each case (a,b) the normal component of the plate vibration is measured with an interferometric optical probe, mounted on a 2D translation stage. (c) Dispersion curves of the first propagating symmetrical Lamb modes in a Duralumin plate computed from the Rayleigh-Lamb equation³⁰.</figcaption></figure>

frequency $f_{ZGV}$ to the cut-off frequency, there is a coexistence between a backward mode, referred to as $S_{2b}$ in the literature, and the forward mode $S_1^{22}$.

In general, at a plate edge, a single Lamb mode is reflected into a large number of modes that can be either propagative, inhomogeneous or evanescent. The key phenomenon occurring with the $S_1$ Lamb mode in the vicinity of the ZGV point is that it is mainly reflected into the $S_{2b}$ mode. An analytical study of this phenomenon⁴³ has shown that the reflection coefficient of the $S_1$ mode into the $S_{2b}$ mode is close to -1, while the conversion into the other modes can be neglected. Furthermore, this peculiar property holds over a large angular spectrum close to the ZGV-point. This phenomenon arises from the equality of $S_1$ and $S_{2b}$ wave numbers at the ZGV-point. This implies that these two modes are associated with similar stress-displacement fields and only differ by their opposite Poynting vectors. As a consequence, the stress-free boundary condition at the edge of the plate can be satisfied with a simple combination of $S_1$ and $S_{2b}$ modes stress field, leading to a reflection coefficient $\rho = -1$. As a result, a free edge acts as a nearly perfect NR mirror for the $S_1$ and $S_{2b}$ modes⁴³. In this paper, we will now show how this peculiar property holds in much more complex geometries and how negative reflection can be taken advantage of to harness waves in chaotic or scattering media.

**Negative reflection in chaotic cavity.** This striking effect is first investigated in an elastic chaotic billiard [see Fig. 2(a)]: a truncated circular plate of diameter 80 mm and thickness $2h = 1.5$ mm. The ZGV-point occurs at $f_{ZGV} = 1.91$ MHz and $k_{ZGV} = 1.06 \text{ mm}^{-1}$. The ultrasound source consists in a 7 mm-diameter piezoelectric transducer placed on top of the plate. The transducer's width is larger than the ZGV wavelength ($\lambda_{ZGV} \sim 6$ mm) so that a dominant excitation of the $S_{2b}$-mode is provided. The out-of-plane component of the induced wave-field is recorded in the time domain over 1 ms by means of laser interferometry (see Methods). A discrete Fourier transform of the recorded signals is first performed over the temporal domain. The $S_{2b}$ and $S_1$ contributions are then separated with low and high pass spatial filters, respectively, with a cut-off at $k_{ZGV}$ (see Supplementary Information S1). Figure 3(a) shows the incident $S_{2b}$ wave field in the vicinity of the ZGV frequency ($f = 1.91$ MHz). Not surprisingly, it displays divergent cylindrical wave fronts centered around the transducer location. Figure 3(b) and (c) show the reflected $S_1$ wave field at $f = 1.92$ MHz and $f = 1.91$ MHz, respectively. The first frequency is well above the ZGV-point displayed in Fig. 2(c). There is a clear mismatch between the $S_1$ and $S_{2b}$ wave numbers: $k_1 = 1.26 \text{ mm}^{-1}$ and $k_{2b} = 0.86 \text{ mm}^{-1}$, respectively (see Supplementary Information S1). As dictated by Snell's law, the reflection angle is, in this case, not opposite to the incident angle [Fig. 3(e)]. This explains the random feature displayed by the reflected wave field in Fig. 3(b), which is characteristic of a chaotic cavity. In the vicinity of the ZGV frequency ($f = 1.91$ MHz), the wave numbers of $S_1$ and $S_{2b}$ almost coincide: $k_1 = 1.145 \text{ mm}^{-1}$ and $k_{2b} = 0.975 \text{ mm}^{-1}$ (see Supplementary Information S1). As a consequence, NR gives rise to a cylindrical wave front that converges back towards the initial source location [Fig. 3(c,f)], thus mimicking a PC experiment. Note that the NR wave field [Fig. 3(e)] is out-of-phase compared to the incident wave front [Fig. 3(a)], thereby confirming the theoretical prediction of a reflection coefficient, $\rho = -1$, between the $S_1$ and $S_{2b}$ modes at the ZGV-point⁴³. As for a TR experiment in a cavity, the size of the focal spot originating from the NR process is limited to a half-wavelength²⁰. This can be interpreted in terms of the diffraction limit and the loss of the evanescent components⁴⁶⁻⁴⁸.

The dynamics of the NR process is revealed by performing a discrete Fourier transform of the recorded signals over apodized time windows. The corresponding wave-fields are displayed in Fig. 4. The top and bottom panels show the temporal evolution of the $S_{2b}$ and $S_1$ modes, respectively, at the central frequency $f = 1.91$ MHz. In the time window [0–25] µs (short time), the $S_{2b}$-mode is largely predominant confirming its selective excitation by the transducer. The incident wave front then propagates across the cavity and finally vanishes when it reaches the cavity boundary. The residual $S_{2b}$ random wave field observed at long times of flight ([200–400] µs) can be explained by the spurious reflections undergone by the NR wave field when it converges back to the transducer. Figure 4(b) displays the time evolution of the $S_1$-mode. As soon as the incident $S_{2b}$-mode reaches the cavity boundary, it is converted into the NR $S_1$-mode. As in a PC or TR experiment, the NR wave field focuses back on the source location. A nice focus is obtained for instance over the time range [0–200] µs. Beyond that time, we observe a
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Negative reflection in the truncated billiard. (a) Incident wave field associated with the excited $S_{2b}$-mode at frequency $f = 1.91$ MHz. (b,c) Reflected wave field associated with the $S_1$-mode at frequency $f = 1.92$ MHz and $f = 1.91$ MHz, respectively. (d,e,f) Ray path trajectories involved in (a,b,c), respectively.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Time evolution of the negative reflection process taking place in the truncated billiard. The $S_{2b}$ and $S_1$ wave fields are displayed in panels (a and b), respectively, over the indicated time windows. Note the different color scales in (a and b).</figcaption></figure>

progressive degradation of the focal spot. Longer NR paths are actually more sensitive on the slight mismatch between the wave numbers $k_1$ and $k_{2b}$ in the vicinity of the ZGV-point (see Supplementary Information S1). Figures 3 and 4 illustrate the remarkable properties of the NR phenomenon in a complex medium. Whatever the chaotic feature of the environment, the incident wave field is re-created, but traveling backward, retracing its passage through the medium and focusing back to the source location. However, unlike TR or PC, a NR mirror is fully passive. No energy is injected into the system through a TR mirror or a nonlinear wave mixing process.

**Negative reflection by a random disordered slab.** Following the demonstration of the NR phenomenon in a cavity, we now address the more challenging case of a scattering layer made of randomly distributed holes. Such finite-size scatterers actually correspond to the small-scale limit of a curved plate edge. As NR occurs up to large angles of incidence$^{43}$, holes drilled in a plate constitute good candidates for finite-size NR mirrors.

The experiment is depicted in Fig. 2(b). The system under study is a plate of thickness $2h \sim 2$ mm for which the ZGV-point occurs at $f_{ZGV} = 1.41$ MHz and $k_{ZGV} = 0.89 \text{ mm}^{-1}$. Disorder is introduced by drilling circular holes with diameters ranging from 4 to 12.5 mm over a thickness of 20 mm from one of the free edge of the plate. The
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Negative reflection at a disordered plate edge. Normal displacement field filtered at 1.42 MHz (left) and 1.41 MHz close to the ZGV point (right). The Fourier analysis of the recorded signals performed on the indicated time windows yields the incident (a,b,d,e) and reflected (c,f) wave-fields. Time windows are different to account for the frequency-dependence of the group velocity for Lamb modes. Ray path trajectories are illustrated for incident (g) and negatively reflected waves with $k_1 \neq k_{2b}(h)$ and $k_1 = k_{2b}(i)$.</figcaption></figure>

holes diameters are comparable with the wavelength $\lambda_{ZGV} \sim 8$ mm. The ultrasound source is a piezoelectric strip of 25 mm-length, 5 mm-width and 1 mm-thickness glued on top of the plate [Fig. 2(b)]. This strip is located at 40 mm from the free edge and tilted 25° with respect to it. The source width ($\sim \lambda_{ZGV}/2$) allows to excite both the $S_1$ and $S_{2b}$ Lamb modes in the vicinity of the ZGV-point. The out-of-plane component of the local vibration of the plate is measured in the time domain during 1 ms over an area of dimensions $60 \times 120 \text{ mm}^2$ (see Methods). In order to isolate the incident and reflected wave fronts in the vicinity of the ZGV-point, a discrete Fourier transform of the recorded signals is performed over distinct time windows. Figure 5(a,b) show the incident wave front at frequency $f = 1.42$ MHz ($f \neq f_{ZGV}$), in the time windows $[0-50]$ µs and $[25-125]$ µs, respectively. Figure 5(d,e) displays the same wave-front but at frequency 1.41 MHz ($f \approx f_{ZGV}$), for the time windows $[0-200]$ µs and $[100-300]$ µs, respectively. The temporal windows considered at these two frequencies are different because the group velocity of Lamb modes strongly varies in the vicinity of the ZGV point. In both cases, the incident wave front is a wave packet of well-resolved momentum impinging the diffusive layer with an angle of incidence equal to 25° [Fig. 5(g)]. This wave packet is then scattered by the holes [Fig. 5(c,f)]. Figure 5(c) shows the reflected wave field in the time window $[100-200]$ µs and at frequency $f = 1.42$ MHz, that is well above the ZGV-point displayed in Fig. 2(c). The wave number mismatch between the $S_1$ and $S_{2b}$ modes results in a randomization of the reflected wave vector induced by the disordered layer [Fig. 5(h)] and a speckle-like feature of the reflected wave field in Fig. 5(c). On the contrary, in the vicinity of the ZGV frequency ($f = 1.41$ MHz), the incident and reflected wave vectors point exactly in opposite directions [Fig. 5(i)]. Figure 5(f) shows the corresponding NR wave field in the time window $[250-450]$ µs and at frequency $f = 1.41$ MHz. Despite disorder, NR gives rise to a well resolved wave packet with a momentum opposed to the incident one. The scattering layer not only retro-reflects the incident wave packet but also preserves its amplitude and phase. To some extent, the scattering layer is thus cloaked by the NR phenomenon. Whatever the position and concentration of scatterers, the NR wave-field always displays the same feature and corresponds to the PC replica of the incident wave-field. This experiment thus illustrates the robustness of the NR phenomenon.

Discussion

In this paper, the NR phenomenon is experimentally investigated in the time domain since it allows a clear discrimination between the incident and NR wave-fields. A first point we would like to discuss is the manifestation of the NR phenomenon in the frequency domain. NR can actually induce striking features in this regime. For sake of simplicity and generality, a scalar model is now considered to discuss about the potential impact of NR. Of course, it does not account for the full elastic problem but our aim is here to grasp the main features of NR for all kind of waves. Let us assume a monochromatic point source in a chaotic cavity [Fig. 2(a)]. The corresponding diverging incident wave-field is given by the free space Green's function $G_0(r|r_0)$. As shown in the Supplementary Information S2, multiple NR events on the cavity boundary have to be taken into account in the monochromatic regime. It yields the following expression for the Green's function $G(r|r_0)$ in the cavity,
$$G(\mathbf{r}|\mathbf{r}_0) = \frac{G_0(\mathbf{r}|\mathbf{r}_0)}{\underset{\text{incident wave-field}}{\overline{\phantom{G_0(\mathbf{r}|\mathbf{r}_0)}}}} + \frac{\frac{\rho}{1-\rho}[G_0(\mathbf{r}|\mathbf{r}_0) - G_0^*(\mathbf{r}|\mathbf{r}_0)]}{\underset{\text{NR wave-field}}{\overline{\phantom{G_0(\mathbf{r}|\mathbf{r}_0)}}}} \quad (1)$$

where the symbol * stands for phase conjugate. The NR wave-field is made of two contributions. The first component is a wave converging towards the initial source at $\mathbf{r}_0$, of opposite phase velocity compared to the incident wave-field. This wave is accounted for by $G_0(\mathbf{r}|\mathbf{r}_0)$ in the second term of Eq. (1). Owing to energy flux conservation$^{47}$, this wave does not stop when it reaches $\mathbf{r}_0$ and gives rise to a diverging wave accounted for by $G_0^*(\mathbf{r}|\mathbf{r}_0)$. As in a TR cavity experiment$^{46–48}$, the superimposition of the converging and diverging waves limits the focal spot size of the NR wave-field to one half wavelength [see Fig. 3(c)]. The pre-factor $\rho/(1-\rho)$ in Eq. (1) accounts for the multiple NR reflection events in the cavity [see Supplementary Information S2].

Depending on the sign of $\rho$, negative reflection can give rise to opposite interference phenomena. A reflection coefficient $\rho$ close to $-1$, as exhibited by the $S_1$ and $S_{2b}$ Lamb modes at a free edge$^{43}$, yields the following expression for $G(\mathbf{r}|\mathbf{r}_0)$:

$$\lim_{\rho \to -1} G(\mathbf{r}|\mathbf{r}_0) = \mathrm{Re}\{G_0(\mathbf{r}|\mathbf{r}_0)\}. \quad (2)$$

The destructive interference between the incident and NR waves gives rise to a stationary wave-field that coincides with the real part of the free-space Green's function. In 2D or 3D configurations, this implies a singularity of the wave-field at the source location and a super-focusing ability [see Supplementary Information S2].

On the contrary, a positive sign of the reflection coefficient yields a stationary wave-field proportional to the imaginary part of the free-space Green's function:

$$\lim_{\rho \to 1} G(\mathbf{r}|\mathbf{r}_0) = \frac{1}{\varepsilon} \mathrm{Im}\{G_0(\mathbf{r}|\mathbf{r}_0)\}, \quad (3)$$

with $\rho=1-\varepsilon$ and $\varepsilon \ll 1$. The constructive interference between the incident and NR waves gives rise to a resonant amplification of the wave-field accounted by the factor $\varepsilon^{-1}$ in the last equation. As in a TR experiment$^{47}$, the focal spot at the initial source location is diffraction-limited since $\mathrm{Im}\{G_0(\mathbf{r}|\mathbf{r}_0)\}$ exhibits a typical width of $\lambda/2$ [see Supplementary Information S2].

These fundamental results of Eqs (2) and (3) can also be extended to a source embedded into a random scattering medium at a depth much larger than the scattering mean free path. NR would induce a stationary wave-field and no energy could escape from the near-field of the source. The randomly disordered slab can therefore act as a band gap material. The observation of such super-focusing phenomena would require an experiment with an extremely fine frequency tuning$^{42}$. Moreover, it should also involve a non-invasive point-like source capable of a selective excitation of the forward or backward mode. Such experimental conditions are not met by the current experimental set up that operates in the time domain for a rigorous proof of the NR concept. An experimental demonstration of the NR super-focusing capabilities in the monochromatic regime thus constitutes a challenging perspective for this work.

From a more practical point-of-view, the experiment displayed in Fig. 5 illustrates how an ensemble of scatterers can become a NR mirror that retro-reflects the incident wave field towards the source itself. NR can thus be extremely useful for non-destructive testing applications. Lamb waves are indeed widely used in structural health monitoring for the detection of cracks in thin sheet materials and tubular products, such as aircraft shells or pipes for instance$^{49,50}$. Usually, an image of the structure is built by means of a pulse echo method with a single transducer. However, a large part of the scattered wave-field is not collected by this transducer because of diffraction. This effect can severely limit the penetration depth of ultrasound imaging techniques because the echoes of interest are drowned into a predominant structure-borne noise. On the contrary, NR can backscatter all the energy of the scattered wave towards the emitting transducer and the signal of interest can then emerge from noise. As a consequence, NR can be taken advantage of for overcoming the imaging depth limit of standard pulse echo methods in structural health monitoring.

In summary, we have demonstrated the NR phenomenon in the time domain with elastic guided waves in plates of arbitrary shape. Whatever the chaotic or random feature of the plate, NR at its boundaries or on the scatterers gives rise to a replica of the incident wave, following the same path but in a reverse way. NR displays focusing abilities similar to a PC mirror while being a fully passive process. On the one hand, NR of elastic waves can be extremely useful in structural health monitoring but may also find applications in the development of new acoustic devices including resonators, lenses, and filters. On the other hand, this phenomenon can also occur in zero-gap systems such as photonic crystals$^{34}$, chiral metamaterials$^{36}$ or graphene$^{32}$. NR can thus become a tool of choice for wave-front shaping and focusing in all fields of wave physics.

## Methods

**Experimental procedure.** The experiment consists in measuring the wave-field induced by a piezo-electric transducer or a piezo-electric strip glued on top of a plate [see Fig. 2(a) and (b)]. In the cavity, a 10 µs chirp signal spanning the frequency range 1.8–2.0 MHz is sent to the transducer, which generates a cylindrical incident wave front in the plate. In the disordered slab, a 15 µs chirp signal spanning the frequency range 1.3–1.5 MHz is sent to the piezoelectric strip, which generates a plane incident wave front. The out-of-plane component of the local vibration is detected on the other side of the plate by a photorefractive interferometer (Bossa Nova Tempo 1D). This probe is sensitive to the phase shift along the path of the optical probe beam. Signals detected by the optical probe are fed into a digital sampling oscilloscope and transferred to a computer. The corresponding wave field is measured over a time length of 1 ms and over a grid of points that map the plate surface with a pitch of 1 mm.
References

1. Fink, M. et al. Time-reversed acoustics. *Rep. Prog. Phys.* **63** (2000).

2. Lerosey, G. et al. Time reversal of electromagnetic waves. *Phys. Rev. Lett.* **92** (2004).

3. Bacot, V., Labousse, M., Eddi, A., Fink, M. & Fort, E. Time reversal and holography with spacetime transformations. *Nature Phys.* **12**, 972–977 (2016).

4. Hellwarth, R. W. Generation of time-reversed wave fronts by nonlinear refraction. *J. Opt. Soc. Am.* **67**, 1–3 (1977).

5. Yariv, A. Phase conjugation. *IEEE J. Quant. Elect.* **14**, 650–660 (1978).

6. Fisher, R. *Optical Phase Conjugation*. (Academic Press, Inc., London, 1983).

7. Mosk, A., Lagendijk, A., Lerosey, G. & Fink, M. Controlling waves in space and time for imaging and focusing in complex media. *Nat. Photon.* **6**, 283–292 (2012).

8. Pendry, J., Holden, A. J., Robbins, D. J. & Stewart, W. J. Magnetism from conductors and enhanced nonlinear phenomena. *IEEE Trans. Microwave Theory Tech.* **47**, 2075–2084 (1999).

9. Smith, D. R., Padilla, W. J., Vier, D. C., Nemat-Nasser, S. C. & Schultz, S. Composite medium with simultaneously negative permeability and permittivity. *Phys. Rev. Lett* **84**, 4184–4187 (2000).

10. Liu, Z. *et al.* Locally resonant sonic materials. *Science* **289**, 1734–1736 (2000).

11. Veselago, V. G. The electrodynamics of substances with simultaneously negative values of. *Sov. Phys. Usp.* **10**, 509 (1968).

12. Pendry, J. B. Negative refraction makes a perfect lens. *Phys. Rev. Lett.* **85**, 3966 (2000).

13. Shelby, R. A., Smith, D. R. & Schultz, S. Experimental verification of a negative index of refraction. *Science* **292**, 77–79 (2001).

14. Pendry, J. B., Schurig, D. & Smith, D. R. Controlling electromagnetic fields. *Science* **312**, 1780 (2006).

15. Leonhardt, U. Optical conformal mapping. *Science* **312**, 1777 (2006).

16. Yu, N. & Capasso, F. Flat optics with designer metasurfaces. *Nature Materials* **13**, 139–150 (2014).

17. Maslovski, S. & Tretyakov, S. Phase conjugation and perfect lensing. *J Appl. Phys.* **94**, 4241 (2003).

18. Pendry, J. Time reversal and negative refraction. *Science* **322**, 71–73 (2008).

19. Derode, A. & Fink, M. Robust acoustic time reversal with high-order multiple scattering. *Phys. Rev. Lett.* **75**, 4206–4210 (1995).

20. Draeger, C. & Fink, M. One-channel time reversal of elastic waves in a chaotic 2d-silicon cavity. *Phys. Rev. Lett.* **79**, 407–410 (1997).

21. Tolstoy, I. & Usdin, E. Wave propagation in elastic plates: low and high mode dispersion. *J. Acoust. Soc. Am.* **29**, 37–42 (1957).

22. Meitzler, A. H. Backward-wave transmission of stress pulses in elastic cylinders and plates. *J. Acoust. Soc. Am.* **38**, 835–842 (1965).

23. Mindlin, R. Waves and vibrations in isotropic elastic plates. *Structural Mechanics* **1960**, 199–232 (1960).

24. Negishi, K. Existence of negative group velocities in Lamb waves. *Jap. J. Appl. Phys.* **26**, 171 (1987).

25. Shuvalov, A. L. & Poncelet, O. On the backward lamb waves near thickness resonances in anisotropic plates. *Int. J. Solids Struct.* **45**, 3430–3448 (2008).

26. Laurent, J., Royer, D., Hussain, T., Ahmad, F. & Prada, C. Laser induced zero-group velocity resonances in transversely isotropic cylinder. *J. Acoust. Soc. Am.* **137**, 3325 (2015).

27. Cross, M. C. & Lifshitz, R. Elastic wave transmission at an abrupt junction in a thin plate with application to heat transport and vibrations in mesoscopic systems. *Phys. Rev. B* **64**, 085324 (2001).

28. Maznev, A. A. & Every, A. G. Existence of backward propagating acoustic waves in supported layers. *Wave Motion* **48**, 401–407 (2011).

29. Malischewsky, P. G., Forbriger, T. & Lomnitz, C. Unusual, equivocal rayleigh-dispersion curves for simple models taking into account the special propagation conditions in the valley of mexico city (cdmx) - preliminary results. *Geofisica Internacional* **56**, 7 (2017).

30. Royer, D. & Dieulesaint, E. *Elastic Waves in Solid I*. (Springer Verlag, Berlin, 2000).

31. Prada, C., Clorennec, D. & Royer, D. Local vibration of an elastic plate and zero-group velocity lamb modes. *J. Acoust. Soc. Am.* **124**, 203 (2008).

32. Cheianov, V. V., Falko, V. & Altshuler, B. L. The focusing of electron flow and a veselago lens in graphene p-n junctions. *Science* **315**, 1252–1255 (2007).

33. Joannopoulos, J. D., Johnson, S. G., Winn, J. N. & Meade, R. D. *Photonic Crystals - Molding the Flow of Light*. (Princeton University, Princeton, NJ, 2008).

34. Sivan, Y. & Pendry, J. B. Time reversal in dynamically tuned zero-gap periodic systems. *Phys. Rev. Lett.* **106**, 193902 (2011).

35. Pendry, J. A chiral route to negative refraction. *Science* **306**, 1353–1355 (2004).

36. Zhang, C. & Cui, T. J. Negative reflections of electromagnetic waves in a strong chiral medium. *Appl. Phys. Lett.* **91**, 194101 (2007).

37. Antoniades, M. A. & Eleftheriades, G. V. Compact linear lead/lag metamaterial phase shifters for broadband applications. *IEEE Antennas Wireless Propag. Lett.* **2**, 103–106 (2003).

38. Bramhavar, S. *et al*. Negative refraction and focusing of elastic Lamb waves at an interface. *Phys. Rev. B.* **83**, 014106 (2011).

39. Philippe, F. D., Murray, T. W. & Prada, C. Focusing on plates: Controlling guided waves using negative refraction. *Scient. Rep.* **5** (2015).

40. Legrand, F., Gérardin, B., Laurent, J., Prada, C. & Aubry, A. Negative refraction of Lamb modes: A theoretical study. *Phys. Rev. B* **98**, 214114 (2018).

41. Germano, M., Alippi, A., Bettucci, A. & Mancuso, G. Anomalous and negative reflection of Lamb waves in mode conversion. *Phys. Rev. B* **85**, 012102 (2012).

42. Veres, I. A., Grëunsteidl, C., Stobbe, D. M. & Murray, T. W. Broad-angle negative reflection and focusing of elastic waves from a plate edge. *Phys. Rev. B* **93**, 174304 (2016).

43. Gérardin, B., Laurent, J., Prada, C. & Aubry, A. Negative reflection of lamb waves at a free edge: Tunable focusing and mimicking phase conjugation. *J. Acoust. Soc. Am.* **140**, 591–600 (2016).

44. Prada, C., Balogun, O. & Murray, T. Laser-based ultrasonic generation and detection of zero-group velocity Lamb waves in thin plates. *Appl. Phys. Lett.* **87**, 194109 (2005).

45. Holland, S. D. & Chimenti, D. E. Air-coupled acoustic imaging with zero-group-velocity Lamb modes. *Appl. Phys. Lett.* **83**, 2704–2706 (2003).

46. Cassereau, D. & Fink, M. Time-reversal of ultrasonic fields. iii. theory of the closed time-reversal cavity. *IEEE Trans Ultrason Ferroelectr Freq Control.* **39**, 579–592 (1992).

47. de Rosny, J. & Fink, M. Overcoming the diffraction limit in wave physics using a time-reversal mirror and a novel acoustic sink. *Phys. Rev. Lett.* **89**, 124301 (2002).

48. Carminati, R., Pierrat, R., de Rosny, J. & Fink, M. Theory of the time reversal cavity for electromagnetic fields. *Opt. Lett.* **32**, 3107–3109 (2007).

49. Boller, C. Next generation structural health monitoring and its integration into aircraft design. *International Journal of Systems Science* **31**, 1333–1349 (2000).

50. Su, Z., Ye, L. & Lu, Y. Guided lamb waves for identification of damage in composite structures: A review. *J. Sound Vibr.* **295**, 753–780 (2006).
**Acknowledgements**

The authors are grateful for funding provided by LABEX WIFI (Laboratory of Excellence within the French Program Investments for the Future, ANR-10-LABX-24 and ANR-10-IDEX-0001-02 PSL*) and by the Agence Nationale de la Recherche (ANR-15-CE24-0014-01, Research Project COPPOLA). B.G. acknowledges financial support from the French "Direction Générale de l'Armement" (DGA).

**Author Contributions**

A.A. and C.P. suggested the experiments; B.G., J.L. and F.L. conducted the experiments; B.G. and A.A. analyzed the experimental results and all authors discussed the results; A.A. wrote the manuscript and all authors reviewed the manuscript.

**Additional Information**

**Supplementary information** accompanies this paper at https://doi.org/10.1038/s41598-019-38480-3.

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_8.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019