二维Helmholtz边界超奇异积分方程解析研究

基于常规边界元法及超奇异边界积分方程复线性耦合的Burton-Miller方法应用于无限域声学问题的最大难点在于处理超奇异积分（二维问题）.目前,此类超奇异积分主要使用各种弱奇异/正则化方法求解,而这些弱奇异/正则化方法具有时间消耗大等弱点.基于围道积分定理,本文给出一种使用常值单元的二维Helmholtz边界超奇异积分的解析表达式.在有限部分积分意义下,所有的奇异和超奇异积分可以解析表达.数值算例表明该解析表达式是有效的.     

A High-Order Discontinuous Galerkin Method for the Two-Dimensional Time-Domain Maxwell's Equations on Curved Mesh

In this paper, a DG (Discontinuous Galerkin) method which has been widely employed in CFD (Computational Fluid Dynamics) is used to solve the two-dimensional time-domain Maxwell's equations for complex geometries on unstructured mesh. The element interfaces on solid boundary are treated in both curved way and straight way. Numerical tests are performed for both benchmark problems and complex cases with varying orders on a series of grids, where the high-order convergence in accuracy can be observed. Both the curved and the straight solid boundary implementation can give accurate RCS (Radar Cross-Section) results with sufficiently small mesh size, but the curved solid boundary implementation can significantly improve the accuracy when using relatively large mesh size. More importantly, this CFD-based high-order DG method for the Maxwell's equations is very suitable for complex geometries.     

Transparent boundary conditions for the Helmholtz equation in some ramified domains with a fractal boundary

The paper addresses a class of boundary value problems in some self-similar ramified domains, with the Laplace or Helmholtz equations. Much stress is placed on transparent boundary conditions which allow the solutions to be computed in subdomains. A self similar finite element method is proposed and tested. It can be used for numerically computing the spectrum of the Laplace operator with Neumann boundary conditions, as well as the eigenmodes. The eigenmodes are normalized by means of a perturbation method and the spectral decomposition of a compactly supported function is carried out. Finally, a numerical method for the wave equation is addressed.     

Near-field performance analysis of locally-conformal perfectly matched absorbers via Monte Carlo simulations

In the numerical solution of some boundary value problems by the finite element method (FEM), the unbounded domain must be truncated by an artificial absorbing boundary or layer to have a bounded computational domain. The perfectly matched layer (PML) approach is based on the truncation of the computational domain by a reflectionless artificial layer which absorbs outgoing waves regardless of their frequency and angle of incidence. In this paper, we present the near-field numerical performance analysis of our new PML approach, which we call as locally-conformal PML, using Monte Carlo simulations. The locally-conformal PML method is an easily implementable conformal PML implementation, to the problem of mesh truncation in the FEM. The most distinguished feature of the method is its simplicity and flexibility to design conformal PMLs over challenging geometries, especially those with curvature discontinuities, in a straightforward way without using artificial absorbers. The method is based on a special complex coordinate transformation which is ‘locally-defined’ for each point inside the PML region. The method can be implemented in an existing FEM software by just replacing the nodal coordinates inside the PML region by their complex counterparts obtained via complex coordinate transformation. We first introduce the analytical derivation of the locally-conformal PML method for the FEM solution of the two-dimensional scalar Helmholtz equation arising in the mathematical modeling of various steady-state (or, time-harmonic) wave phenomena. Then, we carry out its numerical performance analysis by means of some Monte Carlo simulations which consider both the problem of constructing the two-dimensional Green’s function, and some specific cases of electromagnetic scattering.     

Near-field performance analysis of locally-conformal perfectly matched absorbers via Monte Carlo simulations

