Analysis of an augmented mixed‐primal formulation for the stationary Boussinesq problem

In this article, we propose and analyze a new mixed variational formulation for the stationary Boussinesq problem. Our method, which uses a technique previously applied to the Navier–Stokes equations, is based first on the introduction of a modified pseudostress tensor depending nonlinearly on the velocity through the respective convective term. Next, the pressure is eliminated, and an augmented approach for the fluid flow, which incorporates Galerkin‐type terms arising from the constitutive and equilibrium equations, and from the Dirichlet boundary condition, is coupled with a primal‐mixed scheme for the main equation modeling the temperature. In this way, the only unknowns of the resulting formulation are given by the aforementioned nonlinear pseudostress, the velocity, the temperature, and the normal derivative of the latter on the boundary. An equivalent fixed‐point setting is then introduced and the corresponding classical Banach Theorem, combined with the Lax–Milgram Theorem and the Babu?ka–Brezzi theory, are applied to prove the unique solvability of the continuous problem. In turn, the Brouwer and the Banach fixed‐point theorems are used to establish existence and uniqueness of solution, respectively, of the associated Galerkin scheme. In particular, Raviart–Thomas spaces of order k for the pseudostress, continuous piecewise polynomials of degree ≤ k+1 for the velocity and the temperature, and piecewise polynomials of degree ≤ k for the boundary unknown become feasible choices. Finally, we derive optimal a priori error estimates, and provide several numerical results illustrating the good performance of the augmented mixed‐primal finite element method and confirming the theoretical rates of convergence. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 445–478, 2016     

非线性双相滞热传导方程的新混合元超收敛分析

针对非线性双相滞热传导方程,建立了一种自由度少且自然满足B-B条件的新混合元逼近格式.在半离散格式下,基于双线性元的高精度结果,分别导出了原始变量的H~1模及中间变量的L~2模的超逼近性质,进而,借助于插值后处理算子,得到了原始及中间变量比传统误差高一阶的整体超收敛结果.     

Mixed finite element methods for incompressible flow: Stationary Stokes equations

In this article, we develop and analyze a mixed finite element method for the Stokes equations. Our mixed method is based on the pseudostress‐velocity formulation. The pseudostress is approximated by the Raviart‐Thomas (RT) element of index k ≥ 0 and the velocity by piecewise discontinuous polynomials of degree k. It is shown that this pair of finite elements is stable and yields quasi‐optimal accuracy. The indefinite system of linear equations resulting from the discretization is decoupled by the penalty method. The penalized pseudostress system is solved by the H(div) type of multigrid method and the velocity is then calculated explicitly. Alternative preconditioning approaches that do not involve penalizing the system are also discussed. Finally, numerical experiments are presented. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Guaranteed velocity error control for the pseudostress approximation of the Stokes equations

The pseudostress approximation of the Stokes equations rewrites the stationary Stokes equations with pure (but possibly inhomogeneous) Dirichlet boundary conditions as another (equivalent) mixed scheme based on a stress in H(div) and the velocity in L². Any standard mixed finite element function space can be utilized for this mixed formulation, e.g., the Raviart‐Thomas discretization which is related to the Crouzeix‐Raviart nonconforming finite element scheme in the lowest‐order case. The effective and guaranteed a posteriori error control for this nonconforming velocity‐oriented discretization can be generalized to the error control of some piecewise quadratic velocity approximation that is related to the discrete pseudostress. The analysis allows for local inf‐sup constants which can be chosen in a global partition to improve the estimation. Numerical examples provide strong evidence for an effective and guaranteed error control with very small overestimation factors even for domains with large anisotropy.© 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1411–1432, 2016     

A virtual element method-based positivity-preserving conservative scheme for convection–diffusion problems on polygonal meshes

In this paper, we propose a positivity-preserving conservative scheme based on the virtual element method (VEM) to solve convection–diffusion problems on general meshes. As an extension of finite element methods to general polygonal elements, the VEM has many advantages such as substantial mathematical foundations, simplicity in implementation. However, it is neither positivity-preserving nor locally conservative. The purpose of this article is to develop a new scheme, which has the same accuracy as the VEM and preserves the positivity of the numerical solution and local conservation on primary grids. The first step is to calculate the cell-vertex values by the lowest-order VEM. Then, the nonlinear two-point flux approximations are utilized to obtain the nonnegativity of cell-centered values and the local conservation property. The new scheme inherits both advantages of the VEM and the nonlinear two-point flux approximations. Numerical results show that the new scheme can reach the optimal convergence order of the virtual element theory, that is, the second-order accuracy for the solution and the first-order accuracy for its gradient. Moreover, the obtained cell-centered values are nonnegative, which demonstrates the positivity-preserving property of our new scheme.     

Navier－Stokes方程的变网格非协调有限元法

