文档介绍：莁中山大学罿螈本科生实验报告书羇肂羂螈肃院系工学院应用力学与工程系袄专业班级理论与应用力学07级螀实验课程计算流体力学袈姓名丁鹏蒄学号07300129节指导教师詹杰民余凌晖蕿羈袅羄莈肈芆蒂莁***蒃实验四方腔环流问题膄实验目地膀用数值方法计算二维不可压缩无粘流体方腔环流地流函数和势函数,,是二维方形腔体,,,,满足下列无量纲形式地定常方程和边界条件芇芅莄,在腔体四边羂,在AB和CD上莇,在AD和BC上蚆其中,(),肇蒃图表2方腔环流网格葿用二阶精度地差商代替上式中地微商,得薇膃引进松弛因子,方程可化为下列地迭代格式羁芈蚇为了得到涡量地边界条件,令方程在腔体四边也成立,利用Taylor公式,边界条件,边界条件成立薄,在四边蚃,在两侧和底边芁,,可按照下列步骤进行,膁给定,,由于在边界上,因此,+1=n蝿重复步骤3),4)和5)直至袆螇其中,,,C为给定精度,,减少,对于一个地网格,当时,可取,而当时,可取,此时计算结果如图所示羄结果与分析肃一、流函数地数据与作图分析薁Re=10,25*25地网格,迭代精度,流线图和涡量等值线图如下(同时利用surfer和fluent画图)DXDiTa9E3d肆莅螄图表3RE=10流函数图形莀膆图表4RE=10涡量线螅从这个流函数图中可以清楚看到一次涡,可以观察底边两个角处有二次涡不是很明显,=100,50*50地网格,迭代精度,流线图和涡量等值线图如下(同时利用surfer和fluent画图)5PCzVD7HxA膈芆膆图表5RE=100方腔环流地流函数袄膁图表6RE=100涡量线莆从图表RE=100方腔环流地流函数,可以清楚观察底边两个角处有明显地二次涡,=500,100*100地网格,迭代精度,流线图和涡量等值线图如下(只利用surfer)莂羀蒅图表7RE=500方腔环流地流函数蚄此图可以看到二次涡,及其位置,对于流函数涡量法计算流函数时,=1000,50*50地网格,迭代精度,流线图和涡量等值线图如下(同时利用surfer和fluent画图)jLBHrnAILg蝿蝿肅图表8RE=1000方腔环流地流函数薂螂图表9RE=1000涡量线衿Re=2000,150*150地网格,迭代精度,流线图和涡量等值线图如下(同时利用surfer和fluent画图)xHAQX74J0X蒆