Spectral and perturbation analysis for ultrasonic guided waves

The influence of damage on waves propagating in complex geometry waveguides is investigated through a numerical model formulated by combining the Spectral Finite Element Method and Perturbation Techniques. The resulting numerical tool allows efficient computation of the wave propagation response and the analysis of the effects of damages of various extent and location. The dynamic behavior of the damaged waveguides is described through a general higher order model which couples different waves thus allowing the prediction of mode conversion phenomena. Arbitrary cross-section can be considered through Finite Element (FE) discretization according to well-established Semi-Analytical Finite Element (SAFE) procedures. Two types of damages which allow the application of perturbation theory are considered: a small localized reduction of the thickness and a reduction of material stiffness and density. A validation by comparison with a Finite Element Model as well as numerical examples are presented to illustrate the model capabilities.     

Computation of leaky guided waves dispersion spectrum using vibroacoustic analyses and the Matrix Pencil Method: A validation study for immersed rectangular waveguides

The paper aims at validating a recently proposed Semi Analytical Finite Element (SAFE) formulation coupled with a 2.5D Boundary Element Method (2.5D BEM) for the extraction of dispersion data in immersed waveguides of generic cross-section. To this end, three-dimensional vibroacoustic analyses are carried out on two waveguides of square and rectangular cross-section immersed in water using the commercial Finite Element software Abaqus/Explicit. Real wavenumber and attenuation dispersive data are extracted by means of a modified Matrix Pencil Method. It is demonstrated that the results obtained using the two techniques are in very good agreement.     

Elastic wave propagation in sinusoidally corrugated waveguides

The ultrasonic wave propagation in sinusoidally corrugated waveguides is studied in this paper. Periodically corrugated waveguides are gaining popularity in the field of vibration control and for designing structures with desired acoustic band gaps. Currently only numerical method (Boundary Element Method or Finite Element Method) based packages (e.g., PZFlex) are in principle capable of modeling ultrasonic fields in complex structures with rapid change of curvatures at the interfaces and boundaries but no analyses have been reported. However, the packages are very CPU intensive; it requires a huge amount of computation memory and time for its execution. In this paper a new semi-analytical technique called Distributed Point Source Method (DPSM) is used to model the ultrasonic field in sinusoidally corrugated waveguides immersed in water where the interface curvature changes rapidly. DPSM results are compared with analytical solutions. It is found that when a narrow ultrasonic beam hits the corrugation peaks at an angle, the wave propagates in the backward direction in waveguides with high corrugation depth. However, in waveguides with small corrugation the wave propagates in the forward direction. The forward and backward propagation phenomenon is found to be independent of the signal frequency and depends on the degree of corrugation.     

多量子阱平面光波导的有限元分析

本文利用有限元法分析了多量子阱平面光波导的传播特性,给出了任意阱数多量子阱光波导_TE模和TM模的有效折射率和TE模的强度分布.结果表明:γ(势阱与势垒的厚度比)值的改变,影响波导中导模的模式数目、有效折射率、双折射和强度分布;并发现在某些情况下,均方根近似是不适用的.     

An Extension of the Effective Index Method to Analyze Leakage Losses in Rib Type Waveguides

An extension of the effective index method is proposed as a tool to investigate leakage losses in two-dimensional waveguides. The above method, referred to as the Extended Effective Index Method (EEIM), utilizes the Transfer Matrix Technique (TMT) for the calculation of the complex propagation constants. The results show that the EEIM can be applied to conventional rib leaky waveguides as well as to rib ARROW leaky waveguides. The results for rib ARROW waveguides show excellent agreement when compared to those obtained with Finite Element Method.     

一种水平变化波导中声传播问题的耦合模态法

针对介质参数及海底边界水平变化波导中的声传播问题,本文基于多模态导纳法提出一种能量守恒且便于数值稳定求解的耦合模态方法.将声压表示为一组正交完备的本地本征函数之和,对声压满足的Helmholtz方程在本地本征函数上作投影,推导出关于声压模态系数的二阶耦合模态方程组.耦合矩阵直观描述水平变化因素对模态耦合的贡献.为避免直接求解二阶耦合模态方程组可能遇到的数值发散问题,将其重构为两个耦合的一阶演化方程组,引入导纳矩阵并使用Magnus数值积分方法获得稳定的声场解.利用该耦合模态方法数值计算水平变化波导中的声场,并与COMSOL参考解比较,结果表明该耦合模态理论能够精确求解水平变化波导中的点源及分布源传播问题.     

