声波方程吸收边界条件的稳定性分析

引言对于无界区域中波动现象的数值模拟，必需引进人工边界将计算限制在一个有界区域上．为了确定解，需要在人工边界上加适当的边界条件．对于声波和弹性波方程，这样的一组人工边界条件，也叫吸收边界条件，在［1，2］中被系统地构造出来．对于声波方程，这些吸收边界条件恰好是单程波方程的近似．如山中所指出，减少边界反射，便于在计算中应用和稳定性是构造吸收边界条件的三点关键．Ellgqllist和Maid。用模态分析方法15］证明，带有［IJ中构造的吸收边界条件的波动方程初边值问题是适定的，并且估计了人工边界所产生的误差．对于更广…     

抛物型初边值问题的自然积分方程及其数值解法

１．引言数值求解无界区域的偏微分方程,自然的处理方式是削去区域的无界部分,即引入一条适当的人工边界ｒ。,将原问题的求解限制在一个适当的有界区域Ｄ内,这样必须在人工边界上引入适当边界的条件．于是很自然地导致这样一个问题：'是否存在一个人工边界条件,使得在这边界条件下,原问题在区域Ｄ内所求得的数值解与原无界区域的解在Ｄ上的限制是完全一致的？"这里我们的着眼点是寻求与原无界区域问题等价的数学形式,以便于数值求解．因为边界元方法可以将区域内的问题转化到区域的边界上去处理,经典的边界元方法常被应用．七十年代…     

局部扰动半平面上的散射问题

该文讨论半平面上有局部扰动情况下的散射问题．通过位势理论,应用边界积分方程的方法研究了该问题解的存在与唯一性．主要方法是运用对称反射,使该无界区域上的散射问题变成一个有界区域上的散射问题,只是这一有界区域的边界不光滑．通过仔细分析相应的边界积分算子,作者得到了其解的存在与唯一性．     

Webster方程有限元法的超收敛性

<正>1前言由语音生成过程建立的数学模型是一个无界区域上的声压波动方程,引入人工边界可将这个无界区域分成一个内部有界区域和另一个无界区域,再将时间和空间进行变量分离,就可以得到一个如下一维简单的基于内部有界区间的微分方程边值问题     

一个扩散问题的自然边界元法与有限元法组合

本文讨论由Helmholtz方程描述的扩散问题的自然边界元法与有限元法的组合．取一个圆作为公共边界，用Fourier展开建立边界积分方程，将无界区域上的问题化为有界区域上的非局部边值问题．在变分方程中公共边界上的未知量只包含函数本身而不包含其法向导数，从而减少了未知数的数目，并且边界元剐度矩阵只有极少量不同的元素，有利于数值计算．这种组台方法优越于建立在直接边界元法基础上的组合方法．文中证明了变分解的唯一性，数值解的收敛性和误差估计．最后讨论了数值技术并给出一个算倒．     

无界区域上波导问题的有效谱迦略金方法(英文)

旨在求解二维无界区域上的波导问题.通过在人工边界施加精确非反射边界条件(nonreflecting boundary condition,NRBC),可以将无界问题截断为有界问题,然后,给出了在空间上用傅里叶谱迦略金方法、时间上用Newmark格式和Richardson外插技巧的数值格式来求解该有界问题.介绍了一个快速的时间行进格式来处理包含在NRBC中的时间方向的卷积.提出的算法在空间和时间上都可以达到高阶精度.另外,也证明了该截断问题的适定性.     

关于无界区域和非齐次问题的无限元法

<正> 无限相似单元法(简称无限元法)是一种计算有奇点的椭圆型方程的近似解的方法.如果奇点是边界上的一个角点,那末它是有界区域上的无限元法,如果奇点是无穷远点,那末它是无界区域上的无限元法.Babu(s|ˇ)ka 曾经研究过无界区域上的有限元法,由于只有有限个自由度,因此求解区域不能不是有界的,尽管这个区域可以很大.现在,我们用无限元法处理这一类问题.我们证明了在无界区域上各种无限元方法解的存在唯一性.     

Moving mesh method for a reaction-diffusion equation with traveling heat source on unbounded domain

对无界区域上带移动热源的反应扩散方程提出了局部吸收边界条件.移动网格方法对导出的有界区域问题进行了求解.数值例子显示了当热源移动的速度比较慢的时候,方程会在有限时间内发生爆破现象.而当热源移动的速度足够快时,爆破现象不会发生.数值例子验证了新方法的有效性和精确性.     

Dirichlet外问题的概率数值方法

针对定解区域是无界区域的Dirichlet外问题,提出了一种新的有效的概率数值方法,它是从解的随机表达式出发,将无界区域上的问题转化成区域边界上的问题.此时,只要在边界上进行剖分,将问题离散化,然后在无界区域外的有界区域内构作一个辅助球,并且利用布朗运动、漂移布朗运动从球外一点出发,首中球面的位置和时间的分布等,就可以获得Dirichlet外问题的数值解.     

Dirichlet外问题的概率数值方法

针对定解区域是无界区域的Dirichlet外问题,提出了一种新的有效的概率数值方法,它是从解的随机表达式出发,将无界区域上的问题转化成区域边界上的问题.此时,只要在边界上进行剖分,将问题离散化,然后在无界区域外的有界区域内构作一个辅助球,并且利用布朗运动、漂移布朗运动从球外一点出发,首中球面的位置和时间的分布等,就可以获得Dirichlet外问题的数值解.     

A principle of images for absorbing boundary conditions

