The present work describes the application of the method of fundamental solutions （MFS） along with the analog equation method （AEM） and radial basis function （RBF） approximation for solving the 2D isotropic and ...The present work describes the application of the method of fundamental solutions （MFS） along with the analog equation method （AEM） and radial basis function （RBF） approximation for solving the 2D isotropic and anisotropic Helmholtz problems with different wave numbers. The AEM is used to convert the original governing equation into the classical Poisson＇s equation, and the MFS and RBF approximations are used to derive the homogeneous and particular solutions, respectively. Finally, the satisfaction of the solution consisting of the homogeneous and particular parts to the related governing equation and boundary conditions can produce a system of linear equations, which can be solved with the singular value decomposition （SVD） technique. In the computation, such crucial factors related to the MFS-RBF as the location of the virtual boundary, the differential and integrating strategies, and the variation of shape parameters in multi-quadric （MQ） are fully analyzed to provide useful reference.展开更多
A novel panel-free approach based on the method of fundamental solutions (MFS) is proposed to solve the potential flow for predicting ship motion responses in the frequency domain according to strip theory. Compared w...A novel panel-free approach based on the method of fundamental solutions (MFS) is proposed to solve the potential flow for predicting ship motion responses in the frequency domain according to strip theory. Compared with the conventional boundary element method (BEM), MFS is a desingularized, panel-free and integration-free approach. As a result, it is mathematically simple and easy for programming. The velocity potential is described by radial basis function (RBF) approximations and any degree of continuity of the velocity potential gradient can be obtained. Desingularization is achieved through collating singularities on a pseudo boundary outside the real fluid domain. Practical implementation and numerical characteristics of the MFS for solving the potential flow problem concerning ship hydrodynamics are elaborated through the computation of a 2D rectangular section. Then, the current method is further integrated with frequency domain strip theory to predict the heave and pitch responses of a containership and a very large crude carrier (VLCC) in regular head waves. The results of both ships agree well with the 3D frequency domain panel method and experimental data. Thus, the correctness and usefulness of the proposed approach are proved. We hope that this paper will serve as a motivation for other researchers to apply the MFS to various challenging problems in the field of ship hydrodynamics.展开更多
A localized space-time method of fundamental solutions(LSTMFS)is extended for solving three-dimensional transient diffusion problems in this paper.The interval segmentation in temporal direction is developed for the a...A localized space-time method of fundamental solutions(LSTMFS)is extended for solving three-dimensional transient diffusion problems in this paper.The interval segmentation in temporal direction is developed for the accurate simulation of long-time-dependent diffusion problems.In the LSTMFS,the whole space-time domain with nodes arranged i divided into a series of overlapping subdomains with a simple geometry.In each subdomain,the conventional method of fundamental solutions is utilized and the coefficients associated with the considered domain can be explicitly determined.By calculating a combined sparse matrix system,the value at any node inside the space-time domain can be obtained.Numerical experi-ments demonstrate that high accuracy and efficiency can be simultaneously achieved via the LSTMFS,even for the problems defined on a long-time and quite complex computational domain.展开更多
In this paper,some remarks for more efficient analysis of two-dimensional elastostatic problems using the method of fundamental solutions are made.First,the effects of the distance between pseudo and main boundaries o...In this paper,some remarks for more efficient analysis of two-dimensional elastostatic problems using the method of fundamental solutions are made.First,the effects of the distance between pseudo and main boundaries on the solution are investigated and by a numerical study a lower bound for the distance of each source point to the main boundary is suggested.In some cases,the resulting system of equations becomes ill-conditioned for which,the truncated singular value decomposition with a criterion based on the accuracy of the imposition of boundary conditions is used.Moreover,a procedure for normalizing the shear modulus is presented that significantly reduces the condition number of the system of equations.By solving two example problems with stress concentration,the effectiveness of the proposed methods is demonstrated.展开更多
In this work,the localized method of fundamental solution(LMFS)is extended to Signorini problem.Unlike the traditional fundamental solution(MFS),the LMFS approximates the field quantity at each node by using the field...In this work,the localized method of fundamental solution(LMFS)is extended to Signorini problem.Unlike the traditional fundamental solution(MFS),the LMFS approximates the field quantity at each node by using the field quantities at the adjacent nodes.The idea of the LMFS is similar to the localized domain type method.The fictitious boundary nodes are proposed to impose the boundary condition and governing equations at each node to formulate a sparse matrix.The inequality boundary condition of Signorini problem is solved indirectly by introducing nonlinear complementarity problem function(NCP-function).Numerical examples are carried out to validate the reliability and effectiveness of the LMFS in solving Signorini problems.展开更多
The elastic plate vibration model is studied under the external force. The size of the source term by the given mode of the source and some observations from the body of the plate is determined over a time interval, w...The elastic plate vibration model is studied under the external force. The size of the source term by the given mode of the source and some observations from the body of the plate is determined over a time interval, which is referred to be an inverse source problem of a plate equation. The uniqueness theorem for this problem is stated, and the fundamental solution to the plate equation is derived. In the case that the plate is driven by the harmonic load, the fundamental solution method （FSM） and the Tikhonov regularization technique axe used to calculate the source term. Numerical experiments of the Euler-Bernoulli beam and the Kirchhoff-Love plate show that the FSM can work well for practical use, no matter the source term is smooth or piecewise.展开更多
The boundary value problem of plate bending problem on two_parameter foundation was discussed.Using two series of the high_order fundamental solution sequences, namely, the fundamental solution sequences for the multi...The boundary value problem of plate bending problem on two_parameter foundation was discussed.Using two series of the high_order fundamental solution sequences, namely, the fundamental solution sequences for the multi_harmonic operator and Laplace operator, applying the multiple reciprocity method(MRM), the MRM boundary integral equation for plate bending problem was constructed. It proves that the boundary integral equation derived from MRM is essentially identical to the conventional boundary integral equation. Hence the convergence analysis of MRM for plate bending problem can be obtained by the error estimation for the conventional boundary integral equation. In addition, this method can extend to the case of more series of the high_order fundamental solution sequences.展开更多
Fundamental solution of Dirichlet boundary value problem of axisymmetric Helmholtz equation is constructed via modi?ed Bessel function of the second kind, which uni?ed the formulas of fundamental solution of Helmholtz...Fundamental solution of Dirichlet boundary value problem of axisymmetric Helmholtz equation is constructed via modi?ed Bessel function of the second kind, which uni?ed the formulas of fundamental solution of Helmholtz equation, elliptic type Euler-Poisson-Darboux equation and Laplace equation in any dimensional space.展开更多
To obtain the fundamental solution of soil has become the key problem for the semi-analytical and semi-numerical (SASN) method in analyzing plate on layered soil. By applying axisymmetric finite element method (FEM),a...To obtain the fundamental solution of soil has become the key problem for the semi-analytical and semi-numerical (SASN) method in analyzing plate on layered soil. By applying axisymmetric finite element method (FEM),an expression relating the surface settlement and the reaction of the layered soil can be obtained. Such a reaction can be treated as load acting on the applied external load. Having the plate modelled by four-node elements,the governing equation of the plate can be formed and solved. In this case, the fundamental solution can be introduced into the global soil stiffness matrix and five-node or nine-node element soil stiffness matrix.The existing commercial FEM software can be used to solve the fundamental solution of soil, which can bypass the complicated formula derivation and boasts high computational efficiency as well.展开更多
A new three-dimensional fundamental solution to the Stokes flow was proposed by transforming the solid harmonic functions in Lamb's solution into expressions in terms Of the oblate spheroidal coordinates. These fu...A new three-dimensional fundamental solution to the Stokes flow was proposed by transforming the solid harmonic functions in Lamb's solution into expressions in terms Of the oblate spheroidal coordinates. These fundamental solutions are advantageous in treating flows past an arbitrary number of arbitrarily positioned and oriented oblate spheroids. The least squares technique was adopted herein so that the convergence difficulties often encountered in solving three-dimensional problems were completely avoided. The examples demonstrate that present approach is highly accurate, consistently stable and computationally efficient. The oblate spheroid may be used to model a variety of particle shapes between a circular disk and a sphere. For the first time, the effect of various geometric factors on the forces and torques exerted on two oblate spheroids were systematically studied by using the proposed fundamental solutions. The generality of this approach was illustrated by two problems of three spheroids.展开更多
Explicit fomulas for 2-D electroelastic fundamental solutions in general anisotropic piezoelectric media subjected to a line force and a line charge are obtained by using the plane wave decomposition method and a subs...Explicit fomulas for 2-D electroelastic fundamental solutions in general anisotropic piezoelectric media subjected to a line force and a line charge are obtained by using the plane wave decomposition method and a subsequent application of the residue calculus. 'Anisotropic' means that any material symmetry restrictions are not assumed. 'Two dimensional' includes not only in-plane problems but also anti-plane problems and problems in which in-plane and anti-plane deformations couple each other. As a special case, the solutions for transversely isotropic piezoelectric media are given.展开更多
It is expected to solve the piezoelectric problems by using numerical method,e.g.,the boundary element method,since analytical solutions can hardly be obtained due to the deformation,electric field coupling effects of...It is expected to solve the piezoelectric problems by using numerical method,e.g.,the boundary element method,since analytical solutions can hardly be obtained due to the deformation,electric field coupling effects of the material and anisotropy of the piezoelectric properties.However,the appropriate fundamental solutions used in boundary element method have not been found so far.In this paper,the solutions of displacement,potential,stress and dielectric displacement field effected by a unit mechanical force and a unit electric charge applied respectively at a point in a two-dimensional infinite domain are derived based on the coupling theory of electromagnetics and continuum mechanics.Hence,the eight indispensable fundamental solutions in BEM are provided.展开更多
We reduce the initial value problem for the generalized Schroedinger equation with piecewise-constant leading coefficient to the system of Volterra type integral equations and construct new useful integral representat...We reduce the initial value problem for the generalized Schroedinger equation with piecewise-constant leading coefficient to the system of Volterra type integral equations and construct new useful integral representations for the fundamental solutions of the Schroedinger equation. We also investigate some significant properties of the kernels of these integral representations. The integral representations of fundamental solutions enable to obtain the basic integral equations, which are a powerful tool for solving inverse spectral problems.展开更多
A fundamental solution was obtained for an infinite plane bonded by two dissimilar isotropic semi-planes by employing plane elastic complex variable method and theory of boundary value problems for analytic functions....A fundamental solution was obtained for an infinite plane bonded by two dissimilar isotropic semi-planes by employing plane elastic complex variable method and theory of boundary value problems for analytic functions.Fundamental solution was prepared for solving these types of problems with boundary element method.展开更多
A weight double trigonometric series is presented as an approximate fundamental solution for orthotropic plate.Integral equation of orthotropic plate bending is solved by a new method, which only needs one basic bound...A weight double trigonometric series is presented as an approximate fundamental solution for orthotropic plate.Integral equation of orthotropic plate bending is solved by a new method, which only needs one basic boundary integral Eq., puts one fictitious boundary outside plate domain. Examples show that the approximate fundamental solution and solving method proposed in this paper are simple, reliable and quite precise. And they are applicable for various boundary conditions.展开更多
The identification of the traction acting on a portion of the surface of an anisotropic solid is very important in structural health monitoring and optimal design of structures. The traction can be determined using in...The identification of the traction acting on a portion of the surface of an anisotropic solid is very important in structural health monitoring and optimal design of structures. The traction can be determined using inverse methods in which displacement or strain measurements are taken at several points on the body. This paper presents an inverse method based on the method of fundamental solutions for the traction identification problem in two-dimensional anisotropic elasticity. The method of fundamental solutions is an efficient boundary-type meshless method widely used for analyzing various problems. Since the problem is linear, the sensitivity analysis is simply performed by solving the corresponding direct problem several times with different loads. The effects of important parameters such as the number of measurement data, the position of the measurement points, the amount of measurement error, and the type of measurement, i.e., displacement or strain, on the results are also investigated. The results obtained show that the presented inverse method is suitable for the problem of traction identification. It can be concluded from the results that the use of strain measurements in the inverse analysis leads to more accurate results than the use of displacement measurements. It is also found that measurement points closer to the boundary with unknown traction provide more reliable solutions. Additionally, it is found that increasing the number of measurement points increases the accuracy of the inverse solution. However, in cases with a large number of measurement points, further increasing the number of measurement data has little effect on the results.展开更多
In this paper,we develop the method of fundamental solutions(MFS)for solving boundary value problems in the field of optical fluorescence.The governing system of diffusion-absorption equations for the excitation and e...In this paper,we develop the method of fundamental solutions(MFS)for solving boundary value problems in the field of optical fluorescence.The governing system of diffusion-absorption equations for the excitation and emission fluences is transformed into a single fourth-order partial differential equation whose fundamental solution can be expressed as the difference of two fundamental solutions of the complex Helmholtz equation.The numerically obtained results confirm the accuracy of the MFS when compared with an available analytical solution.Numerical results are also provided for a physical application in optical fluorescence.Furthermore,extensions to three dimensions along with numerical verification are performed.展开更多
We describe a numerical method for the solution of acoustic exterior scattering problems based on the time-domain boundary integral representation of the solution.As the spatial discretization of the resulting time-do...We describe a numerical method for the solution of acoustic exterior scattering problems based on the time-domain boundary integral representation of the solution.As the spatial discretization of the resulting time-domain boundary integral equation we use either the method of fundamental solutions(MFS)or the Galerkin boundary element method(BEM).In time we apply either a standard convolution quadrature(CQ)based on an A-stable linear multistep method or a modified CQ scheme.It is well-known that the standard low-order CQ schemes for hyperbolic problems suffer from strong dissipation and dispersion properties.The modified scheme is designed to avoid these properties.We give a careful description of the modified scheme and its implementation with differences due to different spatial discretizations highlighted.Numerous numerical experiments illustrate the effectiveness of the modified scheme and dramatic improvement with errors up to two orders of magnitude smaller in comparison with the standard scheme.展开更多
文摘The present work describes the application of the method of fundamental solutions （MFS） along with the analog equation method （AEM） and radial basis function （RBF） approximation for solving the 2D isotropic and anisotropic Helmholtz problems with different wave numbers. The AEM is used to convert the original governing equation into the classical Poisson＇s equation, and the MFS and RBF approximations are used to derive the homogeneous and particular solutions, respectively. Finally, the satisfaction of the solution consisting of the homogeneous and particular parts to the related governing equation and boundary conditions can produce a system of linear equations, which can be solved with the singular value decomposition （SVD） technique. In the computation, such crucial factors related to the MFS-RBF as the location of the virtual boundary, the differential and integrating strategies, and the variation of shape parameters in multi-quadric （MQ） are fully analyzed to provide useful reference.
基金the Fund of the Minister of Education and Minister of Finance of China (No. ZXZY019)
文摘A novel panel-free approach based on the method of fundamental solutions (MFS) is proposed to solve the potential flow for predicting ship motion responses in the frequency domain according to strip theory. Compared with the conventional boundary element method (BEM), MFS is a desingularized, panel-free and integration-free approach. As a result, it is mathematically simple and easy for programming. The velocity potential is described by radial basis function (RBF) approximations and any degree of continuity of the velocity potential gradient can be obtained. Desingularization is achieved through collating singularities on a pseudo boundary outside the real fluid domain. Practical implementation and numerical characteristics of the MFS for solving the potential flow problem concerning ship hydrodynamics are elaborated through the computation of a 2D rectangular section. Then, the current method is further integrated with frequency domain strip theory to predict the heave and pitch responses of a containership and a very large crude carrier (VLCC) in regular head waves. The results of both ships agree well with the 3D frequency domain panel method and experimental data. Thus, the correctness and usefulness of the proposed approach are proved. We hope that this paper will serve as a motivation for other researchers to apply the MFS to various challenging problems in the field of ship hydrodynamics.
基金the Fundamental Research Funds for the Central Universities(Grants B200203009 and B200202126)the Natural Science Foundation of Jiangsu Province(Grant BK20190073)+2 种基金the State Key Laboratory of Acoustics,Chinese Academy of Sciences(Grant SKLA202001)the State Key Laboratory of Mechanical Behavior and System Safety of Traffic Engineering Structures,Shijiazhuang Tiedao University(Grant KF2020-22)the China Postdoctoral Science Foundation(Grants 2017M611669 and 2018T110430).
文摘A localized space-time method of fundamental solutions(LSTMFS)is extended for solving three-dimensional transient diffusion problems in this paper.The interval segmentation in temporal direction is developed for the accurate simulation of long-time-dependent diffusion problems.In the LSTMFS,the whole space-time domain with nodes arranged i divided into a series of overlapping subdomains with a simple geometry.In each subdomain,the conventional method of fundamental solutions is utilized and the coefficients associated with the considered domain can be explicitly determined.By calculating a combined sparse matrix system,the value at any node inside the space-time domain can be obtained.Numerical experi-ments demonstrate that high accuracy and efficiency can be simultaneously achieved via the LSTMFS,even for the problems defined on a long-time and quite complex computational domain.
文摘In this paper,some remarks for more efficient analysis of two-dimensional elastostatic problems using the method of fundamental solutions are made.First,the effects of the distance between pseudo and main boundaries on the solution are investigated and by a numerical study a lower bound for the distance of each source point to the main boundary is suggested.In some cases,the resulting system of equations becomes ill-conditioned for which,the truncated singular value decomposition with a criterion based on the accuracy of the imposition of boundary conditions is used.Moreover,a procedure for normalizing the shear modulus is presented that significantly reduces the condition number of the system of equations.By solving two example problems with stress concentration,the effectiveness of the proposed methods is demonstrated.
基金supported by the National Science Foundation of China(No.52109089)support of Post Doctor Program(2019M652281)Nature Science Foundation of Jiangxi Province(20192BAB216040).
文摘In this work,the localized method of fundamental solution(LMFS)is extended to Signorini problem.Unlike the traditional fundamental solution(MFS),the LMFS approximates the field quantity at each node by using the field quantities at the adjacent nodes.The idea of the LMFS is similar to the localized domain type method.The fictitious boundary nodes are proposed to impose the boundary condition and governing equations at each node to formulate a sparse matrix.The inequality boundary condition of Signorini problem is solved indirectly by introducing nonlinear complementarity problem function(NCP-function).Numerical examples are carried out to validate the reliability and effectiveness of the LMFS in solving Signorini problems.
文摘The elastic plate vibration model is studied under the external force. The size of the source term by the given mode of the source and some observations from the body of the plate is determined over a time interval, which is referred to be an inverse source problem of a plate equation. The uniqueness theorem for this problem is stated, and the fundamental solution to the plate equation is derived. In the case that the plate is driven by the harmonic load, the fundamental solution method （FSM） and the Tikhonov regularization technique axe used to calculate the source term. Numerical experiments of the Euler-Bernoulli beam and the Kirchhoff-Love plate show that the FSM can work well for practical use, no matter the source term is smooth or piecewise.
文摘The boundary value problem of plate bending problem on two_parameter foundation was discussed.Using two series of the high_order fundamental solution sequences, namely, the fundamental solution sequences for the multi_harmonic operator and Laplace operator, applying the multiple reciprocity method(MRM), the MRM boundary integral equation for plate bending problem was constructed. It proves that the boundary integral equation derived from MRM is essentially identical to the conventional boundary integral equation. Hence the convergence analysis of MRM for plate bending problem can be obtained by the error estimation for the conventional boundary integral equation. In addition, this method can extend to the case of more series of the high_order fundamental solution sequences.
基金The NSF(11326152) of Chinathe NSF(BK20130736) of Jiangsu Province of Chinathe NSF(CKJB201709) of Nanjing Institute of Technology
文摘Fundamental solution of Dirichlet boundary value problem of axisymmetric Helmholtz equation is constructed via modi?ed Bessel function of the second kind, which uni?ed the formulas of fundamental solution of Helmholtz equation, elliptic type Euler-Poisson-Darboux equation and Laplace equation in any dimensional space.
文摘To obtain the fundamental solution of soil has become the key problem for the semi-analytical and semi-numerical (SASN) method in analyzing plate on layered soil. By applying axisymmetric finite element method (FEM),an expression relating the surface settlement and the reaction of the layered soil can be obtained. Such a reaction can be treated as load acting on the applied external load. Having the plate modelled by four-node elements,the governing equation of the plate can be formed and solved. In this case, the fundamental solution can be introduced into the global soil stiffness matrix and five-node or nine-node element soil stiffness matrix.The existing commercial FEM software can be used to solve the fundamental solution of soil, which can bypass the complicated formula derivation and boasts high computational efficiency as well.
文摘A new three-dimensional fundamental solution to the Stokes flow was proposed by transforming the solid harmonic functions in Lamb's solution into expressions in terms Of the oblate spheroidal coordinates. These fundamental solutions are advantageous in treating flows past an arbitrary number of arbitrarily positioned and oriented oblate spheroids. The least squares technique was adopted herein so that the convergence difficulties often encountered in solving three-dimensional problems were completely avoided. The examples demonstrate that present approach is highly accurate, consistently stable and computationally efficient. The oblate spheroid may be used to model a variety of particle shapes between a circular disk and a sphere. For the first time, the effect of various geometric factors on the forces and torques exerted on two oblate spheroids were systematically studied by using the proposed fundamental solutions. The generality of this approach was illustrated by two problems of three spheroids.
文摘Explicit fomulas for 2-D electroelastic fundamental solutions in general anisotropic piezoelectric media subjected to a line force and a line charge are obtained by using the plane wave decomposition method and a subsequent application of the residue calculus. 'Anisotropic' means that any material symmetry restrictions are not assumed. 'Two dimensional' includes not only in-plane problems but also anti-plane problems and problems in which in-plane and anti-plane deformations couple each other. As a special case, the solutions for transversely isotropic piezoelectric media are given.
文摘It is expected to solve the piezoelectric problems by using numerical method,e.g.,the boundary element method,since analytical solutions can hardly be obtained due to the deformation,electric field coupling effects of the material and anisotropy of the piezoelectric properties.However,the appropriate fundamental solutions used in boundary element method have not been found so far.In this paper,the solutions of displacement,potential,stress and dielectric displacement field effected by a unit mechanical force and a unit electric charge applied respectively at a point in a two-dimensional infinite domain are derived based on the coupling theory of electromagnetics and continuum mechanics.Hence,the eight indispensable fundamental solutions in BEM are provided.
文摘We reduce the initial value problem for the generalized Schroedinger equation with piecewise-constant leading coefficient to the system of Volterra type integral equations and construct new useful integral representations for the fundamental solutions of the Schroedinger equation. We also investigate some significant properties of the kernels of these integral representations. The integral representations of fundamental solutions enable to obtain the basic integral equations, which are a powerful tool for solving inverse spectral problems.
文摘A fundamental solution was obtained for an infinite plane bonded by two dissimilar isotropic semi-planes by employing plane elastic complex variable method and theory of boundary value problems for analytic functions.Fundamental solution was prepared for solving these types of problems with boundary element method.
文摘A weight double trigonometric series is presented as an approximate fundamental solution for orthotropic plate.Integral equation of orthotropic plate bending is solved by a new method, which only needs one basic boundary integral Eq., puts one fictitious boundary outside plate domain. Examples show that the approximate fundamental solution and solving method proposed in this paper are simple, reliable and quite precise. And they are applicable for various boundary conditions.
基金funded by Vice Chancellor of Research at Shiraz University(grant 3GFU2M1820).
文摘The identification of the traction acting on a portion of the surface of an anisotropic solid is very important in structural health monitoring and optimal design of structures. The traction can be determined using inverse methods in which displacement or strain measurements are taken at several points on the body. This paper presents an inverse method based on the method of fundamental solutions for the traction identification problem in two-dimensional anisotropic elasticity. The method of fundamental solutions is an efficient boundary-type meshless method widely used for analyzing various problems. Since the problem is linear, the sensitivity analysis is simply performed by solving the corresponding direct problem several times with different loads. The effects of important parameters such as the number of measurement data, the position of the measurement points, the amount of measurement error, and the type of measurement, i.e., displacement or strain, on the results are also investigated. The results obtained show that the presented inverse method is suitable for the problem of traction identification. It can be concluded from the results that the use of strain measurements in the inverse analysis leads to more accurate results than the use of displacement measurements. It is also found that measurement points closer to the boundary with unknown traction provide more reliable solutions. Additionally, it is found that increasing the number of measurement points increases the accuracy of the inverse solution. However, in cases with a large number of measurement points, further increasing the number of measurement data has little effect on the results.
文摘In this paper,we develop the method of fundamental solutions(MFS)for solving boundary value problems in the field of optical fluorescence.The governing system of diffusion-absorption equations for the excitation and emission fluences is transformed into a single fourth-order partial differential equation whose fundamental solution can be expressed as the difference of two fundamental solutions of the complex Helmholtz equation.The numerically obtained results confirm the accuracy of the MFS when compared with an available analytical solution.Numerical results are also provided for a physical application in optical fluorescence.Furthermore,extensions to three dimensions along with numerical verification are performed.
文摘We describe a numerical method for the solution of acoustic exterior scattering problems based on the time-domain boundary integral representation of the solution.As the spatial discretization of the resulting time-domain boundary integral equation we use either the method of fundamental solutions(MFS)or the Galerkin boundary element method(BEM).In time we apply either a standard convolution quadrature(CQ)based on an A-stable linear multistep method or a modified CQ scheme.It is well-known that the standard low-order CQ schemes for hyperbolic problems suffer from strong dissipation and dispersion properties.The modified scheme is designed to avoid these properties.We give a careful description of the modified scheme and its implementation with differences due to different spatial discretizations highlighted.Numerous numerical experiments illustrate the effectiveness of the modified scheme and dramatic improvement with errors up to two orders of magnitude smaller in comparison with the standard scheme.
