Scattering of acoustic waves from a point source over an impedance wedge

In this work we study diffraction of a spherical acoustic wave due to a point source, by an impedance wedge In the exterior of the wedge the acoustic pressure satisfies the stationary wave (Helmholtz) equation and classical impedance boundary conditions on two faces of the wedge, as well as Meixner’s condition at the edge and the radiation conditions at infinity. Solution of the boundary value problem is represented by a Weyl type integral and its asymptotic behavior is discussed. On this way, we derive various components in the far field interpreting them accordingly and discussing their physical meaning.     

Helmholtz边界积分方程中奇异积分间接求解方法

提出了间接求解传统Helmholtz边界积分方程CBIE的强奇异积分和自由项系数,以及Burton-Miller边界积分方程BMBIE中的超强奇异积分的特解法。对于声场的内域问题,给出了满足Helmholtz控制方程的特解,间接求出了CBIE中的强奇异积分和自由项系数。对于声场外域对应的BMBIE中的超强奇异积分,按Guiggiani方法计算其柯西主值积分需要进行泰勒级数展开的高阶近似,公式繁复,实施困难。本文给出了满足Helmholtz控制方程和Sommerfeld散射条件的特解,提出了间接求出超强奇异积分的方法。推导了轴对称结构外场问题的强奇异积分中的柯西主值积分表达式,并通过轴对称问题算例证明了本文方法的高效性。数值结果表明,对于内域问题,采用本文特解法的计算结果优于直接求解强奇异积分和自由项系数的结果,且本文的特解法可避免针对具体几何信息计算自由项系数,因而具有更好的适用性。对于外域问题,两者精度相当,但本文的特解法可避免对核函数进行高阶泰勒级数展开,更易于数值实施。     

Improved non-singular local boundary integral equation method

When the source nodes are on the global boundary in the implementation of local boundary integral equation method (LBIEM),singularities in the local boundary integrals need to be treated specially. In the current paper,local integral equations are adopted for the nodes inside the domain trod moving least square approximation (MLSA) for the nodes on the global boundary,thus singularities will not occur in the new al- gorithm.At the same time,approximation errors of boundary integrals are reduced significantly.As applications and numerical tests,Laplace equation and Helmholtz equa- tion problems are considered and excellent numerical results are obtained.Furthermore, when solving the Hehnholtz problems,the modified basis functions with wave solutions are adapted to replace the usually-used monomial basis functions.Numerical results show that this treatment is simple and effective and its application is promising in solutions for the wave propagation problem with high wave number.     

On Green's function for the reduced wave equation in a spherical annular domain with Dirichlet's boundary conditions

Summary In this study Green's function for the reduced wave equation (Helmholtz equation) in a spherical annular domain with Dirichlet's boundary conditions is derived. The convergence of the series solution representing Green's function is then established. Finally it is shown that Green's function for the Dirichlet problem reduces to Green's function for the exterior of a sphere as given by Franz and Etiènne, when the outer radius is moved towards infinity, and when a special position of the coordinate system is chosen.     

伸缩虚拟边界元法解二维Helmholtz外问题

以位势理论为基础，提出了求解Helmholtz外问题的伸缩虚拟边界元法．给出了该方法在全波数域内获得唯一解的严格数学证明，其核心是通过伸缩虚拟边界使对偶内问题的特征频率(本征值)避开与波数重合，从而保证了解的唯一性，同以往前人提出的几种解法途径相比，该法简单得多；通过诸多边界曲线形状和不同边界量的声辐射算例，从计算精度、稳定性以及克服解的非唯一性等方面，对该方法进行了检验．计算结果表明：对远场或近场辐射声压，该方法都具有非常高的效率和精度．     

Fictitious frequency for the exterior Helmholtz equation subject to the mixed-type boundary condition using BEMs

Boundary integral equations and boundary element methods were employed analytically, semi-analytically and numerically to study the occurrence of fictitious frequency for the exterior Helmholtz equations subject to the mixed-type boundary conditions. A semi-infinite rod and a circular radiator of problems were addressed. Degenerate kernel of the fundamental solution and Fourier series for boundary density were utilized in the null-field integral equation to examine the occurrence of fictitious frequency semi-analytically. The BEM was utilized to solve the solution numerically. The CHIEF technique and Burton and Miller method were adopted to suppress the occurrence of the fictitious frequency. It is emphasized that the occurrence of fictitious frequency depend on the adopted method (singular or hypersingular formulation) no matter what the given type of boundary condition for the problem is. The illustrative examples were tested to verify this finding successfully.     

