流体力学多相流自学作业

流体力学多相流自学作业
流体力学多相流自学作业

多相流及其应用

1.两相与多相流得定义与分类

在物理学中物质有固、液、气与等离子四态或四相。单相物质得流动称为单相流,两种混合均匀得气体或液体得流动也属于单相流。同时存在两种及两种以上相态得物质混合体得流动就就是两相或多相流。在多相流动力学中,所谓得相不仅按物质得状态,而且按化学组成、尺寸与形状等来区分,即不同得化学组成、不同尺寸与不同形状得物质都可能归属不同得相。在两相流研究中,把物质分为连续介质与离散介质。气体与液体属于连续介质,也称连续相或流体相。固体颗粒、液滴与气泡属于离散介质,也称分散相或颗粒相。流体相与颗粒相组成得流动叫做两相流动。

自然界与工业过程中常见得两相及多相流主要有如下几种,其中以两相流最为普遍。(1) 气液两相流

气体与液体物质混合在一起共同流动称为气液两相流。它又可以分单组分工质如水—水蒸气得汽液两相流与双组分工质如空气—水气液两相流两类,前者汽、液两相都具有相同得化学成分,后者则就是两相各具有不同得化学成分。单组分得汽液两相流在流动时根据压力与温度得变化会发生相变。双组分气液两相流则一般在流动中不会发生相变。

自然界中如下雨时得风雨交加,湖面与海面上带雾得上升气流、山区大气中得云遮雾罩。生活中沸腾得水壶中得循环,啤酒及汽水等夹带着气泡从瓶中注人杯子得流动等都属于气液两相流。现代工业设备中广泛应用着气液两相流与传热得原理与技术,如锅炉、核反应堆蒸汽发生器等汽化装置,石油、天然气得管道输送,大量传热传质与化学反应工程设备中得各种蒸发器、冷凝器、反应器、蒸馏塔、汽提塔,各式气液混合器、气液分离器与热交换器等,都广泛存在气液两相流与传热现象。

(2) 气固两相流

气体与固体颗粒混合在一起共同流动称为气固两相流。

空气中夹带灰粒与尘土、沙漠风沙、飞雪、冰雹,在动力、能源、冶金、建材、粮食加工与化工工业中广泛应用得气力输送、气流千燥、煤粉燃烧、石油得催化裂化、矿物得流态化焙烧、气力浮选、流态化等过程或技术,都就是气固两相流得具体实例。

严格地说,固体颗粒没有流动性,不能作流体处理。但当流体中存在大量固体小粒子流时,

如果流体得流动速度足够大,这些固体粒子得特性与普通流体相类似,即可以认为这些固体颗粒为拟流体,在适当得条件下当作流体流动来处理。在流体力学中,尽管流体分子间有间隙,但人们总就是把流体瞧着就是充满整个空间没有间隙得连续介质。由于两相流动研究得不就是单个颗粒得运动特性,而就是大量颗粒得统计平均特性。虽然颗粒得数密度(单位混合物体积中得颗粒数)比单位体积中流体分子数少得多,但当悬浮颗粒较多时,人们仍可设想离散分布于流体中颗粒就是充满整个空间而没有间隙得流体。这就就是常用得拟流体假设。(3) 液固两相流

液体与固体颗粒混合在一些共同流动称液固两相流。

自然界与工业中得典刑实例有夹带泥沙奔流得江河海水,动力、化工,采矿、建筑等工业工程中广泛使用得水力输送,矿浆、纸浆、泥浆、胶桨等浆液流动等。其它像火电厂锅炉得水力除渣管道中得水渣混合物流动,污水处理与排放中得污水管道流动等。

(4) 液液两相流

两种互不相溶得液体混合在一起得流动称液液两相流。

油田开采与地面集输、分离、排污中得油水两相流,化工过程中得乳浊液流动、物质提纯与萃取过程中大量得液液混合物流动均就是液液两相流得工程实例。

(5) 气液液、气液固与液液固多相流

气体、液体与固休颗粒混合在一起得流动称气液固三相流。气体与两种不能均匀混合、互不相溶得液体混合物在一起得共同流动称为气液液三相流;两种不能均匀混合、互不相溶得液体与固体颗粒混合在一起得共同流动称为液液固三相流。

图1 两相流得范围

2.两相流得研究与处理方法

与普通流体动力学类似,研究两相流问题得方法可以分为理论研究与实验研究两方面。从理论分析方法来瞧,仍然存在微观与宏观两种观点。微观分析法就就是从分子运动论出发,利用Boltzman(波尔兹曼)方程与统计平均概念及其理论,建立两相流中各相得基本守恒方程。宏观分析法,就就是以连续介质假设为基础,将两相流中各相都视为连续介质流体,根据每一相得质量、动量与能量宏观守恒方程以及相间相互作用,建立两相流得基本方程组,再利用这些两相流基本方程组去研究分析各种具体得两相流问题。在许多实际问题中,我们所关心得不就是单个物质粒子得运动而就是大量粒子运动所产生得总结果,也就就是所谓得宏观量,如压强、密度、温度、平均流速等。

从宏观观点分析两相流得方法又可以分为3类。

(1) 模型法

即假定相互扩散作用就是连续进行,其基本观点就是:

(1)两相流混合物体中得每一点都同时被两相所占据;

(2)混合物得热力学与输运特性取决于各相得特性与浓度;

(3)各相以自己得质量速度中心移动,相间相互扩散作用反映在模型内。

(2) 容积法

假定过程处于平衡状态,可用平衡方程式进行描述,基本方法就是:

(1)认为流动就是一维得;

(2)对一个有限容积写出质量、动量与能量守衡方程;

(3)守衡方程即可按混合物写出,也可按单独相列出。

(3) 平均法

假定过程处于平衡状态,用平均得守衡方程进行描述,类似低通滤波得方法。

上述3种方法得共同点就就是不考虑局部得与瞬时得特性,仅考虑相界面上流体微粒集中得相互作用,即宏观动力学。

3.两相与多相流得专用术语与基本持性参数

一般用下角标1与2分别表示两相流中得两种相或组分,对液—气两相系统用下角l与g 区分,对流体—固体两相系统用下角标f与s区分。通常选定组分或相2为分散相或为分层流

动中得轻相。

3、1 质量流量、质量流速与质量相含率(相分数)

质量流量就是指单位时间内流过通道总流通截面积得流体质量,用W表示。两相流总得质量流量就是各相质量流量之与,各相得质量流量分别用W1与W2表示,所以有

W= W1+W2(kg/s)

质量流速就是单位流通截面积上得质量流量,用G表示,如流道得总流通截面积为A,各相所占得流通截面积分别为A1与A2,则有

G=W/A G1=W1/A G2=W2/A

各相质量流量与总质量流里之比称为质量相含率或质量相分数。在气液两相流系统中,气体得质量相含率俗称为质量含气率或干度,用x表示;液体得质量相含率俗称为质量含液率。气液两相流得质量含气率与质量含液率之与就是1,所以有

x=G g/G=W g/W

1-x= G l/G=W l/W

3、2 容积流量、容积流速与容积相含率

容积流量就是指单位时间流过通道总流通截面积得流体容积,用Q表示。两相流总得容积流量同样就是各相容积流量Q1与Q2之与,即

Q= Q1+Q2=( W1/ρ1+W2/ρ2)

容积流速就是单位流通截面积上得容积流量,又称折算速度,就是容积流量除以通道总流通面积A,用J表示。

J=Q/A= (Q1+Q2)/A=J1+J2

J1= Q1/A= W1/(ρ1A)

J2= Q2/A= W2/(ρ2A)

各相得折算速度在两相流中就是十分重要并常用得一个术语与参数,它表示两相混合物中任何一相单独流过整个通道截面积时得速度,称为该相得折算速度。

容积相含率就是指各相容积流量与总容积流量之比。在气液两相流系统气相得容积相含率又称为容积含气率,用β表示,液相得容积相含率称为容积含液率。

3、3 各相真实流速

各相容积流量除以流动中各相各自所占流通截面积即为各相得真实流速。用v i表示(i= 1,2,或i= l,g,为气液两相流;i=f,s为流体—固相两相流)。

v i=Q i/A i

3、4 真实相含率或截面相含率

