上海交通大学硕士研究生课程《传热流动的数值分析大作业》.docVIP

“传热流动的数值分析”2015年大作业 2维条件下的无粘、不可压缩流体通过出口和入口流过箱体，具体情况如图所示，求该箱体内的流线情况，所有单位为厘米。 （1）流线方程为： 使用Gauss-Seidel线迭代，，误差，结果输出中，包括在y=0,1,2,3,4,5 处的所有X处对应的流函数值。 （2）设出口处纵向速度V =0，试采用PSOR方法，，计算在x=0,1,2,3,4,5 处的所有Y处对应的流函数值，以及不同的松弛系数和迭代次数的关系曲线（至少三个系数）。 答：（1）该问题为稳态问题，流线方程为椭圆型方程，在求解方程时，首先对方程在计算域内进行离散。计算域为：，在离散时，，因此可以得到流线方程的差分方程为： （1） 整理后可得： （2） 在本题中，采用Gauss-Seidel线迭代方法进行求解，扫描方向选为自左向右，此时有 （3） 由于是线推进，因此在当前线方向求解时，之前的扫描线上的参数已经得到更新，所以方程可改写为： （4） 其中看着当前迭代层中的已知变量。 从方程（4）中可以明显看出，在扫描线方向上，离散方程组的系数矩阵为三对角矩阵，因此该方程组可以采用追赶法进行求解。 从左向右如此推进，在做完了全区内各列的求解后就完成了一轮迭代，具体的程序代码如附录1所示。执行程序请参见文件包。 对于该方法，进出口边界的流函数假设为，对于进口，由于，因此。而对于出口边界，由于，且进口流量相等，因此可以获得。 计算输出在y=0,1,2,3,4,5 处的所有X处对应的流函数值，具体数值请参见文件夹PSOR中执行程序的输出文件RESULTS.TXT，为了方便观察计算结果，本文采用Tecplot作出云图，通过云图可以清晰的观察到计算域内的流线，如图1所示。图中每条网格线的交点即为流函数值。 此方法的迭代次数为645次。 图1 Gauss-Seidel线迭代方法求解结果 （2）本题采用PSOR方法进行计算。在该方法中，需要对进行预计算，本文采用Gauss-Seidel点迭代对该值进行预计算，即： （5） 在本文中，Gauss-Seidel点迭代实施步骤方式为先对同一条X线逐点计算，然后再从左向右进行更新。在程序计算时，方案实施非常简便，在同一轮迭代中，将计算域内所有的流函数值采用同一个矩阵进行存储，当计算更新某个点时，采用矩阵赋值的方法及时进行更新。 PSOR方法的计算方程为： （6） 和题（1）相似，边界条件也采用相同的，此外，出口纵向速度V=0进一步确定。 在x=0,1,2,3,4,5 处的所有Y处对应的流函数值通过云图进行了展示，如图2所示。此外详细的数值结果也可以在该题的程序文件夹PSOR中的RESULTS.TXT中获得，在此不做详细列出。具体程序代码如附录2所示。 图2 PSOR方法求解结果 此外，本文采用了7个不同的松弛系数进行了计算，分布是1.2,1.4,1.5,1.6,1.8,1.9,2.0,其中当松弛系数为2.0时，计算结果不再收敛，其余计算结果如图3所示。随着松弛系数的增大，迭代次数减小。 图3 不同松弛系数下的迭代次数 试使用一阶波动方程，计算波在一个封闭管道内从t=0s到t=0.15s的传播过程。假设声速V=200m/s，在t=0s的初始波形为如下图所示的三角波，请分别用以下方法和步长求解。 一阶上风法 Lax-Wendroff BTCS隐式法 使用以下的不同的步长： （a） （b） （c） 请以0.025s为间隔图示0s到0.15s的波传递情况。 答：一阶波动方程为： （7） 按照题目要求，对上述波动方程进行离散可以获得相应的离散方程并进行求解。 一阶上风法 （8） Lax-Wendroff （9） BTCS隐式法 （10） 从上面的方程中可以看出，方程（8）和（9）