双调和方程混合元的一种新格式

本文介绍了双调和方程混合元的一种新格式,用双二次多项式逼近流函数,双一次多项式逼近涡函数.在拟一致矩形剖分的条件下,证明了此格式具有与C-R格式中分别用双二次多项式逼近相同的收敛阶.     

拟线性Sobolev方程Carey元解的高精度分析

在一种半离散格式下讨论了拟线性Sobolev方程Carey元的超收敛及外推.根据Carey元的构造证明了其有限元解的线性插值与三角形线性元的解相同,再结合线性元的高精度分析和插值后处理技巧导出了超逼近和整体超收敛及后验误差估计.与此同时,根据线性元的误差渐近展开式,构造了一个新的辅助问题,得到了比传统的有限元误差高三阶的外推结果.     

关于线弹性问题混合格式稳定化思想的探讨

首先从混合有限元理论出发,探讨线弹性问题混合变分格式所满足的稳定性条件,从而保证解的存在唯一性.使用连续的分片线性函数和分片常数来分别逼近应力和位移,详细分析了混合格式下稳定化的必要性,有助于更加深入地了解稳定化的基本思想.然后,通过在混合格式中引入位移的跳跃惩罚项,展示了一个无闭锁稳定化混合有限元方法,并证明了此方法是稳定的且是线性收敛的.     

非定常Stokes方程的多级增强Petrov-Galerkin低阶有限元方法

本文对非定常的Stokes方程采用等阶P1/P1元逼近,增强速度有限元空间,应用Petrov-Galerkin途径,根据多级增强空间是否与时间t有关,对时间项采用向后差分,给出了两种全离散的稳定有限元格式.并且证明了这两种格式的稳定性与收敛性.     

对流反应扩散方程的稳定化时间间断时空有限元解的误差估计

在流线迎风Petrov-Galerkin（SUPG）稳定化有限元数值格式的基础上,结合时间方向的变分离散,构造对流反应扩散方程的稳定化时间间断时空有限元格式.该类格式在工程上有一些数值模拟应用,但相关文献没有看到类似数值格式的理论证明.本文以Radau点为节点,构造时间方向的Lagrange插值多项式,证明了稳定化有限元解的稳定性,时间最大模、空间L²（Ω）-模误差估计.文中利用插值多项式和有限元方法相结合的技巧,解耦时空变量,去掉了时空网格的限制条件,提供了时间间断稳定化时空有限元方法的理论证明思路,克服了因时空变量统一导致的实际计算时的复杂性.     

半线性抛物最优控制问题全离散插值系数有限元方法的收敛性分析

本文考虑全离散插值系数有限元方法求解半线性抛物最优控制问题,其中控制变量用分片常数函数逼近,状态变量和对偶状态变量用分片线性函数逼近.对于方程中的半线性项,先用插值系数技巧处理,再用牛顿迭代法求解.通过引入一些辅助变量和投影算子,并利用有限元空间的逼近性质,得到半线性抛物最优控制问题插值系数有限元方法的收敛性结果;数值算例结果验证了理论结果的正确性.     

解两点边值问题的一类修改的三次有限体积元法

构造了求解两点边值问题的一类修改的Lagrange型三次有限体积元法.试探函数空间取以四次Lobatto多项式的零点作为插值节点的Lagrange型三次有限元空间.将插值多项式的导数超收敛点(应力佳点)作为对偶单元的节点,检验函数空间取相应于对偶剖分的分片常数函数空间.证明了新方法具有最优的H1模和L2模收敛阶,讨论了在应力佳点导数的超收敛性,并通过数值实验验证了理论分析结果.     

拟线性双相滞热传导方程的双线性元高精度分析及外推

研究一类拟线性双相滞热传导方程的双线性有限元逼近,利用该元的Ritz投影和插值相结合的技巧,并结合高精度分析和插值后处理技术分别导出了半离散和全离散格式的超逼近和超收敛结果.同时通过构造合适的辅助问题,对半离散格式导出了具有三阶精度的外推解.     

四阶强阻尼非线性波动方程的Hermite型矩形混合有限元分析

讨论了四阶强阻尼非线性波动方程的Hermite型混合有限元方法,并证明了半离散格式下解的存在唯一性.基于该元积分恒等式结果,利用插值与Ritz投影之间的误差估计,可得到半离散格式下O(h~3)阶的超逼近性质,再借助于插值后处理技术导出整体超收敛.进而,通过构造一个新的金离散格式,得到了O(h~3+τ~2)的超逼近和超收敛结果.     

一种推广的对流扩散方程的局部化间断Galerkin方法

研究求解一种产生于径向渗流问题的推广的对流扩散方程的局部化间断Galerkin方法，对一般非线性情形证明了方法的L^2稳定性；对线性情形证明了，当方法取有限元空间为κ次多项式空间时，数值解逼近的L^∞（0，T；L^2)模的误差阶为κ。     

THE DEFECT ITERATION OF THE FINITE ELEMENT FOR ELLIPTIC BOUNDARY VALUE PROBLEMS AND PETROV-GALERKIN APPROXIMATION

1.IntroductionFranketc.of.[l]establishedtheiterateddefectcorrectionschemeforfiniteelemelltofellipticboundaryproblems.FOrlinearellipticboundaryvalueproblem[2--5]havediscllssedtheefficiencyoftheschemebyusillgsuperconvergenceandasymptoticexpansion"lidertheco…     