本文通过所谓的速度－压力型公式讨论了Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的变网格非协调有限元逼近，得到了在确定模意义下的速度、压力误差估计，且在一定条件下，某些误差估计能达到最优。     

Analysis of DG approximations for Stokes problem based on velocity‐pseudostress formulation

In this article, we first discuss the well posedness of a modified LDG scheme of Stokes problem, considering a velocity‐pseudostress formulation. The difficulty here relies on the fact that the application of classical Babu?ka‐Brezzi theory is not easy, so we proceed in a nonstandard way. For uniqueness, we apply a discrete version of Fredholm's alternative theorem, while the a priori error analysis is done introducing suitable projections of exact solution. As a result, we prove that the method is convergent, and under suitable regularity assumptions on the exact solution, the optimal rate of convergence is guaranteed. Next, we explore two stabilizations to the previous scheme, by adding least squares type terms. For these cases, well posedness and the a priori error estimates are proved by the application of standard theory. We end this work with some numerical experiments considering our third scheme, whose results are in agreement with the theoretical properties we deduce.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1540–1564, 2017     

An augmented stress‐based mixed finite element method for the steady state Navier‐Stokes equations with nonlinear viscosity

A new stress‐based mixed variational formulation for the stationary Navier‐Stokes equations with constant density and variable viscosity depending on the magnitude of the strain tensor, is proposed and analyzed in this work. Our approach is a natural extension of a technique applied in a recent paper by some of the authors to the same boundary value problem but with a viscosity that depends nonlinearly on the gradient of velocity instead of the strain tensor. In this case, and besides remarking that the strain‐dependence for the viscosity yields a more physically relevant model, we notice that to handle this nonlinearity we now need to incorporate not only the strain itself but also the vorticity as auxiliary unknowns. Furthermore, similarly as in that previous work, and aiming to deal with a suitable space for the velocity, the variational formulation is augmented with Galerkin‐type terms arising from the constitutive and equilibrium equations, the relations defining the two additional unknowns, and the Dirichlet boundary condition. In this way, and as the resulting augmented scheme can be rewritten as a fixed‐point operator equation, the classical Schauder and Banach theorems together with monotone operators theory are applied to derive the well‐posedness of the continuous and associated discrete schemes. In particular, we show that arbitrary finite element subspaces can be utilized for the latter, and then we derive optimal a priori error estimates along with the corresponding rates of convergence. Next, a reliable and efficient residual‐based a posteriori error estimator on arbitrary polygonal and polyhedral regions is proposed. The main tools used include Raviart‐Thomas and Clément interpolation operators, inverse and discrete inequalities, and the localization technique based on triangle‐bubble and edge‐bubble functions. Finally, several numerical essays illustrating the good performance of the method, confirming the reliability and efficiency of the a posteriori error estimator, and showing the desired behavior of the adaptive algorithm, are reported. © 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1692–1725, 2017     

非定常Stokes方程的全离散稳定化有限元格式

本文研究二维非定常Stokes方程全离散稳定化有限元方法.首先给出关于时间向后一步Euler半离散格式,然后直接从该时间半离散格式出发,构造基于两局部高斯积分的稳定化全离散有限元格式,其中空间用P_1—P_1元逼近,证明有限元解的误差估计.本文的研究方法使得理论证明变得更加简便,也是处理非定常Stokes方程的一种新的途径.     

A NEW NONCONFORMING MIXED FINITE ELEMENT SCHEME FOR THE STATIONARY NAVIER-STOKES EQUATIONS

In this article, a new stable nonconforming mixed finite element scheme is proposed for the stationary Navier-Stokes equations, in which a new low order CrouzeixRaviart type nonconforming rectangular element is taken for approximating space for the velocity and the piecewise constant element for the pressure. The optimal order error estimates for the approximation of both the velocity and the pressure in L2-norm are established, as well as one in broken H1-norm for the velocity. Numerical experiments are given which are consistent with our theoretical analysis.     

A Second Order Nonconforming Triangular Mixed Finite Element Scheme for the Stationary Navier-Stokes Equations

In this paper,a nonconforming triangular mixed finite element scheme with second order convergence behavior is proposed for the stationary Navier-Stokes equations.The new nonconforming triangular element is taken as approximation space for the velocity and the linear element for the pressure.The convergence analysis is presented and optimal error estimates of both broken H1-norm and L2-norm for velocity as well as the L2-norm for the pressure are derived.     

对流扩散方程一类改进的特征线修正有限元方法

１引言在地下水污染,地下渗流驱动,核污染,半导体等问题的数值模拟中,均涉及抛物型对流扩散方程（或方程组）的数值求解问题．这些对流扩散型偏微分方程（或方程组）具有共同的特点：对流的影响远大于扩散的影响,即对流占优性,对流占优性给问题的数值求解带来许多困难,因此对流占优问题的有效数值解法一直是计算数学中重要的研究内容．用通常的差分法或有限元法进行数值求解将出现数值振荡．为了克服数值振荡,提出各种迎风方法和修正的特征方法并在这些问题上得到成功的实际应用、８０年代,Ｄｏｕｇｌａｓ和Ｒｕｓｓｅｌｌ［２］等…     

