计算流体实验报告4.docVIP

中山大学 本科生实验报告书 院 系 工学院应用力学与工程系 专业班级 理论与应用力学07级 实验课程 计算流体力学 姓 名 丁 鹏 学 号 指导教师 詹杰民 余凌晖 实验四 方腔环流问题 实验目的 用数值方法计算二维不可压缩无粘流体方腔环流的流函数和势函数，用相关软件绘制流函数和势函数的图形并作分析。 实验内容 如下图所示，是二维方形腔体，腔体内部充满流体。当顶板眼水平方向被均匀拉动时，腔体内的流体将被带动而做环状运动。这种环流导致了腔体底边出现二次涡。 方程和边界条件 用流函数涡量法，，满足下列无量纲形式的定常方程和边界条件 ，在腔体四边 ，在AB和CD上 ，在AD和BC上 其中，Re为雷诺数 图表 1 方腔环流 2.取正方形网格（如图6.7所示）， 图表 2 方腔环流网格 用二阶精度的差商代替上式中的微商，得 引进松弛因子，方程可化为下列的迭代格式 为了得到涡量的边界条件，令方程在腔体四边也成立，利用Taylor公式，边界条件，边界条件成立 ，在四边 ，在两侧和底边 ，在顶面 3.计算步骤 求解过程，可按照下列步骤进行， 给定， 的初值。可取 令。 由上述方程解出内点的新值和。 由方程解出边界点的新值，由于在边界上，因此，不必计算。 n+1=n 重复步骤3），4）和5）直至 其中，，，C为给定精度，可取 4.计算结果 随着雷诺数的增加，减少，对于一个的网格，当时，可取，而当时，可取，此时计算结果如图所示 结果与分析 一、流函数的数据与作图分析 Re=10，25*25的网格，迭代精度，流线图和涡量等值线图如下（同时利用surfer和fluent画图） 图表 3 RE=10流函数图形 图表 4 RE=10涡量线 从这个流函数图中可以清楚看到一次涡，可以观察底边两个角处有二次涡不是很明显，整个流函数图大致对称。而且利用surfer和fluent进行数据对比。 Re=100，50*50的网格，迭代精度，流线图和涡量等值线图如下（同时利用surfer和fluent画图） 图表 5 RE=100方腔环流的流函数 图表 6 RE=100涡量线 从图表 RE=100方腔环流的流函数，可以清楚观察底边两个角处有明显的二次涡，整个流函数图不再对称。 Re=500，100*100的网格，迭代精度，流线图和涡量等值线图如下（只利用surfer） 图表 7 RE=500方腔环流的流函数 此图可以看到二次涡，及其位置，对于流函数涡量法计算流函数时，当RE增大的时候对应的二次涡也会变大。 Re=1000，50*50的网格，迭代精度，流线图和涡量等值线图如下（同时利用surfer和fluent画图） 图表 8 RE=1000方腔环流的流函数 图表 9RE=1000涡量线 Re=2000，150*150的网格，迭代精度，流线图和涡量等值线图如下（同时利用surfer和fluent画图） 图表 10 RE=2000方腔环流的流函数 图表 11 RE=2000涡量线 Re=3000，150*150的网格，迭代精度，流线图和涡量等值线图如下（同时利用surfer和fluent画图） 图表 12 RE=3000方腔环流的流函数 图表 13 RE=3000涡量线 当用到该种方法计算Re=4000的时候，就会计算的结果又很大偏差，所产生的低端的两个二次涡就会在一起，一次涡整个就会往上提。 Re=4000和5000，利用fluent进行求解，定常，线性粘性（牛顿粘性），求解其等值流函数和等涡线（如下图所示） 图表 16 RE=4000方腔环流 图表 17 RE=4000涡量线 图表 14 RE=5000方腔环流流函数 图表 15 RE=5000涡量线 实验总结 此次实验根据流函数涡量法，计算二维不可压粘性方腔环流问题。在流函数和 1）在自己的操作过程中当自己的雷诺数超过RE=4000的时候，所画的图形就不稳定了。 2）此次在计算过程中，迭代过程中所使用的是流函数涡量法。而在定常的流函数方程和涡量方程都采用的是超松弛迭代法，选用的松弛因子很重要。在程序当中两个方程分别采用不同的松弛因子，这对于迭代的次数会大大降低。 3)本次实验处理数据多，充分利用surf，fluent等软件，在处理不同雷诺数下的流线情况，在自己的程序当中直接设雷诺数re的大小，而且程序中不得超过最大的临界雷诺数，该雷诺数与本身的网格划分有关，还与自己所用的定常方程及超松弛迭代法有关。在FLUENT，利用的是SIMPLE方法，而且在网格划分的很细（100*100）（如图所示），可以算到Re=5000以上的