二维possion方程c语言实现,二维Possion方程两步预估校正格式.pdf

二维Possion 方程的两步预估校正格式

1,2 3 3 3

孙亮 ，马东军，秦丰华，孙德军

(1．中国科学技术大学 地球和空间科学学院 安徽合肥 230026 ；

2 ．中国科学院大气物理研究所，LASG 国家重点实验室，北京 100029 ；

3 ．中国科学技术大学 力学和机械工程系 安徽合肥 230026 )

摘要：本文提出了一种求解Possion 方程的两步预估校正格式，它用两次二阶格式的计算量，达

到四阶精度。与传统的四阶中心差分格式和四阶紧致格式不同，新格式以传统五点二阶中心差

分格式的解为预估解，通过增加对残差的校正计算，使得最终解具有四阶精度。与普通的五点

格式相似，新格式的边界条件处理简单，易于编程实现，并且很容易推广到其它复杂情形。通

过数值实验，结果证明上述格式的确具有编程和计算简单、精度高的优点。预估校正格式也为

改造低精度格式提供了一种很好的思路。

关键词：Possion 方程，预估校正格式，四阶精度

1、引言：

Possion方程的数值求解常常出现在物理、力学等科学计算问题中，特别是在计算不可压缩流

动问题中，总需要在每一步求解压力Possion方程。Possion方程的快速求解和高精度计算成为提

高数值格式整体精度和速度的一个关键因素。传统的五点中心差分格式计算量小，且具有二阶

精度，在科学计算和工程中得到了广泛的应用。近年来，随着科学和工程计算的需要，人们对

[1-3]

计算精度和速度的追求越来越高。这些需求推动着Possion方程求解算法的不断发展 。

一方面为了提高数值精度，使用较少网格点所构造的紧致格式得到了广泛的应用[4-7] 。另一

方面，多重网格方法的出现，极大地提高了Possion方程的求解速度[1,3] 。两者的结合，能有同时

效提高数值精度和计算速度，因而在近年来也得到了大量的应用[2-7] 。虽然紧致格式加多重网格

与传统的五点中心差分格式相比具有明显优势，但是在应用时也存在一些的不便。由于传统的

紧致格式构造比较复杂，而且需要有复杂的边界条件处理以保证达到相应的精度。多重网格的

思想虽然简单，但是实际的插值算子和投影算子构造也很复杂。当多重网格与紧致格式相结合

时，程序编写变得更困难。

本文通过构造新的四阶格式，旨在避免这一困难。我们的方法是将整个过程分为两步：首

先是求解出一个预估值，然后通过对残差的求解，在校正步中求得具有四阶精度的结果。在每

一步求解中我们都采用传统的五点格式，这样可以继承那些求解五点格式Possion 方程的快速算

1 本文受国家自然科学基金(No No.)和LASG 国家重点实验室开放基金资助

通讯作者：孙亮，中国科学技术大学地球与空间科学学院 0551-3600175 Email: sunl@

1

法，从而同时解决速度和精度问题。

2 、算法

2 ．1 格式

u u u

i-1,j+1 i,j+1 i+1,j+1

u