Dispersion free control of hydroelastic waves down to sub-wavelength scale

Dispersion free control of hydroelastic waves down to sub-wavelength scale

L. Domino Laboratoire de Physique et Mécanique des Milieux Hétérogènes (PMMH), UMR CNRS 7636 ; PSL - ESPCI, 10 rue Vauquelin, 75005 Paris, France; Sorbonne Université - UPMC, Univ. Paris 06; Sorbonne Paris Cité - UDD, Univ. Paris 07    M. Fermigier Laboratoire de Physique et Mécanique des Milieux Hétérogènes (PMMH), UMR CNRS 7636 ; PSL - ESPCI, 10 rue Vauquelin, 75005 Paris, France; Sorbonne Université - UPMC, Univ. Paris 06; Sorbonne Paris Cité - UDD, Univ. Paris 07    E. Fort Institut Langevin, ESPCI, CNRS, PSL Research University, 6 rue Jussieu, 75005, Paris, France    A. Eddi Laboratoire de Physique et Mécanique des Milieux Hétérogènes (PMMH), UMR CNRS 7636 ; PSL - ESPCI, 10 rue Vauquelin, 75005 Paris, France; Sorbonne Université - UPMC, Univ. Paris 06; Sorbonne Paris Cité - UDD, Univ. Paris 07

Hydroelastic surface waves propagate at the surface of water covered by a thin elastic sheet and can be directly measured with accurate space and time resolution. We present an experimental approach using hydroelastic waves that allows us to control waves down to the sub-wavelength scale. We tune the wave dispersion relation by varying locally the properties of the elastic cover and we introduce a local index contrast. This index contrast is independent of the frequency leading to a dispersion-free Snell-Descartes law for hydroelastic waves. We then show experimental evidence of broadband focusing, reflection and refraction of the waves. We also investigate the limits of diffraction through the example of a macroscopic analog to optical nanojets, revealing that any sub-wavelength configuration gives access to new features for surface waves.

Gravity-capillary waves have been extensively used as model waves to tackle the issue of wave control at macroscopic scale. Contrary to optics and acoustics, their temporal and spatial typical scales allow for direct and accurate observation of wave propagation inside the medium. The design and fabrication of media with given properties is generally obtained by tuning the local bathymetry, which modifies the wave phase velocity Lighthill (1978). Immersed structures have been used to obtain Anderson localization Belzons et al. (1988) or to create macroscopic metamaterials for wave focusing Hu and Chan (2005); Hu et al. (2013); Wang et al. (2014); Bobinski et al. (2015) and cloaking Farhat et al. (2008); Berraquero et al. (2013). Wave control is achieved for waves with wavelength larger or comparable to the liquid depth (shallow water approximation). In this regime, damping becomes a major issue in laboratory experiments as viscous friction at the bottom dissipates most of the mechanical wave energy. These two limitations narrow the effective bandwidth of the devices created with gravity-capillary waves.

Here, we propose a novel approach based on hydroelastic waves i.e. waves that propagate at the surface of water covered with an elastic sheet. These waves were initially introduced to describe motion in ice sheets located in the marginal ice zone Greenhill (1887); Davys et al. (1985); Schulkes et al. (1987); Squire et al. (1988) and later to study floating structures Watanabe et al. (2004), wakes in the lubrication approximation Arutkin et al. (2016) or cloaking Zareei and Alam (2016). In the limit of thin membranes their dispersion relation writes


where is the pulsation, is the wavenumber, m.s is the acceleration of gravity, is the mechanical tension in the elastic sheet, is the fluid density, is the flexural modulus of the elastic sheet and the fluid depth. The flexural modulus depends on the Young’s modulus of the material, its Poisson modulus and its thickness . Eq. 1 exhibits three distinct regimes depending on the material properties and the wave pulsation : gravity waves, tension waves and flexural waves. So far, very few experiments at the laboratory scale highlighted the flexural regime using whether thin elastic polymer sheets Deike et al. (2013); Montiel et al. (2013) or granular rafts Planchette et al. (2012).

Here, we propose to achieve spatial control of the propagation of hydroelastic waves by modifying the dispersion relation (eq. 1) through local variations of the sheet’s flexural modulus . We first describe our experimental set-up and verify quantitatively the prediction from eq. 1. We introduce a local index contrast using the local phase velocity. This index contrast is independent of the frequency, which allows us to define a dispersion-free Snell-Descartes law for hydroelastic waves. To show the versatility of the system we then implement more complex structures to focus wave energy and probe wave effects due to the finite size of the system.

