数值传热学大作业3gg.docVIP

数值传热学大作业 （第3题） 组长：王晏林 组员：禹国军 王 岩 刘文婷 韩立涛 完成日期：2010年01月03日数值传热学2009-2010学年第一学期大作业3 阶梯形标量场（长宽均为1）的纯对流传递（），控制方程为，上游的边界条件都是第一类的，即给定了的分布，下游按开口边界的方式处理。 图1 示意图 要求： 1、分别利用FUD,CD,SUD,QUICK,CLAM,EULER,MINMOD,MUSCL,OSHER,SMART, STOIC来离散对流项，观察它们的计算结果有何不同。 2、用延迟修正进行求解。 3、写出详细的离散过程和求解方法的分布，并画出时，随的分布。 5、编程采用C/C++或FORTRAN语言。 6、将源程序附于作业之后，程序中要有详细的注释，以反映出思路 (1) 开始离散处理：视u、v为常数，对(1)式进行积分： (2) (3) 已知： (4) (5) (6) (7) 将式(4)、(5)、(6)、(7)分别代入式(2)、(3)，整理可得： (8) 因此可知： 若对流方向沿斜线向上，则u的方向为从左向右，v的方向为从下向上。可将以上各式进行化简，如下所示： 设沿斜线的对流速度为f，则可知： 其中对引入规正变量，进行延迟修正处理，且对于不同的格式有不同的表达式，具体如下： 规正变量： 因此： 对不同格式有不同表达式： FUD：； CD： SUD： QUICK： CLAM: EULER： MINMOD： MUSCL： OSHER： SMART： STOIC: 边界点的处理： 由于，因此，将上边界与右边界做绝热处理，认为边界处值等于内点值；在左上边界处，因此，将边界值可赋值为1，紧邻边界点的内点计算时按划分网格的一半计算系数参数来近似处理。 在下边界处，因此将边界值可赋值为0，紧邻边界点的内点计算时按划分网格的一半计算系数参数来近似处理。 特殊处理： 在计算边界处的时按一阶迎风格式近似处理来求解其值。 具体编程处理： 取网格数为100×100，所以x、y方向空间步长分别为0.01，进行编程计算。引入规正变量，采用延迟修正方法，采用SOR进行迭代（w=0.8），当累积误差err10-5时，迭代终止，输出各网格点数据并绘制图形。 二、程序 见附页。 三、y=0.5时的图形 见附页 四、计算结果分析 （1）计算过程中引入规正变量，采用延迟修正处理，保证编程过程中方程对角占优，从而易于求得稳定解。 （2）从图中可以得知：数值计算在剧烈变化区域（y=0.5处），采用CD、SUD、QUICK格式时，产生越界现象。 （3）计算过程中采用STOIC格式，计算效果最好。 （4）采用不同的格式时，其表达式在CBC线内，则会出现稳定解，否则会出现解得越界现象，越界现象与对流稳定性不同。 （5）编程过程中，采用SOR低松弛迭代方法，所得结果比较理想，计算速度远远Gauss—Seidel方法。在本次作业中，松弛系数取为0.8。 （6）编程过程中，要将边值点带入循环进行计算，否则会导致计算结果出错。进而也证明了，对流现象是具有方向性，其扰动只能沿下游传播。 五、心得体会 通过此次大作业，我们认识到大家合作的重要性，同时也锻炼了我们共同处理问题的能力。此次作业中遇到了很困惑的一个问题，就是一开始采用规正变量后得不到收敛解，原因是计算规正变量时分母会出现为0的情况，一开始的做法是当分母为0时在分母上加上一个很小的数，这样做的话或导致分数值很大，导致最终不收敛。最后我们的做法是，当分母为0时，所计算的边界值采用一阶迎风形式。从这次亲身体会，我们认识到在得到离散的格式之后，不要盲目的去算，而是应该判断其特点，然后再进行程序的编写，并且，在编程过程中要熟知语句的用法，这样不仅会节省很多时间，而且容易发现错误。同时，我们学会了解决问题的一些基本思想，并可以应用这些基本思想解决其它学科中所遇到的问题。 六、每个小组成员承担工作 工作 编程 绘图 报告、ppt 姓名 禹国军 王晏林 王岩 韩立涛 刘文婷 附页 (FUD)