The coupling of boundary integral and finite element methods for the bidimensional exterior steady stokes problem

In this paper, we represent a new numerical method for solving the steady-state Stokes equations in an unbounded plane domain. The technique consists in coupling the boundary integral and the finite element methods. An artificial smooth boundary is introduced separating an interior inhomogeneous region from an exterior one. The solution in the exterior domain is represented by an integral equation over the artificial boundary. This integral equation is incorporated into a velocitypressure formulation for the interior region, and a finite element method is used to approximate the resulting variational problem. This is studied by means of an abstract framework, well adapted to the model problem, in which convergence results and optimal error estimates are derived. Computer results will be discussed in a forthcoming paper.     

Error estimates for the finite element approximation of linear elastic equations in an unbounded domain

In this paper we present error estimates for the finite element approximation of linear elastic equations in an unbounded domain. The finite element approximation is formulated on a bounded computational domain using a nonlocal approximate artificial boundary condition or a local one. In fact there are a family of nonlocal approximate boundary conditions with increasing accuracy (and computational cost) and a family of local ones for a given artificial boundary. Our error estimates show how the errors of the finite element approximations depend on the mesh size, the terms used in the approximate artificial boundary condition, and the location of the artificial boundary. A numerical example for Navier equations outside a circle in the plane is presented. Numerical results demonstrate the performance of our error estimates.

     

The boundary integral method for the linearized rotating Navier-Stokes equations in exterior domain

In this paper, we apply the boundary integral method to the linearized rotating Navier-Stokes equations in exterior domain. Introducing some open ball which decomposes the exterior domain into a finite domain and an infinite domain, we obtain a coupled problem by the linearized rotating Navier-Stokes equations in finite domain and a boundary integral equation without using the artificial boundary condition. For the coupled problem, we show the existence and uniqueness of solution. Finally, we study the finite element approximation for the coupled problem and obtain the error estimate between the solution of the coupled problem and its approximation solution.     

The coupling of boundary integral and finite element methods for the nonstationary exterior flow

In this article, we represent a new numerical method for solving the nonstationary Navier–Stokes equations in an unbounded domain. The technique consists of coupling the boundary integral and the finite element method. The variational formulation and the well-posedness of the coupling method are obtained. The convergence and optimal error estimates for the approximate solution are provided. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 549–565, 1998     

定常不可压Navier—Stokes方程迎风格式

本文利用齐次定解条件对定常不可压Navier—Stokes方程的非线性项进行处理,给出了相应的一种迎风Galerkin有限元算法;针对这种迎风Galerkin有限元算法,在迎风参数满足一定条件下,利用其三项式具有的一些很好性质,更简单地证明了该问题解的存在唯一性。     

三维热传导型半导体问题的交替方向特征有限元方法及理论分析

１．引言 三维热传导型半导体器件瞬态问题的数学模型由四个非线性偏微分方程描述 ［１,２］．工程研究中一般考虑绝流边条件,由于绝流条件可以看作一反射条件来处理、为了数值分析方便,我们在此考虑三维周期问题： 其中,　＝［０,１］３,未知函数是电子位势　;电子,空穴浓度ｅ,ｐ;温度函数Ｔ．方程（１,１）－（１.４）中出现的系数均有正的上下界,且是　周期的． ａ＝Ｑ／ε,Ｑ,ε分别表示电子负荷和介电系数,均为正常数．Ｎ（ｘ）是给定的函数．Ｄｓ（ｘ）为扩散系数,μｓ（ｘ）为迁移率,ｓ＝ｅ,Ｐ．Ｒ（ｅ,ｐ,Ｔ）…     

The coupling of bem and fem for the two-dimensional viscous flow problem

In this paper we propose a numerical scheme for treating the problem of sJow viscous flow past an obstacle in the plane. This scheme is a combination of boundary element and finite element methods. By introducing an auxiliary boundary curve, we divide the region under consideration into two subregions, an inner and an outer region. In the inner region, we employ a finite element method (FEM) for solving a system of simplified field equations with proper natural boundary conditions. In the outer region, the solution is expressed in the form of a simple-layer potential with density function satisfying a system of modified integral equations of the first kind. The latter are solved by a boundary element method (BEM). Both solutions are matched on the common auxiliary boundary curve. Error estimates in suitable function spaces are derived in terms of the mesh widths as well as the small parameters, the Reynolds numbers     

