势问题的径向基函数——局部边界积分方程方法

将基于径向基函数构造的具有插值特性的近似函数和局部边界积分方程方法相结合，建立了求解势问题的径向基函数——局部边界积分方程方法，推导了相应离散方程.与其他边界积分方程的无网格方法相比，本文方法具有数值实现过程简单、计算量小、精度高的优点，并可直接施加边界条件.最后通过算例说明了该方法的有效性. 关键词： 径向基函数 无网格方法 局部边界积分方程 势问题     

An FDTD scheme on a face-centered-cubic (FCC) grid for the solution of the wave equation

A method is proposed to improve the numerical dispersion characteristics for simulations of the scalar wave equation in 3D using the FDTD method. The improvements are realized by choosing a face-centered-cubic (FCC) grid instead of the typical Cartesian (Yee) grid, which exhibits non-physical distortions of the wavefront due to the FD stencil. FCC grids are the logical extension of hexagonal grids in 2D, and have been shown previously to provide optimal sampling of space based on close packing of spheres (highest density). The difference equations are developed for the wave equation on this alternative grid, and the dispersion relationship and stability for grids of equal and non-equal aspect ratios are derived. A comparison is made between FCC and Cartesian formulations, based upon having an equal volume density of gridpoints in each method (i.e. the computational storage requirements of each method would be the same for the same simulated space). The comparison shows that the FCC grid exhibits a much more isotropic dispersion relation than the Cartesian grid of equivalent density. Furthermore, for an equivalent density, the FCC method has a more relaxed stability criterion by a factor of approximately 1.35, resulting in a further reduction in computational resources.     

A method for multi-frequency calculation of boundary integral equation in acoustics based on series expansion

The paper presents a method to solve the problem of multi-frequency calculation of Helmholtz boundary integral equation in acoustics. Based on series expansion, system matrices are independent of wavenumber and become the matrix power series of wavenumber. As a result, all matrices in the matrix power series are only dependent on the structure geometry. In addition, an element transform method to calculate the singular integral and Cauchy singular integral is also discussed because the singular integral need to be solved using the method. The convergence of the series expansion method is also proved in this paper. The effectiveness of the method is confirmed by two numerical examples.     

The novel boundary integral equation method for solving arbitrary cross-section waveguides

A novel method for solving arbitrary cross-sectionn waveguides is presented. The novel method is a modification of the eigen-weighted boundary integral equation method; the EWBIEM is modified by using the eigenfunction of a fictitious regular boundary as weighting function, whose eigenvalue may be the known value, and meanwhile using the domain-bases. To confirm the validity of the novel method, numerical analysis are presented for circular groove guide as an example.The Project Supported by National Natural Science Foundation of P.R. China     

Coupling of Dirichlet-to-Neumann boundary condition and finite difference methods in curvilinear coordinates for multiple scattering

The applicability of the Dirichlet-to-Neumann technique coupled with finite difference methods is enhanced by extending it to multiple scattering from obstacles of arbitrary shape. The original boundary value problem (BVP) for the multiple scattering problem is reformulated as an interface BVP. A heterogenous medium with variable physical properties in the vicinity of the obstacles is considered. A rigorous proof of the equivalence between these two problems for smooth interfaces in two and three dimensions for any finite number of obstacles is given. The problem is written in terms of generalized curvilinear coordinates inside the computational region. Then, novel elliptic grids conforming to complex geometrical configurations of several two-dimensional obstacles are constructed and approximations of the scattered field supported by them are obtained. The numerical method developed is validated by comparing the approximate and exact far-field patterns for the scattering from two circular obstacles. In this case, for a second order finite difference scheme, a second order convergence of the numerical solution to the exact solution is easily verified.     

开口/封闭薄壳体声辐射和散射的统一边界积分方程解法

