首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
非线性发展方程的小模板简化Padé格式   总被引:3,自引:0,他引:3  
在有理逼近的紧致格式的理论基础上,采用特别的统一的Pad啨逼近形式,构造了针对高阶非线性发展方程的、简单小模板的差商格式· 不仅保持了格式的四阶精度,而且还可以采用追赶法求解得到的3对角矩阵,或者采用三阶Runge_Kutta法直接求解积分· 计算效果通过多种算例表明是十分令人满意的· 相对于其他差分格式,此方法具有模板较小而精度保持四阶的优点·  相似文献   

2.
推导了分数阶积分的梯形逼近格式以及Caputo导数的L1逼近格式的四阶展开公式.并利用L1格式的展开式得到了Caputo导数的具有3-α阶精度的三点逼近格式,该逼近格式被应用于数值求解分数阶松弛方程和时间分数阶次扩散方程.  相似文献   

3.
不可压混溶驱替问题的流线扩散──混合元数值模拟   总被引:2,自引:0,他引:2  
采用标准元模拟不可压混溶流问题,当扩散系数矩阵小过剖分参数时,有限元格式仅能给出比最优精度低一阶的逼近解,格式稳定性差并伴有强烈的数值弥散现象.为了克服上述缺陷,本文对压力方程采用混合元,而对浓度方程采用流线扩散格式,在扩散矩阵为线性的假定下,证明了该格式具有较标准元更高的逼近精度(比最优阶低1/2)和更好的稳定性.  相似文献   

4.
以四阶CWENO重构为基础,通过将对流项采用低耗散中心迎风格式离散,扩散项采用四阶中心差分格式离散,对得到的半离散格式采用四阶龙格库塔方法在时间方向上推进,得到一种求解对流扩散方程的高阶有限差分格式.数值结果验证了该格式的四阶精度和基本无振荡特性.  相似文献   

5.
为了降低经典的三阶加权本质无振荡(WENO)格式的数值耗散,提出了一种新的三阶WENO格式的修正模板近似方法.改进了经典WENO-JS格式中各候选模板上数值通量的一阶多项式逼近,通过加入二次项使模板逼近达到三阶精度.计算了相应的候选通量,并且通过引入可调函数φ(x),使得新的格式具有ENO性质.最后给出了一系列数值算例,证明了该方法的有效性.  相似文献   

6.
崔吉田  王同科 《应用数学》2012,25(1):96-104
本文针对常系数和变系数两点混合边值问题提出一种紧有限体积格式,该格式形成的线性代数方程组具有三对角性质,可以使用追赶法求解.证明格式按照H1半范数具有四阶收敛精度.利用节点计算值,给出单元中点值和一阶导数值的高精度后处理计算公式,这两个公式同样具有四阶精度.数值算例验证了理论分析的正确性,并说明了格式的有效性.  相似文献   

7.
任全伟  庄清渠 《计算数学》2013,35(2):125-136
针对研究吊桥模型而建立的四阶微积分方程, 提出Legendre谱逼近法进行求解.构造迭代算法来求解得到的线性系统, 证明了迭代格式的收敛性, 对问题进行了误差分析.数值算例验证了迭代的收敛性和方法的高精度.  相似文献   

8.
色散方程的一类新的并行交替分段隐格式   总被引:14,自引:0,他引:14  
王文洽 《计算数学》2005,27(2):129-140
本文给出了一组逼近色散方程的非对称差分格式,并用这组格式和对称的Crank-Nicolson型格式构造了求解色散方程的并行交替分段差分隐格式.这个格式是无条件稳定的,能直接在并行计算机上使用.数值试验表明,这个格式有很好的精度.  相似文献   

9.
线性正算子只能达到一阶代数精度,通过对线性算子进行组合,可以得到具有高阶代数精度的逼近算子.本文给出保持x3不变的线性逼近算子的具体形式.  相似文献   

10.
郑宁  殷俊锋 《计算数学》2013,35(3):275-285
本文讨论基于不光滑边界的变系数抛物型方程求解的高精度紧格式.首先构造一般变系数抛物型方程的高精度紧格式,并在理论上证明格式具有空间方向四阶精度.然后针对非光滑边界条件,引入局部网格加密技巧在奇异点附近进行不均匀的网格加密.数值实验以期权定价中Black-Scholes偏微分方程的求解为例,验证高精度紧格式用于光滑边界条件的微分方程离散可以达到四阶精度.对于处理非光滑边界条件,网格局部加密技巧能有效的提高数值解精度,使得高精度紧格式用于定价欧式期权可以接近四阶精度.  相似文献   

