A2B分子反映系统经典轨迹计算的有限元方法

本文利用常微分方程的连续有限元法计算了A2B模型的分子的经典轨迹和能量误差，将计算延长到10^-sS，并与辛算法进行了比较．结果表明，在微观反映动力学研究所考虑的时间范围内，有限元法的结果与理论分析一致，能较长时间保持能量守恒和系统整体结构，并在数值计算上探讨了守恒性和近似程度，结果与理论相吻合。     

用龙格-库塔法求解非线性方程组

本文介绍了一种求解非线性方程组的新方法龙格-库塔法。     

有限元方法的超收敛性

有限元方法的理论研究主要围绕着解的收敛性及误差估计.早在六十年代,冯康教授在一般情况下已独立研究过这个问题.到现在,对大多数重要情形已证明m次分片多项式有限元解u_h有如下标准误差估计 ||u-u_h||_s=O(h~(m+1-s)),s=0,1.另一方面人们早在实际计算中注意到:在一定条件下有限元解或其导数在某些特殊点上可能有更高的精度.1972年前后,瑞典V·Thomee,美国C.de Boor及J.Douglas等人从不同角度发现并论证了有限元的这种奇特性质.Douglas称它为超收敛性(Supercon-     

三角形线元与二次元L~2投影的整体超收敛

Based on an orthogonal expansion in a triangle, superconvergence for L^2-projection to linear and quadratic triangular elements is studied.Assume that Ω is a polygonal domain,triangulation is uniform and Th is a set of all vertexesand side midpoints.Then, on Th,the average gradient ^-D(u-uh)=O(h^2) for linear element uh,and u-uh=O(h^4) for quadratic element uh. Under someboundary conditions,these properties upto the boundary are valid.     

非线性问题有限元的超收敛性

对线性椭圆边值问题有限元的超收敛性已有许多研究。它们大致可分为两类:一是直接研究有限元解在某些特殊点上“自然地”具有的超收敛性,另一类是利用有限元解作局部积分平均得到高精度数值。但是对非线性问题研究较少。本文在某些情形下证明了如下重要事实:上述超收敛性对非线性椭圆问题的有限元仍然成立(注。进一步结果可多看作者论文,Superconvergence of finite element approximations tonoulinear elliptic problems,1982年4月19—23日北京,中法有限元方法讨论会)。     

求非线性问题多解的搜索延拓法

基于-Δ的特征系λ_j,φ_j,逼近非线性问题 Δu+f(u)=0(在Ω中), u=0(在¶Ω上)的多重解. 提出了一种新的搜索延拓法(SEM),它由三层子空间上的三种算法组成. 对f(u)=u³,在正方形和L形区域上完成了数值实验.这些结果表明,对应于-Δ的每个k重特征值, 至少存在3^k-1个不同的非零解(猜想1).     

抛物方程时间连续有限元全离散格式的超收敛性

利用张量积分解和时间方向单元正交分解,证明了线性抛物型方程的时间连续全离散有限元在单元节点和内部的特征点的超收敛性.并用连续有限元计算了非线性Schrodinger方程,验证了能量的守恒性.计算结果与理论相吻合.     

一类偏积分微分方程二阶差分全离散格式

本文给出了数值求解一类偏积分微分方程的二阶差分全离散格式．时间方向采用了二阶向后差分格式,积分项的离散利用了Lubich的二阶卷积求积公式,给出了稳定性的证明、误差估计及收敛性的结果,并给出了数值例子．     

三角形二次插值系数有限元法解半线性椭圆问题的超收敛性

基于均匀三角形的剖分求解一类二阶半线性椭圆问题,用插值系数有限元方法比经典有限元法更容易实现,与经典二次有限元一样,二次插值系数有限元方法在对称点处也有四阶超收敛精度,数值计算表明这些结论是正确的.     

二阶常微初值问题的单步格式

讨论二阶常微分方程初值问题utt+au=f,u(0)=u0，ut(0)=v0 的一种单步格式，采用u及v=ut为未知量，计算简单. 证明了此格式的稳定性及对u，v皆有二阶精度. 此格式可用于双曲问题.