解双调和方程边值问题离散化方程组的PEk方法

建立了求解双调和方程边值问题离散化得到的大型块五对角线性代数方程组的PEk方法，对系数矩阵为Hermite正定矩阵的情形，证明了PEk方法的可解性和收敛性，并给出了参数k的选取范围。     

Cauchy奇异积分及积分方程数值方法的讨论

本文给出新的积分公式,利用此积分公式来求解带有Cauchy核奇异积分。在此基础上讨论带有Cauchy核的奇异积分方程的数值解,并且与真解进行比较,验证此积分公式的优越性。     

椭圆类平片裂纹二维有限部积分方程的多项式数值方法

根据均质弹性体中平面裂纹问题的一维Ｃａｕｃｈｙ型主值积分方程的Ｃｈｅｂｙｓｈｅｖ多项式数值求解方法,提出了三维断裂力学问题的椭圆类平片裂纹二维有限部积分方程中未知位移间断用Ｃｈｅｂｙｓｈｅｖ多项式与位移间断基本函数之积来表示的近似数值解法,并导出了与多项式系数相对应的应力强度因子计算公式最后给出了若干不同长短轴半径之比的椭圆平片裂纹应力强度因子计算例计算表明,本文方法的数值结果不但收敛速度快,而且精度也大大高于现有的有限部积分———边界元方法的精度     

变限积分的有限体积法解决对流扩散方程

针对一维对流扩散方程,基于变限积分的有限体积法,提出一种高精度全离散方法。该方法在控制体内对方程进行变限积分,引入变限因子,然后分别对上下限再进行积分,从而将微分方程转化为积分方程,最后运用插值的方法对目标函数进行近似代替。该方法提高了计算精度,结果得到一维的收敛精度。采用Fourier分析法证明了格式为条件稳定。最后给出了非稳态和稳态2种情况下的数值算例,验证了所提出的格式具有高精度和易于编程计算的优点。     

化抛物型方程为边界积分方程的一个方法

给出一个化二阶线性抛物方程为边界积分方程的方法①．     

解一类积分方程的两种方法

本文主要讨论用两种方法解一类常见的积分方程,第一种方法,先对积分方程两端微分,将其转化成与之等价的微分方程,然后利用径向基函数配置法对所得到的微分方程求数值解,也即积分方程的数值解;第二种方法,直接利用径向基函数配置法对积分方程求数值解。最后,分别用两种方法对一个积分方程进行了数值试验,验证了两种方法的可行性。     

双材料裂纹问题的积分方程方法

以双材料单位力基本解为基础,推导得到了单位位移不连续基本解.以裂纹上下表面相对位移为未知量,建立相应的奇异积分方程,为了数值求解积分方程,构建了一种奇异单元,以反映裂纹上下表面相对位移在裂纹周边附近的r~(1/2)奇性(非界面裂纹)和r~(1/2)·r~(iε)Ⅰ-Ⅱ型耦合振荡奇性(界面裂纹).对数值计算方法、特别是奇异和超奇异积分的计算进行了研究,最后通过若干算例考证了本方法的精度.所建立的积分方程及数值方法,是一种求解双材料平面内多裂纹系统(界面和非界面)的一般方法.     

一类积分方程

1 我们研究积分方程 (1.1) integral from -1 to 1(q(ξ)k((ξ-x)/λ)dξ)=πf(x) (|x|≤1,λ∈(0,∞)) 这里q(x)是要求的函数,k(t)和f(x)是已知函数,这个方程的核是 (1.2) k(t)=integral from 0 to ∞([∧_1(u)cosut-∧_2(u)sinut]du/u)     

重调和方程的混合有限元的区域分解算法

主要考虑将区域分解算法应用于混合有限元方法的情形。基于Schwarz交替法,讨论了重调和方程的混合有限元格式的区域分解算法,证明了它的收敛性,利用区域分解技术,给出它的有限离散格式和预处理矩阵。本文表明基于Lions框架的Schwarz算法也适用于混合有限元。     

抛物型偏微分方程的有限解析法

本文提出了抛物型偏微分方程的双层四点格式的有限解析数值格式,证明了该格式的稳定、收敛性、相容性,给出了求解程序。并与文献[1]中双层六点格式,结合算例进行了比较,结果表明,该方法比[1]中的六点格式所用的CPU至少节约一半,最后对双层四点格式进行了讨论。     

用有限元方法解带粘滞项的弹性方程

本文对实际地震勘探仿真中碰到的弹性粘滞波方程用有限元方法对之求数值解作了尝试,给出了有限元计算公式。     

区间有限元静力控制方程的一种迭代解法

当结构的不确定参量可用区间限界时，将区间分析和有限元方法相结合，可建立起区间有限元方法，区间控制方程组的求解是其核心问题，根据所推导出的区间变量的运算特性，提出了求解区间有限元静力控制方程的一种迭代解法。将区间方程组的求解归结为一点值迭代方过程，算例分析表明了中的方法计算量小，易于实施，是实用和可行的。     

二维及三维 N-S 方程的混合有限分析多重网格方法

提出了二维及三维Ｎ－Ｓ方程的快速数值方法—混合有限分析多重网格方法（ＨＦＡ－ＭＧ）,标准算例的计算结果证明ＨＦＡ－ＭＧ在精确地模拟复杂流动特征的前提下,能显著地降低计算工作量。     

二阶椭圆方程自适应最小二乘混合有限元法

研究了二阶椭圆方程的自适应最小二乘混合有限元法,利用二次非协调有限元空间和Raviatr-Thomas有限元空间进行逼近,利用最小二乘函数构造了进行自适应计算的后验误差估计子,并进行了后验误差估计。     

高阶间断有限元方法求解Euler方程时的数值积分问题研究

文章对非结构网格上高阶间断有限元方法求解Euler方程时的数值积分方式进行了研究。首先根据间断有限元方法的基本原理推导出了有限元离散控制方程中线积分和体积分项的数值积分精度要求。然后给出了采用Gauss-Legendre和Gauss-Lobatto积分公式处理线积分项,以及采用Guass积分公式和重构积分方法处理体积分项的情况下,为满足积分精度要求所需使用的最少积分节点数目。最后,通过具体算例对上述积分精度要求进行了验证,并考察了不同数值积分方法对于求解效率和精度的影响。     

二阶椭圆型方程的有限体积解法

给出了有限体积法处理问题的一般原理,在现有文献的基础上对有限体积算法进行了改进,详细推导了一类二阶椭圆型方程的有限体积法列式,并结合具体的算例比较了有限体积法同有限元的计算精度,最后总结了有限体积法的若干优越性。     

Burgers方程的隐-显多步有限元方法

考虑Burgers方程的初边值问题,提出了隐-显多步有限元方法的逼近格式,并证明了该格式的最优阶误差估计。     

求解具有间断系数椭圆型方程的修正有限体积法

对经典的有限体积方法进行了修正,并用于求解具有间断系数的椭圆型方程.通过增加有限体积方法中通量函数泰勒展开式的项数来提高通量函数的计算精度,对与边界相邻的网格点改用中心差分离散,以使得边界处的精度更高.修正后的方法是有限体积方法和有限差分方法的结合形式.得到了格式二阶逐点收敛并且在界面处具有二阶精度.数值实验表明:新方法较经典有限体积方法中求解间断系数问题的算术平均法和调和平均法更具优势.