推导了在二维和三维空间下开口和封闭薄壳体在任意阻抗边界条件下声辐射和散射的统一边界积分方程.相对于以前的求解方法,该方程求解声辐射和散射问题具有相同的影响矩阵,能够同时求解薄壳体气动和振动噪声的辐射和散射现象,以及分析壳体声阻抗对声波传播的影响.推导的方程可以应用于叶轮机械、管道等噪声和消声器消声性能的预测等方面.在此方程基础上,可以进一步考虑运动边界和运动介质对声辐射和散射的影响. 关键词： 薄壳体 声阻抗 积分方程 边界元方法     

基于Burton-Miller方程的轴对称结构声学边界元方法

提出了一种求取轴对称结构任意边界条件下声辐射特性的边界元方法。采用Burton和Miller改进型公式将高阶奇异项转化为弱奇异项之和,保证声辐射参数的唯一性,且计算简单精确。将结构表面声压与振速按照旋转轴角度进行Fourier级数展开,利用级数的正交性建立各项待定系数的求解公式;然后转化格林函数的法向偏导为切向偏导,方便直接计算各项积分,并将面积分公式表示为沿结构边界的线积分和沿旋转角度的积分;进一步采用二次等参单元离散结构边界线,建立声压与振速的关系矩阵,从而确定结构声辐射参数。以脉动球源和横向振动球源为例计算,与解析解和传统边界元法结果作对比,说明该方法的有效精确性。     

一种适用于任意阶空间差分时域有限差分方法的色散介质通用吸收边界条件算法

基于伸展坐标完全匹配层方程和辅助微分方程 方法, 给出了一种在时域有限差分(FDTD) 计算中适用于常见色散介质的通用边界条件算法. 该算法适用于任意阶的FDTD空间差分, 并且由于所采用的D-H方程独立于计算区域, 所以可以被用于截断任意电介质. 数值试验结果表明, 与卷积完全匹配层 算法相比较, 所提出的吸收边界条件算法不仅通用性强、 计算复杂度低、 计算时间短, 并且吸收效果有明显的提高.     

On open boundary conditions for long wave equation in the hodograph plane

We discuss implications of the seaward boundary conditions used in initial-boundary value problem formulation of nonlinear shallow-water wave propagation over a linear slope. We first demonstrate the reflection of wave velocity in the case of Dirichlet condition and that of water elevation in the case of Neumann condition. We then show that linear superposition of the two boundary conditions results in much less reflection at the artificial boundary. We also propose a new boundary condition of mixed type and compare its results with that of the aforementioned conditions.     

A basic study on the sound field analysis of periodic structures by boundary integral equation method

This study deals with the development of the approximate method to analyze the sound field around equally spaced finite obstacles, using the periodic boundary condition. First, on the assumption that the equally spaced finite obstacles are the periodically arranged obstacles, the sound field is analyzed by boundary integral equation method with a Green’s function which satisfies the periodic boundary condition. Furthermore, by comparing these results and the exact solution by using the fundamental solution as Green’s function, the validity of the approximate method is also investigated. Next, in order to evaluate the applicability of the approximate method, the simple formula using some parameters, i.e., the frequency, the period, and the number of obstacles, etc., is proposed. The results of the sound field analysis applied the formula are presented.     

Modified periodic-shell boundary condition for dissipative particle dynamics simulations

We have developed the modified periodic-shell boundary condition (BC) for dissipative particle dynamics (DPD) simulations, which enables us to simulate an outer flow problem around an obstacle using a small simulation region. In order to clarify the validity of this BC, we have treated a uniform flow past a circular cylinder. The present BC has been compared with the ordinary BC such as the uniform flow condition. Also, the present results have been compared with those of the numerical results of the Navier–Stokes equation. The ordinary uniform BC is seen to give rise to significantly distorted flow fields and also to significant disappearance of dissipative particles from the simulation region. In contrast, for the present modified periodic-shell BC, the number density of dissipative particles is kept almost constant during a simulation run, and the flow field is in reasonable agreement with the result, which has been obtained by numerical simulations of the Navier–Stokes equation.     

双曲方程外问题的完全无反射人工边界条件

研究了数值求解双曲方程外问题的人工边界方法.在圆形的人工边界上,得到了三类等价的完全无反射的人工边界条件.给出的数值例子验证了方法的有效性.     

Universal quadratures for boundary integral equations on two-dimensional domains with corners