I Experimental setup

Figure 1: (a) Sketch of the experimental set-up showing the container filled with water and covered with an elastic sheet, the wave generation device and the imaging system. (b) Measured dispersion relation for different film thicknesses m (squares), m (circles), m (diamonds) and m (triangles). For each thickness, the dashed line is the theoretical dispersion relations obtained using equation 1. The plain red line shows the theoretical gravity-capillary dispersion relation for water waves. Inset: Typical wave field measured for a point source vibrating at Hz. The wave travels from left to right.

We use a glass tank ( cm cm cm) filled to a depth cm of water. We cover its surface with a cm cm wide elastic sheet of thickness m made of an optically transparent silicone rubber sheet with Young’s modulus MPa, density kg/m and Poisson’s ratio . This elastic film floats freely at the surface of water so that the mechanical tension inside reduces to the water surface tension mN/m. The waves are generated with a vibration exciter powered with an amplifier controlled with a waveform generator. We work with frequencies ranging from to Hz and with amplitudes to ensure the waves are in the linear regime. In addition, to guarantee that we are in the thin membrane limit we check that , so that equation 1 is valid.

To analyze quantitatively the wave field we use the Free-Surface Synthetic Schlieren optical technique Moisy et al. (2009) based on the apparent displacement of a random dot pattern due to the local slope of the interface. The pattern is located underneath the tank and we observe it from the top using a CCD camera located at m from the fluid surface [Fig. 1(a)]. The sampling frequency of the camera is set to obtain stroboscopic images of the wave propagation with at least 12 images per period. The area we observe with our camera is about cm cm wide ( pixels). A Digital Image Correlation (DIC) algorithm (PIVlab Thielicke (2014); Thielicke and Stamhuis (2014)) is used to compute the displacement field between each recorded image and the reference image. After reconstruction, we obtain 2D elevation fields ( points) and we are able to measure amplitudes down to m. An example of the obtained circular height field generated by a point source is shown in [Fig. 1b] (inset).

Ii Dispersion relation

We first probe the validity of the theoretical dispersion relation predicted by eq. (1). We measure the wave field for circular waves at various frequencies in the range Hz obtained with a point source vibrating on elastic films with thicknesses m. For each field, we perform 2D spatial Fourier transforms to determine the wavenumber associated to each frequency. We plot in figure 1(b) the forcing frequency against the measured wavelength on a log-log scale. Our measurements show that the wavelength decreases with the forcing frequency , typically ranging from cm to cm. For low frequencies ( Hz) we observe a slope of revealing that the gravity term in eq. (1) is dominant. For larger frequencies, two regimes can be observed. For thin elastic sheets ( m) we observe a slope of which corresponds to the tension term . This experimental dispersion relation is in perfect agreement with theory. The transition between the gravity and the tension regime occurs for m. This part of the dispersion relation corresponds to standard water waves, confirming that tension in the film is solely due to the liquid surface tension. For thicker films the behavior is markedly different. The measured slope is showing that flexural term is leading. This is confirmed by the theoretical dispersion relation that is in excellent agreement with our experimental data. This flexural regime is reached when and . In the following we will only consider waves in the flexural regime, i.e. hydroelastic waves.

Iii Effective index and broadband refraction

Figure 2: (a) Wave field showing refraction at the interface between two media, with m (left) and m (right). The wave travels from the left to the right. The scale bar represents cm. (b) Sinus of the measured refracted angle as a function of the sinus of the incident angle for 3 different frequencies. The dashed line corresponds to Snell’s law prediction with . Inset: Schematic drawing of the experiment, showing the wave vectors and the angles of the incident wave and of the transmitted wave . (c) Wave field showing the total reflection at the interface between the two media previously described. The wave travels from the left to the right. The scale bar represents cm. (d) Blue circles: normalized profile of the intensity of the evanescent wave taken along the dashed line shown in (c). Red line: exponential fit.

In this regime the dispersion relation (eq. 1) in the deep water approximation can be simplified as


