高维分数阶cable方程隐式差分逼近

针对高维分数阶cable方程的数值差分逼近问题,采用有限体积方法,构造高维分数阶cable方程一种隐式差分逼近格式.结果表明:该隐式差分格式是无条件稳定和收敛的.利用隐式差分方法求解三维情况的数值例子,将数值解与精确解进行比较,说明隐式差分方法的有效性.此方法可应用于其它类型的高维分数阶微分方程.     

椭圆界面问题的高阶差分格式

构造混合边界条件下椭圆界面问题的一个高阶数值格式.在求解区域内部及界面处采用四阶逼近,边界处采用三阶数值格式,得到一个整体四阶精度的求解格式.数值实验证明了格式的高精度及有效性.     

改进的有限体积法扩散项离散格式

分析了传统的有限体积法求解扩散问题中,当扩散系数在有限控制体界面两侧差异较大时,扩散通量在界面附近出现的数值振荡现象.根据界面两侧扩散通量相等和有限控制体内变量线性变化假设,推导了改进的有限体积法扩散项离散格式.数值算例表明：由2种扩散项离散格式计算得到的变量及其梯度分布很接近,而改进的扩散项离散格式消除了扩散通量的振荡问题,计算结果更加准确,能应用于扩散系数剧烈变化时的扩散问题求解.     

四边形剖分下一类对称有限体格式的构造

本文针对一类自共轭椭圆问题,在四边形网格剖分下,通过选取“verterx-centred”型控制体,利用有限体积元方法,构造了一类保对称的有限体格式,并将之转化为Galerkin变分形式.数值实验表明其误差均达到了饱和阶.由于格式对称,相应离散化系统可采用PCG等方法快速求解,从而提高了运算效率.另外,新格式对flux(流)函数具有超逼近性.     

一种求解变波数Helmholtz方程的高精度紧致差分方法

【目的】进一步研究Helmholtz方程对于大波数和变波数问题的数值计算,数值求解Helmholtz方程具有重要的理论价值和现实意义。【方法】利用泰勒级数展开,并结合混合型紧致格式的思想,推导了数值求解一维和二维Helmholtz方程的六阶精度紧致差分格式。并且格式涉及到未知函数及其一阶和二阶导数值,为保证格式的整体精度,对一阶和二阶导数的计算也采用六阶紧致差分格式。【结果】格式在小波数和变波数的情况下都有六阶精度,在大波数的情况下仍然能保持三阶以上精度。【结论】数值实验验证了格式的精确性和可靠性。     

一维Helmholtz方程的四阶紧致有限体积方法

本文分别提出了一维Helmholtz方程基于Dirichlet和周期边值问题的四阶紧致有限体积方法.对于Dirichlet边值问题,通过Taylor展开给出了方程的四阶紧致有限体积格式,并结合边界处的四阶近似,证明了此问题的离散格式是四阶格式.对于周期边值问题,利用周期边界条件,同样得到了此问题的四阶紧致有限体积格式.数值实验表明给出的两种格式均是四阶格式.     

求解一类多项四阶时间分数阶慢扩散系统的有限差分格式

介绍求解多项四阶时间分数阶慢扩散方程的有限差分方法.利用L1公式逼近时间分数阶导数,用降阶法处理空间四阶导数项,再借助离散能量方法证明差分格式是无条件稳定的且在无穷范数下其收敛阶为O(τ2-β+h2),其中τ和h分别为时间方向和空间方向的步长,β是时间分数导数的最大阶.最后用数值实验验证所提出差分格式的精度和有效性.     

有限体积法定价欧式跳扩散期权模型

考虑有限体积法求解Kou跳扩散期权定价模型.基于线性有限元空间,构造了向后Euler和Crank-Nicolson两种全离散有限体积格式,并结合简单高效的递推公式逼近方程中的积分项.理论分析表明所得的离散矩阵为M-矩阵.数值实验验证了方法的有效性.     

求解一类多项四阶时间分数阶慢扩散系统的有限差分格式

介绍求解多项四阶时间分数阶慢扩散方程的有限差分方法.利用L1公式逼近时间分数阶导数，用降阶法处理空间四阶导数项，再借助离散能量方法证明差分格式是无条件稳定的且在无穷范数下其收敛阶为O(t^(2-B)+h^2)，其中t和h分别为时间方向和空间方向的步长，B是时间分数导数的最大阶.最后用数值实验验证所提出差分格式的精度和有效性.     

四阶线性方程的紧致体积格式及误差分析

本文主要研究四阶线性方程的紧有限体积格式.在四阶问题中引入一个中间函数将其转化为二阶方程组,对其中的方程分别利用四阶紧致体积方法来处理,得到紧致体积格式,并对格式进行了误差分析.数值算例表明该格式具有比较好的计算效果.