半线性抛物方程的时空有限元方法

讨论半线性抛物方程的连续Galerkin时空有限元方法,利用有限元方法和有限差分方法相结合的技巧,证明了弱解的存在唯一性,给出了时间最大模,空间L~2模,即L~∞(L~2)模误差估计.并给出数值算例证明了连续时空有限元方法对于半线性抛物方程的有效性.     

热传导方程一个新的H~1-Galerkin非协调混合有限元格式

对热传导方程提出了一个新的H~1-Galerkin非协调混合有限元格式,其逼近空间不需满足LBB相容性条件,且在不引进传统的Rutz投影的情况下,得到了与以往协调有限元方法相同的L~2-模和H~1-模的误差估计.     

电报方程H~1-Galerkin非协调混合有限元分析

主要研究一类电报方程的H~1-Galerkin非协调混合有限元方法,在任意四边形网格剖分下,其逼近空间分别取为类Wilson元与双线性Q_1元,在不需要满足LBB相容性条件及不采用传统的Ritz投影的情况下,得到了与常规有限元方法相同的L~2-模和H~1-模的误差估计,进一步拓展了H~1-Galerkin混合有限元和类Wilson元的应用范围.     

多孔介质二相驱动问题一种修正特征有限元方法的迭代解法及其收敛性分析

0 引言 多孔介质二相驱动问题的数学模型是由压力方程与浓度方程组成的偏微分方程组的初边值问题.关于该问题的数值解问题,已有大量的文献.为了得到最优的L~2-模误差估计,好多方法用混合元方法解压力方程.我们知道,混合元法得到的方程组系数矩阵是非正定的,从而解混合元比解标准元要困难得多,虽然许多人研究了混合元方法的求解问题,但到目前为止,还没有看到令人满意的好的算法.为了避开对混合元的求解,著名学者T.F.Russell考虑了用标准有限元方法解压力方程,用特征有限元方法解浓度方程的求解方法及其迭代解法,对只有分子扩散的二相驱动问题得到了最优的L~2模误差估计,对有机械弥散的一般二相驱动问题得不到最优的L~2模误差估计,同时在收敛性证明中要求压力有限元空间的指数至少是二.     

对流扩散问题流线扩散法的最优误差估计

在本文中,我们考虑对流占优扩散问题流线扩散双线性有限元方法。原先的文献在ε≤h~2的条件下,得到了L~2-模最优误差估计,而本文则在ε≤h的条件下得到了相同估计。     

一类三维变动区域上的非线性抛物型方程之全离散有限元形式

1 引言 1986年,L.Cermak和M.Zlamal研究了半导体器件中杂质的重新分布,对具有活动边界的二维非线性扩散问题。给出在时间方向上是一阶精度的全离散有限元格式。证明了格式最优的H~1模和次最优的L~2模估计。1989年.P.Lesaint和R.Touzani对一维变动区域上的热传导方程。经过坐标变换,给出了在固定区域上的全离散有限元格式和最优的L~2模估计。1990年,梁国平和陈志明利用时空有限元,给出了变动区域上线性抛物型的方程的全离散变网格有限元格式。证明了最优的L~2收敛性。本文考虑了一类具有活动边界的三维     

非定常Navier-Stokes方程的质量集中各向异性非协调有限元逼近

该文将一个低阶Crouzeix-Raviart型非协调三角形元应用到非定常Navier-Stokes方程,给出了其质量集中有限元逼近格式.在不需要传统Ritz-Volterra投影下,通过引入两个辅助有限元空间对边界进行估计的技巧,在各向异性网格下导出了速度的L~2模和能量模及压力的L~2模的误差估计.     

页岩气藏渗流模型的特征混合有限元数值模拟

将特征有限元方法和混合有限元方法进行耦合,对页岩气藏渗流模型进行了数值模拟,给出了详细的误差分析,得到了最优的L~2模误差估计,并用数值实验验证了方法的有效性.     

二阶特征值问题的非协调元逼近

本文以非协调三角形线性元为例,讨论了二阶特征值问题的非协调有限元逼近,基于二阶变分问题非协调有限元逼近的有关分析结果,不仅得到了特征值逼近解的误差估计,而且得到了特征函数逼近解的最优的L~2-误差估计和拟最优的L~∞-误差估计。     

二阶双曲方程一个新的非协调混合元超收敛性分析

研究了一类二阶双曲型方程在新混合元格式下的非协调混合有限元方法.在抛弃传统有限元分析的必要工具-Ritz投影算子的前提下,直接利用单元的插值性质,运用高精度分析和对时间t的导数转移技巧,借助于插值后处理技术,分别导出了关于原始变量u的H~1-模和通量=-▽u在L~2-模下的O(h~2)阶超逼近性质和整体超收敛结果.进一步,给出了一些数值算例验证了理论分析的正确性.     

A Low Order Nonconforming Mixed Finite Element Method for Non-Stationary Incompressible Magnetohydrodynamics System

A low order nonconforming mixed finite element method (FEM) is established for the fully coupled non-stationary incompressible magnetohydrodynamics (MHD) problem in a bounded domain in 3D. The lowest order finite elements on tetrahedra or hexahedra are chosen to approximate the pressure, the velocity field and the magnetic field, in which the hydrodynamic unknowns are approximated by inf-sup stable finite element pairs and the magnetic field by $H^1(\Omega)$-conforming finite elements, respectively. The existence and uniqueness of the approximate solutions are shown. Optimal order error estimates of $L^2(H^1)$-norm for the velocity field, $L^2(L^2)$-norm for the pressure and the broken $L^2(H^1)$-norm for the magnetic field are derived.     