基于模式分离的兰姆波逆时偏移成像

应用频率域逆时偏移方法实现各向同性和各向异性板中缺陷的兰姆波成像.由于缺陷引起的多模态散射信号会在重建图像中形成伪像,根据基本导波模式振动对称性的差别进行了模式分离预处理.基于多元阵列超声技术,开展了铝板和复合板内缺陷频率域逆时偏移超声成像方法的数值仿真研究.首先,建立有限元模型,采用环形传感器数值采集由缺陷引起的兰姆波散射信号,然后,将采集到的多模式散射信号进行模式分离处理,再将模式分离后的兰姆波散射信号经时间反转后并在相应的接收器处重新激励,在频域中运用格林函数反向传播兰姆波散射信号,获取监测区域的声场信息,与正向传播声场进行互相关,重建缺陷图像.首先对铝板中单缺陷以及复合材料板中相邻的两个相同缺陷进行数值仿真,对比有无模式分离处理的缺陷逆时偏移成像效果,体现出模式分离的重要性.在此基础上,采用逆时偏移方法对复合板材内位置邻近、深度不同的双缺陷进行识别.数值结果表明,模式分离预处理后的缺陷重建图像能够有效去除多模式干扰产生的伪像.文中提出的成像方法对各向同性板和各向异性板内缺陷的检测和成像具有很好的发展潜力,可以准确地探测多个缺陷的形状、尺寸和深度.     

室内粗糙表面声散射体掠入射角度散射估计

针对经典边界元方法对掠入射角度下的散射估计精度较低的问题,发展了一种边界无网格模型.此模型将散射体视为具有无限延伸的形式,避免了经典边界元法中薄板形式所导致的声压差问题,而且更加符合实际房间中散射体的存在形式;模型利用无网格算法实现数值仿真,可对任意表面形状特别是曲面散射体具有更高的仿真精度。利用边界无网格模型计算了不同形状散射体的散射系数及散射声场,并将结果与解析方法、测量实验进行了对比.对比结果表明,边界无网格模型可以准确预测散射体的散射性质,特别是对掠入射角度的估计要优于经典边界元法.研究结论可应用于室内声学散射体特征预测及优化设计,对提高声场扩散及室内音质水平具有重要意义.      

Analysis and design of efficient coupling in photonic crystal circuits

A rigorous analysis and design of efficient coupling from photonic crystal (PhC) waveguides into conventional dielectric waveguides is reported. Closed-form expressions for the reflection and transmission matrices that completely characterize the scattering that occurs at the interface are derived based on an eigenmode expansion technique and a Bloch basis. Analytic expressions are used to analyze the reflection into PhC waveguides. We obtain that negligible reflection can be achieved by choosing a certain interface within a PhC unit cell. Furthermore, analytic expressions are used to design a novel and compact coupler structure in order to achieve high coupling efficiency when broad dielectric waveguides are considered. Thereby, transmission efficiencies near 100 from the fundamental guided Bloch mode into the fundamental waveguide mode are achieved.     

Computation of the scattering matrix of guided acoustical propagation by the Wave Finite Element approach

This work aims to study the guided acoustical propagation. The Wave Finite Element method is applied rather than the classic Finite Element method. Finite Element libraries are used for an elementary portion modelling. The dynamic stiffness matrix is then calculated. Periodicity conditions lead to a simple eigenvalue problem and the wave basis can be extracted. The dynamic stiffness matrix of the coupling elements corresponding to the lined parts can be also calculated considering the acoustical impedance. The use of the coupling’s conditions provides the scattering matrix. Within the framework of this method, the dispersion curves, the evolution of the scattering coefficients and the forced response to pressure excitations for both single and coupled waveguides can be represented, and compared to the FE method. The WFE method, treated in literature, and its main interests are reviewed, while adding the impedance notion into the calculation of the scattering matrix.     

