Computing leaky Lamb waves for waveguides between elastic half-spaces using spectral collocation.

In non-destructive evaluation guided wave inspections, the elastic structure to be inspected is often embedded within other elastic media and the ensuing leaky waves are complex and non-trivial to compute; we consider the canonical example of an elastic waveguide surrounded by other elastic materials that demonstrates the fundamental issues with calculating the leaky waves in such systems. Due to the complex wavenumber solutions required to represent them, leaky waves pose significant challenges to existing numerical methods, with methods that spatially discretise the field to retrieve them suffering from the exponential growth of their amplitude far into the surrounding media. We present a spectral collocation method yielding an accurate and efficient identification of these modes, leaking into elastic half-spaces. We discretise the elastic domains and, depending on the exterior bulk wavespeeds, select appropriate mappings of the discretised domain to complex paths, in which the numerical solution decays and the physics of the problem are preserved. By iterating through all possible radiation cases, the full set of dispersion and attenuation curves are successfully retrieved and validated, where possible, against the commercially available software disperse. As an independent validation, dispersion curves are obtained from finite element simulations of time-dependent waves using Fourier analysis.

Reinforcement learning optimisation for graded metamaterial design using a physical-based constraint on the state representation and action space.

The energy harvesting capability of a graded metamaterial is maximised via reinforcement learning (RL) under realistic excitations at the microscale. The metamaterial consists of a waveguide with a set of beam-like resonators of variable length, with piezoelectric patches, attached to it. The piezo-mechanical system is modelled through equivalent lumped parameters determined via a general impedance analysis. Realistic conditions are mimicked by considering either magnetic loading or random excitations, the latter scenario requiring the enhancement of the harvesting capability for a class of forcing terms with similar but different frequency content. The RL-based optimisation is empowered by using the physical understanding of wave propagation in a such local resonance system to constrain the state representation and the action space. The procedure outcomes are compared against grading rules optimised through genetic algorithms. While genetic algorithms are more effective in the deterministic setting featuring the application of magnetic loading, the proposed RL-based proves superior in the inherently stochastic setting of the random excitation scenario.

Models for resonant acoustic metasurfaces with application to moth wing ultrasound absorption.

Taking as bioinspiration the remarkable acoustic absorption properties of moth wings, we develop a simple analytical model that describes the interaction between acoustic pressure fields, and thin elastic plates incorporating resonant sub-structures. The moth wing is an exemplar of a natural acoustic metamaterial; the wings are deeply subwavelength in thickness at the frequencies of interest, the absorption is broadband and the tiny scales resonate on the moth wing acting in concert. The simplified model incorporates only the essential physics and the scales are idealized to flat rigid rectangular plates coupled via a spring to an elastic plate that forms the wing; all the components are deep-subwavelength at desired frequencies. Based on Fourier analysis, complemented by phenomenological modelling, our theory shows excellent agreement with simulation mimicking the moth-wing structure. Moth wings operate as broadband sound absorbers employing a range of scale sizes. We demonstrate that a random distribution of scale sizes generates a broadband absorption spectrum. To further illustrate the potential of the model, we design a deeply sub-wavelength acoustic counterpart of electromagnetically induced reflectance. This article is part of the theme issue 'Wave generation and transmission in multi-scale complex media and structured metamaterials (part 2)'.

Leaky wave characterisation using spectral methods.

Leaky waves are an important class of waves, particularly for guiding waves along structures embedded within another medium; a mismatch in wavespeeds often leads to leakage of energy from the waveguide, or interface, into the medium, which consequently attenuates the guided wave. The accurate and efficient identification of theoretical solutions for leaky waves is a key requirement for the choices of modes and frequencies required for non-destructive evaluation inspection techniques. We choose a typical situation to study: an elastic waveguide with a fluid on either side. Historically, leaky waves are identified via root-finding methods that have issues with conditioning, or numerical methods that struggle with the exponential growth of solutions at infinity. By building upon a spectral collocation method, we show how it can be adjusted to find exponentially growing solutions, i.e., leaky waves, leading to an accurate, fast, and efficient identification of their dispersion properties. The key concept required is a mapping, in the fluid region, that allows for exponential growth of the physical solution at infinity, whilst the mapped numerical setting decays. We illustrate this by studying leaky Lamb waves in an elastic waveguide immersed between two different fluids and verify this using the commercially available software disperse.

Drying-induced stresses in poroelastic drops on rigid substrates.

