THE ARTIFICIAL BOUNDARY CONDITIONS FOR NUMERICAL SIMULATIONS OF THE COMPLEX AMPLITUDE IN A COUPLED BAY-RIVER SYSTEM

We consider the numerical approximations of the complex amplitude in a coupled bayriver system in this work. One half-circumference is introduced as the artificial boundary in the open sea, and one segment is introduced as the artificial boundary in the river if the river is semi-infinite. On the artificial boundary a sequence of high-order artificial boundary conditions are proposed. Then the original problem is solved in a finite computational domain, which is equivalent to a variational problem. The numerical approximations for the original problem are obtained by solving the variational probiem with the finite element method. The numerical examples show that the artificial boundary conditions given in this work are very effective.     

NUMERICAL SOLUTIONS OF PARABOLIC PROBLEMS ON UNBOUNDED 3-D SPATIAL DOMAIN

In this paper,the numerical solutions of heat equation on 3-D unbounded spatial do-main are considered. n artificial boundary Γ is introduced to finite the computationaldomain.On the artificial boundary Γ,the exact boundary condition and a series of approx-imating boundary conditions are derived,which are called artificial boundary conditions.By the exact or approximating boundary condition on the artificial boundary,the originalproblem is reduced to an initial-boundary value problem on the bounded computationaldomain,which is equivalent or approximating to the original problem.The finite differencemethod and finite element method are used to solve the reduced problems on the finitecomputational domain.The numerical results demonstrate that the method given in thispaper is effective and feasible.     