Virtual element methods for parabolic problems on polygonal meshes

The virtual element method (VEM) is a recent technology that can make use of very general polygonal/polyhedral meshes without the need to integrate complex nonpolynomial functions on the elements and preserving an optimal order of convergence. In this article, we develop for the first time, the VEM for parabolic problems on polygonal meshes, considering time‐dependent diffusion as our model problem. After presenting the scheme, we develop a theoretical analysis and show the practical behavior of the proposed method through a large array of numerical tests. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 2110–2134, 2015     

A family of mixed finite elements for the elasticity problem

Summary A new mixed finite element formulation for the equations of linear elasticity is considered. In the formulation the variables approximated are the displacement, the unsymmetric stress tensor and the rotation. The rotation act as a Lagrange multiplier introduced in order to enforce the symmetry of the stress tensor. Based on this formulation a new family of both two-and three-dimensional mixed methods is defined. Optimal error estimates, which are valid uniformly with respect to the Poisson ratio, are derived. Finally, a new postprocessing scheme for improving the displacement is introduced and analyzed.     

A reduced FE formulation based on POD method for hyperbolic equations

A proper orthogonal decomposition (POD) method was successfully used in the reduced-order modeling of complex systems. In this paper, we extend the applications of POD method, namely, apply POD method to a classical finite element (FE) formulation for second-order hyperbolic equations with real practical applied background, establish a reduced FE formulation with lower dimensions and high enough accuracy, and provide the error estimates between the reduced FE solutions and the classical FE solutions and the implementation of algorithm for solving reduced FE formulation so as to provide scientific theoretic basis for service applications. Some numerical examples illustrate the fact that the results of numerical computation are consistent with theoretical conclusions. Moreover, it is shown that the reduced FE formulation based on POD method is feasible and efficient for solving FE formulation for second-order hyperbolic equations.     

Study of a finite element method for the time-dependent generalized Stokes system associated with viscoelastic flow

A three-field finite element scheme designed for solving systems of partial differential equations governing time-dependent viscoelastic flows is studied. Once a classical backward Euler time discretization is performed, the resulting three-field system of equations allows for a stable approximation of velocity, pressure and extra stress tensor, by means of continuous piecewise linear finite elements, in both two- and three- dimensional space. This is proved to hold for the linearized form of the system. An advantage of the new formulation is the fact that it provides an algorithm for the explicit iterative resolution of system nonlinearities. Convergence in an appropriate sense applying to these three flow fields is demonstrated.     

Stokes方程的压力梯度局部投影间断有限元法

本文对定常的Stokes方程提出了一种新的间断有限元法,通过将通常的间断Galerkin有限元法与压力梯度局部投影相结合,建立了一个稳定的间断有限元格式,对速度和压力的任意分片多项式空间P_l(K),P_m(K)的间断有限元逼近证明了解的存在唯一性,给出了关于速度和压力的L~2范数的最优误差估计.     

Optimal error estimates for the pseudostress formulation of the Navier–Stokes equations

In this paper, we prove optimal a priori error estimates for the pseudostress-velocity mixed finite element formulation of the incompressible Navier–Stokes equations, thus improve the result of Cai et al. (SINUM 2010). This is achieved by applying Petrov–Galerkin type Brezzi–Rappaz–Raviart theory.     

Virtual Elements for the Reissner‐Mindlin plate problem

In this work, we present a virtual element method for the approximation of the plate bending problem in the Reissner‐Mindlin formulation. The proposed method follows the MITC approach of the FEM context. We construct a family of VEM spaces with arbitrary degree of accuracy that satisfies the conditions of the MITC philosophy. We perform some numerical tests which allow us to assess the convergence and the robustness of the method.     

不可压混流驱动问题的最小二乘混合元方法

1引言考虑多孔介质中两相不可压缩可混溶渗流驱动问题,它是由一组非线性耦合的椭园型压力方程和抛物型浓度方程组成：dVV。—一山人V什）gVV却）一q,VEn,（．1）＆,,。＿．、。。—一。x）＿＋u·grade－dlv（D（u）grade）一（1－c）q－,xEn,tEJ,（1．2）＆”－－’”””‘”－”””——－’——,、—’一其中a（）一a（x,c）一是（x）／卢（c）,J一［0,Ti,DcyR‘为水平油藏区域．方程式（1．l）一（1．2）中各物理量的意义如下：广为流体压力,c为流体的浓度,u为流体的Darer速度,叶为源汇项,／一—。x（q,O）,…