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

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.     

THE NONCONFORMING FINITE ELEMENT METHOD FOR SIGNORINI PROBLEM

We present the Crouzeix-Raviart linear nonconforming finite element approximation of the variational inequality resulting from Signorini problem. We show if the displacement field is of H2 regularity, then the convergence rate can be improved from O（h3/4） to quasi-optimal O（h｜log h｜1/4） with respect to the energy norm as that of the continuous linear finite element approximation. If stronger but reasonable regularity is available, the convergence rate can be improved to the optimal O（h） as expected by the linear approximation.     

COUPLING OF FINITE ELEMENT AND BOUNDARY ELEMENT METHODS FOR THE SCATTERING BY PERIODIC CHIRAL STRUCTURES

Consider a time-harmonic electromagnetic plane wave incident on a biperiodic structure in R^3. The periodic structure separates two homogeneous regions. The medium inside the structure is chiral and nonhomogeneous. In this paper, variational formulations coupling finite element methods in the chiral medium with a method of integral equations on the periodic interfaces are studied. The well-posedness of the continuous and discretized problems is established. Uniform convergence for the coupling variational approximations of the model problem is obtained.     

MULTIGRID ALGORITHM FOR THE COUPLING SYSTEM OF NATURAL BOUNDARY ELEMENT METHOD AND FINITE ELEMENT METHOD FOR UNBOUNDED DOMAIN PROBLEMS

In this paper, some V-cycle multigrid algorithms are presented for the coupling system arising from the discretization of the Dirichlet exterior problem by coupling the natural boundary element method and finite element method. The convergence of these multigrid algorithms is obtained even with only one smoothing on all levels. The rate of convergence is found uniformly bounded independent of the number of levels and the mesh sizes of all levels, which indicates that these multigrid algorithms are optimal. Some numerical results are also reported.     

ACCURACY ENHANCEMENT FOR THE SIGNORINI PROBLEM WITH FINITE ELEMENT METHOD

In this paper,we study the accuracy enhancement for the frictionless Signorini problem on a polygonal domain with linear finite elements.Numerical test is given to verify our result.     

A BOUNDARY ELEMENT PROCEDURE FOR THE SIGNORINI PROBLEM IN THREE-DIMENSIONAL ELASTICITY

In this paper, a new numerical method for the Signorini problem in three-dimensional elasticity is presented. The problem is reduced to a boundary variational inequality based on a new representation of the derivative of the doublelayer potential. Furthermore, a boundary element procedure is described for the numerical approximation of its solution and an abstract error estimate is given.     

LEAST-SQUARES MIXED FINITE ELEMENT METHODS FOR NONLINEAR PARABOLIC PROBLEMS

1. IntroductionA large number of physical phenomena are modeled by partial differelltial equations orsystems of parabolic type in an evolutionary or eIliptic type at steady state. It is frequentlythe case that a good approximation of some function of the …     

A BOUNDARY ELEMENT METHOD FOR SIGNORINI PROBLEM IN LINEAR ELASTICITY

A Signorini problem in linear elasticity is reduced to a variational inequality on the boundary, and a boundary element approach is described for the numerical approximation of its solution; an error estimate is also given.     

THE ADAPTIVE IMMERSED INTERFACE FINITE ELEMENT METHOD FOR ELASTICITY INTERFACE PROBLEMS

In this paper, we propose adaptive finite element methods with error control for solving elasticity problems with discontinuous coefficients. The meshes in the methods do not need to fit the interfaces. We establish a residual-based a posteriori error estimate which is $λ$-independent multiplicative constants; the Lamé constant $λ$ steers the incompressibility. The error estimators are then implemented and tested with promising numerical results which will show the competitive behavior of the adaptive algorithm.     

热传导—对流问题的混合有限元分析注记

本文研究流体力学中的热传导－对流问题,分析其混合有限元解的存在性和误差估计,并给出一些混合有限元格式的例子。     

ON THE FINITE ELEMENT APPROXIMATION OF SYSTEMS OF REACTION-DIFFUSION EQUATIONS BY p/hp METHODS

We consider the approximation of systems of reaction-diffusion equations, with the finite element method. The highest derivative in each equation is multiplied by a parameter ε∈ （0, 1], and as ε → 0 the solution of the system will contain boundary layers. We extend the analysis of the corresponding scalar problem from [Melenk, IMA J. Numer. Anal. 17（1997）, pp. 577-601], to construct a finite element scheme which includes elements of size O（εp） near the boundary, where p is the degree of the approximating polynomials. We show that, under the assumption of analytic input data, the method yields exponential rates of convergence, independently of ε, when the error is measured in the energy norm associated with the problem. Numerical computations supporting the theory are also presented, which also show that the method yields robust exponential convergence rates when the error in the maximum norm is used.     

