Richard V. Craster

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 150u2009Hz, 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 ≤10u2009Hz could be designed starting from vertical pillars coupled to the ground with longitudinal resonance ≤10u2009Hz.

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.

Seismic metasurfaces: Sub-wavelength resonators and Rayleigh wave interaction

We consider the canonical problem of an array of rods, which act as resonators, placed on an elastic substrate; the substrate being either a thin elastic plate or an elastic half-space. In both cases the flexural plate, or Rayleigh surface, waves in the substrate interact with the resonators to create interesting effects such as effective band-gaps for surface waves or filters that transform surface waves into bulk waves; these effects have parallels in the field of optics where such sub-wavelength resonators create metamaterials in the bulk and metasurfaces at the free surfaces. n nHere we carefully analyse this canonical problem by extracting the dispersion relations analytically thereby examining the influence of both the flexural and compressional resonances on the propagating wave. For an array of resonators atop an elastic half-space we augment the analysis with numerical simulations. Amongst other effects, we demonstrate the striking effect of a dispersion curve which corresponds to a mode that transitions from Rayleigh wave-like to shear wave-like behaviour and the resultant change in the fields from surface to bulk waves.

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–10u2009Hz, and for low-speed sedimentary soil (vs: 300–500u2009m/s), that the vibration of a structure is reduced by up to 6u2009dB at its resonance frequency.

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.

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–600u2009kHz, 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.

The validity of Kirchhoff theory for scattering of elastic waves from rough surfaces

The Kirchhoff approximation (KA) for elastic wave scattering from two-dimensional (2D) and three-dimensional (3D) rough surfaces is critically examined using finite-element (FE) simulations capable of extracting highly accurate data while retaining a fine-scale rough surface. The FE approach efficiently couples a time domain FE solver with a boundary integration method to compute the scattered signals from specific realizations of rough surfaces. Multiple random rough surfaces whose profiles have Gaussian statistics are studied by both Kirchhoff and FE models and the results are compared; Monte Carlo simulations are used to assess the comparison statistically. The comparison focuses on the averaged peak amplitude of the scattered signals, as it is an important characteristic measured in experiments. Comparisons, in both two dimensions and three dimensions, determine the accuracy of Kirchhoff theory in terms of an empirically estimated parameter σ2/λ0 (σ is the RMS value, and λ0 is the correlation length, of the roughness), being considered accurate when this is less than some upper bound c, (σ2/λ0<c). The incidence and scattering angles also play important roles in the validity of the Kirchhoff theory and it is found that for modest incidence angles of less than 30°, the accuracy of the KA is improved even when σ2/λ0>c. In addition, the evaluation results are compared using 3D isotropic rough surfaces and 2D surfaces with the same surface parameters.

A time-domain finite element boundary integration method for ultrasonic nondestructive evaluation.

A 2-D and 3-D numerical modeling approach for calculating the elastic wave scattering signals from complex stress-free defects is evaluated. In this method, efficient boundary integration across the complex boundary of the defect is coupled with a time-domain finite element (FE) solver. The model is designed to simulate time-domain ultrasonic nondestructive evaluation in bulk media. This approach makes use of the hybrid concept of linking a local numerical model to compute the near-field scattering behavior and theoretical mathematical formulas for postprocessing to calculate the received signals. It minimizes the number of monitoring signals from the FE calculation so that the computation effort in postprocessing decreases significantly. In addition, by neglecting the conventional regular monitoring box, the region for FE calculation can be made smaller. In this paper, the boundary integral method is implemented in a commercial FE code, and it is validated by comparing the scattering signals with results from corresponding full FE models. The coupled method is then implemented in real inspection scenarios in both 2-D and 3-D, and the accuracy and the efficiency are demonstrated. The limitations of the proposed model and future works are also discussed.

Unveiling Extreme Anisotropy in Elastic Structured Media

Periodic structures can be engineered to exhibit unique properties observed at symmetry points, such as zero group velocity, Dirac cones, and saddle points; identifying these and the nature of the associated modes from a direct reading of the dispersion surfaces is not straightforward, especially in three dimensions or at high frequencies when several dispersion surfaces fold back in the Brillouin zone. A recently proposed asymptotic high-frequency homogenization theory is applied to a challenging time-domain experiment with elastic waves in a pinned metallic plate. The prediction of a narrow high-frequency spectral region where the effective medium tensor dramatically switches from positive definite to indefinite is confirmed experimentally; a small frequency shift of the pulse carrier results in two distinct types of highly anisotropic modes. The underlying effective equation mirrors this behavior with a change in form from elliptic to hyperbolic exemplifying the high degree of wave control available and the importance of a simple and effective predictive model.

An asymptotic theory for waves guided by diffraction gratings or along microstructured surfaces

An effective surface equation, that encapsulates the detail of a microstructure, is developed to model microstructured surfaces. The equations deduced accurately reproduce a key feature of surface wave phenomena, created by periodic geometry, that are commonly called Rayleigh–Bloch waves, but which also go under other names, for example, spoof surface plasmon polaritons in photonics. Several illustrative examples are considered and it is shown that the theory extends to similar waves that propagate along gratings. Line source excitation is considered, and an implicit long-scale wavelength is identified and compared with full numerical simulations. We also investigate non-periodic situations where a long-scale geometrical variation in the structure is introduced and show that localized defect states emerge which the asymptotic theory explains.