Primary wave transmission in systems of elastic rods with granular interfaces

We study analytically and numerically primary pulse transmission in one dimensional systems of identical linearly elastic non-dispersive rods separated by identical homogeneous granular layers composed of n beads. The beads interact elastically through a strongly (essentially) nonlinear Hertzian contact law. The main challenge in studying pulse transmission in such strongly nonlinear media is to analyze the ‘basic problem’, namely, the dynamical response of a single intermediate granular layer, confined from both ends by barely touching linear elastic rods subject to impulsive excitation of the left rod. The analysis of the basic problem is carried out under two basic assumptions; namely, of sufficiently small duration of the shock excitation applied to the first layer of the system, and of sufficiently small mass of each bead in the granular interface compared to the mass of each rod. In fact, the smallness of the mass of the bead defines the small parameter in the asymptotic analysis of this problem. Both assumptions are reasonable from the point of view of practical applications. In the analysis we focus only in primary pulse propagation, by neglecting secondary pulse reflections caused by wave scattering at each granular interface and considering only the transmission of the main (primary) pulse across the interface to the neighboring elastic rod. Two types of shock excitations are considered. The first corresponds to fixed time duration (but still much smaller compared to the characteristic time of pulse propagation through the length of each rod), whereas the second type corresponds to a pulse duration that depends on the small parameter of the problem. The influence of the number of beads of the granular interface on the primary wave transmission is studied, and it is shown that at granular interfaces with a relatively low number of beads fast time scale oscillations are excited with increasing amplitudes with increasing number of beads. For a larger number of beads, primary pulse transmission is by means of solitary wave trains resulting from the dispersion of the original shock pulse; in that case fast oscillations result due to interference phenomena caused by the scattering of the main pulse at the boundary of the interface. Considering a periodic system of rods we demonstrate significant reduction of the primary pulse when transmitted through a sequence of granular interfaces. This result highlights the efficacy of applying granular interfaces for passive shock mitigation in layered elastic media.     

Asymptotic approximations for Bloch waves and topological mode steering in a planar array of Neumann scatterers

We study the canonical problem of wave scattering by periodic arrays, either of infinite or finite extent, of Neumann scatterers in the plane; the characteristic lengthscale of the scatterers is considered small relative to the lattice period. We utilise the method of matched asymptotic expansions, together with Fourier series representations, to create an efficient and accurate numerical approach for finding the dispersion curves associated with Floquet–Bloch waves through an infinite array of scatterers. The approach lends itself to direct scattering problems for finite arrays and we illustrate the flexibility of these asymptotic representations on topical examples from topological wave physics.     

Approximation by Multipoles of the Multiple Acoustic Scattering by Small Obstacles in Three Dimensions and Application to the Foldy Theory of Isotropic Scattering

The asymptotic analysis carried out in this paper for the problem of a multiple scattering in three dimensions of a time-harmonic wave by obstacles whose size is small as compared with the wavelength establishes that the effect of the small bodies can be approximated at any order of accuracy by the field radiated by point sources. Among other issues, this asymptotic expansion of the wave furnishes a mathematical justification with optimal error estimates of Foldy’s method that consists in approximating each small obstacle by a point isotropic scatterer. Finally, it is shown how this theory can be further improved by adequately locating the center of phase of the point scatterers and the taking into account of self-interactions. In this way, it is established that the usual Foldy model may lead to an approximation whose asymptotic behavior is the same than that obtained when the multiple scattering effects are completely neglected.     

PROPAGATION OF A LONG WAVE IN A CURVED DUCT (Ⅱ)APPLICATIONS OF MATCHED EXPANSION TO LONG WAVE PROPAGATION THROUGH A HOLE WITH VARIABLE CROSS-SECTION

