线性抛物型积分微分方程的扩展混合体积元方法

1 引言 考虑线性抛物型积分微分方程初边值问题: {pt(x,t)-▽.{A(x,t)▽p(x,t) +∫t0 B(x,t,τ)▽p(x,τ)dτ}=f(x,t),(x,t)∈Ω×(0,T],(1.1) p(x,0):p0(x), x∈Ω, p(x,t)=0, (x,t)∈(a)Ω×(0,T]. 这里x=(x,y),Ω=(a,b)×(c,d),(e)Ω是区域Ω的边界,p为未知函数,A=(aij)2×2为已知的对称正定矩阵,B=(bij)2×2为已知矩阵,而且aij,bij,(aij)t(i,j=1,2)光滑有界,f∈L2(Ω).     

Optimal error analysis of Crank–Nicolson lowest-order Galerkin-mixed finite element method for incompressible miscible flow in porous media

Numerical methods for incompressible miscible flow in porous media have been studied extensively in the last several decades. In practical applications, the lowest-order Galerkin-mixed method is the most popular one, where the linear Lagrange element is used for the concentration and the lowest order Raviart–Thomas mixed element pair is used for the Darcy velocity and pressure. The existing error estimate of the method in L² -norm is in the order in spatial direction, which however is not optimal and valid only under certain extra restrictions on both time step and spatial meshes, excluding the most commonly used mesh h = h_p = h_c . This paper focuses on new and optimal error estimates of a linearized Crank–Nicolson lowest-order Galerkin-mixed finite element method (FEM), where the second-order accuracy for the concentration in both time and spatial directions is established unconditionally. The key to our optimal error analysis is an elliptic quasi-projection. Moreover, we propose a simple one-step recovery technique to obtain a new numerical Darcy velocity and pressure of second-order accuracy. Numerical results for both two and three-dimensional models are provided to confirm our theoretical analysis.     

Error Estimates for Mixed Finite Element Methods for Sobolev Equation

1 IntroductionLet fl be a bounded domain in R2 with Lipschitz continuous boundaxy 0fl. For thed0 < T < co, we consider the fo1lowing initial-boun'lar}-ralue problem for thc Sobolevequation:where ut denotes the time derivative of the function (1. Vu denotes the gradient of thefunction u, and divv denotes the divergence of the vect{Jr tulued function v, a1 b1, f, anduo are known functions.The standard finite element method for (1.1) (1.3) llas received considerable attentionand is well studied…     

On the mixed finite element method with Lagrange multipliers

In this note we analyze a modified mixed finite element method for second‐order elliptic equations in divergence form. As a model we consider the Poisson problem with mixed boundary conditions in a polygonal domain of R ². The Neumann (essential) condition is imposed here in a weak sense, which yields the introduction of a Lagrange multiplier given by the trace of the solution on the corresponding boundary. This approach allows to handle nonhomogeneous Neumann boundary conditions, theoretically and computationally, in an alternative and usually easier way. Then we utilize the classical Babu?ka‐Brezzi theory to show that the resulting mixed variational formulation is well posed. In addition, we use Raviart‐Thomas spaces to define the associated finite element method and, applying some elliptic regularity results, we prove the stability, unique solvability, and convergence of this discrete scheme, under appropriate assumptions on the mesh sizes. Finally, we provide numerical results illustrating the performance of the algorithm for smooth and singular problems. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 192–210, 2003     

半线性Sobolev方程的H~1-Galerkin混合有限元方法

利用H~1-Galerkin混合有限元方法研究了一维半线性Sobolev方程,得到了半离散解的最优阶误差估计,优点是不需验证LBB相容性条件.     

Sobolev 方程的$H^1$-Galerkin混合有限元方法

对Sobolev方程采用H1-Galerkin混合有限元方法进行数值模拟．给出了一维空间中该方法的半离散和全离散格式及其最优误差估计;并将该方法推广到二维和三维空间．与H1-Galerkin有限元方法相比,该方法不仅降低了对有限元空间的连续性要求;而且与传统的混合有限元方法具有相同的收敛阶,但其有限元空间的选取却不需要满足LBB相容条件．数值例子将进一步说明该方法的可行性与有效性．     

A new mixed scheme based on variation of constants for Sobolev equation with nonlinear convection term

A new mixed scheme which combines the variation of constants and the H ¹-Galerkin mixed finite element method is constructed for nonlinear Sobolev equation with nonlinear convection term. Optimal error estimates are derived for both semidiscrete and fully discrete schemes. Finally, some numerical results are given to confirm the theoretical analysis of the proposed method.     

四阶强阻尼波动方程的混合控制体积法

本文利用混合控制体积方法在三角网格剖分下求解四阶强阻尼波动方程.通过使用最低阶Raviart-Thomas混合有限元空间和引入迁移算子把解函数空间映射成试探函数空间,构造了半离散和全离散的混合控制体积格式,得到了最优阶误差估计.     

Optimal error estimates for the hp-version interior penalty discontinuous Galerkin finite element method

We consider the hp-version interior penalty discontinuous Galerkinfinite-element method (hp-DGFEM) for second-order linear reaction–diffusionequations. To the best of our knowledge, the sharpest knownerror bounds for the hp-DGFEM are due to Rivière et al.(1999,Comput. Geosci., 3, 337–360) and Houston et al.(2002,SIAM J. Numer. Anal., 99, 2133–2163). These are optimalwith respect to the meshsize h but suboptimal with respect tothe polynomial degree p by half an order of p. We present improvederror bounds in the energy norm, by introducing a new functionspace framework. More specifically, assuming that the solutionsbelong element-wise to an augmented Sobolev space, we deducefully hp-optimal error bounds.     

A Priori and a Posteriori Error Analysis of a Wavelet-Based Stabilization for the Mixed Finite Element Method

We use Galerkin least-squares terms and biorthogonal wavelet bases to develop a new stabilized dual-mixed finite element method for second-order elliptic equations in divergence form with Neumann boundary conditions. The approach introduces the trace of the solution on the boundary as a new unknown that acts also as a Lagrange multiplier. We show that the resulting stabilized dual-mixed variational formulation and the associated discrete scheme defined with Raviart–Thomas spaces are well-posed and derive the usual a priori error estimates and the corresponding rate of convergence. Furthermore, a reliable and efficient residual-based a posteriori error estimator and a reliable and quasi-efficient one are provided.     

Mixed covolume method for parabolic problems on triangular grids

In this paper, we present a mixed covolume method for parabolic equations on triangular grids. This method use the lowest order Raviart–Thomas (R–T) mixed finite element space as the trial space. We prove the optimal order of convergence for the approximate pressure and velocity in L²-norm. Furthermore, we obtain the quasi-optimal error estimates for the approximate pressure in L^∞-norm.     

Sobolev方程的连续时空有限元方法

本文研究Sobolev方程的连续时空有限元方法.首先建立Sobolev方程的连续时空有限元格式,然后证明了解的存在唯一性和稳定性并给出连续时空有限元解各种范数下的误差估计.最后给出数值算例来验证理论分析的正确性,并进一步说明本文所建立的格式关于时间可以得到比传统有限元方法更高的精度.     

Sobolev方程的混合连续时空有限元解的误差估计

通过引入辅助变量构造Sobolev方程的混合连续时空有限元离散格式,使得该格式既利用混合法将方程降阶,又将时间和空间两个变量同时用有限元方法离散,从而获得时空形式高精度数值模型.证明了Sobolev方程混合时空有限元解的存在唯一性、稳定性,并利用时间和空间投影算子推导出时空数值解的误差估计.     

定常的磁流体力学方程的非线性Galerkin混合元法

提出了定常的磁流体动力学方程的一种非线性Galerkin混合元法,并导出非线性Galerkin混合元解的存在性和误差估计．