11.
不增加基点,仅摄动二阶ENO格式的系数(简记为MCENO),得到一类求解双曲型守恒律方程的三阶MCENO格式.由MCENO格式的构造过程可以看出,MCENO格式保留了ENO格式的许多性质,例如本质无振荡性、TVB性质等,且能提高一阶精度.进一步,利用MCENO格式模拟二维Rayleigh-Taylor(RT)不稳定性和Lax激波管的数值求解问题.数值结果表明,t=2.0时,MCENO格式的密度曲线处于三阶WENO格式和五阶WENO格式之间,是一个高效高精度格式.值得注意的是,三阶MCENO格式,三阶WENO格式和五阶WENO格式的CPU时间之比为0.62:1:2.19.表明相对于原始ENO格式,MCENO格式在光滑区域有较高精度,能提高格式精度.  相似文献   

12.
提出了数值求解三维热传导方程的一个四阶精度的有限差分格式,首先对三个空间方向上的二阶导数项,采用四次样条函数来近似,从而得到半离散的常微分方程.然后利用常微分方程的解析解表达式,时间矩阵利用Padé近似,得到时间和空间均为四阶精度的差分格式.最后利用方法计算了两个数值算例,并与文献中结果进行了对比,从而验证了高精度格式的性能.  相似文献   

13.
提出了数值求解一维非定常对流扩散反应方程的一种高精度紧致隐式差分格式,其截断误差为O(τ~4+τ~2h~2+h~4),即格式整体具有四阶精度.差分方程在每一时间层上只用到了三个网格节点,所形成的代数方程组为三对角型,可采用追赶法进行求解,最后通过数值算例验证了格式的精确性和可靠性.  相似文献   

14.
借助显式紧致格式和隐式紧致格式的思想,基于截断误差余项修正,并结合原方程本身,构造出了一种求解一维定常对流扩散反应方程的高精度混合型紧致差分格式.格式仅用到三个点上的未知函数值及一阶导数值,而一阶导数值利用四阶Pade格式进行计算,格式整体具有四阶精度.数值实验结果验证了格式的精确性和可靠性.  相似文献   

15.
基于复化Simpson公式和复化两点Gauss-Legendre公式,构造了两个求解时间分布阶扩散方程的高阶有限差分格式.不同于以往文献中提出的时间一阶或二阶格式,这两种格式在时间方向都具有三阶精度,而在分布阶和空间方向可达到四阶精度.数值结果表明,两种算法都是稳定且收敛的,从而是有效的.两种格式的收敛速率也通过数值实验进行了验证,并且通过和文献中的算法对比可以得出其更为高效,  相似文献   

16.
张静  王彩华 《应用数学》2022,(2):355-373
本文针对含扰动参数ε的含源反应扩散方程,采用待定系数法,在三点模板的中心点处进行泰勒展开,对泰勒展式中的高阶导数项充分利用原微分方程进行"降阶",然后分别从"横向"和"纵向"两个角度进行修正,得到了两类差分格式,其中横向系列差分格式(HDS)的精度分别达到二阶、四阶和六阶.数值实验与参考格式比对效果较好,且横向差分格式...  相似文献   

17.
基于紧致差分方法,推导出一个时间和空间方向均为四阶精度的三层隐式紧致格式,并采用Fourier分析法给出了格式的稳定性条件.最后,数值例子验证了所给出来的格式的精度和可靠性.  相似文献   

18.
三维泊松方程的高精度多重网格解法   总被引:7,自引:0,他引:7  
利用对称网格点泰勒展开式中各阶导数项明显的对称性,得到了数值求解三维泊松方程的四阶和六阶精度的紧致差分格式,其推导过程简便直接.为了克服传统迭代法在求解高维问题时计算量大、收敛速度慢的缺陷,采用了多重网格加速技术,设计了相应的多重网格算法,求解了三维泊松方程的Dirichlet边值问题.数值实验结果表明,本文所提出的高精度紧致格式达到了期望的精度并且多重网格方法的加速效果是非常显著的.  相似文献   

19.
采用标量辅助变量(scalar auxiliary variable, SAV)方法结合重心插值配点法求解二维Allen-Cahn方程.在时间方向上分别采用Crank-Nicolson格式、二阶向后差分格式离散,空间方向上采用重心Lagrange插值配点法离散,建立了两种无条件能量稳定SAV格式,并给出了重心插值配点格式的逼近性质.数值实验表明:两种SAV配点格式的时间收敛阶为二阶,并满足能量递减规律.与空间采用有限差分法离散对比,重心Lagrange配点格式具有指数收敛的特性.  相似文献   

20.
侯波  葛永斌 《应用数学》2019,32(3):635-642
本文提出数值求解一维对流方程的一种两层隐式紧致差分格式,采用泰勒级数展开法以及对截断误差余项中的三阶导数进行修正的方法对时间和空间导数进行离散.格式的截断误差为O(τ~4+τ~2h~2+h~4),即该格式在时间和空间上均可达到四阶精度.利用von Neumann方法分析得到该格式是无条件稳定的.通过数值实验验证了本文格式的精确性和稳定性.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号