Only the case in which the parameterε=ka《1 is considered in this paper,where k is the wave number and a is the characteris-tic radius of the cross-section of the hole.The general asymp-totic expansion of the complex velocity potential of a long wave propagating in the hole with variable cross-section is obtained by regular perturbation:The methods of matched asymptotic ex-pansion are employed to calculate the reflection coefficients,scattering coefficients and radiation coefficients at the open ends of the hole when a long wave propagates through it,which may be open at both ends or only at one end.Three examples of different kinds of holes are given to show the way to solve such two-dimensional or three-dimensional problems.     

Wave propagation with tunneling in a highly discontinuous layered medium

An impulsive plane wave traverses a stratified medium consisting of a large number N of homogeneous isotropic perfectly elastic layers. The directly transmitted wave is greatly reduced by the cumulative effect of scattering loss at each of the many interfaces. However, close to the arrival of the direct wave is a broad pulse, arising from multiple scattering; this pulse does not decay as rapidly as the direct wave and ultimately appears to diffuse about a moving center. The latter process, which is determined by the medium statistics, leads to time delays, effective anisotropy, and apparent attenuation.

The present work may be regarded as an extension of that described by Burridge, White and Papanicolaou (1988) and Burridge and Chang (1989) to allow for tunneling P waves for S-wave incidence beyond the critical angle.

When the reflection coefficients at the interfaces are scaled as 1/√N while N → ∞, and when time is measured in units of vertical travel time across an average layer, numerical solutions of the exact problem show that the shape of the broad transmitted pulse approaches the limiting form given as the solution of a certain integrodifferential equation in accordance with our asymptotic theory.     

Diffraction of solitary water waves around three-dimensional floating bodies

By using the matched asymptotic expansion method and the idea of edge layer, a mathematic model for describing the interaction between weakly nonlinear shallow-water waves and three-dimensional floating bodies is formed in the paper. As a numerical example, the diffraction of a solitary wave around a vertically floating circular cylinder has been investigated and the results are presented. The present method can further be extended to the study of wave diffraction around floating bodies of general shape. The project is supported by the National Natural Science Foundation of China.     

各向异性平板开孔动应力集中问题的研究

采用各向异性平板弯曲波动理论及摄动方法，对正交各向异性平板开孔弯曲波的散射及动应力集中问题进行了分析研究，得到了此种平板稳态弯曲波动问题的渐近形式的分析解。同时采用保角映射技术，为求解正交各向异性平板开孔弹性波的散射及动应力集中问题提供了一种统一规范的方法。     

Scattering by penetrable acoustic targets

An acoustic target of constant density ?_t and variable index of refraction is imbedded in a surrounding acoustic fluid of constant density ?_a. A time harmonic wave propagating in the surrounding fluid is incident on the target. We consider two limiting cases of the target where the parameter ε ≡ ?_a/?_t → 0 (the nearly rigid target) or ε → ∞ (the nearly soft target). Wh en the frequency of the incident wave is bounded away from the ‘in-vacuo’ resonant frequencies of the target, the resulting scattered field is essentially the field scattered by the rigid target for ε = 0 or the soft target if ε → ∞. However, when the frequency of the incident wave is near a resonant frequency,the target oscillates and its interaction with the surrounding fluid produces peaks in the scattered field amplitude. In this paper we obtain asymptotic expansions of the solutions of the scattering problems for the nearly rigid and the nearly soft targets as ε → 0 or ε → ∞, respectively, that are uniformly valid in the incident frequency. The method of matched asymptotic expansions is used in the analysis. The outer and inner expansions correspond to the incident frequencies being far or near to the resonant frequencies, respectively. We have applied the method only to simple resonant frequencies, but it can be extended to multiple resonant frequencies. The method is applied to the incidence of a plane wave on a nearly rigid sphere of constant index of refraction. The far field expressions for the scattered fields, including the total scattering cross-sections, that are obtained from the asymptotic method and from the partial wave expansion of the solution are in close agreement for sufficiently small values of ε.     

Interactions of propagating waves in a one-dimensional chain of linear oscillators with a strongly nonlinear local attachment