The Helmholtz Equation in a Locally Perturbed Half-Space with Non-Absorbing Boundary

We obtain uniqueness and existence results of an outgoing solution for the Helmholtz equation in a half-space, or in a compact local perturbation of it, with an impedance boundary condition. It is worth noting that these kinds of domains have unbounded boundaries which lead to a non-classical exterior problem. The established radiation condition is somewhat different from the usual Sommerfeld’s one, due to the appearance of surface waves (in the case of a non-absorbing boundary). A half-space Green’s function framework is used to carry out our computations. This is an extended and detailed version of the previous article “The Helmholtz equation with impedance in a half-space,” Duran et al. (CR Acad Sci Paris Ser I 341:561–566, 2005).     

平面Laplace外边值问题

证明平面调和函数的Dirichlet外问题解存在唯一的充要条件,在此基础上,确立外问题的等价边界积分方程,首次给出外域上的极值原理,对第一类Fredholm边界的积分方程的可解性进行了讨论。     

Integral Solution of a Problem with Non-Standard Boundary Conditions in an Enhanced Theory of Bending of Elastic Plates

A modified boundary integral equation method is used to solve a specific type of mixed boundary value problem in an enhanced theory of bending of elastic plates in which the effects of transverse shear deformation and transverse normal strain are taken into account. The problem considered is characterized by the fact that a combination of transverse displacement and bending and twisting moments is prescribed on the curve which bounds the middle surface of the plate. Both interior and exterior problems are formulated and the corresponding existence and uniqueness results derived.     

Invariant imbedding method for wave problems

A new method, which reduces various boundary value problems for a wave equation to initial value problems, is developed. The scalar Helmholtz equation for the one- and three-dimensional cases is considered. The method is extended to the case of nonlinear media. Its applications to the wave equation for different dimensions and various media are described. The source of the wave field may be situated outside or inside the layer occupied by the medium. Governing equations are obtained for cases when one can neglect the backward scattering. The operators that arise are reduced to integral operators. The problem of wave scattering by a weakly rough surface is briefly considered.     

On Green's functions for the reduced wave equation in a circular annular domain with dirichlet, Neumann and radiation type boundary conditions

Summary In this paper Green's functions for the reduced wave equation (Helmholtz equation) in a circular annular domain with the Dirichlet, the radiation, and Neumann boundary conditions are derived. The convergence of the series representing Green's functions is then established. Finally it is shown that these functions reduce to Green's function for the exterior of a circle as given by Franz and Etiènne when the outer radius is moved towards infinity.     

A second-order radiation boundary condition for the shallow water wave equations on two-dimensional unstructured finite element grids

A second-order radiation boundary condition (RBC) is derived for 2D shallow water problems posed in ‘wave equation’ form and is implemented within the Galerkin finite element framework. The RBC is derived by matching the dispersion relation for the interior wave equation with an approximate solution to the exterior problem for outgoing waves. The matching is correct to second order, accounting for curvature of the wave front and the geometry. Implementation is achieved by using the RBC as an evolution equation for the normal gradient on the boundary, coupled through the natural boundary integral of the Galerkin interior problem. The formulation is easily implemented on non-straight, unstructured meshes of simple elements. Test cases show fidelity to solutions obtained on extended meshes and improvement relative to simpler first-order RBCs.     

敷设多孔吸声材料声腔特征值分析的径向积分边界元法

由于Helmholtz方程的基本解是频率的函数,因此传统边界元法在处理声场特征值问题时具有天生的缺陷。本文采用Laplace方程基本解生成积分方程,通过径向积分法将在此过程中产生的域积分项转化为边界积分。此方法克服了传统边界元法系数矩阵对频率的依赖,同时克服了特解积分法对特解的依赖,并通过对表面声导纳的多项式逼近,将敷设多孔吸声材料声腔特征值问题转化为矩阵多项式,从而避免了复杂的非线性求解。通过数值算例验证了算法的有效性。     

