求解二维多区域声波散射问题的边界元方法

对于多散射区域的声波散射问题的外Neumann边值问题,用单层位势来逼近每个散射域上的散射波,再利用位势理论的跳跃关系将问题转换为第二类边界积分方程组的求解问题,然后用Nystrom方法进行了求解.对多个随机散射区域的声波散射问题,数值例子体现了该求解方法的可行性和准确性.     

用Backus-Gilbert方法求解声波散射问题

利用位势理论将散射问题的外边界问题转化为第一类边界积分方程求解,再利用Backus-Gilbert方法给出了二维空间的数值结果,与Tikhonov正则化方法比较,虽然精度稍差一些,但是计算方法和计算机实现比较简单.     

无界区域抛物方程自然边界元方法

本文应用自然边界元方法求解无界区域抛物型初边值问题。首先将控制方程对时间进行离散化,得到关于时间步长离散化的椭圆型问题。通过Fourier展开,导出相应问题的自然积分方程和Poisson积分公式。研究了自然积分算子的性质,并讨论了自然积分方程的数值解法,最后给出数值例子。从而解决了抛物型问题的自然边界归化和自然边界元方法。     

利用正则化方法求解浅海声波散射问题

讨论了浅海中波导的声波散射问题，我们将该问题归结为一个第一类边界积分方程，利用正同化方法求解，并证明了该方法的收敛性，数值例子表明了该方法的简单与有效性。     

反演声波阻尼系数的一个逼近方法

１．引言 考虑在均匀介质中传播的声波,此声波碰到障碍Ｄ发生散射．设 Ｄ Ｒ~３为一单连通区域, Ｄ∈ Ｃ~2．设入射波为平面波ｕ~i（ｘ）＝ ｅｘｐ［ｉｋｘ·α］．其中 ｋ＞ ０为波数, ａ为入射角．记总体场ｕ＝ｕｉ十ｕｓ,ｕｓ满足阻尼边界条件．正散射问题归结为求满足其中。表示单位外法向,称为阻尼系数,（１．３）称为 Ｓｏｍｍｅｒｆｅｌｄ辐射条件．将满足辐射条件的Ｈｅｌｍｈｏｌｔｚ方程的解称为辐射解．散射波ｕｓ具有渐近性质［１］上式在所有方向一致成立．我们称为散射波ｕｓ的远场模式（ｆａｒ ｆｉｌｅｄｐａｔｔｅｒ…     

边界元方法的一些研究进展

本文旨在综述我们小组近二十年来在边界元方法这一领域的一些研究成果,在简要介绍边界元方法的基本思想后,主要介绍了一类非线性界面问题的有限元-边界元耦合方法、求解电磁散射问题的有限元-边界元耦合方法和超奇异积分的一类计算方法.     

利用边界元法求解一类重调和方程

边界元法(BEM)和多重互易法(MRM)相结合求解一类重调和方程.通过重调和基本解序列给出的MRM-方法和BEM, 推导出该类问题的MRM-边界变分方程, 用边界元法求解该变分方程, 从而得到重调和方程的近似解, 并给出了解的存在唯一性证明.通过数值算例说明了MRM-方法具有收敛速度快、计算精度高, 易编程等优点, 为使用边界元法数值求解重调和方程提供了方法和理论依据.适合于工程中的实际运算.     

调和方程自然边界元Shannon 小波方法

1 引言 调和方程无论是在数学上还是在物理学中都占有重要地位,它有很多不同的物理背景,在力学和物理学中研究的许多问题都可归结为调和方程的边值问题,所以对调和方程进行深入研究有重要意义.余德浩教授在[3]中主要对调和方程在典型域(即单位圆,上半平面)上的情形进行了考虑.特别地,对单位圆的情形给出了刚度矩阵系数的计算公式和调和方程解的存在唯一性.本文采用由冯康教授[1]开创的自然边界元方法和Galerkin小波方法相耦合,对上半平面的调和方程Neumann问题进行了研究,得到十分有效的计算结果.     

圆外区域求解Helmholtz方程

1引言许多科学和工程计算问题都可以归结为无界区域上的偏微分方程边值问题．而求解椭圆方程边值问题的常用技术是有限元方法,可是对于无界区域,在用有限元方法求解时,往往遇到困难．最简单的办法显然是直接略去区域的无界部分求解,但这样做或者导致过低的计算精度,或者要付出很高的计算代价．边界归化,即将求解偏微分方程边值问题转化为边界积分方程,是求解某些无界区域问题的强有力的手段．自70年代以来,有限元和     

