Local and Parallel Finite Element Algorithms Based on Two-Grid Discretizations for Nonlinear Problems

In this paper, some local and parallel discretizations and adaptive finite element algorithms are proposed and analyzed for nonlinear elliptic boundary value problems in both two and three dimensions. The main technique is to use a standard finite element discretization on a coarse grid to approximate low frequencies and then to apply some linearized discretization on a fine grid to correct the resulted residual (which contains mostly high frequencies) by some local/parallel procedures. The theoretical tools for analyzing these methods are some local a priori and a posteriori error estimates for finite element solutions on general shape-regular grids that are also obtained in this paper.     

基于二重网格的定常Navier-Stokes方程的局部和并行有限元算法

对二维定常的不可压缩的Navier-Stokes方程的局部和并行算法进行了研究．给出的算法是多重网格和区域分解相结合的算法,它是基于两个有限元空间:粗网格上的函数空间和子区域的细网格上的函数空间．局部算法是在粗网格上求一个非线性问题,然后在细网格上求一个线性问题,并舍掉内部边界附近的误差相对较大的解．最后,基于局部算法,通过有重叠的区域分解而构造了并行算法,并且做了算法的误差分析,得到了比标准有限元方法更好的误差估计,也对算法做了数值试验,数值结果通过比较验证了本算法的高效性和合理性．     

Finite Volume Methods for Eigenvalue Problems

Based on a linear finite element space, in this paper, two symmetric finite volume schemes are proposed for self-adjoint elliptic boundary eigenvalue problems. Both convergence and superconvergence are discussed.This revised version was published online in October 2005 with corrections to the Cover Date.     

非协调元特征值渐近下界

利用有限元收敛速度下界的结果获得某些非协调元方法新的Aubin-Nitsche估计形式,然后再结合非协调元特征值的展开式获得不需要额外条件下非协调元特征值渐近下界的结果.     

特征值问题的Lagrange型各向异性有限元方法

在各向异性网格下首先研究了二阶椭圆特征值问题算子谱逼近的若干抽象结果.然后将这些结果具体应用于线性和双线性Lagrange型协调有限元,得到了与传统有限元网格剖分下相同的最优误差估计,从而拓宽了已有的成果.     

解椭圆边值问题的有限元并行算法

本文研究椭圆边值问题有限元方程的求解，在对限元基函数一种特定的“红黑”排序基础上，构造出具有异步并行计算结构的迭代算法，并证明了算法的收敛性。     

一种迭代格式的有限元并行算法*

本文提出了一种求解有限元方程的迭代格式的并行算法.该方法在线性代数方程迭代解法的基础上,引进并行运算步骤;并且运用加权残数方法,通过选择适当的权函数,推导了该并行算法的有限元基本格式.该方法在西安交通大学BLXSI-6400并行计算机上程序实现.计算结果表明它能有效地提高运算速度,减少计算时间,是一种有效的求解大型结构有限元方程的并行算法.     

最优化问题的并行算法

本文对求解非线性最优化问题的几种主要并行思想，即按变量分裂的并行算法，函数值、梯度值的并行计算，计算步骤并行的算法等，作了简要的综述，并介绍了近几年在这方面取得的进展．     

特征值问题的三角形二次元的展开式及其外推

我们考虑利用三角形二次元来求解特征值问题,并给出特征值的误差展开式,以此为基础进行外推获得高精度.     

关于特征值问题有限元方法研究的一个说法

受林群的微积分哲学公式((相对真理)/(绝对真理)=0.9)的启示,总结了近些年来关于特征值问题有限元方法的研究,并发现其背后同样隐藏着该哲学公式,换句话说,所追求的是特征值问题的有限元数值解和真解的零距离,其实就是追求真解的过程,要经多道(即0.9,0.99,0.999,…),再将比例,即(数值解)/(真解),提到1.     

三维有限元本征值的外推

通过二维和三维积分恒等式,探讨泊松方程本征值问题三角线元和四面体线元Richardson外推的可行性.理论分析表明,如果剖分为均匀一致和拟一致,外推均可将解的精度提高二阶.     

弹性接触问题参数变分原理的有限元并行算法*

本文基于弹性接触问题的参数变分原理的有限元解法,利用并行计算机的特性和并行处理结构,建立了相应的并行算法.该算法从刚度阵的生成和组集,静凝聚过程,求应力过程等多方面实现了并行化.该算法在西安交通大学ELXSI-6400并行计算机上程序实现,计算结果表明能有效地节省计算时间,是一种分析接触问题的有效的并行算法.     

Parallel Algorithms for Orthotropic Problems