In the numerical solution of some boundary value problems by the finite element method (FEM), the unbounded domain must be truncated by an artificial absorbing boundary or layer to have a bounded computational domain. The perfectly matched layer (PML) approach is based on the truncation of the computational domain by a reflectionless artificial layer which absorbs outgoing waves regardless of their frequency and angle of incidence. In this paper, we present the near-field numerical performance analysis of our new PML approach, which we call as locally-conformal PML, using Monte Carlo simulations. The locally-conformal PML method is an easily implementable conformal PML implementation, to the problem of mesh truncation in the FEM. The most distinguished feature of the method is its simplicity and flexibility to design conformal PMLs over challenging geometries, especially those with curvature discontinuities, in a straightforward way without using artificial absorbers. The method is based on a special complex coordinate transformation which is ‘locally-defined’ for each point inside the PML region. The method can be implemented in an existing FEM software by just replacing the nodal coordinates inside the PML region by their complex counterparts obtained via complex coordinate transformation. We first introduce the analytical derivation of the locally-conformal PML method for the FEM solution of the two-dimensional scalar Helmholtz equation arising in the mathematical modeling of various steady-state (or, time-harmonic) wave phenomena. Then, we carry out its numerical performance analysis by means of some Monte Carlo simulations which consider both the problem of constructing the two-dimensional Green’s function, and some specific cases of electromagnetic scattering.     

A kernel-free boundary integral method for elliptic boundary value problems

This paper presents a class of kernel-free boundary integral (KFBI) methods for general elliptic boundary value problems (BVPs). The boundary integral equations reformulated from the BVPs are solved iteratively with the GMRES method. During the iteration, the boundary and volume integrals involving Green’s functions are approximated by structured grid-based numerical solutions, which avoids the need to know the analytical expressions of Green’s functions. The KFBI method assumes that the larger regular domain, which embeds the original complex domain, can be easily partitioned into a hierarchy of structured grids so that fast elliptic solvers such as the fast Fourier transform (FFT) based Poisson/Helmholtz solvers or those based on geometric multigrid iterations are applicable. The structured grid-based solutions are obtained with standard finite difference method (FDM) or finite element method (FEM), where the right hand side of the resulting linear system is appropriately modified at irregular grid nodes to recover the formal accuracy of the underlying numerical scheme. Numerical results demonstrating the efficiency and accuracy of the KFBI methods are presented. It is observed that the number of GMRES iterations used by the method for solving isotropic and moderately anisotropic BVPs is independent of the sizes of the grids that are employed to approximate the boundary and volume integrals. With the standard second-order FEMs and FDMs, the KFBI method shows a second-order convergence rate in accuracy for all of the tested Dirichlet/Neumann BVPs when the anisotropy of the diffusion tensor is not too strong.     

A kernel-free boundary integral method for elliptic boundary value problems

This paper presents a class of kernel-free boundary integral (KFBI) methods for general elliptic boundary value problems (BVPs). The boundary integral equations reformulated from the BVPs are solved iteratively with the GMRES method. During the iteration, the boundary and volume integrals involving Green’s functions are approximated by structured grid-based numerical solutions, which avoids the need to know the analytical expressions of Green’s functions. The KFBI method assumes that the larger regular domain, which embeds the original complex domain, can be easily partitioned into a hierarchy of structured grids so that fast elliptic solvers such as the fast Fourier transform (FFT) based Poisson/Helmholtz solvers or those based on geometric multigrid iterations are applicable. The structured grid-based solutions are obtained with standard finite difference method (FDM) or finite element method (FEM), where the right hand side of the resulting linear system is appropriately modified at irregular grid nodes to recover the formal accuracy of the underlying numerical scheme. Numerical results demonstrating the efficiency and accuracy of the KFBI methods are presented. It is observed that the number of GMRES iterations used by the method for solving isotropic and moderately anisotropic BVPs is independent of the sizes of the grids that are employed to approximate the boundary and volume integrals. With the standard second-order FEMs and FDMs, the KFBI method shows a second-order convergence rate in accuracy for all of the tested Dirichlet/Neumann BVPs when the anisotropy of the diffusion tensor is not too strong.     

二维多流体域耦合声场的光滑有限元解法