We study the interaction of propagating wavetrains in a one-dimensional chain of coupled linear damped oscillators with a strongly nonlinear, lightweight, dissipative local attachment which acts, in essence, as nonlinear energy sink—NES. Both symmetric and highly un-symmetric NES configurations are considered, labelled S-NES and U-NES, respectively, with strong (in fact, non-linearizable or nearly non-linearizable) stiffness nonlinearity. Especially for the case of U-NES we show that it is capable of effectively arresting incoming slowly modulated pulses with a single fast frequency by scattering the energy of the pulse to a range of frequencies, by locally dissipating a major portion of the incoming energy, and then by backscattering residual waves upstream. As a result, the wave transmission past the location of the NES is minimized, and the NES acts, in effect, as passive wave arrestor and reflector. Analytical reduced-order modeling of the dynamics is performed through complexification/averaging. In addition, governing nonlinear dynamics is studied computationally and compared to the analytical predictions. Results from the reduced order model recover the exact computational simulations.     

Equivalent multipolar point-source modeling of small spheres for fast and accurate electromagnetic wave scattering computations

In this paper, we develop reduced models to approximate the solution of the electromagnetic scattering problem in an unbounded domain which contains a small perfectly conducting sphere. Our approach is based on the method of matched asymptotic expansions. This method consists in defining an approximate solution using multi-scale expansions over outer and inner fields related in a matching area. We make explicit the asymptotics up to the second order of approximation for the inner expansion and up to the fifth order for the outer expansion. We validate the results with numerical experiments which illustrate theoretical orders of convergence for the asymptotic models requiring negligible computational cost.     

Wave radiation and diffraction from a small submerged circular cylinder

The linear water wave radiation/diffraction problem for a small submerged cylinder is solved analytically in terms of an asymptotic expansion. The expansion parameter is the radius of the cylinder divided by its depth of submergence. The influence of the cylinder on the wave field is represented by a rotating dipole, which obeys the free-surface conditions and radiation conditions. Our results agree fully with the classical results. We find a basic physical reason why the circular shape is exceptional compared to all other obstacle geometries: it is because of the circular paths of fluid particles in free linear deep-water waves.     

Mass,momentum and total excess energy transported by a weak planeN wave

A weakly nonlinear plane acoustic wave is emitted into an ideal gas of semi-infinite extent from an infinite plate by its sinusoidal motion of single period. The wave develops into anN wave in the far field, as long as the energy dissipation is negligible everywhere except for discontinuous shock fronts. The third-order effects at shock fronts are evaluated, i.e., the generation of reflected acoustic wave as a result of the interaction of shock and expansion wave and the production of entropy by the energy dissipation at shock fronts. Consideration of these effects enables one to estimate the whole mass, momentum and total excess energy (sum of the kinetic energy and excess of internal energy over an initial undisturbed value) transported by theN wave to the accuracy of third order of wave amplitude. It is shown that the mass and total excess energy transported by theN wave increase and the momentum decreases to asymptotic limits as the wave propagates. The result shows good agreement with a numerical result obtained by solving the Euler equations with a high-resolution TVD upwind scheme.     

Perturbation solutions for asymmetric laminar flow in porous channel with expanding and contracting walls

The cases of large Reynolds number and small expansion ratio for the asym- metric laminar flow through a two-dimensional porous expanding channel are considered. The Navier-Stokes equations are reduced to a nonlinear fourth-order ordinary differential equation by introducing a time and space similar transformation. A singular perturbation method is used for the large suction Reynolds case to obtain an asymptotic solution by matching outer and inner solutions. For the case of small expansion ratios, we are able to obtain asymptotic solutions by double parameter expansion in either a small Reynolds number or a small asymmetric parameter. The asymptotic solutions indicate that the Reynolds number and expansion ratio play an important role in the flow behavior. Nu- merical methods are also designed to confirm the correctness of the present asymptotic solutions.     

Reflection and transmission of a compression wave at a tunnel portal

