A novel unsplit perfectly matched layer for the second-order acoustic wave equation

When solving acoustic field equations by using numerical approximation technique, absorbing boundary conditions (ABCs) are widely used to truncate the simulation to a finite space. The perfectly matched layer (PML) technique has exhibited excellent absorbing efficiency as an ABC for the acoustic wave equation formulated as a first-order system. However, as the PML was originally designed for the first-order equation system, it cannot be applied to the second-order equation system directly. In this article, we aim to extend the unsplit PML to the second-order equation system. We developed an efficient unsplit implementation of PML for the second-order acoustic wave equation based on an auxiliary-differential-equation (ADE) scheme. The proposed method can benefit to the use of PML in simulations based on second-order equations. Compared with the existing PMLs, it has simpler implementation and requires less extra storage. Numerical results from finite-difference time-domain models are provided to illustrate the validity of the approach.     

Unsplit perfectly matched layer for acoustic field simulation based on second-order wave equation

The Perfectly Matched Layer(PML) is an effective absorbing boundary and has been widely used in acoustic simulation.In order to develop an absorbing boundary condition for numerical simulation based on the second-order acoustic wave equation,an Unsplit PML algorithm is proposed.Firstly,frequency-domain expression of this method is formulated based on the complex stretching coordinate schema.Then,its time-domain expression is derived by formulating auxiliary differential equations.Finally,relative theoretical analysis and numerical simulations are carried out,the results of which demonstrate that compared with the existing Split PMLs,the proposed method has the same absorbing efficiency and can reduce storage greatly.It can also increase calculation efficiency with easier implementation.     

有耗介质完全匹配层在时域多分辨分析中的应用

 对无耗和有耗介质中的完全匹配问题进行了研究,将广义完全匹配吸收层(GPML)应用到时域多分辨分析（MRTD）中。GPML是在扩展坐标系下由Maxwell方程得到的,在MRTD中实现对GPML的求解,并在频域对GPML进行了有效性分析,得出了GPML在不同空间步长和不同吸收层厚度情况下入射波由真空入射到吸收层的反射系数在频域的分布,并给出了有耗介质中GPML吸收效果的例子。结果表明,GPML的反射系数在一定频率范围内小于-40 dB,且随着计算空间步长的减小和吸收层厚度的增加而减小。GPML既可以用来截断无耗介质也可以用来截断有耗介质,为时域多分辨分析方法提供了一种较为通用的吸收边界。     

常规分裂和非分裂完全匹配层吸收边界比较研究

分别对常规分裂和非分裂完全匹配层吸收边界在完全弹性介质和孔隙介质的应用效果进行了比较研究。比较了两种吸收边界应用于完全弹性介质和孔隙介质所需计算量、存储量大小以及吸收效果。为比较吸收效果,首先给出两种吸收边界下声场快照图的比较,然后利用吸收边界附近的接收点波形,计算和分析了两种吸收边界在固定入射角度处随时间变化的边界反射幅度,以及对不同入射角度声波的边界反射系数。研究结果表明,无论是完全弹性介质还是孔隙介质,常规非分裂完全匹配层吸收边界比分裂完全匹配层吸收边界对垂直入射乃至一般入射角度的声波都有更好的吸收效果,但需要更大的存储量和计算量。     

常规分裂和非分裂完全匹配层吸收边界比较研究

分别对常规分裂和非分裂完全匹配层吸收边界在完全弹性介质和孔隙介质的应用效果进行了比较研究。比较了两种吸收边界应用于完全弹性介质和孔隙介质所需计算量、存储量大小以及吸收效果。为比较吸收效果,首先给出两种吸收边界下声场快照图的比较,然后利用吸收边界附近的接收点波形,计算和分析了两种吸收边界在固定入射角度处随时间变化的边界反射幅度,以及对不同入射角度声波的边界反射系数。研究结果表明,无论是完全弹性介质还是孔隙介质,常规非分裂完全匹配层吸收边界比分裂完全匹配层吸收边界对垂直入射乃至一般入射角度的声波都有更好的吸收效果,但需要更大的存储量和计算量。     