针对声学有限元分析中四节点等参单元计算精度低,对网格质量敏感的问题,将光滑有限元法引入到多流体域耦合声场的数值分析中,提出了二维多流体域耦合声场的光滑有限元解法。该方法在Helmholtz控制方程与多流体域耦合界面的声压/质点法向速度连续条件的基础上,得到二维多流体耦合声场的离散控制方程,并采用光滑有限元的分区光滑技术将声学梯度矩阵形函数导数的域内积分转换形函数的域边界积分,避免了雅克比矩阵的计算。以管道二维多流体域耦合内声场为数值分析算例,研究结果表明,与标准有限元相比,对单元尺寸较大或扭曲严重的四边形网格模型,光滑有限元的计算精度更高。因此光滑有限元能很好地应用于大尺寸单元或扭曲严重的网格模型下二维多流体域耦合声场的预测,具有良好的工程应用前景。      

Singular meshless method using double layer potentials for exterior acoustics

Time-harmonic exterior acoustic problems are solved by using a singular meshless method in this paper. It is well known that the source points cannot be located on the real boundary, when the method of fundamental solutions (MFS) is used due to the singularity of the adopted kernel functions. Hence, if the source points are right on the boundary the diagonal terms of the influence matrices cannot be derived. Herein we present an approach to obtain the diagonal terms of the influence matrices of the MFS for the numerical treatment of exterior acoustics. By using the regularization technique to regularize the singularity and hypersingularity of the proposed kernel functions, the source points can be located on the real boundary and therefore the diagonal terms of influence matrices are determined. We also maintain the prominent features of the MFS, that it is free from mesh, singularity, and numerical integration. The normal derivative of the fundamental solution of the Helmholtz equation is composed of a two-point function, which is one of the radial basis functions. The solution of the problem is expressed in terms of a double-layer potential representation on the physical boundary based on the potential theory. The solutions of three selected examples are used to compare with the results of the exact solution, conventional MFS, boundary element method, and Dirichlet-to-Neumann finite element method. Good numerical performance is demonstrated by close agreement with other solutions.     

小尺度阻尼边界封闭空间声传递函数的有限元素法计算

如何求解阻尼边界封闭空间中声源点到接收点的低频声传递函数已成为目前小尺度封闭空间可听化技术研究的关键技术，能处理任意形状及复杂边界条件的有限元素法可作为求解该问题的适合方法，以室内声声有源Helmholtz方程及其相应边界方程为基础，本文推导出了用于小尺度阻尼边界封闭空间声传递函数的有限元素求解方法，并编制了相应的计算机程序，在算例中，首先通过与模态叠加法计算结果进行比较，验证了该方法的正确性。最后计算了某型车体内腔中任意两点间声传递函数。     

二阶声场波动方程的非分裂完全匹配层算法

在声场仿真中,完全匹配层(Perfectly Matched Layer,PML)是一种十分有效的吸收边界并得到广泛应用。为了解决基于二阶声场波动方程数值仿真中的吸收边界问题,提出了一种非分裂PML算法。首先,基于伸缩坐标变换,推导了PML算法的频域表达式。然后,通过构造辅助微分方程,得到了非分裂PML的时域表达式。最后,进行了相关理论分析和数值仿真,结果表明:相对于已有的声场分裂PML算法,该算法在保持相同的吸收效率的同时,能较大地节约存储空间,提高计算效率,且更易于实现。      

Helmholtz水声换能器弹性壁液腔谐振频率研究

针对传统Helmholtz水声换能器设计中刚性壁假设的局限性,将Helmholtz腔体的弹性计入到液腔谐振频率计算中,实现低频弹性Helmholtz水声换能器液腔谐振频率精确设计.基于细长圆柱壳腔体的低频集中参数模型,导出了腔体弹性引入的附加声阻抗表达式,得到了弹性壁条件下Helmholtz水声换能器等效电路图,给出了考虑了末端修正的弹性壁Helmholtz共振腔液腔谐振频率计算公式.利用ANSYS软件建立了算例模型,仿真分析了不同材质、半径、长度时的Helmholtz共振腔液腔谐振频率.结果对比表明弹性理论值与仿真值符合得很好,相比起传统的刚性壁理论计算结果,本文的弹性壁理论得出的液腔谐振频率值有所降低,与真实情况更加接近.本文的结论可以为精确设计低频弹性Helmholtz水声换能器提供理论支持.     

