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

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

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

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

波动方程单位分解有限元法及收敛性分析

<正>1引言与问题的提出自黄云清和许进超提出基于单位分解技术的重叠型非匹配网格有限元方法以来,这类方法日益引起人们极大的兴趣.这篇文章将这类有限元方法运用到波动方程当中,分别研究了重叠型非匹配网格波动方程半离散和全离散有限元方法,给出了有限元解及其梯度的误差估计.下面先给出模型问题的简单介绍.     

非线性抛物型方程的变网格有限元法

Three linear finite element schemes with moving mesh are established for a class of nonlinear parabolic equations.Optimal L^2-norm and energy norm error estimates are derived under certain conditions.     

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

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

SUPERCONVERGENCE ANALYSIS FOR CUBIC TRIANGULAR ELEMENT OF THE FINITE ELEMENT SUPERCONVERGENCE ANALYSIS FOR CUBIC TRIANGULAR ELEMENT OF THE FINITE ELEMENT

1. IntroductionAlthough we bed proved the superconvergence of quadratic triangular elem6ntsbefore 1985, the superconvergence reseaxch of k(k 2 3)-degree triangulax elemeds Onlyhas a few advances, e.g.,.Lin, Yan and Zhou (see [15]) Prove that the three degreeHerlnite elements possess superconvergence and Walilbin (see [5-7]) obtains a roughresult by using a fine interinr estimation, that is, the placement fUnCtion or its gradestmay have weak superconvergence in the lOcal syUUnetric points,e…     

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

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

一类波动方程有限元投影格式的误差分析

1.引 言 本文考虑如下不含阻尼项的波动方程的有限元逼近: 其中区域Ω Rd（d=2,3）是足够光滑的有界多边形区域,其边界为Γ=　Ω.初始条件为:当t=0时,u=u0,ut=u1.     

高波数问题的超收敛性

 本文考虑求解Helmholtz方程的有限元方法的超逼近性质以及基于PPR后处理方法的超收敛性质.我们首先给出了矩形网格上的p-次元在收敛条件k（kh）^2p+1≤C₀下的有限元解和基于Lobatto点的有限元插值之间的超逼近以及重构的有限元梯度和精确解之间的超收敛分析.然后我们给出了四边形网格上的线性有限元方法的分析.这些估计都给出了与波数k和网格尺寸h的依赖关系.同时我们回顾了三角形网格上的线性有限元的超收敛结果.最后我们给出了数值实验并且结合Richardson外推进一步减少了误差.     

UNIFORMLY-STABLE FINITE ELEMENT METHODS FOR DARCY-STOKES-BRINKMAN MODELS

In this paper, we consider 2D and 3D Darcy-Stokes interface problems. These equations are related to Brinkman model that treats both Darcy＇s law and Stokes equations in a single form of PDE but with strongly discontinuous viscosity coefficient and zerothorder term coefficient. We present three different methods to construct uniformly stable finite element approximations. The first two methods are based on the original weak formulations of Darcy-Stokes-Brinkman equations. In the first method we consider the existing Stokes elements. We show that a stable Stokes element is also uniformly stable with respect to the coefficients and the jumps of Darcy-Stokes-Brinkman equations if and only if the discretely divergence-free velocity implies almost everywhere divergence-free one. In the second method we construct uniformly stable elements by modifying some well-known H（div）-conforming elements. We give some new 2D and 3D elements in a unified way. In the last method we modify the original weak formulation of Darcy-Stokes- Brinkman equations with a stabilization term. We show that all traditional stable Stokes elements are uniformly stable with respect to the coefficients and their jumps under this new formulation.     

GLOBAL SUPERCONVERGENCE ESTIMATES OF FINITE ELEMENT METHOD FOR SCHRODINGER EQUATION

1.IntroductionConsiderthefollowinginitialboundaryvalueproblemofSchr6dingerequationwheren~[0,1]',at~On/Ot,T>0isaconstant.Theequivalentvariationalformof(1.1)is:foralltE[0,T],findu(t)6Hi(n)satisfiesthefollowingvariationalequation:where(w,v)~IwvdxdenotestheinnerproductofL'(fl)anda(u,v)~(Vu,Vv),ibetheimaginaryunit.Weassumethatthefunctionsarecomplex--valuedandHibertspacesarecomplexspaces.LetThbeaquasiuniformrectangulationoffiwithmeshsizeh>0andS'(O)CHi(fi)bethecorrespondingpiecewisebilinearpol…     

