一维双曲守恒方程组的Taylor-Galerkin有限元方法

１．引言 在文章［８］中,利用双曲守恒律的Ｈａｍｉｌｔｏｎ－Ｊａｃｏｂｉ方程形式,应用 Ｇａｌｅｒｋｉｎ有限元给出了求解一维双曲守恒律的计算方法．不同于间断有限元方法［２］、［３］和 Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ有限元方法［１］求解双曲守恒律,文章［８］采用连续 Ｇａｌｅｒｋｉｎ有限元求解双曲守恒律． 在文章［８］中,对差分方法和有限元方法求解双曲守恒律作了较为详细的讨论．同时在文章［８］中,采用积分变换,将双曲守恒律方程变成 Ｈａｍｉｌｔｏｎ－Ｊａｃｏｂｉ方程形式．由于 Ｈａｍｉｌｔｏｎ－Ｊａｃｏ…     

基于Rosenbrock型指数积分的一维间断Galerkin有限元方法

提出基于Rosenbrock型指数积分的一维间断Galerkin有限元方法.该方法在空间上使用间断有限元方法离散,在时间上采用Rosenbrock型指数积分方法.这样不仅可以保持空间离散上的高精度,而且继承了指数时间积分方法具有显式大步长时间推进的优点.数值试验的结果表明,对于一维双曲守恒律问题,这种方法是一种有效的数值算法.     

三维高速冲击动力有限元滑移面算法

本文给出三维高速冲击动力有限元滑移面算法.该算法不但能保证单元结点的动量守恒、动量矩守恒,而且由实例计算表明该算法在处理高速穿、破甲过程中是稳定和有效的.     

一个求解多维守恒律方程组的二阶显式有限元格式

１．引言 近年来,在非结构网格上求解双曲型守恒律的数值方法引起了较为广泛的关注,出现了有限体积方法［１］,间断 Ｇａｌｅｒｋｉｎ方法 ［２］,流线扩散方法［３］,以及 ＮＮＤ格式 ［４］等．我们在［６,７］中提出了一种求解双曲型守恒律方程式的有限元方法,它是在一个求解对流扩散问题的有限元方法 ［５］的基础上发展起来的．它是一个显式有限元方法,因此计算量很小．在这个方法中,我们将任意维的问题归结为在单元棱边上的一维计算,引入了积分因子,因此在单元内部可以容纳边界层．这样,它特别适合于对流占优问题以及双曲…     

一阶双曲问题间断有限元的后验误差分析

一阶双曲问题的有限元后验误差估计至今没有得到很好的解决.本文对d维区域上一阶双曲问题的k次间断有限元逼近提出了一种新的后验误差分析方法, 进而建立了间断有限元解在DG范数下(强于L₂范数)基于误差余量型的后验误差估计. 数值计算验证了本文理论分析的有效性. 本文方法也适用于其他变分问题有限元逼近的后验误差分析.     

二维非线性Schr?dinger方程的两类局部守恒算法

本文针对二维非线性Schr?dinger方程,提出两类局部守恒算法.不需要考虑边界条件,即可保持任意时空区域上相应的局部能量守恒律和局部动量守恒律.在合适的边界条件下,它们能自然地保持电荷、全局能量或全局动量守恒律.本文同时对算法进行了守恒分析和误差分析.在数值实验部分,本文构造了类似的多辛Preissman算法进行比较,数值结果验证了其长时间计算的优势.     

对流反应扩散方程的稳定化时间间断时空有限元解的误差估计

在流线迎风Petrov-Galerkin（SUPG）稳定化有限元数值格式的基础上,结合时间方向的变分离散,构造对流反应扩散方程的稳定化时间间断时空有限元格式.该类格式在工程上有一些数值模拟应用,但相关文献没有看到类似数值格式的理论证明.本文以Radau点为节点,构造时间方向的Lagrange插值多项式,证明了稳定化有限元解的稳定性,时间最大模、空间L²（Ω）-模误差估计.文中利用插值多项式和有限元方法相结合的技巧,解耦时空变量,去掉了时空网格的限制条件,提供了时间间断稳定化时空有限元方法的理论证明思路,克服了因时空变量统一导致的实际计算时的复杂性.     

