非结构网格上一类满足局部极值原理的三阶精度有限体积方法

对二维标量双曲型守恒律方程,发展了一类满足局部极值原理的非结构网格有限体积格式.其构造思想是,以单调数值通量为基础,通过应用基于最小二乘法的二次重构和极值限制器,使数值解满足局部极值原理.为保证数值解在光滑区域达到三阶精度,该格式可结合局部光滑探测器使用.本文从理论上分析了格式的稳定性条件,数值实验验证了格式的精度和对间断的分辨能力.     

非协调有限元的全塑性分析

本文采用满足相容条件的非协调有限元模型以解决全塑性分析中有限元解的数值精度问题.文中讨论了该模型适用于全塑性分析的机理和判据,还设计了一个确定塑性极限载荷的算法.     

求解带有移动界面的线性对流方程的浸入界面有限体积法（英文）

本文为一类带有移动界面的守恒律方程提出了耦合高分辨率格式的数值算法.这种算法是在一致大小的笛卡尔网格上导出而满足标准的双曲型稳定条件.文末列举数值算例研究这种算法的收敛性和数值精度.     

二阶扩张状态观测器的误差分析

本文利用分片光滑的Ｌｙａｐｕｎｏｖ函数来进行扩张状态观测器的误差分析,指出了非线性扩张状态观测器有更好的估计精度能力,给出了为提高估计精度扩张状态观测器的参数所满足的条件．     

双曲型守恒律方程的熵稳定格式的一些讨论

本文讨论双曲型守恒律方程的熵稳定格式.对于给定的熵对,格式所满足的熵条件中的数值熵通量是不唯一的.Tadmor的充分条件可以唯一地确定标量方程的熵守恒通量,但不能唯一确定方程组的熵守恒通量,却可以给出方程组的空间一阶精度的熵守恒格式.也讨论了在熵守恒通量上添加数值粘性得到的显式熵稳定格式需要满足的条件及常见的时间离散对熵守恒和熵稳定的影响.     

板弯曲断裂计算的边界配置解法

采用复变函数理论和边界配置方法,分析计算了Kirchhoff板的弯曲断裂问题.假设了位移及内力的复变函数式,它们能满足一系列的基本方程和支配条件,例如域内的平衡方程、裂纹表面的边界条件、裂纹尖端的应力奇异性质.这样,仅板边界的边界条件需要考虑.它们可用边界配置法和最小二乘法近似满足.对不同边界条件和载荷情形进行了分析计算.数值算例表明,本文方法精度较高,计算量小,是一种有效的半解析、半数值计算方法.     

曲梁单元和它的收敛率

本文给出了拟协调曲梁和扁曲梁单元.数值结果表明,用于近似曲梁的拟协调曲梁和扁曲梁单元较直梁单元具有更好的精度.由位移法构造的曲梁单元不能够满足刚体位移条件,为了近似地满足刚体位移条件必须用很多的单元.本文证明了直梁单元、拟协调曲梁和扁曲梁单元,当单元尺寸无限缩小时,具有相同的收敛速度O(l~2),当使用均匀网格时,其中l是单元的长度.     

非线性不适定问题一种双循环的牛顿型迭代格式

本文研究非线性算子方程F(x)=y的解,结合最速下降法,Newton-Landweber迭代格式及正则化思想,在F满足适当的条件下,构造出新的双循环迭代格式.本文对格式的收敛性进行了严格论证,并估计出迭代格式的收敛精度.     

Cahn-Hilliard方程非条件能量稳定二阶方法的稳定性和收敛性分析

在这项工作中,我们构建了一种非条件能量稳定的高效不变能量四分法(IEQ)来求解Cahn-Hilliard方程.所构建的数值格式是线性的、具有二阶时间精度和非条件能量稳定性.我们仔细分析了数值格式的唯一可解性、稳定性和误差估计.结果表明,所构建的格式满足唯一可解性、非条件能量稳定性和时间方向的二阶收敛性.通过大量的二维和三维数值实验,我们进一步验证了所提出格式的收敛阶、非条件能量稳定性和有效性.     

非协调曲边四边形十二自由度平板弯曲单元*

本文利用精确元法[1],给出一个十二自由度曲边四边形板弯曲单元.该方法不需要变分原理,适用于任意正定和非正定偏微分方程.利用这个方法,单元之间的协调条件很容易满足,仅须位移和内力在单元节点上连续,即可保证所得到的解收敛于精确解.利用本文方法所获得的解,无论是位移还是内力可同时有二阶收敛精度.文末给出数值算例.表明了本文所得到的单元有非常好的精度.     

Numerical analysis for a variable-order nonlinear cable equation

In this paper, a variable-order nonlinear cable equation is considered. A numerical method with first-order temporal accuracy and fourth-order spatial accuracy is proposed. The convergence and stability of the numerical method are analyzed by Fourier analysis. We also propose an improved numerical method with second-order temporal accuracy and fourth-order spatial accuracy. Finally, the results of a numerical example support the theoretical analysis.     

Numerical analysis for a variable-order nonlinear cable equation

