油气多相流试井模型
多相流试井模型
对于复杂、多变的非均质油藏的试井问题以及多相流试井的问题，描述多相流的方程是高度非线性的，在不做较大简化的情况下，没有解析解，因此，对多相流试井问题需要用数值解法。试井问题数据解法称为数值试井。
数值试井作为一个完整的体系，其研究所涉及的内容包括四个部分：空间离散、产量描述、具体问题的一般规律和参数拟合。对于数值试井问题的研究首先必须有一套合适的网格，此网格要既能描述和反映开关井后压力的细微变化，又能兼顾计算效率，在此基础上形成有限差分方程才能更好地解决数值试井问题。
油气两相数值试井模型
油藏注入、采出的不平衡造成了油田注水区块内部的压力不均衡。油藏的非均质性的客观存在，也给稳定油藏压力、控制水驱突进造成了很大的困难。从油田普遍的生产实践看来，都无法完全保证油井生产时的井底压力维持在原油饱和压力以上，这必然造成近井地带原油脱气现象的存在。由此可见，油、气、水三相同产的现象在油田的生产实际中是普遍存在的。
多相流试井分析人们仍在沿用Perrine，Al-Khalifa及Raghavan等人提出的近似解析方法。多相流试井分析近似的方法或多或少地带有单相流试井的影子，因而蒙蔽了多相流试井曲线的一些重要特征。作为数值试井理论的应用，本节将通过对油气两相流动试井问题的分析，介绍油气两相流试井曲线的基本特征。
对于油气两相流动，只需要在第六章离散控制方程及本章1节相应的部分中去掉与水相相关的部分即可。
两相流动时，是否存在饱和度梯度对试井曲线是有很大影响的。这使得多相流动试井具有比单相流动更多的多解性。作为对油气两相试井规律的研究，本节将饱和度梯度的影响剥离开来，不考虑饱和度梯度，来讨论不同因素对曲线特征的影响。
一、油气两相流试井一般曲线特征
均质油藏中油气流动压降复合曲线如图12.5.1，此图显示出了一个有趣的现象,均质油藏油气两相流试井曲线竟然表现出了单相油流双渗介质的特征:在压力导数曲线
（14.2.1－14.2.4式）上的驼峰结束以后，曲线急转而下，出现了一个凹子，然后开始上升，最后维持在一个平稳的水平上。
根据曲线，可以将油气两相流试井可以分为三个阶段：第一阶段为井储影响阶段。此阶段压力梯度尚未建立，从地层向井底的流入非常有限，因此地层中基本保持原来的状态。这一阶段与均质单相流动也是相同的。第二阶段为过渡阶段，此时，随着地层中流体向井底流入，压力梯度逐步形成。压力的降低造成气体迅速膨胀，由于气体的压缩率比原油的压缩率大得多，大大有利于流体得产出，必然会对压力下降造成缓冲，使压力下降速率降低，从而形成凹子。当压力梯度形成后，流动进入第三阶段，是为拟稳定阶段，压力下降速度维持不变。压力恢复复合曲线如图12.5.2。
0.01
0.1
1
10
100
0.001
0.010.11101001000
logtD/CD
l o g P w D & l o g P 'w D
上凸部分
下凹部分
图12.5.1 油气两相流，压降试井，压力和压力导数复合曲线
压恢曲线的导数曲线上在驼峰结束后也有一个“凹子”。曲线从“凹子”上升到一定的阶段后，稍微下降又再度恢复至与前次上升相同的水平，保持此水平不变。曲线同样具有与双渗介质类似的特征。压恢油气两相流也可以分为三个阶段，但其变化机理基本上是一个压缩的过程：第一阶段为井储影响阶段。此时，随流体持续流入，井底开始憋压，地层中整体压力变化不太明显。第二阶段为过渡段，由于地层中压力普遍回升，气体受压缩，同样由于气体的压缩率大大高于原油的压缩率，将产生很大的储集空间，从而使压力恢复速率降低，缓冲了压力回升的过程，从而形成凹子。当由井筒而向外的地层中的压力基本回升后，压力恢复基本稳定，此为流动的第三阶段，拟稳定阶段。
0.01
0.1
1
10
100
0.001
0.010.11101001000
logtD/CD
l o g P w D & l o g P 'w D
图12.5.2 油气两相流，压恢试井，压力和压力导数复合曲线
总之，油气两相流之所以出现类似双重介质的凹子现象，主要是由于参与流动的气体压缩率很大，因而对压力变化极为敏感的缘故。当压力稍微降低时，气体膨胀，可以将原油挤出地层，使压力降低暂时减缓；当压力升高，气体压缩，为原油提供空间，抵消一部分压力升高的因素，从而起到了缓冲压力变化的作用。
图12.5.1、12.5.2生动地说明了油气两相渗流试井与单相流试井乃至油水两相流试井的不同。这是油气两相的物性差异悬殊造成的。在油水两相流情况下之所以与单相流的曲线基本相似完全是因为两者之间的差异还不足以产生任何非单相流的特征。由此可见，在有气相参与流动的情况下，用单相流度简单迭加（多相流分析压力方法）构成多相流的总流度必然是不合适的，毕竟，用解析的方法来研究单相流试井的时候是无法把流体的参数看作是试井的函数的。换句话说，由于均质油藏单相流试井曲线上并没有出现过类似于双重介质的“凹子”，又由于相对渗透率和流体粘度与时间无关（如果考虑相对渗透率和流体粘度与时间有关则用Martin 的推导方法是没有结果的），因此进行总流度和综合压缩系数的校正之后的单相流试井曲线相当于对其压力乘上了一个系数，则曲线只会在横、纵坐标上发生相对位移，曲线形态是不会发生改变的。以上对于油气两相流动的新认识，能为实际的生产试井解释提供一个指示，通过是否出现“凹子”来判断是否有气相参与流动，当有“凹子”出现时，进行解释时就应充分考虑气体的影响。
二、饱和度的影响
饱和度对产量规律的影响在于两个方面：一方面是气体的膨胀与压缩，当气体的饱和度增加后，膨胀与压缩量将大大增加；另一方面，油气两相渗透率受两相饱和度的影响相当大，如果在油藏的局部地区发生了脱气现象，受气相饱和度增加的影响，相对渗
透率可能发生数量级的变化。进而，气、油流度比相对于气相饱和度的关系在半对数坐标上基本上呈一条直线（图12.5.3），说明，气、油流度比与气相饱和度近似呈指数关系。因此，气体饱和度越大，流度差异越大，气体越流现象越严重，近井带越易富集气体，从而更进一步加大气体的膨胀与压缩量，缩短了井储结束时间。同时，当油相流度随油相饱和度减少而减小后，也使在开、关井后的压力降落与恢复的过程中压力变化的幅度提高、速度加快。
110100
1000
0.3
0.4
0.5
0.6
Sg
M o g
图12.5.3 气、油流度比与气相饱和度关系
0.01
0.1
1
10
100
0.001
0.010.1
110100
logtD/CD
l o g P w D & l o g P 'w D
0.7/0.3
0.6/0.4
0.5/0.50.45/0.55
图12.5.4 油气两相压降试井曲线
从压降复合曲线（图12.5.4）可以看到，随气相饱和度增加，压力迅速大幅度的下降掩盖了早期井筒反应的特征，能很清楚地看出来。受饱和度的影响，试井曲线有以下规律：
1. 随气相渗透率的增加、油相渗透率的减小，压力和压力导数曲线对往左上角移动。
2. 但曲线总的形态并未改变。压力曲线，先有一段时期的快速上升，而后曲线又显现
出稳中有升的趋势。压力导数的曲线，首先上升大到一个“驼峰”，随后曲线迅疾下降至“峰谷”（凹子），曲线从“峰谷”以较大的速度上升到一定程度之后，曲线再
呈现出平稳上升的走势。这是它们之间的共性。
0.01
0.1
1
10100
0.001
0.010.1
110100
logtD/CD
l o g P w D & l o g P 'w D
0.7/0.3
0.6/0.
0.5/0.5
0.45/0.55
图12.5.5 油气两相压恢试井曲线对比
3. 对于不同的初始饱和度，压力曲线平缓上升段与压力导数曲线平缓上升段之间的相
对距离基本保持不变。
4. 对于不同的饱和度初始分布，“峰顶”和“峰谷”之间的垂向相对距离是不同的。且
随着油相饱和度的减小，这个距离有加大的趋势。即气相饱和度越大，“凹子”越深。注意到在双渗介质中，定义了一个地层系数比，2
2111
1h k h k h k +=
κ。这个系数表征不同
地层渗透率差异的大小，κ值越大，表明层间差异越大，在双渗介质复合图版上表现出来的“凹子”也就越深。可见，κ值与气相饱和度之间具有一定的对应性。 5. 油相饱和度越大，曲线受井筒储存影响的时间越长；油相饱和度越小，曲线受井筒
储存影响的时间越短，径向拟稳态流动开始的时间越早，曲线开始平稳上升的时间也就越早。
图12.5.5乃压力恢复复合曲线，其受饱和度影响的规律与压降情况下是相同的。
三、油气粘度差异的影响
受油相粘度的增加，气、油流度比为气相饱和度的线性函数，变化不大。重要的是，气体的膨胀与压缩对试井过程的作用程度相同。图12.5.6是受粘度差异影响的油水试井压降复合曲线。
除在过渡段（图12.5.6两虚线之间的曲线部分）曲线存在差异以外，油相粘度的变化对试井曲线没有影响。而在过渡段，当油相粘度较小时，曲线产生的凹子浅，随粘度增加，凹子加深。之所以出现凹子加深，是因为随油粘度增加，近井带气体的越流程
度增加，饱和度相应较大，所能产生的缓冲相应就会大一些。但越流现象毕竟只局限于近井地带，并不影响远井处，因此在拟稳定阶段各曲线基本重合。
0.01
0.1
1
10100
0.001
0.010.1
110100
log tD/CD
l o g P w D & l o g P 'w D
图12.5.6 粘度对压降复合曲线的影响
与油水两相流动又一次不同的是，油气初始饱和度相同时，只有在相同的油气粘度对比下，压力降落与恢复复合曲线才是重合的。这再次显示了多相流动的复杂性。
Fluent多相流模型选择
FLUENT多相流模型 分类 1、气液或液液流动 气泡流动：连续流体中存在离散的气泡或液泡 液滴流动：连续相为气相，其它相为液滴 栓塞（泡状）流动：在连续流体中存在尺寸较大的气泡 分层自由流动：由明显的分界面隔开的非混合流体流动。 2、气固两相流动 粒子负载流动：连续气体流动中有离散的固体粒子 气力输运：流动模式依赖，如固体载荷、雷诺数和例子属性等。最典型的模式有沙子的流动，泥浆流，填充床以及各相同性流 流化床：有一个盛有粒子的竖直圆筒构成，气体从一个分散器进入筒内，从床底不断冲入的气体使得颗粒得以悬浮。 3、液固两相流动 泥浆流：流体中的大量颗粒流动。颗粒的stokes数通常小于1。大于1是成为流化了的液固流动。 水力运输：在连续流体中密布着固体颗粒 沉降运动：在有一定高度的盛有液体的容器内，初始时刻均匀散布着颗粒物质，随后，流体会出现分层。 4、三相流 以上各种情况的组合 多相流动系统的实例 气泡流：抽吸、通风、空气泵、气穴、蒸发、浮选、洗刷。 液滴流：抽吸、喷雾、燃烧室、低温泵、干燥机、蒸发、气冷、洗刷。 栓塞流：管道或容器中有大尺度气泡的流动 分层流：分离器中的晃动、核反应装置沸腾和冷凝 粒子负载流：旋风分离器、空气分类器、洗尘器、环境尘埃流动 气力输运：水泥、谷粒和金属粉末的输运 流化床：流化床反应器、循环流化床 泥浆流：泥浆输运、矿物处理 水力输运：矿物处理、生物医学、物理化学中的流体系统 沉降流动：矿物处理。 多相流模型的选择原则 1、基本原则
1)对于体积分数小于10%的气泡、液滴和粒子负载流动，采用离散相 模型。 2)对于离散相混合物或者单独的离散相体积率超出10%的气泡、液滴 和粒子负载流动，采用混合模型或欧拉模型。 3)对于栓塞流、泡状流，采用VOF模型 4)对于分层/自由面流动，采用VOF模型 5)对于气动输运，均匀流动采用混合模型，粒子流采用欧拉模型。 6)对于流化床，采用欧拉模型 7)泥浆和水力输运，采用混合模型或欧拉模型。 8)沉降采用欧拉模型 9)对于更一般的，同时包含多种多相流模式的情况，应根据最感兴趣 的流动特种，选择合适的流动模型。此时由于模型只是对部分流动特 征采用了较好的模拟，其精度必然低于只包含单个模式的流动。 2、混合模型和欧拉模型的选择原则 VOF模型适合于分层的或自由表面流，而混合模型和欧拉模型适合于流动中有相混合或分离，或者分散相的体积分数超过10%的情况（小于10%可使用离散相模型）。 1)如果分散相有宽广的分布（如颗粒的尺寸分布很宽），最好采用混 合模型，反之使用欧拉模型。 2)如果相间曳力规律一直，欧拉模型通常比混合模型更精确；若相间 曳力规律不明确，最好选用混合模型。 3)如果希望减小计算了，最好选用混合模型，它比欧拉模型少解一部 分方程；如果要求精度而不在意计算量，欧拉模型可能是更好的选择。 但是要注意，复杂的欧拉模型比混合模型的稳定性差，可能会遇到收 敛困难。
(推荐)FLUENT中两相流、多相流中模型的的选择问题
两相流：通常把含有大量固体或液体颗粒的气体或液体流动称为两相流；其中含有多种尺寸组颗粒群为一个“相”，气体或液体为另一“相”，由此就有气—液，气—固，液—固等两相流之分。 两相流的研究：对两相流的研究有两种不同的观点：一是把流体作为连续介质，而把颗粒群作为离散体系；而另一是除了把流体作为连续介质外，还把颗粒群当作拟连续介质或拟流体。 引入两种坐标系：即拉格朗日坐标和欧拉坐标，以变形前的初始坐标为自变量称为拉格朗日Langrangian 坐标或物质坐标；以变形后瞬时坐标为自变量称为欧拉Eulerian 坐标或空间坐标。 一.离散相模型 FLUENT在求解连续相的输运方程的同时，在拉格朗日坐标下模拟流场中离散相的第二相； 离散相模型解决的问题：煤粉燃烧、颗粒分离、喷雾干燥、液体燃料的燃烧等； 应用范围：FLUENT中的离散相模型假定第二相体积分数一般说来要小于10-12%（但颗粒质量承载率可以大于10-12%，即可模拟离散相质量流率等/大于连续相的流动）；不适用于模拟在连续相中无限期悬浮的颗粒流问题，包括：搅拌釜、流化床等； 颗粒-颗粒之间的相互作用、颗粒体积分数对连续相的影响未考虑； 湍流中颗粒处理的两种模型：Stochastic Tracking，应用随机方法来考虑瞬时湍流速度对颗粒轨道的影响；Cloud Tracking，运用统计方法来跟踪颗粒围绕某一平均轨道的湍流扩散。通过计算颗粒的系统平均运动方程得到颗粒的某个“平均轨道” 二.多相流模型 FLUENT中提供的模型： VOF模型(Volume of Fluid Model) 混合模型(Mixture Model) 欧拉模型(Eulerian Model) 1.VOF模型(Volume of Fluid Model) VOF模型用来处理没有相互穿插的多相流问题，在处理两相流中，假设计算的每个控制容积中第一相的体积含量为α1，如果α1=0，表示该控制容积中不含第一相，如果α1=1，则表示该控制容积中只含有第一相，如果0<α1<1，表示该控制容积中有两相交界面； VOF方法是用体积率函数表示流体自由面的位置和流体所占的体积，其方法占内存小，是一种简单而有效的方法。 2.混合模型(Mixture Model) 用混合特性参数描述的两相流场的场方程组称为混合模型； 考虑了界面传递特性以及两相间的扩散作用和脉动作用；使用了滑移速度的概念，允许相以不同的速度运动； 用于模拟各相有不同速度的多相流；也用于模拟有强烈耦合的各向同性多相流和各相以相同速度运动的多相流； 缺点：界面特性包括不全，扩散和脉动特性难于处理。 3.欧拉模型(Eulerian Model) 欧拉模型指的是欧拉—欧拉模型； 把颗粒和气体看成两种流体，空间各点都有这两种流体各自不同的速度、温度和密度，这些流体其存在在同一空间并相互渗透，但各有不同的体积分数，相互间有滑移；
多相流模拟知识讲解
多相流模拟
多相流模拟介绍 自然界和工程问题中会遇到大量的多相流动。物质一般具有气态、液态和固态三相，但是多相流系统中相的概念具有更为广泛的意义。在多项流动中，所谓的“相”可以定义为具有相同类别的物质，该类物质在所处的流动中具有特定的惯性响应并与流场相互作用。比如说，相同材料的固体物质颗粒如果具有不同尺寸，就可以把它们看成不同的相，因为相同尺寸粒子的集合对流场有相似的动力学响应。本章大致介绍一下Fluent中的多相流建模。 多相流动模式 我们可以根据下面的原则对多相流分成四类： ?气-液或者液-液两相流： o气泡流动：连续流体中的气泡或者液泡。 o液滴流动：连续气体中的离散流体液滴。 o活塞流动：在连续流体中的大的气泡 o分层自由面流动：由明显的分界面隔开的非混合流体流动。 ?气-固两相流： o充满粒子的流动：连续气体流动中有离散的固体粒子。 o气动输运：流动模式依赖诸如固体载荷、雷诺数和粒子属性等因素。最典型的模式有沙子的流动，泥浆流，填充床，以及各向同性流。 o流化床：由一个盛有粒子的竖直圆筒构成，气体从一个分散器导入筒内。从床底不断充入的气体使得颗粒得以悬浮。改变气体的流量，就会有气泡不断的出 现并穿过整个容器，从而使得颗粒在床内得到充分混合。 ?液-固两相流
o泥浆流：流体中的颗粒输运。液-固两相流的基本特征不同于液体中固体颗粒的流动。在泥浆流中，Stokes数通常小于1。当Stokes数大于1时，流动成为 流化（fluidization）了的液-固流动。 o水力运输：在连续流体中密布着固体颗粒 o沉降运动：在有一定高度的成有液体的容器内，初始时刻均匀散布着颗粒物质。随后，流体将会分层，在容器底部因为颗粒的不断沉降并堆积形成了淤积 层，在顶部出现了澄清层，里面没有颗粒物质，在中间则是沉降层，那里的粒 子仍然在沉降。在澄清层和沉降层中间，是一个清晰可辨的交界面。 三相流 (上面各种情况的组合) 多相系统的例子 ?气泡流例子：抽吸，通风，空气泵，气穴，蒸发，浮选，洗刷 ?液滴流例子：抽吸，喷雾，燃烧室，低温泵，干燥机，蒸发，气冷，刷洗 ?活塞流例子：管道或容器内有大尺度气泡的流动 ?分层自由面流动例子：分离器中的晃动，核反应装置中的沸腾和冷凝 ?粒子负载流动例子：旋风分离器，空气分类器，洗尘器，环境尘埃流动 ?风力输运例子：水泥、谷粒和金属粉末的输运 ?流化床例子：流化床反应器，循环流化床 ?泥浆流例子: 泥浆输运，矿物处理 ?水力输运例子：矿物处理，生物医学及物理化学中的流体系统 ?沉降例子：矿物处理 多相建模方法 计算流体力学的进展为深入了解多相流动提供了基础。目前有两种数值计算的方法处理多相流：欧拉－拉格朗日方法和欧拉－欧拉方法。 欧拉－拉格朗日方法
多相流模型和离散相模型的区别
多相流模型和离散相模型的区别 2008-03-30 10:18 两相流：通常把含有大量固体或液体颗粒的气体或液体流动称为两相流；其中含有多种尺寸组颗粒群为一个“相”，气体或液体为另一“相”，由此就有气—液，气—固，液—固等两相流之分。 两相流的研究：对两相流的研究有两种不同的观点：一是把流体作为连续介质，而把颗粒群作为离散体系；而另一是除了把流体作为连续介质外，还把颗粒群当作拟连续介质或拟流体。 引入两种坐标系：即拉格朗日坐标和欧拉坐标，以变形前的初始坐标为自变量称为拉格朗日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方法是用体积率函数表示流体自由面的位置和流体所占的体积，其方法占内存小，是一种简单而有效的方法。 混合模型(Mixture Model) 用混合特性参数描述的两相流场的场方程组称为混合模型；← ←考虑了界面传递特性以及两相间的扩散作用和脉动作用；使用了滑移速度的概念，允许相以不同的速度运动； ←用于模拟各相有不同速度的多相流；也用于模拟有强烈耦合的各向同性多相流和各相以相同速度运动的多相流； ←缺点：界面特性包括不全，扩散和脉动特性难于处理。
第20章 通用多相流模型--60页 多相流数据后处理
20.通用多相流模型（General Multiphase Models） 本章讨论了在FLUENT中可用的通用的多相流模型。第18章提供了多相流模型的简要介绍。第19章讨论了Lagrangian离散相模型，第21章讲述了FLUENT中的凝固和熔化模型。20.1选择通用多相流模型（Choosing a General Multiphase Model） 20.2VOF模型(Volume of Fluid(VOF)Model) 20.3混合模型(Mixture Model) 20.4欧拉模型(Eulerian Model) 20.5气穴影响(Cavity Effects) 20.6设置通用多相流问题(Setting Up a General Multiphase Problem) 20.7通用多相流问题求解策略(Solution Strategies for General Multiphase Problems) 20.8通用多相流问题后处理(Postprocessing for General Multiphase Problems) 20.1选择通用的多相流模型（Choosing a General Multiphase Model） 正如在Section 18.4中讨论过的，VOF模型适合于分层的或自由表面流，而mixture和Eulerian 模型适合于流动中有相混合或分离，或者分散相的volume fraction超过10%的情形。（流动中分散相的volume fraction小于或等于10%时可使用第19章讨论过的离散相模型）。 为了在mixture模型和Eulerian模型之间作出选择，除了Section18.4中详细的指导外，你还应考虑以下几点： ★如果分散相有着宽广的分布，mixture模型是最可取的。如果分散相只集中在区域的一部分，你应当使用Eulerian模型。 ★如果应用于你的系统的相间曳力规律是可利用的（either within FLUENT or through a user-defined function），Eulerian模型通常比mixture模型能给出更精确 的结果。如果相间的曳力规律不知道或者它们应用于你的系统是有疑问的， mixture模型可能是更好的选择。 ★如果你想解一个需要计算付出较少的简单的问题，mixture模型可能是更好的选择，因为它比Eulerian模型要少解一部分方程。如果精度比计算付出更重要， Eulerian模型是更好的选择。但是请记住，复杂的Eulerian模型比mixture模型 的计算稳定性要差。 三种模型概要的讲述，包括它们各自的局限，在Sections20.1.1，20.1.2，20.1.3中给出。 三种模型详细的讲述在Sections20.2,20.3和20.4中给出。 20.1.1VOF模型的概述及局限(Overview and Limitations of the VOF Model) 概述（Overview） VOF模型通过求解单独的动量方程和处理穿过区域的每一流体的volume fraction来模拟两种或三种不能混合的流体。典型的应用包括预测， jet breakup、流体中大泡的运动（the motion of large bubbles in a liquid）、the motion of liquid after a dam break和气液界面的稳态和瞬态处理（the steady or transient tracking of any liquid-gas interface）。 局限（limitations） 下面的一些限制应用于FLUENT中的VOF模型： ★你必须使用segregated solver. VOF 模型不能用于coupled solvers. ★所有的控制容积必须充满单一流体相或者相的联合；VOF模型不允许在那些空的区域中没有任何类型的流体存在。 ★只有一相是可压缩的。
Fluent多相流模型选择及设定
1.多相流动模式 我们可以根据下面的原则对多相流分成四类： ?气-液或者液-液两相流： o 气泡流动：连续流体中的气泡或者液泡。 o 液滴流动：连续气体中的离散流体液滴。 o 活塞流动: 在连续流体中的大的气泡 o 分层自由面流动：由明显的分界面隔开的非混合流体流动。 ?气-固两相流： o 充满粒子的流动：连续气体流动中有离散的固体粒子。 o 气动输运：流动模式依赖诸如固体载荷、雷诺数和粒子属性等因素。最典型的模式有沙子的流动，泥浆流，填充床，以及各向同性流。 o 流化床：由一个盛有粒子的竖直圆筒构成，气体从一个分散器导入筒内。从床底不断充入的气体使得颗粒得以悬浮。改变气体的流量，就会有气泡不断 的出现并穿过整个容器，从而使得颗粒在床内得到充分混合。 ?液-固两相流 o 泥浆流：流体中的颗粒输运。液-固两相流的基本特征不同于液体中固体颗粒的流动。在泥浆流中，Stokes 数通常小于1。当Stokes数大于1 时，流动成为流化（fluidization）了的液-固流动。 o 水力运输: 在连续流体中密布着固体颗粒 o 沉降运动: 在有一定高度的成有液体的容器内，初始时刻均匀散布着颗粒物质。随后，流体将会分层，在容器底部因为颗粒的不断沉降并堆积形成了淤 积层，在顶部出现了澄清层，里面没有颗粒物质，在中间则是沉降层，那里 的粒子仍然在沉降。在澄清层和沉降层中间，是一个清晰可辨的交界面。 ?三相流 (上面各种情况的组合) 各流动模式对应的例子如下： ?气泡流例子：抽吸，通风，空气泵，气穴，蒸发，浮选，洗刷 ?液滴流例子：抽吸，喷雾，燃烧室，低温泵，干燥机，蒸发，气冷，刷洗?活塞流例子：管道或容器内有大尺度气泡的流动 ?分层自由面流动例子：分离器中的晃动，核反应装置中的沸腾和冷凝 ?粒子负载流动例子：旋风分离器，空气分类器，洗尘器，环境尘埃流动 ?风力输运例子：水泥、谷粒和金属粉末的输运 ?流化床例子：流化床反应器，循环流化床 ?泥浆流例子: 泥浆输运，矿物处理 ?水力输运例子：矿物处理，生物医学及物理化学中的流体系统 ?沉降例子：矿物处理 2. 多相流模型 FLUENT中描述两相流的两种方法:欧拉一欧拉法和欧拉一拉格朗日法，后面
如何在fluent中设置多相流
3 设置一般的多相流问题（Setting Up a General Multiphase Problem） 设置和求解一般多相流问题的步骤的要点如下，各个子部分详细的讲述在随后的章节中。记住这里给出的仅是与一般多相流计算相关的步骤。有关你使用的其它模型和相关的多相流模型的输入的详细信息，将在这些模型中合适的部分给出。 1)选中你想要使用的多相流模型（VOF, mixture, or Eulerian）并指定相数。 Define Models Multiphase... 2)从材料库中复制描述每相的材料。 Define Materials... 如果你使用的材料在库中没有，应创建一种新材料。 !!如果你的模型中含有微粒（granular）相，你必须在fluid materials category中为它创建新材料（not the solid materials category.） 3)定义相，指定相间的相互作用（interaction）(例如，使用欧拉模型时的drag functions) Define Phases... 4)(仅对欧拉模型)如果流动是紊流，定义多相紊流模型。 Define Models Viscous... 5)如果体积力存在，turn on gravity and specify the gravitational acceleration. Define Operating Conditions... 6)指定边界条件，包括第二相体积份额在流动边界和壁面上的接触角。 Define Boundary Conditions... 7)设置模拟具体的解参数 Solve Controls Solution... 8)初始化解和为第二相设定初始体积份额。 Solve Initialize Patch... 9)计算求解和检查结果 *欧拉多相流模拟的附加指南（Additional Guidelines for Eulerian Multiphase Simulations）一旦你决定了欧拉多相流模型适合你的问题，你应当考虑求解你的多相流问题的需求计算能力。要求的计算能力很强的依赖于所求解的输运方程的个数和耦合程度。对欧拉多相流模型，有大数量的高度耦合的输运方程，计算的耗费将很高，在设置你的问题前，尽可能减少问题的statement到最简化的可能形式。 在你开始第一次求解尝试，取而代之尽力去求解多相流动的所有的复杂方面，你可以以简单近似地开始并且知道问题定义的最终形式。简化多相流问题的一些建议列举如下： 1．使用六面体或四边形网格（而不用四面体或三角形网格）。 2．减少相的数目。 你会发现即使简单的近似也会给你的问题提供有用的信息。 3.2选用多相流模型并指定相数（Enabling the Multiphase Model and Specifying the Number of Phases） 为了选VOF, mixture, Eulerian多相流模型，在Multiphase Model panel下选Volume of Fluid, Mixture, or Eulerian as the Model。 Define Models Multiphase... 如果你选的欧拉模型，输入如下： ?number of phases：为了给多相流计算指定相数，在Number of Phases下输入合适的值。你最多可以指定20相。 ?(optional) cavitation effects：包含气穴影响（Including Cavitation Effects） 对混合的欧拉模型计算，包含气穴影响是可能的。为了选气穴模型，在Multiphase Model panel中Interphase
在fluent中设置多相流
如何在fluent中设置多相流
————————————————————————————————作者：————————————————————————————————日期：
3 设置一般的多相流问题（Setting Up a General Multiphase Problem） 3.1使用一般多相流模型的步骤（Steps for Using the General Multiphase Models） 设置和求解一般多相流问题的步骤的要点如下，各个子部分详细的讲述在随后的章节中。记住这里给出的仅是与一般多相流计算相关的步骤。有关你使用的其它模型和相关的多相流模型的输入的详细信息，将在这些模型中合适的部分给出。 1)选中你想要使用的多相流模型（VOF, mixture, or Eulerian）并指定相数。 Define Models Multiphase... 2)从材料库中复制描述每相的材料。 Define Materials... 如果你使用的材料在库中没有，应创建一种新材料。 !!如果你的模型中含有微粒（granular）相，你必须在fluid materials category中为它创建新材料（not the solid materials category.） 3)定义相，指定相间的相互作用（interaction）(例如，使用欧拉模型时的drag functions) Define Phases... 4)(仅对欧拉模型)如果流动是紊流，定义多相紊流模型。 Define Models Viscous... 5)如果体积力存在，turn on gravity and specify the gravitational acceleration. Define Operating Conditions... 6)指定边界条件，包括第二相体积份额在流动边界和壁面上的接触角。 Define Boundary Conditions... 7)设置模拟具体的解参数 Solve Controls Solution... 8)初始化解和为第二相设定初始体积份额。 Solve Initialize Patch... 9)计算求解和检查结果 *欧拉多相流模拟的附加指南（Additional Guidelines for Eulerian Multiphase Simulations）一旦你决定了欧拉多相流模型适合你的问题，你应当考虑求解你的多相流问题的需求计算能力。要求的计算能力很强的依赖于所求解的输运方程的个数和耦合程度。对欧拉多相流模型，有大数量的高度耦合的输运方程，计算的耗费将很高，在设置你的问题前，尽可能减少问题的statement到最简化的可能形式。 在你开始第一次求解尝试，取而代之尽力去求解多相流动的所有的复杂方面，你可以以简单近似地开始并且知道问题定义的最终形式。简化多相流问题的一些建议列举如下： 1．使用六面体或四边形网格（而不用四面体或三角形网格）。 2．减少相的数目。 你会发现即使简单的近似也会给你的问题提供有用的信息。 3.2选用多相流模型并指定相数（Enabling the Multiphase Model and Specifying the Number of Phases） 为了选VOF, mixture, Eulerian多相流模型，在Multiphase Model panel下选Volume of Fluid, Mixture, or Eulerian as the Model。 Define Models Multiphase... 如果你选的欧拉模型，输入如下： ?number of phases：为了给多相流计算指定相数，在Number of Phases下输入合适的值。你最多可以指定20相。?(optional) cavitation effects：包含气穴影响（Including Cavitation Effects） 对混合的欧拉模型计算，包含气穴影响是可能的。为了选气穴模型，在Multiphase Model panel中Interphase Mass Transfer下打开Cavitation。 由于气穴影响，接下来你应指定在使用传质计算时的两个参数。这些参数的指定应当于调查下的流动特征参数
多相流理论建模
多相流理论在石油化工领域的应用研究综述 摘要:多相管道流动广泛存在于石油化工行业中，但目前国内外研究学者对多相管流机理问题在实验和理论方面尚存在有不一致的结论和认识，通过阅读相关领域文献，本文对多相管流研究进行调研整理,阐述了多相管流发展历程及研究现况，介绍了近些年来的研究热点，最后概述多相管流在石油化工行业的应用。关键词：多相管道，石油化工，研究热点，应用。 1.引言 多相管流是研究相异组份或者相异状态的物质共同存在于管流中的动力学及传热传质规律，广泛存在于石油化工、能源动力、低温制冷及航空航天等学科领域中。在石油化工中,常常需要处理不同组份或者不同相态物质的流动问题,这种体系通常被称作多相体系,相应的流动被称为多相流[1]。两相以上的流动广泛存在于石油化工中,如石油、天然气、水三相流动,甚至油、气、水、沙四相同时流动。 随着工业科技的发展，多相管流的原理研究及其应用俨然是石油化工领域的重要研究课题。但时至目前，研究人员仍不能彻底掌握多相流的原理，对许多问题在实验和理论方面尚存在有不一致的结论和认识,所以也还不能提出精确描述多相流动特征的模型，多相流总体的理论体系框架还很不完善，高精度、适用范围宽、便于实际应用的多相流技术还有待逐步发展,数值计算方法也还有发展空间[2]。 2.多相管流研究历史和现状 2.1研究历史 1949年，Lockhart-Martineli[3]首次提出利用经验关系代数式来求解多相管流的流动特性参数,直到上个世纪80年代,多相管流的学者都在探讨较广泛的流动前提下,广泛得到实验数据和现场数据，并利用统计学等数学的方法,得出相应经验或者半经验的公式。其优势是计算简单,在实验范围内计算精度高，其缺陷是只适合在特定的条件下应用。 此后研究学者们几乎都将科研重心转到了机理研究上,探究多相管流的物理特性,建立相应数学物理模型,用数值方法简化并利用计算器得出其解。Barnea[4]提出的统一流动模型适用任意倾角下的流型判断,并指出各种流型之间变化的规律。Taitel[5]的段塞流动适用于任意倾角，根据质量守恒定理得出单元段塞内的平均持液率，根据动量定理得出液膜区长度和一个段塞单元长度,根据作用力定理得出通过段塞单元的压降。Alves[6]的统一环状流模型在任意倾角均适用,发现环状流和分层流流型有许多共同点,不同点是环状流的气芯在中，使气液相界面间摩阻系数计算方式不同。Gozllez[7]的统一泡状流数理模型能对泡状流中的液体持液率进行精确预测。统一流动模型的优势是其精度较高，可以得到流型的流动原理并且有更宽泛的适用性。 随着多相管流的研究深入，研究者逐渐从研究现象向原理研究发展，因为研究多相流机理才可以真正抓住流型及各种流型变换的本质。B.D.wood[8]等人对段塞流在下倾管中运动原理进行分析,发现下倾管为大倾角时，不会出现段塞流，而仅仅在下倾角较小时，由于长波出现导致流型不稳定才会出现段塞流。Taitel[9]等人对水平管、上倾管中产生段塞流的本质进行深入的分析研究,认为多相管路的扰动、管路倾斜角度的变化以及气液两相流量的变化等情况使管中出现Kelvin-Helmholtz小波的不稳定现象是水平管和上倾管中易出现段塞流的本质因
多相流模拟
多相流模拟介绍 自然界和工程问题中会遇到大量的多相流动。物质一般具有气态、液态和固态三相，但是多相流系统中相的概念具有更为广泛的意义。在多项流动中，所谓的“相”可以定义为具有相同类别的物质，该类物质在所处的流动中具有特定的惯性响应并与流场相互作用。比如说，相同材料的固体物质颗粒如果具有不同尺寸，就可以把它们看成不同的相，因为相同尺寸粒子的集合对流场有相似的动力学响应。本章大致介绍一下Fluent中的多相流建模。 多相流动模式 我们可以根据下面的原则对多相流分成四类： ?气-液或者液-液两相流： o气泡流动：连续流体中的气泡或者液泡。 o液滴流动：连续气体中的离散流体液滴。 o活塞流动：在连续流体中的大的气泡 o分层自由面流动：由明显的分界面隔开的非混合流体流动。 ?气-固两相流： o充满粒子的流动：连续气体流动中有离散的固体粒子。 o气动输运：流动模式依赖诸如固体载荷、雷诺数和粒子属性等因素。最典型的模式有沙子的流动，泥浆流，填充床，以及各向同性流。 o流化床：由一个盛有粒子的竖直圆筒构成，气体从一个分散器导入筒内。从床底不断充入的气体使得颗粒得以悬浮。改变气体的流量，就会有气泡不断的出现并穿过整个容器，从 而使得颗粒在床内得到充分混合。 ?液-固两相流 o泥浆流：流体中的颗粒输运。液-固两相流的基本特征不同于液体中固体颗粒的流动。在泥浆流中，Stokes数通常小于1。当Stokes数大于1时，流动成为流化（fluidization）了的液 -固流动。 o水力运输：在连续流体中密布着固体颗粒 o沉降运动：在有一定高度的成有液体的容器内，初始时刻均匀散布着颗粒物质。随后，流体将会分层，在容器底部因为颗粒的不断沉降并堆积形成了淤积层，在顶部出现了澄清层， 里面没有颗粒物质，在中间则是沉降层，那里的粒子仍然在沉降。在澄清层和沉降层中间， 是一个清晰可辨的交界面。 三相流(上面各种情况的组合) 多相系统的例子 ?气泡流例子：抽吸，通风，空气泵，气穴，蒸发，浮选，洗刷 ?液滴流例子：抽吸，喷雾，燃烧室，低温泵，干燥机，蒸发，气冷，刷洗 ?活塞流例子：管道或容器内有大尺度气泡的流动 ?分层自由面流动例子：分离器中的晃动，核反应装置中的沸腾和冷凝 ?粒子负载流动例子：旋风分离器，空气分类器，洗尘器，环境尘埃流动 ?风力输运例子：水泥、谷粒和金属粉末的输运 ?流化床例子：流化床反应器，循环流化床 ?泥浆流例子: 泥浆输运，矿物处理 ?水力输运例子：矿物处理，生物医学及物理化学中的流体系统 ?沉降例子：矿物处理 多相建模方法 计算流体力学的进展为深入了解多相流动提供了基础。目前有两种数值计算的方法处理多相流：欧拉－拉格朗日方法和欧拉－欧拉方法。 欧拉－拉格朗日方法
FLUENT中两相流多相流中模型的的选择问题
F L U E N T中两相流多相流中模型的的选择问题 GE GROUP system office room 【GEIHUA16H-GEIHUA GEIHUA8Q8-
两相流：通常把含有大量固体或液体颗粒的气体或液体流动称为两相流；其中含有多种尺寸组颗粒群为一个“相”，气体或液体为另一“相”，由此就有气—液，气—固，液—固等两相流之分。 两相流的研究：对两相流的研究有两种不同的观点：一是把流体作为连续介质，而把颗粒群作为离散体系；而另一是除了把流体作为连续介质外，还把颗粒群当作拟连续介质或拟流体。 引入两种坐标系：即拉格朗日坐标和欧拉坐标，以变形前的初始坐标为自变量称为拉格朗日Langrangian 坐标或物质坐标；以变形后瞬时坐标为自变量称为欧拉Eulerian 坐标或空间坐标。 一.离散相模型 FLUENT在求解连续相的输运方程的同时，在拉格朗日坐标下模拟流场中离散相的第二相；← 离散相模型解决的问题：煤粉燃烧、颗粒分离、喷雾干燥、液体燃料的燃烧等；← 应用范←围：FLUENT中的离散相模型假定第二相体积分数一般说来要小于10-12%（但颗粒质量承载率可以大于10-12%，即可模拟离散相质量流率等/大于连续相的流动）；不适用于模拟在连续相中无限期悬浮的颗粒流问题，包括：搅拌釜、流化床等； 颗粒-颗粒之间的相互作用、颗粒体积分数对连续相的影响未考虑；← 湍流中颗粒处理的两种模型：Stochastic Tracking，应用随机方法来考虑瞬时湍流速度对颗粒轨道的影响；Cloud Tracking，运用统计方法来跟踪颗粒围绕某一平均轨道的湍流扩散。通过计算颗粒的系统平均运动方程得到颗粒的某个“平均轨道”← 二.多相流模型 FLUENT中提供的模型： VOF模型(Volume of Fluid Model)←