Navier-Stokes方程的局部投影稳定化方法

将Matthies,Skrzypacz和Tubiska的思想从线性的Oseen方程拓展到了非线性的Navier-Stokes方程,针对不可压缩的定常Navier-Stokes方程,提出了一种局部投影稳定化有限元方法.该方法既克服了对流占优,又绕开了inf-sup条件的限制.给出的局部投影空间既可以定义在两种不同网格上,又可以定义在相同网格上.与其他两级方法相比,定义在同一网格空间上的局部投影稳定化格式更紧凑.在同一网格上,除了给出需要bubble函数来增强的逼近空间外,还特别考虑了两种不需要用bubble函数来增强的新的空间.基于一种特殊的插值技巧,给出了稳定性分析和误差估计.最后,还列举了两个数值算例,进一步验证了理论结果的正确性.     

关于“Bubble-隐定化有限元方法”的两点注记

本文针对Ｓｔｏｋｅｓ－问题给出了［１２］发展的基于局部ｂｕｂｂｌｅ－函数稳定化有限元法与Ｇａｌｓ－稳定化有限元法,ｂｕｂｂｌｅ函数扩充元法的等价关系     

关于“Bbble—隐定化有限元方法”的两点注记

本文针对Ｓｔｏｋｅｓ－问题给出了「１２」发展的基于局部ｂｕｂｂｌｅ－函数稳定化有限元与Ｇａｌｓ－稳定化有限元法，ｂｕｂｂｌｅ函数扩充元法的等价关系。     

Stokes问题基于泡函数的简化的稳定化混合元格式的收敛性

利用泡函数导出Stokes问题的两种新的，简化的稳定化混合有限元格式，并证明这些格式与通常带泡函数的稳定化格式具有相同的收敛性，但是自由度可以大大减少。     

对流占优问题稳定化的扩展混合有限元方法

讨论了对流占优问题稳定化的扩展混合元数值模拟.把稳定化的思想与扩展混合元方法相结合,既可以高精度逼近未知函数,未知函数的梯度及伴随向量函数,又能保证格式的稳定性.理论分析表明,方法是有效的,具有最优L²逼近精度.     

二阶椭圆问题基于泡函数的简化的稳定化二阶混合有限元格式

本文研究二阶椭圆方程基于泡函数的稳定化的二阶混合有限元格式,通过消去泡函数导出一种自由度很少的简化的稳定化的二阶混合有限元格式, 误差分析表明消去泡函数的简化格式与带有泡函数的格式具有相同的精度而可以节省6N_p个自由度(其中N_p三角形剖分中的顶点数目).       

平面弹性力学问题的混合元法的泡函数稳定性及其后验误差估计

该文讨论平面弹性力学问题的混合元法的泡函数稳定性,并导出基于简化的稳定化格式的一种先验误差估计和后验误差估计．这种简化的稳定化格式较通常的格式节省自由度．     

Bubble—函数稳定化混合有限元分析

本文针对混合结构抽象问题，基于「９」的非标准稳定化有限元方法的一般框架研究了ｂｕｂｂｌｅ－函数稳定化方法，该逼近代格式使得Ｂａｂｕｓｋａ－Ｂｒｅｚｚｉ条件是不必要的。     

关于线弹性问题混合格式稳定化思想的探讨

首先从混合有限元理论出发,探讨线弹性问题混合变分格式所满足的稳定性条件,从而保证解的存在唯一性.使用连续的分片线性函数和分片常数来分别逼近应力和位移,详细分析了混合格式下稳定化的必要性,有助于更加深入地了解稳定化的基本思想.然后,通过在混合格式中引入位移的跳跃惩罚项,展示了一个无闭锁稳定化混合有限元方法,并证明了此方法是稳定的且是线性收敛的.     

Navier-Stokes方程的低阶混合有限元方法

1引言 定常N-s方程是流体力学中一类非常重要的方程,而经典的混合有限元方法要求有限元空间组合满足B-B条件.这一条件限制了工程中常用的低阶有限元空间如:P1/P1,P,/Po等.为了去掉LBB条件限制,产生了一种新的方法--稳定化方法(也成CBB方法).1988年,F.Brezzi和J.Douglas.Jr对线性的Stokes方程建立了一种稳定化格式([2]).对于低阶的有限元应用压力投影稳定项构造了一种稳定化格式,并给出了格式的解存在唯一性,且给出了几种有限元的算例.     

伪双曲方程的新混合有限元方法

构造分析一类二阶伪双曲方程的H1-Galerkin扩展混合有限元方法,该方法采用了扩展混合有限元方法和H1-Galerkin混合有限元方法相结合的技巧.新的格式同时保持了扩展混合有限元方法和H1-Galerkin混合有限元方法的优点.该混合格式与标准的混合格式相比能同时逼近三个变量:未知函数、梯度和流量(系数乘以梯度),并且不必满足LBB相容性条件.     