The phase velocity then only depends on the film properties and the wavenumber . From this phase velocity we can define a relative effective refractive index . is spatially tunable by varying locally the value of the flexural coefficient . This can be achieved by changing locally the Young’s modulus or the film thickness . For two domains covered with elastic films with different thicknesses and and same Young’s modulus , the ratio of their refractive indices and writes


Thicker regions (resp. thinner) thus correspond to smaller (resp. higher) refractive indices. In the hydroelastic regime the index ratio is given by the local values of the film thickness and does not depend on nor .

We perform experiments to test eq. 3 through Snell-Descartes law using the refraction of a plane wave at an interface between two media. The interface is obtained using two thicknesses (m and m, respectively) [Fig. 2(a)], and the frequency ranges from Hz. The incident (resp. transmitted) waves have a wave vector (resp. ) that forms an angle (resp. ) with the interface normal [Fig. 2b]. We measure these angles by means of spatial Fourier transforms for varying incidence angle, ranging from to . We plot in figure 2(b) against for three different frequencies. The result is linear which means that hydroelastic waves obey the Snell-Descartes law of refraction: . This result based on translation invariance holds whatever the frequency. The expected slope given by the refractive index ratio is in excellent agreement with our experimental data.

We also study the situation where the angle of incidence is larger than the critical angle, here . Such a wave field is presented in figure 2(c) where . As expected, the wave undergoes a Total Internal Reflection (TIR) with the presence of an evanescent wave in medium 2. No energy is transmitted through the interface but this evanescent wave can be observed: the height of the wave decreases rapidly away from the interface. We can directly measure the amplitude of these waves on our fields, and we represent the profile obtained in the inset of figure 2(d). The amplitude of these waves decreases exponentially, with a typical penetration length . This penetration length is in good agreement with the expected value of obtained with .

Iv Lenses, focalisation

Figure 3: (a) Wave field for a lens with a radius of curvature of cm. Circular waves traveling at Hz from the left to the right. The plain red line denotes the shape of the lens, the dashed line denotes its focal plane. (b) Measurement of as a function of for various lenses and 5 different frequencies: Hz (squares), Hz (diamonds), Hz (up triangles), Hz (down triangles) and Hz (circles). The error on the measure is lower than and is much smaller than the marker size. Plain line denotes the theoretical focal length for a thin lens (see text). (c) Profile of the intensity field along the focal plane for a lens with radius cm and for Hz.

Fine sub-wavelength wave control can be achieved easily by designing engineered shapes to focus and guide waves. Here, as an example, we design 2D lenses by cutting out symmetric circular arcs in the silicon polymer. The obtained shapes [Figure 3] are defined by their radius of curvature . We then deposit them on the first membrane to locally increase the thickness. Note that this shape should create a convergent lens as thicker regions have a lower refractive index . We excite the system using a point source located on the left of the lens. Fig. 3(a) presents a typical wave field for a lens with cm excited with Hz and shows a focal spot on the right side of the lens (see also supplementary videos: Movie1.mov and Movie2.mov). Using the location of this focal spot and that of the source we define the focal length as , where (resp. ) is the distance between the lens and the source (resp. the image). [Fig. 3(b)] shows the measured as a function of for different frequencies ranging from Hz to Hz. We observe that increases linearly with while being independent of as the refractive index ratio only depends on the film thickness ratio in the flexural regime. Here the typical size of the lens compares with and the paraxial approximation is clearly not satisfied. This makes ray optics a poor candidate to model our results. However, its prediction agrees surprisingly well with our experimental data [Fig. 3(b)].

We now characterize the profile of the wave field at the focus. Fig. 3 (c) presents the lateral normalized intensity profile of the focal spot obtained for a lens with a curvature radius of cm at Hz. The profile exhibits a central peak with a Full Width at Half Maximum (FWHM) of , and 2 side lobes. The presence of these secondary peaks is the signature of diffraction: as both the typical width of the focal spot and the typical size of the lens compare with , wave propagation should be described at the wavelength scale.

V Sub-wavelength focusing

Figure 4: (a) Intensity profile for a plane wave propagating at Hz through a thin disk of diameter , here denoted with the dark dashed line. The film is m thick inside the circle, and m thick elsewhere. We take an intensity profile along the red vertical dashed line. (b) Intensity profile taken along the red dashed line plotted in (a).

