双调和方程混合元的一种新格式

本文介绍了双调和方程混合元的一种新格式,用双二次多项式逼近流函数,双一次多项式逼近涡函数.在拟一致矩形剖分的条件下,证明了此格式具有与C-R格式中分别用双二次多项式逼近相同的收敛阶.     

关于Ciarlet-Raviart混合有限元法的一个注记

本文推广解双调合方程的Ciarlet-Raviart混合有限元方案：用二次元逼近流函数φ．一次元逼近涡度-Δφ．在拟一致三角形剖分的条件下,证明了推广方案具有φ和-Δφ都用二次元逼近的标准Ciarlet－Raviart方案同样的精度阶．     

球Couette流的数值模拟

该文利用谱方法对同心旋转球间轴对称Couette流进行数值模拟.给出Navier Stokes方程的流函数涡度形式,利用Stokes流把边界条件齐次化, 选取Stokes算子的特征函数做为逼近子空间的基函数,对同心旋转球间轴对称Couette流进行谱逼近     

不可压磁流体力学方程组的高精度数值解法

为数值求解低雷诺数下不可压流体在电磁场作用下的流动,提出一种四阶紧致差分方法.由二维原始变量的MHD方程组出发,推导出具有较少未知量的电流密度-涡量-流函数形式MHD方程组.建立了求解二维非定常不可压MHD方程组的电流密度-涡量-流函数形式的四阶精度紧致差分格式.为验证本文提出的高精度紧致差分方法的精确性和可靠性,对有...     

求解流函数涡度方程的特征线—差分方法

本文讨论了用特征线方法与有限差分方法相结合的数值方法(特征线—差分方法)求解流函数涡度形式的Navier-Stokes方程的问题.证明了该方法的收敛性,给出了数值例子.     

对Taylor—Gelerkin有限元法的一点改进和它的应用

本文针对Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ有限元法的两个基本假设进行讨论。改进了原假设，仅以一个假设作为出发点。得到了广义的有限元离散公式。对具体流函数－涡量方程的求解进行了改进的Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ有限元分析。提出了组合式的求解方法，使求解过程更为合理。算例计算表明，该方法的效果是很好的。     

近壁湍流脉动的概率分布函数

采用大涡模拟方法,模拟槽道湍流,获得了不同雷诺数情形下的槽道流大涡模拟数据库．在此基础上,获得了流向和垂向脉动速度的概率分布函数,并运用假设检验,分析了其与正态分布的定量差别．进一步计算了流向和垂向脉动速度的偏斜度、平坦度,讨论了二者在粘性子层、过渡区和对数律区的变化．同时,讨论了粘性子层、过渡区和对数律区流向和垂向脉动速度概率分布函数的特点及其与湍流猝发的高速流下扫和低速流喷发事件的关系．最后,分析了雷诺数对流向、垂向脉动速度分布的影响．     

对Taylor-Galerkin有限元法的一点改进和它的应用

本文针对Taylor-Galerkin有限元法的两个基本假设进行讨论.改进了原假设,仅以一个假设作为出发点,得到了广义的有限元离散公式.对具体流函数—涡量方程的求解进行了改进的Taylor-Galerkin有限元分析.提出了组合式的求解方法,使求解过程更为合理.算例计算表明,该方法的效果是很好的.     

液固两相圆柱绕流尾迹内颗粒扩散分布的离散涡数值研究

基于离散涡方法求得的非定常水流场和颗粒的Lagrange运动方程,数值模拟了稀疏液固两相圆柱绕流尾迹内颗粒的扩散分布．获得了流动的涡谱与3种不同St数颗粒（St=0.25,1.0,40）在流场中的分布．通过引入扩散函数来定量表示颗粒在流场中的纵向扩散强度,并计算得到了不同St数颗粒的扩散函数随时间的变化．数值结果揭示出了液固两相圆柱绕流尾迹中的颗粒扩散分布与颗粒的St数和尾涡结构密切相关:1) 中小St数(St=0.25～4.0)颗粒在运动过程中不能进入涡核区,而在旋涡结构的外沿聚集,且颗粒的St数愈大,其越远离涡核区域;2） 在圆柱绕流尾迹区域内,中小St数(St=0.25～4.0)颗粒的纵向扩散强度随其St数的增大而减小．     

环形空腔内自然对流问题的Galerkin方法

本文讨论了环形空腔内自然对流问题所满足的Boussinesq方程组——关于涡度ζ、流函数φ及温度θ的椭圆-抛物非线性耦合方程组 用Galerkin方法对其进行了数值分析,得到了Galerkin逼近(含半离散和全离散)的最优先验误差估计。     

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     