We develop a theory for drying-induced stresses in sessile, poroelastic drops undergoing evaporation on rigid surfaces. Using a lubrication-like approximation, the governing equations of three-dimensional nonlinear poroelasticity are reduced to a single thin-film equation for the drop thickness. We find that thin drops experience compressive elastic stresses but the total in-plane stresses are tensile. The mechanical response of the drop is dictated by the initial profile of the solid skeleton, which controls the in-plane deformation, the dominant components of elastic stress, and sets a limit on the depth of delamination that can potentially occur. Our theory suggests that the alignment of desiccation fractures in colloidal drops is selected by the shape of the drop at the point of gelation. We propose that the emergence of three distinct fracture patterns in dried blood drops is a consequence of a nonmonotonic drop profile at gelation. We also show that depletion fronts, which separate wet and dry solid, can invade the drop from the contact line and localize the generation of mechanical stress during drying. Finally, the finite element method is used to explore the stress profiles in drops with large contact angles.

Octet lattice-based plate for elastic wave control.

Motivated by the importance of lattice structures in multiple fields, we numerically investigate the propagation of flexural waves in a thin reticulated plate augmented with two classes of metastructures for wave mitigation and guiding, namely metabarriers and metalenses. The cellular architecture of this plate invokes the well-known octet topology, while the metadevices rely on novel customized octets either comprising spherical masses added to the midpoint of their struts or variable node thickness. We numerically determine the dispersion curves of a doubly-periodic array of octets, which produce a broad bandgap whose underlying physics is elucidated and leveraged as a design paradigm, allowing the construction of a metabarrier effective for inhibiting the transmission of waves. More sophisticated effects emerge upon parametric analyses of the added masses and node thickness, leading to graded designs that spatially filter waves through an enlarged bandgap via rainbow trapping. Additionally, Luneburg and Maxwell metalenses are realized using the spatial modulation of the tuning parameters and numerically tested. Wavefronts impinging on these structures are progressively curved within the inhomogeneous media and steered toward a focal point. Our results yield new perspectives for the use of octet-like lattices, paving the way for promising applications in vibration isolation and energy focusing.

Tailored elastic surface to body wave Umklapp conversion.

Elastic waves guided along surfaces dominate applications in geophysics, ultrasonic inspection, mechanical vibration, and surface acoustic wave devices; precise manipulation of surface Rayleigh waves and their coupling with polarised body waves presents a challenge that offers to unlock the flexibility in wave transport required for efficient energy harvesting and vibration mitigation devices. We design elastic metasurfaces, consisting of a graded array of rod resonators attached to an elastic substrate that, together with critical insight from Umklapp scattering in phonon-electron systems, allow us to leverage the transfer of crystal momentum; we mode-convert Rayleigh surface waves into bulk waves that form tunable beams. Experiments, theory and simulation verify that these tailored Umklapp mechanisms play a key role in coupling surface Rayleigh waves to reversed bulk shear and compressional waves independently, thereby creating passive self-phased arrays allowing for tunable redirection and wave focusing within the bulk medium.

Fresnel drag in space-time-modulated metamaterials.

A moving medium drags light along with it as measured by Fizeau and explained by Einstein's theory of special relativity. Here we show that the same effect can be obtained in a situation where there is no physical motion of the medium. Modulations of both the permittivity and permeability, phased in space and time in the form of traveling waves, are the basis of our model. Space-time metamaterials are represented by effective bianisotropic parameters, which can in turn be mapped to a moving homogeneous medium. Hence these metamaterials mimic a relativistic effect without the need for any actual material motion. We discuss how both the permittivity and permeability need to be modulated to achieve these effects, and we present an equivalent transmission line model.

Extraordinarily transparent compact metallic metamaterials.

The design of achromatic optical components requires materials with high transparency and low dispersion. We show that although metals are highly opaque, densely packed arrays of metallic nanoparticles can be more transparent to infrared radiation than dielectrics such as germanium, even when the arrays are over 75% metal by volume. Such arrays form effective dielectrics that are virtually dispersion-free over ultra-broadband ranges of wavelengths from microns up to millimeters or more. Furthermore, the local refractive indices may be tuned by altering the size, shape, and spacing of the nanoparticles, allowing the design of gradient-index lenses that guide and focus light on the microscale. The electric field is also strongly concentrated in the gaps between the metallic nanoparticles, and the simultaneous focusing and squeezing of the electric field produces strong 'doubly-enhanced' hotspots which could boost measurements made using infrared spectroscopy and other non-linear processes over a broad range of frequencies.

Moving Contact Lines: Linking Molecular Dynamics and Continuum-Scale Modeling.

Despite decades of research, the modeling of moving contact lines has remained a formidable challenge in fluid dynamics whose resolution will impact numerous industrial, biological, and daily life applications. On the one hand, molecular dynamics (MD) simulation has the ability to provide unique insight into the microscopic details that determine the dynamic behavior of the contact line, which is not possible with either continuum-scale simulations or experiments. On the other hand, continuum-based models provide a link to the macroscopic description of the system. In this Feature Article, we explore the complex range of physical factors, including the presence of surfactants, which governs the contact line motion through MD simulations. We also discuss links between continuum- and molecular-scale modeling and highlight the opportunities for future developments in this area.