EQUIVALENT BOUNDARY INTEGRAL EQUATIONS WITH INDIRECT VARIABLES FOR PLANE ELASTICITY PROBLEMS

IntroductionIthasbeenratheralonghistorythattheBoundaryElementMethod (BEM )isappliedtosolvetheplaneelasticityproblems[1～2 ].However,theEBIE ,whichisequivalenttotheoriginalboundaryvalueproblem ,hasnotbeenfullyappreciatedandsolvedinBEMcommunity .TheconventionalboundaryintegralequationswithindirectvariablesarediscussedthoroughlyanditisshownthatthepreviousresultsarenotEBIE ,i.e .,sometimes,thereexistsnosolutionormorethanonesolutiontothem .Themainkeyliesintheexactformoftheexteriorproblems.I…     

求解Helmholtz方程基于核重构思想的最小二乘配点法

基于核重构思想构造近似函数，将配点法和最小二乘原理相结合对微分方程进行离散， 建立了Helmholtz方程的最小二乘配点格式，并分别研究了Helmholtz方程的波传播问题和 边界层问题. 通过数值算例可以发现，给出的数值计算结果非常接近于精确解，计算精度明显高于SPH 法的数值结果，且随着节点数目的增加，其精确度越来越高，具有良好的收敛性.     

Solutions of the interior and exterior boundary value problems in plane elasticity by using dislocation distribution layer

Based on a previous publication (Savruk, 1981), a dislocation distribution layer method for the solution of interior and exterior boundary value problem (BVP) is studied in more detail. Properties of an integral operator in the resulting integral equation are studied. It is proved theoretically that the tractions applied on the outer boundary should be in equilibrium. In addition, a dislocation distribution layer method for the solution of exterior BVP is also suggested. In the exterior BVP, the tractions applied on the boundary may not be in equilibrium. In the exterior BVP, one must consider the single-valued condition of displacements. The formulation in the exterior BVP is not same as in the interior BVP. In the process of discretization, a technique for balance of the numbers of resulting algebraic equations and unknowns is suggested. Numerical examples prove that the suggested method can give sufficient accurate results.     

Laplace transform method in boundary value problems for the Helmholtz equation

The purpose of this present note is to obtain boundary integral equation formulations of boundary value problems for the two-dimensional Helmholtz equation. The intergral equations are derived using the Laplace transform method.     

A boundary integral equation method for Navier-Stokes equations—application to flow in annulus of eccentric cylinders

A novel Navier-Stokes solver based on the boundary integral equation method is presented. The solver can be used to obtain flow solutions in arbitrary 2D geometries with modest computational effort. The vorticity transport equation is modelled as a modified Helmholtz equation with the wave number dependent on the flow Reynolds number. The non-linear inertial terms partly manifest themselves as volume vorticity sources which are computed iteratively by tracking flow trajectories. The integral equation representations of the Helmholtz equation for vorticity and Poisson equation for streamfunction are solved directly for the unknown vorticity boundary conditions. Rapid computation of the flow and vorticity field in the volume at each iteration level is achieved by precomputing the influence coefficient matrices. The pressure field can be extracted from the converged streamfunction and vorticity fields. The solver is validated by considering flow in a converging channel (Hamel flow). The solver is then applied to flow in the annulus of eccentric cylinders. Results are presented for various Reynolds numbers and compared with the literature.     

A discontinuous Galerkin method with plane waves and Lagrange multipliers for the solution of short wave exterior Helmholtz problems on unstructured meshes

Recently, a discontinuous Galerkin method with plane wave basis functions and Lagrange multiplier degrees of freedom was proposed for the efficient solution of the Helmholtz equation in the mid-frequency regime. This method was fully developed however only for regular meshes, and demonstrated only for interior Helmholtz problems. In this paper, we extend it to irregular meshes and exterior Helmholtz problems in order to expand its scope to practical acoustic scattering problems. We report preliminary results for two-dimensional short wave problems that highlight the superior performance of this discontinuous Galerkin method over the standard finite element method.     

平面弹性外问题的等价的直接变量边界积分方程

用变分法证明平面弹性力学外边值问题的正确提法。在此基础之上，确立外问题的等价的直接变量边界积分方程。对传统的惯用的直接变量边界积分方程进行了深入的讨论，表明它与原边值问题不等价。