The key problem of finite-difference time-domain (FD-TD) method is the skillful application of special conditions on the boundaries of the computational domain. A new technique named Perfectly Matched Layer(PML) yield...The key problem of finite-difference time-domain (FD-TD) method is the skillful application of special conditions on the boundaries of the computational domain. A new technique named Perfectly Matched Layer(PML) yields a robust Absorbing Boundary Condition(ABC) independent of the angle of incidence and the frequency of outgoing waves. In this paper, the principle of the PML technique is briefly presented. Then some problems in the application and their settlements are discussed emphatically. Finally three numerical tests and a measured result are devoted to examine the accuracy and effectiveness of this approach.展开更多
The theory of perfectly matched layer (PML) artificial boundary condition (ABC), which is characterized by absorption any wave motions with arbitrary frequency and arbitrarily incident angle, is introduced. The co...The theory of perfectly matched layer (PML) artificial boundary condition (ABC), which is characterized by absorption any wave motions with arbitrary frequency and arbitrarily incident angle, is introduced. The construction process of PML boundary based on elastodynamic partial differential equation (PDE) system is developed. Combining with velocity-stress hybrid finite element formulation, the applicability of PML boundary is investigated and the numerical reflection of PML boundary is estimated. The reflectivity of PML and multi-transmitting formula (MTF) boundary is then compared based on body wave and surface wave simulations. The results show that although PML boundary yields some reflection, its absorption performance is superior to MTF boundary in the numerical simulations of near-fault wave propagation, especially in comer and large angle grazing incidence situations. The PML boundary does not arise any unstable phenomenon and the stability of PML boundary is better than MTF boundary in hybrid finite element method. For a specified problem and analysis tolerance, the computational efficiency of PML boundary is only a little lower than MTF boundary.展开更多
The perfectly matched layer(PML) is an effective technique for truncating unbounded domains with minimal spurious reflections. A fluid parabolic equation(PE) model applying PML technique was previously used to analyze...The perfectly matched layer(PML) is an effective technique for truncating unbounded domains with minimal spurious reflections. A fluid parabolic equation(PE) model applying PML technique was previously used to analyze the sound propagation problem in a range-dependent waveguide(Lu and Zhu, 2007). However, Lu and Zhu only considered a standard fluid PE to demonstrate the capability of the PML and did not take improved one-way models into consideration. They applied a [1/1] Padé approximant to the parabolic equation. The higher-order PEs are more accurate than standard ones when a very large angle propagation is considered. As for range-dependent problems, the techniques to handle the vertical interface between adjacent regions are mainly energy conserving and single-scattering. In this paper, the PML technique is generalized to the higher order elastic PE, as is to the higher order fluid PE. The correction of energy conserving is used in range-dependent waveguides. Simulation is made in both acoustic cases and seismo-acoustic cases. Range-independent and range-dependent waveguides are both adopted to test the accuracy and efficiency of this method. The numerical results illustrate that a PML is much more effective than an artificial absorbing layer(ABL) both in acoustic and seismo-acoustic sound propagation modeling.展开更多
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 t...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.展开更多
A good artificial boundary treatment in a seismic wave grid-based numerical simulation can reduce the size of the computational region and increase the computational efficiency,which is becoming increasingly important...A good artificial boundary treatment in a seismic wave grid-based numerical simulation can reduce the size of the computational region and increase the computational efficiency,which is becoming increasingly important for seismic migration and waveform inversion tasks requiring hundreds or thousands of simulations.Two artificial boundary techniques are commonly used:perfectly matched layers(PMLs),which exhibit the excellent absorption performance but impose a greater computational burden by using finite layers to gradually reduce wave amplitudes;and absorbing boundary conditions(ABCs),which have the high computational efficiency but are less effective in absorption because they employ the one-way wave equation at the exterior boundary.Naturally,PMLs have been combined with ABCs to reduce the number of PMLs,thus improving the computational efficiency;many studies have proposed such hybrid PMLs.Depending on the equations from which the ABCs are derived,there are two hybrid PML variants:the PML+unstretched ABC(UABC),in which the ABC is derived from a physical equation;or the PML+stretched ABC(SABC),in which the ABC is derived from the PML equation.Even though all the previous studies concluded that hybrid PMLs can improve the absorption performance,none of them quantified how many PMLs can be removed by combining the PML with the ABC compared with the pure PML.In this paper,we systematically study the absorption performance of the two hybrid PML variants.We develop a method to distinguish the artificial reflections from the PML-interior interface and those caused by the PML exterior boundary to accurately approximate the additional absorption achieved by using the UABC and the SABC.The reflection coefficients based on a theoretical derivation and numerical tests both show that the UABC amplifies most reflections and is not recommended in any situation;conversely,the SABC can always diminish reflections,but the additional absorption achieved by the SABC is relatively poor and cannot effectively reduce the number of PMLs.In contrast,we find that simply increasing the damping parameter improves absorption better than the PML+SABC.Our results show that the improvement in absorption achieved by combining the PML with either the SABC or the UABC is not better than that obtained by simply adjusting the damping profile of the PML;thus,combining the PML with the ABC is not recommended in practice.展开更多
An FDTD system associated with uniaxial perfectly matched layer(UPML) for an electromagnetic scattering problem in two-dimensional space in polar coordinates is considered.Particularly the FDTD system of an initial-...An FDTD system associated with uniaxial perfectly matched layer(UPML) for an electromagnetic scattering problem in two-dimensional space in polar coordinates is considered.Particularly the FDTD system of an initial-boundary value problems of the transverse magnetic(TM) mode to Maxwell's equations is obtained by Yee's algorithm,and the open domain of the scattering problem is truncated by a circle with a UPML.Besides,an artificial boundary condition is imposed on the outer boundary of the UPML.Afterwards,stability of the FDTD system on the truncated domain is established through energy estimates by the Gronwall inequality.Numerical experiments are designed to approve the theoretical analysis.展开更多
In this paper,we study the two-dimensional Helmholtz scattering problem by a locally perturbed line with impedance boundary condition.Different from the problem with Dirichlet boundary condition,the Green function of ...In this paper,we study the two-dimensional Helmholtz scattering problem by a locally perturbed line with impedance boundary condition.Different from the problem with Dirichlet boundary condition,the Green function of the Helmholtz equation with impedance boundary condition becomes very complicated and comprises surface waves along the locally perturbed line.A uniaxial perfectly matched layer(UPML)method is proposed to truncate the half plane into a bounded computational domain.The main contribution of this paper is to prove the well-posedness of the PML problem and the exponential convergence of the approximate solution to the exact solution as either the thickness or the medium parameter of PML increases.展开更多
We review time-domain formulations of radiation boundary conditions for Maxwell's equations, focusing on methods which can deliver arbitrary accuracy at acceptable compu-tational cost. Examples include fast evaluatio...We review time-domain formulations of radiation boundary conditions for Maxwell's equations, focusing on methods which can deliver arbitrary accuracy at acceptable compu-tational cost. Examples include fast evaluations of nonlocal conditions on symmetric and general boundaries, methods based on identifying and evaluating equivalent sources, and local approximations such as the perfectly matched layer and sequences of local boundary conditions. Complexity estimates are derived to assess work and storage requirements as a function of wavelength and simulation time.展开更多
从TTI介质一阶应力—速度方程出发,利用旋转交错网格高阶有限差分方法,将非分裂完全匹配层(Non-spliting Perfect Match Layer,简称NPML)边界吸收条件和自由边界条件相结合形成组合边界条件,进行了二维三分量TTI介质弹性波场数值模拟。...从TTI介质一阶应力—速度方程出发,利用旋转交错网格高阶有限差分方法,将非分裂完全匹配层(Non-spliting Perfect Match Layer,简称NPML)边界吸收条件和自由边界条件相结合形成组合边界条件,进行了二维三分量TTI介质弹性波场数值模拟。波场快照和炮记录表明:①采用非分裂式边界条件能较好地消除近地表大角度入射波和瞬逝波;②组合边界条件与NPML边界吸收条件相比,不仅有效地压制了边界反射,同时实现了对自由地表的模拟,获得了丰富的全波场信息,其中在地表产生的PS转换横波作为一种特殊的横波现象,可为近地表结构调查以及多波波场分析等提供有益信息;③自由地表引起的面波以及多次波对偏移结果有着重要影响,因此在实际地震资料处理中应当充分考虑自由地表条件对波场的影响效应。数值模拟结果证实了组合边界条件下二维三分量TTI介质波场数值模拟方法的可行性和正确性。展开更多
基金Supported by the National Natural Science Foundation of China
文摘The key problem of finite-difference time-domain (FD-TD) method is the skillful application of special conditions on the boundaries of the computational domain. A new technique named Perfectly Matched Layer(PML) yields a robust Absorbing Boundary Condition(ABC) independent of the angle of incidence and the frequency of outgoing waves. In this paper, the principle of the PML technique is briefly presented. Then some problems in the application and their settlements are discussed emphatically. Finally three numerical tests and a measured result are devoted to examine the accuracy and effectiveness of this approach.
基金National Natural Science Foundation of China (50608024 and 50538050).
文摘The theory of perfectly matched layer (PML) artificial boundary condition (ABC), which is characterized by absorption any wave motions with arbitrary frequency and arbitrarily incident angle, is introduced. The construction process of PML boundary based on elastodynamic partial differential equation (PDE) system is developed. Combining with velocity-stress hybrid finite element formulation, the applicability of PML boundary is investigated and the numerical reflection of PML boundary is estimated. The reflectivity of PML and multi-transmitting formula (MTF) boundary is then compared based on body wave and surface wave simulations. The results show that although PML boundary yields some reflection, its absorption performance is superior to MTF boundary in the numerical simulations of near-fault wave propagation, especially in comer and large angle grazing incidence situations. The PML boundary does not arise any unstable phenomenon and the stability of PML boundary is better than MTF boundary in hybrid finite element method. For a specified problem and analysis tolerance, the computational efficiency of PML boundary is only a little lower than MTF boundary.
基金supported by the Foundation of State Key Laboratory of Acoustics,Institute of Acoustics,Chinese Academy of Sciences(No.SKLA201303)the National Natural Science Foundation of China(Nos.11104044,11234002,and 11474073)
文摘The perfectly matched layer(PML) is an effective technique for truncating unbounded domains with minimal spurious reflections. A fluid parabolic equation(PE) model applying PML technique was previously used to analyze the sound propagation problem in a range-dependent waveguide(Lu and Zhu, 2007). However, Lu and Zhu only considered a standard fluid PE to demonstrate the capability of the PML and did not take improved one-way models into consideration. They applied a [1/1] Padé approximant to the parabolic equation. The higher-order PEs are more accurate than standard ones when a very large angle propagation is considered. As for range-dependent problems, the techniques to handle the vertical interface between adjacent regions are mainly energy conserving and single-scattering. In this paper, the PML technique is generalized to the higher order elastic PE, as is to the higher order fluid PE. The correction of energy conserving is used in range-dependent waveguides. Simulation is made in both acoustic cases and seismo-acoustic cases. Range-independent and range-dependent waveguides are both adopted to test the accuracy and efficiency of this method. The numerical results illustrate that a PML is much more effective than an artificial absorbing layer(ABL) both in acoustic and seismo-acoustic sound propagation modeling.
基金supported by the National Science and Technology Major Special Sub-project of China(No.2016ZX05024-001-008)the National Natural Science Foundation Joint Fund Prcject of China(No.U1562215).
文摘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.
基金the National Key R&D Program of China(No.2018YFC1504204)National Natural Science Foundation of China(No.U1901602)+3 种基金Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory(Guangzhou)(No.GML2019ZD0203)Shenzhen Science and Technology Program(No.KQTD20170810111725321)Shenzhen Key Laboratory of Deep Offshore Oil and Gas Exploration Technology(No.ZDSYS20190902093007855)the Center for Computational Science and Engineering of Southern University of Science and Technology.
文摘A good artificial boundary treatment in a seismic wave grid-based numerical simulation can reduce the size of the computational region and increase the computational efficiency,which is becoming increasingly important for seismic migration and waveform inversion tasks requiring hundreds or thousands of simulations.Two artificial boundary techniques are commonly used:perfectly matched layers(PMLs),which exhibit the excellent absorption performance but impose a greater computational burden by using finite layers to gradually reduce wave amplitudes;and absorbing boundary conditions(ABCs),which have the high computational efficiency but are less effective in absorption because they employ the one-way wave equation at the exterior boundary.Naturally,PMLs have been combined with ABCs to reduce the number of PMLs,thus improving the computational efficiency;many studies have proposed such hybrid PMLs.Depending on the equations from which the ABCs are derived,there are two hybrid PML variants:the PML+unstretched ABC(UABC),in which the ABC is derived from a physical equation;or the PML+stretched ABC(SABC),in which the ABC is derived from the PML equation.Even though all the previous studies concluded that hybrid PMLs can improve the absorption performance,none of them quantified how many PMLs can be removed by combining the PML with the ABC compared with the pure PML.In this paper,we systematically study the absorption performance of the two hybrid PML variants.We develop a method to distinguish the artificial reflections from the PML-interior interface and those caused by the PML exterior boundary to accurately approximate the additional absorption achieved by using the UABC and the SABC.The reflection coefficients based on a theoretical derivation and numerical tests both show that the UABC amplifies most reflections and is not recommended in any situation;conversely,the SABC can always diminish reflections,but the additional absorption achieved by the SABC is relatively poor and cannot effectively reduce the number of PMLs.In contrast,we find that simply increasing the damping parameter improves absorption better than the PML+SABC.Our results show that the improvement in absorption achieved by combining the PML with either the SABC or the UABC is not better than that obtained by simply adjusting the damping profile of the PML;thus,combining the PML with the ABC is not recommended in practice.
文摘An FDTD system associated with uniaxial perfectly matched layer(UPML) for an electromagnetic scattering problem in two-dimensional space in polar coordinates is considered.Particularly the FDTD system of an initial-boundary value problems of the transverse magnetic(TM) mode to Maxwell's equations is obtained by Yee's algorithm,and the open domain of the scattering problem is truncated by a circle with a UPML.Besides,an artificial boundary condition is imposed on the outer boundary of the UPML.Afterwards,stability of the FDTD system on the truncated domain is established through energy estimates by the Gronwall inequality.Numerical experiments are designed to approve the theoretical analysis.
基金supported by China NSF grants Nos.11771057,11401040 and 11671052.The research of X.J.Li is supported by China NSF grant Nos.11805049 and 11771440 and by the National Magnetic Confinement Fusion Science Program No.2015GB110003.
文摘In this paper,we study the two-dimensional Helmholtz scattering problem by a locally perturbed line with impedance boundary condition.Different from the problem with Dirichlet boundary condition,the Green function of the Helmholtz equation with impedance boundary condition becomes very complicated and comprises surface waves along the locally perturbed line.A uniaxial perfectly matched layer(UPML)method is proposed to truncate the half plane into a bounded computational domain.The main contribution of this paper is to prove the well-posedness of the PML problem and the exponential convergence of the approximate solution to the exact solution as either the thickness or the medium parameter of PML increases.
文摘We review time-domain formulations of radiation boundary conditions for Maxwell's equations, focusing on methods which can deliver arbitrary accuracy at acceptable compu-tational cost. Examples include fast evaluations of nonlocal conditions on symmetric and general boundaries, methods based on identifying and evaluating equivalent sources, and local approximations such as the perfectly matched layer and sequences of local boundary conditions. Complexity estimates are derived to assess work and storage requirements as a function of wavelength and simulation time.
文摘从TTI介质一阶应力—速度方程出发,利用旋转交错网格高阶有限差分方法,将非分裂完全匹配层(Non-spliting Perfect Match Layer,简称NPML)边界吸收条件和自由边界条件相结合形成组合边界条件,进行了二维三分量TTI介质弹性波场数值模拟。波场快照和炮记录表明:①采用非分裂式边界条件能较好地消除近地表大角度入射波和瞬逝波;②组合边界条件与NPML边界吸收条件相比,不仅有效地压制了边界反射,同时实现了对自由地表的模拟,获得了丰富的全波场信息,其中在地表产生的PS转换横波作为一种特殊的横波现象,可为近地表结构调查以及多波波场分析等提供有益信息;③自由地表引起的面波以及多次波对偏移结果有着重要影响,因此在实际地震资料处理中应当充分考虑自由地表条件对波场的影响效应。数值模拟结果证实了组合边界条件下二维三分量TTI介质波场数值模拟方法的可行性和正确性。