某相得流动在任意流通截面上所占通道截面积与总得流通截面积之比称作该相得真实相含率或截面相含率。对气液两相流,气相得真实相含率又称为截面含气率、真实含气率或空隙率,用α表示。而液相所占截面积与总流通截面积之比称为截面含液率,用(1-α}表示。即;

α=A g/A 1-α=A l/A

3、5 滑动比、滑移速度、漂移速度与漂移流率

两相流中各相真实速度得比值称为滑动比。气液两相流得滑动比用S表示,就是气相真实速度与液相真实速度之比

S=v g/v l

滑移速度就是指两相流各相真实速度得差,用v s表示

v s= v g-v l=J g/α-J l/(1-α)

漂移速度就是指轻相(如气相)速度与两相混合物平均速度v H之差,用v D表示

v D= v g-v H

两相混合物平均速度指当滑动比S=1时两相混合物得速度。

漂移流率就是指滑移速度v s得式两边乘以通分后得分母项,消去分母后得等式,用j D表示,有

j D=(v g-v l) α(1-α)= J g(1-α)-J lα

3、6 两相混合物得密度与比容

两相流体得密度有两种表示法:

(1)流动密度。指单位时间内流过截面得两相混合物得质量与容积之比,即

ρ0=W/Q

(2)真实密度。指流动瞬间任一流动截面上两相流混合物得密度,用ρm表示,定义如下:

ρm=αρg+(1-α)ρm

由截面含气率α与容积含气率β及滑移比S各自得定义可推得:

显然当S=1,即v g=v l,两相间无相对速度时,α=β流动密度才等于真实密度。

比容就是密度得倒数,因此,两相混合物得比容为:

流动比容;

真实比容

4.气液两相流得流型与流型图

4、1 垂直上升管中得气液两相流流型及其流型图

气液两相流在垂直管中上升流动时得几种常见流型。

(1)细泡状流型

细泡状流动就是最常见得流型之一。其特征为在液相中带有散布在液体中得细小气泡。直径在1mm以下得气泡就是球形得。直径在1mm以上得气泡外形就是多种多样得。

(2)气弹状流型

气弹状流型由一系列气弹组成。气弹端部呈球形而尾部就是平得。在两气弹之间夹有小气泡而气弹与管壁之间存在液膜。

(3)块块流型

当管内气速增大时,气弹发生分裂形成块状流型。此时大小不一得块状气体在液流中以混乱状态流动。

(4)带纤维得环状流型

在这种流型中,管壁上液膜较厚且含有小气泡。管子核心部分主要就是气体,但在气流中含有由被气体从液膜带走得细小液滴形成得长条纤维。

(5)环状流型

在这种流型中,管壁上有一层液膜,管子核心部分为带有自液膜卷入得细小液滴得气体。环状流型都发生在较高气体流速时。在受热管道中,当管壁温度高到足以使管壁液膜汽化时,流动结构就会发展到壁上无液膜,只有气相中还含有细小液滴得雾状流型。

图3表示有单组分气液两相流体在垂直上升受热管中得流型与管壁热流密度得关系。在图中,温度低于饱与温度得液体以固定流量进入各受热管。各受热管得热流密度依次自左向右逐渐增加。由图3可见,随着热流密度得增大,各管中得沸腾点逐渐移向管子进口,各管中得流型也逐渐由单相液体、细泡状流型、气弹状流型、块状流型、环状流型、雾状流型一直发展到干饱与蒸汽与过热蒸汽流动。

在气液两相流中,在两相流量、流体得物性值(密度、粘度、表面张力等)、管道得几何形状,管道尺寸以及热流密度确定得条件下,要判断管内气液两相流得流型可应用流型图。流型图主要就是根据试验资料总结而成得,因而应用流型图时不应超出获得该流型图得试验范围。

在判断垂直上升得流型图中,图4所示得流型图得到较为广泛得应用。此图适用于空气—水与汽—水两相流,就是在管子内径为31.2mm得管子中用压力为0、14-0、54 MPa

得空气—水混合物为工质得出得,此图与应用压力为3、45 -6、9 MPa得汽水混合物在管径为12.7mm管子中得出得试验数据相符良好,所以也可适用于上述参数得汽水混合物。

图中横坐标为脚ρL J2L纵坐标为ρG J2G可分别按下列两式计算:

4.2垂直下降管中得气液两相流流型及其流型图

在垂直管中气液两相一起往下流动时得流型示于图5。

这些流型就是由空气—水混合物得试验结果得出得。在气液两相做下降流动时得细泡状流型与上升流动时得细泡状流型不同。前者得细泡集中在管子核心部分,而后者则散布在整个管子截面上。如液相流量不变而使气相流量增大,则细泡将聚集成气弹。下降流动时得气弹状流型比上升流动时稳定。下降流动时得环状流动有几种流型,在气相及液相流量小时,有一层液膜沿管壁下流,核心部分为气相,这称为下降液膜流型。当液相流量增大,气相将进人液膜,这称为带气泡得下降液膜流型,当气液两相流量都增大时会出现块状流型。在气相流量较高时能发展为核心部分为雾状流动,壁面有液膜得雾式环状流型。

图6示出得就是作下降流动得气液两相流流型图。该图就是以空气与多种液体混合物作试验得出得,试验管径为25、 4 mm,试验压力为0、17 MPa。

图6选用作横坐标,用作纵坐标。,

Y为液相物性系数,

—液相动力粘度P a·s;—20℃,0、1MPa时水得动力粘度P a·s;

—液相密度kg/m3;—20℃,0、1MPa时水得密度kg/m3;

—液相表面张力N/m;—20℃,0、1MPa时水得表面张力N/m

4.3水平管中得气液两相流流型及其流型图

气液两相流体在水平管中流动时得流型种类比垂直管中得多。这主要就是由于重力得影响使气液两相有分开流动得倾向造成得。

气液两相流体在水平管中流动时流型大致可分为6种,即:细泡状流型、气塞状流型、分层流型、波状分层流型、气弹状流型及环状流型。图7中有这6种流型得示意图。

在细泡状流型中,由于重力得影响,细泡大都位于管子上部。当气体流量增加时,小气泡合并成气塞形成气塞状流型。分层流型发生于气液两相得流量均小时甲此时,气液两相分开流动,两相之间存在一平滑得分界面。当气相流量较高时,两相分界面上出现流动波,形成波状分层流型。气相流量再增大会形成气弹状流型,但此时气弹偏向管子上部。当气相流量很高而液相流量较低时就出现环状流型。

水平管中得气液两相流流型也可按相应得流型图确定。在现有得各种水平管流型图中,Baker流型图建立最早且得到广泛应用,特别在石油工业与冷凝工程设计中。图8所示为坐标经过改进后得Baker流型图。

横坐标为J LФ,纵坐标为J G/λ。其中J G、J L分别为气相及液相得折算质量流速;系数λ与Ф值与流体得物性有关。

5.界面现象

5.1概述

两相与多相流中另一个最基本最显著得特征就是混合物系统中存在明确得相界面,并且伴随着混合物得流动,这些相界面得状态、形状也处于不断得运动或变化发展中。相界面就是指将两种不同相得物质分隔开得区域,在此区域内得物质得特征与性质不同于相邻得区域。界面得邻近区域内物质得密度与浓度均有显著得差异,在接近两相分隔得临界点,密度在与界面相垂直得方向上就是位置得连续函数。也就就是说界面就是一个三维区域,其厚度可能有几个分子直径大小或者更厚一些,这已经为统计力学得计算结果所证实。

所谓界面现象,就是指与相界面上得动量、能量与质量传递相关得所有效应。许多物理、化学作用都发生在表面与界面上,例如沸腾与凝结得状态会因界面上存在即便就是极少量得杂质而发生明显得改变,泡沫得稳定性、湖面得波度以及镜体表面得汽水雾罩等现象都就是发生在界面上得物理、化学变化。

5.2速度边界条件

如果不考虑相变或传质、速度相容性在单组分与双组分得两相流动中都就是相同得。但在有相变时,就可能在界面上存在一有限得速度跳跃。流体流经界面必须满足连续性条件,如果密度有变化则速度也会有变化。囚此,相对于图9所示得相界面,在平行与垂直相界面得两个方向上都必须满足如下得相容性条件:

(1)连续切向速度要求

v1t=v2t

(2)穿透相界面得连续性要求