截断各向异性等离子体的修正各向异性完全匹配层吸收边界

 基于磁化等离子体本构方程,提出了一种截断各向异性色散介质的修正的各向异性完全匹配层(M-UPML)吸收边界算法。通过等效相对介电常数,将UPML推广到截断各向异性等离子体介质的情形,并推导了其时域有限差分方法(FDTD)迭代式。用该方法计算了半空间磁化与非磁化等离子体的反射系数,计算结果与解析解相一致,表明该吸收边界具有良好的吸收效果。     

Drift motion of domain boundaries in garnet ferrites in an acoustic wave field

The drift motion of a 180° domain boundary in garnet ferrites with two nonequivalent sublattices is studied in an elastic stress field induced by an acoustic wave propagating in the domain boundary plane. The dependences of the drift velocity on the amplitude and polarization of the acoustic wave are obtained, and the drift motion conditions for a strip domain structure are determined.     

A non-uniform basis order for the discontinuous Galerkin method of the 3D dissipative wave equation with perfectly matched layer

In this study a discontinuous Galerkin method (DG) for solving the three-dimensional time-dependent dissipative wave equation is investigated. In the case of unbounded problems, the perfectly matching layer (PML) is used to truncate the computational domain. The aim of this work is to investigate a simple selection method for choosing the basis order for elements in the computational mesh in order to obtain a predetermined error level. The selection method studied here relies on the error estimates provided by Ainsworth [M. Ainsworth, Dispersive and dissipative behaviour of high order discontinuous Galerkin finite element methods, Journal of Computational Physics 198(1) (2004) 106–130]. The performance of the non-uniform basis is examined using numerical experiments. In the simulated model problems, a feasible method choosing the basis order for arbitrary sized elements is achieved. In simulations, the effect of dissipation and the choices of the PML parameters on the performance of the DG method are also investigated.     

Frequency/time domain modeling of a direct drive point absorber wave energy converter

Two different methods to model a point absorber wave energy converter(WEC)with direct drive linear power take-off(PTO)are proposed in the present study:the frequency domain(FD)method and the time domain(TD)method.In the FD analysis,the frequency response function(FRF)of the WEC device is obtained via the equation of motion,and the expressions of power capture width in regular and random waves are derived as well.In the TD modeling,based on a state space approximation of the convolution term in the motion equation,both regular wave and random wave simulations are carried out.The regular wave simulation results indicate that the state space approximation is sufficiently accurate and the capture width reaches the maximum in the vicinity of the natural frequency.In the random wave simulations,the effects of buoy size,the PTO damping and wave climate on the power capture width are discussed in detail,which leads to the conclusion that the capture widths are influenced by the natural frequency of the WEC device,peak frequency of the wave spectrum,the amplitude of FRF and PTO damping.Furthermore,the increase of the capture width is at the cost of a relatively large buoy size and PTO damping when control is not included.     

Optical modelling of photonic crystals and VCSELs using eigenmode expansion and perfectly matched layers

We present a modelling approach for photonic crystal structures and vertical-cavity surface-emitting lasers (VCSELs). This method is based on vectorial eigenmode expansion combined with perfectly matched layer (PML) boundary conditions. Compared to other methods, a relatively small computational effort is required, while at the same time accurate results are obtained, even in the presence of strong scattering and diffraction losses.     

Fast analysis of transient acoustic wave scattering from rigid bodies using the multilevel plane wave time domain algorithm

The analysis of transient wave scattering from rigid bodies using integral equation-based techniques is computationally intensive: if carried out using classical schemes, the evaluation of the velocity potential on the surface of a three-dimensional scatterer, represented in terms of Ns spatial basis functions for Nt time steps, requires O(NtNs2) operations. The recently developed plane wave time domain (PWTD) algorithm permits the rapid evaluation of transient fields that are generated by bandlimited source distributions. It has been shown that incorporation of the PWTD algorithm into integral equation-based solvers in a two-level setting reduces the computational complexity of a transient analysis to O(NtNs1.5 logNs). In this paper, it is shown that casting the PWTD scheme into a multilevel framework permits the analysis of transient acoustic surface scattering phenomena in O(NtNslog2Ns) operations using O(NtNs) memory. Numerical examples that demonstrate the efficacy of the multilevel implementation are also presented.     

