第4章 多相流体温度分布计算
多相流基础第四章 可变形颗粒动力学
4、雾化射程 水平方向喷射时，喷雾液滴丧失动能时所能达到的平面与喷 口距离 射程较短：雾化角大，雾化粒度很细； 射程较长：密集喷雾炬，吸入空气量较少。 5、流量密度分布 在单位时间内，通过与喷射方向相垂直的单位横截面上液体 质量（体积）径向分布
燃料分布特性
煤燃烧国家重点实验室 SKLCC
§4－4 液滴的蒸发
直至再分裂成细小液滴；
4）扰动的增长速度取决于韦伯数和雷诺数。
煤燃烧国家重点实验室 SKLCC
3、各种雾化方式及装置中的各种力
1）压力雾化 强制液体通过一个小孔以达到液体的雾化。
主要力：通过小孔的压力，液体的表面张力，液体粘性力；
流量取决于压降； 主要的无量纲参数为韦伯数和雷诺数。
2）离心力压力雾化
T
5）对一个直径为d的大液滴的分裂，取L=d, ρ=ρg, U=UgUL, 则， 若We>12，则液滴将分裂。
煤燃烧国家重点实验室 SKLCC
2、液体系统的外形
1）液体的雾化取决于液体的形状； 2）实用中应以最小的代价使液体在最大的面积上分裂并
散开；
3）一种可改善雾化过程性能和雾化效果的方法是采用扇 形喷嘴以形成液体薄片，在有限扰动时分裂成细条，
CD f
24 2 g 3l k 24 Hc Re 3 g 3l k Re
k为阻滞效率引起的系数，取决于环流型式及Re数
对于非球形液滴，CD 应取实际液滴形状的系数值；
对于大液滴，在下落过程中液滴将发生改变：
大液滴 —— 扁平液滴 —— 分裂成小液滴
煤燃烧国家重点实验室 SKLCC
β= -b2g/σT, θ为点(x,z)处的倾斜角。
β决定液顶的形状；
气液两相流 整理
第一章概论相的概念：相是体系中具有相同化学组成和物理性质的一部分，与体系的其它均匀部分有界面隔开两相流动的处理方法：双流体瞬态模拟方法和精确描述物理现象的稳态机理模型是多相管流研究的主要方法目前研究存在的问题：1、多相流问题未得到解析解；2、油气水三相流的研究不够深入；3、水平井段变质量流动研究较少；4、缺乏向下流动的综合机理模型；5、缺乏专用研究仪器气液两相流的分类：1、细分散体系：细小的液滴或气泡均匀分散在连续相中2、粗分散体系：较大的气泡或液滴分散在连续相中3、混合流动型：两相均非连续相4、分层流动：两相均为连续相气液两相流的基本特征：1、体系中存在相界面：两相之间也存在力的作用，出现质量和能量的交换时伴随着机械能的损失2、两相的分布情况多种多样：两相流动中两相介质的分布称为流型3、两相流动中存在滑脱现象：相间速度的差异称为滑脱，滑脱将产生附加的能量损失4、沿程流体体积流量有很大变化，质量流量不变气液两相流研究方法：1、经验方法：从气液两相流动的物理概念出发，或者使用因次分析法，或者根据流动的基本微分方程式，得到反映某一特定的两相流动过程的一些无因次参数，然后依据实验数据整理出描述这一流动过程的经验关系式。
优点：使用方便，在一定条件下能取得好的结果缺点：使用有局限性，且很难从其中得出更深层次的关系2、半经验方法：根据所研究的气液两相流动过程的特点，采用适当的假设和简化，再从两相流动的基本方程式出发，求得描述这一流动过程的函数关系式，最后用实验方法确定出函数关系式中的经验系数。
优点：有一定的理论基础，应用广泛缺点：存在简化和假设，具有不准确性3、理论分析方法：针对各种流动过程的特点，应用流体力学方法对其流动特性进行分析，进而建立起描述这一流动过程的解析关系式。
优点：以理论分析为基础，可以得到解析关系式缺点：建立关系式困难，求解复杂研究气液两相流应考虑的几个问题：1、不能简单地用层流或紊流来描述气液两相流2、水平或倾斜流动是轴不对称的3、由于相界面的存在增加了研究的复杂性4、总能量方程中应考虑与表面形成的能量问题5、多相流动中各相的温度、组分的浓度都不是均匀的，相之间有传热和传质6、各相流速不同，出现滑脱问题，是多相流研究的核心与重点流动型态：相流动中两相介质的分布状况称为流型或两相流动结构流型图：描述流型变化及其界限的图。
化工原理第四章第二节讲稿
温度梯度:
温度梯度是向量，其方向垂直于等温面，并以温度增加的 方向为正。
傅立叶定律是热传导的基本定律，它指出：单位时间内传导的热量与温度梯度及垂直于热流方向的截面积成正比，即
导热系数表征物质导热能力的大小，是物质的物理性质之一，其值与物质的组成、结构、密度、温度及压强有关。
3、液体的导热系数 在非金属液体中，水的导热系数最大。除水和甘油外， 绝大多数液体的导热系数随温度的升高而略有减小，
纯液体的导热系数比溶液的导热系数大。
3、气体的导热系数 气体的导热系数很小，不利于导热，但有利于保温。 气体的导热系数随温度升高而加大 。 在相当大的压强范围内，气体的导热系数随压强变化极小 注意：在传热过程中，物质内不同位置的温度可能不相同，因而导热系数也不同，在工程计算中常取导热系数的算术平均值。
（4-7）
四、圆筒壁的稳定热传导
1、单层圆筒壁的热传导
仿照平壁热传导公式，通过该圆筒壁的导热速率可以表示为：
分离变量积分：
——圆筒壁的导热热阻
这个式子也可以写成与平壁传导速率方程类似的形式
（4-9）
精品课件！
精品课件！
——圆筒壁的内外表面的对数平均面积，m2
当r2/r1≤2时可用算术平均值代替对数平均值
1、导热系数的定义
在数值上等于单位温度梯度下的热通量 ，是物质的物理性质之一 。 一般，金属的导热系数最大，非金属的固体次之，液体的较小，气体的最小。
2、固体的导热系数 纯金属的导热系数一般随温度的升高而降低， 金属的导热系数大都随纯度的增加而增大。 非金属的建筑材料或绝热材料的导热系数随密度增加而增大，也随温度升高而增大。
1、温度场和温度梯度
一、傅里叶定律
化工流体流动与传热4.3 对流传热概述
换热器任一截 面上热流体的 平均温度
换热器任一截 面上与热流体 相接触一侧的 壁温
17
2. 热边界层
λ dt λ dt ( )w = − ( )w 因此有 α = − T − Tw dy ∆t d y
上式为对流传热系数的另一定义式， 上式为对流传热系数的另一定义式，该式表 对于一定的流体和温度差， 明，对于一定的流体和温度差，只要知道壁面附 近的流体层的温度梯度， 近的流体层的温度梯度，就可由该式求得α。 热边界层的厚薄影响层内的温度分布， 热边界层的厚薄影响层内的温度分布，因而 影响温度梯度。当边界层内、 影响温度梯度。当边界层内、外侧的温度差一定 热边界层愈薄， 愈大， 时，热边界层愈薄，则(dt/dy)w愈大，因而α就 愈大。反之，则相反。 愈大。反之，则相反。
24
4.3.3 保温层的临界直径
dc
图4-15 保温层的临界直径
25
第 4 章 传热
4.1 概述 4.2 热传导 4.3 对流传热概述 4.4 传热过程计算
4.4.1 热量衡算
26
热平衡方程
假设换热器的热损失可忽略， 假设换热器的热损失可忽略 ， 则单位时间 内热流体放出的热量等于冷流体吸收的热量。 内热流体放出的热量等于冷流体吸收的热量。 对于换热器的微元面积d 对于换热器的微元面积 dS ， 其热量衡算式 可表示为
dQ = α i (T − Tw )dS i =
λ
b
(Tw − t w )dS m = α o (t w − t )dS o
或
T − Tw Tw − t w tw − t dQ = = = 1 b 1
α i dSi
λ dS m
α o dS o
dQ = K (T − t )dS
第四章 两(多)相流动4.4
边界条件
y 0, u x 0 dux y=h, dx
2 g uR
2
求解结果
2 g uR g g 2 u x hg y gy 2 2 气流方向向下时，取正号 液体在边界上（y=h）的速度ui为 2 g uR h g 2 ui gh 2 2
4.4 液膜流动
薄层液体沿固体壁面 波动强化了传热、传质
一、液膜流动特性
1．液膜流动的基本状态
Re<20~30 层流 Re>30~50 产生了波动 Re=250~500 层流转变为湍流
2.液膜运动方程 可以近似认为液膜是沿平面运动，横向变 化大于纵向变化。可以沿用二维平板流动的边 界层方程。
滞留膜流量与厚度的计算对每一个区域分别求解流体动力学方程在动态弯月面区域端点处的自由面曲率进行匹配使用上面的方法求得低速抽出时的单位宽度上薄膜的液流流量在自由面较高处薄膜厚度不变它和流量的关系随着雷诺数的加大薄膜出现波动分析有三种途径
第四章 两（多）相流动
郭宁
山东轻工业学院
化学工程学院 2009.5
2 g ug ux y 2
当气体静止时 当y=0时， ux=uy=0
u x 0 y
求解过程
由连续性方程得
u x u y dy x
由液膜微元体积dV=h*1*dx得
d b un u x dy dx 0
假定膜只在流动方向x向改变曲率，则
对每一个区域分别求解流体动力学方程 在动态弯月面区域端点处的自由面曲率进行匹配 使用上面的方法，求得低速抽出时的单位宽度上 薄膜的液流流量
Q 0.944
流体的热传导性质和温度分布计算的统一化方程
流体的热传导性质和温度分布计算的统一化方程在热传导过程中，了解流体的热传导性质和温度分布非常重要。
这些性质和分布可以通过统一化方程来计算和描述。
本文将介绍流体的热传导性质和温度分布的统一化方程，并进行详细的说明。
1. 热传导性质的统一化方程热传导是指热量以分子间碰撞的方式在流体内传递，其速率与温度梯度成正比。
流体的热传导性质可以用热传导系数来描述，该系数表示单位长度和单位面积的传热速率。
热传导性质的统一化方程可以用如下的形式表示：q = -Λ * ∇T其中，q表示单位时间内通过单位面积的传热量，Λ表示热传导系数，∇T表示温度梯度。
这个方程表明热传导量与温度梯度成正比，并且方向是从温度高的地方向温度低的地方传导。
热传导系数Λ是材料的性质之一，不同的材料有不同的热传导系数。
2. 温度分布的统一化方程温度分布是指流体内各点的温度分布情况，通过计算和分析温度分布可以了解流体的传热情况和热平衡状态。
温度分布的统一化方程可以用如下的形式表示：∇·(Λ * ∇T) + Q = 0其中，∇·表示散度运算符，Λ表示热传导系数，∇T表示温度梯度，Q表示单位体积内的热源或热汇。
这个方程表示了温度分布变化的规律，左边的第一项表示热传导对温度分布的影响，右边的第二项表示热源或热汇对温度分布的影响。
3. 统一化方程的应用举例为了更好地理解和应用统一化方程，下面以一些具体的例子来说明。
例1：热传导问题考虑一个均匀材料的棒状结构，两端分别连接固体和冷却器，在固体端施加恒定的温度，冷却器端恒定的冷却温度。
假设棒状结构的热传导系数Λ为常数，温度分布可以通过统一化方程计算得出。
例2：热对流问题考虑一个流体在闭合容器内的对流传热问题，容器内的流体受到恒定的加热。
通过统一化方程可以计算出流体内各点的温度分布，进而分析热对流的特性。
例3：热辐射问题考虑一个辐射热传导问题，辐射热通量通过边界表面的辐射传热。
通过统一化方程可以计算出边界表面的温度分布和热通量分布，进而分析热辐射的传热特性。
第四章 两(多)相流动4.2
界面上的内部压力分布 ps p0 d ga 1 cos 终端速度
2 U ga / c 3
泰勒-戴维斯方程
尾流
液滴为球形re1忽略惯性力sinsinsinursinur为液滴腰线上速度的绝对值表面旋度coscos哈马德赖伯钦斯基公式hr公式sthr阻力曲线及速度曲线三液体中的气泡运动小气泡的运动d001cmre1若气泡在黏度较大的液体中运动甚至对于较大的气泡还可应用hr方程若有表面活性剂存在气泡速度可用stockes公式求取中等尺寸气泡的运动d02cmre700800气体的密度和黏度可以忽略若气泡为球形半径为a则绕流液体的速度势为sincosre48re48终端速度气泡恒速上升时阻力和升力相等上升速度大气泡的运动eo40re12气泡体积大于3cm18cm
d c
表面速度
u , s u0 sin
u0为液滴腰线上速度的绝对值
U U 1 u0 2 c d 2 1
c
表面旋度
U sin 2 3 s 2a 1
阻力
cUa cos 2 3 pc p0 ,c 2 2r 1 5 dUr cos pd p0d 2 a 1
1 c 2 1 d r a, pc 2 c 2 pd 2 d 2 r r sin a r r sin
U 2 r , 2 sin r 2
c
r a, c d 0
1 4 6 c d
0
d sin d d
2
3. 阻力曲线及速度曲线
三液体中的气泡运动
化工原理第四章两流体间传热过程的计算
2019/8/14
一、热交换的基本计算式
1、总传热方程式
QKAtm
式中 Q——传热速率，W； K——总传热系数，W/m2.℃ A——传热面积，m2
T1、T2——热流体的进、出口温度，℃； t1、t2——冷流体的进、出口温度，℃。
2019/8/14
T1 Cp1
qm1
qm2
t2
qm2
t1
Cp2
T2
qm1
对于热流体： Q qm 1C p1T 1T 2 对于冷流体： Q qm 2C p2t2t1
2019/8/14
②热焓法 由热力学得知，在等压条件下，物系与外界交换的
Q q m 1 r 1 C p 1 T s T 2 q m 2 r 2 C p 2 t s t 2
式中 Ts、ts——饱和温度（沸点 ）。
2019/8/14
20℃（l）
Q1qmCpt2t1
Q2 qmr
100℃（l）
100℃（g）
Q Q 1 Q 2 q m C p ( t s t ) q m r q m r C p t s t 2
1.0 0.9
ψ 0.8
0.7
0.6
0.5 0.0
0.1 0.2
0.3
0.4
0.5
0.6
0.7
p
四壳程对数平均温度差的
校正系数εψΔ t值
2019/8/14
0.8
井筒流体温度分布计算方法
井筒流体温度分布计算方法在多相管流压力计算中，需要油藏流体的高压物性数据，而流体的高压物性对压力和温度非常敏感，因而准确预测多相流体的温度是压力梯度计算的基础。
另外，油藏流体沿井筒向地面流动过程中，随着不断散热，其温度将不断降低，油温过低可能导致原油结蜡，因而多相流体温度的准确预测对怎样采取防蜡措施、是否增加井口加热设备等也是很重要的。
国内外对井筒流体温度分布进行了大量的工作。
早在1937年，Schlumberger 等人就提出了考虑井筒温度分布的意义。
五十年代初期，Nowak 和Bird 通过井筒温度分布曲线解释注水和注汽剖面。
Lasem 等人于1957年首先提出了计算井筒温度分布的方法。
Ramey.H.J 于1962年首先用理论模型描述了井筒中流体温度分布于井深和生产时间的关系。
Ramey.H.J 从能量守恒的观点出发，建立了计算井筒温度分布的能量守恒方程JdW dQ J g uduJ g gdZ dH l c c -=++(2-8)Ramey.H.J 利用该模型推导了向井中注入液体和气体时的温度分布公式。
当注入液体时：Azl e b aA t T b aA aZ t Z T --+++-=])([),(0(2-9)当注入气体时：A zl e c a A b t T c a A b aZ t Z T -⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛±+-++⎪⎭⎫ ⎝⎛±-+=7781)(7781),(0(2-10)式中： []Ukr t Uf r k W A c 112)(π+=Eickmeier 等人于1970年在Ramey.H.J 研究的基础上建立了一套关于注液和产液期间液体和井筒周围地层间热交换的有限差分模型。
计算过程中，将油管、套管、水泥环及地层的传热全都考虑在内。
但作者仍然只是研究单相流体的温度分布，传热计算中把流体的物性等都看作是常数。
后来，Satter 对注蒸汽是相态的变化对温度分布的影响进行了研究。
流体的多相流动和多相传热
流体的多相流动和多相传热多相流动和多相传热是流体力学和传热学中的重要研究领域，涉及到两种或多种不同相态的流体在相互作用中的行为及其传热特性。
在工程和科学领域中，多相流动和多相传热的研究对于理解和解决现实生活中的许多问题具有重要意义。
一、多相流动的基本概念多相流动是指在同一空间内同时存在着两种或多种不同相态的流体。
常见的多相流动包括气固、液固和气液两相流动。
在多相流动中，流体与固体或另一种流体之间通过界面相互作用并产生各种现象。
多相流动的特点包括相态转换、相分离、相互作用等。
1. 相态转换相态转换是指流体在不同条件下从一种相态转变为另一种相态的过程。
例如，水从液态转变为气态的过程称为蒸发，气体从气态转变为液态的过程称为冷凝。
相态转换会伴随着能量的吸收或释放，对多相流动和多相传热的研究具有重要影响。
2. 相分离相分离是指在两种或多种相态的流体中，不同相态的物质会产生分离现象。
例如，在液固两相流动中，固体颗粒可能会沉积在底部形成固体堆积。
相分离现象对于流体力学和传热学的研究与应用非常重要。
3. 相互作用相互作用是指在多相流动中不同相态的流体之间通过界面相互作用产生的现象。
例如，在两种相态的气体流动中，气泡的生成和破裂是气液两相流动的典型现象。
相互作用对于多相流动的传热特性具有重要影响。
二、多相流动的数学模型多相流动的数学模型在流体力学和传热学的研究中起着重要作用，它描述了不同相态的流体在空间和时间上的分布及其相互作用规律。
常见的多相流动模型包括欧拉模型、拉格朗日模型和欧拉-拉格朗日混合模型。
1. 欧拉模型欧拉模型假设流体处于连续介质的状态，将流体视为宏观量的集合。
欧拉模型通过质量、动量和能量守恒方程来描述流体的运动及其与固体界面的相互作用。
欧拉模型在多相流动的研究和工程应用中被广泛采用。
2. 拉格朗日模型拉格朗日模型将流体视为由大量微观粒子组成的离散系统。
拉格朗日模型通过分析流体粒子的运动轨迹和相互作用来描述流体的运动规律。
多相流数值计算
膨胀过程喷管的设计可认为其流动为一维、无粘、超声速膨胀流动，过程如图1。
蒸汽在点1处等熵膨胀，在点2处越过饱和线，此时过饱和度S=1，蒸汽为饱和蒸汽。
由入口条件T0和P0决定，点2发生在收缩断（亚音速）或者扩张段（超音速）。
蒸汽在膨胀过程中快速冷却，通常为10^6K/s，但是在点2并没有发生凝结。
当等熵膨胀到点3，均值凝结成核率急剧上升，大量的分子簇以10^19/cm^-3s^-1的速度形成。
汽化潜热的释放使蒸汽的热力学参数偏离等熵线，压力升高1%，该点称为“凝结起始点”。
也是在这一点，光散射的现象发生。
在蒸汽喷管当中，这一点的过饱和度通常在4-15，取决于起始条件。
凝结核的直径在10^-7cm数量级。
继续流动，蒸汽在这些凝结核上凝结，汽化潜热的释放使其热力学参数远离等熵线到达4点。
3点到4点的区域称为“凝结区域”。
在凝结区域最后，蒸汽的热力学状态接近饱和线。
凝结核生长的速度下降，流动继续膨胀和冷却。
通过改变入口条件，即P0和T0值，蒸汽可以沿着不同的等熵线膨胀，以致发生在点3附近的成核过程可以通过三相点温度检查出。
喷管内流动的控制方程：连续性方程：dρρ+dAA+duu=−dg1−g或者：ρAu1−g=m，动量方程：ρudu+dp=0能量方程：d +u22=dq≅Ldg状态方程：p=1−gρRT/μ我们有四个方程组，五个未知数。
求解有两种方法，一是测量一个量比如压力，二是使g为其他参数的函数。
J=C(p)2(2ς)1/2v c exp(−∆G∗/kT)J是成核速率（cm^-3s^-1），∆G∗为Gibbs自由能（∆G∗=4πr∗2ς/3），临界半径r由Gibbs-Thompson-Helmholtz方程给出。
正如很多文献所述，表面张力对成核速率J的影响非常大。
r∗=2ςv c/kT ln(p/p∞)因为凝结核的尺寸远远小于蒸汽分子平均自由程，生长规律简化了。
dr=f(p,T,p D,T D,α,ξ)方程已经闭合，喷管内任何地方所有的热力学参数可以求解，甚至包括液滴尺寸分布。
化工原理第四章对流传热
【解】在确定各物理量时，先确定定性温度。
一般情况下，用进出设备流体的温度的平均值
（算术平均值），即：
t t进＋t出 ＝20＋40＝30℃
2
2
查数据手册，30℃时水的物性数据为：
Cp=4183J/(K.kg) ρ=996kg/m3 μ=8.01×10－4Pa.s λ=0.618W/(m.K)
【注意事项】
（1）定性温度取流体进出温度的算术平均值tm； （2）特征尺寸为管内径d；
（3）流体被加热时，n＝0.4；
流体被冷却时，n＝0.3。
（4）若l/d＜60 ，进行校正：
'
1
d
0.7
l
3/24/2020
（2）圆形直管内的湍流（高粘度流体）
0.027 ( du )0.8 ( c p )0.33 ( )0.14
（1）什么是定性温度 【定义】确定物性参数 数值的温度称为定性温 度。
Re du
T1
t2
Pr c p
T2
t1
3/24/2020
（2）定性温度的取法 ①流体进、出口温度的平均值
②膜温
tm
t1
t2 2
t tm tw 2
th T1
热Φ 流 体
th,w
t2
Φ
冷 流 tc,w 体
式中 tw——壁面上的温度；
bt
Q bt A(tw t) 当流体被冷却时：
Q
bt'
A(T
Tw )
bt’
3/24/2020
4、牛顿冷却定律
令：
bt
Q
bt
A(t w
t)
流体被加热： Q A(tw t)
第四章传热—3(对流传热分析和计算)讲解
流体物理特性参数的温度称为定性温度。一般定性温
度有三种取法：进、出口流体的平均温度，壁面平均
温度，流体和壁面的平均温度（膜温）。
4 、准数是一个无因次数群，其中涉及到的物理量必 须用统一的单位制度。
二、流体无相变时对流传热系数的关联式
1、流体在圆形直管内作强制对流 （1） 低粘度流体
Nu=0.023Re0.8Prn
准数关联式是一种经验公式，在利用关联式求对流传 热系数时，不能超出实验条件范围。 在应用关联式时应注意以下几 点： 1、应用范围 2、特性尺寸 无因次准数 Nu、Re等中所包含的传热面尺寸称为特征 尺寸。通常是选取对流体流动和传热发生主要影响的 尺寸作为特征尺寸。
3、定性温度
流体在对流传热过程中温度是变化的。确定准数中
液膜两侧温差
表面粗糙度
一般换热设备中的冷凝可按膜状冷凝考虑。 冷凝的传热系数一般都很大，如水蒸汽作膜状冷 凝时的传热系数α通常为5000~15000 W/(m2· ℃ )。 因而传热壁的另一侧热阻相对地大，是传热过程 的主要矛盾。 当蒸汽中不凝性气体的含量为1%时，可降低α 60%左右。因此冷凝器应装有放气阀，以便及 时排除不凝性气体。
查25℃时水的物性数据（见附录）如下 ： Cp=4.179×103J/kg· ℃ μ=90.27×10-3N· s/m2
Re=duρ/μ=(0.02×0.4 ×997)/(90.27 ×10-5）=8836 Re在2300~10000之间，为过渡流区
Pr=cpμ/λ=(4.179 ×103 ×90.27 × 10-5)/60.8 × 10-2 =6.2 a可按式
3、流体的运动状况：
传热热阻主要集中在层流底层，湍流使滞流底层
厚度减薄，对流表面传热系数也就随之增大。
求解圆管中流体以及管壁的温度分布
求解圆管中流体以及管壁的温度分布
一，设置参数
pi=acos(-1)
ri=((5/16)-0.083)*0.0254 !管内半径
r0=(5/16)*0.0254 !管外半径
l=12*0.0254 !管长度
nthick=1 !径向单元个数
size=(r0-ri)/nthick !径向单元长度
nlength=nint(l/size) !长度方向单元个数
ktube=100 !管导热系数
kwater=0.6 !流体导热系数
denswatr=1000 !流体密度
cwater=4183 !流体比热
vinput=3 !入口流速
carea=pi*(ri**2) !流通面积
massflow=vinput*carea*denswatr !质量流率
tinlet=20 !流体入口温度
tbulk=100 !管外环境温度
hi=15000 !管内壁与流体的对流系数
h0=10000 !管外壁对流系数
二，定义三类单元
FLUID116, k1=1 , k2=1
SURF151, k3=1, k4=1,k5=1,k8=2
Plane55,k3=1
定义材料
关键点k,1,0,0,0 k,2,0,-l,0 k,3,0,-l-size,0
建直线l,1,2 l,2,3 创建组元waterlin
生成管壁
四，划分网格Lmesh,1，lmesh,2
Amesh,1
创建表面效应单元
五·，施加载荷1，流体入口温度
4，管内流体质量流率
六，求解
七，后处理
各节点温度云图。
多相流数值计算
FLUENT计算两相流相关问题：通常把含有大量固体或液体颗粒的气体或液体流淌称为两相流；两相流的讨论：对两相流的讨论有两种不同的观点：一是把流体作为连续介质而把颗粒群作为离散体系；而另一是除了把流体作为连续介质外，还把颗粒群当作拟连续介质或拟流体。
引入两种坐标系：即拉格朗日坐标和欧拉坐标，以变形前的初始坐标为自变量称为拉格朗日Langrangian坐标或物质坐标；以变形后瞬时坐标为自变量称为欧拉Eulerian坐标或空间坐标。
离散相模型FLUENT在求解连续相的输运方程的同时,在拉格朗日坐标下模拟流场中离散相的其次相；离散相模型解决的问题：煤粉燃烧、颗粒分别、喷雾干燥、液体燃料的燃烧等；应用范围：FLUENT中的离散相模型假定其次相体积分数一般说来要小于10-12% (但颗粒质量承载率可以大于10-12%,即可模拟离散相质量流率等大于连续相的流淌)；不适用于模拟在连续相中无限期悬浮的颗粒流问题，包括：搅拌釜、流化床等；颗粒颗粒之间的相互作用、颗粒体积分数对连续相的影响未考虑；湍流中颗粒处理的两种模型：Stochastic Tracking,应用随机方法来考虑瞬时湍流速度对颗粒轨道的影响；Cloud Tracking,运用统计方法来跟踪颗粒围绕某一平均轨道的湍流集中。
通过计算颗粒的系统平均运动方程得到颗粒的某个平均轨道。
多相流模型FLUENT中供应的模型：VOF 模型(VblUme ofFluid Model)混合模型(MiXtUre Model)欧拉模型(EUlerian Model)1.VOF 模型(VOIilme Of Fluid Model)VOF模型用来处理没有相互穿插的多相流问题，在处理两相流中，假设计算的每个掌握容积中第一相的体积含量为αl,假如αl=0,表示该掌握容积中不含第一相，假如αl=l,则表示该掌握容积中只含有第一相，假如O<αl<l,表示该掌握容积中有两相交界面；VOF方法是用体积率函数表示流体自由面的位置和流体所占的体积，其方法占内存小，是一种简洁而有效的方法。
气液两相流 整理
第一章概论相的概念：相是体系中具有相同化学组成和物理性质的一部分，与体系的其它均匀部分有界面隔开两相流动的处理方法：双流体瞬态模拟方法和精确描述物理现象的稳态机理模型是多相管流研究的主要方法目前研究存在的问题：1、多相流问题未得到解析解；2、油气水三相流的研究不够深入；3、水平井段变质量流动研究较少；4、缺乏向下流动的综合机理模型；5、缺乏专用研究仪器气液两相流的分类：1、细分散体系：细小的液滴或气泡均匀分散在连续相中2、粗分散体系：较大的气泡或液滴分散在连续相中3、混合流动型：两相均非连续相4、分层流动：两相均为连续相气液两相流的基本特征：1、体系中存在相界面：两相之间也存在力的作用，出现质量和能量的交换时伴随着机械能的损失2、两相的分布情况多种多样：两相流动中两相介质的分布称为流型3、两相流动中存在滑脱现象：相间速度的差异称为滑脱，滑脱将产生附加的能量损失4、沿程流体体积流量有很大变化，质量流量不变气液两相流研究方法：1、经验方法：从气液两相流动的物理概念出发，或者使用因次分析法，或者根据流动的基本微分方程式，得到反映某一特定的两相流动过程的一些无因次参数，然后依据实验数据整理出描述这一流动过程的经验关系式。
优点：使用方便，在一定条件下能取得好的结果缺点：使用有局限性，且很难从其中得出更深层次的关系2、半经验方法：根据所研究的气液两相流动过程的特点，采用适当的假设和简化，再从两相流动的基本方程式出发，求得描述这一流动过程的函数关系式，最后用实验方法确定出函数关系式中的经验系数。
优点：有一定的理论基础，应用广泛缺点：存在简化和假设，具有不准确性3、理论分析方法：针对各种流动过程的特点，应用流体力学方法对其流动特性进行分析，进而建立起描述这一流动过程的解析关系式。
优点：以理论分析为基础，可以得到解析关系式缺点：建立关系式困难，求解复杂研究气液两相流应考虑的几个问题：1、不能简单地用层流或紊流来描述气液两相流2、水平或倾斜流动是轴不对称的3、由于相界面的存在增加了研究的复杂性4、总能量方程中应考虑与表面形成的能量问题5、多相流动中各相的温度、组分的浓度都不是均匀的，相之间有传热和传质6、各相流速不同，出现滑脱问题，是多相流研究的核心与重点流动型态：相流动中两相介质的分布状况称为流型或两相流动结构流型图：描述流型变化及其界限的图。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
4 实测值与计算值比较
右图为某油井实测的温度 值与计算值的比较。从图 中看出，实测的井温值与 计算的井温值重合较好， 且与地温梯度曲线近乎平 行两条曲线对比，最大误 差为1.6℃。
实测井温曲线与计算值比较
实际应用中传热系数的处理
t
t1s
h
GC KD
1
e
KDh GC
Kl K D
t
t1s
h
联立上述公式可得 KDt t1s hdh GCdt
求解方程得
t
t1s
h
GC KD
1
KDh
e GC
此式即井筒 内任意深度 点油温计算 公式
t
t1s
h
GC KD
1
e
KDh GC
分析
第一项： t1s h 它反映地温自然变化规律，意谓着油
流静止时，原油本身的温度完全为环境地温所决定，因此
可称“静态温度”。
公式的推导 传热公式 油损失热量为
根据热量守恒
dq K t ts Ddh
K—J/(m2.s. ℃) K—kJ/(m2.h. ℃) K—kcal/(m2.h. ℃)
dq GCdt
C为原油 比热 C＝2.1kJ／kg.℃
C＝0.5kcal／kg.℃
K t ts Ddh GCdt
地层温度与深度的关系 ts t1s h
(2) 地温梯度
一般认为 =0.03℃／m
(3) 距井底高度 h
取某一点至油层中部的距离
(4) 原油质量流量
为油井的实际产量或设计产量
(5) 原油比热C
一般计算时取C＝2.1kJ／kg.℃ 当原油含水时 GC G油C油 G水C水
(6) 总传热系数K
1
K
1
1
i i
1
2
R0
(7)油管外径D
其数据在完井后即可提供
3 假设条件的分析
(1) 气体质量流量 国内多数油田的气油比一般在10～80之间，气量可忽略不 计，假设条件（1）可以成立。但对气油比大于100的井， 可将气量折换成油量。 GC G油C油 G水C水 G气C气
(2) 流动型态
井筒内的流型较复杂，不同流型的传热特性并不相同， 但目前还缺乏反映具体流型传热规律的量化模型，考虑 到只研究宏观整体，故可将流动视为准稳定流，同时忽 略流型的影响。
稠油井筒降粘工艺
热流体循环
电热杆采油工艺 伴热电缆采油工艺
流动保障 Flow Assurance
1、公式推导及分析
设 t1 为井底油温，T为井口油温， t1s 为井底处地层温度， t2s 为井 口处地层温度。
假设条件：
•气体质量忽略不计 •井筒中液体流动为准稳定流，体积 和流型变化的影响忽略不计 •流体对地层放热，其总传热系数K 为常数 •油流在油管中流动时因摩擦而产生 的热量忽略不计 •因天然气析出及膨胀吸热忽略不计
课程回顾
➢多相流计算中为什么要计算油藏流体的高压物性？
➢井筒中某点(温度和压力分别为T 和p)油相和气相实
际体积流量的计算步骤。
Qo×Bo
Zp0T pT0
Qo
(Rp
Rs
)
第四章 多相流体温度分布计算
在多相管流压力计算中，首先要计算油藏流体的高压 物性，而流体的高压物性对压力和温度非常敏感，因 而准确预测多相流体的温度是压力计算的基础。另外， 油藏流体在地层中温度高，沿井筒向地面流动过程中， 随着不断散热，其温度将不断降低，油温过低可能导 致原油结蜡，因而多相流体温度的准确预测对怎样采 取防蜡措施，是否增加井口加热设备等是很重要的。
GC Kl
1
Klh
e GC
1
Kl 1.1573 5.4246exp(
G
)
1000
K—kJ/(m2.h. ℃) K—kcal/(m2.h. ℃)
Kl—kJ/(m.h. ℃) Kl—kcal/(m.h. ℃) G — kg/h Kl —kcal/(m.h. ℃)
第二项意味着因油流运动和地层温度对井筒油温的影响， 从而产生“静态温度”的增量，即“动态温度”。如果总 传热系数的值或套管直径较大，散热情况良好，井筒油温 就低；如果油流量G或液体比热C较大，则井筒油温就高。
2 公式中各项参数取值 (1) 井底油温 t1s
t
t1s
h
GC KD
1
e
KDh GC
井底油温也就是油层温度
当压力低于饱和压力时，有天然气析出。 析出气体需要热量，已析出气体不断膨胀， 又会吸收一部分热量，这两部分热量的计 算比较复杂。
焦耳-汤普逊系数：
h
T p
h
焦耳-汤普逊效应
在井底压力大于饱和压力且气油比小于100时，气体析出 膨胀而引起的温降一般在2℃左右。而且此温降又被油气 在油管中由于摩擦所产生的热量补偿掉一部分。因此，在 一般工程计算中，可忽略此温降。
(3) 总传热系数
井筒外部不同位置的岩层性质 亦不同，井筒内的流体流型又 有变化，所以严格地说，总传 热系数值应该是一个变量，认 为其变化与平均值之间的差值 不大，取实测平均值即可[实测 在21～25kJ／(m2·h·℃)左右]。 因此，计算时可按常数考虑。
1
K
1
1
i i
1
2
R0
(4) 天然气析出和膨胀Байду номын сангаас题