Cattaneo模型的紧致有限差分法

紧致差分格式是一种高精度的有限差分方法.本文给出了Cattaneo模型的四阶紧致差分格式,通过对具体算例进行数值模拟,和二阶差分格式比较,验证了紧致差分方法的精确性和有效性.     

一维Zakharov-Rubenchik方程的有效紧致差分格式

偏微分方程的有限差分法是科学计算中的一种有效方法,采用经典的一阶和二阶有限差分格式对方程进行数值求解,要想得到较高精度的近似解是不容易的,一种合理的方法是设计高阶紧致差分格式.为了研究一维Zakharov-Rubenchik方程的有效紧致差分格式及其数值计算.针对一般形式的Zakharov-Rubenchik方程,提出了一种半隐式紧致有限差分格式,该格式克服了传统差分格式效率低、精确度不足的缺点,并在离散层次上保持了质量和能量的守恒性.最后,通过数值算例验证了该格式的精确程度及守恒性,并对几种不同差分格式的误差和计算耗时进行了比较,数值结果表明了半隐式紧致差分格式的高阶收敛性及有效性.     

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

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

极坐标系下求解一维Helmholtz方程的六阶紧致差分方法

基于极坐标系下四阶紧致差分方法,运用Taylor级数展开构造了一种极坐标系下求解一维Helmholtz方程的六阶紧致差分方法,通过分析所构造差分格式的截断误差确保该方法在理论上可达到六阶精度,最后通过数值算例验证了所构造差分格式的精确性和可靠性.     

三维双调和方程的高阶紧致差分格式及其多重网格方法

提出一类求解三维双调和方程的高精度紧致差分格式.该类格式是以泊松方程的高精度格式为基础的四阶精度19点紧致差分格式和六阶精度27点紧致差分格式.采用多重网格方法求解由高精度紧致差分格式所形成的代数方程组,并与低精度方法进行比较.讨论多重网格方法中不同松驰算子的迭代收敛效果.数值实验结果验证四阶紧致差分格式和六阶紧致差分格式的精度以及多重网格方法的可靠性和高效性.     

分数阶CEV模型下未定权益的紧致差分法

提出一种关于求解分数阶CEV模型下未定权益的紧致差分法.在时间上采用Caputo导数进行离散,在空间上采用4阶紧致差分格式进行离散.针对未定权益,得到一个时间2-α阶,空间4阶精度的紧致差分格式.并且运用傅里叶分析法和数学归纳法验证该方法的稳定性和收敛性.最后,通过数值实验验证该方法的有效性.     

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

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

求解热传导方程的四阶有限容积紧致格式

对非齐次热传导方程提出了一种数值求解的有限容积紧致格式,该格式具有空间上的四阶精度,且与有限差分紧致格式相比,更好地保持了问题的物理守恒性.数值算例表明,在相同的结点下,有限容积紧致格式比有限差分非紧致格式的精度更高.     

2维Maxwell方程的局部1维高阶紧致格式

将算子分裂方法与高阶紧致差分方法相结合,构造了2维Maxwell方程的局部1维紧致时域有限差分格式.该格式在时间方向和空间方向分别具有1阶和4阶收敛精度,并且具有计算效率高、无条件稳定的优点.数值实验表明:新构造的格式是能量守恒、高效率的.     

紧致差分格式求解二维非稳态不可压N-S方程

采用泰勒展式系数匹配的方法构造基于非等距网格的紧致差分格式并得出了它的截断误差.紧致差分格式能够很好的模拟不同时刻流场的变化情况,网格系统的选择对精度的影响很大,基于非等距网格的紧致差分方法是一种比经典差分方法精度更高的求解非稳态纳维斯托克斯方程的有效算法.     

一类四阶非线性抛物方程的紧致差分格式

针对一类四阶非线性抛物方程的初边值问题建立紧致差分格式,利用降阶的思想,通过引入中间变量将原四阶问题转化成二阶非线性方程组.对方程中的时间导数项和空间导数项分别采用Crank-Nicolson格式和四阶紧致差分格式进行离散,对非线性项采用外插的方法进行处理,从而得到原问题的三层线性紧致差分格式,其局部截断误差为■.数值算例表明该格式具有良好的计算效果.基于四阶非线性抛物方程在薄膜理论等问题中的重要作用,对此类方程构造高精度的紧致差分格式,可以使该方程在有关工程计算方面得到更好的应用,因此该研究成果具有重要的理论意义和广泛的应用前景.     

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