非线性抛物型积分微分方程非协调三角形Carey元的收敛性分析

将非协调三角形Carey元应用于二维空间中的非线性抛物型积分微分方程.通过一些新的特殊方法和技巧,给出了有限元解的最优L<'2>模和能量模误差估计.     

具有调和曲率与有限$L^p$曲率的完备流形

Let(M~n, g)(n ≥ 3) be an n-dimensional complete Riemannian manifold with harmonic curvature and positive Yamabe constant. Denote by R and R?m the scalar curvature and the trace-free Riemannian curvature tensor of M, respectively. The main result of this paper states that R?m goes to zero uniformly at infinity if for p ≥ n, the L~p-norm of R?m is finite.As applications, we prove that(M~n, g) is compact if the L~p-norm of R?m is finite and R is positive, and(M~n, g) is scalar flat if(M~n, g) is a complete noncompact manifold with nonnegative scalar curvature and finite L~p-norm of R?m. We prove that(M~n, g) is isometric to a spherical space form if for p ≥n/2, the L~p-norm of R?m is sufficiently small and R is positive.In particular, we prove that(M~n, g) is isometric to a spherical space form if for p ≥ n, R is positive and the L~p-norm of R?m is pinched in [0, C), where C is an explicit positive constant depending only on n, p, R and the Yamabe constant.     

Low Regularity Error Analysis for Weak Galerkin Finite Element Methods for Second Order Elliptic Problems

This paper presents error estimates in both an energy norm and the $L^2$-norm for the weak Galerkin (WG) finite element methods for elliptic problems with low regularity solutions. The error analysis for the continuous Galerkin finite element remains same regardless of regularity. A totally different analysis is needed for discontinuous finite element methods if the elliptic regularity is lower than H-1.5. Numerical results confirm the theoretical analysis.     

Stabilized Nonconforming Mixed Finite Element Method for Linear Elasticity on Rectangular or Cubic Meshes

Based on the primal mixed variational formulation, a stabilized nonconforming mixed finite element method is proposed for the linear elasticity on rectangular and cubic meshes. Two kinds of penalty terms are introduced in the stabilized mixed formulation, which are the jump penalty term for the displacement and the divergence penalty term for the stress. We use the classical nonconforming rectangular and cubic elements for the displacement and the discontinuous piecewise polynomial space for the stress, where the discrete space for stress are carefully chosen to guarantee the well-posedness of discrete formulation. The stabilized mixed method is locking-free. The optimal convergence order is derived in the $L^2$-norm for stress and in the broken $H^1$-norm and $L^2$-norm for displacement. A numerical test is carried out to verify the optimal convergence of the stabilized method.     

Finite difference/$H^1$-Galerkin MFE procedure for a fractional water wave model

In this article, an $H^1$-Galerkin mixed finite element (MFE) method for solving the time fractional water wave model is presented. First-order backward Euler difference method and $L1$ formula are applied to approximate integer derivative and Caputo fractional derivative with order $1/2$, respectively, and $H^1$-Galerkin mixed finite element method is used to approximate the spatial direction. The analysis of stability for fully discrete mixed finite element scheme is made and the optimal space-time orders of convergence for two unknown variables in both $H^1$-norm and $L^2$-norm are derived. Further, some computing results for a priori analysis and numerical figures based on four changed parameters in the studied problem are given to illustrate the effectiveness of the current method     

A superconvergent $L^{\infty}$-error estimate of RT1 mixed methods for elliptic control problems with an integral constraint

In this paper, we investigate the superconvergence property of mixed finite element methods for a linear elliptic control problem with an integral constraint. The state and co-state are approximated by the order $k=1$ Raviart-Thomas mixed finite element spaces and the control variable is approximated by piecewise constant functions. A superconvergent approximation of the control variable $u$ will be constructed by a projection of the discrete adjoint state. It is proved that this approximation have convergence order $h^{2}$ in $L^{\infty}$-norm. Finally, a numerical example is given to demonstrate the theoretical results.     

Unconditionally Optimal Error Analysis of the Second-Order BDF Finite Element Method for the Kuramoto-Tsuzuki Equation

This paper aims to study a second-order semi-implicit BDF finite element scheme for the Kuramoto-Tsuzuki equations in two dimensional and three dimensional spaces. The proposed scheme is stable and the nonlinear term is linearized by the extrapolation technique. Moreover, we prove that the error estimate in $L^2$-norm is unconditionally optimal which means that there has not any restriction on the time step and the mesh size. Finally, numerical results are displayed to illustrate our theoretical analysis.     

A global superconvergent $L^{\infty}$-error estimate of mixed finite element methods for semilinear elliptic optimal control problems

In this paper, we discuss the superconvergence of mixed finite element methods for a semilinear elliptic control problem with an integral constraint. The state and costate are approximated by the order $k=1$ Raviart-Thomas mixed finite element spaces and the control variable is approximated by piecewise constant functions. Approximation of the optimal control of the continuous optimal control problem will be constructed by a projection of the discrete adjoint state. It is proved that this approximation has convergence order $h^{2}$ in $L^{\infty}$-norm. Finally, a numerical example is given to demonstrate the theoretical results.     

A Two-Grid Finite Element Method for Time-Dependent Incompressible Navier-Stokes Equations with Non-Smooth Initial Data

We analyze here, a two-grid finite element method for the two dimensional time-dependent incompressible Navier-Stokes equations with non-smooth initial data. It involves solving the non-linear Navier-Stokes problem on a coarse grid of size $H$ and solving a Stokes problem on a fine grid of size $h, h <相似文献