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

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等迭代方法求解线性方程组，即可得到函数的在网格点处的离散值。同时，计算域为的矩形区域，划分结构网格，均匀步长，方向网格步长，方向网格步长，则方向网格点数为m=， 方向网格点数为n=。 四、计算公式和步骤 4.1