Macro-Hybrid Variational Formulations of Constrained Boundary Value Problems

In the context of convex analysis, macro-hybrid variational formulations of constrained boundary value problems are presented. Monotone mixed variational inclusions are macro-hybridized on the basis of nonoverlapping domain decompositions, and corresponding three-field versions are derived. Then, for regularization purposes, augmented formulations are established via preconditioned exact penalizations and expressed in terms of proximation operators. Optimization interpretations are given for potential problems, recovering the classic two- and three-field augmented Lagrangian formulations. Furthermore, associated parallel two- and three-field proximal-point algorithms are discussed for numerical resolution of finite element discretizations. Applications to dual mixed variational formulations of problems from mechanics illustrate the theory.     

A New Finite Element Space for Expanded Mixed Finite Element Method

In this article, we propose a new finite element space Λ$_h$ for the expanded mixed finite element method (EMFEM) for second-order elliptic problems to guarantee its computing capability and reduce the computation cost. The new finite element space Λ$_h$ is designed in such a way that the strong requirement V$_h\subset$Λ$_h$ in [9] is weakened to {v$_h\in$V$_h$; divv$_h$=0}$\subset$Λ$_h$ so that it needs fewer degrees of freedom than its classical counterpart. Furthermore, the new Λ$_h$ coupled with the Raviart-Thomas space satisfies the inf-sup condition, which is crucial to the computation of mixed methods for its close relation to the behavior of the smallest nonzero eigenvalue of the stiff matrix, and thus the existence, uniqueness and optimal approximate capability of the EMFEM solution are proved for rectangular partitions in $\mathbb{R}^d, d=2,3$ and for triangular partitions in $\mathbb{R}^2$. Also, the solvability of the EMFEM for triangular partition in $\mathbb{R}^3$ can be directly proved without the inf-sup condition. Numerical experiments are conducted to confirm these theoretical findings.     

Navier-Stokes问题的六面体单元的混合有限元法

本文.我们给出三维空间的Navier-Stokes问题的一种新的六面体单元的混合有限元格式.     

三维Stokes问题各向异性混合元分析

本文提出了一个一般的立方体单元格式并将其应用到三维Stokes问题的混合有限元逼近,给出了各向异性插值误差估计,相容误差估计和LBB条件成立的验证,从而证明了其在不满足正则性和拟一致条件下的收敛性．另外我们还得到了其一个特殊收敛性质,即在解(u,p)∈(H3(Ω))3×H2(Ω)时,相容误差阶为O(h2max),比插值误差阶O(hmax)高一阶．     

发展基于RIESZ表示算子的最小二乘混合有限元方法

针对一类广义抽象混合结构变分方程,研究了其解的存在唯一性问题;同时发展基于Riesz表示算子的有限元方法.     

Convergence of the Augmented Lagrangian Method for Nonlinear Optimization Problems over Second-Order Cones

The paper analyzes the rate of local convergence of the augmented Lagrangian method for nonlinear second-order cone optimization problems. Under the constraint nondegeneracy condition and the strong second order sufficient condition, we demonstrate that the sequence of iterate points generated by the augmented Lagrangian method locally converges to a local minimizer at a linear rate, whose ratio constant is proportional to 1/τ with penalty parameter τ not less than a threshold . Importantly and interestingly enough, the analysis does not require the strict complementarity condition. Supported by the National Natural Science Foundation of China under Project 10771026 and by the Scientific Research Foundation for the Returned Overseas Chinese Scholars, State Education Ministry.     

拟线性抛物型积分微分方程的一个新最低阶混合元格式

对一类拟线性抛物型积分微分方程构造了一个新的最低阶三角形协调混合元格式,并直接利用单元插值的性质,给出了相应的收敛性分析和H~1-模及L~2-模意义下的最优误差估计.     

退化椭圆问题的最小二乘混合有限元方法

本文研究了电磁场中关于共振现象的一类退化的椭圆问题 ,提出了最小二乘混合有限元方法 .这一方法的好处是可以去掉传统混合元空间的LBB条件所得到的系数矩阵是对称正定的 ,使得法语解更加方便 .本文得到了最小二乘混合有限元方法的L2 和H1估计 .     

Method of Nonconforming Mixed Finite Element for Second Order Elliptic Problems

In this paper, the method of non-conforming mixed finite element for second order elliptic problems is discussed and a format of real optimal order for the lowest order error estimate.     

Preconditioning Mixed Finite Element Saddle-point Elliptic Problems