Analyses are made of the interaction of the nonlinearly steepened, compression wavefront generated by a high-speed train in a tunnel with the tunnel portal ahead of the train. The ‘micro-pressure’ pulse emitted from the portal can rattle structures in nearby buildings, and the expansion wave reflected back towards the train can cause discomfort to passengers. It is concluded that the usual simplified approximation of one-dimensional propagation within the tunnel provides an adequate representation of interactions of the wave with the portal, and also with ‘windows’ in the tunnel wall near the portal. It is shown how a discrete distribution of windows can be used to produce a reflected expansion wave that varies linearly across the wavefront, and how the thickness of that wavefront can be made many times larger than the thickness of the incident compression wave profile. A detailed analysis of the wave radiated from the portal reveals that cumulative nonlinear effects of propagation over long distances make little or no contribution to the free-space radiation of the micro-pressure wave.     

Spectral method for multiple edge diffraction by a flat strip

A spectral iteration scheme is employed to analyze time-harmonic and transient scattering of an E- or H-polarized incident plane wave by a perfectly conducting plane strip. The scattered field is synthesized by successive interactions between the edges, with each interaction modeled by half-plane diffraction. The plane wave spectrum generating a particular order of diffraction consists, in addition to the incident plane wave excitation, of a portion determined from the previous diffraction. The multiple integral spectral representations constructed in this manner satisfy the edge condition, and they are in a form suitable for inversion into the time domain by the modified Cagniard-deHoop method. Asymptotic reductions for special cases yield agreement with results from other methods, when available. Numerical calculations including up to triple diffraction have been performed for H- and E-polarized impulse and Gaussian pulse scattering. The results are clearly seen to repair the deficiencies of wavefront approximations at longer observation times, and from comparison with data generated independently by eigenfunction expansion, they describe accurately the total scattered response, owing to the high damping rate of higher-order diffractions.     

Designing a Linear Structure with a Local Nonlinear Attachment for Enhanced Energy Pumping

We present a design procedure for enhancing nonlinear energy pumping from a mode of a linear-damped substructure to a weakly coupled, essentially nonlinear oscillator. By this we denote the one way, irreversible passive transfer of vibrational energy from the mode to the nonlinear attachment. The design relies in the asymptotic expansion for large energies of a nonlinear normal mode of the underlying conservative system that provides an analytic estimate of the level of the amplitude reached by the nonlinear attachment in the energy pumping regime. The analytical findings are validated by direct numerical simulations.     

Efficient computation of steady solitary gravity waves

An efficient numerical method to compute solitary wave solutions to the free surface Euler equations is reported. It is based on the conformal mapping technique combined with an efficient Fourier pseudo-spectral method. The resulting nonlinear equation is solved via the Petviashvili iterative scheme. The computational results are compared to some existing approaches, such as Tanaka’s method and Fenton’s high-order asymptotic expansion. Several important integral quantities are computed for a large range of amplitudes. The integral representation of the velocity and acceleration fields in the bulk of the fluid is also provided.     

A justification of nonlinear properly invariant plate theories

A single asymptotic derivation of three classical nonlinear plate theories is presented in a setting which preserves the frame-invariance properties of three-dimensional finite elasticity. By a successive scaling of the external loading on the three-dimensional body, the nonlinear membrane theory, the nonlinear inextensional theory and the von Kármán equations are derived as the leading-order terms in the asymptotic expansion of finite elasticity. The governing equations of the nonlinear inextensional theory are of particular interest where 1) plane-strain kinematics and plane-stress constitutive equations are derived simultaneously from the asymptotic analysis, 2) the theory can be phrased as a minimization problem over the space of isometric deformations of a surface, and 3) the local equilibrium equations are identical to those arising in the one-director Cosserat shell model. Furthermore, it can be concluded that with a regular, single-scale asymptotic expansion it is not possible to obtain a system of plate equations in which finite membrane strain and finite bending strain occur simultaneously in the leading-order term of an asymptotic analysis.