Approximation of weak adjoints by reverse automatic differentiation of BDF methods

We shed light on the relation between the discrete adjoints of multistep backward differentiation formula (BDF) methods and the solution of the adjoint differential equation. To this end, we develop a functional-analytic framework based on a constrained variational problem and introduce the notion of weak adjoint solutions of ordinary differential equations. We devise a Petrov-Galerkin finite element (FE) interpretation of the BDF method and its discrete adjoint scheme obtained by reverse internal numerical differentiation. We show how the FE approximation of the weak adjoint is computed by the discrete adjoint scheme and prove its convergence in the space of normalized functions of bounded variation. We also show convergence of the discrete adjoints to the classical adjoints on the inner time interval. Finally, we give numerical results for non-adaptive and fully adaptive BDF schemes. The presented framework opens the way to carry over techniques on global error estimation from FE methods to BDF methods.     

EXTRAPOLATION AND A-POSTERIORI ERROR ESTIMATORS OF PETROV-GALERKIN METHODS FOR NON-LINEAR VOLTERRA INTEGRO-DIFFERENTIAL EQUATIONS

1. IntroductionThe pmpme of tabs Paper is to show that the ~ardson edrapolation can be used toenhance the nUmerical solutions generated by a cab of Petrov-Gaierkin lhate element methodsfor the nonlinear VOlterra integrO-chrential equation (VIDE):where j = j(t,y): I x R --+ R and k = k(t,8,g): D x R - R (with D:= {(t,8): 0 S & S t ST}) denote given hmctions.Throughout tab paperl it will always be assumed that the problem (1.1) possesses a piquesolution y E C'(I), namely, the given hmc…     

关于九参数双参数元与广义协调元的对称列式

关于九参数双参数元与广义协调元的对称列式石钟慈（中国科学院计算数学与科学工程计算研究所）石东洋（西安交通大学）ＯＮＴＨＥＳＹＭＭＥＴＲＩＣＡＬＦＯＲＭＯＦＴＨＥ９－ＰＡＲＡＭＥＴＥＲＤＯＵＢＬＥＳＥＴＰＡＲＭＥＴＥＲＥＬＥＭＥＮＴＡＮＤＴＨＥＧＥＮＥ...     

An h-p Petrov-Galerkin finite element method for linear Volterra integro-differential equations

We analyze an h-p version Petrov-Galerkin finite element method for linear Volterra integrodifferential equations. We prove optimal a priori error bounds in the L ²- and H ¹-norm that are explicit in the time steps, the approximation orders and in the regularity of the exact solution. Numerical experiments confirm the theoretical results. Moreover, we observe that the numerical scheme superconverges at the nodal points of the time partition.     

Galerkin finite element method and error analysis for the fractional cable equation

The cable equation is one of the most fundamental equations for modeling neuronal dynamics. These equations can be derived from the Nernst-Planck equation for electro-diffusion in smooth homogeneous cylinders. Fractional cable equations are introduced to model electrotonic properties of spiny neuronal dendrites. In this paper, a Galerkin finite element method(GFEM) is presented for the numerical simulation of the fractional cable equation(FCE) involving two integro-differential operators. The proposed method is based on a semi-discrete finite difference approximation in time and Galerkin finite element method in space. We prove that the numerical solution converges to the exact solution with order O(τ+h^l+1) for the lth-order finite element method. Further, a novel Galerkin finite element approximation for improving the order of convergence is also proposed. Finally, some numerical results are given to demonstrate the theoretical analysis. The results show that the numerical solution obtained by the improved Galerkin finite element approximation converges to the exact solution with order O(τ²+h^l+1).     

GLOBAL EXISTENCE, UNIFORM DECAY AND EXPONENTIAL GROWTH FOR A CLASS OF SEMI-LINEAR WAVE EQUATION WITH STRONG DAMPING

In this paper, we consider the nonlinearly damped semi-linear wave equation associated with initial and Dirichlet boundary conditions. We prove the existence of a local weak solution and introduce a family of potential wells and discuss the invariants and vacuum isolating behavior of solutions. Furthermore, we prove the global existence of solutions in both cases which are polynomial and exponential decay in the energy space respectively, and the asymptotic behavior of solutions for the cases of potential well family with 0相似文献   

Stability of semilinear elliptic optimal control problems with pointwise state constraints

A semi-linear elliptic control problems with distributed control and pointwise inequality constraints on the control and the state is considered. The general optimization problem is perturbed by a certain class of perturbations, and we establish convergence of local solutions of the perturbed problems to a local solution of the unperturbed optimal control problem. This class of perturbations include finite element discretization as well as data perturbation such that the theory implies convergence of finite element approximation and stability w.r.t.?noisy data.     

Some remarks on the Balog–Wooley decomposition theorem and quantities <Emphasis Type="Italic">D</Emphasis><Superscript>+</Superscript>, <Emphasis Type="Italic">D</Emphasis><Superscript>×</Superscript>

A finite element with new properties of approximation of higher derivatives is constructed, and a method for the construction of a finite element space in the planar case is proposed. The method is based on Yu.N. Subbotin’s earlier results as well as on the results obtained in this paper. The constructed piecewise polynomial function possesses the continuity property and new approximation properties.