速度—压力的Q2—Q1有限元

本文讨论了Stokes方程的Q_2-Q_1有限元,即速度空间采用双二次分片多项式插值,压力空间采用双一次分片多项式播值.在不满足经典的Babuska-Brezzi条件下,本注记进一步讨论了混合有限元方法和简化积分的罚方法,当解光滑性加强时,分别得到最优阶误差估计式|u-u_h|=O(h_2)及|u-u_h~2|_1=O(h~(2+)),改进了G.F.Carey,J.T.Oden等的结果.     

A penalty‐FEM for navier‐stokes type variational inequality with nonlinear damping term

In this article, we consider a penalty finite element (FE) method for incompressible Navier‐Stokes type variational inequality with nonlinear damping term. First, we establish penalty variational formulation and prove the well‐posedness and convergence of this problem. Then we show the penalty FE scheme and derive some error estimates. Finally, we give some numerical results to verify the theoretical rate of convergence. © 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 918–940, 2017     

Stokes问题低阶非协调混合元的改进加罚算法(英文)

通过对由经典加罚算法得到的两个解进行线性组合,研究Stokes方程低阶非协调混合元的改进加罚算法.该方法利用较大的罚参数能得到同使用较小参数的经典加罚方法一样的收敛阶.此外,基于单元的特性和插值后处理技巧,得到一些超收敛结果,从而改进以往的文献结果.     

Penalty finite element approximations for the Stokes equations by L2 projection

This paper considers the penalty finite element method for the Stokes equations, based on some stable finite elements space pair (X_h, M_h) that do satisfy the discrete inf–sup condition. Theoretical results show that the penalty error converges as fast as one should expect from the order of the elements. Moreover, the penalty finite element method by L² projection can improve the penalty error estimates. Finally, we confirm these results by a series of numerical experiments. Copyright © 2008 John Wiley & Sons, Ltd.     

On the Stability of Biquadratic-Bilinear Velocity-Pressure Finite Elements

In this paper, the stability of biquadratic-bilinear velocity-pressure finite elements are discussed for the stationary Stokes problem. It is proved that there exist constants $c ,c' ＞0$ independent of $h$ such thatHence a question in [1] is answered.     

四阶障碍问题的稳定化混合有限元方法

本文讨论了四阶障碍问题的稳定化混合有限元方法.首先,引入网格依赖范数,通过加罚方法得到了与四阶障碍问题的等价的混合变分形式.随后给出了基于C~0协调有限元空间(W_h,V_h)的混合有限元逼近,例如P_k-P_k三角形有限元.在网格依赖范数下,(W_h,V_h)满足离散的inf-sup条件.最后,我们在不同的假设下,得到了一些误差估计.     

Composite penalty method of a low order anisotropic nonconforming finite element for the Stokes problem

Composite penalty method of a low order anisotropic nonconforming quadrilateral finite element for the Stokes problem is presented. This method with a large penalty parameter can achieve the same accuracy as the stand method with a small penalty parameter and the convergence rate of this method is two times as that of the standard method under the condition of the same order penalty parameter. The superconvergence for velocity is established as well. The results of this paper are also valid to the most of the known nonconforming finite element methods.     

A penalty method for the time-dependent Stokes problem with the slip boundary condition and its finite element approximation

We consider the finite element method for the time-dependent Stokes problem with the slip boundary condition in a smooth domain. To avoid a variational crime of numerical computation, a penalty method is introduced, which also facilitates the numerical implementation. For the continuous problem, the convergence of the penalty method is investigated. Then we study the fully discretized finite element approximations for the penalty method with the P1/P1-stabilization or P1b/P1 element. For the discretization of the penalty term, we propose reduced and non-reduced integration schemes, and obtain an error estimate for velocity and pressure. The theoretical results are verified by numerical experiments.     

高波数Helmholtz方程的有限元方法和连续内罚有限元方法

本文介绍高波数Helmholtz方程的有限元方法和连续内罚有限元方法.将以线性元情形为例,给出方法的明显依赖于波数k的预渐近稳定性和误差分析.我们将介绍三种证明方法.我们还讨论了内罚有限元方法的罚参数的选取以显著减少方法的污染误差·最后还给出数值例子验证理论结果.     

An unfitted interface penalty method for the numerical approximation of contrast problems

We aim to approximate contrast problems by means of a numerical scheme which does not require that the computational mesh conforms with the discontinuity between coefficients. We focus on the approximation of diffusion-reaction equations in the framework of finite elements. In order to improve the unsatisfactory behavior of Lagrangian elements for this particular problem, we resort to an enriched approximation space, which involves elements cut by the interface. Firstly, we analyze the H¹-stability of the finite element space with respect to the position of the interface. This analysis, applied to the conditioning of the discrete system of equations, shows that the scheme may be ill posed for some configurations of the interface. Secondly, we propose a stabilization strategy, based on a scaling technique, which restores the standard properties of a Lagrangian finite element space and results to be very easily implemented. We also address the behavior of the scheme with respect to large contrast problems ending up with a choice of Nitsche?s penalty terms such that the extended finite element scheme with penalty is robust for the worst case among small sub-elements and large contrast problems. The theoretical results are finally illustrated by means of numerical experiments.     

求解定常不可压Stokes方程的两层罚函数方法