We describe the construction of a collection of quadrature formulae suitable for the efficient discretization of certain boundary integral equations on a very general class of two-dimensional domains with corner points. The resulting quadrature rules allow for the rapid high-accuracy solution of Dirichlet boundary value problems for Laplace’s equation and the Helmholtz equation on such domains under a mild assumption on the boundary data. Our approach can be adapted to other boundary value problems and certain aspects of our scheme generalize to the case of surfaces with singularities in three dimensions. The performance of the quadrature rules is illustrated with several numerical examples.     

An assessment of the high frequency boundary element and Rayleigh integral approximations

This paper revisits the popular Rayleigh integral approximation and also considers a second approximation, the high frequency boundary element method, which is similar to the Rayleigh integral. The Rayleigh integral approximation under consideration is enhanced so that only the elements visible to a particular point in the field are used to calculate the sound pressure at that point. It is demonstrated how both the Rayleigh integral and high frequency boundary element method are approximations to the boundary integral equation so that similarities between the two methods are recognized. Several test cases were conducted in order to assess and compare both methods. The first set of test cases was the pulsating and oscillating sphere. Both methods were then compared on more applied examples including a running engine, construction cab, and transmission housing. It was concluded that though both methods can reliably predict the sound power for some problems, the high frequency boundary element method is the more robust.     

Exact boundary conditions for the initial value problem of convex conservation laws

The initial value problem of convex conservation laws, which includes the famous Burgers’ (inviscid) equation, plays an important rule not only in theoretical analysis for conservation laws, but also in numerical computations for various numerical methods. For example, the initial value problem of the Burgers’ equation is one of the most popular benchmarks in testing various numerical methods. But in all the numerical tests the initial data have to be assumed that they are either periodic or having a compact support, so that periodic boundary conditions at the periodic boundaries or two constant boundary conditions at two far apart spatial artificial boundaries can be used in practical computations. In this paper for the initial value problem with any initial data we propose exact boundary conditions at two spatial artificial boundaries, which contain a finite computational domain, by using the Lax’s exact formulas for the convex conservation laws. The well-posedness of the initial-boundary problem is discussed and the finite difference schemes applied to the artificial boundary problems are described. Numerical tests with the proposed artificial boundary conditions are carried out by using the Lax–Friedrichs monotone difference schemes.     

Exact solutions of the Klein—Gordon equation with ring-shaped oscillator potential by using the Laplace integral transform

We present exact solutions for the Klein-Gordon equation with a ring-shaped oscillator potential. The energy eigenvalues and the normalized wave functions are obtained for a particle in the presence of non-central oscillator potential. The angular functions are expressed in terms of the hypergeometric functions. The radial eigenfunctions have been obtained by using the Laplace integral transform. By means of the Laplace transform method, which is efficient and simple, the radial Klein-Gordon equation is reduced to a first-order differential equation.     

Eigenvalue estimates for the Dirac operator with generalized APS boundary condition

Under two boundary conditions, the generalized Atiyah–Patodi–Singer boundary condition and the modified generalized Atiyah–Patodi–Singer boundary condition, we get the lower bounds for the eigenvalues of the fundamental Dirac operator on compact spin manifolds with nonempty boundary.     

Coupling discontinuous Galerkin methods and retarded potentials for transient wave propagation on unbounded domains

This work deals with the numerical simulation of wave propagation on unbounded domains with localized heterogeneities. To do so, we propose to combine a discretization based on a discontinuous Galerkin method in space and explicit finite differences in time on the regions containing heterogeneities with the retarded potential method to account the unbounded nature of the computational domain. The coupling formula enforces a discrete energy identity ensuring the stability under the usual CFL condition in the interior. Moreover, the scheme allows to use a smaller time step in the interior domain yielding to quasi-optimal discretization parameters for both methods. The aliasing phenomena introduced by the local time stepping are reduced by a post-processing by averaging in time obtaining a stable and second order consistent (in time) coupling algorithm. We compute the numerical rate of convergence of the method for an academic problem. The numerical results show the feasibility of the whole discretization procedure.     

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.