This paper is concerned with the Cauchy problem for the modified Helmholtz equation in an infinite strip domain 0<x≤1,y∈R.The Cauchy data at x = 0 is given and the solution is then sought for the interval 0<x...This paper is concerned with the Cauchy problem for the modified Helmholtz equation in an infinite strip domain 0<x≤1,y∈R.The Cauchy data at x = 0 is given and the solution is then sought for the interval 0<x≤1.This problem is highly ill-posed and the solution(if it exists) does not depend continuously on the given data. In this paper,we propose a fourth-order modified method to solve the Cauchy problem. Convergence estimates are presented under the suitable choices of regularization parameters and the a priori assumption on the bounds of the exact solution.Numerical implementation is considered and the numerical examples show that our proposed method is effective and stable.展开更多
In this paper, we consider the reconstruction of the wave field in a bounded domain. By choosing a special family of functions, the Cauchy problem can be transformed into a Fourier moment problem. This problem is ill-...In this paper, we consider the reconstruction of the wave field in a bounded domain. By choosing a special family of functions, the Cauchy problem can be transformed into a Fourier moment problem. This problem is ill-posed. We propose a regularization method for obtaining an approximate solution to the wave field on the unspecified boundary. We also give the convergence analysis and error estimate of the numerical algorithm. Finally, we present some numerical examples to show the effectiveness of this method.展开更多
In this paper, we investigate the mixed spectral method using generalized Laguerre functions for exterior problems of fourth order partial differential equations. A mixed spectral scheme is provided for the stream fun...In this paper, we investigate the mixed spectral method using generalized Laguerre functions for exterior problems of fourth order partial differential equations. A mixed spectral scheme is provided for the stream function form of the Navier-Stokes equations outside a disc. Numerical results demonstrate the spectral accuracy in space.展开更多
There are two cases of the exterior problems of the Helmholtz equation. If λ ≥ 0 the bilinear form is coercive, and if λ < 0 it is the scattering problem. We give a new approach of the infinite element method, w...There are two cases of the exterior problems of the Helmholtz equation. If λ ≥ 0 the bilinear form is coercive, and if λ < 0 it is the scattering problem. We give a new approach of the infinite element method, which enables us to solve these exterior problems as well as corner problems. A numerical example of the scattering problem is given. [ABSTRACT FROM AUTHOR]展开更多
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.展开更多
Considering the fundamental solution of the Laplace equation as the weight function,we give the iterative format for solving the nonhomogeneous Helmholtz equation with variable coefficients. Furthermore, the iteratio...Considering the fundamental solution of the Laplace equation as the weight function,we give the iterative format for solving the nonhomogeneous Helmholtz equation with variable coefficients. Furthermore, the iteration method of BEM for solving the equation mentioned above is obtained. The numerical example is given in this paper. Finally, the iteration method of BUM mentioned above is compared with the coupled method of BEM that was presented before then by authors.展开更多
We apply the fast multipole method (FMM) accelerated boundary element method (BEM) for the three-dimensional (3D) Helmholtz equation, and as a result, large-scale acoustic scattering problems involving 400000 elements...We apply the fast multipole method (FMM) accelerated boundary element method (BEM) for the three-dimensional (3D) Helmholtz equation, and as a result, large-scale acoustic scattering problems involving 400000 elements are solved efficiently. This is an extension of the fast multipole BEM for two-dimensional (2D) acoustic problems developed by authors recently. Some new improvements are obtained. In this new technique, the improved Burton-Miller formulation is employed to over-come non-uniqueness difficulties in the conventional BEM for exterior acoustic problems. The computational efficiency is further improved by adopting the FMM and the block diagonal preconditioner used in the generalized minimum residual method (GMRES) iterative solver to solve the system matrix equation. Numerical results clearly demonstrate the complete reliability and efficiency of the proposed algorithm. It is potentially useful for solving large-scale engineering acoustic scattering problems.展开更多
This paper is concerned with the reconstruction of the radiation wave field in the exterior of a bounded two- or three-dimensional domain from the knowledge of Cauchy data on a part of the boundary of the aforemention...This paper is concerned with the reconstruction of the radiation wave field in the exterior of a bounded two- or three-dimensional domain from the knowledge of Cauchy data on a part of the boundary of the aforementioned domain. It is described by the Cauchy problem for the Helmholtz equation. By using the Dirichlet-to-Neumann map, this problem is transformed into an operator equation with compact operator. We rigorously justify the asymptotic behaviors of singular values of the compact operator. Then a projection method with regularization is applied to solve the operator equation, and the convergence of the regularization method is discussed. Finally, several numerical examples are presented to illustrate the approach. The results demonstrate that the algorithm is effective.Mathematics subject classification: 35R25, 35R30, 78A40.展开更多
From the point of view of energy analysis, the cause that the uniqueness of the boundary integral equation induced from the exterior Helmholtz problem does not hold is investigated in this paper. It is proved that the...From the point of view of energy analysis, the cause that the uniqueness of the boundary integral equation induced from the exterior Helmholtz problem does not hold is investigated in this paper. It is proved that the Sommerfeld's condition at the infinity is changed so that it is suitable not only for the radiative wave but also for the absorptive wave when we use the boundary integral equation to describe the exterior Helmholtz problem. There fore, the total energy of the system is conservative. The mathematical dealings to guarantee the uniqueness are discussed based upon this explanation.展开更多
We study the coupled mKdV equation by the dressing method via local Riemann-Hilbert problem. With the help of the Lax pairs, we obtain the matrix Riemann-Hilbert problem with zeros. The explicit solutions for the coup...We study the coupled mKdV equation by the dressing method via local Riemann-Hilbert problem. With the help of the Lax pairs, we obtain the matrix Riemann-Hilbert problem with zeros. The explicit solutions for the coupled mKdV equation are derived with the aid of the regularization of the Riemann-Hilbert problem.展开更多
This paper describes formulation and implementation of the fast multipole boundary element method (FMBEM) for 2D acoustic problems. The kernel function expansion theory is summarized, and four building blocks of the...This paper describes formulation and implementation of the fast multipole boundary element method (FMBEM) for 2D acoustic problems. The kernel function expansion theory is summarized, and four building blocks of the FMBEM are described in details. They are moment calculation, moment to moment translation, moment to local translation, and local to local translation. A data structure for the quad-tree construction is proposed which can facilitate implementation. An analytical moment expression is derived, which is more accurate, stable, and efficient than direct numerical computation. Numerical examples are presented to demonstrate the accuracy and efficiency of the FMBEM, and radiation of a 2D vibration rail mode is simulated using the FMBEM.展开更多
We extend the pure source transfer domain decomposition method(PSTDDM)to solve the perfectly matched layer approximation of Helmholtz scattering problems in heterogeneous media.We first propose some new source transfe...We extend the pure source transfer domain decomposition method(PSTDDM)to solve the perfectly matched layer approximation of Helmholtz scattering problems in heterogeneous media.We first propose some new source transfer operators,and then introduce the layer-wise and block-wise PSTDDMs based on these operators.In particular,it is proved that the solution obtained by the layer-wise PSTDDM in R2 coincides with the exact solution to the heterogeneous Helmholtz problem in the computational domain.Second,we propose the iterative layer-wise and blockwise PSTDDMs,which are designed by simply iterating the PSTDDM alternatively over two staggered decompositions of the computational domain.Finally,extensive numerical tests in two and three dimensions show that,as the preconditioner for the GMRES method,the iterative PSTDDMs are more robust and efficient than PSTDDMs for solving heterogeneous Helmholtz problems.展开更多
In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variat...In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variational formulation and well posedness of the coupling method are obtained. The convergence and optimal estimates for the approximation solution are provided.展开更多
This paper presents a novel wideband fast multipole boundary element approach to 3D half-space/planesymmetric acoustic wave problems.The half-space fundamental solution is employed in the boundary integral equations s...This paper presents a novel wideband fast multipole boundary element approach to 3D half-space/planesymmetric acoustic wave problems.The half-space fundamental solution is employed in the boundary integral equations so that the tree structure required in the fast multipole algorithm is constructed for the boundary elements in the real domain only.Moreover,a set of symmetric relations between the multipole expansion coefficients of the real and image domains are derived,and the half-space fundamental solution is modified for the purpose of applying such relations to avoid calculating,translating and saving the multipole/local expansion coefficients of the image domain.The wideband adaptive multilevel fast multipole algorithm associated with the iterative solver GMRES is employed so that the present method is accurate and efficient for both lowand high-frequency acoustic wave problems.As for exterior acoustic problems,the Burton-Miller method is adopted to tackle the fictitious eigenfrequency problem involved in the conventional boundary integral equation method.Details on the implementation of the present method are described,and numerical examples are given to demonstrate its accuracy and efficiency.展开更多
In this paper, we give a hybrid method to numerically solve the inverse open cavity scattering problem for cavity shape, given the scattered solution on the opening of the cavity. This method is a hybrid between an it...In this paper, we give a hybrid method to numerically solve the inverse open cavity scattering problem for cavity shape, given the scattered solution on the opening of the cavity. This method is a hybrid between an iterative method and an integral equations method for solving the Cauchy problem. The idea of this hybrid method is simple, the operation is easy, and the computation cost is small. Numerical experiments show the feasibility of this method, even for cases with noise.展开更多
We establish variational formulation and prove the existence and uniqueness of the three dimen- sional axisymmetric Stokes exterior problem in weighted spaces. Error estimates and convergence for P2 - P0 elements with...We establish variational formulation and prove the existence and uniqueness of the three dimen- sional axisymmetric Stokes exterior problem in weighted spaces. Error estimates and convergence for P2 - P0 elements with infinite element methods are also obtained. Numerical experiments are presented to verify the theoretical analysis.展开更多
In this paper, we study a class of boundary value problems for conformable fractional differential equations under a new definition. Firstly, by using the monotone iterative technique and the method of coupled upper a...In this paper, we study a class of boundary value problems for conformable fractional differential equations under a new definition. Firstly, by using the monotone iterative technique and the method of coupled upper and lower solution, the sufficient condition for the existence of the boundary value problem is obtained, and the range of the solution is determined. Then the existence and uniqueness of the solution are proved by the proof by contradiction. Finally, a concrete example is given to illustrate the wide applicability of our main results.展开更多
基金supported by the NSF of China(10571079,10671085)and the program of NCET
文摘This paper is concerned with the Cauchy problem for the modified Helmholtz equation in an infinite strip domain 0<x≤1,y∈R.The Cauchy data at x = 0 is given and the solution is then sought for the interval 0<x≤1.This problem is highly ill-posed and the solution(if it exists) does not depend continuously on the given data. In this paper,we propose a fourth-order modified method to solve the Cauchy problem. Convergence estimates are presented under the suitable choices of regularization parameters and the a priori assumption on the bounds of the exact solution.Numerical implementation is considered and the numerical examples show that our proposed method is effective and stable.
文摘In this paper, we consider the reconstruction of the wave field in a bounded domain. By choosing a special family of functions, the Cauchy problem can be transformed into a Fourier moment problem. This problem is ill-posed. We propose a regularization method for obtaining an approximate solution to the wave field on the unspecified boundary. We also give the convergence analysis and error estimate of the numerical algorithm. Finally, we present some numerical examples to show the effectiveness of this method.
基金supported by the National Natural Science Foundation of China (No.10871131)the Science and Technology Commission of Shanghai Municipality (No.075105118)+1 种基金the Shanghai Leading Academic Discipline Project (No.S30405)the Fund for E-institutes of Shanghai Universities(No.E03004)
文摘In this paper, we investigate the mixed spectral method using generalized Laguerre functions for exterior problems of fourth order partial differential equations. A mixed spectral scheme is provided for the stream function form of the Navier-Stokes equations outside a disc. Numerical results demonstrate the spectral accuracy in space.
基金the China State Major Key Project for Basic Researches and the Science Fund of the Ministry of Education of China.
文摘There are two cases of the exterior problems of the Helmholtz equation. If λ ≥ 0 the bilinear form is coercive, and if λ < 0 it is the scattering problem. We give a new approach of the infinite element method, which enables us to solve these exterior problems as well as corner problems. A numerical example of the scattering problem is given. [ABSTRACT FROM AUTHOR]
基金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.
文摘Considering the fundamental solution of the Laplace equation as the weight function,we give the iterative format for solving the nonhomogeneous Helmholtz equation with variable coefficients. Furthermore, the iteration method of BEM for solving the equation mentioned above is obtained. The numerical example is given in this paper. Finally, the iteration method of BUM mentioned above is compared with the coupled method of BEM that was presented before then by authors.
基金supported by the Fundamental Research Funds for the Central Universities (Grant No. 2010MS080)the Research Fund for Doctoral Program of Higher Education of China (Grant No. 20070487403)
文摘We apply the fast multipole method (FMM) accelerated boundary element method (BEM) for the three-dimensional (3D) Helmholtz equation, and as a result, large-scale acoustic scattering problems involving 400000 elements are solved efficiently. This is an extension of the fast multipole BEM for two-dimensional (2D) acoustic problems developed by authors recently. Some new improvements are obtained. In this new technique, the improved Burton-Miller formulation is employed to over-come non-uniqueness difficulties in the conventional BEM for exterior acoustic problems. The computational efficiency is further improved by adopting the FMM and the block diagonal preconditioner used in the generalized minimum residual method (GMRES) iterative solver to solve the system matrix equation. Numerical results clearly demonstrate the complete reliability and efficiency of the proposed algorithm. It is potentially useful for solving large-scale engineering acoustic scattering problems.
文摘This paper is concerned with the reconstruction of the radiation wave field in the exterior of a bounded two- or three-dimensional domain from the knowledge of Cauchy data on a part of the boundary of the aforementioned domain. It is described by the Cauchy problem for the Helmholtz equation. By using the Dirichlet-to-Neumann map, this problem is transformed into an operator equation with compact operator. We rigorously justify the asymptotic behaviors of singular values of the compact operator. Then a projection method with regularization is applied to solve the operator equation, and the convergence of the regularization method is discussed. Finally, several numerical examples are presented to illustrate the approach. The results demonstrate that the algorithm is effective.Mathematics subject classification: 35R25, 35R30, 78A40.
文摘From the point of view of energy analysis, the cause that the uniqueness of the boundary integral equation induced from the exterior Helmholtz problem does not hold is investigated in this paper. It is proved that the Sommerfeld's condition at the infinity is changed so that it is suitable not only for the radiative wave but also for the absorptive wave when we use the boundary integral equation to describe the exterior Helmholtz problem. There fore, the total energy of the system is conservative. The mathematical dealings to guarantee the uniqueness are discussed based upon this explanation.
文摘We study the coupled mKdV equation by the dressing method via local Riemann-Hilbert problem. With the help of the Lax pairs, we obtain the matrix Riemann-Hilbert problem with zeros. The explicit solutions for the coupled mKdV equation are derived with the aid of the regularization of the Riemann-Hilbert problem.
基金Project supported by the National Natural Science Foundation of China(No.11074170)the State Key Laboratory Foundation of Shanghai Jiao Tong University(No.MSVMS201105)
文摘This paper describes formulation and implementation of the fast multipole boundary element method (FMBEM) for 2D acoustic problems. The kernel function expansion theory is summarized, and four building blocks of the FMBEM are described in details. They are moment calculation, moment to moment translation, moment to local translation, and local to local translation. A data structure for the quad-tree construction is proposed which can facilitate implementation. An analytical moment expression is derived, which is more accurate, stable, and efficient than direct numerical computation. Numerical examples are presented to demonstrate the accuracy and efficiency of the FMBEM, and radiation of a 2D vibration rail mode is simulated using the FMBEM.
基金funded by the Natural Science Foundation of China under grants 12071401,12171238,12261160361,and 11525103the science and technology innovation Program of Hunan Province 2022RC1191.
文摘We extend the pure source transfer domain decomposition method(PSTDDM)to solve the perfectly matched layer approximation of Helmholtz scattering problems in heterogeneous media.We first propose some new source transfer operators,and then introduce the layer-wise and block-wise PSTDDMs based on these operators.In particular,it is proved that the solution obtained by the layer-wise PSTDDM in R2 coincides with the exact solution to the heterogeneous Helmholtz problem in the computational domain.Second,we propose the iterative layer-wise and blockwise PSTDDMs,which are designed by simply iterating the PSTDDM alternatively over two staggered decompositions of the computational domain.Finally,extensive numerical tests in two and three dimensions show that,as the preconditioner for the GMRES method,the iterative PSTDDMs are more robust and efficient than PSTDDMs for solving heterogeneous Helmholtz problems.
文摘In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variational formulation and well posedness of the coupling method are obtained. The convergence and optimal estimates for the approximation solution are provided.
基金supported by the National Natural Science Foundation of China (11172291)the National Science Foundation for Post-doctoral Scientists of China (2012M510162)the Fundamental Research Funds for the Central Universities (KB2090050024)
文摘This paper presents a novel wideband fast multipole boundary element approach to 3D half-space/planesymmetric acoustic wave problems.The half-space fundamental solution is employed in the boundary integral equations so that the tree structure required in the fast multipole algorithm is constructed for the boundary elements in the real domain only.Moreover,a set of symmetric relations between the multipole expansion coefficients of the real and image domains are derived,and the half-space fundamental solution is modified for the purpose of applying such relations to avoid calculating,translating and saving the multipole/local expansion coefficients of the image domain.The wideband adaptive multilevel fast multipole algorithm associated with the iterative solver GMRES is employed so that the present method is accurate and efficient for both lowand high-frequency acoustic wave problems.As for exterior acoustic problems,the Burton-Miller method is adopted to tackle the fictitious eigenfrequency problem involved in the conventional boundary integral equation method.Details on the implementation of the present method are described,and numerical examples are given to demonstrate its accuracy and efficiency.
基金Supported by Major State Research Development Program of China(2005CB321701)National Natural Science Foundation of China(10971083)
文摘In this paper, we give a hybrid method to numerically solve the inverse open cavity scattering problem for cavity shape, given the scattered solution on the opening of the cavity. This method is a hybrid between an iterative method and an integral equations method for solving the Cauchy problem. The idea of this hybrid method is simple, the operation is easy, and the computation cost is small. Numerical experiments show the feasibility of this method, even for cases with noise.
文摘We establish variational formulation and prove the existence and uniqueness of the three dimen- sional axisymmetric Stokes exterior problem in weighted spaces. Error estimates and convergence for P2 - P0 elements with infinite element methods are also obtained. Numerical experiments are presented to verify the theoretical analysis.
文摘In this paper, we study a class of boundary value problems for conformable fractional differential equations under a new definition. Firstly, by using the monotone iterative technique and the method of coupled upper and lower solution, the sufficient condition for the existence of the boundary value problem is obtained, and the range of the solution is determined. Then the existence and uniqueness of the solution are proved by the proof by contradiction. Finally, a concrete example is given to illustrate the wide applicability of our main results.