Properties of symmetrical layers as matching layers in multilayer thin film design were analyzed. A calculation method was presented to derive parameters of desired equivalent refractive index. A harmonic beam splitte...Properties of symmetrical layers as matching layers in multilayer thin film design were analyzed. A calculation method was presented to derive parameters of desired equivalent refractive index. A harmonic beam splitter was designed and fabricated to test this matching method. OCIS codes: 230.1360, 220.0220, 310.6860.展开更多
In numerical simulation of wave propagation,both viscoelastic materials and perfectly matched layers(PMLs)attenuate waves.The wave equations for both the viscoelastic model and the PML contain convolution operators.Ho...In numerical simulation of wave propagation,both viscoelastic materials and perfectly matched layers(PMLs)attenuate waves.The wave equations for both the viscoelastic model and the PML contain convolution operators.However,convolution operator is intractable in finite-difference time-domain(FDTD)method.A great deal of progress has been made in using time stepping instead of convolution in FDTD.To incorporate PML into viscoelastic media,more memory variables need to be introduced,which increases the code complexity and computation costs.By modifying the nonsplitting PML formulation,I propose a viscoelastic model,which can be used as a viscoelastic material and/or a PML just by adjusting the parameters.The proposed viscoelastic model is essentially equivalent to a Maxwell model.Compared with existing PML methods,the proposed method requires less memory and its implementation in existing finite-difference codes is much easier.The attenuation and phase velocity of P-and S-waves are frequency independent in the viscoelastic model if the related quality factors(Q)are greater than 10.The numerical examples show that the method is stable for materials with high absorption(Q=1),and for heterogeneous media with large contrast of acoustic impedance and large contrast of viscosity.展开更多
It is an important issue to numerically solve the time fractional Schrödinger equation on unbounded domains, which models the dynamics of optical solitons propagating via optical fibers. The perfectly matched lay...It is an important issue to numerically solve the time fractional Schrödinger equation on unbounded domains, which models the dynamics of optical solitons propagating via optical fibers. The perfectly matched layer approach is applied to truncate the unbounded physical domain, and obtain an initial boundary value problem on a bounded computational domain, which can be efficiently solved by the finite difference method. The stability of the reduced initial boundary value problem is rigorously analyzed. Some numerical results are presented to illustrate the accuracy and feasibility of the perfectly matched layer approach. According to these examples, the absorption parameters and the width of the absorption layer will affect the absorption effect. The larger the absorption width, the better the absorption effect. There is an optimal absorption parameter, the absorption effect is the best.展开更多
In this work,a numerical study of the effects of soil-structure interaction(SSI)and granular material-structure interaction(GSI)on the nonlinear response and seismic capacity of flat-bottomed storage silos is conducte...In this work,a numerical study of the effects of soil-structure interaction(SSI)and granular material-structure interaction(GSI)on the nonlinear response and seismic capacity of flat-bottomed storage silos is conducted.A series of incremental dynamic analyses(IDA)are performed on a case of large reinforced concrete silo using 10 seismic recordings.The IDA results are given by two average IDA capacity curves,which are represented,as well as the seismic capacity of the studied structure,with and without a consideration of the SSI while accounting for the effect of GSI.These curves are used to quantify and evaluate the damage of the studied silo by utilizing two damage indices,one based on dissipated energy and the other on displacement and dissipated energy.The cumulative energy dissipation curves obtained by the average IDA capacity curves with and without SSI are presented as a function of the base shear,and these curves allow one to obtain the two critical points and the different limit states of the structure.It is observed that the SSI and GSI significantly influence the seismic response and capacity of the studied structure,particularly at higher levels of PGA.Moreover,the effect of the SSI reduces the damage index of the studied structure by 4%.展开更多
Large calculation error can be formed by directly employing the conventional Yee’s grid to curve surfaces.In order to alleviate such condition,unconditionally stable CrankNicolson Douglas-Gunn(CNDG)algorithm with is ...Large calculation error can be formed by directly employing the conventional Yee’s grid to curve surfaces.In order to alleviate such condition,unconditionally stable CrankNicolson Douglas-Gunn(CNDG)algorithm with is proposed for rotationally symmetric multi-scale problems in anisotropic magnetized plasma.Within the CNDG algorithm,an alternative scheme for the simulation of anisotropic plasma is proposed in body-of-revolution domains.Convolutional perfectly matched layer(CPML)formulation is proposed to efficiently solve the open region problems.Numerical example is carried out for the illustration of effectiveness including the efficiency,resources,and absorption.Through the results,it can be concluded that the proposed scheme shows considerable performance during the simulation.展开更多
We develop a uniaxial optimal perfectly matched layer (opt PML) method for solving the time-harmonic scattering problems by choosing a particular absorbing function with unbounded integral in a rectangular domain. W...We develop a uniaxial optimal perfectly matched layer (opt PML) method for solving the time-harmonic scattering problems by choosing a particular absorbing function with unbounded integral in a rectangular domain. With this choice, the solution of the optimal PML problem not only converges exponentially to the solution of the original scatting problem, but also is insensitive to the thickness of the PML layer for sufficiently small parameter ε0. Numerical experiments are included to illustrate the competitive behavior of the proposed optimal method.展开更多
Alternating direction implicit finite difference time domain (ADI-FDTD) method is unconditionally stable and the maximum time step is not limited by the Courant stability condition, but rather by numerical error. Co...Alternating direction implicit finite difference time domain (ADI-FDTD) method is unconditionally stable and the maximum time step is not limited by the Courant stability condition, but rather by numerical error. Compared with the conventional FDTD method, the time step of ADI-FDTD can be enlarged arbitrarily and the CPU cost can be reduced. 2D perfectly matched layer (PML) absorbing boundary condition is proposed to truncate computation space for ADI-FDTD in dispersive media using recursive convolution(RC) method and the 2D PML formulations for dispersive media are derived. ADI-FDTD formulations for dispersive media can be obtained from the simplified PML formulations. The scattering of target in dispersive soil is simulated under sine wave and Gaussian pulse excitations and numerical results of ADI-FDTD with PML are compared with FDTD. Good agreement is observed. At the same time the CPU cost for ADI-FDTD is obviously reduced.展开更多
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.展开更多
Despite of the limitation in modeling infinite space, the finite element method(FEM) is one of the most used tools to numerically study the geotechnical problems regarding the capacity of simulating different geometri...Despite of the limitation in modeling infinite space, the finite element method(FEM) is one of the most used tools to numerically study the geotechnical problems regarding the capacity of simulating different geometries, conditions and material behaviors. A kind of absorbing layer named perfectly matched layer(PML) has been applied to modeling the radiation damping using FEM, which makes the dynamic analysis of soil-structure interaction more accurate. The PML is capable of absorbing incident waves under any angle and frequency, ensuring them to pass through the model boundaries without reflection.In this context, a new FEM program has been written and the PML formula has been implemented by rewriting the dynamic equation of motion and deriving new properties for the quadrilateral elements.The analysis of soil-foundation interaction by applying the PML is validated by the evaluation of impedance/compliance functions for different ground conditions. The results obtained from the PML model match the extended mesh results, even though the domain is small enough that other types of absorbing boundaries can reflect waves back to the foundation. The mechanism of the wave propagation in the region shows that the forced vibrations can be fully absorbed and damped by the boundaries surrounded by PMLs which is the role of radiation damping in FEM modeling.展开更多
The characterization of finite length Surface Acoustic Wave (SAW) and Bulk acoustic Wave (BAW) resonators is addressed here. The Finite Element Analysis (FEA) induces artificial wave reflections at the edges of the me...The characterization of finite length Surface Acoustic Wave (SAW) and Bulk acoustic Wave (BAW) resonators is addressed here. The Finite Element Analysis (FEA) induces artificial wave reflections at the edges of the mesh. In fact, these ones do not contribute in practice to the corresponding experimental response. The Perfectly Matched Layer (PML) method, allows to suppress the boundary reflections. In this work, we first demonstrate the basis of PML adapted to FEA formalism. Next, the results of such a method are depicted allowing a discussion on the behavior of finite acoustic resonators.展开更多
A higher-order finite-difference time-domain(HO-FDTD) in the spherical coordinate is presented in this paper. The stability and dispersion properties of the proposed scheme are investigated and an air-filled spheric...A higher-order finite-difference time-domain(HO-FDTD) in the spherical coordinate is presented in this paper. The stability and dispersion properties of the proposed scheme are investigated and an air-filled spherical resonator is modeled in order to demonstrate the advantage of this scheme over the finite-difference time-domain(FDTD) and the multiresolution time-domain(MRTD) schemes with respect to memory requirements and CPU time. Moreover, the Berenger's perfectly matched layer(PML) is derived for the spherical HO-FDTD grids, and the numerical results validate the efficiency of the PML.展开更多
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.展开更多
A new computational algorithm is introduced for solving scattering problem in periodic structure. The PML technique is used to deal with the difficulty on truncating the unbounded domain while the DSC algorithm is uti...A new computational algorithm is introduced for solving scattering problem in periodic structure. The PML technique is used to deal with the difficulty on truncating the unbounded domain while the DSC algorithm is utilized for the spatial discretization. The present study reveals that the method is efficient for solving the problem.展开更多
A novel 212DFDTD method aiming to deal with dispersive media is developed. It is named as 212D(FD)2TD method. The main idea of it is to represent the constitutive equation by a convolution in time domain. Consequently...A novel 212DFDTD method aiming to deal with dispersive media is developed. It is named as 212D(FD)2TD method. The main idea of it is to represent the constitutive equation by a convolution in time domain. Consequently, this newly developed method is used to analyze some of the characteristics of an ideal line antenna in pure water. A good agreement with the former experimented results is obtained.展开更多
Composite electromagnetic scattering from a two-dimensional (2D) ship-like target on a one-dimensional sea surface is investigated by using the finite-difference time-domain (FDTD) method. A uniaxial perfectly mat...Composite electromagnetic scattering from a two-dimensional (2D) ship-like target on a one-dimensional sea surface is investigated by using the finite-difference time-domain (FDTD) method. A uniaxial perfectly matched layer is adopted for truncation of FDTD lattices.The FDTD updated equations can be used for the total computation domain by choosing the uniaxial parameters properly. To validate the proposed numerical technique,a 2D infinitely long cylinder over the sea surface is taken into account first.The variation of angular distribution of the scattering changing with incident angle is calculated. The results show good agreement with the conventional moment method. Finally,the influence of the incident angle,the polarization,and the size of the ship-like target on the composite scattering coefficient is discussed in detail.展开更多
In this paper, we consider electromagnetic scattering problems for two-dimensional overfilled cavities. A half ringy absorbing perfectly matched layer (PML) is introduced to enclose the cavity, and the PML formulati...In this paper, we consider electromagnetic scattering problems for two-dimensional overfilled cavities. A half ringy absorbing perfectly matched layer (PML) is introduced to enclose the cavity, and the PML formulations for both TM and TE polarizations are presented. Existence, uniqueness and convergence of the PML solutions are considered. Numerical experiments demonstrate that the PML method is efficient and accurate for solving cavity scattering problems.展开更多
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.展开更多
To efficiently simulate and calculate the radar cross section(RCS) related electromagnetic problems by employing the finite-difference time-domain(FDTD) algorithm, an efficient stretched coordinate perfectly matched l...To efficiently simulate and calculate the radar cross section(RCS) related electromagnetic problems by employing the finite-difference time-domain(FDTD) algorithm, an efficient stretched coordinate perfectly matched layer(ESC-PML) based upon the exponential time differencing(ETD) method is proposed.The proposed implementation can not only reduce the number of auxiliary variables in the SC-PML regions but also maintain the ability of the original SC-PML in terms of the absorbing performance. Compared with the other existed algorithms, the ETDFDTD method shows the least memory consumption resulting in the computational efficiency. The effectiveness and efficiency of the proposed ESC-PML scheme is verified through the RCS relevant problems including the perfect E conductor(PEC) sphere model and the patch antenna model. The results indicate that the proposed scheme has the advantages of the ETD-FDTD method and ESC-PML scheme in terms of high computational efficiency and considerable computational accuracy.展开更多
Spectral element method(SEM) for elastic media is well known for its great flexibility and high accuracy in solving problems with complex geometries.It is an advanced choice for wave simulations.Due to anelasticity ...Spectral element method(SEM) for elastic media is well known for its great flexibility and high accuracy in solving problems with complex geometries.It is an advanced choice for wave simulations.Due to anelasticity of earth media,SEM for elastic media is no longer appropriate.On fundamental of the second-order elastic SEM,this work takes the viscoelastic wave equations and the vertical transversely isotropic(VTI) media into consideration,and establishes the second-order SEM for wave modeling in viscoelastic VTI media.The second-order perfectly matched layer for viscoelastic VTI media is also introduced.The problem of handling the overlapped absorbed corners is solved.A comparison with the analytical solution in a twodimensional viscoelastic homogeneous medium shows that the method is accurate in the wave-field modeling.Furtherly,numerical validation also presents its great flexibility in solving wave propagation problems in complex heterogeneous media.This second-order SEM with perfectly matched layer for viscoelastic VTI media can be easily applied in wave modeling in a limited region.展开更多
In order to simulate metamaterial rotational symmetric open region problems,unconditionally stable perfectly match layer(PML)implementation is proposed in the body of revolution(BOR)finite-difference time-domain(FDTD)...In order to simulate metamaterial rotational symmetric open region problems,unconditionally stable perfectly match layer(PML)implementation is proposed in the body of revolution(BOR)finite-difference time-domain(FDTD)lattice.More precisely,the proposed algorithm is implemented by the Crank-Nicolson(CN)Douglas-Gunn(DG)procedure for BOR metamaterial simulation.The constitutive relationship of metamaterial can be expressed by the Drude model and calculated by the piecewise linear recursive convolution(PLRC)approach.The effectiveness including absorption,efficiency,and accuracy is demonstrated through the numerical example.It can be concluded that the proposed implementation is to take the advantages of the CNDG-PML procedure,PLRC approach,and BORFDTD algorithm in terms of considerable accuracy,enhanced absorption and remarkable efficiency.Meanwhile,it can be demonstrated that the proposed scheme can maintain its unconditional stability when the time step exceeds the CourantFriedrichs-Levy(CFL)condition.展开更多
文摘Properties of symmetrical layers as matching layers in multilayer thin film design were analyzed. A calculation method was presented to derive parameters of desired equivalent refractive index. A harmonic beam splitter was designed and fabricated to test this matching method. OCIS codes: 230.1360, 220.0220, 310.6860.
文摘In numerical simulation of wave propagation,both viscoelastic materials and perfectly matched layers(PMLs)attenuate waves.The wave equations for both the viscoelastic model and the PML contain convolution operators.However,convolution operator is intractable in finite-difference time-domain(FDTD)method.A great deal of progress has been made in using time stepping instead of convolution in FDTD.To incorporate PML into viscoelastic media,more memory variables need to be introduced,which increases the code complexity and computation costs.By modifying the nonsplitting PML formulation,I propose a viscoelastic model,which can be used as a viscoelastic material and/or a PML just by adjusting the parameters.The proposed viscoelastic model is essentially equivalent to a Maxwell model.Compared with existing PML methods,the proposed method requires less memory and its implementation in existing finite-difference codes is much easier.The attenuation and phase velocity of P-and S-waves are frequency independent in the viscoelastic model if the related quality factors(Q)are greater than 10.The numerical examples show that the method is stable for materials with high absorption(Q=1),and for heterogeneous media with large contrast of acoustic impedance and large contrast of viscosity.
文摘It is an important issue to numerically solve the time fractional Schrödinger equation on unbounded domains, which models the dynamics of optical solitons propagating via optical fibers. The perfectly matched layer approach is applied to truncate the unbounded physical domain, and obtain an initial boundary value problem on a bounded computational domain, which can be efficiently solved by the finite difference method. The stability of the reduced initial boundary value problem is rigorously analyzed. Some numerical results are presented to illustrate the accuracy and feasibility of the perfectly matched layer approach. According to these examples, the absorption parameters and the width of the absorption layer will affect the absorption effect. The larger the absorption width, the better the absorption effect. There is an optimal absorption parameter, the absorption effect is the best.
文摘In this work,a numerical study of the effects of soil-structure interaction(SSI)and granular material-structure interaction(GSI)on the nonlinear response and seismic capacity of flat-bottomed storage silos is conducted.A series of incremental dynamic analyses(IDA)are performed on a case of large reinforced concrete silo using 10 seismic recordings.The IDA results are given by two average IDA capacity curves,which are represented,as well as the seismic capacity of the studied structure,with and without a consideration of the SSI while accounting for the effect of GSI.These curves are used to quantify and evaluate the damage of the studied silo by utilizing two damage indices,one based on dissipated energy and the other on displacement and dissipated energy.The cumulative energy dissipation curves obtained by the average IDA capacity curves with and without SSI are presented as a function of the base shear,and these curves allow one to obtain the two critical points and the different limit states of the structure.It is observed that the SSI and GSI significantly influence the seismic response and capacity of the studied structure,particularly at higher levels of PGA.Moreover,the effect of the SSI reduces the damage index of the studied structure by 4%.
文摘Large calculation error can be formed by directly employing the conventional Yee’s grid to curve surfaces.In order to alleviate such condition,unconditionally stable CrankNicolson Douglas-Gunn(CNDG)algorithm with is proposed for rotationally symmetric multi-scale problems in anisotropic magnetized plasma.Within the CNDG algorithm,an alternative scheme for the simulation of anisotropic plasma is proposed in body-of-revolution domains.Convolutional perfectly matched layer(CPML)formulation is proposed to efficiently solve the open region problems.Numerical example is carried out for the illustration of effectiveness including the efficiency,resources,and absorption.Through the results,it can be concluded that the proposed scheme shows considerable performance during the simulation.
基金The Major State Research Development Program (2005CB321701) of Chinathe NSF(10801063) of China
文摘We develop a uniaxial optimal perfectly matched layer (opt PML) method for solving the time-harmonic scattering problems by choosing a particular absorbing function with unbounded integral in a rectangular domain. With this choice, the solution of the optimal PML problem not only converges exponentially to the solution of the original scatting problem, but also is insensitive to the thickness of the PML layer for sufficiently small parameter ε0. Numerical experiments are included to illustrate the competitive behavior of the proposed optimal method.
文摘Alternating direction implicit finite difference time domain (ADI-FDTD) method is unconditionally stable and the maximum time step is not limited by the Courant stability condition, but rather by numerical error. Compared with the conventional FDTD method, the time step of ADI-FDTD can be enlarged arbitrarily and the CPU cost can be reduced. 2D perfectly matched layer (PML) absorbing boundary condition is proposed to truncate computation space for ADI-FDTD in dispersive media using recursive convolution(RC) method and the 2D PML formulations for dispersive media are derived. ADI-FDTD formulations for dispersive media can be obtained from the simplified PML formulations. The scattering of target in dispersive soil is simulated under sine wave and Gaussian pulse excitations and numerical results of ADI-FDTD with PML are compared with FDTD. Good agreement is observed. At the same time the CPU cost for ADI-FDTD is obviously reduced.
基金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.
文摘Despite of the limitation in modeling infinite space, the finite element method(FEM) is one of the most used tools to numerically study the geotechnical problems regarding the capacity of simulating different geometries, conditions and material behaviors. A kind of absorbing layer named perfectly matched layer(PML) has been applied to modeling the radiation damping using FEM, which makes the dynamic analysis of soil-structure interaction more accurate. The PML is capable of absorbing incident waves under any angle and frequency, ensuring them to pass through the model boundaries without reflection.In this context, a new FEM program has been written and the PML formula has been implemented by rewriting the dynamic equation of motion and deriving new properties for the quadrilateral elements.The analysis of soil-foundation interaction by applying the PML is validated by the evaluation of impedance/compliance functions for different ground conditions. The results obtained from the PML model match the extended mesh results, even though the domain is small enough that other types of absorbing boundaries can reflect waves back to the foundation. The mechanism of the wave propagation in the region shows that the forced vibrations can be fully absorbed and damped by the boundaries surrounded by PMLs which is the role of radiation damping in FEM modeling.
文摘The characterization of finite length Surface Acoustic Wave (SAW) and Bulk acoustic Wave (BAW) resonators is addressed here. The Finite Element Analysis (FEA) induces artificial wave reflections at the edges of the mesh. In fact, these ones do not contribute in practice to the corresponding experimental response. The Perfectly Matched Layer (PML) method, allows to suppress the boundary reflections. In this work, we first demonstrate the basis of PML adapted to FEA formalism. Next, the results of such a method are depicted allowing a discussion on the behavior of finite acoustic resonators.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.61301063 and 41305017)
文摘A higher-order finite-difference time-domain(HO-FDTD) in the spherical coordinate is presented in this paper. The stability and dispersion properties of the proposed scheme are investigated and an air-filled spherical resonator is modeled in order to demonstrate the advantage of this scheme over the finite-difference time-domain(FDTD) and the multiresolution time-domain(MRTD) schemes with respect to memory requirements and CPU time. Moreover, the Berenger's perfectly matched layer(PML) is derived for the spherical HO-FDTD grids, and the numerical results validate the efficiency of the PML.
基金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 NNSF of China(10626017)the Science Foundation of the Education Committee of Heilongjiang Province(11511276)the Foundation of Heilongjiang Province(LBH-Q05114).
文摘A new computational algorithm is introduced for solving scattering problem in periodic structure. The PML technique is used to deal with the difficulty on truncating the unbounded domain while the DSC algorithm is utilized for the spatial discretization. The present study reveals that the method is efficient for solving the problem.
文摘A novel 212DFDTD method aiming to deal with dispersive media is developed. It is named as 212D(FD)2TD method. The main idea of it is to represent the constitutive equation by a convolution in time domain. Consequently, this newly developed method is used to analyze some of the characteristics of an ideal line antenna in pure water. A good agreement with the former experimented results is obtained.
基金Project supported by the National Natural Science Foundation of China (Grant No 60571058)the Specialized Research Fund for the Doctoral Program of Higher Education,China (Grant No 20070701010)
文摘Composite electromagnetic scattering from a two-dimensional (2D) ship-like target on a one-dimensional sea surface is investigated by using the finite-difference time-domain (FDTD) method. A uniaxial perfectly matched layer is adopted for truncation of FDTD lattices.The FDTD updated equations can be used for the total computation domain by choosing the uniaxial parameters properly. To validate the proposed numerical technique,a 2D infinitely long cylinder over the sea surface is taken into account first.The variation of angular distribution of the scattering changing with incident angle is calculated. The results show good agreement with the conventional moment method. Finally,the influence of the incident angle,the polarization,and the size of the ship-like target on the composite scattering coefficient is discussed in detail.
文摘In this paper, we consider electromagnetic scattering problems for two-dimensional overfilled cavities. A half ringy absorbing perfectly matched layer (PML) is introduced to enclose the cavity, and the PML formulations for both TM and TE polarizations are presented. Existence, uniqueness and convergence of the PML solutions are considered. Numerical experiments demonstrate that the PML method is efficient and accurate for solving cavity scattering problems.
文摘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 the National Natural Science Foundation of China(61571022611971022)。
文摘To efficiently simulate and calculate the radar cross section(RCS) related electromagnetic problems by employing the finite-difference time-domain(FDTD) algorithm, an efficient stretched coordinate perfectly matched layer(ESC-PML) based upon the exponential time differencing(ETD) method is proposed.The proposed implementation can not only reduce the number of auxiliary variables in the SC-PML regions but also maintain the ability of the original SC-PML in terms of the absorbing performance. Compared with the other existed algorithms, the ETDFDTD method shows the least memory consumption resulting in the computational efficiency. The effectiveness and efficiency of the proposed ESC-PML scheme is verified through the RCS relevant problems including the perfect E conductor(PEC) sphere model and the patch antenna model. The results indicate that the proposed scheme has the advantages of the ETD-FDTD method and ESC-PML scheme in terms of high computational efficiency and considerable computational accuracy.
基金financially supported by the National Natural Science Foundation of China (Grant No.41304077)Postdoctoral Science Foundation of China (Grant No.2013M531744,2014T70740)+1 种基金Key Laboratory of Geospace Environment and Geodesy (Grant No.12-02-03)Subsurface Multi-scale Imaging Laboratory (Grant No.SMIL-2014-01)
文摘Spectral element method(SEM) for elastic media is well known for its great flexibility and high accuracy in solving problems with complex geometries.It is an advanced choice for wave simulations.Due to anelasticity of earth media,SEM for elastic media is no longer appropriate.On fundamental of the second-order elastic SEM,this work takes the viscoelastic wave equations and the vertical transversely isotropic(VTI) media into consideration,and establishes the second-order SEM for wave modeling in viscoelastic VTI media.The second-order perfectly matched layer for viscoelastic VTI media is also introduced.The problem of handling the overlapped absorbed corners is solved.A comparison with the analytical solution in a twodimensional viscoelastic homogeneous medium shows that the method is accurate in the wave-field modeling.Furtherly,numerical validation also presents its great flexibility in solving wave propagation problems in complex heterogeneous media.This second-order SEM with perfectly matched layer for viscoelastic VTI media can be easily applied in wave modeling in a limited region.
基金supported by the National Key Laboratory of Science and Technology on Space Microwave(6142411032201)the National Key Research and Development Program of China(2020YFB1807400)+2 种基金the National Natural Science Foundation of China(6157102261971022)the National Key Laboratory Foundation of China(61424020305)。
文摘In order to simulate metamaterial rotational symmetric open region problems,unconditionally stable perfectly match layer(PML)implementation is proposed in the body of revolution(BOR)finite-difference time-domain(FDTD)lattice.More precisely,the proposed algorithm is implemented by the Crank-Nicolson(CN)Douglas-Gunn(DG)procedure for BOR metamaterial simulation.The constitutive relationship of metamaterial can be expressed by the Drude model and calculated by the piecewise linear recursive convolution(PLRC)approach.The effectiveness including absorption,efficiency,and accuracy is demonstrated through the numerical example.It can be concluded that the proposed implementation is to take the advantages of the CNDG-PML procedure,PLRC approach,and BORFDTD algorithm in terms of considerable accuracy,enhanced absorption and remarkable efficiency.Meanwhile,it can be demonstrated that the proposed scheme can maintain its unconditional stability when the time step exceeds the CourantFriedrichs-Levy(CFL)condition.