Acoustic Far-Field Hypersonic Surface Wave Detection with Single Plasmonic Nanoantennas.

The optical properties of small metallic particles allow us to bridge the gap between the myriad of subdiffraction local phenomena and macroscopic optical elements. The optomechanical coupling between mechanical vibrations of Au nanoparticles and their optical response due to collective electronic oscillations leads to the emission and the detection of surface acoustic waves (SAWs) by single metallic nanoantennas. We take two Au nanoparticles, one acting as a source and the other as a receptor of SAWs and, even though these antennas are separated by distances orders of magnitude larger than the characteristic subnanometric displacements of vibrations, we probe the frequency content, wave speed, and amplitude decay of SAWs originating from the damping of coherent mechanical modes of the source. Two-color pump-probe experiments and numerical methods reveal the characteristic Rayleigh wave behavior of emitted SAWs, and show that the SAW-induced optical modulation of the receptor antenna allows us to accurately probe the frequency of the source, even when the eigenmodes of source and receptor are detuned.

Simulation of immiscible liquid-liquid flows in complex microchannel geometries using a front-tracking scheme.

The three-dimensional two-phase flow dynamics inside a microfluidic device of complex geometry is simulated using a parallel, hybrid front-tracking/level-set solver. The numerical framework employed circumvents numerous meshing issues normally associated with constructing complex geometries within typical computational fluid dynamics packages. The device considered in the present work is constructed via a module that defines solid objects by means of a static distance function. The construction combines primitive objects, such as a cylinder, a plane, and a torus, for instance, using simple geometrical operations. The numerical solutions predicted encompass dripping and jetting, and transitions in flow patterns are observed featuring the formation of drops, 'pancakes', plugs, and jets, over a wide range of flow rate ratios. We demonstrate the fact that vortex formation accompanies the development of certain flow patterns, and elucidate its role in their underlying mechanisms. Experimental visualisation with a high-speed imaging are also carried out. The numerical predictions are in excellent agreement with the experimental data.

Enhanced sensing and conversion of ultrasonic Rayleigh waves by elastic metasurfaces.

Recent years have heralded the introduction of metasurfaces that advantageously combine the vision of sub-wavelength wave manipulation, with the design, fabrication and size advantages associated with surface excitation. An important topic within metasurfaces is the tailored rainbow trapping and selective spatial frequency separation of electromagnetic and acoustic waves using graded metasurfaces. This frequency dependent trapping and spatial frequency segregation has implications for energy concentrators and associated energy harvesting, sensing and wave filtering techniques. Different demonstrations of acoustic and electromagnetic rainbow devices have been performed, however not for deep elastic substrates that support both shear and compressional waves, together with surface Rayleigh waves; these allow not only for Rayleigh wave rainbow effects to exist but also for mode conversion from surface into shear waves. Here we demonstrate experimentally not only elastic Rayleigh wave rainbow trapping, by taking advantage of a stop-band for surface waves, but also selective mode conversion of surface Rayleigh waves to shear waves. These experiments performed at ultrasonic frequencies, in the range of 400-600 kHz, are complemented by time domain numerical simulations. The metasurfaces we design are not limited to guided ultrasonic waves and are a general phenomenon in elastic waves that can be translated across scales.

Physical insights into the blood-brain barrier translocation mechanisms.

The number of individuals suffering from diseases of the central nervous system (CNS) is growing with an aging population. While candidate drugs for many of these diseases are available, most of these pharmaceutical agents cannot reach the brain rendering most of the drug therapies that target the CNS inefficient. The reason is the blood-brain barrier (BBB), a complex and dynamic interface that controls the influx and efflux of substances through a number of different translocation mechanisms. Here, we present these mechanisms providing, also, the necessary background related to the morphology and various characteristics of the BBB. Moreover, we discuss various numerical and simulation approaches used to study the BBB, and possible future directions based on multi-scale methods. We anticipate that this review will motivate multi-disciplinary research on the BBB aiming at the design of effective drug therapies.

High-frequency homogenization for travelling waves in periodic media.