On the numerical analysis of finite element and dirichlet-to-neumann methods for nonlinear exterior transmission problems

We provide the numerical analysis of the combination of finite elements and Dirichlet-to-Neumann mappings (based on boundary integral operators) for a class of nonlinear exterior transmission problems whose weak formulations reduce to Lipschitz-continuous and strongly monotone operator equations. As a model we consider a nonlinear second order elliptic equation in divergence form in a bounded inner region of the plane, coupled with the Laplace equation in the corresponding unbounded exterior part. A discrete Galerkin scheme is presented by using linear finite elements on a triangulation of the domain, and then applying numerical quadrature and analytical formulae to evaluate all the linear, bilinear and semilinear forms involved. We prove the unique solvability of the discrete equations, and show the strong convergence of the approximate solutions. Furthermore, assuming additional regularity on the solution of the continuous operator equation, the asymptotic rate of convergence O(h) is also derived. Finally, numerical experiments are presented, which confirm the convergence results.     

A System of Plane Elasticity Canonical Integral Equations and Its Application

In this paper, we obtain a new system of canonical integral equations for the plane elasticity problem over an exterior circular domain, and give its numerical solution. Coupling with the classical finite element method, it can be used for solving general plane elasticity exterior boundary value problems. This system of highly singular equations is also an exact boundary condition on the artificial boundary. It can be approximated by a series of nonsingular integral boundary conditions.     

The coupling of boundary integral and finite element methods for nonmonotone nonlinear problems ∗

In this paper, we apply the coupling of the boundary integral and finite element methods to study the weak solvability of certain nonmonotone nonlinear exterior boundary value problems. In order to convert the original exterior problem into an equivalent nonlocal boundary value problem on a finite region, we employ two different approaches based on the use of one and two integral equations on the coupling boundary. Existence of a solution for the associated weak formulation, and convergence properties of the corresponding Galerkin approximations are deduced from fundamental results in nonlinear functional analysis. Indeed, the main arguments of our proofs are based on a surjectivity theorem for mappings of type (S) and on the Fredholm alternative for nonlinear A-proper mappings.     

The h-p version of the coupling of finite element and boundary element methods for transmission problems in polyhedral domains

Summary. This paper analyzes the rate of convergence of the h-p version of the coupling of the finite element and boundary element method for transmission problems with a linear differential operator with variable coefficients in a bounded polyhedral domain and with constant coefficients in the exterior domain . This procedure uses the variational formulation of the differential equation in and involves integral operators on the interface between and . The finite elements are used to obtain approximate solutions of the differential equation in and the boundary elements are used to obtain approximate solutions of the integral equations. For given piecewise analytic data we show that the Galerkin solution of this coupling procedure converges exponentially fast in the energy norm if the h-p version is used both for finite elements and boundary elements. Received February 10, 1996 / Revised version received April 4, 1997     

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

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

ON THE COUPLING OF BOUNDARY INTEGRAL AND FINITE ELEMENT METHODS FOR SIGNORINI PROBLEMS

1.IntroductionPartialdifferentialequationssubjecttounilateralboundaryconditionsareusuallycalledSignoriniproblemsintheliterature.TheseproblemshavebeenstudiedbymanyauthodssincetheappearenceofthehistoricalpaperbyA.Signoriniin1933[25].Signoriniproblemsaroseinmanyareasofapplicationse.g.,theelasticitywithunilateralconditions[lo],thefluidmechnicsproblemsinmediawithsemipermeableboundaries[8,12],theelectropaintprocess[1]etc.Fortheexistence,uniquenessandregularityresultsforSignorinitypeproblemswerefer…     

On a Class of Coupled Variational Formulations for a Nonlinear Parabolic Equation

