A novel design scheme for acoustic cloaking of complex shape based on region partitioning and multi-origin coordinate transformation

Acoustic cloaking is an important application of acoustic metamaterials. This article proposes a novel design scheme for acoustic cloaking based on the region partitioning and multi-origin coordinate transformation. The cloaked region is partitioned into multiple narrow strips. For each strip, a local coordinate system is established with the local origin located at the strip center, and a coordinate transformation in the local coordinate system is conducted to squeeze the material along the strip length direction to form the cloaked region. To facilitate the implementation of the acoustic cloak, the multilayer effective medium is used to approximate the non-uniform anisotropic material parameters. The effectiveness of the proposed coordinate transformation method is verified by comparing the results from our method with those in the literature. Firstly, the results of a circular acoustic cloak in the literature are reproduced by using our finite element (FE) simulations for validation. Then, a comparison is made between the traditional coordinate transformation scheme and our new scheme for simulating an elliptical acoustic cloak. The results indicate that the proposed multi-origin coordinate transformation method has a better cloaking effect on the incident wave along the ellipse minor axis direction than the traditional method. This means that for the same object, an appropriate transformation scheme can be selected for different incident wave directions to achieve the optimal control effect. The validated scheme is further used to design an arch-shaped cloak composed of an upper semicircular area and a lower rectangular area, by combining the traditional single-centered coordinate transformation method for the semicircular area and the proposed multi-origin method for the rectangular area. The results show that the designed cloak can effectively control the wave propagation with significantly reduced acoustic pressure level. This work provides a flexible acoustic cloak design method applicable for arbitrary shapes and different wave incident directions, enriching the theory of acoustic cloaking based on coordinate transformation.     

Simulation-based conceptual design of an acoustic metamaterial with full band gap using an air-based 1-3 piezoelectric composite for ultrasonic noise control

This paper aims at proposing a novel type of acoustic metamaterials with complete band gap composed of piezoelectric rods with square array as inclusions embedded in an air background (matrix). A modified plane wave expansion method accompanied with the principles of the Bloch–Floquet method with electromechanical coupling effect and also impedance spectra are used to get a band frequency and to investigate the passband for the selected cut of piezoelectric rods. We investigate both the electromechanical coupling coefficient and mechanical quality factor and their dependency to passband and bandwidth, which depends on both the density and the wave impedance of the matrix and the inclusions (rods). The ratio of the volume of inclusion to the matrix is used to define the fill factor or the so-called inclusion ratio, to introduce the bandwidth as a function of that. Furthermore, the fabrication method is presented in this paper. The results make a suitable foundation for design purposes and may develop an inherently passive ultrasonic noise control. In addition, the results provide the required guidance for a simulation-based design of elastic wave filters or wave guide that might be useful in high-precision mechanical systems operated in certain frequency ranges and switches made of piezoelectric materials; they also propose a novel type of elastic metamaterials, which is independent of the wave direction and has an equal sensitivity in all directions in which it reacts omnidirectionally and mitigates the occupational noise exposure.     

A THEORY OF EFFECTIVE THERMAL CONDUCTIVITY FOR MATRIX-INCLUSION- MICROCRACK THREE-PHASE HETEROGENEOUS MATERIALS BASED ON MICROMECHANICS

The effective thermal conductivity of matrix-inclusion-microcrack three-phase heterogeneous materials is investigated with a self-consistent micromechanical method (SCM) and a random microstructure finite element method(RMFEM). In the SCM, microcracks are assumed to be randomly distributed and penny-shaped and inclusions to be spherical, the crack effect is accounted for by introducing a crack density parameter, the effective thermal conductivity is derived which relates the macroscopic behavior to the crack density parameter. In the RMFEM, the highly irregular microstructure of the heterogeneous media is accurately described, the interaction among the matrix-inclusion-microcracks is exactly treated, the inclusion shape effect and crack size effect are considered. A Ni/ZrO₂ particulate composite material containing randomly distributed, penny-shaped cracks is examined as an example. The main results obtained are: (1) the effective thermal conductivity is sensitive to the crack density and exhibits essentially a linear relationship with the density parameter; (2) the inclusion shape has a significant effect on the effective thermal conductivity and a polygon-shaped inclusion is more effective in increasing or decreasing the effective thermal conductivity than a sphere-shaped one; and (3) the SCM and RMFEM are compared and the two methods give the same effective property in the case in which the matrix thermal conductivity λ₁ is greater than the inclusion one λ₂. In the inverse case of λ₁ < λ₂, the two methods agree as the inclusion volume fraction and crack density are low and differ as they are high. A reasonable explanation for the agreement and deviation between the two methods in the case of λ₁ < λ₂ is made. This work was supported by the National Natural Science Foundation of China and Chnese “863” High-Tech, Program.     