带滑动边界条件的Stokes方程的边界元逼近

前言 带滑动边界条件的Stokes方程,在诸如具有自由表面或具有大攻角的流体模型中起着重要作用。在电镀或容器壁可与流体起化学反应等流动问题中,经典的Stokes问题的不滑动边界条件不再成立,而滑动边界条件才是适当的物理模型。关于这一类实际问题,已有一些数值结果,但仅有很少的工作是就一般问题进行系统的分析。在文献[9]中,R.Verfuth就带滑动边界条件的定常Navier-Stokes方程给出了一种混合有限     

A symmetric boundary element method for the Stokes problem in multiple connected domains

We consider a symmetric Galerkin boundary element method for the Stokes problem with general boundary conditions including slip conditions. The boundary value problem is reformulated as Steklov–Poincaré boundary integral equation which is then solved by a standard approximation scheme. An essential tool in our approach is the invertibility of the single layer potential which requires the definition of appropriate factor spaces due to the topology of the domain. Here we describe a modified boundary element approach to solve Dirichlet boundary value problems in multiple connected domains. A suitable extension of the standard single layer potential leads to an operator which is elliptic on the original function space. Copyright © 2003 John Wiley & Sons, Ltd.     

A system of boundary integral equations for the transmission problem in acoustics

In this paper, we reduce the classical two-dimensional transmission problem in acoustic scattering to a system of coupled boundary integral equations (BIEs), and consider the weak formulation of the resulting equations. Uniqueness and existence results for the weak solution of corresponding variational equations are established. In contrast to the coupled system in Costabel and Stephan (1985) [4], we need to take into account exceptional frequencies to obtain the unique solvability. Boundary element methods (BEM) based on both the standard and a two-level fast multipole Galerkin schemes are employed to compute the solution of the variational equation. Numerical results are presented to verify the efficiency and accuracy of the numerical methods.     

A convolution quadrature Galerkin boundary element method for the exterior Neumann problem of the wave equation

The numerical solution of the Neumann problem of the wave equation on unbounded three‐dimensional domains is calculated using the convolution quadrature method for the time discretization and a Galerkin boundary element method for the spatial discretization. The mathematical analysis that has been built up for the Dirichlet problem is extended and developed for the Neumann problem, which is important for many modelling applications. Numerical examples are then presented for one of these applications, modelling transient acoustic radiation. Copyright © 2009 John Wiley & Sons, Ltd.     

A boundary integral equation for the transmission eigenvalue problem for Maxwell equation

We propose a new integral equation formulation to characterize and compute transmission eigenvalues in electromagnetic scattering. As opposed to the approach that was recently developed by Cakoni, Haddar and Meng (2015) which relies on a two‐by‐two system of boundary integral equations, our analysis is based on only one integral equation in terms of the electric‐to‐magnetic boundary trace operator that results in a simplification of the theory and in a considerable reduction of computational costs. We establish Fredholm properties of the integral operators and their analytic dependence on the wave number. Further, we use the numerical algorithm for analytic nonlinear eigenvalue problems that was recently proposed by Beyn (2012) for the numerical computation of the transmission eigenvalues via this new integral equation.     

A boundary integral equation method for the transmission eigenvalue problem

We propose a new integral equation formulation to characterize and compute transmission eigenvalues for constant refractive index that play an important role in inverse scattering problems for penetrable media. As opposed to the recently developed approach by Cossonnière and Haddar [1,2] which relies on a two by two system of boundary integral equations our analysis is based on only one integral equation in terms of Dirichlet-to-Neumann or Robin-to-Dirichlet operators which results in a noticeable reduction of computational costs. We establish Fredholm properties of the integral operators and their analytic dependence on the wave number. Further we employ the numerical algorithm for analytic non-linear eigenvalue problems that was recently proposed by Beyn [3] for the numerical computation of transmission eigenvalues via this new integral equation.     

A degenerate diffusion problem with dynamical boundary conditions Diffusion problem with dynamical boundary conditions

The purpose of this paper is to study the existence, the uniqueness and the limit in , as of solutions of general initial-boundary-value problems of the form and in a bounded domain with dynamical boundary conditions of the form Received: 5 December 2000 / Revised version: 20 November 2001 / Published online: 4 April 2002