51.Introducti0nSince198O)stheoriesandapplicationsofboundaryelementmethods(BEM)orboundaryintegralmethods(BIM)havemadegreatsuccessesfortheparaboliclnit1alboundaryvalueproblems(seeL1-12j),andtheapproachhasbeenappliedtonumericalsolutionsofinitialboundaryva1ueproblemssuccessfully(seeL1-5j'L8j).Thepropertiesofboundaryelementoperatorshavebeenstudiedbyboundaryintegralmethodsbymanyauthors(see.[4j,L6J'[7j'L12J).Theseresultsprovideabasisforconvergencesanderrorestimatesfornumericalapproximationofbou…     

区域分裂法求非线性抛物方程的扩张混合元解

针对非线性抛物方程,给出了全离散的扩张混合元格式,利用一个建立在非重叠型区域分裂技巧上的并行迭代法求解了最后的非线性代数方程组,证明了迭代法的收敛性并给出了最优阶的误差估计.     

On efficient numerical methods for an initial-boundary value problem with nonlocal boundary conditions

Many physical phenomena are modeled by nonclassical hyperbolic boundary value problems with nonlocal boundary conditions. In this paper, the problem of solving the one-dimensional wave equation subject to given initial and non-local boundary conditions is considered. These non-local conditions arise mainly when the data on the boundary cannot be measured directly. Several finite difference methods with low order have been proposed in other papers for the numerical solution of this one dimensional non-classic boundary value problem. Here, we derive a new family of efficient three-level algorithms with higher order to solve the wave equation and also use a Simpson formula with higher order to approximate the integral conditions. Additionally, the fourth-order formula is also adapted to nonlinear equations, in particular to the well-known nonlinear Klein–Gordon equations which many physical problems can be modeled with. Numerical results are presented and are compared with some existing methods showing the efficiency of the new algorithms.     

求解Brinkman方程的修正弱有限元方法

本文针对Brinkman方程引入了一种修正弱Galerkin(MWG)有限元方法.我们通过具有两个离散弱梯度算子的变分形式来逼近模型. 在MWG方法中, 分别用次数为$k$和$k-1$的不连续分段多项式来近似速度函数$u$和压力函数$p$. MWG方法的主要思想是用内部函数的平均值代替边界函数. 因此, 与WG方法相比, MWG方法在不降低准确性的同时, 具有更少的自由度, 对于任意次数不超过$k-1$ 的多项式,MWG方法均可以满足稳定性条件. MWG 方法具有高度的灵活性, 它允许在具有一定形状正则性的任意多边形或多面体上使用不连续函数. 针对$H^1$和$L^22$范数下的速度和压力近似解, 建立了最优阶误差估计. 数值算例表明了该方法的准确性, 收敛性和稳定性.     

Weak Galerkin finite element method for valuation of American options

We introduce a weak Galerkin finite element method for the valuation of American options governed by the Black-Scholes equation. In order to implement, we need to solve the optimal exercise boundary and then introduce an artificial boundary to make the computational domain bounded. For the optimal exercise boundary, which satisfies a nonlinear Volterra integral equation, it is resolved by a higher-order collocation method based on graded meshes. With the computed optimal exercise boundary, the front-fixing technique is employed to transform the free boundary problem to a one- dimensional parabolic problem in a half infinite area. For the other spatial domain boundary, a perfectly matched layer is used to truncate the unbounded domain and carry out the computation. Finally, the resulting initial-boundary value problems are solved by weak Galerkin finite element method, and numerical examples are provided to illustrate the efficiency of the method.     

THE APPROXIMATIONS OF THE EXACT BOUNDARY CONDITION AT AN ARTIFICIAL BOUNDARY FOR LINEARIZED INCOMPRESSIBLE VISCOUS FLOWS

1.IntroductionManyproblemsarisinginfluidmechanicsaregiveninanunboundeddomain,suchasfluidflowaroundobstacles.Whencomputingthenumericalsolutionsoftheseproblems,oneoftenintroducesartificialboundariesandsetsupaxtificialboundaryconditionsonthem.Thentheoriginal…     

Error bounds for the finite element approximation of an incompressible material in an unbounded domain

Summary. In this paper we design high-order local artificial boundary conditions and present error bounds for the finite element approximation of an incompressible elastic material in an unbounded domain. The finite element approximation is formulated in a bounded computational domain using a nonlocal approximate artificial boundary condition or a local one. In fact there are a family of nonlocal approximate artificial boundary conditions with increasing accuracy (and computational cost) and a family of local ones for a given artificial boundary. Our error bounds indicate how the errors of the finite element approximations depend on the mesh size, the terms used in the approximate artificial boundary condition and the location of the artificial boundary. Numerical examples of an incompressible elastic material outside a circle in the plane is presented. Numerical results demonstrate the performance of our error bounds. Received August 31, 1998 / Revised version received November 6, 2001 / Published online March 8, 2002