Consider the scattering of a time-harmonic elastic plane wave by a periodic rigid surface. The elastic wave propagation is governed by the two-dimensional Navier equation. Based on a Dirichlet-to-Neumann (DtN) map, a transparent boundary condition (TBC) is introduced to reduce the scattering problem into a boundary value problem in a bounded domain. By using the finite element method, the discrete problem is considered, where the TBC is replaced by the truncated DtN map. A new duality argument is developed to derive the a posteriori error estimate, which contains both the finite element approximation error and the DtN truncation error. An a posteriori error estimate based adaptive finite element algorithm is developed to solve the elastic surface scattering problem. Numerical experiments are presented to demonstrate the effectiveness of the proposed method.
more »
« less
Elastic scattering from rough surfaces in three dimensions
Consider the elastic scattering of a plane or point incident wave by an unbounded and rigid rough surface. The angular spectrum representation (ASR) for the time-harmonic Navier equation is derived in three dimensions. The ASR is utilized as a radiation condition to the elastic rough surface scattering problem. The uniqueness is proved through a Rellich-type identity for surfaces given by uniformly Lipschitz functions. In the case of flat surfaces with local perturbations, an equivalent variational formulation is deduced in a truncated bounded domain and the existence of solutions are shown for general incoming waves. The main ingredient of the proof is the radiating behavior of the Green tensor to the first boundary value problem of the Navier equation in a half-space.
more »
« less
- Award ID(s):
- 1912704
- NSF-PAR ID:
- 10182368
- Date Published:
- Journal Name:
- Journal of differential equations
- Volume:
- 269
- ISSN:
- 0022-0396
- Page Range / eLocation ID:
- 4045-4078
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The existing concepts of non-reciprocity in propagation of acoustic or elastic waves are based either on nonlinear effects, or on local circulation of linear elastic fluid that leads to red or blue Doppler shift, depending on the direction of sound wave. The same concepts exist for electromagnetic non-reciprocity, where external magnetic field may produce the effect similar to local rotation of the medium. These two concepts originate from two known methods of breaking a time-reversal symmetry (T-symmetry), that is necessary for observation of nonreciprocal wave propagation. Both concepts require additional electrical or mechanical devices to be installed with their own power sources. Here we propose to explore viscosity of fluid as a natural factor of the T-symmetry breaking through energy dissipation. We report experimental observation of the nonreciprocal transmission of ultrasound through a water-submerged phononic crystals consisting of several layers of aluminum rods arranged in a square lattice. While viscous losses break the T-symmetry, making the wave propagation thermodynamically irreversible, the transmission remains reciprocal if the scatterers are symmetrical. To generate different energy losses for opposite directions of propagation, the P-symmetry of the crystal is broken by using asymmetric scatterers. Due to asymmetry, two sound waves propagating in the opposite directions produce different distributions of velocity and pressure that leads to different local absorption. Dissipation of acoustic energy occurs mostly near the surface of the scatterers and it strongly depends on surface roughnesses. Using two phononic crystal with smooth and rough aluminum rods we demonstrate low (2-5 dB) and high (10-15 dB) level of non-reciprocity within a wide range of frequencies, 300-600 kHz. Experimental results are in agreement with numerical simulations based on the Navier-Stokes equation. This nonreciprocal linear device is very cheap, robust and does not require energy source.more » « less
-
The paper studies the equilibrium configurations of inextensible elastic membranes exhibiting lateral fluidity. Using a continuum description of the membrane's motions based on the surface Navier–Stokes equations with bending forces, the paper derives differential equations governing the mechanical equilibrium. The equilibrium conditions are found to be independent of lateral viscosity and relate tension, pressure, and tangential velocity of the fluid. These conditions suggest that either the lateral fluid motion ceases or non-decaying stationary flow of mass can only be supported by surfaces with Killing vector fields, such as axisymmetric shapes. A shape equation is derived that extends the classical Helfrich model with an area constraint to membranes of non-negligible mass. Furthermore, the paper suggests a simple numerical method to compute solutions of the shape equation. Numerical experiments conducted reveal a diverse family of equilibrium configurations. The stability of equilibrium states involving lateral flow of mass remains an unresolved question.more » « less
-
more » « less
Abstract We present new methodologies for the numerical solution of problems of elastic scattering by open arcs in two dimensions. The algorithms utilize weighted versions of the classical elastic integral operators associated with Dirichlet and Neumann boundary conditions, where the integral weight accounts for (and regularizes) the singularity of the integral‐equation solutions at the open‐arc endpoints. Crucially, the method also incorporates a certain “open‐arc elastic Calderón relation” introduced in this article, whose validity is demonstrated on the basis of numerical experiments, but whose rigorous mathematical proof is left for future work. (In fact, the aforementioned open‐arc elastic Calderón relation generalizes a corresponding elastic Calderón relation for closed surfaces, which is also introduced in this article, and for which a rigorous proof is included.) Using the open‐surface Calderón relation in conjunction with spectrally accurate quadrature rules and the Krylov‐subspace linear algebra solver GMRES, the proposed overall open‐arc elastic solver produces results of high accuracy in small number of iterations, for both low and high frequencies. A variety of numerical examples in this article demonstrate the accuracy and efficiency of the proposed methodology.
-
ABSTRACT The northern clingfish (Gobiesox maeandricus) has a suction-based adhesive disc that can stick to incredibly rough surfaces, a challenge for stiff commercial suction cups. Both clingfish discs and bioinspired suction cups have stiff cores but flexible edges that can deform to overcome surface irregularities. Compliant surfaces are common in nature and technical settings, but performance data for fish and commercial cups are gathered from stiff surfaces. We quantified the interaction between substrate compliance, surface roughness and suction performance for the northern clingfish, commercial suction cups and three biomimetic suction cups with disc rims of varying compliance. We found that all cups stick better on stiffer substrates and worse on more compliant ones, as indicated by peak stress values. On compliant substrates, surface roughness had little effect on adhesion, even for commercial cups that normally fail on hard, rough surfaces. We propose that suction performance on compliant substrates can be explained in part by effective elastic modulus, the combined elastic modulus from a cup–substrate interaction. Of all the tested cups, the biomimetic cups performed the best on compliant surfaces, highlighting their potential to be used in medical and marine geotechnical fields. Lastly, we discuss the overmolding technique used to generate the bioinspired cups and how it is an important tool for studying biology.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1016/j.jde.2020.03.022
