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

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

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

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

封闭方腔内自然对流问题的高精度紧致差分格式

提出数值求解二维非定常不可压涡量-流函数Navier-Stokes/Boussinesq方程组的高精度紧致差分格式,格式空间为四阶精度,时间为二阶精度,并且是无条件稳定的.为了验证高精度紧致差分格式的精确性和可靠性,对有解析解的二维非定常不可压Navier-Stokes/Boussinesq方程组的Dirichlet问题和典型的封闭方腔自然对流问题进行数值模拟.     

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

提出了三维Helmholtz方程等距网格上的一种四阶精度19点紧致差分格式。结合多重网格V循环算法和红黑高斯-塞德尔迭代法进行求解,并与二阶中心差分格式进行了比较。计算结果验证了本文方法的精确性和有效性。     

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

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