Finite element meshes and node-numberings suitable for parallel solution with equally loaded processors are presented for linear orthotropic elliptic partial differential equations. These problems are of great importance, for instance in the oil and airfoil industries. The linear systems of equations are solved by the conjugate gradient method preconditioned by modified incomplete factorization, MIC. The basic method presented, is based on fronts of uncoupled nodes and unlike earlier methods it has the advantage of no requirement of a specific orientation of the mesh. This method is however, in general, restricted to small degree of anisotropy in the differential equation. Another method, which does not suffer from this limitation, uses rotation of the differential equation and spectral equivalence. The rotation is made in such a way that in the new co-ordinate system, the basic method is applicable. The spectral equivalence property is used for estimation of the condition number of the preconditioned system. Both methods are suitable for implementation on parallel computers. The computer architecture could be single instruction multiple data (SIMD) as well as multiple instruction multiple data (MIMD) with shared or distributed memory. Implementation of the basic method on a shared memory parallel computer shows a significant improvement by use of the MIC method compared with the diagonal scaling preconditioning method.     

Parallel Synchronous and Asynchronous Space-Decomposition Algorithms for Large-Scale Minimization Problems

Three parallel space-decomposition minimization (PSDM) algorithms, based on the parallel variable transformation (PVT) and the parallel gradient distribution (PGD) algorithms (O.L. Mangasarian, SIMA Journal on Control and Optimization, vol. 33, no. 6, pp. 1916–1925.), are presented for solving convex or nonconvex unconstrained minimization problems. The PSDM algorithms decompose the variable space into subspaces and distribute these decomposed subproblems among parallel processors. It is shown that if all decomposed subproblems are uncoupled of each other, they can be solved independently. Otherwise, the parallel algorithms presented in this paper can be used. Numerical experiments show that these parallel algorithms can save processor time, particularly for medium and large-scale problems. Up to six parallel processors are connected by Ethernet networks to solve four large-scale minimization problems. The results are compared with those obtained by using sequential algorithms run on a single processor. An application of the PSDM algorithms to the training of multilayer Adaptive Linear Neurons (Madaline) and a new parallel architecture for such parallel training are also presented.     

非定常对流扩散方程的局部和并行有限元算法

对基于两重网格的非定常对流扩散方程的局部和并行有限元算法进行了研究.算法的理论依据是两重网格的思想,解的低频分量可以用一个整体的粗网格空间来逼近,高频分量可以用局部和并行的细网格空间来逼近.因此,这种局部和并行算法仅仅涉及一个粗网格上的整体逼近和细网格上的局部校正.得到了算法的误差估计,一些数值例子验证了算法的有效性.     

各向异性网格下特征值问题的有限元分析

主要目的是在各向异性网格下研究二阶椭圆特征值问题的两类非协调有限元—类Wilson矩形元和Carey三角形元—的收敛性分析.通过新的技巧和方法,得到了与传统有限元网格剖分下相同的特征对的最优误差估计.推广了已有的结果.     

Parallel Finite Element Software Development and Performance Analysis in an Object-Oriented Programming Framework

The use of object-oriented programming techniques in the development of parallel, finite element analysis software enhances code reuse and increases efficiency during application development. In this paper, an object-oriented programming framework developed by the authors is utilized in the implementation of parallel finite element software for modeling of the resin transfer molding manufacturing process. The motivation for choosing the resin transfer molding finite element application and implementing it with the object-oriented framework is that it was originally developed and parallelized in a functional programming paradigm thus offering the possibility of direct comparisons. Discussion of the software development effort and performance results are presented and analyzed.Mathematics Subject Classifications (2000) 65M60, 65Y05.     

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

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

一类优化控制问题的有限元解的误差估计和超收敛

对积分微分方程的优化控制问题进行了介绍.讨论了积分微分方程的优化控制问题的混合有限元逼近,给出了优化控制问题的有限元逼近解的误差估计和超收敛性质.     

Extrapolation of Mixed Finite Element Approximations for the Maxwell Eigenvalue Problem

In this paper,a general method to derive asymptotic error expansion formulas for the mixed finite element approximations of the Maxwell eigenvalue problem is established.Abstract lemmas for the error of the eigenvalue approximations are obtained.Based on the asymptotic error expansion formulas,the Richardson extrapolation method is employed to improve the accuracy of the approximations for the eigenvalues of the Maxwell system from θ(h2) to θ(h4) when applying the lowest order Nédé1ec mixed finite element and a nonconforming mixed finite element.To our best knowledge,this is the first superconvergence result of the Maxwell eigenvalue problem by the extrapolation of the mixed finite element approximation.Numerical experiments are provided to demonstrate the theoretical results.