Graded Quasiperiodic Metamaterials Perform Fractal Rainbow Trapping.

The rainbow trapping phenomenon of graded metamaterials can be combined with the fractal spectra of quasiperiodic waveguides to give a metamaterial that performs fractal rainbow trapping. This is achieved through a graded cut-and-project algorithm that yields a geometry for which the effective projection angle is graded along its length. As a result, the fractal structure of local band gaps varies with position, leading to broadband "fractal" rainbow trapping. We demonstrate this principle by designing an acoustic waveguide, which is characterised using theory, simulation and experiments.

Erratum: Elastic Orbital Angular Momentum [Phys. Rev. Lett. 128, 064301 (2022)].

This corrects the article DOI: 10.1103/PhysRevLett.128.064301.

Elastic Orbital Angular Momentum.

We identify that flexural guided elastic waves in elastic pipes carry a well-defined orbital angular momentum associated with the compressional dilatational potential. This enables the transfer of elastic orbital angular momentum, that we numerically demonstrate, through the coupling of the compressional potential in a pipe to the acoustic pressure field in a surrounding fluid in contact with the pipe.

Low-frequency wave-energy amplification in graded two-dimensional resonator arrays.

Energy amplification in square-lattice arrays of C-shaped low-frequency resonators, where the resonator radii are graded with distance, is investigated in the two-dimensional linear acoustics setting for both infinite (in one dimension) and finite arrays. Large amplifications of the incident energy are shown in certain array locations. The phenomenon is analysed using: (i) band diagrams for doubly-periodic arrays; (ii) numerical simulations for infinite and finite arrays; and (iii) eigenvalue analysis of transfer matrices operating over individual columns of the array. It is shown that the locations of the large amplifications are predicted by propagation cut-offs in the modes associated with the transfer-matrix eigenvalues. For the infinite array, the eigenvalues form a countable set, and for the low frequencies considered, only a single propagating mode exists for a given incident wave, which cuts off within the array, leading to predictive capabilities for the amplification location. For the finite array, it is shown that (in addition to a continuous spectrum of modes) multiple discrete propagating modes can be excited, with the grading generating new modes, as well as cutting others off, leading to complicated amplification patterns. The numerical simulations reveal that the largest amplifications are achieved for a single row array, with amplifications an order of magnitude smaller for the corresponding infinite array. This article is part of the theme issue 'Modelling of dynamic phenomena and localization in structured media (part 1)'.

Tailoring drug release rates in hydrogel-based therapeutic delivery applications using graphene oxide.

Graphene oxide (GO) is increasingly used for controlling mass diffusion in hydrogel-based drug delivery applications. On the macro-scale, the density of GO in the hydrogel is a critical parameter for modulating drug release. Here, we investigate the diffusion of a peptide drug through a network of GO membranes and GO-embedded hydrogels, modelled as porous matrices resembling both laminated and 'house of cards' structures. Our experiments use a therapeutic peptide and show a tunable nonlinear dependence of the peptide concentration upon time. We establish models using numerical simulations with a diffusion equation accounting for the photo-thermal degradation of fluorophores and an effective percolation model to simulate the experimental data. The modelling yields an interpretation of the control of drug diffusion through GO membranes, which is extended to the diffusion of the peptide in GO-embedded agarose hydrogels. Varying the density of micron-sized GO flakes allows for fine control of the drug diffusion. We further show that both GO density and size influence the drug release rate. The ability to tune the density of hydrogel-like GO membranes to control drug release rates has exciting implications to offer guidelines for tailoring drug release rates in hydrogel-based therapeutic delivery applications.

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.

The symmetry and coupling properties of solutions in general anisotropic multilayer waveguides.

Multilayered plate and shell structures play an important role in many engineering settings where, for instance, coated pipes are commonplace such as in the petrochemical, aerospace, and power generation industries. There are numerous demands, and indeed requirements, on nondestructive evaluation (NDE) to detect defects or to measure material properties using guided waves; to choose the most suitable inspection approach, it is essential to know the properties of the guided wave solutions for any given multilayered system and this requires dispersion curves computed reliably, robustly, and accurately. Here, the circumstances are elucidated, and possible layer combinations, under which guided wave solutions, in multilayered systems composed of generally anisotropic layers in flat and cylindrical geometries, have specific properties of coupling and parity; the partial wave decomposition of the wave field is utilised to unravel the behaviour. A classification into five families is introduced and the authors claim that this is the fundamental way to approach generally anisotropic waveguides. This coupling and parity provides information to be used in the design of more efficient and robust dispersion curve tracing algorithms. A critical benefit is that the analysis enables the separation of solutions into categories for which dispersion curves do not cross; this allows the curves to be calculated simply and without ambiguity.

