求解一维扩散反应方程的隐式高精度紧致差分格式

提出了一维扩散反应方程的一种隐式高精度紧致差分格式,空间二阶导数采用四阶紧致差分格式进行离散,时间导数采用四阶向后欧拉公式进行离散,格式截断误差为O（τ⁴+h⁴）,即时间和空间都可以达到四阶精度,最后通过数值实验验证了本文方法的精确性和可靠性.     

求解一维对流扩散反应方程的一种隐式差分格式

提出了数值求解一维非稳态对流扩散反应方程的一种隐式差分格式。首先将模型方程利用指数函数转化为对流扩散方程,构造它的差分格式,然后对差分方程的系数进行相应处理,并进行回代,得到对流扩散反应方程的隐式差分格式,其截断误差为O(τ2+h2),采用von Neumann方法证明了格式是无条件稳定的,并且由于每一时间层上只用到了3个网格点,所以可直接采用追赶法求解差分方程,数值结果显示了算法的有效性。     

二维非定常对流扩散方程的隐式多重网格方法

提出了数值求解二维非定常对流扩散方程的一种新的加权平均隐格式，利用Fourier分析方法证明了该格式是无条件稳定的．利用多重网格加速技术。提出了基于时间修正的多重网格的全近似格式(FAS)，从而克服了传统迭代法在求解隐格式时收敛速度慢的缺陷，大大加快了迭代收敛速度。提商了问题的求解效率．数值计算结果表明，修正的多重网格FAS格式较传统的多重网格粗网格校正格式(CS)具有更好的收敛效率．并且随着σ和s的增大，它可以将传统迭代法的收敛速度提商几十倍，甚至几百倍．     

扩散反应方程的三次样条高阶差分格式

基于发展的三次样条是差分公式，提出了两种数值求解含源汇扩散反应方程的二层三结点高，精度差分格式，格式推导过程简便，精度在时间和空间方向上分别达到二阶和四阶，且均为无条件稳定的，最后通过数值算例检验了格式的优良性态。     

求解二维扩散方程的加权平均隐式多重网格方法

提出了数值求解二维扩散方程两种精度分别为O(τ^2 h^2)和O(τ^2 h^4)的无条件稳定的加权平均隐格式，并采用多重网格方法进行求解，从而克服了传统迭代法在求解隐格式时收敛速度慢的缺陷，提高了求解效率．数值实验验证了该方法的精确性和可靠性．     

求解扩散方程的一种高精度隐式差分方法

利用一阶微商和二阶微商的四阶紧致差分逼近公式,推导出了数值求解一维扩散方程的两种新的高精度隐式紧致差分格式,其截断误差分别为O(τ^2 h^4)和O(τ^4 h^4)．通过Fourier分析方法证明了格式O(τ^2 h^4)是无条件稳定的,而格式O(r^4 h^4)是无条件不稳定的．并且由于每一时间层上只用到了3个网格点,所以差分方程可采用追赶法直接进行求解．     

扩散方程的高精度加权差分格式

基于已发展的二阶微商三次样条四阶逼近公式,提出了数值求解含源汇项扩散方程的二层三点且在空间方向上达到四阶精度的加权差分格式．通过Ｆｏｕｒｉｅｒ方法讨论了文中格式的稳定性．证明了当１／２θ１时,格式是无条件稳定的,而当０θ＜１／２时,只有０＜ｒ１／［３（１－２θ）］,格式是稳定的．其中θ是权参量,ｒ＝Ｄτ／ｈ２为Ｆｏｕｒｉｅｒ数,Ｄ为扩散系数,而τ,ｈ分别为时间和空间网格长度     

扩散方程高精度加权差分格式的MATLAB实现

对扩散方程混合问题,利用二阶微商三次样条四阶逼近公式构造其四阶加权差分格式.使用MATLAB软件编程,将问题的解用图像表示出来,通过数值结果验证了该方法的可行性和稳定性.     

二维波动方程的高精度隐格式及其多重网格算法

提出了数值求解二维波动方程的两种高精度三层紧致隐格式．利用Fourier分析方法证明了格式均是无条件稳定的．并在此基础上提出了求解该问题的多重网格算法．从而克服了传统迭代法在求解隐格式时收敛速度慢的缺陷，大大加快了迭代收敛速度．提高了求解效率．数值实验结果验证了方法的精确性和可靠性．     