流固耦合地震波动问题的显式谱元模拟方法

流固耦合地震波动问题主要研究由流体和固体构成的复杂系统中地震波传播特性及其规律. 传统模拟方法中一般以声波方程、弹性波方程的数值解分别描述理想流体和弹性固体中的波动, 并实时地处理两种不同性质介质之间的相互耦合作用, 数值格式复杂且限制数值模拟精度与计算效率. 本文采用谱元法结合多次透射公式人工边界条件实现了一种流固耦合地震波动问题的高阶显式数值计算方法. 该方法利用了流固耦合问题统一计算框架,可将饱和多孔介质的Biot波动方程分别退化为理想流体的声波方程和弹性固体的弹性波方程. 通过P波垂直入射的水平成层理想流体-饱和多孔介质-弹性固体场地模型、P波斜入射的不规则层状界面以及任意形状界面的理想流体-饱和多孔介质-弹性固体场地模型等三个算例, 与传递函数法解析解以及集中质量有限元法计算结果进行对比分析, 证明了本文方法的正确性与有效性. 数值模拟结果表明, 本文方法相较传统有限元法可以少得多的节点数量获得更高的数值精度, 并且在较宽的频率范围内都能可靠地模拟出流固耦合系统的动力响应, 充分体现出本文方法兼顾高精度、计算效率和复杂场地建模灵活的特点.      

High accuracy wave simulation – Revised derivation,numerical analysis and testing of a nearly analytic integration discrete method for solving acoustic wave equation

The nearly analytic integration discrete (NAID) method for solving the two-dimensional acoustic wave equation has been fully mathematically revised, analyzed and tested. The NAID method is an alternative numerical modeling method for generating synthetic seismograms. The acoustic wave equation is first transformed into a system of first-order ordinary differential equations (ODEs) with respect to time variable t, and then directly integrated at a small time interval of [t_n, t_n+1] to obtain semi-discrete ordinary differential equations. The integral kernel is expanded into a truncated Taylor series, to which the integration operator is explicitly applied. The high-order temporal derivatives involved in the integral kernel are replaced by high-order spatial derivatives, which then are approximately calculated as a weighted linear combination of the displacement, the particle-velocity, and their spatial gradients. In this article, we investigate the theoretical properties of the revised NAID method, including the discrete error and the stability criteria. Numerical results for constant and layered velocity models show that, comparing to the Lax–Wendroff correction (LWC) scheme and the staggered-grid finite difference method, the NAID method can effectively suppress the numerical dispersion and source-noises caused by the discretization of the acoustic wave equation with too-coarse spatial grids or when models have strong velocity contrasts between adjacent layers. The proposed NAID method has been applied in computing the acoustic wavefields for two heterogeneous models – the corner edge model and the Marmousi model. Promising numerical results illustrate that the NAID method provides an encouraging tool for large-scale and complex wave simulation and inversion problems based on the acoustic equation.     

Direct numerical modeling of time-reversal acoustic subwavelength focusing

