二维方腔流的涡量流函数法数值模拟.docxVIP

粘性不可压缩二维方腔剪切流的涡量流函数法数值模拟 令狐烈 浙江大学 摘要：采用了涡量-流函数法对粘性不可压缩流体的二维瞬态流动进行模拟计算，并使用Tecplot 360软件对计算结果进行可视化处理分析，研究了不同雷若数下方腔流流场结构的变化。 关键词：二维方腔剪切流，涡量流函数方法，数值模拟。 数学模型 1.1问题描述 粘性不可压缩二维方腔流问题表述如下： 有一宽度为 L 的方腔，充满粘性不可压缩流体，初始时刻流体静止不动，如图1所示。下壁面以常速值速度v运动，而其余三壁面固定不动。由于是粘性流体，将会使整个方腔内的流体运动起来，一定时间后达到稳定状态。将方腔划分为 100×100 的网格。 图一 粘性不可压缩二维方腔流问题示意图（左图）与网格示意图（右图） 1.2.无量纲化的控制方程 该问题的控制方程为二维不可压缩N-S方程，使用速度V和长度L进行无量纲化得到的涡量流函数方程为： 其中为雷诺数，涡量定义为。 对涡量速运方程（2）使用FTCS差分得到离散方程： 即： 对方程（1）用中心差分得到离散方程 即： 用显式格式求解： 其中Wp为松弛因子 1.3网格划分 对于本问题，采用等间距网格，将方腔 区域划分为n*n的网格，则有 故： 1.4边界条件 流动采用无滑移边界条件，壁面处速度为零。 流函数边界条件： 其中n表示网格划分份数。 涡量边界条件： 结果分析