In this paper, a variable-order nonlinear cable equation is considered. A numerical method with first-order temporal accuracy and fourth-order spatial accuracy is proposed. The convergence and stability of the numerical method are analyzed by Fourier analysis. We also propose an improved numerical method with second-order temporal accuracy and fourth-order spatial accuracy. Finally, the results of a numerical example support the theoretical analysis.     

一种改进的光电经纬仪两站交会测量方法

提出一种新的光电经纬仪2站交会定位方法。该方法是在常用定位公式的基础上,实时计算弹道参数的精度,根据解算的弹道参数精度变化,动态确定加权系数。仿真结果表明：在不等精度光电经纬仪交会定位情况下,该方法能够有效提高两站定位的精度。     

非饱和土壤水流问题的CN有限元格式

首先给出二维非饱和土壤水流问题基于Crank-Nicolson(CN)方法的具有时间二阶精度的半离散化格式,然后直接从CN时间半离散化格式出发,建立具有时间二阶精度的全离散化CN有限元格式,并给出误差估计,最后用数值例子说明全离散化CN有限元格式的优越性.这种方法可以绕开关于空间变量的半离散化格式的讨论,提高时间离散的精度,极大地减少时间方向的迭代步,从而减少实际计算中截断误差的积累,提高计算精度和计算效率.     

Convergence analysis of generalized nonlinear inexact Uzawa algorithm for stabilized saddle point problems

This paper deals with a modified nonlinear inexact Uzawa (MNIU) method for solving the stabilized saddle point problem. The modified Uzawa method is an inexact inner-outer iteration with a variable relaxation parameter and has been discussed in the literature for uniform inner accuracy. This paper focuses on the general case when the accuracy of inner iteration can be variable and the convergence of MNIU with variable inner accuracy, based on a simple energy norm. Sufficient conditions for the convergence of MNIU are proposed. The convergence analysis not only greatly improves the existing convergence results for uniform inner accuracy in the literature, but also extends the convergence to the variable inner accuracy that has not been touched in literature. Numerical experiments are given to show the efficiency of the MNIU algorithm.     

抛物型方程的一种高精度区域分解有限差分算法

1引言 近年来,区域分解算法以可以将大型问题分解为一系列小型问题以减少计算规模及算法可高度并行实现等特点受到了人们的广泛关注.前人也做了很多很好的工作:参考文献[1]中C.N.Dawson等人提出了显一隐格式的区域分解算法,在时间层不分层的内边界点采用大步长向前-中心差分显格式及在内点采用古典隐格式,取得的精度为O(△t+h2+H3).参考文献[2]中给出了[1]中区域分解算法对于内边界点为等距分布的多子区域时的新的误差估计,使含H3误差项的系数比[1]中缩小了一倍.还将采用大步长日的saul'yev的非对称差分格式应用于内边界点,并给出了两个子区域和多个子区域情形下差分解的先验误差估计.     

非饱和土壤水流方程的CN广义差分法

首先给出二维非饱和土壤水流方程时间二阶精度的Crank-Nicolson(CN)时间半离散化格式,然后直接从CN时间半离散化格式出发,建立具有时间二阶精度的全离散化CN广义差分格式,并给出误差分析,最后用数值例子验证全离散化CN广义差分格式的优越性.这种方法能提高时间离散的精度,极大地减少时间方向的迭代步,从而减少实际计算中截断误差的积累,提高计算精度和计算效率.而且该方法可以绕开对空间变量的半离散化广义差分格式的讨论,使得理论研究更简便.     

Numerical simulation for the variable-order Galilei invariant advection diffusion equation with a nonlinear source term

In this paper, we consider the variable-order Galilei advection diffusion equation with a nonlinear source term. A numerical scheme with first order temporal accuracy and second order spatial accuracy is developed to simulate the equation. The stability and convergence of the numerical scheme are analyzed. Besides, another numerical scheme for improving temporal accuracy is also developed. Finally, some numerical examples are given and the results demonstrate the effectiveness of theoretical analysis.     

Sixth-Order Compact Extended Trapezoidal Rules for 2D Schrödinger Equation

Based on high-order linear multistep methods (LMMs), we use the class of extended trapezoidal rules (ETRs) to solve boundary value problems of ordinary differential equations (ODEs), whose numerical solutions can be approximated by boundary value methods (BVMs). Then we combine this technique with fourth-order Padé compact approximation to discrete 2D Schrödinger equation. We propose a scheme with sixth-order accuracy in time and fourth-order accuracy in space. It is unconditionally stable due to the favourable property of BVMs and ETRs. Furthermore, with Richardson extrapolation, we can increase the scheme to order 6 accuracy both in time and space. Numerical results are presented to illustrate the accuracy of our scheme.     

基于精度的属性约简新方法

自Pawlak提出粗糙集概念以来,人们就一直对粗糙集的近似精度很感兴趣,出现了不少有关近似精度的文献.在粗糙集理论中,精度是量化由粗糙集边界引起的不精确性的一种重要数字特征.在分析传统精度和基于等价关系图的过剩熵的近似精度的基础上,提出了一种新的精度定义.比较发现,新定义的精度更具有合理性.同时把这个新定义的精度运用到了属性约简上,通过实例比较发现,本文提出的属性约简更具有可行性.