ρ1v1n=ρ2v2n =m

式中:m表示垂直通过相界面得流体质量流量。通过相变过程中得平均连续速度被假设为垂直于相界面得。如果界面还正在移动、就将界面得移动速度简单地叠加到上述速度上去。

5.3应力边界条件

表面张力就是界面上每单位面积得自由能,即形成一单位表面积所需外界做得功,因为将物质分子从体相内移动到表面区就是要克服分子间得吸力而做功得。也可以把表面张力想象为表面作用在每单位长度截线上得力。界面稳定存在得必要条件就是表面张力大于0,如果小于0,那么偶然得涨落就可导致表面区不断扩大,最后使一种物质完全分散到另一种物质中,两种气体之间混合得情况就就是这样。

通过界面得应力场得连续性因表面张力效应而改变。如果表面张力就是均匀得,则通过界面得剪切应力也就是连续得。但杂质、灰尘或表面活性剂得存在、或则沿界面得温度梯度等都会引起表面张力梯度得形成。这种情况下,通过界面得剪切应力存在一等于表面张力梯度得突增量。

6.两相流基本数理模型

两相流中各相在空间与时间上随机扩散,同时存在动态得相互作用。对于这种复杂得二维两相瞬态问题,完整得解析解目前还无法导出。人们先后提出了多种数理模型,从最简单得均相流模型一直到最复杂得双流体模型。

6、1均相流动模型(Homogeneous Flow Model)

均相流动模型把气液两相混合物瞧作就是一种均匀介质,相间没有相对速度,流动参数取两相相应参数得平均值。在此基础土,可将两相流视为具有平均流体特性得单相流对待。在均相模型中采用了两个基本假设:

(1)两相间处于热力学平衡状态,即两相具有相同得温度并且都处于饱与状态;

(2)气液两相得流速相等,即为均匀流。其数学表达式及守恒方程如下:

u G=u L=u H

S= u G/u L=1

在通道微元段上得动量守恒方程为:

动量增加率=动量流出率-动量流入率+动量存积率=作用于控制体得力之与

Pp为通道圆周

对于等截面通道得稳定流动,这个方程简化为

能量方程并不直接用于压降得计算,主要用来计算绝热流动中得局部干度,能量守衡方程为:

能量流人率=能量流出率+能量存积率

化简后e—单位质量流体得对流传热热能

实际上气液两相得流速并不相等,只有在高含气量或很小含气量时两相流速才近似相等,因此这一模型实际上只适用于泡状流与雾状流。从前人得研究中可以瞧出,均相模型得计算结果同实验值有较大出人,偏差随质量流速得减小而增大。这种偏差得存在就是因为均相模型假设两相之间没有速度差异。当质量流速较小时,浮力效应显著,引起两相速度间相当大得差异;而质量流速较大时,液相湍动得结果使得两相得混合更加均匀,因此质量流速

增大时偏差减小。

6、2分相流动模型( Separated Flow M0del )

分相流动模型就是将气、液两相都当作连续流体分别来处理,并考虑了两相之间得相互

作用。其基本假设就是:

(1)两相间保持热力学平衡;

(2)气液两相得速度为常量,但不一定相等。假定气液两相都以一定得平均速度在流道中流动。分相流动模型在一定程度上考虑了两相间得相互作用,计算结果比均相模型理想。当两相平均流动速度相等时,分相流动即可转化为均相模型。因此可将均相模型视为分相模型得一个特殊情况。分相流动模型适用于两相间存在微弱耦合得场合,如分层流与环状流。分相流动得基本方程推导如下:

先写出每一相得各个守恒方程,然后将每一组得守恒方程进行合并,成为混合物得总平衡方程这也就是在确定压降时被人们广泛采用得分析方法。

连续性方程,也就就是质量守恒方程得数学式为

()()()()()()01111=-??+--+-??+-z A t

A u z G A u z z A u L L L e L L L L δαραρδαρδαρ Ge 就是单位长度上液体转化为气体得转化率。

分相流动模型得动量守衡方程同样可以从各相得动量方程导出,混合物动量方程

()()()()

22111sin G G L L G G L L p W tp u A u z A u u t A p g z p αραραρρτθρ+-??++?-??=--??- 分相流动模型得两相能量守恒方程式为:

()()()()()()()

()t p A x x G t gGA x x A q z A q qP xr h GA z h h t A G L G L v p L G H L L ??+???? ?????? ??+--??--???

? ?????? ??+--??-+=+??++-??αραρθαραραραυ222222223311sin 1116、3漂移模型(Drift-Flux Model)

这一模型主要由Zuber 等人提出。它就是在热力学平衡得假设下,建立在两相平均速度场基础上得一种模型。漂移模型提出了一个漂移速度得概念,当两相流以某一混合速度流动时,气相相对于这个混合速度有一个漂移速度,液体则有一个反向得漂移速度以保持流动得连续性。在守恒方程组中将相间相对速度以漂移速度来考虑,通过附加得气相连续方程来描写气液两相流动。

6、4 两流体模型(Two Fluid Model)

由于各相得动力学性质不完全相同或浓度分布不均匀,气、液两相得运动存在相当大得差异。上述几种模型中,均相模型完全没有考虑两相得差异。分相模型与漂移模型在一定程

度上引入了气液两相得相互作用,但仍过于简单而无法精确描述两相得运动与空间分布。

目前一致公认最为完善可靠得模型就是两流体模型。它可以用欧拉方法、拉格朗日方法来描述,每种方法都有其固有得优缺点。工程实际中采用得方法一般基于积分法或时均紊流模型,往往不能提供瞬态得流动结构。而这对于离散相得输运恰恰又就是非常关键得。只有拉格朗日方法能够考虑离散相与瞬态流动结构之间复杂得相互作用,但其表达式尚存在某种程度得不确定性,需要很多理想化得假设。

两流体模型将每一种流体都瞧作就是充满整个流场得连续介质,针对两相分别写出质量、动量与能量守恒方程,通过相界面间得相互作用(动量、能量与质量得交换)将两组方程祸合在一起。这种方法只需假设每相在局部范围内都就是连续介质,不必引人其它人为假设,而且对两相流得种类与流型没有任何限制,适用于可当做连续介质研究得任何二元混合物,所建立得两相流方程就是目前最全面完整得,求得得解中包含得信息丰富完全。但两流体模型包含得变量多、方程复杂,因此求解比较困难。

7.多相流应用

7、1 高效旋风式分离器

其办法就是把圆筒壁面做成粗糙面,使湍流增强,碰撞壁面机会增多。旋风分离器可用于浮渣分离、豆奶脱臭、粉尘分离、喷漆分离等有强刺激气体得场合。

7、2高浓度多相流泵

常用得泵多按单相流理论设计。对于高浓度多相流,由于流体得密度、粘度都发生了很大得变化,且随着各相浓度得不同,固相或气相均有不同变化即使旋转速度不变化。即由于粘度变化,雷诺数相似律不存在了因此。高粘流体泵得设计应按流体得流变参数考虑。

我国蔡保元教授设汁得两相流泵只考虑固、液两相得相对运动,把固相作为动得边界条

件,山于两者得惯性力不同,便会发生“相对阻塞”与“相对抽吸”作用。在抽浑水时、两相流泵得效率比清水泵高但该两相流泵没考虑粘度得影响,而在高粘度下这种影响就是不可忽视得。

7、3 三相流垂直提升一泵

为了开采海底矿产,国外已应用三相流垂直提升技术采集锰结核,可把5000米深得海底得锰块提升上来。这一技术也可用于深井开采、井底清淤、河道疏浚、污染底泥处理等方面。

参考文献

[1]郭烈锦、两相与多相流动力学[M]、西安:西安交通大学出版社、

工程流体力学习题全解

