计算流体力学编程作业matlab,中科大计算流体力学CFD之大作业二.doc

中科大计算流体力学CFD之大作业二

CFD实验报告二

姓名：

学号：

一、题目

求解Poisson方程

，

， ，

， ，

， ，

描出等值线：, 0.2, 0.5, 0.75, 1.

要求所用方法：(1) Jacobi, G-S选一；SOR，线SOR，块SOR选一。迭代法要求误差；

(2) CG方法，MG方法选一。

二、报告要求

1)简述问题的性质、求解原则；

2)列出全部计算公式和步骤；

3)表列出程序中各主要符号和数组意义；

4)计算结果与精确解比较

5)结果分析(方案选择、比较、讨论、体会、建议等)；

6)附源程序。

三、问题简要分析

3.1问题性质

从题目给出的问题可以看出，该方程是一个二阶线性非齐次偏微分方程，并且给出了四个边界条件，更具定解条件，所以该方程是可以求解的。

3.2求解原则

题中是关于x和y的二阶导数，因此可以用二阶中心差分离散，即用正五点格式离散泊松方程。将差分方程整理成以五对角矩阵为系数的线性方程组，用Jacobi或者CG或者SOR等迭代方法求解线性方程组，即可得到函数的在网格点处的离散值。同时，计算域为的矩形区域，划分结构网格，均匀步长，方向网格步长，方向网格步长，则方向网格