Rigorous and approximate methods for modeling wave scattering from a locally perturbed perfectly conducting surface

Rigorous and approximate methods are considered for solving the problem of harmonic plane wave scattering from a plane surface arbitrarily perturbed along one dimension on a finite interval. This problem is treated using the Fredholm integral equations of the second kind and the Kirchhoff and Rayleigh approximations. The estimates of the computational efficiency of the integral equation method and the Rayleigh approximation are compared by calculating fields scattered from random rough surfaces in the resonance region (i.e., when the roughness height is comparable to or smaller than the incident wavelength) for an arbitrary incidence of a plane wave. Scattering patterns calculated using the integral equations and the Kirchhoff approximation are discussed in the case of large-scale random rough surface scattering. Particular attention is paid to scattering at near-grazing incidence.     

Modeling and simulation of experimentally fabricated QDSSC using ZnS as light absorbing and blocking layer

Abstract—Two main factors which limit the power conversion efficiency of solar cells are light absorption and recombination processes. In photovoltaic (PV) devices, low energy photons cannot be absorbed and excite electrons from valance band to conduction band, hence do not contribute to the current. On the other hand, high energy photons cannot be efficiently used due to a poor match to the energy gap. Existence of charge recombination in PV devices causes the low conversion performance, which is indicated by the low open-circuit voltage (V _OC ). Using a blocking layer in system could effectively reduce the recombination of charge carriers. In this study, we simulated a solar cell with ITO/ZnO/P3HT&PCBM/Ag structure. To prevent the charge recombination, a ZnS QD layer was used which acts as a light absorbing and a recombination blocking layer in the ITO/ZnO film/ZnS QD/P3HT&PCBM/Ag structure. The simulated J–V characteristics of solar cells showed a close match with the experimental results. Simulate data showed an increase of conversion efficiency in ZnS QDSSC from 1.71 to 3.10%, which is relatively 81.28% increase.     

High-performance modeling acoustic and elastic waves using the parallel Dichotomy Algorithm

A high-performance parallel algorithm is proposed for modeling the propagation of acoustic and elastic waves in inhomogeneous media. An initial boundary-value problem is replaced by a series of boundary-value problems for a constant elliptic operator and different right-hand sides via the integral Laguerre transform. It is proposed to solve difference equations by the conjugate gradient method for acoustic equations and by the GMRES(k) method for modeling elastic waves. A preconditioning operator was the Laplace operator that is inverted using the variable separation method. The novelty of the proposed algorithm is using the Dichotomy Algorithm [26], which was designed for solving a series of tridiagonal systems of linear equations, in the context of the preconditioning operator inversion. Via considering analytical solutions, it is shown that modeling wave processes for long instants of time requires high-resolution meshes. The proposed parallel fine-mesh algorithm enabled to solve real application seismic problems in acceptable time and with high accuracy. By solving model problems, it is demonstrated that the considered parallel algorithm possesses high performance and efficiency over a wide range of the number of processors (from 2 to 8192).     

Evaluation of the acoustic and non-acoustic properties of sound absorbing materials using a three-microphone impedance tube

This paper presents a straightforward application of an indirect method based on a three-microphone impedance tube setup to determine the non-acoustic properties of a sound absorbing porous material. First, a three-microphone impedance tube technique is used to measure some acoustic properties of the material (i.e., sound absorption coefficient, sound transmission loss, effective density and effective bulk modulus) regarded here as an equivalent fluid. Second, an indirect characterization allows one to extract its non-acoustic properties (i.e., static airflow resistivity, tortuosity, viscous and thermal characteristic lengths) from the measured effective properties and the material open porosity. The procedure is applied to four different sound absorbing materials and results of the characterization are compared with existing direct and inverse methods. Predictions of the acoustic behavior using an equivalent fluid model and the found non-acoustic properties are in good agreement with impedance tube measurements.     

A comparison of statistically optimized near field acoustic holography using single layer pressure-velocity measurements and using double layer pressure measurements

