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

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

基于二维泊松方程六阶紧致格式的多重网格方法

利用六阶紧致差分格式、结合多重网格V循环算法求解了二维泊松方程的Dirichlet边值问题，并用不同的松驰算子与四阶精度格式的多重网格方法进行了比较，计算结果表明，该方法在不明显增加计算量的前提下较四阶精度格式的多重网格方法具有更好的精确度和收敛阶，且ZLGS迭代不论对四阶精度还是对六阶精度格式的多重网格算法，都是一种较其他松弛算子更加有效的“光滑剂”。     

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

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

四阶紧致差分格式对三维网格计算特性的影响

为了研究四阶紧致差分格式对三维网格计算特性的影响,引进了一种推导斜压地转适应过程频散关系的通用方法,采用二阶中央差和四阶紧致差分格式从频率和群速度两个方面比较在各种三维网格上描述斜压地转适应过程产生误差的情况.结果表明,采用高精度的四阶紧致差分格式仅能提高三维网格C/LTS和EL/LTS的计算特性,而对其他几种三维网格的影响则不是正面的.四阶紧致差分格式在大气或海洋数值模式中应慎用.     

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

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

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

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

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

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

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

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

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

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

二维泊松方程非均匀网格上的高精度紧致差分格式

提出了数值求解二维泊松方程基于非均匀网格的高阶紧致差分格式,通过选取合适的网格分布参数求解具有边界层的数值算例,空间可以达到四阶精度.并与均匀网格上的计算结果进行比较,充分验证了本文非均匀网格高精度紧致格式的精确性和优越性.