求解定常不可压Stokes方程的两层罚函数方法

借助于两套有限元网格空间提出了一种求解定常不可压Stokes方程的两层罚函数方法.该方法只需要求解粗网格空间上的Stokes方程和细网格空间上的两个易于求解的罚参数方程（离散后的线性方程组具有相同的对称正定系数矩阵）.收敛性分析表明粗网格空间相对于细网格空间可以选择很小，并且罚参数的选取只与粗网格步长和问题的正则性有关.因此罚参数不必选择很小仍能够得到最优解.最后通过数值算例验证了上述理论结果，并且数值对比可知两层罚函数方法对于求解定常不可压Stokes方程具有很好的效果.     

并行间断有限元算法求解Navier-Stokes方程

间断Galerkin有限元方法非常适合在非结构网格上高精度求解Navier-Stokes方程,然而其十分耗费计算资源.为了提高计算效率,提出了高效的MIMD并行算法.采用隐式时间离散GMRES+LU SGS格式,结合多重网格方法,当地时间步长加速算法收敛.为了保证各处理器间负载平衡,采用区域分解二级图方法划分网格,实现内存合理分配,数据只在相邻处理器间传递.数值模拟了RAE2822翼型和M6黏性绕流,加速比基本呈线性变化且接近理想值.结果表明了该算法能有效减少计算时间、合理分配内存,具有较高的加速比和并行效率,适合于MIMD粗粒度科学计算.     

Stokes方程的投影／方向分裂谱方法

我们提出和分析了一种求解Stokes方程的数值方法．新方法基于空间上的Legendre谱离散,时间上则采用投影／方向分裂格式．更确切地说,时间离散的出发点是旋度形式的压力校正投影法,在此基础上进一步应用方向分裂法,把速度和压力方程分裂为一系列一维的椭圆型子问题．然后生成的这些一维子问题用Legendre谱方法进行空间离散．另外,我们证明了全离散格式的稳定性．一些数值实验验证了收敛性和方法的有效性．     

关于非定常不可压Navier-Stokes方程的时间高精度隐式差分方法

The incompressible Navier-Stokes equations,upon spatial discretization,become a system of differential algebraic equations,formally of index2.But due to the special forms of the discrete gradient and disrete divergence,its index can be regarded as 1.Thus,in this paper,a systematic approach following the ODE theory and methods is presented for the construction of high-order time-accurate implicit schemes for the incompressible Navier-Stokes equations,with projection methods for efficiency of numerical solution.The 3rd order 3-step BDF with componentconsistent pressure-correction projection method is a first attempt in this direction;the related iterative solution of the auxiliary velocyty,the boundary conditions and the stability of the algorithm are discussed.Results of numerical tests on the incompressible Navier-Stokes equations with an exact solution are presented,confirming the accureacy,stability and component-consistency of the proposed method.     

定常Navier—Stokes方程的Petrov—Galerkin方法

研究了定常Navier－Stokes方程的四种Petrov－Galerkin有限元方法：PG1，PG2，SD和GLS．它们都是稳定的，避免了经典混合方法中必要的Babuska－Brezzi条件．给出了各种方法有限元解的存在性、唯一性和唯一解的误差估计．     

求解非线性算子方程的加速迭代格式

研究非线性算子方程的近似求解方法.首先对通常的求解非线性方程加速迭代格式进行推广,得到高阶收敛速度的加速迭代格式,最后把这种加速迭代格式推广到非线性算子方程的求解中去,利用非线性算子的渐进展开,证明了这种加速格式具有三阶的收敛速度.     

非定常Stokes方程的全离散稳定化有限元格式

本文研究二维非定常Stokes方程全离散稳定化有限元方法.首先给出关于时间向后一步Euler半离散格式,然后直接从该时间半离散格式出发,构造基于两局部高斯积分的稳定化全离散有限元格式,其中空间用P_1—P_1元逼近,证明有限元解的误差估计.本文的研究方法使得理论证明变得更加简便,也是处理非定常Stokes方程的一种新的途径.     

求Lyapunov矩阵方程的双对称解的迭代算法

本文研究了Lyapunov矩阵方程.利用共轭梯度法,建立了求该矩阵方程双对称解的迭代算法.同时,也能给出指定矩阵的最佳逼近双对称矩阵.     

连续Sylvester矩阵方程求解的分裂迭代算法

有效求解连续的Sylvester矩阵方程对于科学和工程计算有着重要的应用价值,因此该文提出了一种可行的分裂迭代算法.该算法的核心思想是外迭代将连续Sylvester矩阵方程的系数矩阵分裂为对称矩阵和反对称矩阵,内迭代求解复对称矩阵方程.相较于传统的分裂算法,该文所提出的分裂迭代算法有效地避免了最优迭代参数的选取,并利用了复对称方程组高效求解的特点,进而提高了算法的易实现性、易操作性.此外,从理论层面进一步证明了该分裂迭代算法的收敛性.最后,通过数值算例表明分裂迭代算法具有良好的收敛性和鲁棒性,同时也证实了分裂迭代算法的收敛性很大程度依赖于内迭代格式的选取.     

