含开边界二维Stokes问题的Galerkin边界元解法

本文推导了含有开边界的二维有限域上Stokes问题的边界积分方程, 得出基于单层位势的第一类间接边界积分方程.对与之等价的边界变分方程用Galerkin边界元求解以得出单层位势的向量密度. 对于含有开边界端点的边界单元,采用特别的插值函数, 以模拟其固有的奇异性.论文用若干数值算例模拟了含有开边界的有限区域上不可压缩粘性流体的绕流.       

三维Helmholtz方程外问题的自然边界元与有限元耦合法

１．引言 设Γ０是空间闭曲面,Ω是Γ０外部的无界区域,考虑三维Ｈｅｌｍｈｏｌｔｚ 方程外Ｎｅｕｍａｎｎ问题其中 是波数,ｗ是频率,ｃ０为波在均匀介质中的传播速度,ｖ是区域Ω的边界Γ０的外法线方向,即指向由 Γ０包围的内部区域, 　　　　　　　为 Γ０上的已知函数．为了保证问题（１．１）和（１．２）的解的存在唯一性,必须附加上无穷远边界条件,即所谓的 Ｓｏｍｍｅｒｆｅｌｄ辐射条件其中ｉ是虚数单位,　　　　　． 许多数学物理问题,例如时间调和声波对不可穿透的障碍物的散射,海洋水下声波的传播,电磁波的绕射与辐…     

解一类半线性外边值问题的自然边界元与有限元耦合法

In this paper, we combine a finite element approach with the natural boundary element method to stduy the weak solvability and Galerkin approximations of a class of semilinear exterior boundary value problems. Our analysis is mainly based on the variational formulation with constraints. We discuss the error estimate of the finite element solution and obtain the asymptotic rate of convergence O（h^n) Finally, we also give two numerical examples.     

无界区域瞬时涡流问题有限元-边界元耦合的A-φ法的误差分析

针对三维无界区域带有凸多边形导体的瞬时涡流问题,本文提出了一种基于势场的有限元-边界元耦合的方法,从理论上讨论了其能量模误差估计.虽然电场被分解为电矢势A与磁标势φ的梯度之和后增加了方程与未知量的个数,但这种分解可以很好地处理不同介质间的间断.与传统的A-φ法不同,本文讨论了一种全离散的A-φ解耦形式,这样不仅可以避免传统格式所产生的鞍点问题的求解,又可以减少计算量.     

带裂纹方形截面杆扭转问题的自然边界元与有限元耦合法

根据基于区域分解的自然边界元与有限元的耦合法,研究了带裂纹的方形截面杆的扭转问题,编制了耦合法计算程序,计算了几种尺寸截面的抗扭刚度、截面各点的应力及裂纹的应力强度因子,并绘出了裂纹尖端的应力分布图.计算中,还探索了松弛因子对迭代收敛速度的影响.从实践上证实了自然边界元与有限元的耦合法所具有的优点.     

曲边有限元和自然边界元求解外边值问题的耦合法

1　引　　言边界元与有限元耦合法在科学和工程计算中有着独特的作用 .由于区域的无限性给人们常用的有限元方法带来困难 ,边界元方法又难以独立处理非线性和非均质的问题以及具有不规则边界的区域上的问题 ,而两者相结合却可以克服各自的缺点 ,故边界元与有限元耦合法在处理一般区域问题特别是无界区域问题时便得到科学与工程界的青睐 ,获得了比较广泛的应用 .自然边界元方法并不引入新的变量 ,属于直接边界元方法[2 ] [8] .它保持能量不变和原边值问题的许多有用性质 ,例如双线性型的对称性和强制性 ,从而自然积分方程的解的存在唯一性及…     

无穷凹角区域各向异性问题的自然边界元与有限元耦合法

本文研究无穷凹角区域上一类各向异性问题的自然边界元与有限元耦合法.利用自然边界归化原理,获得圆弧或椭圆弧人工边界上的自然积分方程,给出了耦合的变分形式及其数值方法,以及逼近解的收敛性和误差估计,最后给出了数值例子,以示方法的可行性和有效性.     

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

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

弹性动力问题的有限元与特解边界元耦合的数值方法

导出了特解边界元法与有限元法的耦合方程。并应用自由度缩减技术,使耦合方程的自由度缩减到有限元域及其和边界元域的耦合边界上。这样得到的耦合方程不增加原有限元方程的带宽和阶数。耦合方程的求解可以引用求解有限元方程的所有方法,易于程序实现。数值算例结果表明,本文所提出的方法是正确的,是一种较为理想的耦合方法。     

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

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

椭圆外区域上Helmholtz问题的耦合法