第1章 绪论 选择题 【1.1】 按连续介质的概念,流体质点是指:(a )流体的分子;(b )流体内的固体颗粒; (c )几何的点;(d )几何尺寸同流动空间相比是极小量,又含有大量分子的微元体。 解:流体质点是指体积小到可以看作一个几何点,但它又含有大量的分子,且具有 诸如速度、密度及压强等物理量的流体微团。 (d ) 【1.2】 与牛顿内摩擦定律直接相关的因素是:(a )切应力和压强;(b )切应力和剪切变 形速度;(c )切应力和剪切变形;(d )切应力和流速。 解:牛顿内摩擦定律是 d d v y τμ =,而且速度梯度d d v y 是流体微团的剪切变形速度 d d t γ,故d d t γ τμ=。 (b ) 【1.3】 流体运动黏度υ的国际单位是:(a )m 2/s ;(b )N/m 2;(c )kg/m ;(d )N·s/m 2。 解:流体的运动黏度υ的国际单位是/s m 2 。 (a ) 【1.4】 理想流体的特征是:(a )黏度是常数;(b )不可压缩;(c )无黏性;(d )符合RT p =ρ 。 解:不考虑黏性的流体称为理想流体。 (c ) 【1.5】 当水的压强增加一个 大气压时,水的密度增大约为:(a )1/20 000;(b )1/1 000;(c )1/4 000;(d )1/2 000。 解:当水的压强增加一个大气压时,其密度增大约 95d 1 d 0.51011020 000k p ρ ρ -==???= 。 (a ) 【1.6】 从力学的角度分析,一般流体和固体的区别在于流体:(a )能承受拉力,平衡时 不能承受切应力;(b )不能承受拉力,平衡时能承受切应力;(c )不能承受拉力,平衡时不能承受切应力;(d )能承受拉力,平衡时也能承受切应力。 解:流体的特性是既不能承受拉力,同时具有很大的流动性,即平衡时不能承受切应力。 (c ) 【1.7】 下列流体哪个属牛顿 流体:(a )汽油;(b )纸浆;(c )血液;(d )沥青。 解:满足牛顿内摩擦定律的流体称为牛顿流体。 (a ) 【1.8】 15C o 时空气和水的运动黏度6215.210m /s υ-=?空气,62 1.14610m /s υ-=?水 ,这 说明:在运动中(a )空气比水的黏性力大;(b )空气比水的黏性力小;(c )空气 与水的黏性力接近;(d )不能直接比较。 解:空气的运动黏度比水大近10倍,但由于水的密度是空气的近800倍,因此水的黏度反而比空气大近50倍,而黏性力除了同流体的黏度有关,还和速度梯度有 关,因此它们不能直接比较。 (d ) 【1.9】 液体的黏性主要来自于液体:(a )分子热运动;(b )分子间内聚力;(c )易变形

流体力学大作业

《计算流体力学》课程大作业 作业内容: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中的实验结果进行对比。

流体力学课后习题第四章作业答案

第四章作业答案 4-3水在变直径竖管中流动,已知粗管直径 d 1=300mm ,流速v 1=6m/s 。两断面相距3m,为使两断面的压力表读值相同。试求细管直径(水头损失不计)。 解: 221122122222 112222p v p v Z Z g 2g g 2g p v p v v 6 300 3 4.837m v 9.74m/s g 2g g 2g 2g 2g l h ρρρρ++=+++++=+++=+=?= 2 2 2211 21v d v d d 300235.5mm ==== 4—4变直径管段AB ,d A =0.2m,d B =0.4m ,高差△h=1.5m,测得p A =30kPa ,p B =40kPa ,B 点处断面平均流速v B =1.5m/s ,试判断水在管中的流动方向。 解: 222 2222 0.43061.5()6m/s 0 4.900.229.8240 1.51.5 5.69m 29.819.6 B A A A B A A A B B B B d p H z m d g g g p H Z g g υυυρυρ==?==++=++==++=++= H B >H A , 水由B 流向A; 水头损失5.69-4.90=0.79m 4—5用水银压差计测量水管中的点流速u ,如读值 △h=60mm ,(1)求该点流速;(2)若管中流体是30.8/kg m ρ=的油,△h 不变,不计水头损失,则该点的流速是多少? 解: (1) 3.85m/s u === (2) 4.34m/s u === 4—6 利用文丘里管的喉管处负压抽吸基坑中的积水,已经知道管道直径1100d mm =, 喉管直径2 50d mm =,2h m =,能量损失忽略不计。试求管道中流量至少为多大, 才能抽出基坑中的积水? 解:由题意知,只有当12 12()()p p z z h g g ρρ+-+=时,刚好才能把水吸上来,由文丘里流 量计原理有Q =, 其中211 d k π=, 代入数据,有12.7Q l s =。

工程流体力学思考题章模板

第一章绪论 1、什么叫流体? 流体与固体的区别? 流体是指能够流动的物质, 包括气体和液体。 与固体相比, 流体分子间引力较小, 分子运动剧烈, 分子排列松散, 这就决定了流体不能保持一定的形状, 具有较大流动性。 2、流体中气体和液体的主要区别有哪些? (1)气体有很大的压缩性, 而液体的压缩性非常小; (2)容器内的气体将充满整个容器, 而液体则有可能存在自由液面。 3、什么是连续介质假设? 引入的意义是什么? 流体充满着一个空间时是不留任何空隙的, 即把流体看作是自由介质。 意义: 不必研究大量分子的瞬间运动状态, 而只要描述流体宏观状态物理量, 如密度、质量等。 4、何谓流体的压缩性和膨胀性? 如何度量? 压缩性: 温度不变的条件下, 流体体积随压力变化而变化的性质。用体积压缩系数βp表示, 单位Pa-1。 膨胀性: 压力不变的条件下, 流体体积随温度变化而变化的性质。用体积膨胀系数βt表示, 单位K-1。 5、何谓流体的粘性, 如何度量粘性大小, 与温度关系? 流体所具有的阻碍流体流动, 即阻碍流体质点间相对运动的性质称为粘滞性, 简称粘性。用粘度μ来表示, 单位N·S/m2或Pa·S。

液体粘度随温度的升高而减小, 气体粘度随温度升高而增大。 6、作用在流体上的力怎样分类, 如何表示? (1) 质量力: 采用单位流体质量所受到的质量力f 表示; (2) 表面力: 常见单位面积上的表面力Pn 表示, 单位Pa 。 7、什么情况下粘性应力为零? ( 1) 静止流体 ( 2) 理想流体 第二章 流体静力学 1、 流体静压力有哪些特性? 怎样证明? ( 1) 静压力沿作用面内法线方向, 即垂直指向作用面。 证明: ○1流体静止时只有法向力没有切向力, 静压力只能沿法线方向; ○2流体不能承受拉力, 只能承受压力; 因此, 静压力唯一可能的方向就是内法线方向。 ( 2) 静止流体中任何一点上各个方向静压力大小相等, 与作用方向无关。 证明: 2、 静力学基本方程式的意义和使用范围? 静力学基本方程式:Z+ g P ρ=C 或 Z 1+g P ρ1=Z 2+g P ρ2 (1) 几何意义: 静止流体中测压管水头为常数 物理意义: 静止流体中总比能为常数

流体力学 大作业

一.选择题 1.牛顿内摩擦定律适用于()。 A.任何流体B.牛顿流体C.非牛顿流体 2.液体不具有的性质是()。 A.易流动性B.压缩性C.抗拉性D.粘滞性 3连续介质假定认为流体()连续。 A.在宏观上B.在微观上C.分子间D.原子间 4.在国际单位制中流体力学基本量纲不包括()。 A.时间B.质量C.长度D.力. 5.在静水中取一六面体,作用在该六面体上的力有() A.切向力、正压力B.正压力C.正压力、重力D.正压力、切向力、重力 6.下述哪些力属于质量力( ) A.惯性力B.粘性力C.弹性力D.表面张力E.重力 7.某点存在真空时,()() A.该点的绝对压强为正值B.该点的相对压强为正值c.该点的绝对压强为负值D.该点的相对压强为负值 8.流体静压强的()。 A.方向与受压面有关B.大小与受压面积有关B.大小与受压面方位无关 9.流体静压强的全微分式为()。 A.B.C. 10.压强单位为时,采用了哪种表示法()。 A.应力单位B.大气压倍数C.液柱高度 11.密封容器内液面压强小于大气压强,其任一点的测压管液面()。A.高于容器内液面B.低于容器内液面C.等于容器内液面 12.流体运动的连续性方程是根据()原理导出的。 A.动量守恒 B. 质量守恒 C.能量守恒 D. 力的平衡 13. 流线和迹线重合的条件为()。

