对流扩散反应方程基于坐标变换的高阶紧致差分格式

基于非均匀网格,提出了一种求解一维定常对流扩散反应方程的高精度紧致差分格式。首先采用坐标变换方法将原方程由物理空间的非均匀网格转换为计算空间的均匀网格,然后给出一阶导数和二阶导数在均匀网格上的中心差分逼近式,并结合变换后的方程,得到了定常对流扩散反应方程具有四阶精度的紧致差分格式。最后,通过数值算例验证了该方法的精确性和高分辨率的特点。数值实验结果表明,对于所研究问题,该方法较不进行坐标变换而直接在物理域上建立的非均匀网格上的高阶紧致格式具有更高精度。     

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

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

二维Helmholtz方程的高阶紧致差分方法

基于二阶导数的四阶Padé型紧致差分逼近式,并结合原方程本身,得到了二维Helm-holtz一种四阶精度的紧致差分格式.该格式在每个空间方向上只涉及到三个点处的未知量及其二阶导数值,边界处对于二阶导数利用四阶显式偏心格式.然后,利用Richardson外推法、算子插值法及二阶导数在边界点处的六阶显式偏心格式,将本文构造的二维Helmholtz方程四阶紧致差分格式的精度提高到六阶.最后,通过数值实验验证了本文方法的精确性和可靠性.     

求解二维扩散方程的一种高精度紧致差分格式

扩散方程通常用来描述扩散现象中的物质密度的变化或者与扩散相类似的现象,针对二维扩散方程提出了一种高精度紧致差分格式,该格式基于四次样条函数对空间变量进行离散,对时间导数采用(2,2)Padé逼近,从而得到了时间和空间均为四阶精度的紧致差分格式.然后证明了该格式是无条件稳定的.最后通过数值实验,验证方法的精确性和稳定性.     

一维定常对流扩散反应方程的高精度紧致差分格式

针对一维定常对流扩散反应方程,提出了一种四阶精度的有理型紧致差分格式,其局部截断误差为O(h⁴);然后通过Richardson外推技术和算子插值法将本文格式的精度提高到六阶.因为格式仅涉及到3个网格基架点,所以对于Dirichlet边值问题,由差分格式可得三对角线性方程组,可采用追赶法进行求解.最后通过数值算例验证了本文方法的精确性和可靠性.     

三维对流扩散方程的高精度紧致差分格式

以一维定常对流扩散方程的高精度差分格式为基础,构造了三维非定常对流扩散方程的高精度紧致差分格式.该格式为两层格式,时间具有二阶精度,空间具有四阶精度.具体算例说明了上述格式的精确性和可靠性.     

污染物对流扩散方程的预测校正紧差分格式

结合预报校正线性多步法与高阶紧致差分格式方法的优点,空间导数采用四阶紧致差分格式进行离散之后,对得到的空间半离散格式采用改进的预报校正的线性多步法进行时间推进,得到一种时空方向均为四阶精度的求解非线性对流扩散方程的高精度方法。数值试验表明该格式可以有效求解非线性对流扩散方程,验证了格式的良好性能。     

定常对流扩散反应方程的指数型高阶差分格式

提出了一种数值求解一维定常对流扩散反应方程的指数型四阶差分格式．首先,由常数变易法求得模型方程的通解并应用到xi-1,xi,xi＋13点上,得到模型方程的指数型差分格式．然后,利用源项f（x）在点xi处的二阶泰勒展开,得到定常对流扩散反应方程的指数型四阶差分格式．最后,用数值算例验证了该格式的高阶精度和可靠性．     

一维四阶双曲方程的紧致差分格式

本文主要研究一维四阶双曲方程初边值问题.首先通过引入一个中间函数将其转化为二阶方程组,然后对方程中的空间导数项采用四阶紧致差分格式离散,时间导数项采用二阶中心差分格式离散,构造出问题的隐式紧致差分格式.数值算例表明该格式具有较好的计算效果.     

一维非定常对流扩散反应方程的高精度紧致差分格式

针对一维非定常对流扩散反应方程,首先推导了一种新的2层高精度紧致差分隐格式,其截断误差为O(τ~2+τh~2+h~4),即当τ=O(h~2)时,格式空间具有四阶精度;然后采用Fourier分析方法分析了格式的稳定性;最后通过数值算例验证了本文格式的精确性和可靠性.     

一维对流扩散方程的四阶精度有限差分法

文章基于Ｈｅｒｍｉｔｅ插值多项式的构造思路，推导出了一维含源扩散方程、扩散反应方程和无源项对流扩散方程的高精度紧致差分格式，并在导出三类特殊方程差分格式的基础上，建立了统一的含源项一维对流扩散方程的差分格式。本文方法是精确的，稳定的，且易应用到其它问题中。数值算例证明了所建立差分格式的有效性。     

一维对流扩散方程的四阶精度有限差方法

文章基于Ｈｅｒｍｉｔｅ插值多项式的构造思路，推导出了一维含源扩散方程、扩散反应方程和无湖泊项对流扩散方程的高精度紧致差分格式，并在导出三类特殊主程差分格式的基础上，建立了统一的含源项一维对流扩散方程的差分格式。本文方法是精确的。稳定的，且易应用到其它问题中，数值算例证明了所建立差分格式的有效性。     