Burton–Miller-type singular boundary method for acoustic radiation and scattering

This paper proposes the singular boundary method (SBM) in conjunction with Burton and Miller?s formulation for acoustic radiation and scattering. The SBM is a strong-form collocation boundary discretization technique using the singular fundamental solutions, which is mathematically simple, easy-to-program, meshless and introduces the concept of source intensity factors (SIFs) to eliminate the singularities of the fundamental solutions. Therefore, it avoids singular numerical integrals in the boundary element method (BEM) and circumvents the troublesome placement of the fictitious boundary in the method of fundamental solutions (MFS). In the present method, we derive the SIFs of exterior Helmholtz equation by means of the SIFs of exterior Laplace equation owing to the same order of singularities between the Laplace and Helmholtz fundamental solutions. In conjunction with the Burton–Miller formulation, the SBM enhances the quality of the solution, particularly in the vicinity of the corresponding interior eigenfrequencies. Numerical illustrations demonstrate efficiency and accuracy of the present scheme on some benchmark examples under 2D and 3D unbounded domains in comparison with the analytical solutions, the boundary element solutions and Dirichlet-to-Neumann finite element solutions.     

二维声学多层快速多极子边界元及其应用

与模型自由度的平方成正比的存储量和计算量,使传统边界元无法应用到大型模型的计算。为此,发展了一种二维声学多层快速多极子边界元算法。通过二维Helmholtz核函数展开理论的简要介绍,推导了源点矩计算、源点矩转移、源点矩至本地展开转移、本地展开转移公式,并详细描述了二维声学快速多极子边界元算法的具体实现步骤。使用快速傅里叶插值进行源点矩和本地展开系数的多层传递。采用对角左预处理方法,改善边界方程的条件数,减少迭代求解次数。最后通过数值算例,验证了所发展的二维声学快速多极子算法的正确性和高效性。      

Full time-domain nonlinear coupled dynamic analysis of a truss spar and its mooring/riser system in irregular wave

A new full time-domain nonlinear coupled method has been established and then applied to predict the responses of a Truss Spar in irregular wave.For the coupled analysis,a second-order time-domain approach is developed to calculate the wave forces,and a finite element model based on rod theory is established in three dimensions in a global coordinate system.In numerical implementation,the higher-order boundary element method(HOBEM)is employed to solve the velocity potential,and the 4th-order Adams-Bashforth-Moultn scheme is used to update the second-order wave surface.In deriving convergent solutions,the hull displacements and mooring tensions are kept consistent at the fairlead and the motion equations of platform and mooring-lines/risers are solved simultaneously using Newmark-integration scheme including Newton-Raphson iteration.Both the coupled quasi-static analysis and the coupled dynamic analysis are performed.The numerical simulation results are also compared with the model test results,and they coincide very well as a whole.The slow-drift responses can be clearly observed in the time histories of displacements and mooring tensions.Some important characteristics of the coupled responses are concluded.     

Combined Helmholtz equation-least squares method for reconstructing acoustic radiation from arbitrarily shaped objects

A combined Helmholtz equation-least squares (CHELS) method is developed for reconstructing acoustic radiation from an arbitrary object. This method combines the advantages of both the HELS method and the Helmholtz integral theory based near-field acoustic holography (NAH). As such it allows for reconstruction of the acoustic field radiated from an arbitrary object with relatively few measurements, thus significantly enhancing the reconstruction efficiency. The first step in the CHELS method is to establish the HELS formulations based on a finite number of acoustic pressure measurements taken on or beyond a hypothetical spherical surface that encloses the object under consideration. Next enough field acoustic pressures are generated using the HELS formulations and taken as the input to the Helmholtz integral formulations implemented through the boundary element method (BEM). The acoustic pressure and normal component of the velocity at the discretized nodes on the surface are then determined by solving two matrix equations using singular value decomposition (SVD) and regularization techniques. Also presented are in-depth analyses of the advantages and limitations of the CHELS method. Examples of reconstructing acoustic radiation from separable and nonseparable surfaces are demonstrated.     

