An innovative local artificial boundary condition is proposed to numerically solve the Cauchy problem of the Klein-Gordon equation in an unbounded domain.Initially,the equation is considered as the axial wave prop-aga...An innovative local artificial boundary condition is proposed to numerically solve the Cauchy problem of the Klein-Gordon equation in an unbounded domain.Initially,the equation is considered as the axial wave prop-agation in a bar supported on a spring foundation.The numerical model is then truncated by replacing the half-infinitely long bar with an equivalent mechanical structure.The effective frequency-dependent stiffness of the half-infinitely long bar is expressed as the sum of rational terms using Pade approximation.For each term,a corresponding substructure composed of dampers and masses is constructed.Finally,the equivalent mechan-ical structure is obtained by parallelly connecting these substructures.The proposed approach can be easily implemented within a standard finite element framework by incorporating additional mass points and damper elements.Numerical examples show that with just a few extra degrees of freedom,the proposed approach effec-tively suppresses artificial reflections at the truncation boundary and exhibits first-order convergence.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.11832001 and 11702046).
文摘An innovative local artificial boundary condition is proposed to numerically solve the Cauchy problem of the Klein-Gordon equation in an unbounded domain.Initially,the equation is considered as the axial wave prop-agation in a bar supported on a spring foundation.The numerical model is then truncated by replacing the half-infinitely long bar with an equivalent mechanical structure.The effective frequency-dependent stiffness of the half-infinitely long bar is expressed as the sum of rational terms using Pade approximation.For each term,a corresponding substructure composed of dampers and masses is constructed.Finally,the equivalent mechan-ical structure is obtained by parallelly connecting these substructures.The proposed approach can be easily implemented within a standard finite element framework by incorporating additional mass points and damper elements.Numerical examples show that with just a few extra degrees of freedom,the proposed approach effec-tively suppresses artificial reflections at the truncation boundary and exhibits first-order convergence.