半线性抛物型积分微分方程的间断时空有限元方法

本文讨论了一类半线性抛物型积分微分方程的间断时空有限元方法．利用有限元和有限差分方法相结合的技巧,在时间离散区间内,利用Radau点处Lagrange插值多项式的特性,去掉间断时空有限元的传统证明过程中对时空网格的限制条件,并给出了时间最大模、空间L_2模,即L_∞(L_2)模的误差估计．     

对流扩散方程的间断时空有限元方法的误差估计

研究了一类线性对流扩散方程的间断时空有限元方法,即空间连续,时间允许间断的时空有限元方法.将有限元方法和有限差分方法相结合,在每一时间层上充分利用Lagrange插值多项式在Radau点处的特性,给出了有限元解的最优阶L∞(L2)模误差估计.     

多介质流体计算的Lax-Friedrichs格式

针对当前多介质流体计算中出现的速度和压强在介质界面处产生伪振荡的问题，我们设计了一种基于非平衡态的Lax-Friedrichs格式．我们构造的算法保证了多个守恒性质：总质量和分质量守恒，总动量和总能量守恒，它还可以保证分质量非负．更重要的是它消除了速度和压强在介质界面处的伪振荡．数值例子表明这一算法是有效的．     

A framework for obtaining guaranteed error bounds for finite element approximations

We give an overview of our recent progress in developing a framework for the derivation of fully computable guaranteed posteriori error bounds for finite element approximation including conforming, non-conforming, mixed and discontinuous finite element schemes. Whilst the details of the actual estimator are rather different for each particular scheme, there is nonetheless a common underlying structure at work in all cases. We aim to illustrate this structure by treating conforming, non-conforming and discontinuous finite element schemes in a single framework. In taking a rather general viewpoint, some of the finer details of the analysis that rely on the specific properties of each particular scheme are obscured but, in return, we hope to allow the reader to ‘see the wood despite the trees’.     

A coupled continuous-discontinuous FEM approach for convection diffusion equations

In this article, we introduce a coupled approach of local discontinuous Galerkin and standard finite element method for solving convection diffusion problems. The whole domain is divided into two disjoint subdomains. The discontinuous Galerkin method is adopted in the subdomain where the solution varies rapidly, while the standard finite element method is used in the other subdomain due to its lower computational cost. The stability and a priori error estimate are established. We prove that the coupled method has O((ε1 / 2 + h 1 / 2 )h k ) convergence rate in an associated norm, where ε is the diffusion coefficient, h is the mesh size and k is the degree of polynomial. The numerical results verify our theoretical results. Moreover, 2k-order superconvergence of the numerical traces at the nodes, and the optimal convergence of the errors under L 2 norm are observed numerically on the uniform mesh. The numerical results also indicate that the coupled method has the same convergence order and almost the same errors as the purely LDG method.     

MIXED FINITE ELEMENT METHODS FOR FRACTIONAL NAVIER-STOKES EQUATIONS

This paper gives the detailed numerical analysis of mixed finite element method for fractional Navier-Stokes equations.The proposed method is based on the mixed finite element method in space and a finite difference scheme in time.The stability analyses of semi-discretization scheme and fully discrete scheme are discussed in detail.Furthermore,We give the convergence analysis for both semidiscrete and flly discrete schemes and then prove that the numerical solution converges the exact one with order O(h2+k),where h and k:respectively denote the space step size and the time step size.Finally,numerical examples are presented to demonstrate the effectiveness of our numerical methods.     

非稳态奇异系数微分方程的时间间断时空有限元方法

针对一类二维单轴奇异系数非稳态问题构造了一种时间间断时空有限元格式,利用以Radau点为节点的Lagrange插值多项式的特性,结合有限差分法和有限元法的技巧证明了格式的稳定性和有限元解的时间最大模、空间加权L²（?）-模误差估计.最后列举了一些数值试验结果,验证了理论结果和格式的可行性.     