To further show the versatility of hydroelastic wave to control waves at the subwavelength scale, we use a geometry known in optics as “nanojet”. They were first introduced in optics using small cylindrical (or spherical) structures () with a strong index contrast Chen et al. (2004); Ferrand et al. (2008); Geints et al. (2010). These structures are traditionally manufactured or simulated using a dielectric sphere that has a higher refractive index than the surrounding medium, like glass, water or latex. The focal spot is created by the combination of evanescent and propagating waves in the shadow-side of the sphere.

We transpose this object to 2D hydroelastic waves by creating a thinner circular area with diameter in the elastic sheet [Fig. 4(a)]. As for classical microspheres in optics, the refractive index in this region is higher than in the outside medium. [Fig. 4(a)] shows the measured intensity field inside and outside the disk for an incident wave with Hz and mm (see also supplementary videos Movie1.mov and Movie2.mov). The wave propagates from the left to the right, and the thinner region is denoted with the dark dashed circle. We observe that the circular patch distorts the incoming plane wave and that a strong focal spot emerges out of the circular area. We show in Fig. 4(b) the lateral intensity profile of this focal spot that presents a very narrow peak as well as small side lobes on both sides. The FWHM of this focal spot is , which is smaller than the classical diffraction limit at , confirming that sub-wavelength focusing can be achieved using this simple design at the wavelength scale.

Vi Conclusion

We have achieved control of hydroelastic waves propagation in a model experiment. We first confirmed that the waves can be accurately described by the dispersion relation (Eq. 1) at the laboratory scale, and that the elastic sheet’s properties have a crucial incidence on the wave propagation in the flexural waves regime. Indeed, the flexion modulus can be tuned spatially by locally modifying the elastic film’s thickness. We build a local index contrast that only depends on the film’s flexion modulus and is therefore dispersion free despite the dispersive nature of the waves. Using this we first show that Fermat principle applies for hydroelastic waves, as refraction of an incident plane wave on a straigth interface obeys Snell-Descartes law independently of the incident wave frequency. With this feature we implement lenses with tunable focalisation properties. Nevertheless, the index variations occur on a typical scale that compares with the hydroelastic wavelength leading to subtle wave effects. This is particularly revealed by the construction of a macroscopic equivalent of nanojets. These simple circular structures allow to overcome the diffraction limit leading to a focal spot as small as .