A.恒定流 B.非恒定流 C.非恒定均匀流 14.总流伯努利方程适用于()。 A.恒定流 B.非恒定流 C.可压缩流体 15. 总水头线与测压管水头线的基本规律是:()、() A.总水头线总是沿程下降的。 B.总水头线总是在测压管水头线的上方。 C.测压管水头线沿程可升可降。 D.测压管水头线总是沿程下降的。 16 管道中液体的雷诺数与()无关。 A. 温度 B. 管径 C. 流速 D. 管长 17.. 某圆管直径d=30mm,其中液体平均流速为20cm/s。液体粘滞系数为0.0114cm3/s,则此管中液体流态为()。 A. 层流 B. 层流向紊流过渡 C.紊流 18.等直径圆管中紊流的过流断面流速分布是()A呈抛物线分布B. 呈对数线分布 C.呈椭圆曲线分布 D. 呈双曲线分布 19.等直径圆管中的层流,其过流断面平均流速是圆管中最大流速的() A 1.0倍B.1/3倍C. 1/4倍D. 1/2倍 20.圆管中的层流的沿程损失与管中平均流速的()成正比. A. 一次方 B. 二次方 C. 三次方 D. 四次方 21..圆管的水力半径是( ) A. d/2 B. d/3 C. d/4 D. d/5. 22谢才公式中谢才系数的单位是()A. 无量纲B. C. D. . 23. 判断层流和紊流的临界雷诺数是() A.上临界雷诺数 B.下临界雷诺数 C.上下临界雷诺数代数平均 D.上下临界雷诺数几何平均 24.. 对于管道无压流,当充满度分别为()时,其流量和速度分别达到最大。 A. 0.5, 0.5 B. 0.95, 0.81 C. 0.81, 081 D. 1.0, 1.0 25.对于a, b, c三种水面线,下列哪些说法是错误()() A.所有a、c型曲线都是壅水曲线,即,水深沿程增大。B.所有b型曲线都是壅水曲线,即,水深沿程增大。C.所有a、c型曲线都是降水曲线,即,水深沿程减小。C.所有b型曲线都是降水曲线,即,水深沿程减

流体力学作业题参考题答案

流体力学网上作业题参考答案 第一章:绪论(56) 一、名词解释 1. 流体:在任何微小剪切力的持续作用下能够连续不断变形的物质/液体和气体统称为流体。 2. 流体质点:质点亦称为流体微团,其尺度在微观上足够大,大到能包含大量的分子,使得在统计平均 后能得到其物理量的确定值;而宏观行又足够小,远小于所研究问题的特征尺度,使其平均物理量可看成是连续的。 3. 惯性:惯性是物体维持原有运动状态的能力的性质。 4. 均质流体:任意点上密度相同的流体,称为均质流体。 5. 非均质流体:各点密度不完全相同的流体,称为非均质流体。 6. 粘性:流体在运动状态下具有抵抗剪切变形能力的性质,成为粘性或者粘滞性。 7. 内摩擦力:流体在流动时,对相邻两层流体间发生的相对运动,会产生阻碍其相对运动的力,这种力称为内摩擦力。 8. 流体的压缩性:流体的宏观体积随着作用压强的增大而减小的性质,称为流体的压缩性。 9. 不可压缩流体:流体密度随压强变化很小,流体的密度可视为常数的流体。 10. 可压缩流体:流体密度随压强变化不能忽略的流体。 11. 表面张力:表面张力是液体自由表面在分子作用半径范围内,由于分子引力大于斥力而在表层沿表面方向产生的拉力。 12. 表面力:作用在所研究流体外表面上,与表面积大小成正比的力。 13.质量力:作用在液体每一个质点上,其大小与液体质量成正比。 14.牛顿流体:把在作剪切运动时满足牛顿内摩擦定律的流体称为牛顿流体。 二、选择题 1. 理想流体是指可以忽略( C )的流体。A 密度 B 密度变化 C 粘度 D 粘度变化 2. 不可压缩流体是指忽略( B )的流体。A 密度 B 密度变化 C 粘度 D 粘度变化 3. 关于流体粘性的说法中,不正确的是( D )。 A 粘性是流体的固有属性 B 流体的粘性具有传递运动和阻碍运动的双重性 C 粘性是运动流体产生机械能损失的根源 D 流体的粘性随着温度的增加而降低 4.下列各种流体中,属于牛顿流体的是( A ) A 水、汽油、酒精 B 水、新拌混凝土、新拌建筑砂浆 C 泥石流、泥浆、血浆 D 水、泥石流、天然气 5. 单位质量力的量纲与( B )的量纲相同。A 速度 B 加速度 C 粘性力 D 重力 6、在水力学中,单位质量力是指( C ) a、单位面积液体受到的质量力; b、单位体积液体受到的质量力; c、单位质量液体受到的质量力; d、单位重量液体受到的质量力。 7、运动粘度与动力粘度的关系是( A ) A μ υ ρ = B g μ υ= C υ μ ρ = D g υ μ=

工程流体力学_思考题__1~4章

向期末进发!!! 第一章绪论 1、什么叫流体?流体与固体的区别? 流体是指可以流动的物质,包括气体和液体。 与固体相比,流体分子间引力较小,分子运动剧烈,分子排列松散,这就决定了流体不能保持一定的形状,具有较大流动性。 2、流体中气体和液体的主要区别有哪些? (1)气体有很大的压缩性,而液体的压缩性非常小; (2)容器内的气体将充满整个容器,而液体则有可能存在自由液面。 3、什么是连续介质假设?引入的意义是什么? 流体充满着一个空间时是不留任何空隙的,即把流体看作是自由介质。 意义:不必研究大量分子的瞬间运动状态,而只要描述流体宏观状态物理量,如密度、质量等。 4、何谓流体的压缩性和膨胀性?如何度量? 压缩性:温度不变的条件下,流体体积随压力变化而变化的性质。用体积压缩系数βp表示,单位Pa-1。 膨胀性:压力不变的条件下,流体体积随温度变化而变化的性质。用体积膨胀系数βt表示,单位K-1。 5、何谓流体的粘性,如何度量粘性大小,与温度关系? 流体所具有的阻碍流体流动,即阻碍流体质点间相对运动的性质称为粘滞性,简称粘性。用粘度μ来表示,单位N·S/m2或Pa·S。 液体粘度随温度的升高而减小,气体粘度随温度升高而增大。 6、作用在流体上的力怎样分类,如何表示? (1)质量力:采用单位流体质量所受到的质量力f表示; (2)表面力:常用单位面积上的表面力Pn表示,单位Pa。 7、什么情况下粘性应力为零? (1)静止流体(2)理想流体 第二章流体静力学 1、流体静压力有哪些特性?怎样证明? (1)静压力沿作用面内法线方向,即垂直指向作用面。 证明:○1流体静止时只有法向力没有切向力,静压力只能沿法线方向; ○2流体不能承受拉力,只能承受压力; 所以,静压力唯一可能的方向就是内法线方向。 (2)静止流体中任何一点上各个方向静压力大小相等,与作用方向无关。 证明:

第二章计算流体力学的基本知识

第二章计算流体力学的基本知识 流体流动现象大量存在于自然界及多种工程领域中,所有这些工程都受质量守恒、动量守恒和能量守恒等基本物理定律的支配。这章将首先介绍流体动力学的发展和流体力学中几个重要守恒定律及其数学表达式,最后介绍几种常用的商业软件。 2.1计算流体力学简介 2.1.1计算流体力学的发展 流体力学的基本方程组非常复杂,在考虑粘性作用时更是如此,如果不靠计算机,就只能对比较简单的情形或简化后的欧拉方程或N-S方程进行计算。20 世纪30~40 年代,对于复杂而又特别重要的流体力学问题,曾组织过人力用几个月甚至几年的时间做数值计算,比如圆锥做超声速飞行时周围的无粘流场就从1943 年一直算到1947 年。 数学的发展,计算机的不断进步,以及流体力学各种计算方法的发明,使许多原来无法用理论分析求解的复杂流体力学问题有了求得数值解的可能性,这又促进了流体力学计算方法的发展,并形成了"计算流体力学" 。 从20 世纪60 年代起,在飞行器和其他涉及流体运动的课题中,经常采用电子计算机做数值模拟,这可以和物理实验相辅相成。数值模拟和实验模拟相互配合,使科学技术的研究和工程设计的速度加快,并节省开支。数值计算方法最近发展很快,其重要性与日俱增。 自然界存在着大量复杂的流动现象,随着人类认识的深入,人们开始利用流动规律来改造自然界。最典型的例子是人类利用空气对运动中的机翼产生升力的机理发明了飞机。航空技术的发展强烈推动了流体力学的迅速发展。 流体运动的规律由一组控制方程描述。计算机没有发明前,流体力学家们在对方程经过大量简化后能够得到一些线形问题解读解。但实际的流动问题大都是复杂的强非线形问题,无法求得精确的解读解。计算机的出现以及计算技术的迅速发展使人们直接求解控制方程组的梦想逐步得到实现,从而催生了计算流体力