粘弹性方程一种新的分裂正定混合元法

本文用分裂正定混合有限元方法研究二阶粘弹性方程. 首先构造一种新的分裂正定混合变分形式和基于这种分裂正定混合变分形式关于时间的半离散格式, 然后绕开关于空间变量的半离散化格式, 直接从时间半离散出发构造出全离散化的分裂正定混合有限元格式, 并给出这种分裂正定混合有限元解的误差估计. 这种研究思路使得理论论证变得更简单,这是处理二阶粘弹性方程的一种新的尝试.     

一类单边约束问题的数值方法

本文分别基于原始变分形式与对偶混合变分形式，对一类单边约束问题进行了数值求解，提出了求解离散对偶混合变分问题的Uzawa型算法，并用数值例子验证了算法的有效性．     

A stabilized mixed formulation for finite strain deformation for low-order tetrahedral solid elements

A stabilized mixed finite element formulation for four-noded tetrahedral elements is introduced for robustly solving small and large deformation problems. The uniqueness of the formulation lies within the fact that it is general in that it can be applied to any type of material model without requiring specialized geometric or material parameters. To overcome the problem of volumetric locking, a mixed element formulation that utilizes linear displacement and pressure fields was implemented. The stabilization is provided by enhancing the rate of deformation tensor with a term derived using a bubble function approach. The element was implemented through a user-programmable element of the commercial finite element software ANSYS. Using the ANSYS platform, the performance of the element was evaluated by comparing the predicted results with those obtained using mixed quadratic tetrahedral elements and hexahedral elements with a B-bar formulation. Based on the quality of the results, the new element formulation shows significant potential for use in simulating complex engineering processes.     

Une résolution par les elements finis mixtes à une inconnue par maille

To solve the groundwater flow equations, we show how to produce a scheme with one unblown per element shirting from a mixed formulation discretized with the Raviarl-Thomas triangular elements of lowest order. We study the new formulation in the elliptic case with sink/source terms in order to use mixed finite elements with less unknowns without any numerical intégration. The last part of the paper is aimed to study the positive definiteness of the matrix obtained with this new formulation.     

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.     

Coupling of mixed finite element and stabilized boundary element methods for a fluid–solid interaction problem in 3D

We introduce and analyze the coupling of a mixed finite element and a boundary element for a three‐dimensional time‐harmonic fluid–solid interaction problem. We consider a formulation in which the Cauchy stress tensor and the rotation are the main variables in the elastic structure and use the usual pressure formulation in the acoustic fluid. The mixed variational formulation in the solid is completed with boundary integral equations relating the Cauchy data of the acoustic problem on the coupling interface. A crucial point in our formulation is the stabilization technique introduced by Hiptmair and coworkers to avoid the well‐known instability issue appearing in the boundary element method treatment of the exterior Helmholtz problem. The main novelty of this formulation, with respect to a previous approach, consists in reducing the computational domain to the solid media and providing a more accurate treatment of the far field effect. We show that the continuous problem is well‐posed and propose a conforming Galerkin method based on the lowest‐order Arnold–Falk–Winther mixed finite element. Finally, we prove that the numerical scheme is convergent with optimal order.Copyright © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1211–1233, 2014     

对偶变量块体混合元及其位移元的收敛性和精度分析

弹性力学Hamilton正则方程和Hamilton混合元的等效刚度系数矩阵,均具有直观的辛特性.基于H R变分原理和弹性力学保辛理论建立的对偶变量块体混合元,其等效刚度系数矩阵同样具有直观的辛特性.根据对偶变量块体混合元列式,可直接建立问题的控制方程,进行混合法求解.同时,通过对偶变量块体混合元列式可以导出对偶变量块体位移元列式,建立问题的控制方程后,可先求位移的解.数值实例表明:线性8结点对偶变量块体位移减缩积分元的各力学量的收敛速度均衡、收敛过程稳定、结果精度高,其应力变量的收敛速度与传统的20结点位移协调减缩积分元接近.对偶变量块体位移元具有普适性.     

Stabilized Nonconforming Mixed Finite Element Method for Linear Elasticity on Rectangular or Cubic Meshes

Based on the primal mixed variational formulation, a stabilized nonconforming mixed finite element method is proposed for the linear elasticity on rectangular and cubic meshes. Two kinds of penalty terms are introduced in the stabilized mixed formulation, which are the jump penalty term for the displacement and the divergence penalty term for the stress. We use the classical nonconforming rectangular and cubic elements for the displacement and the discontinuous piecewise polynomial space for the stress, where the discrete space for stress are carefully chosen to guarantee the well-posedness of discrete formulation. The stabilized mixed method is locking-free. The optimal convergence order is derived in the $L^2$-norm for stress and in the broken $H^1$-norm and $L^2$-norm for displacement. A numerical test is carried out to verify the optimal convergence of the stabilized method.