Focusing waves back to their original source position is possible both experimentally and numerically thanks to time reversal mirrors (TRM). For a TRM placed in the far field of the source, the focusing spot of the reversed wavefield is subject to the diffraction limit and cannot be smaller than half the minimum wavelength, even for a very small source. Yet, numerous time reversal experiments in resonating media have shown subwavelength focusing. In this work, we show that it is possible to model these subwavelength focusing observations with simple physics, only the 2-D standard acoustic wave equation, and with specific fine scale heterogeneity. Our work is based on the spectral element method to solve the wave equation and to model time reversal experiments. Such a method makes it possible to propagate very long time series in complex and strongly discontinuous media with high accuracy. The acoustic wave equations are solved at the fine scale in media with one or more split rings of size much smaller than the wavelength. Such split rings produce a Helmholtz resonance effect as well as propagation band-gaps. We show that, in such media, even with a single split ring resonator, subwavelength focusing down to 1/13th of the minimum wavelength can be observed.     

Scattering coefficient reconstruction in plates using focused acoustic beams

A method has been developed and proven using highly focused acoustic beams that allows for the rapid reconstruction of scattering coefficients of a thin anisotropic plate immersed in liquid. In a single bistatic coordinate scan, nearly the entire range of wave number within the spatial and temporal frequency bandwidth of the transducer can be reconstructed. This paper also reports the development of a multiple-source complex transducer point model that includes all extrinsic factors and permits prediction of the wave number-frequency (k–f) domain result obtained from a scan of focused transducers in a pitch-catch reflection or transmission arrangement. Extensive experiments have been performed to test the method and the model and to demonstrate transducer beam effects on the k–f domain results, leading to a very efficient method for mapping major portions of the guided wave dispersion spectrum in thin-plate media. As a demonstration of the technique, an estimate of material elastic properties in an isotropic and a transversely-isotropic plate is reported, making only minimal use of the highly redundant dispersion data. Acoustic velocities inferred from these experiments show a disparity of less than 3% from contact acoustic estimates of the same parameters in either plate.     

Interaction of a plane harmonic wave with a thin rigid inclusion of the shape of a cylindrical shell

The paper presents the solution of the problem of determining the stress state in an elastic matrix containing a rigid inclusion of the shape of a thin cylindrical shell. It is assumed that harmonic vibrations occur in the matrix under the conditions of axial symmetry (the symmetry axis is the inclusion axis) and the conditions of full adhesion between the inclusion and the matrix are satisfied. The vibrations are caused by the propagation of a plane wave whose front is perpendicular to the inclusion axis. The solution method is based on representing the displacements in the matrix as discontinuous solutions of the equations of axisymmetric oscillations of an elastic medium with unknown stress jumps on the inclusion surface. The realization of the boundary conditions for these jumps leads to a system of integral equations. Its solution is constructed numerically by the mechanical quadrature method with the use of special quadrature formulas for specific integrals. It is numerically investigated how the ratio of the inclusion geometric dimensions and the propagating wave frequency affect the stress concentration near the inclusion.     

Excitation by sound of Tollmien-Schlichting waves in a supersonic boundary layer

The interaction of sound with a supersonic boundary layer is considered. Because of the dependence of the main flow on the longitudinal coordinate, a sound wave generates unstable oscillations within the boundary layer. Calculations made for Mach number M = 2.0 and dimensionless frequency 2πfv_e/U_e ² = 0.91·10^?4 showed that near the lower branch of the curve of neutral stability a Tollmien—Schlichting wave can be excited with an intensity 2–3 times greater than that of the external acoustic wave.     

水下声波传播过程的时域间断Galerkin有限元方法

本文构建了声压波动方程的改进时域间断Galerkin有限元方法.传统时域连续有限元方法在计算高梯度、强间断特征水中声波传播问题时往往会出现虚假数值振荡现象,这些数值振荡会影响正常波动的计算精度.为了解决这一问题,本文通过引入人工阻尼的方式构建了改进的时域间断Galerkin有限元方法,并针对具有高梯度、强间断特征的多障...     

Intrinsic transfer matrix method and split quaternion formalism for multilayer media