计算流体力学大作业

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 f t x ??+=?? (0.1) 分量方程: 连续性方程、动量方程和能量方程分别是: 2 22,,p u ρ

() ()()()2 000u t x u u p t x x u E p E t x ρρρρ???+ =?????????++=? ??????+?????+ =????? (0.2) 其中 22v u E c T ρ?? =+ ?? ? 对于完全气体,在量纲为一的形式下,状态方程为: ()2 p T Ma ργ∞ = (0.3) 在量纲为一的定义下,定容热容v c 为: () 21 1v 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 λλλ==+=-,依据如下算法将其分裂成正负特征值: () 12 222 k k k λλελ±±+= (0.6) 3.2 流通矢量的分裂 这里对流通矢量的分裂选用Steger-Warming 分裂法,分裂后的流通矢量为 ()()()()()()()12312322232121212122f u u c u c u u c u c w γλλλργλλλγλλγλ?? ? -++ ?=-+-++ ? ? ? -+-+++ ??? +++++++ ++ ++ (0.7)

流体力学作业答案

流体力学作业答案-标准化文件发布号:(9456-EUATWK-MWUB-WUNN-INNUL-DDQTY-KII

例题1:如下图所示,一圆锥体绕自身轴线等速旋转,锥体与固定壁面间的距离为K ,空隙全部被动力粘滞系数为 μ的牛顿流体所充满。当旋转角速度为ω,锥体底部半径 为R ,高为H ,求作用于圆锥的阻力矩。 解:M= ????====K dh r K dA r K dAr K u dAr 3 22cos 2πμω απμωωμμτ= H K R α πμωcos 23 而 22cos R H H +=α;故:M=223 2R H K R +?πμω 例题2:涵洞进口处,装有与水平线成600倾角而边长为1m 的正方形平板闸门(AB=1m ),求闸门所受静水总压力的大小及作用点。

解:坐标只能建在水面上。 a A kp p 807.91807.9=?= a B kp p 300.18)2 3 1(807.9=+?=

KN p p P B A 050.14112 =??+= h h A y I y y C C C C C D 6.160sin 433.1121 60sin 433.1160 sin 121160sin 00 30=+=??+=+= 0=D x 矩形和圆形的 C y 和C I 值 矩形: 2h y C = 123 bh I C =

圆形: r y C = 4 4 r I C π= 例题3:一直立矩形闸门,用三根工字梁支撑,门高及上游水深H 均为3m,把此闸门所受静水压强分布图分为三等份,每根工字梁分别设在这三等份的重心,求三个工字梁的位置? 解:设静水压力分布图的面积为A ,则每一等份为A/3

工程流体力学思考题章样本模板

向期末进发! ! ! 第一章绪论 1、什么叫流体? 流体与固体的区别? 流体是指能够流动的物质, 包括气体和液体。 与固体相比, 流体分子间引力较小, 分子运动剧烈, 分子排列松散, 这就决定了流体不能保持一定的形状, 具有较大流动性。 2、流体中气体和液体的主要区别有哪些? (1)气体有很大的压缩性, 而液体的压缩性非常小; (2)容器内的气体将充满整个容器, 而液体则有可能存在自由液面。 3、什么是连续介质假设? 引入的意义是什么? 流体充满着一个空间时是不留任何空隙的, 即把流体看作是自由介质。 意义: 不必研究大量分子的瞬间运动状态, 而只要描述流体宏观状态物理量, 如密度、质量等。 4、何谓流体的压缩性和膨胀性? 如何度量? 压缩性: 温度不变的条件下, 流体体积随压力变化而变化的性质。用体积压缩系数βp表示, 单位Pa-1。 膨胀性: 压力不变的条件下, 流体体积随温度变化而变化的性质。用体积膨胀系数βt表示, 单位K-1。 5、何谓流体的粘性, 如何度量粘性大小, 与温度关系? 流体所具有的阻碍流体流动, 即阻碍流体质点间相对运动的性

质称为粘滞性, 简称粘性。用粘度μ来表示, 单位N ·S/m 2或Pa ·S 。 液体粘度随温度的升高而减小, 气体粘度随温度升高而增大。 6、作用在流体上的力怎样分类, 如何表示? (1) 质量力: 采用单位流体质量所受到的质量力 f 表示; (2) 表面力: 常见单位面积上的表面力 Pn 表示, 单位Pa 。 7、什么情况下粘性应力为零? ( 1) 静止流体 ( 2) 理想流体 第二章 流体静力学 1、 流体静压力有哪些特性? 怎样证明? ( 1) 静压力沿作用面内法线方向, 即垂直指向作用面。 证明: ○1流体静止时只有法向力没有切向力, 静压力只能沿法线方向; ○2流体不能承受拉力, 只能承受压力; 因此, 静压力唯一可能的方向就是内法线方向。 ( 2) 静止流体中任何一点上各个方向静压力大小相等, 与作用方向无关。 证明: 2、 静力学基本方程式的意义和使用范围? 静力学基本方程式:Z+g P ρ=C 或 Z 1+g P ρ1=Z 2+g P ρ2 (1) 几何意义: 静止流体中测压管水头为常数

计算流体力学过渡到编程的傻瓜入门教程

借宝地写几个小短文,介绍CFD的一些实际的入门知识。主要是因为这里支持Latex,写起来比较方便。 CFD,计算流体力学,是一个挺难的学科,涉及流体力学、数值分析和计算机算法,还有计算机图形学的一些知识。尤其是有关偏微分方程数值分析的东西,不是那么容易入门。大多数图书,片中数学原理而不重实际动手,因为作者都把读者当做已经掌握基础知识的科班学生了。所以数学基础不那么好的读者往往看得很吃力,看了还不知道怎么实现。本人当年虽说是学航天工程的,但是那时本科教育已经退步,基础的流体力学课被砍得只剩下一维气体动力学了,因此自学CFD的时候也是头晕眼花。不知道怎么实现,也很难找到教学代码——那时候网络还不发达,只在教研室的故纸堆里搜罗到一些完全没有注释,编程风格也不好的冗长代码,硬着头皮分析。后来网上淘到一些代码研读,结合书籍论文才慢慢入门。可以说中间没有老师教,后来赌博士为了混学分上过CFD专门课程,不过那时候我已经都掌握课堂上那些了。 回想自己入门艰辛,不免有一个想法——写点通俗易懂的CFD入门短文给师弟师妹们。本人不打算搞得很系统,而是希望能结合实际,阐明一些最基本的概念和手段,其中一些复杂的道理只是点到为止。目前也没有具体的计划,想到哪里写到哪里,因此可能会很零散。但是我争取让初学CFD 的人能够了解一些基本的东西,看过之后,会知道一个CFD代码怎么炼成的(这“炼”字好像很流行啊)。欢迎大家提出意见,这样我尽可能的可以追加一些修改和解释。

言归正传,第一部分,我打算介绍一个最基本的算例,一维激波管问题。说白了就是一根两端封闭的管子,中间有个隔板,隔板左边和右边的气体状态(密度、速度、压力)不一样,突然把隔板抽去,管子内面的气体怎么运动。这是个一维问题,被称作黎曼间断问题,好像是黎曼最初研究双曲微分方程的时候提出的一个问题,用一维无粘可压缩Euler方程就可以描述了。 这里 这个方程就是描述的气体密度、动量和能量随时间的变化()与它们各自的流量(密度流量,动量流量,能量流量 )随空间变化()的关系。 在CFD中通常把这个方程写成矢量形式 这里 进一步可以写成散度形式

计算流体力学课程大作业

《计算流体力学》课程大作业 ——基于涡量-流函数法的不可压缩方腔驱动流问题数值模拟 张伊哲 航博101 1、 引言和综述 2、 问题的提出,怎样使用涡量-流函数方法建立差分格式 3、 程序说明 4、 计算结果和讨论 5、 结论 1引言 虽然不可压缩流动的控制方程从形式上看更为简单,但实际上,目前不可压缩流动的数值方法远远不如可压缩流动的数值方法成熟。 考虑不可压缩流动的N-S 方程: 01()P t νρ??=? ? ??+??=-?+???? U U UU f U (1.1) 其中ν是运动粘性系数,认为是常数。将方程组写成无量纲的形式: 01()Re P t ??=?? ??+??=-?+????U U UU f U (1.2) 其中Re 是雷诺数。 从数学角度看,不可压缩流动的控制方程中不含有密度对时间的偏导数项,方程表现出椭圆-抛物组合型的特点;从物理意义上看,在不可压缩流动中,压力这一物理量的波动具有无穷大的传播速度,它瞬间传遍全场,以使不可压缩条件在任何时间、任何位置满足,这就是椭圆型方程的物理意义。这就造成不可压缩的N-S 方程不能使用比较成熟的发展型...偏微分方程的数值求解理论和方法。 如果将动量方程和连续性方程完全耦合求解,即使使用显示的离散格式,也将会得到一个刚性很强的、庞大的稀疏线性方程组,计算量巨大,更重要的问题是不易收敛。因此,实际应用中,通常都必须将连续方程和动量方程在一定程度上解耦。 目前,求解不可压缩流动的方法主要有涡量-流函数法,SIMPLE 法及其衍生的改进方法,有限元法,谱方法等,这些方法各有优缺点。其中涡量-流函数法是解决二维不可压缩流动的有效方法。作者本学期学习了研究生计算流体课程,为了熟悉计算流体的基本方法,选择使用涡量-流函数法计算不可压缩方腔驱动流问题,并且对于不同雷诺数下的解进行比较和分析,得出一些结论。 本文接下来的内容安排为:第2节提出不可压缩方腔驱动流问题,并分析该问题怎样使用涡量-流函数方法建立差分格式、选择边界条件。第3节介绍程序的结构。第4节对于不同雷诺数下的计算结果进行分析,并且与U.GHIA 等人【1】的经典结论进行对比,评述本

流体力学标准化作业答案

流体力学标准化作业(三) ——流体动力学 本次作业知识点总结 1.描述流体运动的两种方法 (1)拉格朗日法;(2)欧拉法。 2.流体流动的加速度、质点导数 流场的速度分布与空间坐标(,,)x y z 和时间t 有关,即 流体质点的加速度等于速度对时间的变化率,即 投影式为 或 ()du u a u u dt t ?==+??? 在欧拉法中质点的加速度du dt 由两部分组成, u t ??为固定空间点,由时间变化引起的加速度,称为当地加速度或时变加速度,由流场的不恒定性引起。()u u ??为同一时刻,由流场的空间位置变化引起的加速度,称为迁移加速度或位变加速度,由流场的不均匀性引起。 欧拉法描述流体运动,质点的物理量不论矢量还是标 量,对时间的变化率称为该物理量的质点导数或随体导数。例如不可压缩流体,密度的随体导数 3.流体流动的分类 (1)恒定流和非恒定流 (2)一维、二维和三维流动 (3)均匀流和非均匀流 4.流体流动的基本概念 (1)流线和迹线 流线微分方程 迹线微分方程 (2)流管、流束与总流 (3)过流断面、流量及断面平均流速

体积流量 3(/)A Q udA m s =? 质量流量 (/)m A Q udA kg s ρ=? 断面平均流速 A udA Q v A A ==? (4)渐变流与急变流 5. 连续性方程 (1)不可压缩流体连续性微分方程 (2)元流的连续性方程 (3)总流的连续性方程 6. 运动微分方程 (1)理想流体的运动微分方程(欧拉运动微分方程) 矢量表示式 (2)粘性流体运动微分方程(N-S 方程) 矢量表示式 21()u f p u u u t νρ?+?+?=+??? 7.理想流体的伯努利方 (1)理想流体元流的伯努利方程 (2)理想流体总流的伯努利方程 8.实际流体的伯努利方程 (1)实际流体元流的伯努利方程 (2)实际流体总流的伯努利方程 10.恒定总流的动量方程 投影分量形式

杜广生工程流体力学思考题答案

牛顿流体 作用在流体上的切向应力与它所引起的角变形速度之间的关系符合牛顿内摩擦定 律的流体, 1-2: 什么是连续介质模型?为什么要建立? 1) 将流体作为由无穷多稠密、没有间隙的流体质点构成的连续介质,于是可将流体视为 在时间和空间连续分布的函数。 2) ①可以不考虑流体复杂的微观粒子运动,只考虑在外力作用下的微观运动; ②可以用连续函数的解析方法等数学工具去研究流体的平衡和运动规律。 1-3:流体密度、相对密度概念,它们之间的关系? 1) 密度:单位体内流体所具有的质量,表征流体的质量在空间的密集程度。 相对密度:在标准大气压下流体的密度与4℃时纯水的密度的比值。 关系: 1-4:什么是流体的压缩性和膨胀性? 1) 压缩性:在一定的温度下,单位压强增量引起的体积变化率定义为流体的压缩性系 数,其值越大,流体越容易压缩,反之,不容易压缩。 2) 膨胀性:当压强一定时,流体温度变化体积改变的性质称为流体的膨胀性 1-5:举例说明怎样确定流体是可压缩还是不可压缩的? 气体和液体都是可压缩的,通常将气体时为可压缩流体,液体视为不可压缩流体。 水下爆炸:水也要时为可压缩流体;当气体流速比较低时也可以视为不可压缩流体。 1-6:什么是流体的黏性?静止流体是否有黏性? 1) 流体流动时产生内摩擦力的性质程为流体的黏性 2) 黏性是流体的本身属性,永远存在。 3) 形成黏性的原因:1流体分子间的引力,2流体分子间的热运动 1-7:作用在流体上的力有哪些? 质量力、表面力。 表面张力,是液体表面层由于分子引力不均衡而产生的沿表面作用于任一界线上的张力 毛细现象:由于内聚力和附着力的差别使得微小间隙的液面上升和下降的现象。 上升和下降的高度与流体的种类,管子的材料、液体接触的气体种类和温度有关 2-1: 流体静压强有哪些特性 ? 1) 特性一:流体静压强的作用方向沿作用面的内法线方向 特性二:静压强与作用面在空间的方位无关,只是坐标点的连续可微函数 2-2:流体平衡微分方程的物理意义是什么? 在静止流体内的任一点上,作用在单位质量流体上的质量力与静压强的合力相平衡 2-3:什么是等压面?等压面的方程是什么?有什么重要性质? 1) 在流体中压强相等的点组成的面。 2) 0=++dz f dy f dx f z y x 3) 性质:在静止流体中,作用于任意点的质量力垂直于经过该点的等压面。 水头:单位重量流体所具有的能量也可以用液柱高来表示 帕斯卡原理:施于在重力作用下不可压缩流体表面上的压强,将以同样大小传到液体内部任 意点上 2-4:写出流体静力学基本方程的几种表达式。说明流体静力学基本方程的适用范围以及物 理意义、几何意义。 w f d ρρ=