非稳态奇异系数方程的有限元方法

如果求解区域与数据都满足轴对称条件,利用柱坐标变换可将三维Poisson方程-△u=f的第一边值问题化为具有奇异系数的二维稳态问题     

二维抛物型积分微分方程动边界问题的有限元方法

１引言抛物型积微分方程,可广泛用于描述具有记忆的材料的热传导、气体扩散、松散介质中的压力等实际问题中的现象,具有重要研究意义．关于固定空间区域上该类方程的研究,可见文献［１］,［２］;关于动边界抛物型方程,梁国平等已有重要工作［３］,［４］;作者在文［５］中,研究了一维动边界抛物型积微分方程的数值方法．本文研究二维空间区域变动情形下此类方程初边值问题的全离散、半离散有限元逼近格式及有关数值分析．主要特点在于对动边界和时间积分项（Ｖｏｌｔｅｒｒａ项）的处理．对于前者,通过空间变量代换,将问题化为定…     

伪双曲型积分-微分方程的H~1-Galerkin混合元法误差估计

<正>1引言考虑如下一类具有Lipschitz连续边界(?)Ω的凸有界区域Ω上的伪双曲型积分微分方程其中Ω(?)R~d,(d=1,2,3)J=(0,T],对于固定的T,0T∞,函数0a_0≤a(x,t)≤     

Sine-Gordon方程的一类低阶非协调有限元分析

本文讨论了Sine-Gordon方程的一类低阶非协调有限元一般逼近格式,直接利用插值技巧和单元的特殊性质导出了相应未知量的最优误差估计.     

Burgers方程的混合元分析及其数值模拟

1．引言混合有限元法在高阶偏微分方程和含有两个战者两个以上）的未知国数的偏微分方程的数值解的研究中起着重要的作用．但是,到目前为止,混合有限元法主要是用于2n阶或一阶偏微分方程（组）,如二阶椭圆型方程、平面弹性力学方程、双调和方程、Stokes和Navier－stokes方程、抛物型方程以及电磁场方程修见＞到以及当中的参考文献）．然而,R前混合有限元法还没有被用于对非线性的Burgers方程作数值研究．而过去对Burgers方程的数值研究主要采用标准有限元法、差分方法和谱方法修见【IO－12］以及当中的参考文献）．本文的目的是用混…     

A TWO-SCALE HIGHER-ORDER FINITE ELEMENT DISCRETIZATION FOR SCHROEDINGER EQUATION

In this paper, a two-scale higher-order finite element discretization scheme is proposed and analyzed for a Schroedinger equation on tensor product domains. With the scheme, the solution of the eigenvalue problem on a fine grid can be reduced to an eigenvalue problem on a much coarser grid together with some eigenvalue problems on partially fine grids. It is shown theoretically and numerically that the proposed two-scale higher-order scheme not only significantly reduces the number of degrees of freedom but also produces very accurate approximations.     

IMPLICIT-EXPLICIT MULTISTEP FINITE ELEMENT-MIXED FINITE ELEMENT METHODS FOR THE TRANSIENT BEHAVIOR OF A SEMICONDUCTOR DEVICE

The transient behavior of a semiconductor device consists of a Poisson equa-tion for the electric potential and of two nonlinear parabolic equations for the electrondensity and hole density.The electric potential equation is discretized by a mixed finiteelement method. The electron and hole density equations are treated by implicit-explicitmultistep finite element methods. The schemes are very efficient. The optimal order errorestimates both in time and space are derived.     

SUPERCONVERGENCE ANALYSIS OF FINITE ELEMENT METHODS FOR OPTIMAL CONTROL PROBLEMS OF THE STATIONARY BENARD TYPE

In this paper, we consider the finite element approximation of the distributed optimal control problems of the stationary Benard type under the pointwise control constraint. The states and the co-states are approximated by polynomial functions of lowest-order mixed finite element space or piecewise linear functions and the control is approximated by piecewise constant functions. We give the superconvergence analysis for the control; it is proved that the approximation has a second-order rate of convergence. We further give the superconvergence analysis for the states and the co-states. Then we derive error estimates in L^∞-norm and optimal error estimates in L^2-norm.