本文研究了带有阻尼项的四阶梁振动方程初边值问题,基于紧致差分方法,给出了数值求解该问题的四种高阶紧致差分格式.对方程中的一阶和二阶时间导数项采用中心差分离散,对四阶空间导数项分别采用五点、七点和带紧致的五点、七点四种方法进行离散,得到四种高阶紧致差分格式,这四种格式均在时间方向达到二阶精度,在空间方向分别达到二阶、四阶、四阶和六阶精度.最后利用数值算例验证了四种格式的精度阶与理论结果一致.本文相对于之前的研究,对弹性梁的振动增加了阻尼因素,因此也更加适合对实际问题的数值计算.     

弹性波方程的紧致差分方法

在对弹性波方程进行数值模拟时 ,低阶差分格式往往产生严重的数值频散 ,高阶显示差分格式需要用较多的网格点 ,不利于边界的处理。而紧致差分格式吸收了它们的优点 ,弥补了它们的不足。为此该文应用紧致差分格式的思想 ,发展了二维情况下弹性波方程初值问题的紧致差分方法 ,研究了它的稳定性 ,并用 Fourier方法分析了显示差分格式和紧致差分格式的相速度误差 ,最后利用紧致差分方法在粗网格条件下对地震波传播进行了数值模拟 ,并同五点四阶中心差分方法的计算结果进行了对比。结果表明 ,求解弹性波方程的紧致差分方法有效 ,且具有比同网格点差分格式更高的计算精度和较小的数值频散。     

求解二维热传导方程的高精度紧致差分方法

基于Richardson外推法提出了一种数值求解二维热传导方程的高阶紧致差分方法.该方法首先利用时间二阶、空间四阶精度的紧致交替方向隐式(ADI)差分格式在不同尺寸的网格上对原方程进行求解,然后利用Richardson外推技术外推一次,最终得到了二维热传导方程时间四阶、空间六阶精度的数值解,数值实验验证了该方法的高阶精度及有效性.     

求解一类非线性偏微分方程的高精度紧致差分方法

针对一类非线性偏微分方程,提出了一种新的高精度紧致差分方法.首先对内部网格节点处的空间一阶和二阶导数项采用四阶精度的Padé紧致差分格式进行离散,然后对时间导数项采用泰勒级数展开并使用截断误差余项修正法进行离散,最终得到了求解该非线性方程的一种三层隐式高精度紧致差分格式,其截断误差为O(τ2+τh2+h4),即当τ=O(h2)时,该格式在空间上具有四阶精度.最后通过对广义Burgers-Fisher方程和广义Burgers-Huxley方程的数值求解,验证了本文方法的精确性和可靠性.     

二维Helmholtz方程的改进六阶紧致差分法

为得到求解二维Helmholtz方程的高精度差分法, 构造了一种改进六阶紧致差分格式: 首先, 给出一种带优化参数的六阶紧致差分格式的截断误差; 然后, 对此截断误差的部分项进行二阶紧致逼近, 得到一种改进紧致差分格式; 其次, 对该格式进行了收敛性分析, 证明其为六阶收敛的; 最后, 基于极小化数值频散的思想, 给出该格式优化参数的加细选取策略。与带优化参数的六阶紧致差分格式相比, 数值实验说明改进六阶紧致差分格式的数值精度有了显著提高, 且其误差对波数k的依赖性更低。     

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

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

非线性薛定谔方程的几种差分格式

在满足一定的初值、边值条件下,结合不同的差分格式对非线性薛定谔(NLS)方程进行数值求解.分别利用经典的向前差分算子、二阶中心差分算子、Crank-Nicolson方法和紧致差分算子构造向前Euler格式、Crank-Nicolson格式和紧致差分格式,并证明Crank-Nicolson格式和紧致差分格式精确保持离散质量守恒和能量守恒.利用数学软件MATLAB进行实验计算,结果表明:所构造的3种格式具有合理性及有效性.     

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

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

三维泊松方程基于Richardson外推的高精度多重网格解法

基于三维泊松方程的四阶紧致差分格式,利用Richardson外推法、算子插值法和多重网格算法,使已有四阶紧致差分格式的计算精度整体提高二阶,精度达到六阶.数值实验验证六阶格式的精确性和多重网格方法的有效性,并与四阶紧致差分格式多重网格方法的计算结果进行比较.