数值求解一维波动方程的四阶紧致差分方法

针对一维波动方程提出了一种有限差分方法.首先,采用泰勒级数展开公式和原方程代入的方法推导出了第一个时间层未知函数值的四阶紧致差分格式.然后,用四阶紧致差分公式近似空间导数项,采用中心差分格式截断误差余项修正的方法处理时间导数项,推导出了第二个时间层以后未知函数的四阶紧致差分格式.该方法时间和空间具有整体四阶精度.利用Fourier方法分析了所提格式的稳定性.由于本文格式在未知时间层仅涉及3个网格点,因此可采用追赶法求解离散化后所得到的线性方程组.最后,用数值算例验证了本文格式的精确性和稳定性.     

二维扩散反应方程的三次样条高精度加权隐式差分格式及其多重网格方法

利用二阶微商的三次样条四阶紧致差分逼近公式,推导出两种数值求解二维扩散反应方程的两层9点加权隐式紧致差分格式.当θ=1/2时,该格式在时间和空间方向上分别达到二阶和四阶精度.通过Fourier方法讨论知,当1/2≤θ≤1时,格式是无条件稳定的;当0≤θ<1/2时,格式是条件稳定的.为了克服传统迭代法在求解隐格式方面的困难,差分方程采用多重网格方法进行求解并将本文格式的结果与P-R格式及C-N格式下的结果进行比较.数值实验结果验证本文方法的精确性和可靠性及多重网格方法的效率.     

一维对流扩散方程Crank-Nicolson特征紧致差分格式

本文针对一维线性对流扩散方程进行离散,在空间方向采用四阶紧致差分格式,对双曲部分采用时间二阶的Crank-Nicolson型特征差分格式,并在其中使用三次周期样条插值.数值算例表明该格式具有比较好的计算效果.     

椭圆型方程两点边值问题的混合型高精度紧致差分格式

【目的】针对一维椭圆型两点边值问题,发展一种六阶混合型高精度紧致差分格式。【方法】主要利用泰勒级数展开和组合紧致差分格式(Combined compact difference,CCD)的思想,将未知函数和它的一阶导数、二阶导数作为未知变量,利用函数和各阶导数之间的固定关系,将原方程对一阶导数泰勒级数展开式中产生的三阶导数项进行替换,同时也利用了一阶导数和二阶导数的六阶组合紧致格式。它的特点是显式紧致差分格式和隐式紧致差分格式混合在一起。【结果】最终使得混合型紧致差分格式整体达到了六阶精度。此外,提出的格式还具有推导简便,易实现编程,且能直接推广到高维问题的优点。尽管格式是六阶精度,但与四阶精度格式一样,空间方向仅仅需要3个网格点,因此由格式生成的方程组可采用追赶法进行高效求解。【结论】最后通过对具有精确解的4个算例进行数值实验,数值结果验证了该格式的精确性和可靠性。     

求解扩散方程的高精度显式紧致差分格式

首先针对一维扩散方程,空间方向采用二阶导数的四阶紧致差分公式进行离散,时间方向采用泰勒级数展开的方法进行离散,推导出了一种高精度显式紧致差分格式;然后通过Fourier分析方法给出了格式的稳定性条件为λ≤1/2(λ为网格比);最后通过数值实验验证了格式的精确性和可靠性.     

阻尼梁振动方程的高阶数值格式

本文研究了带有阻尼项的四阶梁振动方程初边值问题,基于紧致差分方法,给出了数值求解该问题的四种高阶紧致差分格式.对方程中的一阶和二阶时间导数项采用中心差分离散,对四阶空间导数项分别采用五点、七点和带紧致的五点、七点四种方法进行离散,得到四种高阶紧致差分格式,这四种格式均在时间方向达到二阶精度,在空间方向分别达到二阶、四阶...     

一维抛物型方程的ETF-FDS-四阶紧致差分-MG算法

将ETF-FDS格式和四阶紧致差分格式应用于一维抛物型方程,提出了ETF-FDS-四阶紧致差分-MG格式,用傅里叶方法证明该格式是无条件稳定的,并使用了多重网格法.最后用数值试验验证了方法的精确性与可靠性.     

构造对流扩散方程高精度非均匀网格格式的指数变换方法

提出了一种基于非均匀差分网格,构造求解对流扩散方程的高精度格式的指数变换方法.引入指数函数,将对流扩散方程变换为扩散反应方程,消除了数值求解中较难处理的对流项.采用优化差分方法推导出扩散反应方程基于非均匀网格的高精度差分格式,进而通过逆变换得到对流扩散方程的高精度格式.理论分析表明,该方法具有3至4阶精度,当计算区域为均匀网格时取得4阶精度.数值实例表明,在相同的非均匀网格系统中,此方法的计算精度明显优于传统的隐式差分方法.在水环境的实际模拟计算中,根据物理量的变化规律灵活地调整非均匀网格的间距,不仅能增强高精度差分方法的实用性,而且可以取得比均匀网格方法更为精确的计算结果.