An improved convolution perfectly matched layer for elastic second-order wave equation

A convolution perfectly matched layer(CPML)can efficiently absorb boundary reflection in numerical simulation.However,the CPML is suitable for the first-order elastic wave equation and is difficult to apply directly to the second-order elastic wave equation.In view of this,based on the first-order CPML absorbing boundary condition,we propose a new CPML(NCPML)boundary which can be directly applied to the second-order wave equation.We first systematically extend the first-order CPML technique into second-order wave equations,neglecting the space-varying characteristics of the partial damping coefficient in the complex-frequency domain,avoiding the generation of convolution in the time domain.We then transform the technique back to the time domain through the inverse Fourier transform.Numerical simulation indicates that the space-varying characteristics of the attenuation factor have little influence on the absorption effect and increase the memory at the same time.A number of numerical examples show that the NCPML proposed in this study is effective in simulating elastic wave propagation,and this algorithm is more efficient and requires less memory allocation than the conventional PML absorbing boundary.

Borehole-GPR numerical simulation of full wave field based on convolutional perfect matched layer boundary

The absorbing boundary is the key in numerical simulation of borehole radar.Perfect match layer(PML) was chosen as the absorbing boundary in numerical simulation of GPR.But CPML(convolutional perfect match layer) approach that we have chosen has the advantage of being media independent.Beginning with the Maxwell equations in a two-dimensional structure,numerical formulas of finite-difference time-domain(FDTD) method with CPML boundary condition for transverse electric(TE) or transverse magnetic(TM) wave are presented in details.Also,there are three models for borehole-GPR simulation.By analyzing the simulation results,the features of targets in GPR are obtained,which can provide a better interpretation of real radar data.The results show that CPML is well suited for the simulation of borehole-GPR.

Fast division-free parallel structure for convolution perfectly matched layer in finite difference time domain method

Parallel acceleration of convolution perfectly matched layer （CPML） algorithm suffers from massive division operation which is widely accepted as one of the most expensive operations for the equipment such as graphic processing unit （GPU）, field programmable gate array （FPGA） etc. In pursuit of higher efficiency and lower power consumption, this article revisited the CPML theory and proposed a new fast division-free parallel CPML structure. By optimally rearranging the CPML inner iteration process, all the division operators can be eliminated and replaced by recalculating the related field updating coefficients offline. Experiments show that the proposed division-free structure can save more than 50% arithmetic instructions and 25% execution time of the traditional parallel CPML structure without any accuracy loss.

Application of CPML to truncate the open boundaries of cylindrical waveguides in 2.5-dimensional problems

In order to solve the problem of truncating the open boundaries of cylindrical waveguides used in the simulation of high power microwave （HPM） sources, this paper studies the convolutional PML （CPML） in the cylindrical coordinate system. The electromagnetic field＇s FDTD equations and the expressions of axis boundary conditions are presented. Numerical experiments are conducted to validate the equations and axis boundary conditions. The performance of CPML is simulated when it is used to truncate the cylindrical waveguides excited by the sources with different frequencies and modes in the 2.5-dimensional problems. Numerical results show that the maximum relative errors are all less than -90 dB. The CPML method is introduced in the 2.5-dimensional electromagnetic PIC software, and the relativistic backward wave oscillator is simulated by using this method. The results show that the property of CPML is much better than that of the Mur-type absorbing boundary condition when they are used to truncate the open boundaries of waveguides. The CPML is especially suitable for truncating the open boundaries of the dispersive waveguide devices in the simulation of HPM sources.