流体力学大作业

流体力学-大作业

————————————————————————————————作者: ————————————————————————————————日期:

一.选择题 1.牛顿内摩擦定律适用于()。 A.任何流体B.牛顿流体 C.非牛顿流体 2.液体不具有的性质是()。 A.易流动性B.压缩性C.抗拉性 D.粘滞性 3连续介质假定认为流体()连续。 A.在宏观上 B.在微观上 C.分子间D.原子间 4.在国际单位制中流体力学基本量纲不包括()。 A.时间 B.质量 C.长度D.力. 5.在静水中取一六面体,作用在该六面体上的力有() A.切向力、正压力B.正压力C.正压力、重力D.正压力、切向力、重力 6. 下述哪些力属于质量力() A.惯性力B.粘性力 C.弹性力D.表面张力E.重力 7.某点存在真空时,( )() A.该点的绝对压强为正值 B.该点的相对压强为正值c.该点的绝对压强为负值 D.该点的相对压强为负值 8.流体静压强的( )。 A.方向与受压面有关B.大小与受压面积有关B.大小与受压面方位无关 9.流体静压强的全微分式为()。 A. B. C. 10.压强单位为时,采用了哪种表示法()。 A.应力单位B.大气压倍数C.液柱高度 11.密封容器内液面压强小于大气压强,其任一点的测压管液面( )。 A.高于容器内液面B.低于容器内液面 C.等于容器内液面 12.流体运动的连续性方程是根据( )原理导出的。 A.动量守恒 B. 质量守恒 C.能量守恒 D. 力的平衡 13.流线和迹线重合的条件为()。

