非定常对流扩散问题差分-流线扩散法的线性三角元解的最优误差估计

讨论了差分-流线扩散法(FDSD)求解线性对流占优扩散问题解的精度,利用插值后处理技术,使该格式解的空间精间达到最优.     

非线性对流扩散问题的差分-流线扩散法

1．引言流线扩散法（简称SD方法）是由Huzhes和Brooks在1980年前后提出的一种数值求解对流占优扩散问题的新型有限元算法．随后，Johnson和N8vert将SD方法推广到发展型对流扩散问题（［1］，［2］，［3］）．熟知，对于对流扩散问题，标准有限元法虽具有高阶精度，但常产生数值振荡；古典人工粘性Galerkin法更具有较好的稳定性，但仅具有一阶精度．而（SD方法兼具良好的数值稳定性和高阶精度，因此得到了越来越多的重视，对于发展型对流扩散问题，传统的SD方法均采用时空有限元．这样做，虽然可使时间和空间方向上的精度很好的协调起…     

一类非定常对流占优扩散问题差分-流线扩散法的后处理

讨论了一类非定常对流占优扩散方程的差分-流线扩散格式(FDSD),利用插值后处理技术,提高了特殊网格下该FDSD格式在双线性元空间的精度,从而按L∞(L2(Ω) 模达到最优．     

二维发展型对流占优扩散方程的FD-SD法的后验误差估计

引言 对流占优扩散问题是流体力学中一个典型的模型问题,对其数值求解始终是众多学者相当关心的课题．［１１］中指出,即使对于线性问题,通常其解在外流边界附近也会产生剧烈变化．倘若在内流边界上所给出的边值函数存在不连续点时,则在沿过此不连续点的特征线（流线）附近会出现断层．因此在数值求解对流占优扩散问题时,尽管标准有限元法具有高阶精度,但常产生数值剧烈振荡Ｓ而古典人工粘性Ｇａｌｅｒｋｉｎ法虽具有较好的稳定性,但仅具有一阶精度．流线扩散法（Ｓｔｒｅａｍｌｉｎｅ　 Ｄｉｆｆｕｓｉｏｎ Ｍｅｔｈｏｄ,简称 ＳＤ…     

对流扩散问题流线扩散法的最优误差估计

在本文中,我们考虑对流占优扩散问题流线扩散双线性有限元方法。原先的文献在ε≤h~2的条件下,得到了L~2-模最优误差估计,而本文则在ε≤h的条件下得到了相同估计。     

对流-扩散问题的6节点三角形单元流线迎风有限元法和自适应网格重分技术

提出了使用6节点三角形单元的流线迎风有限元法．该方法沿局部流线，直接用于输运控制方程的对流项．采用多个对流一扩散实例来评价该方法的有效性，结果显示该方法是单调的，并且不产生任何振荡．另外，自适应网格技术和该方法相结合后，进一步提高了解的精度，又减少了计算时间和对计算机内存的需求．     

对流扩散问题的流线扩散有限元分析

对流扩散问题的流线扩散有限元分析孙澈，张阳，魏强（南开大学）ＴＨＥＳＴＲＥＡＭＬＩＮＥ－ＤＩＦＦＵＳＩＯＮＦ．Ｅ．ＡＮＡＬＹＳＩＳＦＯＲＣＯＮＶＥＣＴＩＯＮＤＩＦＦＵＳＩＯＮＰＲＯＢＬＥＭＳ￥ＳｕｎＣｈｅ；ＺｈａｎｇＹａｎｇ；ＷｅｉＱｉａｎｇ（Ｎａｎ...     

非线性对流扩散问题的预测校正型FDSD格式

1 引言 流线扩散法(简称SD方法)是由Hughes和Brooks提出,并经Johnson等人发展的一种数值求解对流占优扩散问题的一种有效的数值方法。然而,传统的SD方法利用时-空有限元求解发展型问题,导致对高维问题工作量过于庞大,其编程实现较复杂,对非线性问题也不便进行线性化处理,为此,孙澈提出了仅对空间域作有限元离散,而对时间域作差分离散的差分-流线扩散法(以下简称FDSD方法),主要讨论了拟线性     

对流扩散问题的交替方向差分-流线扩散格式

1.引 言 差分-流线扩散法（Finite Difference-Streamline Diffusion Method,简称FDSD方法）于1998年由文[1]提出并对线性对流占优扩散问题给出分析,随后文[2],[3]就非线性问题的FDSD格式及FDSD预测-校正格式,分别作出了分析,文[4]讨论了FDSD方法的后验估计及自适应技术,[5],[6]则分别讨论了FDSD方法的某些重要应用.与基于时-空有限元的传统流线扩散法相比,FDSD方法的计算工作量已有成数量级的减少,且较易于推广到非线性问题,然而,对于高维问题,在每一时间层,仍然需要求解一大型线性或非线性方程组,工作量仍然很大.参照J.Douglas与T.Dupont关于抛物问题交替方向     

对流扩散方程的稳定化流线扩散法最小二乘非协调有限元分析

将最小二乘法和稳定化的流线扩散法相结合,研究了对流扩散方程的非协调有限元格式,用矩形EQ_1~(rot)元和零阶R-T元分别来逼近位移和应力,利用单元本身的特殊性质,证明了离散格式解的存在惟一性,得到了位移H~1-模和应力H(div)-模的最优误差估计.     

一种用于索结构分析的悬链线单元

根据弹性悬链线的理论解析解推导出适于索结构有限元分析的悬链线单元.与常用的三节点、五节点曲线单元相比,采用该单元编制的软件具有输入数据少、计算机时省、计算精度高的特点.     

非正则条件下类Wilson元的构造及其应用

本文在非正则性条件下，研究了窄四边形上的类Wilson元。通过参考元上类Wilson元的构造，证明了由此产生的有限元对任意窄四边形剖分通过Irons分片检查，得到了二阶问题的误差估计。结果表明，该单元的收敛性质与Wilson元的类似。     

GLOBAL FINITE ELEMENT NONLINEAR GALERKIN METHOD FOR THE PENALIZED NAVIER-STOKES EQUATIONS GLOBAL FINITE ELEMENT NONLINEAR GALERKIN METHOD FOR THE PENALIZED NAVIER-STOKES EQUATIONS

1. IntroductionIn the numerical simulation of the Navier-Stokes equations one encounters three seriousdifficulties in the case of large Reynolds numbers f the treatment of the incomPressibility con-dition divu = 0, the treatment of the noIilinear terms and the large time integration. For thetreatment of the incoInPressibility condition, one use the penalty method in the case of finiteelemellts [1--2l and for the treatmen of the noulinar terms and the large tfor integration, oneuse the nonlin…     

Rail-bridge coupling element of unequal lengths for analysing train-track-bridge interaction systems

This paper presents a rail-bridge coupling element of unequal lengths, in which the length of a bridge element is longer than that of a rail element, to investigate the dynamic problem of train-track-bridge interaction systems. The equation of motion in matrix form is given for a train-track-bridge interaction system with the proposed element. The first two numerical examples with two types of bridge models are chosen to illustrate the application of the proposed element. The results show that, for the same length of rail element, (1) the dynamic responses of train, track and bridge obtained by the proposed element are almost identical to those obtained by the rail-bridge coupling element of equal length, and (2) compared with the rail-bridge coupling element of equal length, the proposed element can help to save computer time. Furthermore, the influence of the length of rail element on the dynamic responses of rail is significant. However, the influence of the length of rail element on the dynamic responses of bridge is insignificant. Therefore, the proposed element with a shorter rail element and a longer bridge element may be adopted to study the dynamic responses of a train-track-bridge interaction system. The last numerical example is to investigate the effects of two types of track models on the dynamic responses of vehicle, rail and bridge. The results show that: (1) there are differences of the dynamic responses of vehicle, rail and bridge based on the single-layer and double-layer track models, (2) the maximum differences increase with the increase of the mass of sleeper, (3) the double-layer track model is more accurate.     

三维Poisson方程特征值的四种有限元解及比较

应用三维EQ1rot元、三维Crouzeix-Raviart元、八节点等参数元、四面体线性元计算三维Poisson方程的近似特征值.计算结果表明:三维EQ1rot元和三维Crouzeix-Raviart元特征值下逼近准确特征值,八节点等参数元、四面体线性元特征值上逼近准确特征值,三维EQr1ot元和三维Crouzeix-Raviart元外推特征值下逼近准确特征值.计算结果还表明三维Crouzeix-Raviart元是一种计算效率较高的非协调元.     

Infinite Groups with an Anticentral Element

An element of a group is called anticentral if the conjugacy class of that element is equal to the coset of the commutator subgroup containing that element. A group is called Camina group if every element outside the commutator subgroup is anticentral. In this paper, we investigate the structure of locally finite groups with an anticentral element. Moreover, we construct some non-periodic examples of Camina groups, which are not locally solvable.     

半导体器件瞬态模拟的对称正定混合元方法

提出具有对称正定特性的混合元格式求解非稳态半导体器件瞬态模拟问题。提出一个最小二乘混合元方法、一个新的具有分裂和对称正定性质的混合元格式和一个解经典混合元方程的对称正定失窃工格式求解电场位势和电场强度方程;提出一个最小二乘混合元格式求解关于电子与空穴浓度的非稳态对流扩散方程,浓度函数和流函数被同时求解;采用标准的有限元方法求解热传导方程。建立了误差分析理论。     

一族新的三维体元——混合刚度有限元

对于三维连续介质的有限元分析,一个通用的二次有限单元体是所谓20节点等参元.尽管这个元素已富有成效地被普遍应用,但是,它需要的节点自由度太多与其达到的二次多项式的逼近精度却十分不相称,显得计算效率很低.基于混合刚度有限元法的一     

A postprocessed flux conserving finite element solution

We propose a local postprocessing method to get a new finite element solution whose flux is conservative element‐wise. First, we use the so‐called polynomial preserving recovery (postprocessing) technique to obtain a higher order flux which is continuous across the element boundary. Then, we use special bubble functions, which have a nonzero flux only on one face‐edge or face‐triangle of each element, to correct the finite element solution element by element, guided by the above super‐convergent flux and the element mass. The new finite element solution preserves mass element‐wise and retains the quasioptimality in approximation. The method produces a conservative flux, of high‐order accuracy, satisfying the constitutive law. Numerical tests in 2D and 3D are presented.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1859–1883, 2017     

Fourteen-Node Mixed Stiffness Element and Its Computational Comparisons with Twenty-Node Isoparametric Element

In this paper, a family of 3-dimensional elements different from isoparametric serendipity is developed according to the variational principle and the convergence criteria of the mixed stiffness finite element method. For the new family, which is named mixed stiffness elements, the number of nodes on the quadratic element is not 20 but 14. Theoretical analysis and various computational comparisons have found the mixed stiffness element superior over the isoparametric serendipity element, especially a substantial improvement in computational efficiency can be achieved by replacing the 20 node-isoparametric element with the 14-node mixed stiffness element.