A Solver for Helmholtz System Generated by the Discretization of Wave Shape Functions

An interesting discretization method for Helmholtz equations was introduced in B. Després [1]. This method is based on the ultra weak variational formulation (UWVF) and the wave shape functions, which are exact solutions of the governing Helmholtz equation. In this paper we are concerned with fast solver for the system generated by the method in [1]. We propose a new preconditioner for such system, which can be viewed as a combination between a coarse solver and the block diagonal preconditioner introduced in [13]. In our numerical experiments, this preconditioner is applied to solve both two-dimensional and three-dimensional Helmholtz equations, and the numerical results illustrate that the new preconditioner is much more efficient than the original block diagonal preconditioner.     

A Finite Volume Method for the Approximation of Diffusion Operators on Distorted Meshes

A new finite volume method is presented for discretizing general linear or nonlinear elliptic second-order partial-differential equations with mixed boundary conditions. The advantage of this method is that arbitrary distorted meshes can be used without the numerical results being altered. The resulting algorithm has more unknowns than standard methods like finite difference or finite element methods. However, the matrices that need to be inverted are positive definite, so the most powerful linear solvers can be applied. The method has been tested on a few elliptic and parabolic equations, either linear, as in the case of the standard heat diffusion equation, or nonlinear, as in the case of the radiation diffusion equation and the resistive diffusion equation with Hall term.     

Preconditioning based on Calderon’s formulae for periodic fast multipole methods for Helmholtz’ equation

Solution of periodic boundary value problems is of interest in various branches of science and engineering such as optics, electromagnetics and mechanics. In our previous studies we have developed a periodic fast multipole method (FMM) as a fast solver of wave problems in periodic domains. It has been found, however, that the convergence of the iterative solvers for linear equations slows down when the solutions show anomalies related to the periodicity of the problems. In this paper, we propose preconditioning schemes based on Calderon’s formulae to accelerate convergence of iterative solvers in the periodic FMM for Helmholtz’ equations. The proposed preconditioners can be implemented more easily than conventional ones. We present several numerical examples to test the performance of the proposed preconditioners. We show that the effectiveness of these preconditioners is definite even near anomalies.     

Hybrid discrete ordinates—spherical harmonics solution to the Boltzmann Transport Equation for phonons for non-equilibrium heat conduction

The Boltzmann Transport Equation (BTE) for phonons has found prolific use for the prediction of non-equilibrium heat conduction phenomena in semiconductor materials. This article presents a new hybrid formulation and associated numerical procedures for solution of the BTE for phonons. In this formulation, the phonon intensity is first split into two components: ballistic and diffusive. The governing equation for the ballistic component is solved using two different established methods that are appropriate for use in complex geometries, namely the discrete ordinates method (DOM), and the control angle discrete ordinates method (CADOM). The diffusive component, on the other hand, is determined by invoking the first-order spherical harmonics (or P₁) approximation, which results in a Helmholtz equation with Robin boundary conditions. Both governing equations, referred to commonly as the ballistic-diffusive equations (BDE), are solved using the unstructured finite-volume procedure. Results of the hybrid method are compared against benchmark Monte Carlo results, as well as solutions of the BTE using standalone DOM and CADOM for two two-dimensional transient heat conduction problems at various Knudsen numbers. Subsequently, the method is explored for a large-scale three-dimensional geometry in order to assess convergence and computational cost. It is found that the proposed hybrid method is accurate at all Knudsen numbers. From an efficiency standpoint, the hybrid method is found to be superior to direct solution of the BTE both for steady state as well as for unsteady non-equilibrium heat conduction calculations with the computational gains increasing with increase in problem size.