一类抛物型方程组的特征修正区域分解有限元方法

在实际生产和科学研究中,有许多物理问题的数学模型为抛物型方程组问题,如可压缩核废料污染问题,地下水资源问题,杨青提出了差分格式和有限元格式,应用先验估计得到了最优的l^2和L^2模误差估计,江城顺等利用交替方向有限元方法得到了H^1模和L^2模误差估计．杨国强等采用显式可解的三层差分格式求解二维方程组得到了H^1模误差估计．     

混合边界条件下的三维定常 旋转Navier-Stokes方程的原始变量有限元计算

本文利用原始变量有限元法求解混合边界条件下的三维定常旋转Navier-Stokes方程,证明了离散问题解的存在唯一性,得到了有限元解的最优误差估计.给出了求解原始变量有限元逼近解的简单迭代算法,并证明了算法的收敛性.针对三维情况下计算资源的限制,采用压缩的行存储格式存储刚度矩阵的非零元素,并利用不完全的LU分解作预处理的GMRES方法求解线性方程组.最后分析了简单迭代和牛顿迭代的优劣对比,数值算例表明在同样精度下简单迭代更节约计算时间.     

可压缩流动的Fourier谱-有限元解法

本文考虑n维(n=2,3)可压缩流动的带有单向周期边值条件问题的数值解.我们在周期方向采用Fourier谱方法,在非周期方向采用有限元方法,从而构造了一类谱-有限元格式.文中严格分析了计算误差,得到了收敛阶的估计.     

抛物型方程基于POD方法的有限元格式

将特征正交分解(Proper Orthogonal Decomposition, 简记为POD)方法应用于抛物型方程通常的有限元格式, 简化其为一个计算量很少但 具有足够高精度的POD有限元格式, 并给出简化的POD有限元解的误差分析. 数值例子表明在简化的POD有限元解和通常的有限元解之间的误差足够小的情形下, POD有限元格式比通常的有限元格式大大地节省计算量, 从而验证POD方法的有效性.     

外问题涡团法的收敛性

本文考虑欧拉方程初边值外问题涡团法的格式，其中泊松方程用等参有限元求解，这种方法被认为可以节省计算量．同时对半离散和全离散的格式都得到了丰满的误差估计．     

非定常不可压Navier-Stokes方程两重网格方法的稳定性和L2误差估计

讨论了二维非定常不可压Navier-Stokes方程的两重网格方法．此方法包括在粗网格上求解一个非线性问题,在细网格上求解一个Stokes问题．采用一种新的全离散(时间离散用Crank-Nicolson格式,空间离散用混合有限元方法)格式数值求解N-S方程．证明了该全离散格式的稳定性．给出了L2误差估计．对比标准有限元方法,在保持同样精度的前提下,TGM能节省大量的计算量．     

含弥散平面两相渗流驱动问题有限元方法的误差分析

研究了平面两相渗流可压缩问题含弥散情形的矩形有限元格式．引进一类插值算子,通过插值函数证明了有限元解的最优误差估计．     

二阶Dirichlet问题的三维等参元逼近及其数值积分格式

本文研究了求解三维二阶Drichlet问题的四面体等参元逼近格式.基于对等参变换的细致估计,证明了格式的收敛性,并得到了H~1-误差的最优估计.为简化等参元的复杂计算,本文还构造了一种简单的数值积分格式,并证明了其收敛性.     

无界区域瞬时涡流问题有限元-边界元耦合的A-φ法的误差分析

针对三维无界区域带有凸多边形导体的瞬时涡流问题,本文提出了一种基于势场的有限元-边界元耦合的方法,从理论上讨论了其能量模误差估计.虽然电场被分解为电矢势A与磁标势φ的梯度之和后增加了方程与未知量的个数,但这种分解可以很好地处理不同介质间的间断.与传统的A-φ法不同,本文讨论了一种全离散的A-φ解耦形式,这样不仅可以避免传统格式所产生的鞍点问题的求解,又可以减少计算量.     

非定常不可压缩Navier-Stokes方程的旋度形式压力校正格式分析

将求解不可压缩流动的旋度形式压力校正格式从Stokes方程延拓到非定常不可压缩Navier-Stokes方程.在第一步需要求解一个线性的对流-扩散方程,在第二步求解一个Stokes问题.首先给出格式的稳定性分析,然后采用P_2—P_1元分别使用标准形式的压力校正格式和旋度形式的压力校正格式进行了数值模拟,模拟结果表明,对于速度的L~2,H~1误差,两种格式几乎一样,对于压力的L~2误差,旋度形式的压力校正格式略有改进.     

Evaluations of Parameters for the Optimization of S.S.O.R. and A.D.I. Preconditioning

The main goal of this paper is to give two ways to estimate the needed parameters in order to obtain the condition number of S.S.O.R. preconditioned matrices, namely, the algebraic matricial formulation of convexity Riesz theorem and the tridiagonal Fourier analysis. The improvement with respect to Axelsson's approach is explicitly given. Estimations of the condition number in the case of A.D.I. preconditioning is also considered.     

Optimization of the Relaxation Parameter for S.S.O.R. and A.D.I. Preconditioning

We present and compare several approaches for the optimization of the relaxation parameter both for A.D.I. and S.S.O.R. basic iteration and preconditioning conjugate gradient method. For each kind of preconditioning a detailed link between estimates of the spectral radius of the iteration matrix and of the condition number resulting from preconditioning is proposed. It allows to choose the best approach in order to obtain the optimal relaxation parameter and the corresponding optimal estimates either of the spectral radius of the iteration matrix and of the resulting condition mumber of the S.S.O.R. and A.D.I. preconditioning.     

O.R.C.S. 1981 Results

Journal of the Operational Research Society -     

p.p. rings and generalized p.p. rings

This paper concerns two conditions, called right p.p. and generalized right p.p., which are generalizations of Baer rings and von Neumann regular rings. We study the subrings and extensions of them, adding proper examples and counterexamples to some situations and questions that occur naturally in the process of this paper.     

Two-point b.v.p. for multivalued equations with weakly regular r.h.s.

A two-point boundary value problem associated to a semilinear multivalued evolution equation is investigated, in reflexive and separable Banach spaces. To this aim, an original method is proposed based on the use of weak topologies and on a suitable continuation principle in Fréchet spaces. Lyapunov-like functions are introduced, for proving the required transversality condition. The linear part can also depend on the state variable x and the discussion comprises the cases of a nonlinearity with sublinear growth in x or of a noncompact valued one. Some applications are given, to the study of periodic and Floquet boundary value problems of partial integro-differential equations and inclusions appearing in dispersal population models. Comparisons are included, with recent related achievements.