喷管网格生成范文.docVIP

喷管网格生成报告 一、 PDE方法建立坐标转换 假定物理平面为，计算平面为，物理平面和计算平面之间的Jacob矩阵 为，为了满足两平面之间一一对应的关系，应该是非平凡的。 利用微分计算的链式法则 ，可以得到 将其代入Laplace方程中，并做和运算，容易得到如下逆变换方程， 其中， ， ，此为非线性椭圆形方程。 二、 逆变换PDE的差分方程 根据差分知识，对于椭圆形方程，其离散方法宜采用中心型差分格式，因此差分方程是隐式的。x，y的各阶导数的中心差分格式为， 将上述式子代入到逆变换方程的系数中， 代入逆变换方程中，并令得逆变换PDE的差分方程， 三、 边界条件 几何边界 收敛圆弧： 收敛錐： 喉部圆弧： 扩散錐： 直线入口边界： 直线出口边界： 圆弧入口边界： 圆弧出口边界： 四边界节点边值的确定 (1) 直线进出口边界时，喷管壁及喷管轴上的逆变换为： 其中为喷管壁入口处的轴向坐标，为喷管壁出口处的轴向坐标，A为调节参数。 进出口边界等距分割，即等参变换： (2) 圆弧进出口边界时，以喉部为界，设喷管壁在喉部上游最临近喉部的节点轴向坐标为，圆弧入口、出口边界的轴向坐标为和，则喉部上游的喷管轴上的逆转换为： 喉部下游的喷管轴上的逆转换为： 喷管壁上的逆变换为： 左右进出口边界等弧分割： 进口圆弧： 出口圆弧： 四、计算方法 简单迭代法（Jacobi法） Gauss-Seidel迭代法 迭代时，列按从左至右，行按从下到上的顺序。当计算时，总是启用前面最新计算出的，，，，计算时也如此， 逐次松弛迭代（SOR法） SOR迭代格式是第k+1次迭代的结果看成是第k次迭代的结果加上一个校正值乘上一个松弛因子，使得改进后的迭代方法收敛速度较快。 这里将Gauss-Seidel迭代值作为，再将与加权平均，得 这里采用超松弛，取 线（块）迭代法——每次用直接法（解三对角矩阵方程）求解一行或一列未 知数 (1) Jacob法 逆变换PDE差分离散方程经过移项可化为三对角矩阵方程如下（这里是一行一行的求解）： 即 用TDMA（追赶法）求解上面三对角矩阵方程 (2) Gauss-Seidel线迭代 即 (3) SOR法 将Gauss-Seidel线迭代值作为，再将与加权平均，得 交替方向隐式迭代（ADI法＋SOR法） 第一步横扫 然后采用的超松弛 第二步竖扫 然后采用的超松弛 迭代初值选择 采用轴对称形式 五、程序框图 程序结构见下页，原程序见附件的F90文件。 .TRUE. .FALSE. .TRUE. .FALSE. .TRUE. .FALSE. .TRUE. .FALSE. ︰ ︰ .TRUE. .FALSE. 六、计算结果及结果分析讨论 本次结果分析取Imax=41，Jmax=11，收敛阀值为0.0001。 E-005 E-005 迭代次数 占用CPU时间（s） E-001 A＝0.5 直线进出口边界 点迭代Jocobi法 3.3704087618579.951135831842 862 3.90625 点迭代G-S法 9.81609022832 13.16945903386 446 2.18750 点迭代SOR法 9.990048881532.992059048478 180 0.93750 线迭代Jocobi法 9.816501737797 2.559