基于对偶混合变分形式的Uzawa型算法

基于弹性接触问题的三变量（应力，位移，接触边界位移）对偶混合变分形式，对混合有限元离散化的单边约束问题，提出了一种Uzawa型算法。首先证明了迭代算法的收敛性，然后用数值例子验证了迭代算法的有效性。     

弹性接触问题的对偶混合有限元分析

１．引言用混合有限元方法求解弹性力学问题,其优点在于可同时求解位移和应力．力学问题的混合变分形式是混合有限元方法的基础．对于弹性接触问题,文献［６］给出了一种混合变分形式,以及相应的混合有限元分析（也可见［１０］）．其混合变分形式是直接从位移交分方程和Ｈｏｏｋ方程导出的,获得了应力ａ在Ｌ２（Ω）而位移、在Ｈ１（Ω）的一个闭凸子集上求解的混合变分问题．本文在［９］中提出的混合变分形式的基础上,再引入另一个Ｌａｐｒａｎｇｅ乘子,获得了三重组混合变分形式．它能同时求解物体内点的应力,位移和接触边界上的…     

一类线性约束变分不等式问题的幂罚函数法

本文研究了一类线性约束变分不等式(Ⅵ)的幂罚函数法求解问题.利用Ⅵ的KKT条件,将Ⅵ转化为等价的混合互补问题和一个新的Ⅵ问题,并在一定条件下分析了解的存在性和唯一性.利用度理论证明了幂罚方程组解的存在性与唯一性.由以上结果最终证明了幂罚函数法的收敛性,即幂罚方程组的解收敛于Ⅵ问题的解.     

弹性力学中单边接触问题的混合有限元法

给出了弹性力学中单边接触问题模型产生的变分不等式的一种新的混合有限元逼近. 用分片连续的P₂-P₁元来分别逼近位移场和接触域上的法应力, 在合理的正则性假设下得到最优收敛阶. 数值算例验证了该理论结果.     

一类非线性变分不等式及其数值逼近

本文讨论了一类带有两个变函数的非线性变分不等式的问题,以及它们的有限元逼近。所考虑的区域ΩR~n真不必为凸区域,只需边界分片光滑即可。文中证明了有限元逼近的收敛性,并给出了逼近的误差估计,特别是最大模估计。     

一类新的求解双边障碍问题的数值方法

针对双边障碍问题的离散互补形式,提出了一类新的格式将其等价转化为方程组的形式,并采用牛顿迭代法进行求解.实验结果显示所提算法能快速,有效地计算出数值解和接触集.     

美式期权定价问题的数值方法

本文研究美式股票看跌期权定价问题的数值方法。通过将问题转化为等价的变分不等式方程，分别建立了半离散和全离散有限元逼近格式。并给出了有限元解的收敛性和稳定性分析。数值实验表明本文算法是一个高效和收敛的算法。     

椭圆变分不等式离散问题的数值解法

本文对线性及非线性椭圆变分不等式的离散问题的数值解法作一简略的综述,并指出了各解法的若干优缺点及若干待研究的问题。     

不等式约束优化问题的一个势函数

基于Carroll(1961)建立的罚函数，本文给出了不等式约束优化问题的一个势函数，并且讨论了该函数的性质．最后证明了在此基础上建立的对偶算法具有Q-线性收敛性．     

具不等式约束变分不等式的信赖域算法

1　引　　言令Ｘ是Ｒｎ 中的非空闭凸集 ,Ｆ :Ｘ→Ｒｎ 是连续映射 ,〈· ,·〉表示Ｒｎ 中的内积 有限维变分不等式问题 (以下简称变分不等式问题 ,记为ＶＩＰ或ＶＩ(Ｘ ,Ｆ) ) :就是求ｘ ∈Ｒｎ,使ｘ ∈Ｘ且 ｘ ∈Ｘ ,〈Ｆ(ｘ ) ,ｘ -ｘ 〉≥ 0 . ( 1 )在Ｘ =Ｒｎ+ 的特殊情形下 ,( 1 )变为非线性互补问题 (记为ＮＣＰ或ＮＣＰ(Ｆ) ) :就是求ｘ ∈Ｒｎ,使ｘ ≥ 0 ,Ｆ(ｘ ) ≥ 0 ,且〈ｘ ,Ｆ(ｘ )〉 =0 . ( 2 )　　变分不等式长期以来一直用于阐述和研究经济学、控制论、交通运输等领域中出现的各种平衡模型 近二十年来 ,变分不等式及其…     