Statistically optimized near field acoustic holography (SONAH) is usually based on the assumption that all sources are on one side of the measurement plane whereas the other side is source free. An extension of the SONAH procedure based on measurement with an array of pressure-velocity probes has recently been suggested. An alternative method uses a double layer array of pressure transducers. Both methods make it possible to distinguish between sources on the two sides of the array and thus suppress the influence of extraneous noise and reflections coming from the "wrong" side. This letter compares the two methods.     

Forward and backward wave propagation in multilayer planar waveguide using metamaterials layer

In this study, an efficient method is discussed to analyze the multilayer planar waveguides with double negative guided (DNG) and double positive material as guiding film. Here, among various lossless multilayer planar waveguide structures, only three layer and five layer structures are discussed. For these structures, guided dispersion characteristic, along with electric field distribution of TM modes are numerically analyzed and compared. This analysis enables an effective comparison of guided modal properties of various modes and results in a better understanding of the multilayer planar waveguide with DNG as guiding film.     

Measurement of material nonlinearity using surface acoustic wave parametric interaction and laser ultrasonics

A dual frequency mixing technique has been developed for measuring velocity changes caused by material nonlinearity. The technique is based on the parametric interaction between two surface acoustic waves (SAWs): The low frequency pump SAW generated by a transducer and the high frequency probe SAW generated and detected using laser ultrasonics. The pump SAW stresses the material under the probe SAW. The stress (typically <5 MPa) is controlled by varying the timing between the pump and probe waves. The nonlinear interaction is measured as a phase modulation of the probe SAW and equated to a velocity change. The velocity-stress relationship is used as a measure of material nonlinearity. Experiments were conducted to observe the pump-probe interaction by changing the pump frequency and compare the nonlinear response of aluminum and fused silica. Experiments showed these two materials had opposite nonlinear responses, consistent with previously published data. The technique could be applied to life-time predictions of engineered components by measuring changes in nonlinear response caused by fatigue.     

Development of an ultrasonic inspection robot using an electromagnetic acoustic transducer for a Lamb wave and an SH-plate wave

For inspection of a storage tank and pipeline in service, the application of an automatic inspection system (nondestructive inspection robot) is desirable, because manual inspection is difficult to perfectly and exactly perform due to the enormous amount of inspection needed. However, an ultrasonic nondestructive inspection robot with a piezoelectric oscillator needs to touch only the material surface to be directly inspected using a coupling medium. That is, the material surface and the sensor must always be held by constant pressure in the vertical direction on the material side. Actually, it is difficult to overcome these problems; thus an ultrasonic inspection robot could not be widely applied. We then tried to develop an ultrasonic inspection robot with an electromagnetic acoustic transducer (EMAT) which did not require a coupling medium to inspect the circumferential pipe parts. We developed a special EMAT that could transmit and receive alternately a Lamb wave with high sensitivity and a SH-plate wave without influence by the welded part. The method by which the inspection robot turned around the direction of the steel pipe surroundings was executed by observing the tape pasted in the direction of the steel pipe surroundings with an installed CCD camera. In this report, the basic mechanism of this inspection robot and an examination of results are described.     

Design and modeling of inversion layer ultrasonic transducers using LiNbO3 single crystal

Using inversion domain engineering controlled by heating temperature, the LiNbO(3) (LNO) piezoelectric plate with both odd and even-order thickness-extensional modes can be excited simultaneously. Therefore, the inversion layer ultrasound transducer is expected to be capable of operating over a wider frequency range. In this paper, the electrical impedance and the acoustic characteristics of LiNbO(3) (LNO) inversion layer transducer have been studied by finite element modeling (FEM). The transducer designed for this study uses a 36 degrees rotated Y-cut LiNbO(3) thin plate with an active element thickness of approximately 100 microm. First the electrical and elastic properties of the 36 degrees rotated Y-cut LNO were obtained by transforming a basic piezoelectric matrix for Z-cut LNO. In order to validate the FEM using the transformed properties several pieces of pure and 50% inversion layer LNO were tested on the electrical impedance analyzer. The modeled impedance characteristics were consistent with the measured data. Next the model was used to design 50-60 MHz transducers using pure and 30% inversion LNO. Two lambda/4 matching layers and a Tungsten loaded epoxy backing were used in these designs. The modeled results show that an over 90% bandwidth transducer can be made with proper matching and 30% inversion layer.