同轴波导内金属支撑杆的散射特性研究

 对同轴波导内传输TEM模时金属支撑杆的散射特性进行了研究。通过将同轴波导等效为窄边是磁壁的矩形波导,建立了金属支撑杆的散射场计算模型,用积分方程法得到了散射场的计算公式,其计算结果与有限元方法仿真结果吻合良好。分析计算和实验表明：当同轴波导传输TEM模时,n根角向均匀分布的支撑杆在同轴波导内激励起一系列TE(mn)1模式;当n大于同轴波导截面平均周长与波长之比时,金属杆所激励的高阶模截止,输出和反射模式仅为TEM模,其反射系数随支撑杆根数或支撑杆半径的增加而增大,随频率的升高而减小。     

Full-vector mode solver for bending waveguides based on the finite-difference frequency-domain method in cylindrical coordinate systems

A full-vector mode solver for bending waveguides is described based on the finite-difference frequency-domain method in a local cylindrical coordinate system in which the formulas are directly derived from Maxwell's equations with the help of Yee's mesh. The perfectly matched layer absorbing boundary conditions are incorporated into the present approach in order to effectively demonstrate the leaky nature of bending waveguides. A typical bending rib waveguide is considered as a numerical example to show the effectiveness of the established method.     

Scattering from extended targets in range-dependent fluctuating ocean-waveguides with clutter from theory and experiments

Bistatic, long-range measurements of acoustic scattered returns from vertically extended, air-filled tubular targets were made during three distinct field experiments in fluctuating continental shelf waveguides. It is shown that Sonar Equation estimates of mean target-scattered intensity lead to large errors, differing by an order of magnitude from both the measurements and waveguide scattering theory. The use of the Ingenito scattering model is also shown to lead to significant errors in estimating mean target-scattered intensity in the field experiments because they were conducted in range-dependent ocean environments with large variations in sound speed structure over the depth of the targets, scenarios that violate basic assumptions of the Ingenito model. Green's theorem based full-field modeling that describes scattering from vertically extended tubular targets in range-dependent ocean waveguides by taking into account nonuniform sound speed structure over the target's depth extent is shown to accurately describe the statistics of the targets' scattered field in all three field experiments. Returns from the man-made targets are also shown to have a very different spectral dependence from the natural target-like clutter of the dominant fish schools observed, suggesting that judicious multi-frequency sensing may often provide a useful means of distinguishing fish from man-made targets.     

Numerical investigation of dispersion relations for helical waveguides using the Scaled Boundary Finite Element method

In this paper, the dispersion properties of elastic waves in helical waveguides are investigated. The formulation is based on the Scaled Boundary Finite Element method (SBEFM). With a set of orthogonal unit basis introduced as the contravariant basis, the helical coordinate is firstly considered, where components of tensor retain the dimension of original quantity. Based on the strain–displacement relation, the eigenvalue matrix is obtained about wavenumbers and frequencies. The cross section of the waveguides is discretized by using high-order spectral elements. Moreover, the formulated linear matrix is utilized to design efficient and accurate algorithms to compute the eigenvalues of helical waveguides. Compared with the Pochhammer–Chree curves, the convergence and accuracy of the SBFEM are discussed. Finally, we give some dispersion curves for a wide range of lay angles and analyze in detail properties of cut-off frequency, mode separation and mode transition for elastic wave propagation in the helical waveguides.     

Seismic scattering from a spherical inclusion eccentrically located within a homogeneous, spherical host: theoretical derivation

The scattered fields from a spherical body eccentrically located within an otherwise homogeneous host sphere are derived by satisfying the boundary conditions at both interfaces simultaneously. The source, which may be composed of any linear combination of S and P waves, is also located arbitrarily within the host sphere. The scattering system may have applications in seismic scattering, by cavities or dense bodies located near the Earth surface or to the Slichter mode.     

The scattering of guided waves in partly embedded cylindrical structures