Nonconforming mixed finite elements for linear elasticity on simplicial grids

This paper introduces a new family of nonconforming mixed finite elements for solving the linear elasticity equations on simplicial grids. Besides, this paper describes the construction of the lowest order basis functions. The construction only involves simple computations due to the new explicit stress shape function spaces and the procedure applies for high order cases. Numerical experiments for four benchmark problems in mechanics indicate the robust locking‐free behavior and show that the lowest order nonconforming mixed method leads to smaller stress errors than the first and second order standard Galerkin methods for the nearly incompressible case.     

OPTIMAL ERROR ESTIMATES FOR NEDELEC EDGE ELEMENTS FOR TIME-HARMONIC MAXWELL'S EQUATIONS

In this paper, we obtain optimal error estimates in both L^2-norm and H（curl）-norm for the Nedelec edge finite element approximation of the time-harmonic Maxwell＇s equations on a general Lipschitz domain discretized on quasi-uniform meshes. One key to our proof is to transform the L^2 error estimates into the L^2 estimate of a discrete divergence-free function which belongs to the edge finite element spaces, and then use the approximation of the discrete divergence-free function by the continuous divergence-free function and a duality argument for the continuous divergence-free function. For Nedelec＇s second type elements, we present an optimal convergence estimate which improves the best results available in the literature.     

稳定化有限元方法中逆估计常数的确定

0．弓I言稳定化有限元方法［2］［4］［12］［15］ro［18］在固体和流体力学的数值计算中构造有效的格式发挥着很大的作用．从理论分析的角度看，该方法（Galerkin一局部最小二乘方法）是完备的、确定的，但是在实际计算中稳定化参数。E（0，CI）的如何选取直接影响到逼近解的质量．数值实验【9川叫'川'到表明。取得太小会造成。伪l。压力模式rI．因此，对of的估计是一个值得注意的重要问题．文［8］虽然估训一了一些逆常数，但其未能给出确定逆估计常数的一般公式，而且技巧性太强，过于依赖区域剖分的性质，使得逆常数的计算复杂化，不…     

一类抛物型方程的强超逼近性质

本研究一类一维抛物型方程初边值问题半离散有限元解，对iku-kh分别得到函数及导数O(h^k 2.5)阶及O（h^k 2）阶的强超逼近估计。     

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.     

受限制多项式插值及在构造形函数空间中的应用

1引言 G.Strang指出：有限元法的新思想在于试控函数的选择，目标是选择这样的分片多项式，它们被少数几个方面的节点值确定，并仍具有我们需要的连续性和逼近度。受限制多项式插值空间就是这样一类空间，在P.G.Ciarlet~[4]的书中有较多的介绍，采用的方法是通过约束条件来决定试验空间，但正如[1]中指出的，这样约束条件欠直观，且容易产生一些不确     

Approximation by quadrilateral finite elements

We consider the approximation properties of finite element spaces on quadrilateral meshes. The finite element spaces are constructed starting with a given finite dimensional space of functions on a square reference element, which is then transformed to a space of functions on each convex quadrilateral element via a bilinear isomorphism of the square onto the element. It is known that for affine isomorphisms, a necessary and sufficient condition for approximation of order in and order in is that the given space of functions on the reference element contain all polynomial functions of total degree at most . In the case of bilinear isomorphisms, it is known that the same estimates hold if the function space contains all polynomial functions of separate degree . We show, by means of a counterexample, that this latter condition is also necessary. As applications, we demonstrate degradation of the convergence order on quadrilateral meshes as compared to rectangular meshes for serendipity finite elements and for various mixed and nonconforming finite elements.

     

Mapped Finite Elements on Hexahedra. Necessary and Sufficient Conditions for Optimal Interpolation Errors

This paper considers finite elements which are defined on hexahedral cells via a reference transformation which is in general trilinear. For affine reference mappings, the necessary and sufficient condition for an interpolation order O(h ^k+1) in the L ²-norm and O(h ^k ) in the H ¹-norm is that the finite dimensional function space on the reference cell contains all polynomials of degree less than or equal to k. The situation changes in the case of a general trilinear reference transformation. We will show that on general meshes the necessary and sufficient condition for an optimal order for the interpolation error is that the space of polynomials of degree less than or equal to k in each variable separately is contained in the function space on the reference cell. Furthermore, we will show that this condition can be weakened on special families of meshes. These families which are obtained by applying usual refinement techniques can be characterized by the asymptotic behaviour of the semi-norms of the reference mapping.