Acoustic fields with impedance boundary conditions have high engineering applications, such as noise control and evaluation of sound insulation materials, and can be approximated by three-dimensional Helmholtz boundar...Acoustic fields with impedance boundary conditions have high engineering applications, such as noise control and evaluation of sound insulation materials, and can be approximated by three-dimensional Helmholtz boundary value problems. Finite difference method is widely applied to solving these problems due to its ease of use. However, when the wave number is large, the pollution effects are still a major difficulty in obtaining accurate numerical solutions. We develop a fast algorithm for solving three-dimensional Helmholtz boundary problems with large wave numbers. The boundary of computational domain is discrete based on high-order compact difference scheme. Using the properties of the tensor product and the discrete Fourier sine transform method, the original problem is solved by splitting it into independent small tridiagonal subsystems. Numerical examples with impedance boundary conditions are used to verify the feasibility and accuracy of the proposed algorithm. Results demonstrate that the algorithm has a fourth- order convergence in and -norms, and costs less CPU calculation time and random access memory.展开更多
In this paper,Chebyshev interpolation nodes and barycentric Lagrange interpolation basis function are used to deduce the scheme for solving the Helmholtz equation.First of all,the interpolation basis function is appli...In this paper,Chebyshev interpolation nodes and barycentric Lagrange interpolation basis function are used to deduce the scheme for solving the Helmholtz equation.First of all,the interpolation basis function is applied to treat the spatial variables and their partial derivatives,and the collocation method for solving the second order differential equations is established.Secondly,the differential matrix is used to simplify the given differential equations on a given test node.Finally,based on three kinds of test nodes,numerical experiments show that the present scheme can not only calculate the high wave numbers problems,but also calculate the variable wave numbers problems.In addition,the algorithm has the advantages of high calculation accuracy,good numerical stability and less time consuming.展开更多
We study the error analysis of the weak Galerkin finite element method in[24,38](WG-FEM)for the Helmholtz problem with large wave number in two and three dimensions.Using a modified duality argument proposed by Zhu an...We study the error analysis of the weak Galerkin finite element method in[24,38](WG-FEM)for the Helmholtz problem with large wave number in two and three dimensions.Using a modified duality argument proposed by Zhu and Wu,we obtain the pre-asymptotic error estimates of the WG-FEM.In particular,the error estimates with explicit dependence on the wave number k are derived.This shows that the pollution error in the broken H1-norm is bounded by O(k(kh)^(2p))under mesh condition k^(7/2)h^(2)≤C0 or(kh)^(2)+k(kh)^(p+1)≤C_(0),which coincides with the phase error of the finite element method obtained by existent dispersion analyses.Here h is the mesh size,p is the order of the approximation space and C_(0) is a constant independent of k and h.Furthermore,numerical tests are provided to verify the theoretical findings and to illustrate the great capability of the WG-FEM in reducing the pollution effect.展开更多
Conventional finite-difference(FD)methods cannot model acoustic wave propagation beyond Courant-Friedrichs-Lewy(CFL)numbers 0.707 and 0.577 for two-dimensional(2D)and three-dimensional(3D)equal spacing cases,respectiv...Conventional finite-difference(FD)methods cannot model acoustic wave propagation beyond Courant-Friedrichs-Lewy(CFL)numbers 0.707 and 0.577 for two-dimensional(2D)and three-dimensional(3D)equal spacing cases,respectively,thereby limiting time step selection.Based on the definition of temporal and spatial FD operators,we propose a variable-length temporal and spatial operator strategy to model wave propagation beyond those CFL numbers while preserving accuracy.First,to simulate wave propagation beyond the conventional CFL stability limit,the lengths of the temporal operators are modified to exceed the lengths of the spatial operators for high-velocity zones.Second,to preserve the modeling accuracy,the velocity-dependent lengths of the temporal and spatial operators are adaptively varied.The maximum CFL numbers for the proposed method can reach 1.25 and 1.0 in high velocity contrast 2D and 3D simulation examples,respectively.We demonstrate the effectiveness of our method by modeling wave propagation in simple and complex media.展开更多
A weak Galerkin(WG)method is introduced and numerically tested for the Helmholtz equation.This method is flexible by using discontinuous piecewise polynomials and retains the mass conservation property.At the same tim...A weak Galerkin(WG)method is introduced and numerically tested for the Helmholtz equation.This method is flexible by using discontinuous piecewise polynomials and retains the mass conservation property.At the same time,the WG finite element formulation is symmetric and parameter free.Several test scenarios are designed for a numerical investigation on the accuracy,convergence,and robustness of the WG method in both inhomogeneous and homogeneous media over convex and non-convex domains.Challenging problems with high wave numbers are also examined.Our numerical experiments indicate that the weak Galerkin is a finite element technique that is easy to implement,and provides very accurate and robust numerical solutions for the Helmholtz problem with high wave numbers.展开更多
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.展开更多
The scattering of the open cavity filled with the inhomogeneous media is studied.The problem is discretized with a fourth order finite difference scheme and the immersed interfacemethod,resulting in a linear system of...The scattering of the open cavity filled with the inhomogeneous media is studied.The problem is discretized with a fourth order finite difference scheme and the immersed interfacemethod,resulting in a linear system of equations with the high order accurate solutions in the whole computational domain.To solve the system of equations,we design an efficient iterative solver,which is based on the fast Fourier transformation,and provides an ideal preconditioner for Krylov subspace method.Numerical experiments demonstrate the capability of the proposed fast high order iterative solver.展开更多
文摘Acoustic fields with impedance boundary conditions have high engineering applications, such as noise control and evaluation of sound insulation materials, and can be approximated by three-dimensional Helmholtz boundary value problems. Finite difference method is widely applied to solving these problems due to its ease of use. However, when the wave number is large, the pollution effects are still a major difficulty in obtaining accurate numerical solutions. We develop a fast algorithm for solving three-dimensional Helmholtz boundary problems with large wave numbers. The boundary of computational domain is discrete based on high-order compact difference scheme. Using the properties of the tensor product and the discrete Fourier sine transform method, the original problem is solved by splitting it into independent small tridiagonal subsystems. Numerical examples with impedance boundary conditions are used to verify the feasibility and accuracy of the proposed algorithm. Results demonstrate that the algorithm has a fourth- order convergence in and -norms, and costs less CPU calculation time and random access memory.
基金partially supported by National Natural Science Foundation of China(11772165,11961054,11902170)Key Research and Development Program of Ningxia(2018BEE03007)+1 种基金National Natural Science Foundation of Ningxia(2018AAC02003,2020AAC03059)Major Innovation Projects for Building First-class Universities in China’s Western Region(Grant No.ZKZD2017009).
文摘In this paper,Chebyshev interpolation nodes and barycentric Lagrange interpolation basis function are used to deduce the scheme for solving the Helmholtz equation.First of all,the interpolation basis function is applied to treat the spatial variables and their partial derivatives,and the collocation method for solving the second order differential equations is established.Secondly,the differential matrix is used to simplify the given differential equations on a given test node.Finally,based on three kinds of test nodes,numerical experiments show that the present scheme can not only calculate the high wave numbers problems,but also calculate the variable wave numbers problems.In addition,the algorithm has the advantages of high calculation accuracy,good numerical stability and less time consuming.
基金The work was supported in part by the National Natural Science Foundation grants 11471031,91430216,and 11601026NSAF U1530401+1 种基金the U.S.National Science Foundation grant DMS1419040and the China Postdoctoral Science Foundation grant 2016M591053.
文摘We study the error analysis of the weak Galerkin finite element method in[24,38](WG-FEM)for the Helmholtz problem with large wave number in two and three dimensions.Using a modified duality argument proposed by Zhu and Wu,we obtain the pre-asymptotic error estimates of the WG-FEM.In particular,the error estimates with explicit dependence on the wave number k are derived.This shows that the pollution error in the broken H1-norm is bounded by O(k(kh)^(2p))under mesh condition k^(7/2)h^(2)≤C0 or(kh)^(2)+k(kh)^(p+1)≤C_(0),which coincides with the phase error of the finite element method obtained by existent dispersion analyses.Here h is the mesh size,p is the order of the approximation space and C_(0) is a constant independent of k and h.Furthermore,numerical tests are provided to verify the theoretical findings and to illustrate the great capability of the WG-FEM in reducing the pollution effect.
基金the National Natural Science Foundation of China(No.41874144)the Research Foundation of China University of PetroleumBeijing at Karamay(RCYJ2018A-01-001).
文摘Conventional finite-difference(FD)methods cannot model acoustic wave propagation beyond Courant-Friedrichs-Lewy(CFL)numbers 0.707 and 0.577 for two-dimensional(2D)and three-dimensional(3D)equal spacing cases,respectively,thereby limiting time step selection.Based on the definition of temporal and spatial FD operators,we propose a variable-length temporal and spatial operator strategy to model wave propagation beyond those CFL numbers while preserving accuracy.First,to simulate wave propagation beyond the conventional CFL stability limit,the lengths of the temporal operators are modified to exceed the lengths of the spatial operators for high-velocity zones.Second,to preserve the modeling accuracy,the velocity-dependent lengths of the temporal and spatial operators are adaptively varied.The maximum CFL numbers for the proposed method can reach 1.25 and 1.0 in high velocity contrast 2D and 3D simulation examples,respectively.We demonstrate the effectiveness of our method by modeling wave propagation in simple and complex media.
基金supported in part by National Science Foundation Grant DMS-1115097supported in part by National Science Foundation Grants DMS-1016579 and DMS-1318898.
文摘A weak Galerkin(WG)method is introduced and numerically tested for the Helmholtz equation.This method is flexible by using discontinuous piecewise polynomials and retains the mass conservation property.At the same time,the WG finite element formulation is symmetric and parameter free.Several test scenarios are designed for a numerical investigation on the accuracy,convergence,and robustness of the WG method in both inhomogeneous and homogeneous media over convex and non-convex domains.Challenging problems with high wave numbers are also examined.Our numerical experiments indicate that the weak Galerkin is a finite element technique that is easy to implement,and provides very accurate and robust numerical solutions for the Helmholtz problem with high wave numbers.
基金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.
基金The author is grateful for Professor Tao Tang and Dr.Zhonghua Qiao for many helpful and fruitful discussions,and would like to thank Professor Weiwei Sun for constructive suggestions。
文摘The scattering of the open cavity filled with the inhomogeneous media is studied.The problem is discretized with a fourth order finite difference scheme and the immersed interfacemethod,resulting in a linear system of equations with the high order accurate solutions in the whole computational domain.To solve the system of equations,we design an efficient iterative solver,which is based on the fast Fourier transformation,and provides an ideal preconditioner for Krylov subspace method.Numerical experiments demonstrate the capability of the proposed fast high order iterative solver.