ARTIFICIAL BOUNDARY METHOD FOR BURGERS'''' EQUATION USING NONLINEAR BOUNDARY CONDITIONS

This paper discusses the numerical solution of Burgers' equation on unbounded do-mains. Two artificial boundaries are introduced and boundary conditions are obtained onthe artificial boundaries, which are in nonlinear forms. Then the original problem is re-duced to an equivalent problem on a bounded domain. Finite difference method is appliedto the reduced problem, and some numerical examples are given to show the effectivenessof the new approach.     

The artificial boundary conditions for incompressible materials on an unbounded domain

Summary. In this paper we consider the numerical simulations of the incompressible materials on an unbounded domain in . A series of artificial boundary conditions at a circular artificial boundary for solving incompressible materials on an unbounded domain is given. Then the original problem is reduced to a problem on a bounded domain, which be solved numerically by a mixed finite element method. The numerical example shows that our artificial boundary conditions are very effective. ReceivedJune 7, 1995 / Revised version received August 19, 1996     

ARTIFICIAL BOUNDARY METHOD FOR THE THREE-DIMENSIONAL EXTERIOR PROBLEM OF ELASTICITY

The exact boundary condition on a spherical artificial boundary is derived for thethree-dimensional exterior problem of linear elasticity in this paper. After this bound-ary condition is imposed on the artificial boundary, a reduced problem only defined in abounded domain is obtained. A series of approximate problems with increasing accuracycan be derived if one truncates the series term in the variational formulation, which isequivalent to the reduced problem. An error estimate is presented to show how the errordepends on the finite element discretization and the accuracy of the approximate problem.In the end, a numerical example is given to demonstrate the performance of the proposedmethod.     

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

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

凹角区域双曲外问题的精确人工边界条件

研究一类凹角区域双曲型外问题的数值方法.先用Newmark方法对时间进行离散化,在每个时间步求解一个椭圆外问题.然后引入人工边界,并获得精确的人工边界条件.给出半离散化问题的变分问题,证明了变分问题的适定性,并给出了误差估计.最后给出数值例子,以示该方法的可行性与有效性.     

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

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

Numerical solution of the heat equation with nonlinear boundary conditions in unbounded domains

The numerical solution of the heat equation in unbounded domains (for a 1D problem‐semi‐infinite line and for a 2D one semi‐infinite strip) is considered. The artificial boundaries are introduced and the exact artificial boundary conditions are derived. The original problems are transformed into problems on finite domains. The space semi‐discretization by finite element method and the full approximation by the implicit‐explicit Euler's method are presented. The solvability of the full discretization schemes is analyzed. Computational examples demonstrate the accuracy and the efficiency of the algorithms. Also, the behavior of blowing up solutions is examined numerically. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 23: 379–399, 2007     

AN ARTIFICIAL BOUNDARY CONDITION FOR THE INCOMPRESSIBLE VISCOUS FLOWS IN A NO-SLIP CHANNEL

1.IntroductionWhencomputingthenumericals0luti0nsofviscousfluidfl0wproblemsinallun-boundedd0main,0neoftenintroducesartificialboundaries,andsetsupanartificialbopundarycondition0nthem;thenthe0riginalproblemisreducedtoaproblemonab0undedc0mputationald0main.InordertoIimitthecomputatio11alcost,theseboundariesmustnotbet00farfromthedomainofinterest.Theref0re,theartificialboundaryc0nditi0nsmustbegoodapprotimationt0the"exact"boundaryconditions(sothatthes0lutionoftheproblemintheboundeddonlainisequaltothes…     

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…     

带吸收边界条件的声波方程显式差分格式的稳定性分析

１．引言 在进行无界或半无界区域上各种波动方程的数值求解时,需引进入工边界以限制计算范围,在这些边界上应加相应的人工边界条件．这种边界条件应保证所求得的有界区域上的解很好地逼近原来无界区域上的解．对波动方程来说,就是在边界上人工反射应尽可能地小,使之对区域内部解的影响在允许的误差范围以内.因而它们被称为无反射边界条件或吸收边界条件．这种边界条件还应保证所形成的有界区域上的微分方程定解问题是适定的．这也是各种数值方法稳定的必要条件。 近二十多年来,已发展了声波方程的各种类型的吸收边界条件,其中以Ｃｌ…     

三维Poisson方程外问题的高阶局部人工边界条件

１引言假设Ｒ３是一分片光滑的闭曲面．是以为边界的无界区域,=Ｒ3是以为边界的有界区域,并且存在球Ｂ0＝ｘｘＲ0我们考虑下面Ｐｏｉｓｓｏｎ方程的外问题：这里f(x）,ｇ（x）是,上的已知函数,f(x)的支集是紧的,即存在一个球面=x·x=R1,使得x＝xｘＲ１,有ｆx＝０．令＝,则ｆ（ｘ）的支集包含在中,令＝xx＝,表示u在上的外法向微商．用流量为零的条件代替无限远处条件（３）,则我们得到一个新的外问题：我们将分别讨论问题（１）－（３）和（４）－（７）的数值解．由于求解区域的无界性,给数值计算带来了本质性的困难．克服此…     

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.     

Convergence of a difference scheme for the heat equation in a long strip by artificial boundary conditions

The numerical solution of the heat equation on a strip in two dimensions is considered. An artificial boundary is introduced to make the computational domain finite. On the artificial boundary, an exact boundary condition is proposed to reduce the original problem to an initial‐boundary value problem in a finite computational domain. A difference scheme is constructed by the method of reduction of order to solve the problem in the finite computational domain. It is proved that the difference scheme is uniquely solvable, unconditionally stable and convergent with the convergence order 2 in space and order 3/2 in time in an energy norm. A numerical example demonstrates the theoretical results.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2007     

Convergence of a finite element scheme for the two-dimensional time-dependent Schrödinger equation in a long strip

This paper addresses the finite element method for the two-dimensional time-dependent Schrödinger equation on an infinite strip by using artificial boundary conditions. We first reduce the original problem into an initial-boundary value problem in a bounded domain by introducing a transparent boundary condition, then fully discretize this reduced problem by applying the Crank-Nicolson scheme in time and a bilinear or quadratic finite element approximation in space. This scheme, by a rigorous analysis, has been proved to be unconditionally stable and convergent, and its convergence order has also been obtained. Finally, two numerical examples are given to verify the accuracy of the scheme.     

Modelling of the forced motions of an elastic beam using the method of integrodifferential relations

A variational approach to the numerical modelling of forced lateral motions of an Euler–Bernoulli elastic beam is developed for a number of linear boundary conditions using the method of integrodifferential relations. A class of linear boundary actions is considered. A family of quadratic functionals, connecting the displacement field of points of the beam with the bending-moment functions in the cross section and the momentum density is proposed. Variational formulations of the original initial-boundary value problem on the motion of the beam are given and the necessary conditions for the functionals introduced to be stationary are analysed. The integral and local quality characteristics of the admissible approximate solutions are determined. The relation between the variational problems, formulated for the beam model, with the classical Hamilton–Ostrogradskii variational principles is demonstrated. An algorithm for constructing approximate systems of ordinary differential equations is developed, the solution of which yields stationary (minimum) values of the functionals introduced on a specified set of displacement fields, moments and momenta. Examples of calculations of the displacements for an elastic beam and an analysis of the quality of the numerical solutions obtained are presented.     

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

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

A NON-OVERLAPPING DOMAIN DECOMPOSITION ALGORITHM BASED ON THE NATURAL BOUNDARY REDUCTION FOR WAVE EQUATIONS IN AN UNBOUNDED DOMAIN

In this paper, a new domain decomposition method based on the natural boundary reduction, which solves wave problems over an unbounded domain, is suggestted. An circular artificial boundary is introduced. The original unbounded domain is divided into two subdomains, an internal bounded region and external unbounded region outside the artificial boundary. A Dirichlet-Neumann(D-N) alternating iteration algorithm is constructed. We prove that the algorithm is equavilent to preconditional Richardson iteration method. Numerical studies are performed by finite element method. The numerical results show that the convergence rate of the discrete D-N iteration is independent of the finite element mesh size.     

解任意四边形板弯曲问题的样条有限元法

关于用样条函数解板的弯曲问题,[1]在1979年讨论了矩形板和菱形板的弯曲;[2]在1981年对简支边界条件的矩形板,用振动梁函数和B样条函数组合作为插值函数,得到了效率更高的算式;[3]在1984年对[2]作了补充,采用拉格朗日乘子法,得到了在各种边界条件下平板弯曲的近似解,但所讨论的仍然是矩形板.