求解对流扩散反应方程的一种高精度紧致差分方法

首先,针对一维对流扩散反应方程,借助截断误差余项修正的方法,将中心差分格式余项中未知函数的三阶和四阶导数项利用一阶导数的表达式来代替,从而提出一种新的紧致差分格式,具有四阶精度.然后,为了简化计算,对格式常系数形式的耗散误差和色散误差进行分析,证实该格式的低耗散性.接着,将该方法推广到二维,运用降维的思想转化成2个一维形式的定常对流扩散反应方程,并用求解一维方程的方法,离散后相加即得二维对流扩散反应方程的紧致差分格式.最后,通过数值实验验证本文格式的精确性和可靠性.     

二维波动方程的加权平均隐格式及多重网格算法

提出了数值求解二维波动方程的一种加权平均隐式差分格式，理论分析结果表明其为无条件稳定的，为了克服传统迭代法在求解隐格式方面的困难，采用了多重网格算法，大大加快了迭代收敛速度，提高了求解效率，数值实验结果验证了方法的有效性和可靠性。     

求解热传导方程的一族三层隐式差分格式

提出求解热传导方程的一族高精度三层九点隐式格式,格式的截断误差为O(τ2+h4).利用Fourier方法证明差分格式是绝对稳定的.并通过数值试验,比较差分格式的解和精确解,说明差分格式的有效性.     

二维分数阶扩散方程交替差分格式及其一致性

研究二维有限域上的空间分数阶扩散方程的数值解法,通过移位的Grunwald公式对空间分数阶导数进行离散,得到Euler隐式差分格式。利用傅里叶变换理论证明了交替差分格式的一致性。     

具有非局部边值条件的波动方程的加权隐式差分格式

许多物理现象是由具有非局部条件的双曲型方程描述的.具有非局部条件的双曲型方程的数值解法是一个重要研究领域,在现代科学与技术科学有广泛应用.本文讨论了一类具有非局部边值条件的双曲型方程的数值解.通过引入新的未知函数将一类具有非局部边值条件的波动方程定解问题变为Dirichlet和Neumann边值问题,作者给出了该问题的加权隐式差分格式,证明了该差分格式的唯一可解性,利用Fourier方法给出了上述差分格式的稳定性条件.给出的数值例子用以说明差分格式稳定性和收敛性.     

泊松方程的高精度三次样条差分方法

给出一种求解泊松方程的新数值方法：以二维泊松方程为例，首先将其转化成一维方程，然后将根据由三次样条插值公式导出的四阶精度三次样条差分公式，应用到一维方程之中，最终建立起二维泊松方程矩形网格下九结点差分格式，并给出了误差估计和数值结果。     

求解三维非定常对流扩散方程的隐式差分方法

提出了一种数值求解三维非定常变系数对流扩散方程,对角占优、空间为二阶精度的隐格式,利用Fourier分析方法证明了该格式是无条件稳定的,并且由于格式具有对角占优性,因此适合于大梯度(高雷诺数)问题的数值求解.另外,为了克服传统迭代法在求解隐格式时收敛速度慢的缺陷,采用了多重网格加速技术,大大加快了迭代收敛速度,提高了求解效率.数值实验结果证明了该方法的精确性、稳定性和对高网格雷诺数问题的强适应性.     

有限元有限差分法二维波动逆时偏移初探

提出了采用有限元有限差分实现二维波动方程的逆时偏移算法。该方法在空间上 ,联合采用有限元法和有限差分法 ;对于地表 (水平 )方向 ,使用有限元法进行离散 ,将原方程转化为一个一维 (深度和时间 )问题的方程组 ;在深度和时间方向上 ,采用有限差分法来求解。介绍了算法的基本原理 ,给出了计算实例并与使用F K(频率波数 )域相移法、频率空间域有限差分法的结果进行了比较。与采用有限元的偏移方法相比 ,本方法可以节省大量内存 ;与采用有限差分的偏移方法相比 ,可以在一定程度上提高计算精度。本算法有可能在地震勘探数据处理中发挥一定的作用     

解四阶杆振动方程的三层高精度差分格式

对四阶杆振动方程构造含参数高精度三层差分格式,当参数满足一定条件时,差分格式稳定,局部截断误差阶数最高可达O(τ4+h8).数值例子说明该方法对稳定性的分析是正确的.