计算热物理大作业报告.docVIP

计算热物理大作业报告

S陆振能

PDE方法建立坐标转换的根底及控制方程推演

设代表物理平面，代表计算平面。网格生成方程为物理平面上的Laplace方程，即

，(1)

为了进行离散化数值求解，将方程(1)变换到计算平面上为

式中

， ，

进行离散时，上式二阶导数采用二阶中心差分，一阶导数，交错导数，采用两步中心差分得

将上式代入椭圆型PDE

可得

考虑当时的特殊情形，那么有

同理可得

这样我们就得到了离散了的Laplace形式。

定解条件〔边界条件〕确实定——几何分析

喷管型面为：收敛圆弧－收敛锥－圆弧－圆弧－扩散锥；入口边界和出口边界特征：

均为垂直于轴向X轴的直线；

均为一段圆弧。

为使喷管喉部附近网格较密，喷管壁上的转换坐标由下式确定：

即逆变换为：

其中：为喷管壁入口处的轴向坐标，为喷管壁出口处的轴向坐标，A为调节参数，用于调节网格疏密，A越大，喉部网格越密，分别取0.5，1.0，1.5。

喷管轴线上，以喉部为界，要保证喉部上下游节点数与喷管壁喉部上下游的节点数一致。为此，对垂直于X轴的进出口边界，轴上坐标转换与〔5〕式相同，对圆弧形的出入口边界，以喉部为界，按变化了的〔5〕式分两段来计算：设喷管壁在喉部上游最接近喉部的节点轴线坐标为,圆弧入口、出口边界的轴向坐标为和，那么喉部上游的喷管轴上转换为：

而喉部下游的喷管轴上转换为：

进出口边界按其规定的两种特征，分别取等距分割和等弧分割。

当进出口边界为垂直于轴向的直线时，取等距分割，入口边界各点的轴向坐标为；径向坐标为入口上缘径向坐标按径向网格数均分得到。

当进出口边界为一段圆弧时，取等弧分割，入口边界上各点的轴向坐标为，其中为入口边界上缘径向坐标，d为圆弧段轴向长度，dx为该点到入口边界下缘点的轴向距离，可由入口段圆弧的轴向和径向长度求得圆心角及半径后得到；径向坐标也可由入口段圆弧的轴向和径向长度求得圆心角及半径后得到。

根据喷管壁的形状在其上取6个特征点，即入口边界上缘、收敛圆弧终点、收敛锥终点、喉部上缘、扩散锥起点、扩散锥终点。

计算方法——迭代法及初值选择

3.1JACOBI点迭代

式中，。

3.2Gauss-Seidel点迭代

这种方法与jocobi迭代类似，由上图可以看出，如按从下到上逐行扫描，那么求x(i,j)时，x(i-1,j-1),x(i-1,j),x(i-1,j+1),x(i,j-1)都已求出。y也是同样情况。为了及时利用算得的新值，得G-S迭代格式如下：

这里已结合题目条件，令

3.3逐次松驰迭代点SOR

SOR法分为两步：

利用G-S迭代，求出，；

将G-S迭代结果与上次迭代结果进行加权平均。

3.4JACOBI线迭代

对于x

式中

，，

对于，只要将上式中的相应地变为就可以了。

式中

，，

3.5Gauss-Seidel线迭代

对于x

式中

，，

也就是但凡下标中出现j-1的值均用新值代替。

对于，只要将上式中的相应地变为就可以了。

式中

，，

3.6逐次松驰SOR线迭代

SOR法分为两步：

用G-S迭代，求出，；

G-S迭代结果与上次迭代结果进行加权平均。

3.7交替方向隐式迭代

先逐列，再逐行扫描，同时在G-S迭代中，除xi,j外都乘以。

假设按照j从小到大的次序逐列迭代，按照i从小到大的次序逐行迭代那么有

对于x

1〕逐列迭代

式中

，，

2〕逐行迭代

式中

，，

对于y，只要将相对应处的x换成y，就可以了。

3.8迭代边值的选取

喷管壁上的对应的转换坐标由下式确定

(1)

即逆变换为

(2)

由上式可以得到上边界上的x值，上边界的y=f(x)，这里f(x)是描述上边界几何形状的一分段函数。

喷管轴线上，以喉部为界，要保证喉部上、下游节点数与喷管壁喉部上、下游的节点数一致。为此要对(1)式作些变化，将其分为两段来计算：设喷管壁在喉部上游最邻近喉部的节点轴向坐标为，圆弧入口、出口边界的轴向坐标为和，那么喉部上游的喷管轴上转换为

(3)

而喉部下游的喷管轴上转换为

(4)

由(3)与(4)的逆转换即可得下边界的x,y值。

图1为物理平面图。

图1物理平面图

进、出口边界如为垂直于轴线的取等距分割，如为圆弧进出口条件那么取等弧分割。

图2为计算平面图，间距x，y方向多取1，x方向取41个网格点，y方向取11个网格点。

图2计算平面图

3.9迭代初值的选取

考虑到喷管是轴对称的，故初值取如下值：

程序框图

开始

开始

对x,y的边值进行赋值

确定x,y的初值

迭代求值

满足精度吗？

结果输出

结束