The intrinsic transfer matrix method for 1D longitudinal elastic wave propagation through multilayer media is used to obtain an equivalent split quaternion formalism. Periodic media are analysed in this framework and the presence of a defect is considered. A simple one layer defect and an inversion defect are analysed. A commutative type of split quaternion is identified, which corresponds to defects of the periodic structure that can be placed in any position, the overall acoustic properties of the medium being conserved. Also, a nonperiodic medium composed of such commutative elements has a behaviour independent of the order of elements. Several possible applications in sound wave measurement and processing are outlined. The proposed split quaternion formalism is compact and can make analytical and numerical computations easier.     

Active transient elasto-acoustic response damping of a thick-walled liquid-coupled piezolaminated cylindrical vessel

The linear 3D piezoelasticity theory in conjunction with the versatile transfer matrix approach and the wave equation for the internal acoustic domain are employed for active non-stationary vibroacoustic response control of an arbitrarily thick, tri-laminate, fluid-filled, simply supported, piezocomposite cylindrical tank, excited by arbitrary (non-axisymmetric) time-dependent on-surface mechanical loads. The smart structure is composed of a supporting core layer of functionally graded orthotropic material perfectly bonded to inner and outer spatially distributed radially polarized functionally graded piezoceramic sensor and uniform force actuator (FGPM) layers. Active vibration damping is implemented by transferring the accumulated voltage on the sensor layer to the piezoelectric actuator layer in context of proportional and derivative control laws. Durbin's numerical inverse Laplace transform scheme is utilized to calculate the time response histories of the relevant interface displacement/stress components, center-point acoustic pressure, and actuator voltage, for selected loading configurations (i.e., concentrated step, impulse, and moving external loads). Numerical simulations demonstrate the effectiveness of the adopted distributed sensing/actuation configuration together with the active damping control strategy in suppressing the vibroacoustic response of a three-layered (Ba₂NaNb₅O₁₅/Al/PZT4) water-filled piezoelastic cylindrical tank. Limiting cases are considered and the validity of results is established by comparison with the available data as well as with the aid of a commercial finite element package.     

Fast Fourier homogenization for elastic wave propagation in complex media

In the context of acoustic or elastic wave propagation, the non-periodic asymptotic homogenization method allows one to determine a smooth effective medium and equations associated with the wave propagation in a given complex elastic or acoustic medium down to a given minimum wavelength. By smoothing all discontinuities and fine scales of the original medium, the homogenization technique considerably reduces meshing difficulties as well as the numerical cost associated with the wave equation solver, while producing the same waveform as for the original medium (up to the desired accuracy). Nevertheless, finding the effective medium requires one to solve the so-called “cell problem”, which corresponds to an elasto-static equation with a finite set of distinct loadings. For general elastic or acoustic media, the cell problem is a large problem that has to be solved on the whole domain and its resolution implies the use of a finite element solver and a mesh of the fine scale medium. Even if solving the cell problem is simpler than solving the wave equation in the original medium (because it is time and source independent, based on simple tetrahedral meshes and embarrassingly parallel) it is still a challenge. In this work, we present an alternative method to the finite element approach for solving the cell problem. It is based on a well-known method designed by H. Moulinec and P. Suquet in 1998 in structural mechanics. This iterative technique relies on Green functions of a simple reference medium and extensively uses Fast Fourier Transforms. It is easy to implement, very efficient and relies on a simple regular gridding of the medium. Through examples we show that the method gives excellent results, even, under some conditions, for discontinuous media.     

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 ε.     

基于扭曲Kagome点阵结构的一模材料设计

零能模式超材料指弹性矩阵的特征值中有若干为零的弹性材料,根据零特征值的个数可将其分类为一模至五模材料。当前,针对五模材料已有较深入研究,并在水声和弹性波调控方面获得重要应用,而对其他类型零能模式材料的研究尚未展开。本文对扭曲Kagome周期桁架这样一类欠约束点阵材料的有效弹性性质进行了研究,结果表明通过调节点阵材料的微观几何构型和杆件刚度,该类结构能够涵盖一系列一模材料谱系。针对给定一模弹性张量,发展了软-硬模式分离的微结构逆向优化设计策略。通过特定一模材料中的波传播现象对有效性质预测和微结构设计进行了数值验证。     

