ABSTRACT
In this work, the concept of high-frequency homogenization is extended to the case of one-dimensional periodic media with imperfect interfaces of the spring-mass type. In other words, when considering the propagation of elastic waves in such media, displacement and stress discontinuities are allowed across the borders of the periodic cell. As is customary in high-frequency homogenization, the homogenization is carried out about the periodic and antiperiodic solutions corresponding to the edges of the Brillouin zone. Asymptotic approximations are provided for both the higher branches of the dispersion diagram (second-order) and the resulting wave field (leading-order). The special case of two branches of the dispersion diagram intersecting with a non-zero slope at an edge of the Brillouin zone (occurrence of a so-called Dirac point) is also considered in detail, resulting in an approximation of the dispersion diagram (first-order) and the wave field (zeroth-order) near these points. Finally, a uniform approximation valid for both Dirac and non-Dirac points is provided. Numerical comparisons are made with the exact solutions obtained by the Bloch-Floquet approach for the particular examples of monolayered and bilayered materials. In these two cases, convergence measurements are carried out to validate the approach, and we show that the uniform approximation remains a very good approximation even far from the edges of the Brillouin zone.
ABSTRACT
A continued fraction expansion to the immittances defining viscothermal wave propagation in a cylindrical tube has been presented recently in this journal, intended as a starting point for time domain numerical method design. Though the approximation has the great benefit of passivity, or positive realness under truncation, its convergence is slow leading to approximations of high order in practice. Other passive structures, when combined with optimisation methods, can lead to good accuracy over a wide frequency range, and for relatively low order.
ABSTRACT
An explicit finite-difference scheme is presented for solving the two-dimensional Biot equations of poroelasticity across the full range of frequencies. The key difficulty is to discretize the Johnson-Koplik-Dashen (JKD) model which describes the viscous dissipations in the pores. Indeed, the time-domain version of Biot-JKD model involves order 1/2 fractional derivatives which amount to a time convolution product. To avoid storing the past values of the solution, a diffusive representation of fractional derivatives is used: The convolution kernel is replaced by a finite number of memory variables that satisfy local-in-time ordinary differential equations. The coefficients of the diffusive representation follow from an optimization procedure of the dispersion relation. Then, various methods of scientific computing are applied: The propagative part of the equations is discretized using a fourth-order finite-difference scheme, whereas the diffusive part is solved exactly. An immersed interface method is implemented to discretize the geometry on a Cartesian grid, and also to discretize the jump conditions at interfaces. Numerical experiments are proposed in various realistic configurations.