On a mixed finite element formulation of a second‐order quasilinear problem in the plane

We analyze a mixed finite element discretization of a second‐order quasilinear problem based on the Raviart‐Thomas space. We prove that the discrete problem is solvable and provide a local uniqueness result for the solution. We also obtain optimal order L²‐error estimates for both the scalar variable and the associated flux. The main feature of our method is that it is free from the boundness conditions required in previous works on the coefficients of the quasilinear operator. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 20: 90–103, 2004.     

Mixed finite element analysis of semi-coercive unilateral contact problems with given friction

A unilateral contact 2D-problem is considered provided one of two elastic bodies can shift in a given direction as a rigid body. Using Lagrange multipliers for both normal and tangential constraints on the contact interface, we introduce a saddle point problem and prove its unique solvability. We discretize the problem by a standard finite element method and prove a convergence of approximations. We propose a numerical realization on the basis of an auxiliary “ bolted” problem and the algorithm of Uzawa.     

An expanded mixed covolume method for sobolev equation with convection term on triangular grids

The expanded mixed covolume method for the two‐dimensional Sobolev equation with convection term is developed and studied. This method uses the lowest‐order Raviart‐Thomas mixed finite element space as the trial function space. By introducing a transfer operator γ_h which maps the trial function space into the test function space and combining expanded mixed finite element with mixed covolume method, the continuous‐in‐time, discrete‐in‐time expanded mixed covolume schemes are constructed, and optimal error estimates for these schemes are obtained. Numerical results are given to examine the validity and effectiveness of the proposed schemes.© 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

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.     

An extension of Karmarkar's projective algorithm for convex quadratic programming

We present an extension of Karmarkar's linear programming algorithm for solving a more general group of optimization problems: convex quadratic programs. This extension is based on the iterated application of the objective augmentation and the projective transformation, followed by optimization over an inscribing ellipsoid centered at the current solution. It creates a sequence of interior feasible points that converge to the optimal feasible solution in O(Ln) iterations; each iteration can be computed in O(Ln ³) arithmetic operations, wheren is the number of variables andL is the number of bits in the input. In this paper, we emphasize its convergence property, practical efficiency, and relation to the ellipsoid method.     

关于可靠性设施布局问题的近似算法

设施布局问题的研究始于20世纪60年代,主要研究选择修建设施的位置和数量,以及与需要得到服务的城市之间的分配关系,使得设施的修建费用和设施与城市之间的连接费用之和达到最小.现实生活中, 受自然灾害、工人罢工、恐怖袭击等因素的影响,修建的设施可能会出现故障, 故连接到它的城市无法得到供应,这就直接影响到了整个系统的可靠性.针对如何以相对较小的代价换取设施布局可靠性的提升,研究人员提出了可靠性设施布局问题.参考经典设施布局问题的贪婪算法、原始对偶算法和容错性问题中分阶段分层次处理的思想,设计了可靠性设施布局问题的一个组合算法.该算法不仅在理论上具有很好的常数近似度,而且还具有运算复杂性低的优点.这对于之前的可靠性设施布局问题只有数值实验算法, 是一个很大的进步.     

各向异性网格下二阶椭圆问题的一个混合元形式

提出了二阶椭圆问题的一个混合变分形式,同时证明了Rariart-Thomas元的各向异性插值性质,并给出了单元的对二阶问题的最优误差估计。     

CVD全离散化的混合有限元解的数值模拟

在已有的对CVD化学方程半离散化和全离散化混合有限元解的存在性及其误差分析的基础上,对其全离散化混合有限元解进行了数值模拟,结果进一步表明了混合有限元解的高精度、易于计算的良好性质.