以椭圆外区域上Helmholtz方程为例,研究一种带有椭圆人工边界的自然边界元与有限元耦合法,给出了耦合变分问题的适定性及误差分析并给出数值例子.理论分析及数值结果表明,用方法求解椭圆外问题是十分有效的.为求解具有长条型内边界外Helmholtz问题提供了一种很好的数值方法.     

Coupling of FEM and BEM in Shape Optimization

In the present paper we consider the numerical solution of shape optimization problems which arise from shape functionals of integral type over a compact region of the unknown shape, especially L ²-tracking type functionals. The underlying state equation is assumed to satisfy a Poisson equation with Dirichlet boundary conditions. We proof that the shape Hessian is not strictly H ^1/2-coercive at the optimal domain which implies ill-posedness of the optimization problem under consideration. Since the adjoint state depends directly on the state, we propose a coupling of finite element methods (FEM) and boundary element methods (BEM) to realize an efficient first order shape optimization algorithm. FEM is applied in the compact region while the rest is treated by BEM. The coupling of FEM and BEM essentially retains all the structural and computational advantages of treating the free boundary by boundary integral equations.This research has been carried out when the second author stayed at the Department of Mathematics, Utrecht University, The Netherlands, supported by the EU-IHP project Nonlinear Approximation and Adaptivity: Breaking Complexity in Numerical Modelling and Data Representation     

Mixed Problem for an Ultraparabolic Equation in Unbounded Domain

We investigate a mixed problem for a nonlinear ultraparabolic equation in a certain domain Q unbounded in the space variables. This equation degenerates on a part of the lateral surface on which boundary conditions are given. We establish conditions for the existence and uniqueness of a solution of the mixed problem for the ultraparabolic equation; these conditions do not depend on the behavior of the solution at infinity. The problem is investigated in generalized Lebesgue spaces.     

Coupling Canonical Boundary Element Method with FEM to Solve Harmonic Problem over Cracked Domain

Using the canonical boundary reduction, suggested by Feng Kang, coupled with the finite element method, this paper gives the numerical solutions of the harmonic boundary-value problem over the domain with crack or concave angle. When the coupling is conforming, convergence and error estimates are obtained. This coupling removes the limitation of the canonical boundary reduction to some typical domains, and avoids the shortcoming of the classical finite element method, because of which the accuracy is damaged seriously and the approximate solution does not reflect the behaviour of the solution near the singularity. Numerical calculations have verified those conclusions.     

A Mixed Problem for One Pseudoparabolic System in an Unbounded Domain

We prove the existence and uniqueness of a solution of a mixed problem for a system of pseudoparabolic equations in an unbounded (with respect to space variables) domain.     

Coupling of FEM and BEM for a nonlinear interface problem: The h–p version

This article presents some numerical examples for coupling the finite element method (FEM) and the boundary element method (BEM) as analyzed in [11]. This coupling procedure combines the advantages of boundary elements (problems in unbounded regions) and of finite elements (nonlinear problems with inhomogeneous data). In [28], experimental rates of convergence for the h version are presented, where the accuracy of the Galerkin approximation is achieved by refining the mesh. In this article we treat the h–p version, combining an increase of the degree of the piecewise polynomials with a certain mesh refinement. In our model examples, we obtain theoretically and numerically exponential convergence, which indicates a great efficiency in particular if singularities appear. © 1995 John Wiley & Sons, Inc.     

Convergence of the Neumann Series in BEM for the Neumann Problem of the Stokes System

A weak solution of the Neumann problem for the Stokes system in Sobolev space is studied in a bounded Lipschitz domain with connected boundary. A solution is looked for in the form of a hydrodynamical single layer potential. It leads to an integral equation on the boundary of the domain. Necessary and sufficient conditions for the solvability of the problem are given. Moreover, it is shown that we can obtain a solution of this integral equation using the successive approximation method. Then the consequences for the direct boundary integral equation method are treated. A solution of the Neumann problem for the Stokes system is the sum of the hydrodynamical single layer potential corresponding to the boundary condition and the hydrodynamical double layer potential corresponding to the trace of the velocity part of the solution. Using boundary behavior of potentials we get an integral equation on the boundary of the domain where the trace of the velocity part of the solution is unknown. It is shown that we can obtain a solution of this integral equation using the successive approximation method.     

The Density Problem for Unbounded Bergman Operators

The unbounded Bergman operator, the operator of multiplication by on an unbounded open subset of the plane, is considered. We give a complete answer regarding the density problem of unbounded Bergman operators in terms of its equivalence to the problem of bounded point evaluations for the Bergman spaces. Using this equivalence and the notion of Wiener capacity, we obtain simple geometric conditions that classify almost those open subsets of the plane for which the corresponding Bergman operators are densely defined. With the aid of an analytic approach, we are also able to give condition for a large collection of open subsets of the plane for which all the positive integer powers of the corresponding Bergman operators are densely defined. Submitted: December 14, 2001? Revised: January 14, 2001.