Extinction of elastic wave energy due to scattering in a viscoelastic medium

Forward scattering theorem for elastic longitudinal and shear wave scatterings by an arbitrary-shaped three-dimensional object embedded in a viscoelastic medium is derived. It is shown that the formulae for extinction cross-sections of an object in an energy-absorbing medium are formally the same with those of the object in the lossless elastic medium. Numerical calculations are executed for the longitudinal wave scattering in an epoxy matrix by a spherical inclusion with different material properties. The condition of negative extinction is examined with the causality constraint on the viscoelastic medium taken into account. It is found that the negative extinction occurs in the Rayleigh limit when the attenuation of the medium is sufficiently high and, more restrictedly, the wave speed in the object is larger than that in the medium, while it occurs less likely in the high frequency range considered in this paper (0<ka<100).     

Determination of elastic moduli of composite medium containing bimaterial matrix and non-uniform inclusion concentrations

Reservoir porous rocks usually consist of more than two types of matrix materials,forming a randomly heterogeneous material.The determination of the bulk modulus of such a medium is critical to the elastic wave dispersion and attenuation.The elastic moduli for a simple matrix-inclusion model are theoretically analyzed.Most of the efforts assume a uniform inclusion concentration throughout the whole single-material matrix.However,the assumption is too strict in real-world rocks.A model is developed to estimate the moduli of a heterogeneous bimaterial skeleton,i.e.,the host matrix and the patchy matrix.The elastic moduli,density,and permeability of the patchy matrix differ from those of the surrounding host matrix material.Both the matrices contain dispersed particle inclusions with different concentrations.By setting the elastic constant and density of the particles to be zero,a double-porosity medium is obtained.The bulk moduli for the whole system are derived with a multi-level effective modulus method based on Hashin's work.The proposed model improves the elastic modulus calculation of reservoir rocks,and is used to predict the kerogen content based on the wave velocity measured in laboratory.The results show pretty good consistency between the inversed total organic carbon and the measured total organic carbon for two sets of rock samples.     

Boundary element modeling for defect characterization potential in a wave guide

Wave scattering analysis implemented by boundary element methods (BEM) and the normal mode expansion technique is used to study the sizing potential of two-dimensional shaped defects in a wave guide. Surface breaking half-elliptical shaped defects of three opening lengths (0.3, 6.35 and 12.7 mm) and through-wall depths of 10–90% on a 10 mm thick steel plate were considered. The reflection and transmission coefficients of both Lamb and shear horizontal (SH) waves over a frequency range 0.05–2 MHz were studied. A powerfully practical result was obtained whereby the numerical results for the S₀ mode Lamb wave and n₀ mode SH wave at low frequencies showed a monotonic increase in signal amplitude with an increase in the defect through-wall depth. At high frequency (usually above the cut-off frequency of the A₁ mode for Lamb waves and the n₁ mode for SH waves, respectively), the monotonic trend does not hold in general due to the energy redistribution to the higher order wave modes. Guided waves impinging onto an internal stringer-like an inclusion were also studied. Both the Lamb and SH waves were shown to be insensitive to the stringer internal inclusions at low frequency. Experiments with piezoelectric Lamb wave transducers and non-contact SH wave electro-magnetic acoustic transducers (EMAT) verified some of the theoretical results.     

Lattice Boltzmann investigation for enhancing the thermal conductivity of ice using Al2O3 porous matrix

An enthalpy-based Lattice Boltzmann method (LBM) with double-distribution function (DDF) model is used to investigate numerically the effects of inserting a porous matrix on the heat transfer performance of the phase change material (PCM). Simulations are carried out for melting of ice in saturated Al₂O₃ porous matrix encapsulated in a concentric annulus. The process is considered as a conduction/convection controlled phase change problem at a representative elementary volume (REV) scale. The present results are validated by previous published numerical simulations of melting with and without porous media. In this research paper, the effects of decreasing the porosity on the temperature contours, flow patterns within the melt zone, complete melting time of the PCM and average Nusselt number are discussed qualitatively and quantitatively.