西安交大流体力学上机实验
西安交通大学流体力学边界层测定实验报告2
实验数据处理及计算：
150mm截面数据
250mm截面数据
数据结果分析：
1.由雷诺数判断流态：
临界雷诺数Re=3∗10
流态判断：150mm截面为层流流动
250mm截面为层流向紊流过度区域2.根据边界层的速度分布判断流态：
由绘制的曲线分析，实测曲线均与紊流理论曲线吻合较好。
判断结果为：150mm截面、250mm截面均为紊流流态
3.根据边界层厚度判断流态：
层流：δ=5∗√υx
V
紊流：δ=0.37*υ1
5*x
由以上数据判断结果为：150mm截面流态为：紊流250mm截面流态为：紊流
实验总结与思考：
通过如上三种方法认为，通过实验，该平板模型在实验流场中，150mm截面处与250mm界面处均为紊流流态。
原因可能是风洞中流速过快，以及1截面选择过于靠后，因而测不到或测得层流流态。
建议下次试验对100mm截面进行测试。
思考题：
1.流体的流动状态受到哪些因素的影响？
答：Re=LVρ/μ,影响因素有特征长度L，,流速u,流体密度ρ,流体粘度μ.而温度会影响流体的粘度和密度。
2.为何层流和紊流呈现不同的速度分布规律？
答：两种状态和流动的雷诺数雷诺数小,意味着流体流动时各质点间的粘性力占主要地位,流体各质点平行于管路内壁有规则地流动,呈层流流动状态.雷诺数大,意味着惯性力占主要地位,流体呈紊流流动状态。
流体力学上机实验报告
流体力学上机实验报告一、 实验问题描述该算例是一个二维的导流管分析，分析一个雷诺数为400的层流情况，所用二、 分析方法及假定用FLUID141单元来作二维分析，本算例作了如下分析：雷诺数为400的假想流的层流分析，分析时假定进口速度均匀，并且垂直于进口流场方向上的流体速度为零。
在所有壁面上施加无滑移边界条件（即所有速度分量都为零）；假定流体不可压缩，并且其性质为恒值，在这种情况下，压力就可只考虑相对值，因此在出口处施加的压力边界条件是相对压力为零。
分析时，流场为层流，着可以通过雷诺数来判定，其公式为：R e =ρvD ℎμ三、 几何尺寸及流体性质进口段长度 4 m 进口段高度 2 m 过渡段长度 2m 出口段高度 2 m 出口段长度 4 m 假设流体密度 1 Kg/m3 假设流体粘性 0.01Kg/m -s ； 空气密度 1.205 Kg/m3 空气粘性 1.8135*10-5 Kg/m -s 进口速度 2.0 m/s 出口压力 0 nt/m2四、 分析过程1. 进入ANSYS参见ANSYS Operation Guide进口 出口过渡段 流体流动方向外壁面内壁面2. 设置分析选择1）进入Main Menu>Preference2）点取FLOTRAN CFD项3）点取OK3.定义单元类型1）进入Main Menu>Preprocessor> Element Type>Add/Edit/Delete2）点取Add3）在弹出菜单的左框中点取FLOTRAN CFD，右框中点取2D FLOTRAN 1414）点取OK5）点取Close4.生成分析区域的几何面1）生成进口段，进入Main Menu>Preprocessor>-Modeling-Create>-Areas- Rectangle>By Dimensions2）在弹出菜单中的相应区域输入以下值：X1处输入0X2处输入4Y1处输入0Y2处输入23）点取Apply4）生成出口段，再在上面弹出菜单中输入以下值：X1处输入6X2处输入12Y1处输入0Y2处输入25）点取OK6）在工具栏（Toolbar）窗口中点取SAVE_DB7）进入Main Menu>Preprocessor>-Modeling-Create>Lines>Tan to 2 Lines 8）点取左侧矩形的右面一条线作为第一条切线，再在点取菜单中点取OK 9）点取该线的上端点作为第一切点，再在点取菜单中点取OK10）点取右侧矩形的左面一条线作为第二条切线，再在点取菜单中点取OK11）点取该线的上端点作为第二切点，再在点取菜单中点取OK12）在点取菜单中点取Cancel。
《流体力学》课程实验(上机)指导书及实验报告格式
《流体力学》课程实验指导书袁守利编汽车工程学院2005年9月前言1.实验总体目标、任务与要求1)学生在学习了《流体力学》基本理论的基础上，通过伯努利方程实验、雷诺实验、阻力综合实验和动量方程实验，实现对基本理论的验证。
2）通过实验，使学生对水柱（水银柱）、U型压差计、毕托管、孔板流量计、文丘里流量计等流体力学常用的测压、测流量装置的结构、原理和使用有基本认识。
2.适用专业热能与动力工程3.先修课程《流体力学》相关章节。
4.实验项目与学时分配5. 实验改革与特色根据实验内容和现有实验条件，在实验过程中，采取学生自己动手和教师演示相结合的方法，力求达到较好的实验效果。
实验一阻力综合实验一、实验目的1．观察和测试流体稳定地在等直管道中流动及通过阀门时的能量损失情况；2．掌握管道沿程阻力系数和局部阻力系数的测定方法；3．熟悉流量的测量和测定文丘里及孔板流量计的流量系数；4．熟悉毕托管的使用。
二、实验条件阻力综合实验台三、实验原理1．实验装置：图一阻力综合实验台结构示意图1.水泵电机2.水泵3.循环储水箱4.计量水箱5.孔板及比托管实验管段进水阀6.阀门阻力实验管段进水阀7. D=14mm沿程阻力实验管段进水阀8.D=14mm沿程阻力实验管段9. 阀门阻力实验管段10.孔板流量计11. 比托管12. 测阻阀门13.测压管及测压管固定板14. D=14mm沿程阻力实验管段出水阀15阀门阻力实验管段出水阀16. 孔板及比托管实验管段出水阀17.文丘里实验管段出水阀18. D=10mm沿程阻力实验管段出水阀19.管支架20. D=10mm沿程阻力实验管段21. 文丘里流量计22排水阀门2.工作原理阻力综合实验台为多用途实验装置，利用这种实验台可进行下列实验：A 、阻力实验。
1）. 两种不同直径管路的沿程阻力实验。
2）.阀门局部阻力实验。
B 、孔板流量计流量系数和文丘里流量计流量系数的测定方法。
C 、皮托管测流速和流量的方法。
西安交大流体力学题与答案
一、 选择题(略) 二、 判断题(略) 三、 简答题1.等压面是水平面的条件是什么？:①连续介质 ② 同一介质 ③ 单一重力作用下.2. 同一容器中装两种液体，且21ρρ〈，在容器侧壁装了两根测压管。
试问：图中所标明的测压管中液面位置对吗？为什么？C (c) 盛有不同种类溶液的连通器DC D水油BB (b) 连通器被隔断AA(a) 连通容器解：不对，（右测压管液面要低一些，从点压强的大小分析） 3. 图中三种不同情况，试问：A-A 、B-B 、C-C 、D-D 中哪个是等压面？哪个不是等压面？为什么？：（ a ）A-A 是 （b ）B-B 不是 （c ）C-C 不是， D-D 是。
四、作图题(略)五、计算题(解题思路与答案)1. 已知某点绝对压强为80kN/m 2，当地大气压强p a =98kN/m 2。
试将该点绝对压强、相对压强和真空压强用水柱及水银柱表示。
解： 用水柱高表示（1）该点绝对压强：8.16mH 2o （2）该点相对压强：-1.84mH 2o（3）该点真空压强：1.84mH 2o用水银柱高表示（1）该点绝对压强：599.1mm H g（2）该点相对压强：-135.4 mm H g （3）该点真空压强：135.4 mm H g2. 一封闭水箱自由表面上气体压强p 0=25kN/m 2，h 1=5m ，h 2＝2m 。
求A 、B 两点的静水压强。
解：由压强基本公式ghp p ρ+=0求解A p = 7.551 mH 2o (74 kN/m 2) B p = 4.551 mH 2o (44.6 kN/m 2)3 如图所示为一复式水银测压计，已知m 3.21=∇，m 2.12=∇，m5.23=∇，m 4.14=∇，m5.15=∇(改为3.5m)。
试求水箱液面上的绝对压强0p ＝？解：①找已知点压强(复式水银测压计管右上端)②找出等压面③计算点压强,逐步推求水箱液面上的压强0p.: 0p=273.04 kN/m24 某压差计如图所示，已知H A=H B=1m，ΔH=0.5m。
计算流体力学上机实验报告
《计算流体力学》上机实验报告班级：姓名：学号：北京航空航天大学流体力学研究所上机实验名称两平行平板间不可压缩流体绕物体的平面无旋流动一、实验目的通过具体算例，熟悉和掌握使用CFD方法获取给定流场的流动参数。
二、实验内容、方法及步骤1.流动问题描述如下图所示，考虑在平行放置的两平板之间流过的理想不可压缩流体绕方形物体的平面无旋流动。
2.求解区域H；绕流物体是边长为2的正方形，设两平行平板之间的距离为6L。
根据流动的对称性，可取流上游来流入口位置与物体中心的距离为3动区域的四分之一作为求解区域，如下图所示。
3. 控制方程对于不可压缩流体的平面无旋流动，流函数 在区域 内满足Laplace 方程22220xy4. 边界条件（1）OABC 是一条流线，规定0OABC；（2）对 OE 上任意一点 0,P y ，有Py ；（3） ED 也是一条流线，所以2EDH ； （4）根据对称性，在 CD 上有0CDx。
5. 定解问题对于这里考虑的流动，可用下述定解问题来描述22220 , 0 , , , 20 , x yOABC y OE HED CDx在 内在上在上在上在 上6. 求解区域的离散化 － 计算网格将单位长度等分成n 份，记1h n ，于是求解区域沿x 方向可划分成M L n 个网格，用0,1,2,3,,j M 来标记；沿y 方向可划分成2HNn 个网格，用0,1,2,3,,k N 来标记。
这些网格点可分成三类：（1）当 01j L n 且 0k N ，或者当 1L njM 且 n k N 时，网格点落在流场内部，称为内点。
这些网格点上的流函数需通过求解方程组来计算； （2）当 1Ln j M 且 0k n 时，网格点落在正方形物体内部，网格点上不存在流场，无需计算；（3）其余的网格点落在流场的边界上，称为边界点。
这些网格点上的流函数直接由边界条件给定，也无需计算。
7. 定解问题的离散化 － 差分格式Laplace 方程22220xy的差分近似为1,,1,,1,,122220j kj k j kj k j k j k hh边界条件0x的差分近似为,1,0j kj kh8. 内点上流函数的计算－ 迭代算法在实际的计算中，内点的数量非常多，计算流函数需要求解大型的代数方程组。
西安交通大学传热学上机实验报告
φ1 − φ2 E= （φ1 + φ2）2
三、计算过程
用 MATLAB 编写计算程序，取网格步长 ∆x = ∆y = 0.1m 。 1、第一类边界条件 （1）运行程序 1（见附录 1） ，得到等温边界条件下计算墙角温度分布图：
图 4 等温边界条件下计算等温线分布（左图中每两条线间隔为三摄氏度） 运行程序 2（见附录 2） ，得到等温边界条件下实测墙角温度分布图：
s1=0; for i=2:11 s1=s1+(30-T(i,2))*0.53; end for j=2:15 s1=s1+(30-T(11,j))*0.53; end s1=s1+(30-T(1,2))*0.53/2+(30-T(11,16))*0.53/2
%墙角外侧换热量
s2=0; for i=2:6 s2=s2+T(i,5)*0.53; end for j=7:15 s2=s2+T(8,j)*0.53; end s2=s2+T(1,5)*0.53/2+T(8,16)*0.53/2+T(7,5)*0.53/2+T(8,6)*0.53/2 %墙角内侧换热量 s=2*(s1+s2) %单位长度墙壁的总换热量 e=abs(s1-s2)/((s1+s2)/2)
图3
内节点和绝热边界
图 3 所示的内节点和绝热边界节点方程如下： 内节点：
⎡(t −t )∆x (t −t )∆x (t −t )∆y (t −t )∆y⎤ ΦN +ΦS +ΦE +ΦW = λ⋅1⋅ ⎢ i, j+1 i, j + i, j−1 i, j + i+1, j i, j + i−1, j i, j ⎥ = 0 ∆y ∆y ∆x ∆x ⎣ ⎦
西安交通大学流体力学上机实验
流体力学上机实验报告题6-27强度为 24 m2/s 的点源位于坐标原点，与速度为 10m/s沿x正向的均匀流动叠加。
求复合流动滞止点到坐标原点的距离：该流动可以表示的绕流物体的形状及在θ=π/2 和3π/2 处物体外廓线之间的距离；θ=π/2 处物体外廓线上的速度。
分析：复合流动的速度势函数与流函数分别为Φ=10r cosθ+242πln rΨ=10r sinθ+242πθ流线方程为C=10r sinθ+242πθ势流速度场为V r=10r cosθ+242πrVθ=−10sinθ当θ=π2时V r=6.366m/s, Vθ=−10m/s则当Vr=0时，X=r=−0.382m外轮廓线所对应方程为12=10r sinθ+242πθ则在θ=π2，3π2之间时：外轮廓线距离为1.2m下面采用EXCEL工具绘制轮廓线：图中A代表θ，B代表X值，C代表Y值B= (1.2-0.38217*A)/TAN(A)C==1.2-0.38217*A所绘图像如下题7-36设有一虹吸管，其装置如图所示，设管径为 150mm.试求通过该管的流量。
已知水温为10 ℃，Δ=1mm,ξ弯头= 0.4 ,ξ进口= 1.0 。
分析：根据题意列出伯努利方程：p1ρg +z1+V122g=p2ρg+z2+V222g+ℎLT根据题意代入数据得：ℎLT=3−V22 2g考虑沿程和局部水力损失：ℎLT=(f ld+1+2∗0.4)V222g整理得：58.86=（70f+2.8）V22雷诺数Re有：Re=ρV2D μ当温度是10摄氏度时：ρ=999.7kg/m3,μ=1.307*10-3 N*s/m2 则整理上述方程代入数据可得：Re=114730√58.8670f+2.8根据科尔布鲁克公式1√f =−2.0lg (∆D⁄3.7+2.51Re√f)由上式即可构成迭代循环。
其中，对于循环初值，由,在大雷诺数条件下可在莫迪图上查出f 的初值为0.033.最后由Q=AV=πD＾2V/4可解得流量值。
其他系统西安交通大学---流体力学所有答案
其他系统西安交通大学---流体力学所有答案应用总流伯努利方程解题时，两个断面间一定是缓变流，方程才成立。
答案是：错误液体在重力场中作加速直线运动时，其自由面与重力和惯性力的合力处处正交。
答案是：错误液滴内的压强比大气压小。
（答案是：错误相邻两流线的函数值之差，是此两流线间的单宽流量。
（）答案是：正确相对静止状态的等压面一定也是水平面。
答案是：错误相对静止状态的等压面可以是斜面或曲面。
答案是：正确温度升高液体的表面张力系数增大。
答案是：错误温度升高液体的表面张力系数增大答案是：错误通过一个曲面上的体积流量与法向速度有关。
答案是：正确声音传播过程是一个等熵过程。
答案是：正确如果流场中若干流体微团无绕自身轴线旋转运动，刚称为无旋流动。
答案是：错误气体的粘性随温度的升高而增大。
答案是：正确普朗特混合长度理论建立了脉动速度与时均速度梯度之间的关系。
答案是：正确平面无旋流动既存在流函数又存在势函数。
平面流只存在流函数，无旋流动存在势函数。
答案是：正确牛顿流体的粘性切应力与速度梯度，即角变形速率成正比。
答案是：正确马赫线是超音速流动中被扰动和未扰动区域的分界线。
答案是：正确流体质点只有质量没有大小答案是：错误流体静止或相对静止状态的等压面一定是水平面。
答案是：错误流体的粘滞性随温度变化而变化，温度升高粘滞性减少；温度降低粘滞性增大。
（）答案是：错误流体的静压是指流体的点静压。
答案是：正确理想流体定常流动，流线与等势线重合。
答案是：错误静止流体中任意点的压强与方向无关。
答案是：正确静止的流体中任意一点的各个方向的压强值均相等。
答案是：正确恒定流是迁移加速度为零。
答案是：错误附面层外的流体流动时理想无旋流动。
（）答案是：正确附面层内的流体流动是粘性有旋流动。
答案是：正确附面层分离只能发生在增压减速区。
对于无旋流动，其速度势函数满足拉普拉斯方程的必要条件是流动定常。
答案是：错误对于不可压缩流体的平面流动，存在着流函数。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
题6-27
强度为 24 m2/s 的点源位于坐标原点，与速度为 10m/s沿x正向的均匀流动叠加。
求复合流动滞止点到坐标原点的距离：该流动可以表示的绕流物体的形状及在θ=π/2 和3π/2 处物体外廓线之间的距离；θ=π/2 处物体外廓线上的速度。
分析：复合流动的速度势函数与流函数分别为
Φ=10r cosθ+24
2πln rΨ=10r sinθ+24
2π
θ
流线方程为C=10r sinθ+24
2π
θ
势流速度场为V r=10r cosθ+24
2πr
Vθ=−10sinθ
当θ=π
2
时V r=6.366m/s, Vθ=−10m/s
则当Vr=0时，X=r=−0.382m
外轮廓线所对应方程为12=10r sinθ+24
2π
θ
则在θ=π
2，3π
2
之间时：外轮廓线距离为1.2m
下面采用EXCEL工具绘制轮廓线：
图中A代表θ，B代表X值，C代表Y值B= (1.2-0.38217*A)/TAN(A)
C==1.2-0.38217*A
所绘图像如下
题7-36
设有一虹吸管，其装置如图所示，设管径为 150mm.试求通过该管的流量。
已知水温为10 ℃，Δ=1mm,ξ弯头= 0.4 ,ξ进口= 1.0 。
分析：根据题意列出伯努利方程：
p1ρg +z1+
V12
2g
=
p2
ρg
+z2+
V22
2g
+ℎLT
根据题意代入数据得：
ℎLT=3−V22 2g
考虑沿程和局部水力损失：
ℎLT=(f l
d
+1+2∗0.4)
V22
2g
整理得：
58.86=（70f+2.8）V22雷诺数Re有：
Re=ρV2D μ
当温度是10摄氏度时：ρ=999.7kg/m3,
μ=1.307*10-3 N*s/m2 则整理上述方程代入数据可得：
Re=114730√58.86
70f+2.8根据科尔布鲁克公式
1√f =−2.0lg (
∆
D⁄
3.7
+
2.51
Re√f
)
由上式即可构成迭代循环。
其中，对于循环初值，由,在大雷诺数条件下可在莫迪图上查出f 的初值为0.033.
最后由Q=AV=πD＾2V/4可解得流量值。
迭代程序（fortran90语言)及运行结果如下。
PROGRAM GUANDAO
IMPLICIT NONE
REAL V,RE,F,F0,A,Q
REAL,PARAMETER::PI=3.141592653
F0=0.033
A=1
PRINT*," V F Q"
DO WHILE(ABS(A)>0.0000000000000001)
V=((2*9.8*3)/((10.5/0.15)*F0+2.8))**0.5
RE=(999.7*V*0.15)/(1.307/1000)
F=(-2*LOG10(1/(150*3.7)+2.51/(RE*F0**0.5)))**-2
A=F-F0
F0=F
Q=V*PI*0.15**2/4
PRINT*,V,F,Q
END DO
Q=V*PI*0.15**2/4
PRINT*,"---------------------------------------------" PRINT*," Q=",Q
END PROGRAM。