We consider saddle-point problems that typically arise from the mixed finite element discretization of second-order elliptic problems. By proper equivalent algebraic operations the considered saddle-point problem is transformed to another saddle-point problem. The resulting problem can then be efficiently preconditioned by a block-diagonal matrix or by a factored block-matrix (the blocks correspond to the velocity and pressure, respectively). Both preconditioners have a block on the main diagonal that corresponds to the bilinear form (δ is a positive parameter) and a second block that is equal to a constant times the identity operator. We derive uniform bounds for the negative and positive eigenvalues of the preconditioned operator. Then any known preconditioner for the above bilinear form can be applied. We also show some numerical experiments that illustrate the convergence properties of the proposed technique.     

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

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

Finite Element Solution of Conical Diffraction Problems

This paper is devoted to the numerical study of diffraction by periodic structures of plane waves under oblique incidence. For this situation Maxwell's equations can be reduced to a system of two Helmholtz equations in R ² coupled via quasiperiodic transmission conditions on the piecewise smooth interfaces between different materials. The numerical analysis is based on a strongly elliptic variational formulation of the differential problem in a bounded periodic cell involving nonlocal boundary operators. We obtain existence and uniqueness results for discrete solutions and provide the corresponding error analysis.     

Extrapolation of Mixed Finite Element Approximations for the Maxwell Eigenvalue Problem

In this paper,a general method to derive asymptotic error expansion formulas for the mixed finite element approximations of the Maxwell eigenvalue problem is established.Abstract lemmas for the error of the eigenvalue approximations are obtained.Based on the asymptotic error expansion formulas,the Richardson extrapolation method is employed to improve the accuracy of the approximations for the eigenvalues of the Maxwell system from θ(h2) to θ(h4) when applying the lowest order Nédé1ec mixed finite element and a nonconforming mixed finite element.To our best knowledge,this is the first superconvergence result of the Maxwell eigenvalue problem by the extrapolation of the mixed finite element approximation.Numerical experiments are provided to demonstrate the theoretical results.     

Stokes问题非协调混合有限元超收敛分析

本文通过引入全新的技巧,研究了Stokes问题的非协调混合有限元方法,得到了关于速度与压力的超逼近性质.进一步地通过构造一个恰当的插值后处理算子,还得到了关于速度的整体超收敛结果.     

A Priori Error Estimates for Least-Squares Mixed Finite Element Approximation of Elliptic Optimal Control Problems

In this paper, a constrained distributed optimal control problem governed by a first-order elliptic system is considered. Least-squares mixed finite element methods, which are not subject to the Ladyzhenkaya-Babuska-Brezzi consistency condition, are used for solving the elliptic system with two unknown state variables. By adopting the Lagrange multiplier approach, continuous and discrete optimality systems including a primal state equation, an adjoint state equation, and a variational inequality for the optimal control are derived, respectively. Both the discrete state equation and discrete adjoint state equation yield a symmetric and positive definite linear algebraic system. Thus, the popular solvers such as preconditioned conjugate gradient (PCG) and algebraic multi-grid (AMG) can be used for rapid solution. Optimal a priori error estimates are obtained, respectively, for the control function in $L^2(Ω)$-norm, for the original state and adjoint state in $H^1(Ω)$-norm, and for the flux state and adjoint flux state in $H$(div; $Ω$)-norm. Finally, we use one numerical example to validate the theoretical findings.     

Stokes方程非协调混合元的特征值下界

通过利用Crouzeix-Raviart元({1,x,y}),旋转元({1,x,y,x~2-y~2}),拓广旋转元({1,x,y,x~2,y~2})以及拓广Crouzeix-Raviart元({1,x,y,x~2+y~2})这四种混合有限元(参看正文中示图)来提供求Stokes特征值下界的方法.并找到恰当的理论框架,重要的是证明不仅统一,而且出奇的短,仅需几行.最后给出相关的数值结果来验证本文的理论分析.     

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

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

Stokes问题Q_2-P_1混合元外推方法

考虑Stokes问题的有限元解与精确解插值的Q2-P1混合元的渐近误差展开和分裂外推.首先利用积分恒等式技巧确定了微分方程精确解与有限元插值之间积分式的主项,其次再借助插值后处理和分裂外推技术,得到了比通常的误差估计提高两阶的收敛速度.     

蒸汽沉淀化学反应方程混合元法的离散格式研究

蒸汽沉淀化学反应过程有着极其广泛的应用,其数学模型归结为一个包含流速场,温度场,压力场和气体溶质场的非线性偏微分方程组．用混合有限元方法研究蒸汽沉淀化学反应方程组,导出其半离散化和全离散化的混合元格式,并证明这些格式的解的存在性和收敛性(误差估计)．用混合元法处理究蒸汽沉淀化学反应方程组,可以同时求出流速场,温度场,压力场和气体溶质场的数值解. 因此该研究既具有重要的理论意义,又具有广泛的应用前景．