We believe our macroscopic experiment can be used as a model experiment to study the physics of waves with new features so far unachieved. Direct observation of the waves combined with the ability to tune the medium’s properties down to sub-wavelength scale opens promising perspectives to probe wave propagation in structured Krauss et al. (1996); Vasseur et al. (2001); Pendry et al. (2006); Brûlé et al. (2014) or random media Duvall Jr et al. (1993); Vellekoop et al. (2010); Mosk et al. (2012). In particular this system allows for precise spatio-temporal control of wave sources as well as a precise monitoring of dynamical effects. Inspired by recent work on time-reversal of gravito-capillary waves Bacot et al. (2016), we now aim to implement macroscopic dynamical spatial structures Chumak et al. (2010); Sivan and Pendry (2011); Yuan et al. (2016).


  • Lighthill (1978) J. Lighthill, Cambridge, Cambrige  (1978).
  • Belzons et al. (1988) M. Belzons, E. Guazzelli,  and O. Parodi, Journal of Fluid Mechanics 186, 539 (1988).
  • Hu and Chan (2005) X. Hu and C. Chan, Physical review letters 95, 154501 (2005).
  • Hu et al. (2013) X. Hu, J. Yang, J. Zi, C. T. Chan,  and K.-M. Ho, Scientific reports 3, 1916 (2013).
  • Wang et al. (2014) Z. Wang, P. Zhang, X. Nie,  and Y. Zhang, EPL (Europhysics Letters) 108, 24003 (2014).
  • Bobinski et al. (2015) T. Bobinski, A. Eddi, P. Petitjeans, A. Maurel,  and V. Pagneux, Applied Physics Letters 107, 014101 (2015).
  • Farhat et al. (2008) M. Farhat, S. Enoch, S. Guenneau,  and A. Movchan, Physical review letters 101, 134501 (2008).
  • Berraquero et al. (2013) C. Berraquero, A. Maurel, P. Petitjeans,  and V. Pagneux, Physical review E 88, 051002 (2013).
  • Greenhill (1887) A. Greenhill, American Journal of Mathematics , 97 (1887).
  • Davys et al. (1985) J. Davys, R. Hosking,  and A. Sneyd, Journal of Fluid Mechanics 158, 269 (1985).
  • Schulkes et al. (1987) R. M. S. M. Schulkes, R. J. Hosking,  and A. D. Sneyd, Journal of Fluid Mechanics 180, 297 (1987).
  • Squire et al. (1988) V. A. Squire, W. H. Robinson, P. J. Langhorne,  and T. G. Haskell, Nature 333, 159 (1988).
  • Watanabe et al. (2004) E. Watanabe, T. Utsunomiya,  and C. Wang, Engineering structures 26, 245 (2004).
  • Arutkin et al. (2016) M. Arutkin, R. Ledesma-Alonso, T. Salez,  and É. Raphaël, arXiv preprint arXiv:1607.07816  (2016).
  • Zareei and Alam (2016) A. Zareei and R. Alam, in Proc. 31st Int. Workshop on Water Waves and Floating Bodies, Michigan, USA (2016) pp. 197–200.
  • Deike et al. (2013) L. Deike, J.-C. Bacri,  and E. Falcon, Journal of Fluid Mechanics 733, 394 (2013).
  • Montiel et al. (2013) F. Montiel, F. Bonnefoy, P. Ferrant, L. G. Bennetts, V. A. Squire,  and P. Marsault, Journal of Fluid Mechanics 723, 604 (2013).
  • Planchette et al. (2012) C. Planchette, E. Lorenceau,  and A.-L. Biance, Soft Matter 8, 2444 (2012).
  • Moisy et al. (2009) F. Moisy, M. Rabaud,  and K. Salsac, Experiments in Fluids 46, 1021 (2009).
  • Thielicke (2014) W. Thielicke, The flapping flight of birds: Analysis and application ([S.n.], 2014).
  • Thielicke and Stamhuis (2014) W. Thielicke and E. Stamhuis, Journal of Open Research Software 2 (2014).
  • Chen et al. (2004) Z. Chen, A. Taflove,  and V. Backman, Opt. Express 12, 1214 (2004).
  • Ferrand et al. (2008) P. Ferrand, J. Wenger, A. Devilez, M. Pianta, B. Stout, N. Bonod, E. Popov,  and H. Rigneault, Opt. Express 16, 6930 (2008).
  • Geints et al. (2010) Y. Geints, E. Panina,  and A. Zemlyanov, Optics Communications 283, 4775 (2010).
  • Krauss et al. (1996) T. F. Krauss, M. Richard,  and S. Brand, Nature 383, 699 (1996).
  • Vasseur et al. (2001) J. Vasseur, P. A. Deymier, B. Chenni, B. Djafari-Rouhani, L. Dobrzynski,  and D. Prevost, Physical Review Letters 86, 3012 (2001).
  • Pendry et al. (2006) J. B. Pendry, D. Schurig,  and D. R. Smith, science 312, 1780 (2006).
  • Brûlé et al. (2014) S. Brûlé, E. H. Javelaud, S. Enoch,  and S. Guenneau, Phys. Rev. Lett. 112, 133901 (2014).
  • Duvall Jr et al. (1993) T. Duvall Jr, S. Jefferies, J. Harvey,  and M. Pomerantz, Nature 362, 430 (1993).
  • Vellekoop et al. (2010) I. Vellekoop, A. Lagendijk,  and A. Mosk, Nature photonics 4, 320 (2010).
  • Mosk et al. (2012) A. P. Mosk, A. Lagendijk, G. Lerosey,  and M. Fink, Nature photonics 6, 283 (2012).
  • Bacot et al. (2016) V. Bacot, M. Labousse, A. Eddi, M. Fink,  and E. Fort, Nat Phys 12, 972 (2016).
  • Chumak et al. (2010) A. V. Chumak, V. S. Tiberkevich, A. D. Karenowska, A. A. Serga, J. F. Gregg, A. N. Slavin,  and B. Hillebrands, Nature Communications 1, 141 EP (2010).
  • Sivan and Pendry (2011) Y. Sivan and J. B. Pendry, Opt. Express 19, 14502 (2011).
  • Yuan et al. (2016) L. Yuan, M. Xiao,  and S. Fan, Phys. Rev. B 94, 140303 (2016).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description