流体力学大作业
西工大-计算流体力学大作业
计算流体力学大作业学号： 姓名：1、不可压平面流通过二维容器（如图）。
采用 简单迭代、超松弛迭代 求解 势流方程获得容器内的速势和速度分布 。
边界条件按照课本中给，即流经 A 、B 的体积流量为1。
要求： 1）推导差分方程的迭代公式；2）编写计算机程序 ； 3）绘制计算结果曲线 。
答：1）迭代公式推导对于容器中的定常流场，其支配方程为22220x yφφ∂∂+=∂∂ 求解域为下图所示矩形区域则支配方程由有限差分形式代换，得1,,1,,1,,122220()()i j i j i ji j i j i j x y φφφφφφ+-+--+-++=∆∆具有22()()x y ∆+∆的截断误差对于正方形网格，有22()()x y h ∆=∆=，则上式可改写为n=17,1,1,,1,11()4i j i j i j i j i j φφφφφ+-+-=+++若采用简单迭代公式，即Liebmann 公式，则有(1)()(1)()(1),1,1,,1,11()4n n n n n i j i j i j i j i j φφφφφ++++-+-=+++若采用超松弛迭代，即SOR 公式，则有(1)()()(1)()(1),,1,1,,1,1(1)()4n n n n n n i j i j i j i j i j i j ωφωφφφφφ++++-+-=-++++其中松弛因子12ω<<。
ω最佳值opt ω为opt ω=式中cos(/)cos(/)m n αππ=+，m ，n 分别表示在网格系统中垂直线和水平线的总数。
2）计算机程序本程序采用C 语言编写。
程序源代码如下： #include<stdio.h> #include<math.h> void main() { int m=25,n=17,ilast[17],jlast[25]; int step1,step2; double h=0.25; double psi_j[25][17],psiprv_j,vel_j[25][17],velx_j[25][17],vely_j[25][17]; double psi_c[25][17],psiprv_c,vel_c[25][17],velx_c[25][17],vely_c[25][17]; double Pi,Alpha,Omega,Error; int i,j; for(i=0;i<17;i++) jlast[i]=17; for(i=17;i<m;i++) jlast[i]=17-(i-16); for(j=0;j<9;j++) ilast[j]=25; for(j=9;j<n;j++) ilast[j]=25-(j-8); //数据初始化 for(j=0;j<n;j++) { psi_j[0][j]=1.0; psi_c[0][j]=1.0;}for(i=1;i<m;i++){psi_j[i][jlast[i]-1]=1.0;psi_c[i][jlast[i]-1]=1.0; }for(j=0;j<8;j++){psi_j[m-1][j]=1.0;psi_c[m-1][j]=1.0;}for(i=1;i<m-1;i++){if(i>6 && i<21){psi_j[i][0]=0.0;psi_c[i][0]=0.0;}else{psi_j[i][0]=1.0;psi_c[i][0]=1.0;}}for(i=1;i<m-1;i++){for(j=1;j<jlast[i]-1;j++){psi_j[i][j]=0.5;psi_c[i][j]=0.5;}}//处理右上角数据for(i=0;i<m;i++){for(j=0;j<n;j++){if(j>jlast[i]-1){psi_j[i][j]=0;vel_j[i][j]=3;psi_c[i][j]=0;vel_c[i][j]=3;}}}Pi=4.0*atan(1.0);Alpha=cos(Pi/m)+cos(Pi/n);Omega=(8.0-4*sqrt(4-pow(Alpha,2)))/pow(Alpha,2);//计算速势step1=0;step2=0;//简单迭代while(1){Error=0.0;for(i=1;i<m-1;i++){for(j=1;j<jlast[i]-1;j++){psiprv_j=psi_j[i][j];psi_j[i][j]=(psi_j[i-1][j]+psi_j[i+1][j]+psi_j[i][j-1]+psi_j[i][j+1])/4.0;Error=Error+fabs(psi_j[i][j]-psiprv_j);}}step1++;if(step1>1000)break;if(Error<=0.001)break;}//超松弛迭代while(1){Error=0.0;for(i=1;i<m-1;i++){for(j=1;j<jlast[i]-1;j++){psiprv_c=psi_c[i][j];psi_c[i][j]=(1-Omega)*psi_c[i][j]+Omega*(psi_c[i-1][j]+psi_c[i+1][j]+psi_c[i][j-1]+psi_c[i][j+1])/4.0;Error=Error+fabs(psi_c[i][j]-psiprv_c);}}step2++;if(step2>1000)break;if(Error<=0.001)break;}//计算速度for(i=0;i<m;i++){for(j=0;j<jlast[i];j++){if(j==0){vely_j[i][j]=(-3*psi_j[i][j]+4*psi_j[i][j+1]-psi_j[i][j+2])/2/h;vely_c[i][j]=(-3*psi_c[i][j]+4*psi_c[i][j+1]-psi_c[i][j+2])/2/h;}else if(j==jlast[i]-1){vely_j[i][j]=(psi_j[i][j-2]-4*psi_j[i][j-1]+3*psi_j[i][j])/2/h;vely_c[i][j]=(psi_c[i][j-2]-4*psi_c[i][j-1]+3*psi_c[i][j])/2/h;}else{vely_j[i][j]=(psi_j[i][j+1]-psi_j[i][j-1])/2/h;vely_c[i][j]=(psi_c[i][j+1]-psi_c[i][j-1])/2/h;}}}for(j=0;j<n;j++){for(i=0;i<ilast[j];i++){if(i==0){velx_j[i][j]=(-3*psi_j[i][j]+4*psi_j[i+1][j]-psi_j[i+2][j])/2/h;velx_c[i][j]=(-3*psi_c[i][j]+4*psi_c[i+1][j]-psi_c[i+2][j])/2/h;}else if(i==ilast[j]-1){velx_j[i][j]=(psi_j[i-2][j]-4*psi_j[i-1][j]+3*psi_j[i][j])/2/h;velx_c[i][j]=(psi_c[i-2][j]-4*psi_c[i-1][j]+3*psi_c[i][j])/2/h;}else{velx_j[i][j]=(psi_j[i+1][j]-psi_j[i-1][j])/2/h;velx_c[i][j]=(psi_c[i+1][j]-psi_c[i-1][j])/2/h;}}}for(i=0;i<m;i++){for(j=0;j<jlast[i];j++){vel_j[i][j]=sqrt(pow(velx_j[i][j],2)+pow(vely_j[i][j],2));vel_c[i][j]=sqrt(pow(velx_c[i][j],2)+pow(vely_c[i][j],2));}}//输出结果分布FILE *fp;fp=fopen("f:\\ESL\\YFresult.txt","w");fprintf(fp,"简单迭代结果\n");fprintf(fp,"速度势分布\n");for(j=n-1;j>=0;j--){for(i=0;i<ilast[j];i++){fprintf(fp,"%-10.6f\n",psi_j[i][j]);}}fprintf(fp,"速度分布\n");for(j=n-1;j>=0;j--){for(i=0;i<ilast[j];i++){fprintf(fp,"%-10.6f\n",vel_j[i][j]);}}fprintf(fp,"超松弛迭代结果\n");fprintf(fp,"速度势分布\n");for(j=n-1;j>=0;j--){for(i=0;i<ilast[j];i++){fprintf(fp,"%-10.6f\n",psi_c[i][j]);}}fprintf(fp,"速度分布\n");for(j=n-1;j>=0;j--){for(i=0;i<ilast[j];i++){fprintf(fp,"%-10.6f\n",vel_c[i][j]);}}fclose(fp);//输出tecplot数据FILE *fp1;fp1=fopen("f:\\ESL\\TECPLOT-result.txt","w");fprintf(fp1,"title=erwei grid\n");fprintf(fp1,"variables=x, y, psi_easy, velocity_easy, psi_SOR\n, velocity_SOR\n");fprintf(fp1,"zone t=grid,i=25,j=17,f=point\n");for(j=0;j<n;j++){for(i=0;i<m;i++){fprintf(fp1,"%-10.6f,%-10.6f,%-10.6f,%-10.6f,%-10.6f,%-10.6f\n",i*h,j*h,psi_j[i][j],vel_j[i][j],p si_c[i][j],vel_c[i][j]);}}fclose(fp1);}3）计算结果采用简单迭代，容器内的速势和速度分布速势分布（简单迭代）速度分布（简单迭代）采用超松弛迭代，容器内的速势和速度分布速势分布（SOR ） 速度分布（SOR ）2、用点源（汇）分布在对称轴的源汇模拟流体绕过NACA0012旋称体的二维轴对称势流解。
(完整版)流体力学作业试题库及答案
第一章 绪论思考题1-1 何谓流体连续介质模型？含有气泡的液体是否适用连续介质模型？ 答：所谓流体的连续介质模型，即把流体视为没有间隙地由流体质点充满它所占据的整个空间的一种连续介质其物理性质和物理量也是连续的。
若气泡相对于液体而言可以看作孤立的点的话，则含有气泡的液体可以适用连续介质模型。
习题11-3 如题图所示，设平行板间隙为0.5mm ，中间充满液体，上板以U ＝0.25m/s 的速度平移，施于单位面积的力为2Pa ，试求液体的粘度为多少？解：YU dy du A F μμτ===液体粘度s Pa AU FY ⋅⨯=⨯⨯==--3310425.0105.02μ1-4 求题图所示的轴与轴套之间的流体粘度。
解：s Pa dLU FY dLA Y U dy du A F ⋅=⨯⨯⨯⨯⨯⨯==⇒====--0648.0493.010)140120(14.3102.034.863πμπμμτ第二章 流体静力学习题22-5 用多管水银测压计测压，，题图中标高的单位为m ，试求水面的压强p 0。
解：Pam g m g p pap m m g p p m m p p m m g p p m m g p p D D CC B B A A 5001065.29.298002.21334169.22.20)2.13.2()2.15.2(g )4.15.2()4.10.3(⨯=⨯-⨯=⨯-⨯=⇒⎪⎪⎪⎩⎪⎪⎪⎨⎧=-+=--=-+=-+=水汞汞水汞水ρρρρρρ2-9 一盛水的敞口容器作加速运动，试求下列两种情况下容器内静压强的分布规律：（1）自由降落；（2）以等加速度a 向上运动。
解：h a g p p )sin (0αρ++=（1）0,900=∴=︒-=p p 相对压强α（2）)(,900a g h p p p p a a ++=∴=︒=ρα绝对压强2-12 试求开启题图所示水闸闸门所需的单宽拉力F 。
不计闸门自重及转轴摩擦力。
流体力学大作业实验报告
END DO PRINT*,"摩擦系数为:",S END PROGRAM
运行结果：
∴ = 0.033379640， ∴ V = 58.8 ÷ (1.8 + 70 ) /s = 3.77m/s 经检验， = = 433001 > 4000，为紊流，符合假设。 ∴流量Q = AV = × V = 0.0666 /s
S=(1/(-2.0*LOG10((A/D)/3.7+2.51/(SQRT(58.8/(1.8+70*T))*D/U*SQRT(T)))))**2
DO WHILE(ABS(S-T)>E) PRINT*,"摩擦系数为:",S T=S S=(1/(-2.0*LOG10((A/D)/3.7+2.51/(SQRT(58.8/(1.8+70*T))*D/U*SQRT(T)))))**2
∴ V = 58.8 ÷ (1.8 + 70 ) /s
由科尔布鲁克公式 = −2.0 ( . + . )，其中 =
查表得， = 10℃时，运动粘性系数 = 1.306 × 10 /s D=150mm, Δ = 1 将V的表达式代入科尔布鲁克公式得：
= −2.0 ( + . )=−2.0 ( +
.
.
.
==6−.31606mm//ss，
流体力学大作业
7-36
题目：设有一虹吸管，其装置如图所示，设管径为 150mm。
试求通过该管的流量。已知水温为 10℃，∆=
，
弯头 = . ,， 进口 = . 。
推导过程：
选取截面 1-1,2-2 如右图所示， = 0 , = −3 , = = 由总流伯努利方程
流体力学大作业
“水流动力学基本原理的应用”大作业
姓名：
学号：
专业班级：
成绩：
教师评语：
年月日
1、对水流流向问题有如下一些说法：“水一定从高处向低处流”,“水一定从压强大的地方向压强小的地方流”，“水一定从流速大的地方向流速小的地方流”。
这些说法是否正确？为什么？
2、在写总流能量方程221112221222p p z z h g g g g
ωαυαυρρ++=+++ 时，过水断面上的计
算点、基准面、压强标准是否可以任意选取？为什么？
3、液流通过如图所示管道流入大气中，已知：U 形测压管中水银柱高差m h Hg 2.0=∆，
10.72h m =水柱高，管径m d 1.01=，管嘴出口直径m d 05.02=，不计管中水头损失，试
求：管中流量Q 。
4、如图所示为一水平面上的渐变弯管，已知：
断面1-1处的压强2
3
1/1098m N p ⨯=，流速14/m s υ=，管径mm d 2001=，管径mm d 1002=，转角
45=α，不计弯管的水头损失。
试求：水流作用在弯管上的力
5、如图所示为闸下底板上的消力墩，已知：跃前断面水深h 1=0.6m ，流速v 1=15m/s ，跃后断面水深h 2=4m ，墩宽b=1.6m ，试求：水流对消力墩的作用力。
题2.9图
6、学过恒定总流能量方程及动量方程及其应用这部分内容以后，你觉得有些什么收获？有什么疑惑或者模糊的地方？。
计算流体力学课程大作业
《计算流体力学》课程大作业——基于涡量-流函数法的不可压缩方腔驱动流问题数值模拟张伊哲 航博1011、 引言和综述2、 问题的提出，怎样使用涡量-流函数方法建立差分格式3、 程序说明4、 计算结果和讨论5、 结论1引言虽然不可压缩流动的控制方程从形式上看更为简单，但实际上，目前不可压缩流动的数值方法远远不如可压缩流动的数值方法成熟。
考虑不可压缩流动的N-S 方程：01()P t νρ∇⋅=⎧⎪∂⎨+∇⋅=-∇+∆⎪∂⎩U UUU f U (1.1)其中ν是运动粘性系数，认为是常数。
将方程组写成无量纲的形式：01()Re P t∇⋅=⎧⎪∂⎨+∇⋅=-∇+∆⎪∂⎩U UUU f U (1.2) 其中Re 是雷诺数。
从数学角度看，不可压缩流动的控制方程中不含有密度对时间的偏导数项，方程表现出椭圆-抛物组合型的特点；从物理意义上看，在不可压缩流动中，压力这一物理量的波动具有无穷大的传播速度，它瞬间传遍全场，以使不可压缩条件在任何时间、任何位置满足，这就是椭圆型方程的物理意义。
这就造成不可压缩的N-S 方程不能使用比较成熟的发展型．．．偏微分方程的数值求解理论和方法。
如果将动量方程和连续性方程完全耦合求解，即使使用显示的离散格式，也将会得到一个刚性很强的、庞大的稀疏线性方程组，计算量巨大，更重要的问题是不易收敛。
因此，实际应用中，通常都必须将连续方程和动量方程在一定程度上解耦。
目前，求解不可压缩流动的方法主要有涡量-流函数法，SIMPLE 法及其衍生的改进方法，有限元法，谱方法等，这些方法各有优缺点。
其中涡量-流函数法是解决二维不可压缩流动的有效方法。
作者本学期学习了研究生计算流体课程，为了熟悉计算流体的基本方法，选择使用涡量-流函数法计算不可压缩方腔驱动流问题，并且对于不同雷诺数下的解进行比较和分析，得出一些结论。
本文接下来的内容安排为：第2节提出不可压缩方腔驱动流问题，并分析该问题怎样使用涡量-流函数方法建立差分格式、选择边界条件。
流体力学大作业
《计算流体力学》课程大作业作业内容：3-4人为小组完成数值模拟，在第8次课上每组进行成果展示，并在课程结束后每组上交一份纸质版报告。
数值模拟实现形式：自编程或者使用任意的开源、商业模型。
成果展示要求：口头讲述和幻灯片结合的方式，每组限时10分钟（8分钟讲述，2分钟提问和讨论）。
报告要求：按照期刊论文的思路和格式进行撰写（包括但不限于如下内容：摘要、绪论\引言、数值模型简介、数值结果分析\讨论、结论、参考文献）。
（以下题目二选一）题目一：固定单方柱扰流问题根据文章《Interactions of tandem square cylinders at low Reynolds numbers》中的实验进行数值模拟，完成但不局限于如下工作：（1）根据Fig. 2 中的雷诺数和方柱排列形式，进行相同雷诺数不同间距比情况下的方柱绕流数值模拟，并做出流线图和Fig.2中的结果对比。
（2）根据Fig. 3 中的雷诺数和方柱排列形式，进行相同雷诺数后柱不同转角情况下的方柱绕流数值模拟，并做出流线图和Fig.3中的结果对比。
（3）根据Fig. 12, 13 中的雷诺数和方柱间距比的设置进行数值模拟，作出频率、斯特劳哈尔数、阻力系数随雷诺数变化的折线并与图中对应的折线画在同一坐标系下比较。
（中共有4条折线，对应4种不同的方柱排列形式下的物理参数随雷诺数变化的规律，仅需选取单柱模型和其中一种双柱模型进行数值模拟，共计16个工况）。
题目二：溃坝问题根据文章《Experimental investigation of dynamic pressure loads during dam break》中的实验进行数值模拟，完成但不局限于如下工作：（1）分别完成二维、三维的溃坝的数值建模，讨论二维、三维模型的区别。
（2）分别将二维、三维溃坝的数值模拟结果和Fig. 7,10中各时刻的自由面形态进行对比，并分别观测溃坝前端水舌的位置随时间的变化，其结果和Fig. 12 种的各试验结果放在同一坐标系下进行对比。
理想流体力学大作业(Hess-Smith方法求附加质量)
理想流体大作业应用hess—smith方法求附加质量一物理背景当船舶在海洋中运动或潜艇在深海中运动时，可以看成是物体在无界流中的运动，其所受到的力和力矩可以用物体的速度，加速度和附加质量决定，因此，附加质量是物体最重要的水动力特性。
对于附加质量的求解就显得非常重要。
由于附加惯性力的作用，物体在理想流体中的变速运动相当于物体自身质量上增加了一个附加质量而在真空中运动，换句话说，理想流体增大了物体的惯性，是物体很难加速也难减速。
附加质量有物体上的单位运动绝对速度势来决定，与流体密度，物体形状和运动方向有关，而与时间无关。
二理论依据物体的附加质量可由下时表示（2.1）(2.2)其中，分别是总速度势的分量，表示物法线方向在各个方向分量，表示物体在i方向运动引起其在j方向的附加质量,，那么，只要求出在物体表面速度势在法线方向的各个分量与法向量的分量，物体的附加质量便可求出。
在无界流中，物体的总速度势可以分解为定常速度是和扰动速度势。
其在物面上满足定解条件(2.3)应用分布源来模拟物面上的扰动势，则有(2.4)所以，在物体表面上有速度沿法线方向的分量表达式(2.5)由物面条件得到(2.6)所以只要求出分布源密度，就可以求出物体表面的速度势。
为了计算附加质量，要求单位速度势在物面上的值,其满足式（2.6）。
从而得到分布源的单位速度势密度满足的积分方程。
只要求解出单位速度势密度，便可以求得附加质量。
三数值模型主要通过在物面上的分布源来模拟物体在无界流中所产生的扰动势，应用数值离散的方法列代数方程组求解积分方程，即用离散量来表示连续变量的方法。
首先，将物体表面离散，划分成由n个四边形单元近似组成。
物面上各个单元的法向可由节点求出。
其次，可以列出求解的线性代数方程组其中(j/=i),从而，速度势密度可解，单位速度势可求，附加质量可以计算得到。
四几何模型对球体，椭球和圆柱进行网格的划分，应用不同数量的四边形网格来近似物体的表面，通过编写程序对其进行网格的划分，并按顺时针输出个节点的坐标值，相应的程序如图五计算参数主要对半径r=10的球壳，短轴长度为10，长短轴比为5的椭球，直径为20，柱体长和截面直径比为5的圆柱进行计算，对其进行网格的划分并求其附加质量系数。
计算流体力学大作业
1 提出问题[问题描述]Sod 激波管问题是典型的一类Riemann 问题。
如图所示，一管道左侧为高温高压气体，右侧为低温低压气体，中间用薄膜隔开。
t=0 时刻，突然撤去薄膜，试分析其他的运动。
Sod 模型问题：在一维激波管的左侧初始分布为：0 ,1 ,1111===u p ρ，右侧分布为：0 ,1.0 ,125.0222===u p ρ，两种状态之间有一隔膜位于5.0=x 处。
隔膜突然去掉，试给出在14.0=t 时刻Euler 方程的准确解，并给出在区间10≤≤x 这一时刻u p , ,ρ的分布图。
2 一维Euler 方程组分析可知，一维激波管流体流动符合一维Euler 方程，具体方程如下： 矢量方程：0U ft x∂∂+=∂∂ (0.1)分量方程：连续性方程、动量方程和能量方程分别是：222,,p u ρ()()()()2000u tx u u pt x x u E p E tx ρρρρ∂⎧∂+=⎪∂∂⎪⎪∂∂∂⎪++=⎨∂∂∂⎪⎪∂+⎡⎤∂⎣⎦+=⎪∂∂⎪⎩ (0.2)其中 22v u E c T ρ⎛⎫=+ ⎪⎝⎭对于完全气体，在量纲为一的形式下，状态方程为：()2p T Ma ργ∞=(0.3)在量纲为一的定义下，定容热容v c 为：()211v c Ma γγ∞=- (0.4)联立(1.2)，(1.3)，(1.4)消去温度T 和定容比热v c ，得到气体压力公式为：()2112p E u γρ⎛⎫=-- ⎪⎝⎭(0.5)上式中γ为气体常数，对于理想气体4.1=γ。
3 Euler 方程组的离散3.1 Jacibian 矩阵特征值的分裂Jacibian 矩阵A 的三个特征值分别是123;;u u c u c λλλ==+=-，依据如下算法将其分裂成正负特征值：()12222k k k λλελ±±+=(0.6)3.2 流通矢量的分裂这里对流通矢量的分裂选用Steger-Warming 分裂法，分裂后的流通矢量为()()()()()()()12312322232121212122f u u c u c u u c u c w γλλλργλλλγλλγλ⎛⎫⎪-++ ⎪=-+-++ ⎪ ⎪ ⎪-+-+++ ⎪⎝⎭＋＋＋＋＋＋＋＋＋＋＋(0.7)()()()()()()()12312322232121212122f u u c u c u u c u c w γλλλργλλλγλλγλ⎛⎫⎪-++ ⎪=-+-++ ⎪ ⎪ ⎪-+-+++ ⎪⎝⎭－－－－－－－－－－－(0.8)其中：()()()223321c w γλλγ±±±-+=- c 为量纲为一的声速：22Tc Ma ∞=(0.9)联立(1.3)，(1.9)式，消去来流马赫数得：ργp c =3.3 一阶迎风显示格式离散Euler 方程组 10n n i i x i x i U U f f t xδδ+-++--++=∆∆ (0.10)得到()()n+1nj j 11U =U j j j j t f f f f x++---+∆⎡⎤--+-⎣⎦∆ 算法如下：① 已知初始时刻t=0的速度、压力及密度分布000,,j j j u P ρ，则可得到特征值分裂值0k λ±，从而求出流通矢量0j f ±；② 应用一阶迎风显示格式可以计算出1t t =∆时刻的组合变量1j U ，从而得到1t t =∆时刻的速度、压力及密度分布111,,j j j u P ρ；③ 利用1t t =∆时刻的速度、压力及密度分布111,,j j j u P ρ可得特征值分裂值1k λ±，从而求出流通矢量1j f ±；④ 按照步骤2的方法即可得到2t t =∆时刻的速度、压力及密度分布222,,j j j u P ρ；⑤ 循环以上过程即可得到()1t n t =+∆时刻的速度、压力及密度分布n+1n+1n+1,,j j j u P ρ。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
《计算流体力学》课程大作业
作业内容：3-4人为小组完成数值模拟，在第8次课上每组进行成果展示，并在课程结束后每组上交一份纸质版报告。
数值模拟实现形式：自编程或者使用任意的开源、商业模型。
成果展示要求：口头讲述和幻灯片结合的方式，每组限时10分钟（8分钟讲述，2分钟提问和讨论）。
报告要求：按照期刊论文的思路和格式进行撰写（包括但不限于如下内容：摘要、绪论\引言、数值模型简介、数值结果分析\讨论、结论、参考文献）。
（以下题目二选一）
题目一：固定单方柱扰流问题
根据文章《Interactions of tandem square cylinders at low Reynolds numbers》中的实验进行数值模拟，完成但不局限于如下工作：
（1）根据Fig. 2 中的雷诺数和方柱排列形式，进行相同雷诺数不同间距比情况下的方柱绕流数值模拟，并做出流线图和Fig.2中的结果对比。
（2）根据Fig. 3 中的雷诺数和方柱排列形式，进行相同雷诺数后柱不同转角情况下的方柱绕流数值模拟，并做出流线图和Fig.3中的结果对比。
（3）根据Fig. 12, 13 中的雷诺数和方柱间距比的设置进行数值模拟，作出频率、斯特劳哈尔数、阻力系数随雷诺数变化的折线并与图中对应的折线画在同一坐标系下比较。
（中共有4条折线，对应4种不同的方柱排列形式下的物理参数随雷诺数变化的规律，仅需选取单柱模型和其中一种双柱模型进行数值模拟，共计16个工况）。
题目二：溃坝问题
根据文章《Experimental investigation of dynamic pressure loads during dam break》中的实验进行数值模拟，完成但不局限于如下工作：
（1）分别完成二维、三维的溃坝的数值建模，讨论二维、三维模型的区别。
（2）分别将二维、三维溃坝的数值模拟结果和Fig. 7,10中各时刻的自由面形态进行对比，并分别观测溃坝前端水舌的位置随时间的变化，其结果和Fig. 12 种的各试验结果放在同一坐标系下进行对比。
（3）根据实验设置数值观测点，分别观测与实验测点相对应的数值观测点上的水体高度、压力随时间的变化曲线，并和Fig.16, 18,21,30,31,32,33,35中的实验结果进行对比。