奇异非线性抛物方程的时空有限元方法

时空有限元的思想早期出现在Oden和Nickell＆Sackman等人的论文中,它通过统一时空变量,克服了一般有限元方法对时间作差分离散时引起的时间上的低精度,得到了一种解决时间依赖问题的有效方法．之后,在其基础上又发展起来了流线扩散法和特征流线扩散法．Gurtin于1964年提出了一种变分原理,为人们构造时空有限元提供了一个新的途径．1973年Reed和Hill提出间断Galerkin有限元方法．     

Stokes方程的稳定化间断有限元法

本文对定常的Stokes方程提出了一种新的间断有限元法,通过对通常的间断Galerkin有限元法应用稳定化思想,建立了一个相容的稳定间断有限元格式,对速度和压力的任意分片多项式空间Pl(K),Pm(K)的间断有限元逼近证明了解的存在唯一性,给出了关于速度和压力的L2 范数的最优误差估计．     

SPECTRAL AND SPECTRAL ELEMENT METHODS FOR HIGH ORDER PROBLEMS WITH MIXED BOUNDARY CONDITIONS

In this paper, we investigate numerical methods for high order differential equations. We propose new spectral and spectral element methods for high order problems with mixed inhomogeneous boundary conditions, and prove their spectral accuracy by using the recent results on the Jacobi quasi-orthogonal approximation. Numerical results demonstrate the high accuracy of suggested algorithm, which also works well even for oscillating solutions.     

STABILIZED FINITE ELEMENT METHODS FOR THE REISSNER-MINDLIN PLATE

A new stabilized finite element method which is different from Hughes and Franco's (1988) is presented for the Reissner-Mindlin plate model. The least square mesh-dependent residual form of the shear constitute equation is added to the Partial Projection scheme to enhance the stability. Using piecewise polynomials of order k≥1 for the rotations, of order k+1 for the displacement and of order k-1 for the shear, the kth order error-estimates are obtained. Besides, our computing scheme can be also applied to some lower order elements. All error-estimates are obtained independent of the plate thickness, and the stability parameter is an arbitrary positive constant.     

杆有限元模型的特征值反问题

1引言非均匀杆的轴向自由振动方程为其中L为杆的长度,A(x)为杆的横截面积,E(x)为弹性模量,ρ(x)是杆的密度．杆的主振动的一般形式为u(x,t)=u(x)sin(wt φ), (2)其中u(x)满足d/dx[E(x)A(x)du/dx] λρ(x)A(x)u=0,λ=ω2．(3)     

THE COUPLING OF NATURAL BOUNDARY ELEMENT AND FINITE ELEMENT METHOD FOR 2D HYPERBOLIC EQUATIONS

In this paper, we investigate the coupling of natural boundary element and finite element methods of exterior initial boundary value problems for hyperbolic equations. The governing equation is first discretized in time, leading to a time-step scheme, where an exterior elliptic problem has to be solved in each time step. Second, a circular artificial boundary FR consisting of a circle of radius R is introduced, the original problem in an unbounded domain is transformed into the nonlocal boundary value problem in abounded subdomain. And the natural integral equation and the Poisson integral formula are obtained in the infinite domainΩ2 outside circle of radius R. The coupled variational formulation is given. Only the function itself, not its normal derivative at artificial boundary ΓR, appears in the variational equation, so that the unknown numbers are reducedand the boundary element stiffness matrix has a few different elements. Such a coupled method is superior to the one based on direct boundary element method. This paper discusses finite element discretization for variational problem and its corresponding numerical technique, and the convergence for the numerical solutions. Finally, the numerical example is presented to illustrate feasibility and efficiency of this method.     

抛物型初边值问题的有限元与边界积分耦合的离散化及其误差分析

1．引言边界元方法是近二十几年来迅速发展起来的一类新的偏微分方程的数值方法．它的独特之处是将空间的维数降低一维,从而倍受工程技术人员的青睐,并在工程技术与计算数学领域得到越来越广泛的重视和研究．对椭圆型问题,边界元方法的理论与应用研究已取得丰硕成果;对发展型问题,近年来在理论方面的研究也已取得重要进展［6－11］．但边界元方法难以处理非均质问题,而有限元对各类问题及各种区域具有较好的适应性,将两者结合起来可充分发挥各自的优点．文山提出了一种抛物方程初边值问题的有限元与边界积分的耦合方法,其主要思想是…     

边界层问题的小波—有限元解

本文将小波分析与有限元法结合起来，建立了一种小波－有限元计算格式，并用该算法计算了一个典型的边界层问题，探讨了寻找边界层位置的过程以及计算边界层区的内部解及外部解的步骤。计算结果表明，用该法寻找的边界层位置以及所求得的内部解与真实结果完全符合。