Material-Force-Based Refinement Indicators in Adaptive Strategies

A material-force-based refinement indicator for adaptive finite element strategies for finite elasto-plasticity is proposed. Starting from the local format of the spatial balance of linear momentum, a dual material counterpart in terms of Eshelby's energy-momentum tensor is derived. For inelastic problems, this material balance law depends on the material gradient of the internal variables. In a global format the material balance equation coincides with an equilibrium condition of material forces. For a homogeneous body, this condition corresponds to vanishing discrete material nodal forces. However, due to insufficient discretization, spurious material forces occur at the interior nodes of the finite element mesh. These nodal forces are used as an indicator for mesh refinement. Assigning the ideas of elasticity, where material forces have a clear energetic meaning, the magnitude of the discrete nodal forces is used to define a relative global criterion governing the decision on mesh refinement. Following the same reasoning, in a second step a criterion on the element level is computed which governs the local h-adaptive refinement procedure. The mesh refinement is documented for a representative numerical example of finite elasto-plasticity. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A monotone finite element scheme for convection-diffusion equations

A simple technique is given in this paper for the construction and analysis of a class of finite element discretizations for convection-diffusion problems in any spatial dimension by properly averaging the PDE coefficients on element edges. The resulting finite element stiffness matrix is an -matrix under some mild assumption for the underlying (generally unstructured) finite element grids. As a consequence the proposed edge-averaged finite element scheme is particularly interesting for the discretization of convection dominated problems. This scheme admits a simple variational formulation, it is easy to analyze, and it is also suitable for problems with a relatively smooth flux variable. Some simple numerical examples are given to demonstrate its effectiveness for convection dominated problems.

     

On the use of physical spline finite element method for acoustic scattering

This paper presents the comparison of physical spline finite element method (PSFEM), in which differential equations are incorporated into interpolations of basic elements, with least-squares finite element method (LSFEM) and mixed Galerkin finite element method (MGFEM) on the numerical solution of one dimensional Helmholtz equation applied to an acoustic scattering problem. Firstly, all three methods are explained in detail and then it is shown that PSFEM reaches higher precision in a shorter time with fewer nodes than the other methods. It is also observed that this method is well suited for high frequency acoustic problems. Consequently, the results of PSFEM point out better efficiency in terms of number of unknowns and accuracy level.     

高波数问题的超收敛性

 本文考虑求解Helmholtz方程的有限元方法的超逼近性质以及基于PPR后处理方法的超收敛性质.我们首先给出了矩形网格上的p-次元在收敛条件k（kh）^2p+1≤C₀下的有限元解和基于Lobatto点的有限元插值之间的超逼近以及重构的有限元梯度和精确解之间的超收敛分析.然后我们给出了四边形网格上的线性有限元方法的分析.这些估计都给出了与波数k和网格尺寸h的依赖关系.同时我们回顾了三角形网格上的线性有限元的超收敛结果.最后我们给出了数值实验并且结合Richardson外推进一步减少了误差.     

求解流固耦合问题的一种四步分裂有限元算法

基于arbitrary Lagrangian Eulerian (ALE) 有限元方法,发展了一种求解流固耦合问题的弱耦合算法．将半隐式四步分裂有限元格式推广至求解ALE描述下的Navier-Stokes（N-S）方程,并在动量方程中引入迎风流线（streamline upwind/Petrov-Galerkin, SUPG）稳定项以消除对流引发的速度场数值振荡;采用Newmark-β法对结构方程进行时间离散;运用经典的Galerkin有限元法求解修正的Laplace方程以实现网格更新,每个计算步施加网格总变形量防止结构长时间、大位移运动时的网格质量恶化．运用上述算法对弹性支撑刚性圆柱体的流致振动问题进行了数值模拟,计算结果与已有结果相吻合,初步验证了该算法的正确性和有效性．