两相流的一些介绍讲解
气液两相流
热物理量测试技术1 概述两相流广泛应用于热能动力工程、核能工程、低温工程以及航天领域等许多领域。
所谓两相流，广义上讲是指一种物质或两种物质在不同状态下的流动，其中气体和液体一起流动称为气液两相流。
对于两相流中的气液混合物，它们可以是同一种物质，即汽—液(如水和水蒸气)，也可以是两种不同的物质，即气—液(如水和空气混合物)。
气液两相流是一个相当复杂的问题，。
在单相流中，经过一段距离之后，就会建立一个稳定的速度场。
但对于两相流，例如蒸汽和水，则很难建立一个稳定的流动，因为在管道流动中有压降产生，由于此压降作用会产生液体的蒸发，所以在研究气液两相流时必须考虑两相间的传热与传质问题。
两相流学科还处于半经验半理论阶段，对于两相流的流动和传热规律进行研究时，除了依靠各种数学物理模型外，还要依靠实验，这就需要两者相结合从而更好地进行研究。
2 两相流压降测量[1]压降，即两相流通过系统时产生的压力变化，是两相流体流动过程中的一个重要参数。
保持两相流体流动所需的动力以及动力系统的容量和功率就取决于压降的大小。
一般说来，两相流体流动时产生的压降一般由三部分组成，即摩擦阻力压降、重位压降、加速压降，管道系统出现阀门、孔板等管件时，还需测量局部压降。
目前，常用差压计或传感器来测量两相流压降。
2.1 利用差压计测量压降应用差压计测量气液两相流压降的测量原理图如图1所示。
所测压降为下部抽头的压力与上部抽头压力之差。
在差压计的Z1截面上可列出压力平衡式如下：P1+(Z2−Z1)ρC g=P2+(Z4−Z3)ρC g+(Z3−Z1)ρM g（2.1）式中，ρC为取压管中的流体密度；ρM为差压计的流体密度。
由（2.1）可得：P1−P2=(Z3−Z1)g(ρM−ρC)+(Z4−Z2)ρC g（2.2）由上式可知，要算出压降P1−P2的值，必须知道取压管中的流体密度ρC和差压计读数Z3−Z1。
当管中流体不流动时：P1−P2=gρm(Z4−Z2)（2.3）式中，ρm为两相混合物平均密度。
两相流、多相流
两相流的概念及类型两相物质(至少一相为流体)所组成的流动系统。
若流动系统中物质的相态多于两个，则称为多相流，两相或多相流是化工生产中为完成相际传质和反应过程所涉及的最普遍的粘性流体流动。
通常根据构成系统的相态分为气液系、液液系、液固系、气固系等。
气相和液相可以以连续相形式出现，如气体-液膜系统；也可以以离散的形式出现,如气泡-液体系统，液滴-液体系统。
固相通常以颗粒或团块的形式处于两相流中。
两相流的流动形态有多种。
除了同单相流动那样区分为层流和湍流外，还可以依据两相相对含量（常称为相比）、相界面的分布特性、运动速度、流场几何条件（管内、多孔板上、沿壁面等）划分流动形态。
对于管内气液系统，随两相速度的变化，可产生气泡流、塞状流、层状流、波状流、冲击流、环状流、雾状流等形态；对于多孔板上气液系可以产生自由分散的气泡、蜂窝状泡沫、活动泡沫、喷雾等形态。
两相流研究的一个基本课题是判断流动形态及其相互转变。
流动形态不同，则热量传递和质量传递的机理和影响因素也不同。
例如多孔板上气液两相处于鼓泡状态时，正系统混合物（浓度增加时表面张力减低）的板效率(见级效率)高于负系统混合物（浓度增加时表面张力增加）；而喷射状态下恰好相反。
两相流研究的另一个基本课题，是关于分散相在连续相中的运动规律及其对传递和反应过程的影响。
当分散相液滴或气泡时，有很多特点。
例如液滴和气泡在运动中会变形，在液滴或气泡内出现环流，界面上有波动，表面张力梯度会造成复杂的表面运动等。
这些都会影响传质通量，进而影响设备的性能。
两相流研究的课题，还有两相流系统的摩擦阻力，系统的振荡和稳定性等。
两相流研究模型两相流的理论分析比单相流困难得多，描述两相流的通用微分方程组至今尚未建立。
大量理论工作采用的是两类简化模型：①均相模型。
将两相介质看成是一种混合得非常均匀的混合物，假定处理单相流动的概念和方法仍然适用于两相流，但须对它的物理性质及传递性质作合理的假定；②分相模型。
两相流及几种模型介绍
两相流：通常把含有大量固体或液体颗粒的气体或液体流动称为两相流；其中含有多种尺寸组颗粒群为一个“相”，气体或液体为另一“相”，由此就有气—液，气—固，液—固等两相流之分。
两相流的研究：对两相流的研究有两种不同的观点：一是把流体作为连续介质，而把颗粒群作为离散体系；而另一是除了把流体作为连续介质外，还把颗粒群当作拟连续介质或拟流体。
引入两种坐标系：即拉格朗日坐标和欧拉坐标，以变形前的初始坐标为自变量称为拉格朗日Langrangian 坐标或物质坐标；以变形后瞬时坐标为自变量称为欧拉Eulerian 坐标或空间坐标。
离散相模型FLUENT在求解连续相的输运方程的同时，在拉格朗日坐标下模拟流场中离散相的第二相；←离散相模型解决的问题：煤粉燃烧、颗粒分离、喷雾干燥、液体燃料的燃烧等；←应用范围：FLUENT中的离散相模型假定第二相体积分数一般说来要小于10-12%（但颗粒质量承载率可以大于10-12%，即可模拟离散相质量流率等/大于连续相的流动）；不适用于模拟在连续相中无限期悬浮的颗粒流问题，包括：搅拌釜、流化床等；←颗粒-颗粒之间的相互作用、颗粒体积分数对连续相的影响未考虑；湍流中颗粒处理的两种模型：Stochastic←Tracking，应用随机方法来考虑瞬时湍流速度对颗粒轨道的影响；Cloud Tracking，运用统计方法来跟踪颗粒围绕某一平均轨道的湍流扩散。
通过计算颗粒的系统平均运动方程得到颗粒的某个“平均轨道”多相流模型FLUENT中提供的模型：VOF模型(Volume of Fluid Model)←混合模型(Mixture Model)←←欧拉模型(Eulerian Model)VOF模型(Volume of Fluid Model)← VOF模型用来处理没有相互穿插的多相流问题，在处理两相流中，假设计算的每个控制容积中第一相的体积含量为α1，如果α1=0，表示该控制容积中不含第一相，如果α1=1，则表示该控制容积中只含有第一相，如果0<α1<1，表示该控制容积中有两相交界面；← VOF方法是用体积率函数表示流体自由面的位置和流体所占的体积，其方法占内存小，是一种简单而有效的方法。
气液两相流 整理
第一章概论相的概念：相是体系中具有相同化学组成和物理性质的一部分，与体系的其它均匀部分有界面隔开两相流动的处理方法：双流体瞬态模拟方法和精确描述物理现象的稳态机理模型是多相管流研究的主要方法目前研究存在的问题：1、多相流问题未得到解析解；2、油气水三相流的研究不够深入；3、水平井段变质量流动研究较少；4、缺乏向下流动的综合机理模型；5、缺乏专用研究仪器气液两相流的分类：1、细分散体系：细小的液滴或气泡均匀分散在连续相中2、粗分散体系：较大的气泡或液滴分散在连续相中3、混合流动型：两相均非连续相4、分层流动：两相均为连续相气液两相流的基本特征：1、体系中存在相界面：两相之间也存在力的作用，出现质量和能量的交换时伴随着机械能的损失2、两相的分布情况多种多样：两相流动中两相介质的分布称为流型3、两相流动中存在滑脱现象：相间速度的差异称为滑脱，滑脱将产生附加的能量损失4、沿程流体体积流量有很大变化，质量流量不变气液两相流研究方法：1、经验方法：从气液两相流动的物理概念出发，或者使用因次分析法，或者根据流动的基本微分方程式，得到反映某一特定的两相流动过程的一些无因次参数，然后依据实验数据整理出描述这一流动过程的经验关系式。
优点：使用方便，在一定条件下能取得好的结果缺点：使用有局限性，且很难从其中得出更深层次的关系2、半经验方法：根据所研究的气液两相流动过程的特点，采用适当的假设和简化，再从两相流动的基本方程式出发，求得描述这一流动过程的函数关系式，最后用实验方法确定出函数关系式中的经验系数。
优点：有一定的理论基础，应用广泛缺点：存在简化和假设，具有不准确性3、理论分析方法：针对各种流动过程的特点，应用流体力学方法对其流动特性进行分析，进而建立起描述这一流动过程的解析关系式。
优点：以理论分析为基础，可以得到解析关系式缺点：建立关系式困难，求解复杂研究气液两相流应考虑的几个问题：1、不能简单地用层流或紊流来描述气液两相流2、水平或倾斜流动是轴不对称的3、由于相界面的存在增加了研究的复杂性4、总能量方程中应考虑与表面形成的能量问题5、多相流动中各相的温度、组分的浓度都不是均匀的，相之间有传热和传质6、各相流速不同，出现滑脱问题，是多相流研究的核心与重点流动型态：相流动中两相介质的分布状况称为流型或两相流动结构流型图：描述流型变化及其界限的图。
两相流、多相流讲课讲稿
两相流、多相流两相流的概念及类型两相物质(至少一相为流体)所组成的流动系统。
若流动系统中物质的相态多于两个，则称为多相流，两相或多相流是化工生产中为完成相际传质和反应过程所涉及的最普遍的粘性流体流动。
通常根据构成系统的相态分为气液系、液液系、液固系、气固系等。
气相和液相可以以连续相形式出现，如气体-液膜系统；也可以以离散的形式出现,如气泡-液体系统，液滴-液体系统。
固相通常以颗粒或团块的形式处于两相流中。
两相流的流动形态有多种。
除了同单相流动那样区分为层流和湍流外，还可以依据两相相对含量（常称为相比）、相界面的分布特性、运动速度、流场几何条件（管内、多孔板上、沿壁面等）划分流动形态。
对于管内气液系统，随两相速度的变化，可产生气泡流、塞状流、层状流、波状流、冲击流、环状流、雾状流等形态；对于多孔板上气液系可以产生自由分散的气泡、蜂窝状泡沫、活动泡沫、喷雾等形态。
两相流研究的一个基本课题是判断流动形态及其相互转变。
流动形态不同，则热量传递和质量传递的机理和影响因素也不同。
例如多孔板上气液两相处于鼓泡状态时，正系统混合物（浓度增加时表面张力减低）的板效率(见级效率)高于负系统混合物（浓度增加时表面张力增加）；而喷射状态下恰好相反。
两相流研究的另一个基本课题，是关于分散相在连续相中的运动规律及其对传递和反应过程的影响。
当分散相液滴或气泡时，有很多特点。
例如液滴和气泡在运动中会变形，在液滴或气泡内出现环流，界面上有波动，表面张力梯度会造成复杂的表面运动等。
这些都会影响传质通量，进而影响设备的性能。
两相流研究的课题，还有两相流系统的摩擦阻力，系统的振荡和稳定性等。
两相流研究模型两相流的理论分析比单相流困难得多，描述两相流的通用微分方程组至今尚未建立。
大量理论工作采用的是两类简化模型：①均相模型。
将两相介质看成是一种混合得非常均匀的混合物，假定处理单相流动的概念和方法仍然适用于两相流，但须对它的物理性质及传递性质作合理的假定；②分相模型。
两相流的一些介绍
三.结果分析
因为网格形状及其质量对两相流的计算收敛性 影响非常大，所以在计算时，我们采用了两套 网格进行计算即六面体网格和四面体网格，然 后比较计算结果。在两相流计算过程中有主相 和第二相的设置区别，为了考察两种设置的关 系，我进行了水是第二相及空气是第二相的两 种设定计算。
Models -> solve.. Base---基于压力求 解器 2. 隐形格式 ---Implicit 3.非定常流动-Unsteady
设置：1Pressure
Your site here
2.多相流模型
步骤：Define -> Models > Multiphase...
设置：混合模型---Mixture 原因：VOF 模型适合于分层
四面体-空气为第 二相-水的密度
Your site here
四面体网格各相面体-水为第二相-水的体积分数
四面体-空气为第二相-空气的体积分数
四面体-空气为第二相-水的体积分数
Your site here
六面体网格的各相密度
六面体-水为第 二相-混合密度
速度对比
四面体网格 水为第二相
四面体网格 空气为第二相
六面体网格 水为第二相
六面体网格 空气为第二相
Your site here
四面体网格各相密度
四面体-水为第 二相-混合密度
四面体-水为第 二相-空气密度
四面体-水为第 二相-水的密度
四面体-空气为第 二相-混合密度
四面体-空气为第 二相-空气密度
两相之间没有滑移速度，因为这个 过程是两相流是静止不动的，物体 一定速度冲入两相中。
的或自由表面流，mixture 和Eulerian 模型适合于流 动中有相混合或分离，或者 分散相的volume fraction 超过 10%的情形。（流动 中分散相的volume fraction 小于或等于 10% 时可使用离散相模型）
热工水力学(第八讲)两相流的描述
32
4.1.1 质量方程
∂ ∂ ∫∫Az ρdAz + ∂z ∫∫Az ρυ z dAz = 0 ∂t
或
∂ ∂ { ρ vα + ρl (1 − α )} Az + ∂z { ρ vαυv,z + ρl (1 − α )υl, z } Az = 0 ∂t
33
4.1.2 动量方程
p p+ z dz
3
1 两相流算子
• 算子(operator)表示某种函数运算的符号。 (operator) • 两相流算子，包括定相函数、面积平均算 子、体积平均算子、时间平均算子。 • 为什么需要三种不同的平均： 两相流不同于单相流，某一相可能在空间 上不连续，在时间上有脉动。
4
1.1 定相函数
• 如果某一时刻空间某一点r处是k相状态，k 表示汽（v）或液（l）。那么定相函数定义 定相函数 为：
{ jv } z =
qm,v,z
ρ v Az
=
Gv,z
ρv
=
Gm,z χ z
ρv
25
• 混合物：
{ j} z
1− χz χz 1 qm ,l, z qm ,v, z = ( + ) = Gm , z ( + ) Az ρl ρv ρl ρv
体积流密度具有线速度的量纲，在单相流 情况下实际上就是平均速度vm。
k
= ck
k
{c}
c =
k
、
{ } c 、 {{c}} = {c}
c
实际上没有意义
ck = ck
α k , {ck } = {ck }k {α k }
13
• 体积平均量若随时间波动，则瞬时量可表 示为时间平均量与脉动量之和：
2-2.井筒气液两相流
如果忽略气体的密度，则
m
fl f
l
液相的流动断面增大将引起混合物密度的增加。
假定：存在和不存在滑脱两种情况下液、气体积流量不
变。有滑脱时，气体流速大，液体流速小，为保持体
积流量不变，气体过流断面将减小为 f g ，液体的过流
断面将增加为 f1。考虑滑脱后分相过流断面的变化：
f fl fl ( fg fg )
• 泡雷诺数
vs C1C2 gD
Nb
vs DL L
-
图 2 21
摩 擦 阻 力 系 数 曲 线
图2-22 C1 ～Nb曲线
• C2是根据泡雷诺数及雷诺数 N R e来确定：
N
Re
vt DL L
•
段塞流的摩擦梯度：
f
fLvt2 ( qL vs Ap )
2D qt vs Ap
泡流的特点：气体是分散相，液体是连续相；气体主
要影响混合物密度，对摩擦阻力的影响不大；滑脱现 象比较严重。
③段塞流：当混合物继续向上流动，压力逐渐降低，气
体不断膨胀，小气泡将合并成大气泡，直到能够占据 整个油管断面时，井筒内将形成一段油一段气的结构。
段塞流的特点：气体呈分散相，液体呈连续相，炮弹
状的大气泡托着油柱向上流动，象一个破漏的活塞向 上推油。油、气间的相对运动要比泡流小，滑脱也小。
(4)雾流
• 雾流混合物密度计算公式与泡流相同：
m H L L H g g (1 H g )L H g g
• 由于雾流的气液无相对运动速度，即滑脱速度接近于雾，
基本上没有滑脱:
Hg
qg qL qg
• 摩擦梯度:
5.3两相流
两相流的均匀流模型
两相流的均匀流模型是把气（汽）－液两相混 合物看作一种均匀介质，其流动特性参量取两Leabharlann 相的相应参量的平均值。均匀流模型假定：
(1) 气（汽）相和液相的速度相等，即S＝1；
(2) 两相之间处于热力学平衡态，即TL＝TG； (3) 摩擦系数使用单相流的公式计算。
均匀流模型把复杂的两相流动作为单相流动处理，掩 盖了两相流中的一些复杂的问题。因此这种模型比较 粗糙，适用于两相质量流密度较大和压力较高的情况， 符合泡状流型，沫状和雾状流型。
L 1 x ( 1 ) G
两相流的局部压降
孔板 管道 uG uL
孔板
uG A uL A0 AC
pTP 1.26 pL pG
2
1 L 1) 1 x( 2 L A2 A1 A2 G
2
pE ,TP
m2 1 1 L 1) 1 x( 2 L A1 A2 G
两相流的局部压降-截面突然缩小
说明：两相流在流经突然 缩小的截面时，有一个 静压力降。
1 vf 1 x 1 S x v g 1 1 x 1 x
xe
H H
fs
H fg
描述两相流的特征参量2
折算速度：假定两 相流中的某一相介 质单独流过该通道 截面积A时的速度。 两相混合物速度
两相流基本方程二：动量守恒方程
dp dp dp dp dz dz F dz A dz G
摩擦压 力梯度
加速速度 压力梯度
重力压 力梯度
两相流基本方程三：能量守恒方程
两相流_中文
流体流动模型: 两相流两相流简介以下例子分析多孔介质中的两相流. 描述流体如何不发生混合同时通过一个孔隙可以回答很多环境和工业问题. 不幸的是, 多相分析很复杂需要求解一些未知量以及多个相关变量. 多相彼此之间的液体属性与每个液相压力和饱和级有关.这个问题验证两相流按照美国环境保护署的实验结果(参考文献1). 这个试验表明实验和数值计算的两相流结果相符. 通过这些实验, 科研人员估计流体对(空气-水, 空气-石油, 和石油-水)变化的流动, 然后使实验结果和计算机使用保持和渗透解析表达式数值计算的结果相匹配. 这个讨论根据Lincoln土壤和使用Mualem (参考文献2)和van Genuchten(参考文献3)给出的流体属性解析公式展开.图2-42: 在Hopmans和其他(参考文献1)两相流实验几何结构.这个例子分多步. 第一部分建立水和空气的两相流模型; 方程求解压力. 饱和度随解而变化. 需要强调假设在整个分析时间里土体中存在一些残余空气和水. 模型追踪空气前锋移动过湿液体通过观察饱和度而不是假设离散界面. 第二部分修改空气－水属性对空气－石油和空气－水系统进行仿真.模型定义在空气和水的实验中, 插入空气到充满水和沙的实验柱体中. 进入的空气 (对该流体对为非湿相)强迫水(湿相)向柱体基础出口流动. 在入口, 空气压力随时间逐步增加, 并且在柱体顶端没有水存在. 在向出口流动, 水通过一个空气流动不可渗透的圆盘. 空气和水都不能通过柱体垂向壁面. 出口处的水压随时间变化, 对应着滴定管中液体高度增加. 柱体总体高度8.34 cm, 半径6cm, 并且圆盘厚度0.74 cm. 实验时间170 hours.以下讨论给出两相流仿真计算的步骤. 首先观察控制方程和定义流体保持和渗透的解析关系(参考文献2和参考文献3). 然后给出一些详细的执行过程和模型数据表. 后面给出结果. 最后在图形用户界面逐步建立模型.控制方程和边界条件多孔介质两相流满足湿(w)和非湿(nw)流体的独立方程:(2-4)(2-5)其中θs 为总气隙或者饱和体积率; Se 为有效饱和函数; t 为时间(s); κint 为多孔介质固有渗透率(m 2); k r 为给定液体的相对渗透函数; η为流体动力粘度(kg/(m·s)); p 为压力(kg/(m·s 2)); ρ为流体密度(kg/m 3); g 为重力加速度; 和D 为垂向坐标(例如, x , y , 或者z )(m).如果流体分布为连续的, 流体在整个分析时间内未完全填充土壤, 给定湿相容积率, θw , 和非湿项, θnw . 对湿相, θ从0或者很小的残余值θr 到总体孔隙率, θs . 有效饱和度, Se , 源自比例系数θ关于θs 和θr 从0到1变化.θ和Se 为系统中所有液体压力的函数. 定义毛细压力:(2-6).气孔可以在给定时间内被一种气体完全填充: (2-7).饱和度随毛细压力如何变化, 如下所示(2-8)其中C 为指定容积, 下表“p ”表示单位压力.使用方程2-6, 方程2-7, 和方程2-8在方程2-4和方程2-5简化数值模型. 控制方程变为:(2-9)(2-10).可以同时求解系统方程得到p w 和p nw . 这个例子中, 两种流体不可压缩, 但是无需如此假设.初始, 柱体中的水和空气服从静力分布. 边界条件允许水仅存在于土壤柱体的基础中. 对于湿相, 边界条件为(2-11)其中n为边界法向的单位向量.因为空气从柱体顶端进入但是没有出来, 所以非湿项的边界条件为(2-12)保持和渗透关系可以从试验数据, 任意数学公式, 和模型中的其他方程内插建立两相流分析并定义系数θ, C, Se, k r, 和p c如何同时变化. 已有模型从参考文献2和参考文献3中的保持和渗透关系表示系数θ, C, Se, 和k r作为p c函数的变化. 因为p c较大并且由于θ, C, Se,和k r的变化较小, 这些表达式转化为与相同水柱等高的毛细压力或者毛细压力头如H c = p c/(ρwater g). 水压属性相对湿流体为(2-13)其中α, n, m, 和L为van Genuchten参数表示土壤属性. 注意两相流, van Genuchten-Mualem公式与H c 有关.对于非湿流体, 属性为(2-14)对湿相的定义同理可得.不同流体对在空气-水, 空气-油, 和油-水之间转换试验, 可以根据Leverett (参考文献4)灵活使用界面张力比例. 不考虑流体对Leverett比例调整柱体顶端的非湿相压力生成相同体积的湿流体从柱体底端流出. 使用With Leverett比例, 在流体对之间转换要求使用正确的流体对的流体属性ρ和η并根据以下公式调整边界和初始压力在这些方程中, σ表示不同流体之间的界面张力, 并且下标表示流体对. 这些值出现在这部分内容末尾的表格中. 例如, σao /σaw等于0.373, 和σwo/σaw等于0.534; 此外, 对空气－水系统第一项非湿相压力头(水以米为单位)为0.4 m, 空气－油系统为0.1 m, 和水－油系统为 0.2 m.因为多孔介质的相对渗透率和保持属性与流体流过该介质有关, 所以要求改变流体对的时候也改变模型的保持和渗透属性. 这个要求可以通过插入新试验数据或者调整数学公式实现. 这个模型中, 作者通过对解析公式拟合得到渗透率和保持率参数. 通过调整参数α, n, κs, 和θr为每种流体得到最佳拟合. 以下数据表显示不同流体对对应的比值α上面给出的σ比值近似相等.执行过程: 数值微分估计C这个例子使用解析表达式估计指定湿度, C. 因为C为q和H c的曲线斜率, 也可以通过使用COMSOL Multiphysics微分运算符计算得到C例如Cp,w(p w) = pdiff(theta_w,Hc)ρwater-1g-1因此可以书写任意表达式或者使用数据θ. 这个方程对流动体积率的微分考虑到毛细压力头. 此处除以水的重量将表达式转化为与压力变化相关. 相同类型的符号diff可进行微分运算. 两个运算符的差别在于diff可识别空间和时间的导数并且已经定义在COMSOL Multiphysics中(例如, p x, p y, p t), 然而pdiff为纯符号不能应用到与变量相关的链式法则中.执行过程: 逐步改变边界以下部分逐步减少通过使用内插函数定义入口处时变步进非湿相压力. 在COMSOL Multiphysics中内插直接可选. 打开选项菜单对话框, 建立时间和相应压力头表, 为插值函数定义名称, 并在其他需要使用函数的地方通过函数名调用, 在求解域设定对话框. 为了激活创建的函数, 输入带变量的函数名(例如, p_nw), 变量即括号中的时间t. 命令如下p_nw(t) = H_nw_t(t)*rho_water*g_w方程中出现水的密度因为参考文献1利用水柱高度定义边界压力.数据模型中使用的数据对应参考文献1给出的Lincoln 沙中空气－水的试验结果:不同流体对的van Genuchten 参数如下所示:根据下表入口处的空气压力头随时间增加.变量 描述 单位表达式g rAcceleration due to gravity m/s 2 9.82 ρfw Fluid density, water kg/m 31000ηw Dynamic viscosity, water Pa·s1·10-3ρfg Fluid density, gas kg/m 31.28ηg Dynamic viscosity, gasPa·s1.81·10-5κint Intrinsic permeability, column m 2 2.48·10-12 κsPermeability, discm 21.33·10-14 q s,w Saturated volume fraction, column0.32 q s,wSaturated volume fraction, disc0.5 p nw,top Initial nonwetting phasepressure head, inletm water0.2变量 描述 单位 空气-水 空气-油 油-水q r,w Residualvolume fraction0.021 0.00001 0.0072α alpha parameter m-1 1.89 5.29 3.58n n parameter, column 2.811 3.002 3.1365 LL parameter, column 0.5 0.5 0.5 κsPermeability, discm -22.48·10-121.09·10-120.94·10-12压力头(m 水) 起始时间(小时)0.4 0 0.6 21.25 0.8 45.251.0 691.5 932 122.54 155在水的出口处, 滴定管的流体级随时间从0 m到0.1 m线性增加.结果和讨论图2-43给出试验柱体从COMSOL Multiphysics两相流解的早期图像. 阴影描述非湿相的有效饱和(空气), 同时箭头给出湿相(水)速度.图2-43: 两相流模型在0.1小时的解: 非湿相饱和(表面图), 湿相速度(箭头). 结果对应从US EPA (参考文献1)中Lincoln土壤的空气－水试验结果.图像给出非湿相流体进入土壤柱体并取代湿相流体. 因为多步压力变化压迫非湿相进入.图2-44给出沿着柱体不同高度的毛细压力在入口边界处的步进压力头. 指定后处理中输入时指定观察点. 两相流问题的解与参考文献1的结果完美一致.毛细压力头标定空气入口压力如图2-44轨迹一致说明试验成功. 为了获得渗透率和保持属性的高分辨率解, 作者在参考文献1设置压力步长足够大使土体中的冲击变为瞬时. 如图2-45所示, 渗透率在整个土体中瞬时变化.图2-44: 在Lincoln土壤(源自US EPA, 参考文献1)中空气－水流动的入口压力头(实线)和毛细压力头(虚线).图2-45: Lincoln土壤在 x = 0.03 m (源自US EPA, 参考文献1)处水的渗透函数(实线)和空气的渗透函数(虚线). 空气－油和油－水系统的两相流动解如图2-46和图2-47所示.图2-46: Lincoln土壤(源自US EPA, 参考文献1)空气－油的入口空气压力头(实线)和毛细压力头(虚线).图2-47: Lincoln土壤 (源自US EPA, 参考文献1)油－水的入口空气压力头(实线)和毛细压力头(虚线).使用COMSOL Multiphysics计算的空气－油和油－水两相流问题的结果证明与参考文献1的结果完美一致. 通过Leverett比例系数可以设置入口压力以便空气－油和油－水系统从空气－水试验生成体积流出率. 与空气－水系统相同, 空气－油试验的毛细压力头和空气入口压力为瞬态的. 然而, 对水－油系统, 非湿相和湿相之间有滞后.参考文献1. J.W. Hopmans, M.E. Grismer, J. Chen, and Y.P. Liu, Parameter estimation of two-fluid capillary pressure saturation and permeability functions, U.S. Environmental Protection Agency EPA/600/R-98/046, Cincinnati, Ohio, 1998.2. Y. Mualem, “A new model for predicting the hydraulic permeability of unsaturated porous media,” Water Res. Research, vol. 12, 1976, pp. 513–522.3. M.Th. van Genuchten, “A closed-form equation for predicting the hydraulic of conductivity of unsaturated soils,” Soil Sci. Soc. Am. J., vol. 44, 1980, pp. 892–898.4. M.C. Leverett, “Capillary behavior in porous solids,” Trans. AIME, vol. 142, 1941, pp. 152–169.模型库路径: Earth_Science_Module/Fluid_Flow/twophase_flow_air_water 使用图形用户界面建模: 空气-水系统模型导航视窗1 打开模型导航视窗下的新增页面, 在空间维度列表选择2D.2 从应用模式列表选择地球科学模块>流体流动>达西定律>压力分析>瞬态分析.3 在因变量编辑区输入p_w, 并在应用模式名编辑区输入w.4 单击多物理场按钮, 然后单击新增.5 再次在应用模式列表选择地球科学模块>流体流动>达西定律>压力分析>瞬态分析.6 重复以下步骤, 改变因变量为p_nw和应用模式名为nw. 单击多物理场按钮, 然后单击新增.7 单击确定.几何建模绘制两个矩形创建几何结构.1 从绘图菜单选择指定对象>矩形.2 输入以下设置;设置完成, 单击确定.参数值Width 0.06Height 0.0074x 0z 03 除了以下设置重复上面的过程:参数值Width 0.060.0834选项Height x 0 z4 在主工具条单击缩放至视窗大小按钮.1 选择选项>常数.2 输入以下名称, 表达式, 和描述(任意的); 设置完成, 单击确定.名称表达式描述SPH 1[h/s] Seconds per hour rho_water 1000[kg/m^3]Density, waterrho_w 1000[kg/m^3]Density, wetting fluidetaw 0.001[Pa*s]*SPH Dynamic viscosity, wetting fluid rho_nw1.28[kg/m^3]Density, nonwetting fluid etanw 0.0000181[Pa*s]*SPHDynamic viscosity, nonwetting fluid3 从选项菜单选择表达式>标量表达式. 输入以下名称和表达式(都在一条线上); 设定完成, 单击确定.名称表达式p_nw_top 0.2[m]*rho_water*g_nwp_nw_init p_nw_top+(8.34[cm]-y)*rho_nw*g_nw p_w_init -rho_w*g_w*yp_w0 0.1[m]*rho_water*g_w*(t/170[s]) Hc (p_nw-p_w)/(rho_water*g_w)Se_w (1+abs(alpha*Hc)^N)^(-M)*(Hc>0)+1*(Hc<=0) theta_w (theta_r+Se_w*(theta_s-theta_r))*(Hc>0)+theta_s*(Hc<=0) kr_w ((Se_w^L*(1-(1-Se_w^(1/M))^M)^2)+eps)*(Hc>0)+1*(Hc<=0) Cp ((alpha*M/(1-M)*(theta_s-theta_r)*Se_w^(1/M)*(1-Se_w^(1/M))^M))/(rho_water*g_w)*(Hc>0) Se_nw 1-Se_w theta_nw theta_s-theta_wkr_nw ((1-Se_w)^L*(1-Se_w^(1/M))^(2*M))*(Hc>0)+eps4 从选项菜单选择表达式>求解域表达式. 输入以下名称和表达式为两个求解域; 设置完成, 单击确定.名称求解域1求解域2theta_s 0.5 0.32 theta_r 0.0211 kap_s 1.34e-14 2.48e-12alpha1.89N 2.811M 1-1/NL 0.5现在使用内插方法为非湿相边界建立步进压力.5 从选项菜单, 选择函数.6 单击新增按钮.7 在函数名编辑区输入Hp_nw_t并选择内插法选项.8 单击确定.9 输入下表中的值; 完成输入, 单击确定.X Y0 0.421.20 0.421.25 0.645.20 0.645.25 0.868.95 0.869 1.092.95 1.093 1.5122.45 1.5122.5 2154.95 2155 4200 4物理量下面部分首先为每一个相建立材料属性, 初始条件, 和边界条件, 然后将两个方程建立联系.应用标量变量建立模型垂向和重力常数, 转到物理量菜单并选择标量变量. 按照如下参数修改默认值; 完成设定, 单击确定.名称表达式D_w yg_w 9.82[m/s^2]*SPH^2D_nw yg_nw 9.82[m/s^2]*SPH^2这步将模型基本时间单位从秒转化为小时. 因此, 任何时候用户交互界面出现的时间单位为秒—既有显式, 如m/s, 或隐式的, 如Pa (= kg/(m·s2))—代替小时.求解域设定—湿相1 从多物理场菜单选择达西定律(w).边界条件—湿相从物理量选择边界设定, 然后设置以下边界条件; 完成设置, 单击确定.求解域设定—非湿相边界条件—非湿相方程系统2 从物理量菜单选择求解域设定. 输入以下设定, 然后单击应用. 设定 求解域1 求解域2Storage term User defined User definedS0 Cp k Skap_s kap_s*kr_w ρfrho_w rho_w η etaw etaw3 单击初始选项卡. 同时选择求解域1和2, 并在p_w(t 0)编辑区为压力输入p_w_init . 单击确定.设定 边界2 边界1, 3, 5–7Boundary conditionPressure Zero flux/Symmetry p 0 p_w01 从多物理场菜单选择达西定律(nw).2 从物理量菜单选择求解域设定. 在求解域选择列表选择1, 然后清除启动本域的复选框.3 选择求解域2. 输入以下设定.设定 求解域2Storage term User definedSCp κSkap_s*kr_nw ρfrho_nw η etanw4 单击初始选项卡. 同时选择求解域1和2, 并在p_nw(t 0)编辑区为压力输入p_nw_init . 单击确定.1 从物理量菜单, 选择边界设定.2 设置以下边界条件; 完成设置, 单击确定.设置 边界3, 4, 7 边界5BoundaryconditionZero flux/Symmetry Pressure p 0 Hp_nw_t(t)*rho_water*g_nw1 从物理量菜单, 选择方程式系统>求解域设定.在对话框信息中, COMSOL Multiphysics转化输入值为系数求解有限元模型. 框中包含每种模型类型: 2D问题中的求解域, 边界, 和点.2 单击d a选项卡. 此处其中系数为压力变化率乘以时间. 4-元素矩阵原本对角线为: 湿相压力, p w, 对湿相方程; 和非湿相压力, p nw, 对非湿相方程. 因为两相问题有储存项对非湿相和湿相的方程, 可以给矩阵填充合适的项. 修改矩阵如下所示:PW PNWCp+eps -Cp+eps-Cp+eps Cp+eps注意eps为一个非常小数.3 单击确定.网格生成1 从网格菜单, 选择自由网格参数.2 单击边界选项卡. 使用Ctrl键选择边界4和5. 在最大单元尺寸编辑区输入0.001.3 选择边界2, 并在最大单元尺寸编辑区输入0.005.4 单击重化网格按钮, 然后单击确定.求解1 从求解菜单, 打开求解器参数对话框. 在求解器列表选择稳态(如果还未选中).2 在时间编辑区, 输入0,0:0.01:0.1,0.1:0.1:1,1:170. 单击确定.3 在主工具条单击求解按钮.后处理和可视化按照如下步骤生成图2-43:1 从后处理菜单, 打开绘图参数对话框.2 在绘图类型区, 选择表面和箭头复选框. 在不同时间的解列表选择0.1.3 单击表面选贤卡, 并在表达式编辑区输入Se_nw, 其表示名称为nw的应用模式计算的有效饱和(Se).4 单击箭头选项卡. 在内建物理参数列表选择Darcy定律(w)>速度场. 转到箭头位置区. 在点数列, x点和y点编辑区输入25. 转到箭头参数区. 单击颜色按钮, 改变颜色为黑色, 并单击确定. 清除自动复选框, 并在比例刻度编辑区输入0.5.5 单击确定.继续如下步骤生成图2-44:1 从后处理菜单, 打开剖面绘图参数对话框.2 单击点选项卡. 在表达式编辑区, 输入Hc. 在坐标区, x编辑区输入0.03 0.03 0.03, y编辑区输入0.02 0.04 0.06.3 单击线设定按钮. 设定线样式为虚线, 然后单击确定.4 单击应用.5 在表达式编辑区, 输入p_nw/(rho_water*g_nw). 在坐标区, x编辑区输入0.03, 在 y编辑区输入0.0834.6 单击线设定按钮. 设定线样式为实线, 然后单击确定.7 单击通用选项卡. 选择保留目前的图复选框, 然后单击应用.继续下面步骤生成图2-45:1 仍然在通用页面, 从绘图于列表选择新图.2 单击点选项卡. 在表达式编辑区输入kr_nw.3 在坐标区, x编辑区输入0.03 0.03 0.03, y编辑区输入0.02 0.04 0.06.4 单击线设定按钮. 设定线样式为虚线, 然后单击确定.5 单击应用.6 在表达式编辑区, 输入kr_w.7 单击线设定按钮. 设定线样式为实线, 然后单击确定.8 单击通用选项卡. 选择保留目前的图复选框, 然后单击确定.两相流: 改变流体对这是两相流动例子的第二部分. 这里可以修改例子在第一部分创建的空气－水模型文件以仿真计算空气－油和油－水两相系统的流动特征, 参见参考文献1. 在空气－油系统中, 油为湿相而空气为非湿相. 在油－水系统中, 水为湿相.将空气－水系统转变为空气－油或者油－水系统需要改变流体属性和几个多孔介质参数. 为了获得不同流体对生成湿相流出率其仿真过程与空气－水系统类似, 如前面讨论的在入口处界面张力随压力比例增加模型库路径: Earth_Science_Module/Fluid_Flow/twophase_flow_air_oil模型库路径: Earth_Science_Module/Fluid_Flow/twophase_flow_oil_water 使用图形用户界面建模: 空气－油和油－水以下部分打开文件twophase_flow_air_water.mph并修改为空气－油和油－水系统. 下面给出的表和说明描述了如何改变参数以满足两相系统仿真的要求.模型导航视窗1 打开模型导航视窗并单击模型数据库选项卡.2 在模型数据库列表中选择地球科学模块>流体流动>twophase flow air water.3 单击确定.选项物理量—非湿从多物理场菜单, 选择达西定律(nw).边界条件—非湿相1 从选项菜单, 打开常数对话框. 按下表所示修改密度, 粘度和界面张力; 完成设定, 单击确定. 约束 空气-油 油-水rho_w 800[kg/m^3]1000[kg/m^3] etaw 0.00392[Pa*s]*SPH 0.001[Pa*s]*SPH rho_nw 1.28[kg/m^3]800[kg/m^3] etanw 0.0000181[Pa*s]*SPH 0.00392[Pa*s]*SPH sigma_ao 0.0259sigma_aw 0.06810.0681 sigma_ow 0.0364 2 选择选项>表达式>标量表达式.3 修改以下表达式; 设定完成, 单击确定.项 空气-油 油-水p_nw_top 0.2[m]*rho_water*g_nw*sigma_ao/sigma_aw0.2[m]*rho_water*g_nw* sigma_ow/sigma_aw p_w0 0.1[m]*rho_w*g_w*(sigma_ao/ sigma_aw)*(t/170[s])0.1[m]*rho_w*g_w*(sigma_ow/sigma_aw)*(t/170[s])4 选择选项>表达式>求解域表达式.5 在求解域选择列表, 选择2.6 改变如下表达式的值; 完成设定, 单击确定.名称 表达式 kap_s 0.94e-12theta_r 0.0072alpha 3.58 N 3.13651 从物理量菜单, 选择边界设定.2 在边界选择列表选择5.3 设置如下条件; 完成设定, 单击确定.边界条件 变量 表达式Pressure(for Air/Oil p 0 Hp_nw_t(t)*rho_water*g_nw*sigma_ao/sigma_aw求解后处理和可视化重复生成图2-44的步骤生成图2-46和图2-47.本文档由lhpvibration 整理 model)Pressure(for Oil/Watermodel)p 0 Hp_nw_t(t)*rho_water* g_nw*sigma_ow/sigma_aw1 从求解菜单, 打开求解器参数对话框.2 在求解器列表, 选择稳态(如果还未选中).3 转到时间步长区. 在时间编辑区输入0,0:0.01:0.1,0.1:0.1:1,1:170. 单击确定.4 在主工具条单击求解按钮.。
两相流、多相流
两相流的概念及类型两相物质(至少一相为流体)所组成的流动系统。
若流动系统中物质的相态多于两个，则称为多相流，两相或多相流是化工生产中为完成相际传质和反应过程所涉及的最普遍的粘性流体流动。
通常根据构成系统的相态分为气液系、液液系、液固系、气固系等。
气相和液相可以以连续相形式出现，如气体-液膜系统；也可以以离散的形式出现,如气泡-液体系统，液滴-液体系统。
固相通常以颗粒或团块的形式处于两相流中。
两相流的流动形态有多种。
除了同单相流动那样区分为层流和湍流外，还可以依据两相相对含量（常称为相比）、相界面的分布特性、运动速度、流场几何条件（管内、多孔板上、沿壁面等）划分流动形态。
对于管内气液系统，随两相速度的变化，可产生气泡流、塞状流、层状流、波状流、冲击流、环状流、雾状流等形态；对于多孔板上气液系可以产生自由分散的气泡、蜂窝状泡沫、活动泡沫、喷雾等形态。
两相流研究的一个基本课题是判断流动形态及其相互转变。
流动形态不同，则热量传递和质量传递的机理和影响因素也不同。
例如多孔板上气液两相处于鼓泡状态时，正系统混合物（浓度增加时表面张力减低）的板效率(见级效率)高于负系统混合物（浓度增加时表面张力增加）；而喷射状态下恰好相反。
两相流研究的另一个基本课题，是关于分散相在连续相中的运动规律及其对传递和反应过程的影响。
当分散相液滴或气泡时，有很多特点。
例如液滴和气泡在运动中会变形，在液滴或气泡内出现环流，界面上有波动，表面张力梯度会造成复杂的表面运动等。
这些都会影响传质通量，进而影响设备的性能。
两相流研究的课题，还有两相流系统的摩擦阻力，系统的振荡和稳定性等。
两相流研究模型两相流的理论分析比单相流困难得多，描述两相流的通用微分方程组至今尚未建立。
大量理论工作采用的是两类简化模型：①均相模型。
将两相介质看成是一种混合得非常均匀的混合物，假定处理单相流动的概念和方法仍然适用于两相流，但须对它的物理性质及传递性质作合理的假定；②分相模型。
01-两相流概述和基础
气液两相流与沸腾换热
绪论－什么是两相流与沸腾？
• 沸腾：沸腾是通过液气相变将工质由液态转换到汽态的
一种剧烈汽化过程，也是伴随大量气泡的形成、成长和运 动的热量传递过程。沸腾传热和气液两相流是由本质上十 分复杂的沸腾和两相流动两种物理现象耦合在一起的一种 热流体流动过程．在两相界面上存在着热量的传递，并且 有相变发生。
气液两相流与沸腾换热
绪论－两相流与沸腾传热的应用（8）
庞力平
华北电力大学能源动力与机械工程学院
气液两相流与沸腾换热
绪论－两相流与沸腾传热的应用（8）
庞力平
华北电力大学能源动力与机械工程学院
气液两相流与沸腾换热
绪论－两相流与沸腾传热的应用（8）
庞力平
华北电力大学能源动力与机械工程学院
气液两相流与沸腾换热
之为相。因此，由任意两种存在分界面的独立物质组 成的物体或系统都称之为两相物体或两相系统。例 如，水和冰的混合物为一种两相物体，因为水和冰都 是存在分界面的独立物质。但是，盐水溶液是一种单 相物体，因为在此溶液中盐和水之间无分界面，盐和 水不是两种独立存在的物质。
庞力平
华北电力大学能源动力与机械工程学院
庞力平
华北电力大学能源动力与机械工程学院
气液两相流与沸腾换热
管内气液两相流的基本参数
• 气液两相流平均密度ρm ρm = ρG α + ρL (1 − α ) 分流平均密度 ρ = W / Q = ρ β + ρ (1 − β) 混合平均密度 • 气液两相流平均速度um u m = u G α + u L (1 − α) 分流平均速度 u m = Q / A = G /[x / ρG + (1 − x) / ρL ] 混合平均速度
土力学中的两相流推导过程
土力学中的两相流推导过程土力学中的两相流是指在土体中同时存在两种不同相态的流体，通常是土体中的水和空气。
土力学中的两相流推导过程是指对两相流的相关参数进行计算和推导的过程。
我们需要了解土力学中的两相流的基本概念和特点。
在土体中存在两相流时，水和空气会以不同的形式和比例存在。
水以液态形式存在于土体的孔隙中，而空气以气态形式存在于土体的孔隙中。
两相流的存在会对土体的力学性质和行为产生重要影响，因此对两相流的计算和推导非常重要。
我们需要确定两相流的相关参数，例如饱和度、孔隙比、渗透率等。
饱和度是指土体中水的含水量与孔隙容积的比值，孔隙比是指土体中孔隙的容积与土体总体积的比值，渗透率是指土体中水的渗透能力。
通过测量和计算这些参数，我们可以获得土体中两相流的具体情况。
然后，我们可以利用土力学中的一些基本理论和经验公式来推导两相流的相关参数。
例如，达西定律可以用来计算土体中水的渗透速率，根据达西定律，渗透速率与渗透率和水头的差值成正比。
通过测量渗透速率和水头的差值，我们可以求得土体中水的渗透率。
我们还可以利用土力学中的两相流模型来推导两相流的相关参数。
例如，罗斯模型是常用的两相流模型之一，它考虑了土体中水和空气的相互作用和运动规律。
通过求解罗斯模型的方程，我们可以得到土体中两相流的速度场和压力场等参数。
我们还可以利用实验数据和数值模拟来验证和优化推导过程中得到的参数和结果。
实验数据可以通过在实验室或现场进行试验来获得，数值模拟可以通过计算机模拟来实现。
通过与实验数据和数值模拟结果的比较，我们可以验证和优化推导过程中得到的参数和结果，以提高推导的准确性和可靠性。
土力学中的两相流推导过程是通过计算、推导和实验验证来确定土体中两相流的相关参数和结果的过程。
通过对两相流的计算和推导，我们可以更好地理解土体中两相流的行为和特性，从而为土力学和土木工程等领域的研究和应用提供支持和指导。
两相流原理
两相流原理两相流原理是指在一定条件下，两种不同物质（例如液体和气体）同时存在于同一空间中并进行相互作用的现象。
在工程和科学领域中，对于两相流的研究具有重要意义，因为它在许多领域中都有广泛的应用，如核能工程、化工工艺、石油开采等。
两相流的形成通常是由于液体的蒸发或气体的凝结引起的。
在一定的温度和压力条件下，液体和气体可以同时存在于同一空间中。
液体和气体之间存在着物质的传递和能量的交换。
这种相互作用可以通过质量传递、能量传递和动量传递来描述。
在两相流中，液体和气体的相互作用形成了许多有趣的现象。
例如，在汽车发动机中，汽油的喷雾与气体的混合是发动机燃烧过程中关键的一步。
在核电站中，蒸汽与冷却剂之间的相互作用决定了核反应堆的热效率。
在化工工艺中，液体和气体的传递现象对反应速率和生产效率有着重要影响。
两相流的研究方法可以通过实验和数值模拟两种途径来进行。
实验方法通常是通过设计实验装置，模拟两相流现象，并通过观察和测量来获得数据。
数值模拟方法则是通过建立数学模型和计算流体力学方法来模拟和预测两相流的行为。
这两种方法各有优势和限制，需要根据具体情况选择合适的方法。
在实际工程中，两相流现象的研究和控制具有重要意义。
例如，在核电站中，需要确保蒸汽与冷却剂之间的传递过程稳定可靠，以保证核反应堆的安全运行。
在石油开采中，需要研究和控制油水两相流的行为，以提高开采效率和减少环境污染。
为了更好地理解两相流现象，研究人员还提出了许多物理模型和数学模型来描述和预测两相流的行为。
这些模型基于质量守恒、能量守恒和动量守恒等基本原理，并结合实验数据进行验证和修正。
通过不断改进和完善这些模型，可以更准确地预测和控制两相流的行为。
两相流原理是液体和气体在一定条件下相互作用的结果。
研究和控制两相流的行为对于许多工程和科学领域具有重要意义。
通过实验和数值模拟等方法，可以更深入地了解两相流的特性，并应用于实际工程中。
随着科学技术的不断发展，对于两相流的研究将会在各个领域中得到更广泛的应用。
8 两相流动
2
浮力 Fb 与流体运动无关
2010/4/26
流体对颗粒的形体曳力 Fp 正比于流速 u — — 形体曳力(Form drag)
5
曳力与曳力系数（Drag and drag coefficient） 流体流动对颗粒表面的总曳力为摩擦曳力与形体曳力之和
u p u ut
u = 0，up = ut 流体静止，颗粒向下运动； up = 0，u = ut ，颗粒静止地悬浮在流体中； u > ut ， up > 0， 颗粒向上运动； u < ut ， up < 0，颗粒向下运动。
2010/4/26 13
沉降速度的求法：
求沉降速度通常采用试差法。 ① 假设流体流动类型； ② 计算沉降速度； ③ 计算Re，验证与假设是否相符； ④ 如果不相符，则转①。如果相符，OK !
2010/4/26
2
— — 表面曳力 (Wall drag)4
曳力与曳力系数（Drag and drag coefficient） 流体静压强对整个球体表面的作用力在流动方向上的分量为
Fn d cos p
2
0
0
r R
R
2
sin d
3 u d p0 gR cos cos R 2 sin cos d 2 R 0 0
A
ap a
非球形颗粒4个几何参数之间的关系
d eV A d eA
d ea d eV
2
工程上多采用可以测量的等体积当量直径 deV 和具有直观意 2010/4/26 15 义的形状系数A。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
速度迹线
四面体-水是第二相-速度迹线
四面体-空气第二相-速度迹线
六面体-水是第二相-速度迹线
六面体-空气是第二相-速度迹线
Your site here
压力迹线
四面体-水为第二相-压力迹线
四面体-空气为第二相-压力迹线
六面体-水为第二相-压力迹线
六面体-空气为第二相-压力迹线
Your site here
3.湍流模型
步骤：Define -> Models -
>Viscous
设置：K-ε 模型
原因：Spalart-Allmaras 机
翼上的超音速、跨音速流动 ，边界层流动等 k-ε模型应用广泛用于 可压缩， 浮力，燃烧等。多相流动动 量方程中所模拟的项数是非 常大的，这使得多相流模拟 中的紊流模型非常复杂。这 一模型可以满足。 k-ω模型对于有压力梯度的大范 围边界层流动航天和涡轮机 械领域
原因：VOF 模型适合于分层
的或自由表面流，mixture
和Eulerian 模型适合于流
动中有相混合或分离，或者
分散相的volume fraction
超滑移速度，因为这个 过程是两相流是静止不动的，物体 一定速度冲入两相中。
中分散相的volume fraction 小于或等于 10% 时可使用离散相模型）
格式，是专门用于多相流计算的格式，适合大体积分数变化的流场 计算。其他解算器为保证计算精确，设为二阶迎风格式
Your site here
8.初始化设置
在初始化时，Y轴向速 度为863m/s，这个速度 值是流场在刚开始计算 时的冲击速度，随着计 算的深入，流场冲击速 度会逐渐增大，最后达 到预设的1700m/s。
Your site here
4.定义材料
Your site here
5.定义第二相
步骤：define->phase... 设置：1.主相--air
2.第二相--water_liquid
Your site here
6.边界条件设置与第二相体积分数设定
步骤： Define->Boundary Conditions...
Your site here
网格划分
四面体网格
六面面体网格
蓝色面为压力入口，红色面为压力出口，黄色面 为对称面，黑色的地面和RVE框架为WALL
Your site here
压力对比
四面体网格水为第二相
四面体网格 空气为第二相
六面体网格水为第二相
六面体网格 空气为第二相
Your site here
速度对比
四面体网格 水为第二相
四面体网格 空气为第二相
六面体网格 水为第二相
六面体网格 空气为第二相
Your site here
四面体网格各相密度
四面体-水为第 二相-混合密度
四面体-水为第 二相-空气密度
四面体-水为第 二相-水的密度
四面体-空气为第 二相-混合密度
四面体-空气为第 二相-空气密度
四面体-空气为第 二相-水的密度
Your site here
四面体网格各相体积分数
四面体-水为第二相-空气的体积分数
四面体-水为第二相-水的体积分数
四面体-空气为第二相-空气的体积分数
四面体-空气为第二相-水的体积分数 Your site here
六面体网格的各相密度
六面体-水为第 二相-混合密度
Your site here
质量守恒对比
Inlet
四面体网格 166770.31 水为第二相
四面体网格 空气为第二 相
六面体网格 水为第二相
166770.63 166757.72
六面体网格
空气为第二 相
166757.8
outlet
166770.34 166770.78
166757.41 166757.02
入口边界：pressure_inlet ---压力入口
压力计算：
p0 ps 1/ 2 | v |2
p0 总压......ps - - - 静压
- -混合密度.....v - -速度
出口边界：pressure-outlet压力出口
Your site here
6.边界条件设置与第二相体积分数设定
Your site here
9.putch设置
水的体积分数
初始化补充，此项设置为补充两相流的初始化，不经 过此项补充在计算时会默认为单项。
Your site here
三.结果分析
因为网格形状及其质量对两相流的计算收敛性 影响非常大，所以在计算时，我们采用了两套 网格进行计算即六面体网格和四面体网格，然 后比较计算结果。在两相流计算过程中有主相 和第二相的设置区别，为了考察两种设置的关 系，我进行了水是第二相及空气是第二相的两 种设定计算。
多相流模型
王新慧
2011.10.20
一.多相流概念与分类
定义：两种或两种以上不同相的流体混合在一
起的流动
第 二 相
分类：
1气-液或者液-液两
2.气-固两相流
3.液-固两相流 4.三相流
主 相
Your site here
二.问题设置和基本流程
Your site here
1.求解器的设置
步骤：Define ->
Your site here
6.边界条件设置与第二相体积分数设定
在边界条件界面中，选择第二 相设定窗口，可任意设定空气 或水的体积分数
第二相体检分数的设定
Your site here
7.解算器设置
步骤:solve->control->soultion 设置：discretion中体积分数--QUICK格式。QUICK体积分数解算器
net
0.03125 -0.15625
Models -> solve..
设置：1Pressure
Base---基于压力求 解器
2. 隐形格式 ---Implicit
3.非定常流动-Unsteady
Your site here
2.多相流模型
步骤：Define -> Models -
> Multiphase...
设置：混合模型---Mixture
六面体-水为第 二相-空气密度
六面体-水为第 二相-水的密度
六面体-空气第 二相-混合密度
六面体-空气第 二相-空气密度
六面体-空气第 二相-水的密度Your site here
六面体网格的各相体积分数
六面体-水为第二相-空气的体积分数
六面体-水为第二相-水的体积分数
六面体-空气为第二相-空气的体积分数 六面体-空气为第二相-空气的体积分数