A Langevin model for fluctuating contact angle behaviour parametrised using molecular dynamics.

Molecular dynamics simulations are employed to develop a theoretical model to predict the fluid-solid contact angle as a function of wall-sliding speed incorporating thermal fluctuations. A liquid bridge between counter-sliding walls is studied, with liquid-vapour interface-tracking, to explore the impact of wall-sliding speed on contact angle. The behaviour of the macroscopic contact angle varies linearly over a range of capillary numbers beyond which the liquid bridge pinches off, a behaviour supported by experimental results. Nonetheless, the liquid bridge provides an ideal test case to study molecular scale thermal fluctuations, which are shown to be well described by Gaussian distributions. A Langevin model for contact angle is parametrised to incorporate the mean, fluctuation and auto-correlations over a range of sliding speeds and temperatures. The resulting equations can be used as a proxy for the fully-detailed molecular dynamics simulation allowing them to be integrated within a continuum-scale solver.

Whispering Bloch modes.

We investigate eigenvalue problems for the planar Helmholtz equation in open systems with a high order of rotational symmetry. The resulting solutions have similarities with the whispering gallery modes exploited in photonic micro-resonators and elsewhere, but unlike these do not necessarily require a surrounding material boundary, with confinement instead resulting from the geometry of a series of inclusions arranged in a ring. The corresponding fields exhibit angular quasi-periodicity reminiscent of Bloch waves, and hence we refer to them as whispering Bloch modes (WBMs). We show that if the geometry of the system is slightly perturbed such that the rotational symmetry is broken, modes with asymmetric field patterns can be observed, resulting in field enhancement and other potentially desirable effects. We investigate the WBMs of two specific geometries first using expansion methods and then by applying a two-scale asymptotic scheme.

Dynamic interfacial trapping of flexural waves in structured plates.

The paper presents new results on the localization and transmission of flexural waves in a structured plate containing a semi-infinite two-dimensional array of rigid pins. In particular, localized waves are identified and studied at the interface boundary between the homogeneous part of the flexural plate and the part occupied by rigid pins. A formal connection has been made with the dispersion properties of flexural Bloch waves in an infinite doubly periodic array of rigid pins. Special attention is given to regimes corresponding to standing waves of different types as well as Dirac-like points that may occur on the dispersion surfaces. A single half-grating problem, hitherto unreported in the literature, is also shown to bring interesting solutions.

The effect of adsorption kinetics on the rate of surfactant-enhanced spreading.

A comparison of the kinetics of spreading of aqueous solutions of two different surfactants on an identical substrate and their short time adsorption kinetics at the water/air interface has shown that the surfactant which adsorbs slower provides a higher spreading rate. This observation indicates that Marangoni flow should be an important part of the spreading mechanism enabling surfactant solutions to spread much faster than pure liquids with comparable viscosities and surface tensions.

Infinite product expansion of the Fokker-Planck equation with steady-state solution.

We present an analytical technique for solving Fokker-Planck equations that have a steady-state solution by representing the solution as an infinite product rather than, as usual, an infinite sum. This method has many advantages: automatically ensuring positivity of the resulting approximation, and by design exactly matching both the short- and long-term behaviour. The efficacy of the technique is demonstrated via comparisons with computations of typical examples.

Modeling guided elastic waves in generally anisotropic media using a spectral collocation method.

Guided waves are now well established for some applications in the non-destructive evaluation of structures and offer potential for deployment in a vast array of other cases. For their development, it is important to have reliable and accurate information about the modes that propagate for particular waveguide structures. Essential information that informs choices of mode transducer, operating frequencies, and interpretation of signals, among other issues, is provided by the dispersion curves of different modes within various combinations of geometries and materials. In this paper a spectral collocation method is successfully used to handle the more complicated and realistic waveguide problems that are required in non-destructive evaluation; many pitfalls and limitations found in root-finding routines based on the partial wave method are overcome by using this approach. The general cases presented cover anisotropic homogeneous perfectly elastic materials in flat and cylindrical geometry. Non-destructive evaluation applications include complex waveguide structures, such as single or multi-layered fiber composites, lined, bonded and buried structures. For this reason, arbitrarily multi-layered systems with both solid and fluid layers are also addressed as well as the implementation of interface models of imperfect boundary conditions between layers.

Rayleigh-Bloch waves along elastic diffraction gratings.