When one uses high-order finite difference schemes for the wave equation, for instance fourth order schemes, the treatment of boundary conditions poses a real difficulty since one needs several additional equations (for the nodes close to the boundary), while one single scalar boundary condition is available. In the case of perfectly reflecting boundary conditions, namely the homogeneous Neumann or Dirichlet conditions, this difficulty can be overcomed by the use of the well-known image principle, which permits the extension of the equation outside of the domain of calculation by an appropriate symmetrization of the data. We propose in this article a generalization of this principle to the absorbing boundary conditions. Through a symmetrization process, we are led to introduce a damped wave equation with a damping term supported by the boundary. The treatment of the boundary condition is then replaced by the approximation of this new damped wave equation in the whole space. The theoretical justification of our approach is based on new energy estimates for the wave equation (when high-order absorbing boundary conditions are used), and constitutes an alternative to the use of the well-known Kreiss criterion to prove the stability of the associated initial boundary value problems. © 1994 John Wiley & Sons, Inc.     

STABILITY ANALYSIS OF FINITE ELEMENT METHODS FOR THE ACOUSTIC WAVE EQUATION WITH ABSORBINGBOUNDARY CONDITIONS (PART I)

1.IntroductionInthenumericalsimulationofwavepropagationinunboundedorsemi-unboundedmediumitisnecessarytointroduceartificialboundariestoobtainfinitecomputational'regions.Thensomeboundaryconditionshavetobeimposedontheseboundaries,whichshouldeliminatethe...     

Algebraic derivation of discrete absorbing boundary conditions for the wave equation

Summary. We introduce a new algebraic framework to derive discrete absorbing boundary conditions for the wave equation in the one-dimensional case. The idea is to factor directly the discrete wave operator and then use one of the factors as boundary condition. We also analyse the stability of the schemes obtained this way and perform numerical simulations to estimate their practical value. Received June 14, 1997 / Revised version received September 15, 1997     

Absorbing boundary conditions for electromagnetic wave propagation

In this paper, the theoretical perfectly absorbing boundary condition on the boundary of a half-space domain is developed for the Maxwell system by considering the system as a whole instead of considering each component of the electromagnetic fields individually. This boundary condition allows any wave motion generated within the domain to pass through the boundary of the domain without generating any reflections back into the interior. By approximating this theoretical boundary condition a class of local absorbing boundary conditions for the Maxwell system can be constructed. Well-posedness in the sense of Kreiss of the Maxwell system with each of these local absorbing boundary conditions is established, and the reflection coefficients are computed as a plane wave strikes the artificial boundary. Numerical experiments are also provided to show the performance of these local absorbing boundary conditions

     

Mathematical analysis of absorbing boundary conditions for the wave equation: the corner problem

Our goal in this work is to establish the existence and the uniqueness of a smooth solution to what we call in this paper the corner problem, that is to say, the wave equation together with absorbing conditions at two orthogonal boundaries. First we set the existence of a very smooth solution to this initial boundary value problem. Then we show the decay in time of energies of high order--higher than the order of the boundary conditions. This result shows that the corner problem is strongly well-posed in spaces smaller than in the half-plane case. Finally, specific corner conditions are derived to select the smooth solution among less regular solutions. These conditions are required to derive complete numerical schemes.

     

Estudio comparativo de distintos modelos de condiciones de contorno de impedancia en problemas acústicos

In this paper, different implementations of numerical locally reacting boundary conditions are studied for acoustic problems. In this comparative study we analyze two types of equations, the Euler equations and the wave equation. We also analyze both finite-differences time-domain (FDTD) algorithms, and pseudo-spectral time domain (PSTD) numerical schemes. We compare different numerical implementations existing in the literature by means of exhaustive numerical experiments. These numerical experiments allow for the study of the absorbing properties of the different schemes as a function of the frequency and the angle of the incident sound waves. This novel comparative study will help the acoustic engineer in order to choose the proper numerical scheme for his/her simulations.     

Remarks on the Huygens absorbing boundary conditions for electromagnetics

The Huygens absorbing boundary conditions (ABCs) are promising new implementations of operator ABCs. They have certain advantageous features which are lacked in other operator ABCs. Under certain conditions and with the Huygens ABCs, the transmitted wave depends solely upon the second derivative with respect to time or upon the double integral of the incident wave. For such cases and for problems with a Dirichlet boundary condition, the overall reflection is not unique. Two new examples of the Huygens ABCs are given for such cases. For each example and with a FDTD scheme the newly derived reflection is less than that which has been studied by Bérenger [J.-P. Bérenger, On the Huygens absorbing boundary conditions for electromagnetics, J. Comput. Phys. 226 (2007) 354–378].     

Absorbing boundary conditions for the wave equation and parallel computing

Absorbing boundary conditions have been developed for various types of problems to truncate infinite domains in order to perform computations. But absorbing boundary conditions have a second, recent and important application: parallel computing. We show that absorbing boundary conditions are essential for a good performance of the Schwarz waveform relaxation algorithm applied to the wave equation. In turn this application gives the idea of introducing a layer close to the truncation boundary which leads to a new way of optimizing absorbing boundary conditions for truncating domains. We optimize the conditions in the case of straight boundaries and illustrate our analysis with numerical experiments both for truncating domains and the Schwarz waveform relaxation algorithm.

     

Absorbing Interface Conditions for the Simulation of Wave Propagation on Non-Uniform Meshes

We proposed absorbing interface conditions for the simulation of linear wave propagation on non-uniform meshes. Based on the superposition principle of second-order linear wave equations, we decompose the interface condition problem into two subproblems around the interface: for the first one the conventional artificial absorbing boundary conditions is applied, while for the second one, the local analytic solutions can be derived. The proposed interface conditions permit a two-way transmission of low-frequency waves across mesh interfaces which can be supported by both coarse and fine meshes, and perform a one-way absorption of high-frequency waves which can only be supported by fine meshes when they travel from fine mesh regions to coarse ones. Numerical examples are presented to illustrate the efficiency of the proposed absorbing interface conditions.