方腔流动.pdfVIP

计算流体力学作业 题目 方腔流指顶部平板以恒定速度驱动规则区域内封闭的不可压流体（例如水）的流动，在方 腔流的流动中可以观察到几乎所有可能发生在不可压流体中的流动现象，如图 1 所示方腔流计 算模型图。 图1 方腔流动示意图 流函数-涡量法以流函数和涡量为未知量，可以消去控制方程中的压力项。根据差分法编写 程序计算方腔流动。 控制方程 边界条件 流函数边界条件：根据已知条件，在四个壁上流函数均为0，上边界平板速度1，根据不可 滑移条件确定流函数偏导数条件。 涡量边界条件，采用（1）： （1）可以用Taylor 展开建立一般形式的Thom 公式。假设某壁面切向速度 ，沿其内法 v  向n 有一节点，距离壁面距离为h, 此点上的流函数为 ，如果壁面上的流函数值为 ，那么  0 1 (  v h) 1 0     0 2h2 时间导数采用向前Euler 法，空间导数项可以用中心差分格式离散 计算步骤 1、计算n + 1 时刻内点的涡量，需要考虑时间步长和空间步长和粘性项的关系。 2 、计算n + 1 时刻的流函数n1 ，超松弛迭代法，松弛因子为-1.8。 3、计算n + 1 时刻边界上的涡量 初始条件如下： M=50+1; % y 方向网格点 N=50+1; % x 方向网格点 Nu=1e-2 ; % 粘性系数 V=1; % 上边界速度 sc=-1.8 ; %松弛因子 tol=1e-4; %迭代精度 dt=0.005; % 计算时间步长 t0=0; % 计算开始时间 t1=4; % 计算截止时间 tN=floor((t1-t0)/dt); % 计算时间步 x0=0; x1=1; dx=(x1-x0)/(N-1); %空间步长 y0=0; y1=1; dy=(y1-y0)/(M-1); %空间步长 Thom 公式 1 0.8 0.6 0.4 0.2 0 -0.2 0 0.2 0.4 0.6 0.8 1 1.2