计算流体力学上机实验报告

1、计算流体力学上机实验报告班级： 姓名： 学号：北京航空航天大学流体力学研究所上机实验名称两平行平板间不可压缩流体绕物体的平面无旋流动一、 实验目的通过具体算例，熟悉和掌握使用CFD方法获取给定流场的流动参数。二、 实验内容、方法及步骤1. 流动问题描述如下图所示，考虑在平行放置的两平板之间流过的理想不可压缩流体绕方形物体的平面无旋流动。2. 求解区域设两平行平板之间的距离为 ；绕流物体是边长为 的正方形，上游来流入口位置与物体中心的距离为 。根据流动的对称性，可取流动区域的四分之一作为求解区域 ，如下图所示。3. 控制方程对于不可压缩流体的平面无旋流动，流函数 在区域 内满足Laplace方程

2、4. 边界条件（1）OABC 是一条流线，规定 ；（2）对 OE 上任意一点 ，有 ；（3） ED 也是一条流线，所以 ；（4）根据对称性，在 CD 上有。5. 定解问题对于这里考虑的流动，可用下述定解问题来描述6. 求解区域的离散化 计算网格将单位长度等分成份，记，于是求解区域沿方向可划分成个网格，用来标记；沿方向可划分成个网格，用来标记。这些网格点可分成三类：（1）当 且 ，或者当 且 时，网格点落在流场 内部，称为内点。这些网格点上的流函数需通过求解方程组来计算；（2）当 且 时，网格点落在正方形物体内部，网格点上不存在流场，无需计算；（3）其余的网格点落在流场 的边界上，称为边界点。这

3、些网格点上的流函数直接由边界条件给定，也无需计算。7. 定解问题的离散化 差分格式Laplace方程 的差分近似为边界条件 的差分近似为8. 内点上流函数的计算 迭代算法在实际的计算中，内点的数量非常多，计算流函数需要求解大型的代数方程组。将Laplace方程的差分近似改写成使用高斯-赛德尔迭代算法计算： （ ） 给定迭代的初始近似 和迭代的收敛精度 。 对 ： 在流场内的每一个内点上用迭代公式进行计算。等式右边如果有边界点，可直接将边界条件指定的流函数值代入。 在边界 CD上的每一个网格点（不包括端点C、D）处，用公式 计算。 对所有内点，计算 。若 ，迭代收敛。否则将 加 ，返回步骤39.

4、 计算源程序（matlab 编制）function liuxian,vx,vy=jisuanliuti(x,y)%this is for Computational Fluid Dynamics;%使用方式：liuxian,vx,vy=jisuanliuti(x,y)；%输入x:x轴方向单位长度上划分的网格（x=3）%输入y:y方向上单位长度上划分的网格(y=3)%输出liuxian:输出各个网格点的流函数值%输出vx：输出x方向的速度u%输出vy：输出y方向的速度v%if (x3)|(y3) error(number of grids can not be divided too littl

5、e);end %参数定义%n=x;%x方向单位长度划分n个h=y;%y方向单位长度划分h个H=3;%外边界的流函数值fai=ones(3*n+1,3*h+1);%流函数%边界条件%fai(:,1)=0;fai(2*n+1:3*n+1,1:h+1)=0;fai(1,:)=(1:3*h+1)-1)*3/(3*h);fai(:,end)=H;fai(end,:)=fai(end-1,:);%计算流函数，用的是高斯-赛德尔迭代%fai1=fai;for i=1:1000 for j=2:1:3*n for k=2:1:3*h t=(j=2*n+1)&(k=h+2); if t=1 fai1(j,k)=

6、1/4*(fai1(j+1,k)+fai1(j-1,k)+fai1(j,k+1)+fai1(j,k-1); end end end c=max(max(abs(fai1-fai); if c1e-5 fai=fai1; else break; end fai1(end,:)=fai1(end-1,:);endliuxian=fai;%画等高线图%figure;contour(0:1/n:3,0:1/h:3,liuxian,8, LineWidth,2);title(流线图,fontsize,14);%改变流函数矩阵格式%for i=1:length(liuxian(:,1) liuxian1(

7、i,:)=liuxian(length(liuxian(:,1)+1-i,:);endliuxian=liuxian1;%速度的计算%x方向速度%vx=(liuxian(1:end-1,:)-liuxian(2:end,:)*h;%y方向速度%vy=-(liuxian(:,2:end)-liuxian(:,1:end-1)*n;end三、 实验结果 边界CD上的速度u的值如下表：U1.67651.6561.6151.5581.502Y11.1251.3751.6251.875U1.4561.4231.4011.391.3845Y2.1252.3752.6252.8753四、 讨论1. x,y方向划分不同的网格计算精度有所不同，划分网格越多，精度越高。下图是分别按（3,3）和（5,5）划分网格时的流线图。Figure 2网格划分为（3，3）Figure 1网格划分为（5，5）从两幅图处可以清晰地看到，网格