含高次逆幂的矩阵方程对称解的双迭代算法

本文研究了在控制理论和随机滤波等领域中遇到的一类含高次逆幂的矩阵方程的等价矩阵方程对称解的数值计算问题.采用牛顿算法求等价矩阵方程的对称解,并采用修正共轭梯度法求由牛顿算法每一步迭代计算导出的线性矩阵方程的对称解或者对称最小二乘解,建立了求这类矩阵方程对称解的双迭代算法,数值算例验证了双迭代算法是有效的.     

The Calahan Method for Parabolic Equations with Time-Dependent Coefficients

A modified Calahan method for parabolic equations with time-dependent coefficients is presented. It is shown that the convergence order is $O(h^{r+1}+k^3)$ while the convergence order obtained in [1] for a standard Calahan method is only $O(h^{r+1}+k^2)$.     

Band-Toeplitz Preconditioned GMRES Iterations for Time-Dependent PDEs

Nonsymmetric linear systems of algebraic equations which are small rank perturbations of block band-Toeplitz matrices from discretization of time-dependent PDEs are considered. With a combination of analytical and experimental results, we examine the convergence characteristics of the GMRES method with circulant-like block preconditioning for solving these systems.This revised version was published online in October 2005 with corrections to the Cover Date.     

预处理Uzawa算法及其在求解Stokes问题上的应用

在求解鞍点问题的经典Uzawa算法收敛性的基础上,对预处理Uzawa算法收敛性做出进行进一步的研究,得到其收敛的充要条件及误差传播矩阵的谱半径;并将其应用到Mini元离散求解Stokes问题中,通过数值计算验证所得结论的正确性.     

A Semi-Lagrangian Runge-Kutta Method for Time-Dependent Partial Differential Equations

In this paper, a Semi-Lagrangian Runge-Kutta method is proposed to com-pute the numerical solution of time-dependent partial di®erential equations.The method is based on Lagrangian trajectory or the integration from the de-parture points to the arrival points (regular nodes). The departure points aretraced back from the arrival points along the trajectory of the path. The highorder interpolation is needed to compute the approximations of the solutionson the departure points, which most likely are not the regular nodes. On thetrajectory of the path, the similar techniques of Runge-Kutta are applied to theequations to generate the high order Semi-Lagrangian Runge-Kutta method.The numerical examples show that this method works very effient for thetime-dependent partial di®erential equations.     

用摄动配置方法求解含时薛定谔方程

将摄动配置方法应用到含时薛定谔方程,在计算实现的基础上结合摄动配置的特征提出了一类新的数值积分方法,并给出了一个2级2阶和一个3级4阶的辛摄动配置方法对含时薛定谔方程的数值算例.为了检验新的数值积分方法,我们还给出了与两个辛摄动配置格式在理论上等价的辛龙格-库塔方法以及同阶的非辛方法的数值模拟.展示了一些数值结果,并给出了一些分析.     

An Analysis of Penalty-Nonconforming Finite Element Method for Stokes Equations

In this paper, the penalty-nonconforming finite element method for Stokes equations is considered. An abstract error estimate is given. For Crouzeix-Raviart nonconforming triangular elements, in particular, the analysis shows that the "reduced integration" technique is not necessary in the integration of the penalty term on each element. It means that a loss of precision is avoided in this penalty method.     

Nonconforming Mixed Finite Element Method for Time-Dependent Maxwell's Equations with ABC

In this paper, a nonconforming mixed finite element method (FEM) is presented to approximate time-dependent Maxwell's equations in a three-dimensional bounded domain with absorbing boundary conditions (ABC). By employing traditional variational formula, instead of adding penalty terms, we show that the discrete scheme is robust. Meanwhile, with the help of the element's typical properties and derivative transfer skills, the convergence analysis and error estimates for semi-discrete and backward Euler fully-discrete schemes are given, respectively. Numerical tests show the validity of the proposed method.     

The Dual Singular Function Method for 2D Boundary Integral Equations

The accuracy of standard boundary element methods for elliptic boundary value problems deteriorates if the boundary of the domain contains corners or if the boundary conditions change along the boundary. Here we first investigate the convergence behaviour of standard spline Galerkin approximation on quasi-uniform meshes for boundary integral equations on polygonal domains. It turns out, that the order of convergence depends on some constant describing the singular behaviour of solutions near corner points of the boundary. In order to recover the full order of convergence for the Galerkin approximation we propose the dual singular function method which is often used for improving the accuracy of finite element methods. The theoretical convergence results are confirmed and illustrated by a numerical example.     

非定常线性化Navier-Stokes方程的非协调流线扩散有限元法分析

对非定常线性化Navier-Stokes方程提出了非协调流线扩散有限元方法．用向后Euler格式离散时间,用流线扩散法处理扩散项带来的非稳定性．速度采用不连续的分片线性逼近,压力采用分片常数逼近．得到了离散解的存在唯一性以及在一定范数意义下离散解的稳定性和误差估计．