Rayleigh-Bloch (RB) waves in elasticity, in contrast to those in scalar wave systems, appear to have had little attention. Despite the importance of RB waves in applications, their connections to trapped modes and the ubiquitous nature of diffraction gratings, there has been no investigation of whether such waves occur within elastic diffraction gratings for the in-plane vector elastic system. We identify boundary conditions that support such waves and numerical simulations confirm their presence. An asymptotic technique is also developed to generate effective medium homogenized equations for the grating that allows us to replace the detailed microstructure by a continuum representation. Further numerical simulations confirm that the asymptotic scheme captures the essential features of these waves.

The speed of an inclined ruck.

Steady rucks in an elastic beam can roll at constant speed down an inclined plane. We examine the dynamics of these travelling-wave structures and argue that their speed can be dictated by a combination of the physical conditions arising in the vicinity of the 'contact points' where the beam is peeled off the underlying plane and stuck back down. We provide three detailed models for the contact dynamics: viscoelastic fracture, a thermodynamic model for bond formation and detachment and adhesion mediated by a thin liquid film. The results are compared with experiments.

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.

Asymptotics for metamaterials and photonic crystals.

Metamaterial and photonic crystal structures are central to modern optics and are typically created from multiple elementary repeating cells. We demonstrate how one replaces such structures asymptotically by a continuum, and therefore by a set of equations, that captures the behaviour of potentially high-frequency waves propagating through a periodic medium. The high-frequency homogenization that we use recovers the classical homogenization coefficients in the low-frequency long-wavelength limit. The theory is specifically developed in electromagnetics for two-dimensional square lattices where every cell contains an arbitrary hole with Neumann boundary conditions at its surface and implemented numerically for cylinders and split-ring resonators. Illustrative numerical examples include lensing via all-angle negative refraction, as well as omni-directive antenna, endoscope and cloaking effects. We also highlight the importance of choosing the correct Brillouin zone and the potential of missing interesting physical effects depending upon the path chosen.

Nonequilibrium hysteresis and Wien effect water dissociation at a bipolar membrane.

As in electrochemical cyclic voltammetry, time-periodic reverse voltage bias across a bipolar membrane is shown to exhibit hysteresis due to transient effects. This is due to the incomplete depletion of mobile ions, at the junction between the membranes, within two adjoining polarized layers; the layer thickness depends on the applied voltage and the surface charge densities. Experiments show that the hysteresis consists of an Ohmic linear rise in the total current with respect to the voltage, followed by a decay of the current. A limiting current is established for a long period when all the mobile ions are depleted from the polarized layer. If the resulting high field within the two polarized layers is sufficiently large, water dissociation occurs to produce proton and hydroxyl traveling wave fronts which contribute to another large jump in the current. We use numerical simulation and asymptotic analysis to interpret the experimental results and to estimate the amplitude of the transient hysteresis and the water-dissociation current.

Dynamics of a climbing surfactant-laden film--I: base-state flow.

The dynamics of a surfactant-laden film climbing up an inclined plane is investigated through a two-dimensional (2-D), nonlinear evolution equation for the interface coupled to convective-diffusion equations for the surfactant, derived using lubrication theory. One-dimensional (1-D) solutions, representing the base-state flow, are investigated for constant flux and constant volume configurations; these flows are parameterised by capillarity, gravity, convection-diffusion ratios (represented by Péclét numbers at the surface and bulk), a solubility parameter, sorption kinetics constants, the number of surfactant monomers in a micelle, and the nonlinearity of the surfactant equation of state. In both configurations studied, a front develops spreading up the substrate against the direction of gravity whereby the leading edge of the front follows a power-law as a function of time. The effect of system parameters on the base-state flow is explored through an extensive parametric study, while the stability of the above-mentioned system to spanwise perturbations is the focus of Part II.

Dynamics of a climbing surfactant-laden film II: stability.

The linear and nonlinear stability of a spreading film of constant flux and a drop of constant volume, discussed in [1], are examined here. A linear stability analysis (LSA) is carried out to investigate the stability to spanwise perturbations, by linearisation of the two-dimensional (2-D) evolution equations derived in [1] for the film thickness and surfactant concentration fields. The latter correspond to convective-diffusion equations for the surfactant, existing in the form of monomers (present at the free surface and in the bulk) and micelles (present in the bulk). The results of the LSA indicate that the thinning region, present upstream of the leading front in the constant flux case, and the leading ridge in the constant volume case, are unstable to spanwise perturbations. Numerical simulations of the 2-D system of equations demonstrate that the above-mentioned regions exhibit finger formation; the effect of selected system parameters on the fingering patterns is discussed.