The scattering of ultrasonic guided waves at a point where a free cylindrical waveguide enters an embedding material is investigated. A modal solution that is valid when the guided waves are incident from the free section of the waveguide is developed. It is shown that in this case it is valid to consider only the modal fields over the cross section of the waveguide, neglecting the fields in the embedding material. As an application, the scattering of the lowest-order longitudinal mode in a cylindrical waveguide, L(0,1), is examined in detail. As well as considering epoxy resin as an embedding material, the case where the embedding material is replaced by a perfectly rigid boundary is discussed. The latter gives some insight into the role of nonpropagating and inhomogeneous waveguide modes in the scattering process. The results from the modal solution are validated using Finite Element modeling, very good agreement being obtained.     

Modal analysis of the range evolution of broadband wavefields in the North Pacific Ocean: low mode numbers

The results of mode-processing measurements of broadband acoustic wavefields made in the fall of 2004 as part of the Long-Range Ocean Acoustic Propagation Experiment (LOAPEX) in the eastern North Pacific Ocean are reported here. Transient wavefields in the 50-90 Hz band that were recorded on a 1400-m long 40 element vertical array centered near the sound channel axis are analyzed. This array was designed to resolve low-order modes. The wavefields were excited by a ship-suspended source at seven ranges, between approximately 50 and 3200 km, from the receiving array. The range evolution of broadband modal arrival patterns corresponding to fixed mode numbers ("modal group arrivals") is analyzed with an emphasis on the second (variance) and third (skewness) moments. A theory of modal group time spreads is described, emphasizing complexities associated with energy scattering among low-order modes. The temporal structure of measured modal group arrivals is compared to theoretical predictions and numerical simulations. Theory, simulations, and observations generally agree. In cases where disagreement is observed, the reasons for the disagreement are discussed in terms of the underlying physical processes and data limitations.     

The use of an orthogonality relation for reducing the size of finite element models for 3D guided waves scattering problems

The scattering of guided waves by complex shaped defects in three-dimensional (3D) waveguides is considered. For such problems, analytical solutions do not exist, and modal decomposition techniques based on the establishment of the displacement and stress fields in the vicinity of the scatterer are quite heavy and complicated to perform. On the other hand, finite elements (FE)-based methods constitute a powerful way to obtain solutions, but they are known to be very memory consuming. This paper proposes a post-processing technique, based on a 3D orthogonality relation, to decompose a complex acoustic field produced by a scatterer and predicted by a 3D FE model, into plane waves, the amplitudes of which are quantified in the far field. This technique allows important reductions in the size of the FE models to be made. Two applications are presented to demonstrate the potential of this method. The first one concerns the scattering of the S₀ Lamb wave incident on a flat bottom circular hole. In this example, the amplitude of each mode is calculated via the orthogonality relation-based method, and compared to that obtained by simply monitoring the displacements at appropriate through-thickness positions. In the second application, the incident S₀ Lamb mode is converted into five modes scattered by a defect of complex geometry.     

Waveguide mode solver based on Neumann-to-Dirichlet operators and boundary integral equations

For optical waveguides with high index-contrast and sharp corners, existing full-vectorial mode solvers including those based on boundary integral equations typically have only second or third order of accuracy. In this paper, a new full-vectorial waveguide mode solver is developed based on a new formulation of boundary integral equations and the so-called Neumann-to-Dirichlet operators for sub-domains of constant refractive index. The method uses the normal derivatives of the two transverse magnetic field components as the basic unknown functions, and it offers higher order of accuracy where the order depends on a parameter used in a graded mesh for handling the corners. The method relies on a standard Nyström method for discretizing integral operators and it does not require analytic properties of the electromagnetic field (which are singular) at the corners.     

Interference effect in scattering loss of high-index-contrast planar waveguides caused by boundary reflections

We present theoretical and experimental results on an interference effect caused by boundary reflections on the optical scattering loss in high-index-contrast planar waveguides. Analytical expressions for the polarization-dependent scattering loss are derived using a surface Green's function. For high-index-contrast waveguides of submicrometer dimensions a significant deviation from accepted theory arises, including scattering loss suppression owing to a thin-film interference effect. Our theoretical predictions are confirmed by loss measurements on silicon-on-insulator channel waveguides.