一元流体动力学基础
一元流体动力学基础
1.直径为150mm 的给水管道，输水量为h kg /7.980，试求断面平均流速。
解：由流量公式vA Q ρ= 注意：()vA Q s kg h kg ρ=?→//
A
Q
v ρ=
得：s m v /0154.0= 2.断面为300mm ×400mm 的矩形风道,风量为2700m 3
/h,求平均流速.如风道出口处断面收缩为150mm ×400mm,求该断面的平均流速
解：由流量公式vA Q = 得：A
Q v =
由连续性方程知2211A v A v = 得：s m v /5.122=
3.水从水箱流经直径d 1=10cm,d 2=5cm,d 3=2.5cm 的管道流入大气中. 当出口流速10m/ 时,求(1)容积流量及质量流量;(2)1d 及2d 管段的流速 解
：
(1)
由
s m A v Q /0049.0333==
质量流量s kg Q /9.4=ρ (2)由连续性方程：
33223311,A v A v A v A v ==
得：s m v s m v /5.2,/625.021==
4.设计输水量为h kg /1.2942的给水管道，流速限制在9.0∽s m /4.1之间。试确定管道直径，根据所选直径求流速。直径应是mm
50
的倍数。
解：vA Q ρ= 将9.0=v ∽s m /4.1代入得343.0=d ∽m 275.0 ∵直径是mm 50的倍数，所以取m d 3.0= 代入vA Q ρ= 得m v 18.1=
5.圆形风道，流量是10000m 3
/h,，流速不超过20 m/s 。试设计直径，根据所定直径求流速。直径规定为50 mm 的倍数。 解：vA Q = 将s m v /20≤代入得：mm d 5.420≥ 取mm d 450= 代入vA Q = 得：s m v /5.17=
6.在直径为d 圆形风道断面上，用下法选定五个点，以测局部风速。设想用和管轴同心但不同半径的圆周，将全部断面分为中间是圆，其他是圆环的五个面积相等的部分。测点即位于等分此部分面积的圆周上，这样测得的流速代表相应断面的平均流速。（1）试计算各测点到管心的距离，表为直径的倍数。（2）若各点流速为54321u u u u u ，，，，，空气密度为ρ，求质量流量G 。
解：（1）由题设得测点到管心的距离依次为1r ……5r
∵
10
310222
1S r S
r =
=
ππ 4
2
d S π=
∴
d
r d r 10
2310221=
=
同理d
r 10
253=
d r 10
274=
d
r 10
295=
（2））
（5125
14
u u d v S G ????????+==πρρ 7.某蒸汽管干管的始端蒸汽流速为25 m/s ，密度为2.62 kg/
m 3
.干管前段直径为50 mm ，接出直径40 mm 支管后，干管后段直径改为45 mm 。如果支管末端密度降为2.30 kg/m 3，干管后段末端密度降为2.24 kg/m 3，但两管质量流量相同，求两管终端流速。 解：由题意可得支
干终干始支
干）（）（）（vA vA vA Q Q ρρρ===2
1
得：???==s m v s
m v /8.22/18支
干终
8.空气流速由超音速过渡到亚超音速时，要经过冲击波。如果在冲击波前，风道中流速为=v 660 m/s ，密度为=ρ 1 kg/m 3。冲击波后速度降至=v 250 m/s 。求冲击波后的密度。 解：2211Q Q ρρ= 又面积相等
32
1
12/64.2m kg v v ==
ρρ
9.管道由不同直径的两管前后相连接组成，小管直径A d =0.2 m ，大管直径B d =0.4 m 。水在管中流动时，A 点压强A p =70kN/ m 2
，B 点压强B p =40kN/ m 2
。B 点流速B v =1 m/s 。试判断水在管中流动方向。并计算水流经过两断面间的水头损失。 解：设水流方向B A → 由连续性方程知：
B B A A A v A v =
得：s m v A /4= 由能量方程知：
12222220h Z g
v p g v p B
B A A
+++=++γγ
得：0824.212>=m h ∴水流方向B A →
10.油沿管线流动,A 断面流速为2 m/s,不记损失,求开口C 管中的液面高度。
解：由连续性方程知：2211A v A v = 得：s m v /5.42= 由能量方程得：
g
v p g v p 2022.122
2211
++=++γγ
其中：m p 5.11
=γ
代入数据解得：m p 86.12
=γ
11.水沿管线下流,若压力计的读数相同,求需要的小管直径
0d ，不计损失。
解：由连续性方程：0011A v A v =
由能量方程得
g
v 221+3=
g
v 220
得面积m d A 12.000=?
12.用水银比压计量测管中水流,过流断面中点流速u 。如图，测得A 点的比压计读数h ?=60mm 汞柱。(1)求该点的流速u ，(2)若管中流体密度为0.8g/cm 3
的油,h ?不变,该点流速为若干,不计损失。
解：设水银容重为‘γ （1）1u =()g h
ρ
ρρ?-’2=3.85m/s
（2）2u =
()
g
h
2
22ρρρ?-‘=4.31m/s
13.水由图中喷嘴流出,管嘴出口d =75 mm ，不考虑损失，计算
H 值以m 计，p 值2/m kN 计。
解：gH
v 23=
由连续性方程得：332211A v A v A v == 由1—2
断面列能量方程：g
v P g v P Z 222
2
2211
+
=++γγ
由断面压强公式：22211175.0)175.0(Z P Z Z P 水汞水γγγ+?+=+++
列水箱水面至喷口处方程：g
v H 223
=
得：m H 8.11= 列压力表至喷口处方程：
g
v g v P 22232
2
2
=+γ 得：2/79m kN P =
14.计算管线流量,管出口d=50mm,求出A,B,C,D 各点的压强,不计水头损失
解：对出口D ，4222?==g h g v d ? 由连续性方程知c b a v v v ==
又d d a
a A v A v = 得：a d v v 9=
由A →D 列能量方程
g
v g v p d
a a
2032022+
+=++γ
得：2/68m kN p a =
同理可得：2
/48.0m kN p b -=2/1.20m kN p c -=0=d p
15.水由管中铅直流出,求流量及测压计读数.水流无损失 解：设水容重为γ，水银容重为1γ 由连续性方程b b a a A v A v =
2
05.04
）
（π
=a A 001.06.0??==πδπd A b
由能量方程知g
v g v b
a 2002322+
+=+
解得：=Q s m A v a a /00815.03=
列管口到水银测压计出口的方程：
h g
v a
12
25.13γγγ=++）（
得：mm h 395=汞柱
16.同一水箱上下两孔口出流,求证:在射流交点处, h 1y 1=h 2y 2
解：1h =
g v 22
1，g
v h 222
2=
21121gt y =
，2
2
22
1gt y = ∵2211t v t v = ∴2211y h y h =
17.一压缩空气罐与文丘里式的引射管连接h d d ，，21均为已知,
问气罐压强0p 多大方才能将B 池水抽空出。
解：设水的容重为γ，g v p 2220γ= g
v g v p 222
2
211=
+γ
h p γ-=1 2211A v A v =
得： 14
1
2
0-???
? ??≥
d d h
p γ
18.如图，闸门关闭时的压力表读数为49 kN/ 2m ，闸门打开后，压力表读数为0.98 kN/ m 2
，由管进口到闸门的水头损
失为1 m ，求管中的平均流速。
流体力学龙天渝课后答案第三章一元流体动力学基础
第三章 一元流体动力学基础 1.直径为150mm 的给水管道，输水量为h kN /7.980，试求断面平均流速。 解：由流量公式vA Q ρ= 注意：()vA Q s kg h kN ρ=?→// A Q v ρ= 得：s m v /57.1= 2.断面为300mm ×400mm 的矩形风道,风量为2700m 3/h,求平均流速.如风道出口处断面收缩为150mm ×400mm,求该断面的平均流速 解：由流量公式vA Q = 得：A Q v = 由连续性方程知2211A v A v = 得：s m v /5.122= 3.水从水箱流经直径d 1=10cm,d 2=5cm,d 3=2.5cm 的管道流入大气中. 当出口流速10m/ 时,求 (1)容积流量及质量流量;(2)1d 及2d 管段的流速 解：(1)由s m A v Q /0049.0333== 质量流量s kg Q /9.4=ρ (2)由连续性方程： 33223311,A v A v A v A v == 得：s m v s m v /5.2,/625.021== 4.设计输水量为h kg /294210的给水管道，流速限制在9.0∽s m /4.1之间。试确定管道直径，根据所选直径求流速。直径应是mm 50的倍数。 解：vA Q ρ= 将9.0=v ∽s m /4.1代入得343.0=d ∽m 275.0 ∵直径是mm 50的倍数，所以取m d 3.0= 代入vA Q ρ= 得m v 18.1= 5.圆形风道，流量是10000m 3/h,，流速不超过20 m/s 。试设计直径，根据所定直径求流速。直径规定为50 mm 的倍数。 解：vA Q = 将s m v /20≤代入得：mm d 5.420≥ 取mm d 450= 代入vA Q = 得：s m v /5.17= 6.在直径为d 圆形风道断面上，用下法选定五个点，以测局部风速。设想用和管轴同心但不同半径的圆周，将全部断面分为中间是圆，其他是圆环的五个面积相等的部分。测点即位于等分此部分面积的圆周上，这样测得的流速代表相应断面的平均流速。（1）试计算各测点到管心的距离，表为直径的倍数。（2）若各点流速为54321u u u u u ，，，，，空气密度为ρ，求质量流量G 。
流体力学复习要点(计算公式)
D D y S x e P gh2 gh1 h2 h1 b L y C C D D y x P hc 第一章 绪论 单位质量力： m F f B m = 密度值： 3 m kg 1000=水ρ， 3 m kg 13600=水银ρ， 3 m kg 29.1=空气ρ 牛顿内摩擦定律：剪切力： dy du μ τ=， 内摩擦力：dy du A T μ= 动力粘度： ρυ μ= 完全气体状态方程：RT P =ρ 压缩系数： dp d 1dp dV 1ρρκ= -=V （N m 2 ） 膨胀系数：T T V V V d d 1d d 1ρρα - == (1/C ?或1/K) 第二章 流体静力学+ 流体平衡微分方程： 01;01;01=??-=??-=??- z p z y p Y x p X ρρρ 液体平衡全微分方程：)(zdz ydy xdx dp ++=ρ 液体静力学基本方程：C =+ +=g p z gh p p 0ρρ或 绝对压强、相对压强与真空度：a abs P P P +=；v a abs P P P P -=-= 压强单位换算：水银柱水柱mm 73610/9800012 ===m m N at 2/101325 1m N atm = 注： h g P P →→ρ ； P N at →→2m /98000乘以 2/98000m N P a = 平面上的静水总压力：（1）图算法 Sb P = 作用点e h y D +=α sin 1 ) () 2(32121h h h h L e ++= ρ 若01 =h ，则压强为三角形分布，3 2L e y D == ρ 注：①图算法适合于矩形平面；②计算静水压力首先绘制压强分布图， α 且用相对压强绘制。 （2）解析法 A gh A p P c c ρ== 作用点A y I y y C xc C D + = 矩形12 3 bL I xc = 圆形 64 4 d I xc π= 曲面上的静水总压力： x c x c x A gh A p P ρ==；gV P z ρ= 总压力z x P P P += 与水平面的夹角 x z P P arct an =θ 潜体和浮体的总压力： 0=x P 排浮gV F P z ρ== 第三章 流体动力学基础 质点加速度的表达式??? ? ? ? ??? ??+??+??+??=??+??+??+??=??+??+??+??=z u u y u u x u u t u a z u u y u u x u u t u a z u u y u u x u u t u a z z z y z x z z y z y y y x y y x z x y x x x x A Q V Q Q Q Q Q G A = === ? 断面平均流速重量流量质量流量体积流量g udA m ρρ 流体的运动微分方程： t z t y t x d du z p z d du y p Y d du x p X = ??-=??-=??- ρρρ1;1;1 不可压缩流体的连续性微分方程 ： 0z u y u x u z y x =??+??+?? 恒定元流的连续性方程： dQ A A ==2211d u d u 恒定总流的连续性方程：Q A A ==2211νν 无粘性流体元流伯努利方程：g 2u g p z g 2u g p z 2 2 222 111++=++ρρ 粘性流体元流伯努利方程： w 2 2222111'h g 2u g p z g 2u g p z +++=++ρρ
工程流体力学(水力学)闻德第五章-实际流体动力学基础课后答案
工程流体力学闻德课后习题答案 第五章 实际流体动力学基础 5—1设在流场中的速度分布为u x =2ax ，u y =-2ay ，a 为实数，且a >0。试求切应力τxy 、τyx 和附加压应力p ′x 、p ′y 以及压应力p x 、p y 。 解：0y x xy yx u u x y ττμ??? ?==+= ????? 24x x u p a x μμ?'=-=-?，24y y u p a y μμ?'=-=?， 4x x p p p p a μ'=+=-，4y y p p p p a μ'=+=+ 5－2 设例５－1中的下平板固定不动，上平板以速度v 沿x 轴方向作等速运动（如图 所示），由于上平板运动而引起的这种流动，称柯埃梯（Couette ）流动。试求在这种流动情况下，两平板间的速度分布。（请将 d 0d p x =时的这一流动与在第一章中讨论流体粘性时的流动相比较） 解：将坐标系ox 轴移至下平板，则边界条件为 y ＝０，0X u u ==；y h =，u v =。 由例５－1中的（11）式可得 2d (1)2d h y p y y u v h x h h μ=- - （１） 当d 0d p x =时，y u v h =，速度ｕ为直线分布，这种特殊情况的流动称简单柯埃梯流动或简单剪切流动。它只是由于平板运动，由于流体的粘滞性带动流体发生的流动。 当 d 0d p x ≠时，即为一般的柯埃梯流动，它是由简单柯埃梯流动和泊萧叶流动叠加而成，速度分布为 (1)u y y y p v h h h =-- （２） 式中2d ()2d h p p v x μ= - （３） 当p ＞０时，沿着流动方向压强减小，速度在整个断面上的分布均为正值；当p ＜０时，沿流动方向压强增加，则可能在静止壁面附近产生倒流，这主要发生p ＜－１的情况． 5－3 设明渠二维均匀（层流）流动，如图所示。若忽略空气阻力，试用纳维—斯托克斯方程和连续性方程，证明过流断面上的速度分布为2sin (2)2 x g u zh z ，单宽流量 3 sin 3 gh q 。
流体力学计算公式
C3.6.2 达西摩擦因子 为了确定λ与Re 的关系,人们作了大量实验和理论研究，下面介绍有代表性的结果。 1．尼古拉兹实验 尼古拉兹（J.Nikuradse，1932）分析了达西的圆管沿程阻力实验数据后，发现壁面粗糙度对λ的影响很大，决定用人工粗糙度方法实现对粗糙度的控制。他用当地黄砂砂粒经筛选后分类均匀粘贴在管内壁上，相对粗糙度ε/d 从1/30—1/1014分6种，测得λ与Re 的关系，得到尼古拉兹图（图C3.6.1）。 2． 常用计算公式 从尼古拉兹图中看到在不同Re 数和ε/d 值的区域，λ有不同的变化规律。 图C3.6.1
（1）层流区 由泊肃叶定律推导的沿程水头损失（C3.4.10）式可得 代入达西公式（C3.6.3）式，可得层流区λ的解析式 上式表明层流区λ与管壁粗糙度无关，写成常用对数形式为 上式在双对数坐标系中是一条直线，与尼古拉兹图吻合。 （2）过渡区 该区是层流向湍流的转捩区（2000 1、单位质量力：m F f B B = 2、流体的运动粘度：ρ μ=v （μ[动力]粘度，ρ密度） 3、压缩系数：dp d dp dV V ρρκ?=?-=11（κ的单位是N m 2）体积模量为压缩系数的倒数 4、体积膨胀系数：dT d dT dV V v ρρα?-=?=11（v α的单位是C K ?1,1） 5、牛顿内摩擦定律：为液体厚）为运动速度，以应力表示为y u dy du dy du A T (,μτμ== 6、静止液体某点压强：为该点到液面的距离）h gh p z z g p p ()(000ρρ+=-+= 7、静水总压力： )h (为受压面积，为受压面形心淹没深度为静水总压力，A p ghA A p p c ρ== 8、元流伯努利方程;'2221112w h g p z g u g p z ++=++ρρ('w h 为粘性流体元流单位重量流体由过流断面1-1运动至过流断面2-2的机械能损失，z 为某点的位置高度或位置水头，g p ρ为测压管高度或压强水头，g u ρ2是单位流体具有的动能，u gh g p p g u 22'=-=ρ，u gh C g p p g C u 22'=-=ρC 是修正系数，数值接近于1） 9、总流伯努利方程:w h g v g p z g v g p z +++=++222 221221111αραρ(α为修正系数通常取1） 10、文丘里流量计测管道流量：)21)(41()()(42 122211g d d d k h k g p z g p z k Q -=?=+-+=πμρρμ 11、沿程水头损失一般表达式：g v d l h f 22 λ=（l 为管长，d 为管径，v 为断面平均流速，g 工程流体力学公式总结 第二章流体得主要物理性质 ?流体得可压缩性计算、牛顿内摩擦定律得计算、粘度得三种表示方法。１.密度ρ= m／V 2.重度γ= G /V 3．流体得密度与重度有以下得关系:γ= ρg或ρ= γ/ g ４.密度得倒数称为比体积,以υ表示υ= １/ ρ= V/m 5.流体得相对密度:d = γ流/γ水= ρ流/ρ水 6．热膨胀性 ７.压缩性、体积压缩率κ 8．体积模量 ９.流体层接触面上得内摩擦力 10．单位面积上得内摩擦力(切应力)(牛顿内摩擦定律) 1１.、动力粘度μ： 12.运动粘度ν:ν=μ/ρ 1３．恩氏粘度°E:°E = t 1 /t 2 第三章流体静力学 ?重点：流体静压强特性、欧拉平衡微分方程式、等压面方程及其、流体静力学基本方程意义及其计算、压强关系换算、相对静止状态流体得压强计算、流体静压力得计算（压力体)。 1.常见得质量力: 重力ΔW = Δmg、 直线运动惯性力ΔFI ＝Δm·ａ 离心惯性力ΔＦR ＝Δm·rω２、 ２.质量力为F。:F= ｍ·aｍ= m(ｆｘi+f yj＋ｆzk) am ＝F/m = f xi＋f yj＋fｚk为单位质量力,在数值上就等于加速度 实例:重力场中得流体只受到地球引力得作用，取z轴铅垂向上,xoy为水平面，则单位质量力在x、y、z轴上得分量为 fｘ= ０,ｆy＝0 , fz=-mg/ｍ= －ｇ式中负号表示重力加速度g与坐标轴z方向相反 ３流体静压强不就是矢量,而就是标量,仅就是坐标得连续函数。即:p=p(x,y,z),由此得静压强得全微分为： 4.欧拉平衡微分方程式 单位质量流体得力平衡方程为: 思考题及答案 一、选择 (1) 二、例题 (2) 三、问答 (14) 一、选择 问题：恒定流是： A、流动随时间按一定规律变化； B、流场中任意空间点的运动要素不随时间变化； C、各过流断面的速度分布相同； D、各过流断面的压强相同。 问题：非恒定流是： A、； B、； C、； D、。 问题：一元流动是： A、均匀流； B、速度分布按直线变化； C、运动参数是一个空间坐标和时间变量的函数； D、限于直线流动。 问题：均匀流是： A、当地加速度为零； B、迁移加速度为零； C、向心加速度为零； D、合加速度为零。 问题1：流速势函数存在的必要与充分条件是: A、平面无旋流动； B、理想流体平面流动； C、不可压缩流体平面流动； D、无旋流动。 问题2：设流速势函数j=xyz，则点B（1，2，1）处的速度u 为： B A、5； B、1； C、3； D、2。 判断：公式（3－14）与公式（3－16）两式形式完全相同，因此其应用条件也相同。 你的回答：对错 判断：土坝渗流中的流网网格一定是直线正方形网格。 你的回答：对错 二、例题 例1如图3-7，已知流速场为，其中C为常数，求流 线方程。 解：由式得 图3-7 积分得： 则： 此外，由得： 因此，流线为Oxy平面上的一簇通过原点的直线，这种流动称为平面点源流动（C＞0时）或平 面点汇流动（C＜0时） 例2已知平面流动 试求：（1）t=0时，过点M（-1，-1）的流线。 （2）求在t=0时刻位于x=-1，y=-1点处流体质点的迹线。解：（1）由式 （2）由式 得 得 得： 由t=0时，x=-1，y=-1得C 1=0, C 2 =0，则有： 将：t=0，x=-1，y=-1 代入得瞬时流线 xy=1 最后可得迹线为： 即流线是双曲线。 例3已知流动速度场为 第五章 实际流体动力学基础 5—1设在流场中的速度分布为u x =2ax ，u y =-2ay ，a 为实数，且a >0。试求切应力τxy 、τyx 和附加压应力p ′x 、p ′y 以及压应力p x 、p y 。 解：0y x xy yx u u x y ττμ??? ?==+= ????? 24x x u p a x μμ?'=-=-?，24y y u p a y μμ?'=-=?， 4x x p p p p a μ'=+=-，4y y p p p p a μ'=+=+ 5－2 设例５－1中的下平板固定不动，上平板以速度v 沿x 轴方向作等速运动（如图所示），由于上平板运动而引起的这种流动，称柯埃梯（Couette ）流动。试求在这种流动情况下，两平板间的速度分布。（请将 d 0d p x =时的这一流动与在第一章中讨 论流体粘性时的流动相比较） 解：将坐标系ox 轴移至下平板，则边界条件为 y ＝０，0X u u ==；y h =，u v =。 由例５－1中的（11）式可得 2d (1)2d h y p y y u v h x h h μ=- - （１） 当d 0d p x =时，y u v h =，速度ｕ为直线分布，这种特殊情况的流动称简单柯埃梯流动或简单剪切流动。它只是由于平板运动，由于流体的粘滞性带动流体发生的流动。 当d 0d p x ≠时，即为一般的柯埃梯流动，它是由简单柯埃梯流动和泊萧叶流动叠加而成，速度分布为 (1)u y y y p v h h h =-- （２） 式中2d ()2d h p p v x μ= - （３） 当p ＞０时，沿着流动方向压强减小，速度在整个断面上的分布均为正值；当p ＜０时，沿流动方向压强增加，则可能在静止壁面附近产生倒流，这主要发生p ＜－１的情况． 5－3 设明渠二维均匀（层流）流动，如图所示。若忽略空气阻力，试用纳维—斯托克斯方程和连 续性方程，证明过流断面上的速度分布为2 sin (2)2x g u zh z r q m = -，单宽流量3sin 3gh q r q m =。 解：（1）因是恒定 二维流动， 0y x z u u u t t t ???===抖?，u u x =，0y u =， 0z u =，由纳维——斯托克 斯方程和连续性方程可 得 计算流体力学常用数值方法简介 李志印 熊小辉 吴家鸣 (华南理工大学交通学院) 关键词 计算流体力学 数值计算 一 前 言 任何流体运动的动力学特征都是由质量守恒、动量守恒和能量守恒定律所确定的,这些基本定律可以由流体流动的控制方程组来描述。利用数值方法通过计算机求解描述流体运动的控制方程,揭示流体运动的物理规律,研究流体运动的时一空物理特征,这样的学科称为计算流体力学。 计算流体力学是一门由多领域交叉而形成的一门应用基础学科,它涉及流体力学理论、计算机技术、偏微分方程的数学理论、数值方法等学科。一般认为计算流体力学是从20世纪60年代中后期逐步发展起来的,大致经历了四个发展阶段:无粘性线性、无粘性非线性、雷诺平均的N-S方程以及完全的N-S方程。随着计算机技术、网络技术、计算方法和后处理技术的迅速发展,利用计算流体力学解决流动问题的能力越来越高,现在许多复杂的流动问题可以通过数值计算手段进行分析并给出相应的结果。 经过40年来的发展,计算流体力学己经成为一种有力的数值实验与设计手段,在许多工业领域如航天航空、汽车、船舶等部门解决了大量的工程设计实际问题,其中在航天航空领域所取得的成绩尤为显著。现在人们已经可以利用计算流体力学方法来设计飞机的外形,确定其气动载荷,从而有效地提高了设计效率,减少了风洞试验次数,大大地降低了设计成本。此外,计算流体力学也己经大量应用于大气、生态环境、车辆工程、船舶工程、传热以及工业中的化学反应等各个领域,显示了计算流体力学强大的生命力。 随着计算机技术的发展和所需要解决的工程问题的复杂性的增加,计算流体力学也己经发展成为以数值手段求解流体力学物理模型、分析其流动机理为主线,包括计算机技术、计算方法、网格技术和可视化后处理技术等多种技术的综合体。目前计算流体力学主要向二个方向发展:一方面是研究流动非定常稳定性以及湍流流动机理,开展高精度、高分辩率的计算方法和并行算法等的流动机理与算法研究;另一方面是将计算流体力学直接应用于模拟各种实际流动,解决工业生产中的各种问题。 二 计算流体力学常用数值方法 流体力学数值方法有很多种,其数学原理各不相同,但有二点是所有方法都具备的,即离散化和代数化。总的来说其基本思想是:将原来连续的求解区域划分成网格或单元子区 第三章 流体动力学基础 习 题 一、单选题 1、在稳定流动中，在任一点处速度矢量是恒定不变的，那么流体质点是 （ ） A ．加速运动 B ．减速运动 C ．匀速运动 D ．不能确定 2、血管中血液流动的流量受血管内径影响很大。如果血管内径减少一半，其血液的流量将变为原来的（ ）倍。 A ．21 B ．41 C ．81 D ．161 3、人在静息状态时，整个心动周期内主动脉血流平均速度为0.2 m/s ，其内径d =2×10-2 m ，已知血液的粘度η =×10-3 Pa·S，密度ρ=×103 kg/m 3 ，则此时主动脉中血液的流动形态处于（ ）状态。 A ．层流 B ．湍流 C ．层流或湍流 D ．无法确定 4、正常情况下，人的小动脉半径约为3mm ，血液的平均速度为20cm/s ，若小动脉某部分被一硬斑阻塞使之变窄，半径变为2mm ，则此段的平均流速为（ ）m/s 。 A ．30 B ．40 C ．45 D ．60 5、有水在同一水平管道中流动，已知A 处的横截面积为S A =10cm 2 ，B 处的横截面积为 S B =5cm 2，A 、B 两点压强差为1500Pa ，则A 处的流速为（ ）。 A ．1m/s B ．2m/s C ．3 m/s D ．4 m/s 6、有水在一水平管道中流动，已知A 处的横截面积为S A =10cm 2 ，B 处的横截面积为S B =5cm 2 ，A 、B 两点压强之差为1500Pa ，则管道中的体积流量为（ ）。 A ．1×10-3 m 3 /s B ．2×10-3 m 3 /s C ．1×10-4 m 3 /s D ．2×10-4 m 3 /s 7、通常情况下，人的小动脉内径约为6mm ，血流的平均流速为20cm/s ，若小动脉某处被一硬斑阻塞而变窄，测得此处血流的平均流速为80cm/s ，则小动脉此处的内径应为（ ）mm 。 A ．4 B ．3 C ．2 D ．1 8、正常情况下，人的血液密度为×103 kg/m 3 ，血液在内径为6mm 的小动脉中流动的平均速度为20cm/s ，若小动脉某处被一硬斑阻塞而变窄，此处内径为4mm ，则小动脉宽处与窄处压强之差（ ）Pa 。 二、判断题 第三章 一元流体动力学基础 1.直径为150mm 的给水管道，输水量为h kN /7.980，试求断面平均流速。 解：由流量公式vA Q ρ= 注意：()vA Q s kg h kN ρ=?→// A Q v ρ= 得：s m v /57.1= 2.断面为300mm ×400mm 的矩形风道,风量为2700m 3/h,求平均流速.如风道出口处断面收缩为150mm ×400mm,求该断面的平均流速 解：由流量公式vA Q = 得：A Q v = 由连续性方程知2211A v A v = 得：s m v /5.122= 3.水从水箱流经直径d 1=10cm,d 2=5cm,d 3=2.5cm 的管道流入大气中. 当出口流速10m/ 时,求 (1)容积流量及质量流量;(2)1d 及2d 管段的流速 解：(1)由s m A v Q /0049.0333== 质量流量s kg Q /9.4=ρ (2)由连续性方程： 33223311,A v A v A v A v == 得：s m v s m v /5.2,/625.021== 4.设计输水量为h kg /294210的给水管道，流速限制在9.0∽s m /4.1之间。试确定管道直径，根据所选直径求流速。直径应是mm 50的倍数。 解：vA Q ρ= 将9.0=v ∽s m /4.1代入得343.0=d ∽m 275.0 ∵直径是mm 50的倍数，所以取m d 3.0= 代入vA Q ρ= 得m v 18.1= 5.圆形风道，流量是10000m 3/h,，流速不超过20 m/s 。试设计直径，根据所定直径求流速。直径规定为50 mm 的倍数。 解：vA Q = 将s m v /20≤代入得：mm d 5.420≥ 取mm d 450= 代入vA Q = 得：s m v /5.17= 6.在直径为d 圆形风道断面上，用下法选定五个点，以测局部风速。设想用和管轴同心但不同半径的圆周，将全部断面分为中间是圆，其他是圆环的五个面积相等的部分。测点即位于等分此部分面积的圆周上，这样测得的流速代表相应断面的平均流速。（1）试计算各测点到管心的距离，表为直径的倍数。（2）若各点流速为54321u u u u u ，，，，，空气密度为ρ，求质量流量G 。 流体力学的计算公式 众所周知，锅炉之类的热力设备可以采用水封管作为限压装置。设备工作时，水封管内的水将设备与大气隔绝，使设备保持一定的压力。水封管内水的多少，即水封管内的实际水位的高低与设备汽压大小有关。汽压大，水位就高；汽压小，水位就低。当设备压力升至额定蒸汽压力时，水封管内的水位应该达到最高水位；而在设备压力稍微大于额定蒸汽压力时，水封管内的水应能立即被冲除掉，使设备内的蒸汽能够迅速得到排放。但是，水封管的内径必须足够大；以保证它的排汽能力大于设备的最大蒸发量，从而防止设备发生超压事故。《蒸汽锅炉安全技术监察规程》第132条规定：“对于额定蒸汽压力小于等于0.1MPa的锅炉可以采用静重式安全阀或水封式安全装置。水封装置的水封管内径不应小于25mm”。《小型和常压热水锅炉安全监察规定》第二十八条要求：“水封管的内径应当根据锅炉的额定容量和压力确定，且内径不得小于25毫米”。 但是，水封管的内径究竟应该根据什么样的公式去计算？这里没有具体加以规定。水封管的内径为何不能小于25mm？这里没有加以说明。 1 确定水封管内径的一种错误观点 有人说，可以按安全阀的排放量计算公式去间接确定水封管的内径。笔者认为这是不对的。理由如下： 其一，按流体力学，安全阀的局部阻碍因为其进出口呈直角型式，边壁是突变的，主流与边壁之间形成大尺度旋涡，蒸汽排放时能量损失很大；而水封管的边壁是渐变的，又不出现减速增压现象的部位，故蒸汽排放时的能量损失很小。其二，安全阀的出口之外肯定还有排汽管，它的排放量计算公式应该考虑这一额外的因素；而水封管却只相当于它的排汽管。其三，安全阀在排汽时，汽流还自始至终受到弹簧或者重锤施加的反方向作用力，能量损失非常之大；而水封管在排汽时却不会受到这样的反作用力，无此项能量损失。显然，安全阀排汽与水封管排汽区别太大，不能按安全阀排放量计算公式去确定水封管的内径。 2 用流体力学理论建立水封管内径计算方法 为了正确确定水封管的应有内径，我们有必要按流体力学的理论来分析一下水封管的流动阻力和能量损失。蒸汽流过水封管时的阻力损失ΔP等于水封管进口压力P1与出口压力P2之差。管道的阻力损失 ΔP＝Pa（1） 式中λ—管道摩擦系数； d —管道内径，m； ρ—流体的密度，k g／m3； u —流体在管内的流速，m／s； L—平直管段的总长度，m； Σξ—管道中各种局部阻力之和。 管道的摩擦系数λ值取决于流体流动的雷诺数Re 和管壁的相对粗糙度。管壁的相对粗糙度等于管壁的绝对粗糙度（即管壁内凸起高度）K与管内径d之比。雷诺数Re =(2) 式中ν－流体运动粘滞系数，m2／s； 其他符号说明同上。 由于水封管进出口压差通常大于40000Pa，排汽时的流速都大于100m／s ，雷诺数Re 都大于80000000，所以蒸汽的流动状态为紊流流动，而且位于阻力平方区。在此区域内，摩擦阻力系数λ值仅取决于管壁的相对粗糙度，可以用尼古拉兹公式进行计算：λ＝（3） 水 水银 题1图 1 2 3 题型一：曲面上静水总压力的计算问题（注：千万注意方向，绘出压力体） 1、AB 曲面为一圆柱形的四分之一，半径R=0.2m ，宽度（垂直纸面）B=0.8m ，水深H=1.2m ，液体密度3 /850m kg =ρ，AB 曲面左侧受到液体压力。求作用在AB 曲面上的水平分力和铅直分力。（10分） 解：（1）水平分力： RB R H g A h P z c x ?-==)2 (ργ…….(3分) N 1.14668.02.0)2 2 .02.1(8.9850=??- ??=，方向向右（2分） 。 （2）铅直分力：绘如图所示的压力体，则 B R R R H g V P z ??? ? ????+-==4)(2πργ……….（3分） 1.1542 8.042.014.32.0)2.02.1(8.98502=???? ? ?????+?-??=，方向向下（2分）。 2．有一圆滚门，长度l=10m ，直径D=4.2m ，上游水深H1=4.2m ，下游水深H2=2.1m ，求作用于圆滚门上的水平和铅直分压力。 解题思路：（1）水平分力： l H H p p p x )(2 1 222121-= -=γ 方向水平向右。 （2）作压力体，如图，则 l D Al V p z 4 432 πγγγ? === 方向垂直向上。 3.如图示，一半球形闸门，已知球门的半径m R 1= ，上下游水位差m H 1= ，试求闸门受到的水平分力和竖直分力的 大小和方向。 解： (1)水平分力： ()2R R H A h P c πγγ?+===左,2R R A h P c πγγ?=' =右 右左P P P x -= kN R H 79.30114.31807.92=???=?=πγ， 方向水平向右。 （２）垂直分力： V P z γ=，由于左、右两侧液体对曲面所形成的压力体均为半球面，且两侧方向相反，因而垂直方向总的压力为0。 4、密闭盛水容器，已知h 1=60cm,h 2=100cm ，水银测压计读值cm h 25=?。试求半径R=0.5m 的半球盖AB 所受总压力的水平分力和铅垂分力。 第一章 绪论 单位质量力： m F f B m = 密度值： 3 m kg 1000=水ρ， 3 m kg 13600=水银ρ，3 m kg 29.1=空气 ρ 牛顿内摩擦定律：剪切力： dy du μ τ=， 内摩擦力：dy du A T μ= 动力粘度：ρυμ= 完全气体状态方程：RT P =ρ 压缩系数： dp d 1dp dV 1ρρκ= -=V （N m 2 ） 膨胀系数：T T V V V d d 1d d 1ρρα - == (1/C ?或1/K) 第二章 流体静力学+ 流体平衡微分方程： 01;01;01=??-=??-=??- z p z y p Y x p X ρρρ 液体平衡全微分方程：)(zdz ydy xdx dp ++=ρ 液体静力学基本方程：C =+ +=g p z gh p p 0ρρ或 绝对压强、相对压强与真空度：a abs P P P +=；v a abs P P P P -=-= 压强单位换算：水银柱水柱m m 73610/9800012===m m N at 2/1013251m N atm = 注： h g P P →→ρ ； P N at →→2m /98000乘以 2/98000m N P a = 平面上的静水总压力：（1）图算法 Sb P = 作用点e h y D += 1 ) () 2(32121h h h h L e ++= 3 2L e y D = = （2）解析法 A gh A p P c c ρ== 作用点A y I y y C xc C D + = 矩形 12 3 bL I xc = 圆形 64 4 d I xc π= 曲面上的静水总压力： x c x c x A gh A p P ρ==；gV P z ρ= 总压力 z x P P P += 与水平面的夹角 x z P P arctan =θ 潜体和浮体的总压力： 0=x P 排浮gV F P z ρ== 第三章 流体动力学基础 质点加速度的表达式??? ? ?? ?????+??+??+??=??+??+??+??= ??+??+??+??=z u u y u u x u u t u a z u u y u u x u u t u a z u u y u u x u u t u a z z z y z x z z y z y y y x y y x z x y x x x x A Q V Q Q Q Q Q G A = === ? 断面平均流速重量流量质量流量体积流量g udA m ρρ 流体的运动微分方程： t z t y t x d du z p z d du y p Y d du x p X = ??-=??-=??- ρρρ1;1;1 不可压缩流体的连续性微分方程 ： 0z u y u x u z y x =??+??+?? 恒定元流的连续性方程： dQ A A ==2211d u d u 恒定总流的连续性方程：Q A A ==2211νν 无粘性流体元流伯努利方程：g 2u g p z g 2u g p z 2 2 222 111++=++ρρ 粘性流体元流伯努利方程： w 2 2222111'h g 2u g p z g 2u g p z +++=++ρρ 以下为 program.cs的代码 using System; using System.Collections.Generic; using System.Linq; //using System.Threading.Tasks; using System.Windows.Forms; namespace水力计算优化版 { staticclass Program { /// 1738年瑞士数学家：伯努利在名著《流体动力学》中提出了伯努利方程。 1755年欧拉在名著《流体运动的一般原理》中提出理想流体概念，并建立了理想流体基本方程和连续方程，从而提出了流体运动的解析方法，同时提出了速度势的概念。 1781年拉格朗日首先引进了流函数的概念。 1826年法国工程师纳维，1845年英国数学家、物理学家斯托克思提出了著名的N-S方程。 1876年雷诺发现了流体流动的两种流态：层流和紊流。 1858年亥姆霍兹指出了理想流体中旋涡的许多基本性质及旋涡运动理论，并于1887年提出了脱体绕流理论。 19世纪末，相似理论提出，实验和理论分析相结合。 1904年普朗特提出了边界层理论。 20世纪60年代以后，计算流体力学得到了迅速的发展。流体力学内涵不断地得到了充实与提高。 理想势流伯努利方程 （3-14） 或（3-15） 物理意义:在同一恒定不可压缩流体重力势流中，理想流体各点的总比能相等即在整个势流场中，伯努利常数C 均相等。 （应用条件：“”所示） 符号说明 二、沿流线的积分 1.只有重力作用的不可压缩恒定流，有 2.恒定流中流线与迹线重合： 沿流线（或元流）的能量方程： （3-16） 注意：积分常数C，在非粘性、不可压缩恒定流流动中，沿同一流线保持不变。一般不同流线各不相同（有旋流）。 （应用条件：“”所示，可以是有旋流） 流速势函数（势函数）观看录像>> ?存在条件：不可压缩无旋流，即或 必要条件存在全微分d 直角坐标 （3-19） 式中：——无旋运动的流速势函数，简称势函数。 ?势函数的拉普拉斯方程形式 对于不可压缩的平面流体流动中，将（3-19）式代入连续性微分方程（3-18），有： 或（3-20） 适用条件：不可压缩流体的有势流动。 点击这里练习一下 极坐标 （3-21） 流函数 第一章 绪论 单位质量力：m F f B m = 密度值：3 m kg 1000=水ρ，3 m kg 13600=水银 ρ ，3 m kg 29.1=空气 ρ 牛顿内摩擦定律：剪切力：dy du μτ=， 内摩擦力：dy du A T μ= 动力粘度：ρυμ= 完全气体状态方程：RT P =ρ 压缩系数： dp d 1dp dV 1ρρκ= -=V （N m 2 ） 膨胀系数：T T V V V d d 1d d 1ρρα-==(1/C ?或1/K) 第二章 流体静力学+ 流体平衡微分方程：01;01;01=??-=??-=??-z p z y p Y x p X ρρρ 液体平衡全微分方程：)(zdz ydy xdx dp ++=ρ 液体静力学基本方程：C =++=g p z gh p p 0 ρρ或 绝对压强、相对压强与真空度：a abs P P P +=；v a abs P P P P -=-= 压强单位换算：水银柱水柱mm 73610/9800012===m m N at 2/101325 1m N atm = 注： h g P P →→ρ ； P N at →→2m /98000乘以 2/98000m N P a = 平面上的静水总压力：（1）图算法 Sb P = 作用点e h y D +=α sin 1 ) ()2(32121h h h h L e ++= 若01 =h ，则压强为三角形分布，3 2L e y D = = 注：①图算法适合于矩形平面；②计算静水压力首先绘制压强分布图， 且用相对压强绘制。 （2）解析法 A gh A p P c c ρ== 作用点A y I y y C xc C D + = 矩形123bL I xc = 圆形64 4d I xc π= 曲面上的静水总压力： x c x c x A gh A p P ρ==；gV P z ρ= 总压力z x P P P += 与水平面的夹角 x z P P arct an =θ 潜体和浮体的总压力：0=x P 排浮gV F P z ρ== 第三章 流体动力学基础 质点加速度的表达式??? ? ? ? ??? ??+??+??+??=??+??+??+??= ??+??+??+??=z u u y u u x u u t u a z u u y u u x u u t u a z u u y u u x u u t u a z z z y z x z z y z y y y x y y x z x y x x x x A Q V Q Q Q Q Q G A = === ? 断面平均流速重量流量质量流量体积流量g udA m ρρ 流体的运动微分方程： t z t y t x d du z p z d du y p Y d du x p X = ??-=??-=??- ρρρ1;1;1 不可压缩流体的连续性微分方程 ：0z u y u x u z y x =??+ ??+?? 恒定元流的连续性方程：dQ A A ==2211d u d u 恒定总流的连续性方程：Q A A ==2211νν 无粘性流体元流伯努利方程：g 2u g p z g 2u g p z 2 2222 111++=++ρρ 粘性流体元流伯努利方程：w 22222111'h g 2u g p z g 2u g p z +++=++ρρ 一元流体动力学基础 1.直径为150mm 的给水管道，输水量为h kg /7.980，试求断面平均流速。 解：由流量公式vA Q ρ= 注意：()vA Q s kg h kg ρ=?→// A Q v ρ= 得：s m v /0154.0= 2.断面为300mm ×400mm 的矩形风道,风量为2700m 3 /h,求平均流速.如风道出口处断面收缩为150mm ×400mm,求该断面的平均流速 解：由流量公式vA Q = 得：A Q v = 由连续性方程知2211A v A v = 得：s m v /5.122= 3.水从水箱流经直径d 1=10cm,d 2=5cm,d 3=2.5cm 的管道流入大气中. 当出口流速10m/ 时,求(1)容积流量及质量流量;(2)1d 及2d 管段的流速 解 ： (1) 由 s m A v Q /0049.0333== 质量流量s kg Q /9.4=ρ (2)由连续性方程： 33223311,A v A v A v A v == 得：s m v s m v /5.2,/625.021== 4.设计输水量为h kg /1.2942的给水管道，流速限制在9.0∽s m /4.1之间。试确定管道直径，根据所选直径求流速。直径应是mm 50 的倍数。 解：vA Q ρ= 将9.0=v ∽s m /4.1代入得343.0=d ∽m 275.0 ∵直径是mm 50的倍数，所以取m d 3.0= 代入vA Q ρ= 得m v 18.1= 5.圆形风道，流量是10000m 3 /h,，流速不超过20 m/s 。试设计直径，根据所定直径求流速。直径规定为50 mm 的倍数。 解：vA Q = 将s m v /20≤代入得：mm d 5.420≥ 取mm d 450= 代入vA Q = 得：s m v /5.17= 6.在直径为d 圆形风道断面上，用下法选定五个点，以测局部风速。设想用和管轴同心但不同半径的圆周，将全部断面分为中间是圆，其他是圆环的五个面积相等的部分。测点即位于等分此部分面积的圆周上，这样测得的流速代表相应断面的平均流速。（1）试计算各测点到管心的距离，表为直径的倍数。（2）若各点流速为54321u u u u u ，，，，，空气密度为ρ，求质量流量G 。 解：（1）由题设得测点到管心的距离依次为1r ……5r ∵ 10 310222 1S r S r = = ππ 4 2 d S π= ∴ d r d r 10 2310221= = 同理d r 10 253= d r 10 274= d r 10 295= （2）） （5125 14 u u d v S G ????????+==πρρ 7.某蒸汽管干管的始端蒸汽流速为25 m/s ，密度为2.62 kg/流体力学计算公式
流体力学公式总结
3 流体动力学基础
实际流体动力学基础
计算流体力学常用数值方法简介[1]
第三章 流体动力学基础
最新流体力学龙天渝课后答案第三章一元流体动力学基础
流体力学的计算公式
流体力学计算题..
流体力学复习要点(计算公式)
流体力学计算器 源代码
流体主要计算公式
流体力学复习要点(计算公式)
一元流体动力学基础
- 最新三章一元流体动力学基础
- 第三章 一元流体动力学基础(第一次).
- 第二章一元流体动力学
- 流体力学龙天渝课后答案第三章一元流体动力学基础(供参考)
- 第三章一元流体动力学基础ppt
- 第三章 一元流体动力学基础
- 一元流体动力学基础
- 一元流体动力学基础.
- 流体力学泵与风机-第3章 一元流体动力学基础
- 流体力学龙天渝课后答案第三章一元流体动力学基础
- 流体力学3 一元流体动力学基础
- 《流体力学》第三讲 一元流体动力学基本3.11-3.12
- 流体力学_03一元流体动力学基础讲解
- 一元流体动力学基础
- 流体力学龙天渝课后答案第三章一元流体动力学基础
- 第三章 一元流体动力学基础(第一次)解读
- 一元流体动力学基础
- 一元流体动力学
- 流体力学龙天渝课后答案第三章一元流体动力学基础
- 第三章 一元流体动力学基础