借助于两套有限元网格空间提出了一种求解定常不可压Stokes方程的两层罚函数方法.该方法只需要求解粗网格空间上的Stokes方程和细网格空间上的两个易于求解的罚参数方程（离散后的线性方程组具有相同的对称正定系数矩阵）.收敛性分析表明粗网格空间相对于细网格空间可以选择很小，并且罚参数的选取只与粗网格步长和问题的正则性有关.因此罚参数不必选择很小仍能够得到最优解.最后通过数值算例验证了上述理论结果，并且数值对比可知两层罚函数方法对于求解定常不可压Stokes方程具有很好的效果.     

A partially penalty immersed interface finite element method for anisotropic elliptic interface problems

In this article, we develop a partially penalty immersed interface finite element (PIFE) method for a kind of anisotropy diffusion models governed by the elliptic interface problems with discontinuous tensor‐coefficients. This method is based on linear immersed interface finite elements (IIFE) and applies the discontinuous Galerkin formulation around the interface. We add two penalty terms to the general IIFE formulation along the sides intersected with the interface. The flux jump condition is weakly enforced on the smooth interface. By proving that the piecewise linear function on an interface element is uniquely determined by its values at the three vertices under some conditions, we construct the finite element spaces. Therefore, a PIFE procedure is proposed, which is based on the symmetric, nonsymmetric or incomplete interior penalty discontinuous Galerkin formulation. Then we prove the consistency and the solvability of the procedure. Theoretical analysis and numerical experiments show that the PIFE solution possesses optimal‐order error estimates in the energy norm and norm.© 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1984–2028, 2014     

Optimal Error Estimate of the Penalty Finite Element Method for the Micropolar Fluid Equations

We present an optimal error estimate of the numerical velocity, pressure, and angular velocity for the fully discrete penalty finite element method of the micropolar equations when the parameters ?, Δ t, and h are sufficiently small. In order to obtain this estimate, we present the time discretization of the penalty micropolar equation that is based on the backward Euler scheme; the spatial discretization of the time discretized penalty micropolar equation is based on a finite elements space pair (X _h , M _h ) that satisfies some approximations properties.     

Domain Decomposition Preconditioners for Discontinuous Galerkin Discretizations of Compressible Fluid Flows

In this article we consider the application of Schwarz-type domain decomposition preconditioners to the discontinuous Galerkin finite element approximation of the compressible Navier-Stokes equations. To discretize this system of conservation laws, we exploit the (adjoint consistent) symmetric version of the interior penalty discontinuous Galerkin finite element method. To define the necessary coarse-level solver required for the definition of the proposed preconditioner, we exploit ideas from composite finite element methods, which allow for the definition of finite element schemes on general meshes consisting of polygonal (agglomerated) elements. The practical performance of the proposed preconditioner is demonstrated for a series of viscous test cases in both two- and three-dimensions.     

Morley-Wang-Xu element methods with penalty for a fourth order elliptic singular perturbation problem

Two Morley-Wang-Xu element methods with penalty for the fourth order elliptic singular perturbation problem are proposed in this paper, including the interior penalty Morley-Wang-Xu element method and the super penalty Morley-Wang-Xu element method. The key idea in designing these two methods is combining the Morley-Wang-Xu element and penalty formulation for the Laplace operator. Robust a priori error estimates are derived under minimal regularity assumptions on the exact solution by means of some established a posteriori error estimates. Finally, we present some numerical results to demonstrate the theoretical estimates.     

Convergence of a penalty-finite element approximation for an obstacle problem

Summary This study establishes an error estimate for a penalty-finite element approximation of the variational inequality obtained by a class of obstacle problems. By special identification of the penalty term, we first show that the penalty solution converges to the solution of a mixed formulation of the variational inequality. The rate of convergence of the penalization is where is the penalty parameter. To obtain the error of finite element approximation, we apply the results obtained by Brezzi, Hager and Raviart for the mixed finite element method to the variational inequality.     

Convergence and Quasi-Optimality of an Adaptive Multi-Penalty Discontinuous Galerkin Method

An adaptive multi-penalty discontinuous Galerkin method (AMPDG) for the diffusion problem is considered. Convergence and quasi-optimality of the AMPDG are proved. Compared with the analyses for the adaptive finite element method or the adaptive interior penalty discontinuous Galerkin method, extra works is done to overcome the difficulties caused by the additional penalty terms.     

Error estimate of the P 1 nonconforming finite element method for the penalized unsteady Navier-Stokes equations

We consider a finite element method for the penalty formulation of the time dependent Navier-Stokes equations. Usually the improper choice of the finite element space will lead that the error estimate (inversely) depends on the penalty parameter e{\epsilon}. We use the classical P ₁ nonconforming finite element space for the spatial discretization. Optimal e{\epsilon}-uniform error estimations for both velocity and pressure are obtained.     

Analysis of the Fictitious Domain Method with an L 2-Penalty for Elliptic Problems

The L ²-penalty fictitious domain method is based on a reformulation of the original problem in a larger simple-shaped domain by introducing a discontinuous reaction term with a penalty parameter ε > 0. We first derive regularity results and some a priori estimates and then prove several error estimates. We also give several error estimates for discretization problems by the finite element and finite volume methods.     

惩罚函数法在三维等参杂交元的应用

提出在等参杂交元中用惩罚函数法引入平衡约束条件,具体讨论了惩罚函数法在三维等参杂交元中的运用,并提出采用分项罚数的方法,建立最佳的罚平衡杂交元模型．罚平衡法可以在不增加自由度的前提下,有效地扼制寄生应力．数值实验表明,新建立的单元,可以有效地抑制单元畸变对计算精度的影响,从而大幅度提高畸变网格下的计算精度,方法带有普遍性．