We consider high-frequency homogenization in periodic media for travelling waves of several different equations: the wave equation for scalar-valued waves such as acoustics; the wave equation for vector-valued waves such as electromagnetism and elasticity; and a system that encompasses the Schrödinger equation. This homogenization applies when the wavelength is of the order of the size of the medium periodicity cell. The travelling wave is assumed to be the sum of two waves: a modulated Bloch carrier wave having crystal wavevector [Formula: see text] and frequency ω1 plus a modulated Bloch carrier wave having crystal wavevector [Formula: see text] and frequency ω2. We derive effective equations for the modulating functions, and then prove that there is no coupling in the effective equations between the two different waves both in the scalar and the system cases. To be precise, we prove that there is no coupling unless ω1=ω2 and [Formula: see text] where Λ=(λ1λ2…λ d ) is the periodicity cell of the medium and for any two vectors [Formula: see text] the product a⊙b is defined to be the vector (a1b1,a2b2,…,adbd ). This last condition forces the carrier waves to be equivalent Bloch waves meaning that the coupling constants in the system of effective equations vanish. We use two-scale analysis and some new weak-convergence type lemmas. The analysis is not at the same level of rigour as that of Allaire and co-workers who use two-scale convergence theory to treat the problem, but has the advantage of simplicity which will allow it to be easily extended to the case where there is degeneracy of the Bloch eigenvalue.

Surface plasmon resonances of arbitrarily shaped nanometallic structures in the small-screening-length limit.

According to the hydrodynamic Drude model, surface plasmon resonances of metallic nanostructures blueshift owing to the non-local response of the metal's electron gas. The screening length characterizing the non-local effect is often small relative to the overall dimensions of the metallic structure, which enables us to derive a coarse-grained non-local description using matched asymptotic expansions; a perturbation theory for the blueshifts of arbitrary-shaped nanometallic structures is then developed. The effect of non-locality is not always a perturbation and we present a detailed analysis of the 'bonding' modes of a dimer of nearly touching nanowires where the leading-order eigenfrequencies and eigenmode distributions are shown to be a renormalization of those predicted assuming a local metal permittivity.

A seismic metamaterial: The resonant metawedge.

Critical concepts from three different fields, elasticity, plasmonics and metamaterials, are brought together to design a metasurface at the geophysical scale, the resonant metawedge, to control seismic Rayleigh waves. Made of spatially graded vertical subwavelength resonators on an elastic substrate, the metawedge can either mode convert incident surface Rayleigh waves into bulk elastic shear waves or reflect the Rayleigh waves creating a "seismic rainbow" effect analogous to the optical rainbow for electromagnetic metasurfaces. Time-domain spectral element simulations demonstrate the broadband efficacy of the metawedge in mode conversion while an analytical model is developed to accurately describe and predict the seismic rainbow effect; allowing the metawedge to be designed without the need for extensive parametric studies and simulations. The efficiency of the resonant metawedge shows that large-scale mechanical metamaterials are feasible, will have application, and that the time is ripe for considering many optical devices in the seismic and geophysical context.

Transformation seismology: composite soil lenses for steering surface elastic Rayleigh waves.

Metamaterials are artificially structured media that exibit properties beyond those usually encountered in nature. Typically they are developed for electromagnetic waves at millimetric down to nanometric scales, or for acoustics, at centimeter scales. By applying ideas from transformation optics we can steer Rayleigh-surface waves that are solutions of the vector Navier equations of elastodynamics. As a paradigm of the conformal geophysics that we are creating, we design a square arrangement of Luneburg lenses to reroute Rayleigh waves around a building with the dual aim of protection and minimizing the effect on the wavefront (cloaking). To show that this is practically realisable we deliberately choose to use material parameters readily available and this metalens consists of a composite soil structured with buried pillars made of softer material. The regular lattice of inclusions is homogenized to give an effective material with a radially varying velocity profile and hence varying the refractive index of the lens. We develop the theory and then use full 3D numerical simulations to conclusively demonstrate, at frequencies of seismological relevance 3-10 Hz, and for low-speed sedimentary soil (vs: 300-500 m/s), that the vibration of a structure is reduced by up to 6 dB at its resonance frequency.

Forests as a natural seismic metamaterial: Rayleigh wave bandgaps induced by local resonances.

We explore the thesis that resonances in trees result in forests acting as locally resonant metamaterials for Rayleigh surface waves in the geophysics context. A geophysical experiment demonstrates that a Rayleigh wave, propagating in soft sedimentary soil at frequencies lower than 150 Hz, experiences strong attenuation, when interacting with a forest, over two separate large frequency bands. This experiment is interpreted using finite element simulations that demonstrate the observed attenuation is due to bandgaps when the trees are arranged at the sub-wavelength scale with respect to the incident Rayleigh wave. The repetitive bandgaps are generated by the coupling of the successive longitudinal resonances of trees with the vertical component of the Rayleigh wave. For wavelengths down to 5 meters, the resulting bandgaps are remarkably large and strongly attenuating when the acoustic impedance of the trees matches the impedance of the soil. Since longitudinal resonances of a vertical resonator are inversely proportional to its length, a man-made engineered array of resonators that attenuates Rayleigh waves at frequency ≤10 Hz could be designed starting from vertical pillars coupled to the ground with longitudinal resonance ≤10 Hz.