A.恒定流B.非恒定流C.非恒定均匀流 14.总流伯努利方程适用于()。 A.恒定流 B.非恒定流C.可压缩流体 15. 总水头线与测压管水头线的基本规律是:( )、( ) A.总水头线总是沿程下降的。B.总水头线总是在测压管水头线的上方。 C.测压管水头线沿程可升可降。 D.测压管水头线总是沿程下降的。 16 管道中液体的雷诺数与()无关。 A.温度B.管径C. 流速D.管长 17.. 某圆管直径d=30mm,其中液体平均流速为20cm/s。液体粘滞系数为0.0114cm3/s,则此管中液体流态为( )。 A. 层流 B. 层流向紊流过渡C.紊流 18.等直径圆管中紊流的过流断面流速分布是()A呈抛物线分布B.呈对数线分布 C.呈椭圆曲线分布D.呈双曲线分布19.等直径圆管中的层流,其过流断面平均流速是圆管中最大流速的() A 1.0倍 B.1/3倍C.1/4倍D. 1/2倍 20.圆管中的层流的沿程损失与管中平均流速的()成正比. A. 一次方 B.二次方 C. 三次方D. 四次方 21..圆管的水力半径是() A. d/2B.d/3 C. d/4D. d/5. 22谢才公式中谢才系数的单位是()A.无量纲B.C.D.. 23.判断层流和紊流的临界雷诺数是() A.上临界雷诺数 B.下临界雷诺数 C.上下临界雷诺数代数平均 D.上下临界雷诺数几何平均 24..对于管道无压流,当充满度分别为( )时,其流量和速度分别达到最大。A.0.5,0.5B.0.95,0.81 C.0.81, 081 D. 1.0,1.0 25.对于a, b,c三种水面线,下列哪些说法是错误( )() A.所有a、c型曲线都是壅水曲线,即,水深沿程增大。B.所有

工程流体力学课后练习题答案

工程流体力学练习题 第一章 1-1解:设:柴油的密度为ρ,重度为γ;40 C 水的密度为ρ0,重度为γ0。则在同一地点的相对密度和比重为: 0ρρ=d ,0 γγ=c 30/830100083.0m kg d =?=?=ρρ 30/81348.9100083.0m N c =??=?=γγ 1-2解:336/1260101026.1m kg =??=-ρ 3/123488.91260m N g =?==ργ 1-3解:269/106.191096.101.0m N E V V V V p p V V p p p ?=??=?-=?-=????-=ββ 1-4解:N m p V V p /105.210 41010002956 --?=?=??-=β 299/104.0105.211m N E p p ?=?==-β 1-5解:1)求体积膨涨量和桶内压强 受温度增加的影响,200升汽油的体积膨涨量为: ()l T V V T T 4.2202000006.00=??=?=?β 由于容器封闭,体积不变,从而因体积膨涨量使容器内压强升高,体积压缩量等于体积膨涨量。故: 26400/1027.16108.9140004 .22004.2m N E V V V V V V p p T T p T T ?=???+=?+?-=?+?-=?β 2)在保证液面压强增量个大气压下,求桶内最大能装的汽油质量。设装的汽油体积为V ,那么:体积膨涨量为:

T V V T T ?=?β 体积压缩量为: ()()T V E p V V E p V T p T p p ?+?=?+?=?β1 因此,温度升高和压强升高联合作用的结果,应满足: ()()???? ? ??-?+=?-?+=p T p T E p T V V T V V 1110ββ ()())(63.197108.9140001018.01200006.012001145 0l E p T V V p T =???? ?????-??+=???? ???-?+=β ()kg V m 34.1381063.19710007.03=???==-ρ 1-6解:石油的动力粘度:s pa .028.01.010028=?= μ 石油的运动粘度:s m /1011.39 .01000028.025-?=?==ρμν 1-7解:石油的运动粘度:s m St /1044.0100 4025-?===ν 石油的动力粘度:s pa .0356.010 4100089.05=???==-ρνμ 1-8解:2/1147001 .01147.1m N u =?==δμ τ 1-9解:()()2/5.1621196.012.02 15.0065.021m N d D u u =-?=-==μδμτ N L d F 54.85.16214.01196.014.3=???=???=τπ 第二章 2-4解:设:测压管中空气的压强为p 2,水银的密度为1ρ,水的密度为2ρ。在水银面建立等压面1-1,在测压管与容器连接处建立等压面2-2。根据等压面理论,有 21p gh p a +=ρ (1) gz p z H g p 2221)(ρρ+=++(2) 由式(1)解出p 2后代入(2),整理得:

最新杜广生工程流体力学思考题答案

牛顿流体作用在流体上的切向应力与它所引起的角变形速度之间的关系符合牛顿内摩擦定律的流体, 1-2: 什么是连续介质模型?为什么要建立? 1)将流体作为由无穷多稠密、没有间隙的流体质点构成的连续介质,于是可将流体视为在 时间和空间连续分布的函数。 2)① 可以不考虑流体复杂的微观粒子运动,只考虑在外力作用下的微观运动; ②可以用连续函数的解析方法等数学工具去研究流体的平衡和运动规律。 1-3:流体密度、相对密度概念,它们之间的关系? 1)密度:单位体内流体所具有的质量,表征流体的质量在空间的密集程度。相对密度:在 标准大气压下流体的密度与4℃时纯水的密度的比值。 关系:d= f w 1-4:什么是流体的压缩性和膨胀性? 1)压缩性:在一定的温度下,单位压强增量引起的体积变化率定义为流体的压缩性系 数,其值越大,流体越容易压缩,反之,不容易压缩。 2)膨胀性:当压强一定时,流体温度变化体积改变的性质称为流体的膨胀性 1-5:举例说明怎样确定流体是可压缩还是不可压缩的?气体和液体都是可压缩的,通常将气体时为可压缩流体,液体视为不可压缩流体。水下爆炸:水也要时为可压缩流体;当气体流速比较低时也可以视为不可压缩流体。 1-6:什么是流体的黏性?静止流体是否有黏性? 1)流体流动时产生内摩擦力的性质程为流体的黏性 2)黏性是流体的本身属性,永远存在。 3)形成黏性的原因:1 流体分子间的引力,2 流体分子间的热运动 1-7:作用在流体上的力有哪些? 质量力、表面力。 表面张力,是液体表面层由于分子引力不均衡而产生的沿表面作用于任一界线上的张力毛细现象:由于内聚力和附着力的差别使得微小间隙的液面上升和下降的现象。 上升和下降的高度与流体的种类,管子的材料、液体接触的气体种类和温度有关 2-1: 流体静压强有哪些特性? 1)特性一:流体静压强的作用方向沿作用面的内法线方向特性二:静压强与作用面在空 间的方位无关,只是坐标点的连续可微函数 2-2:流体平衡微分方程的物理意义是什么?在静止流体内的任一点上,作用在单位质量流体上的质量力与静压强的合力相平衡 2-3:什么是等压面?等压面的方程是什么?有什么重要性质? 1)在流体中压强相等的点组成的面。 2) f x dx f y dy f z dz 0 3)性质:在静止流体中,作用于任意点的质量力垂直于经过该点的等压面。 水头:单位重量流体所具有的能量也可以用液柱高来表示帕斯卡原理:施于在重力作用下不可压缩流体表面上的压强,将以同样大小传到液体内部任意点上 2-4:写出流体静力学基本方程的几种表达式。说明流体静力学基本方程的适用范围以及物理意义、几何意义。 z p c1;z1 p1z2 p2适用于不可压缩重力流体的平衡状态; 物理意义:当连续不可压缩的重力流体处于平衡状态时,在流体中的任意点上,单位重

相关文档
最新文档