Fluent多相流模型选择及设定.docx

1、Fluent多相流模型选择及设定1.多相流动模式我们可以根据下面的原如此对多相流分成四类：气-液或者液-液两相流：o 气泡流动：连续流体中的气泡或者液泡。o 液滴流动：连续气体中的离散流体液滴。o 活塞流动: 在连续流体中的大的气泡o 分层自由面流动：由明显的分界面隔开的非混合流体流动。气-固两相流：o 充满粒子的流动：连续气体流动中有离散的固体粒子。o 气动输运：流动模式依赖诸如固体载荷、雷诺数和粒子属性等因素。最典型的模式有沙子的流动，泥浆流，填充床，以与各向同性流。o 流化床：由一个盛有粒子的竖直圆筒构成，气体从一个分散器导入筒。从床底不断充入的气体使得颗粒得以悬浮。改变气体的流量，就会

2、有气泡不断的出现并穿过整个容器，从而使得颗粒在床得到充分混合。液-固两相流o 泥浆流：流体中的颗粒输运。液-固两相流的根本特征不同于液体中固体颗粒的流动。在泥浆流中，Stokes 数通常小于1。当Stokes数大于1 时，流动成为流化fluidization了的液-固流动。o 水力运输: 在连续流体中密布着固体颗粒o 沉降运动: 在有一定高度的成有液体的容器，初始时刻均匀散布着颗粒物质。随后，流体将会分层，在容器底部因为颗粒的不断沉降并堆积形成了淤积层，在顶部出现了澄清层，里面没有颗粒物质，在中间如此是沉降层，那里的粒子仍然在沉降。在澄清层和沉降层中间，是一个清晰可辨的交界面。三相流(上面各种

3、情况的组合)各流动模式对应的例子如下：气泡流例子：抽吸，通风，空气泵，气穴，蒸发，浮选，洗刷液滴流例子：抽吸，喷雾，燃烧室，低温泵，枯燥机，蒸发，气冷，刷洗活塞流例子：管道或容器有大尺度气泡的流动分层自由面流动例子：别离器中的晃动，核反响装置中的沸腾和冷凝粒子负载流动例子：旋风别离器，空气分类器，洗尘器，环境尘埃流动风力输运例子：水泥、谷粒和金属粉末的输运流化床例子：流化床反响器，循环流化床泥浆流例子: 泥浆输运，矿物处理水力输运例子：矿物处理，生物医学与物理化学中的流体系统沉降例子：矿物处理2. 多相流模型FLUENT中描述两相流的两种方法:欧拉一欧拉法和欧拉一拉格朗日法，后面分别简称欧拉法

4、和拉格朗日法。欧拉法即为两相流模型，拉格朗日法即为离散相模型欧拉法着眼于空间的点，根本思想是考察空间一个点上的物理量与其变化。在欧拉方法中，FLUENT将不同的相被处理成互相贯穿的连续介质。各相的体积率是时间和空间的连续函数，其体积分率之等于1。欧拉法中两相流模型包括:VOF(the volumeoffluid)模型，混合模型和欧拉一欧拉模型VOF模型(Volume of Fluid Model) 混合模型(Mixture Model) 欧拉模型(Eulerian Model) 2.1 VOF模型(Volume of Fluid Model) 􀂑 VOF模型用来处理没有相互穿

5、插的多相流问题，在处理两相流中，假设计算的每个控制容积中第一相的体积含量为1，如果1=0，表示该控制容积中不含第一相，如果1=1，如此表示该控制容积中只含有第一相，如果011，表示该控制容积中有两相交界面；􀂑 VOF方法是用体积率函数表示流体自由面的位置和流体所占的体积，其方法占存小，是一种简单而有效的方法。VOF模型在应用的过程中存在某些局限性:(l)在利用该模型进展模拟时要求所有的控制体积必须被任何一种流体相或混合相所填满，即不能存在无流体流动的区域;(2)只允许一相流体是可压缩的;(3)很难对具有混合物料和反响存在的流动进展模拟;(4)相间存在较大速度差时，界面的速度精

6、度会受到很大的影响。2.2 混合模型(Mixture Model) 混合模型(Mixture Model)是一种简化的两(多)相流模型，它使用单流体方用于模拟各相有不同速度的两(多)相流，但是假定了在短空间尺度上局部的，相之间的耦合很强。同时也用于模拟有强烈藕合的各向同性相流和各相以相度运动的两(多)相流。混合模型可以通过求解混合相的动量、连续性和能量，第二相的体积分率方程，以与相对速度的代数表达式模拟多相(fluldorculate)。典型的应用包括低负载的粒子负载流，沉降，旋风别离器以与气相容很低的泡状流。混合物模型也可用于没有离散相相对速度的均匀多相流。􀂑 用混合特性参

7、数描述的两相流场的场方程组称为混合模型； 􀂑 考虑了界面传递特性以与两相间的扩散作用和脉动作用；使用了滑移速度的概念，允许相以不同的速度运动； 􀂑 用于模拟各相有不同速度的多相流；也用于模拟有强烈耦合的各向同性多相流和各相以一样速度运动的多相流；􀂑 缺点：界面特性包括不全，扩散和脉动特性难于处理。 2.3 欧拉模型(Eulerian Model) 欧拉一欧拉模型(Euler-Euler Model)是两(多)相流中最复杂的两(多)相流模型，也称为双流体模型。连续相与分散相被视为连续的一体。欧拉一欧拉模型对每一相都建立动量方程和连续性方程，通

8、过压力和相间交换系数的藕合来计算求解。欧拉模型的应用包括气泡柱、颗粒悬浮以与流化床的模拟。有人将其成功地应用欧拉-模型模拟了鼓泡塔中两(多)相流的模拟与气泡聚并和破碎的影响。􀂑 欧拉模型指的是欧拉欧拉模型；􀂑 把颗粒和气体看成两种流体，空间各点都有这两种流体各自不同的速度、温度和密度，这些流体其存在同一空间并相互渗透，但各有不同的体积分数，相互间有滑移； 􀂑 颗粒群与气体有相互作用，并且颗粒与颗粒之间相互作用，颗粒群紊流输运取决于与气相间的相互作用而不是颗粒间的相互作用； 􀂑 各颗粒相在空间中有连续的速度、温度与体积分数分

9、布。 几种多相流模型的选择􀂑 VOF模型适合于分层流动或自由外表流；􀂑 Mixture和Eulerian模型适合于流动中有混合或别离，或者离散相的体积份额超过10%12的情况。Mixture模型和Eulerian模型区别􀂑 如果离散相在计算域分布较广，采用 Mixture模型；如果离散相只集中在一局部，使用Eulerian模型；􀂑 当考虑计算域的interphase drag laws 时，Eulerian模型通常比Mixture模型能给出更准确的结果；􀂑 从计算时间和计算精度上考虑。拉格朗日法着眼于流体

10、的质点，根本思想是跟踪每个流体质点在流动过程中的运动全过程，记录每个质点在每一时刻、每一位置的各个物理量与变化。在拉格朗日方法中，FLUENT将主体相视为连续相，稀疏相视为离散颗粒，主体相用欧拉法，而离散相利用拉格朗日法进展粒子跟踪，这就是所谓的欧拉一拉格朗日模型。此模型中需要离散相体积含量不超过15%，离散相和主体相都有自己的压力、粘度与湍流扩散稀疏参数，并在拉格朗日坐标系中考察离散相颗粒的运动轨迹。该模型能详细地分析粒子/液滴间的作用力以与流体间复杂的作用力，防止了应用大量的经验关系，又防止了离散相数值解的扩散问题，虽然计算量庞大，但是相对欧拉模型来讲，精度要更高一些。比拟了各种模型，认为

11、离散相模型能更准确地模拟气固两相流动，能更好的跟踪固体颗粒、气泡、液滴在连续相中运动轨迹。通常，你一旦决定了采用何种模式最能符合实际的流动，那么就可以根据以下的原如此来挑选最优的模型。对于体积率小于10的气泡、液滴和粒子负载流动，采用离散相模型。对于离散相混合物或者单独的离散相体积率超出10的气泡、液滴和粒子负载流动，采用混合物模型或者欧拉模型。对于活塞流，采用VOF 模型。对于分层/自由面流动，采用VOF 模型。对于气动输运，如果是均匀流动，如此采用混合物模型；如果是粒子流，如此采用欧拉模型。对于流化床，采用欧拉模型模拟粒子流。对于泥浆流和水力输运，采用混合物模型或欧拉模型。对于沉降，采用欧

12、拉模型。对于更加一般的，同时包含假如干种多相流模式的情况，应根据最感兴趣的流动特征，选择适宜的流动模型。此时由于模型只是对局部流动特征做了较好模拟，其精度必然低于只包含单个模式的流动。Fluent软件中对喷雾这类气液两相流问题的模拟主要采用其自带的离散相模型(DPMDiscrete Phase Model)。此模型是以欧拉拉格朗日方法为根底建立的。它把流体作为连续介质，在欧拉坐标系加以描述，对此连续相求解输送方程，而把雾滴颗粒群作为离散体系，通过积分拉氏坐标系下的颗粒作用力微分方程来求解离散相颗粒的轨道，可以计算出这些颗粒的轨道以与由颗粒引起的热量/质量传递。同时，在计算中，相间耦合以与耦合结

13、果对离散相轨道、连续相流动的影响均可考虑进去。当计算颗粒的轨道时，Fluent跟踪计算颗粒沿轨道的热量、质量、动量的得到与损失，这些物理量可作用于随后的连续相的计算中去。于是，在连续相影响离散相的同时，用户也可以考虑离散相对连续相的作用。交替求解离散相与连续相的控制方程，直到二者均收敛二者计算解不再变化为止，这样，就实现了双向耦合计算。在采用FLUENT中的离散相模型时，需要定义每个粒子尺寸以与温度。这些初始条件以与有关离散相物理性质的输入量/质量计算的必要条件。轨迹以与热量/质量传递的计算是粒子的对流或辐射传热、质量传递以与粒子在流场运动时的。而预测所得的轨迹以与相关的质量、热量传递可以通过

14、1稳态问题建立与求解程序纲要建立和求解稳态离散相问题的一般程序如下所示:(l)求解连续相流动;(2)生成离散相的入射;(3)根据需要选择是否连续相与离散相关联求解;(4)用画图或者提取数据来跟踪离散相入射。2非稳态问题建立与求解程序纲要建立和求解非稳态离散相问题的一般程序如下所示:(l)生成离散相入射;(2)初始化流场;(3)取适宜的时间步长数目进展求解。随着求解的进展，粒子的位置将会被更新。利用Fluent自带的空气雾化喷嘴模型预测雾化颗粒的颗粒行为。首先假设不带颗粒的空气为连续相，对其进展单相模拟。之后，假设雾化喷嘴喷出的甲烷颗粒为离散相，进展了气液两相耦合模拟。单相稳态模拟的根底上打开D

15、PM模型Discrete Phase Model 参加离散相甲烷雾滴进展两相耦合模拟，重点介绍了DPM中参数的设定。1 打开DPM模型利用Define/Models/Discrete Phase Model打开DPM，本文截取了Discrete Phase Model设置面板的一局部，对其中参数的设定进展详细的分析，如图1所示。图1 Discrete Phase Model面板当模拟两相耦合过程时，用户应该首先计算得到收敛或局部收敛的连续相流场，然后再创建喷射源进展耦合计算。在每一轮离散相的计算，FLUENT 计算颗粒/液滴轨迹并且更新每一个流体计算单元的相间动量、热量以与质量交换项。然后，这

16、些交换项就会作用到随后的连续相的计算。耦合计算时FLUENT 在连续相迭代计算的过程中，按照一定的迭代步数间隔来计算离散相迭代。直到连续相的流场计算结果不再随着迭代步数加大而发生变化即，达到了所有的收敛标准，耦合计算才会停止。当达到收敛时，离散相的轨迹也不再发生变化假如离散相轨迹发生变化将会导致连续相流场的变化。耦合计算的设定步骤如下：1. 计算连续相流场；2. 在Discrete Phase Model panel 面板中，激活Interaction with Continuous Phase 选项；3. 在Number Of Continuous Phase Iterations Per

17、DPM Iteration 文本框中设定颗粒轨迹的计算频率即连续相迭代多少步，就进展一轮离散相的计算。假如用户设定此参数为5，即意味着在连续相进展了五步迭代之后，就开始离散相的迭代计算。两个离散相计算中间应该间隔多少连续相的迭代步，要视用户问题的物理意义而定。需要注意的是，【*假如此参数设定为0，那么FLUENT 将不进展离散相的计算。】另外，图1中绿色圈的2个参数是最大计算步数Max. Number Of Steps和积分尺度Length Scale。最大计算步数Max. Number Of Steps是用积分方程(1)，(2) 求解颗粒轨道时，允许的最大时间步数。当某个颗粒轨道计算达到此时

18、间步数时，FLUENT 就自动中止了此颗粒的轨道计算，输出时，此颗粒被标记为“inplete。对最大时间步数的规定消除了对某些在流场中不停循环的颗粒的无休止的计算。但是，对于缺省的500 步的最大时间步数，很多问题的计算都不止这么多。这种情况下，当颗粒信息在输出时被标记未完成，而实际颗粒并不是在流场中无休止的打转，那么，用户可以增加最大时间步数注值得注意的是：设定上述各个参数的一个简便方法是，假如用户希望颗粒穿越长度为D的计算域，那么用长度标尺乘以最大积分时间步数，其结果应该大致等于D，即等于所设定的Number Of Continuous Phase Iterations Per DPM I

19、teration的值。2 创建injection 通过Define/injection/create进入创建injection面板，如如下图所示：在Injection Type中选择射流源类型，本文选定空气雾化喷嘴(air-blast-atomizer)。在Particle Type中选择颗粒类型，本文选择Droplet液滴是一种存在于连续相气流中的液体颗粒。它服从力的平衡并受到加热冷却的影响由定律1 确定。此外，他还由定律2 和3 确定自身的蒸发与沸腾请参阅Users Guide中的。只有传热选项被激活并且至少两种化学组份在计算中是被激活的，或者已经选择了非预混燃烧或局部预混燃烧模型，液滴类

20、型才是可选的。当选择了液滴类型之后，用户应该使用理想气体定律来定义气相密度在空气辅助雾化模型里，用户应直接设定液膜厚度，如图3所示。在Point Properties面板上，设定喷口处液膜的外半径，即液膜的厚度。另外，用户还必须设定液膜与空气间的最大相对速度差和喷射角度，如图4所示。液膜离开喷口之后，它的初始轨道沿着设定的喷射角。注意：如果初始液膜的轨道指向中心线，那么，喷射角度为负值。3 离散相边界条件的设定在Discrete Phase Model Conditions 属性框下的Boundary Cond. Type 下拉框中选择reflect，trap，或escape 边界条件在面板中

21、，需要点击DPM 才能激活Discrete Phase Model conditions。如图5所示。FLUENT 中的离散相缺省边界条件为：1.壁面wall、对称面symmetry、轴对称的轴线axis均为reflect边界条件，且恢复系数均为；2.所有的流动类型边界压力入口-pressure inlets、速度入口-velocity inlets、压力出口-pressure outlets 等，均为escape边界条件；3.所有的部区域边界辐射体- radiator、多孔介质连续面- porous jump均为边界条件；4.有对壁面边界wall才可以修改恢复系数。注意：在Boundary

22、Conditions 面板打开的面板中可以设定离散相边界条件。当设定完一个以上的喷射源之后，离散相边界条件的输入项就会出现在相应的面板中。4 模拟结果与后处理颗粒轨道的输出时，颗粒的可能的结果如下：1Escaped：逃逸意味着颗粒在已经设定了逃逸边界条件的流动边界终止了轨迹的计算。2Inplete：未完成：意味着颗粒轨迹的计算时间步长已经达到设定的最大步数在Discrete Phase Model panel 面板中的Max. Number Of Steps 文本框中设定，3Trapped：捕获：意味着颗粒在已经设定了捕集边界条件的流动边界终止了轨迹的计算。4Evaporated：蒸发：意味着

23、颗粒在计算域中被完全蒸发掉了。5Aborted：忽略：意味着颗粒由于舍入误差原因而不能进展计算。用户可以修改长度标尺或设定不同的初始条件来重新计算颗粒轨迹。需要注意的是，除了用连续相的变量值来着色颗粒轨迹外，也可以使用离散相的各种变量值来进展着色。这些变量值包括：颗粒已停留时间、颗粒速度、颗粒直径、颗粒密度、颗粒质量、颗粒温度、颗粒所使用的定律、颗粒积分时间步长、颗粒雷诺数。在Color By类目框下的Particle Variables.下拉框中列出了所有可选的着色颗粒变量。为了显示计算域的最大/最小值，可以点击Update Min/Max 按钮更新。DPM 模型的根本操作和须知事项1 DP

24、M 模型概述DPM 模型可以用来模拟流场中的离散相，它的特点是使用方便，模拟思路清晰，计算中可以对颗粒运动轨迹进展跟踪，结果直观；其缺点是，计算结果无法得到离散相各种场图，为结果分析造成很大不便。FLUENT 提供了如下的离散相模型选项：1. 使用Lagrangian 坐标下的公式计算颗粒的轨迹。这些公式涉与了稳态与非稳态条件下离散相的惯性力、曳力和重力。2. 连续相中的漩涡对于离散相扩散产生的扰动进展预测。3. 离散相的加热与冷却。4. 液滴的蒸发和沸腾。5. 提供对颗粒燃烧的模拟，可以通过对挥发份析出和焦炭燃烧来模拟煤粉的燃烧。6. 可以选择是否进展连续相与离散相的耦合计算。7. 液滴的破

25、碎与合并。这些模型时的FLUENT 可以用来对许多种离散相的问题进展模拟，包括颗粒的别离与分级，喷雾枯燥，烟雾的扩散，液体中气泡的搅浑，液体燃料和煤的燃烧。当需要在FLUENT 的模型中参加离散相时，可以通过定义颗粒的初始位置、速度、粒径、温度等参数实现，具体的操作过程在“Discrete PhaseModel面板中完成。以上的参数再加上颗粒的物理属性，就可以作为计算颗粒轨迹和颗粒热、质传递的初始化条件。下面就使用DPM 模型的根本步骤归纳如下：对于稳态问题，可采用以下步骤求解：1. 求解连续相流动；2. 添加离散相；3. 如果需要的话可以求解耦合流动；4. 对计算结果进展后处理对于非稳态问题

26、，可通过以下步骤求解；1. 添加离散相；2. 初始化流场；3. 设定时间步长。对于非耦合问题，FLUENT 会在每个时间步长的最后更新离散相的位置；对于耦合问题，在每次相间耦合计算中离散相的位置都回更新。2 应用DPM 模型需要注意的一些问题在Fluent 中应用DPM 模型进展计算时，需要注意DPM 模型忽略了两相流中颗粒之间的相互作用，以与颗粒相对连续相流动产生的影响。这就决定了两相流中颗粒相的体积分数不能太高，通常情况下这一体积分数要小于10%20%。但是，这并不意味着在应用DPM 模型时颗粒相的质量分数也要小于10%20%，实际上，我们可以使用DPM 模型来模拟离散相质量分数等于或超过

27、连续相质量分数的流动。*【1如果颗粒是以喷射的形式进入连续相的，而且流场中有明确的入口和出口，这种情况下可以使用稳态的DPM 模型来计算；2如果颗粒相在连续相中处于一种无限期的悬浮状态，这种情况下稳态的Lagrangian 模型就不再适用了，对于这样的工况可以考虑使用非稳态的DPM 模型来进展求解。换句话说，对于搅拌器、混和器、流化床这一类容器如果应用DPM 模型来模拟其流场，应该在非稳态的前提下进展。】一旦应用DPM 模型来对流动进展模拟后，Fluent 中的某些功能将不能再被使用。具体如下：1. 周期性的边界条件；2. 可调的时间步长；3. 使用非预混燃烧模型时，颗粒不能参加反响；4. 当

28、使用动网格或变形网格时，颗粒喷射的外表便不能随网格一起运动；5. 如果使用了复合参考系，在参考系下颗粒轨道失去了原有的意义，同理，相间耦合计算也失去了意义。解决这个问题的方法就是采用绝对速度来对颗粒进展跟踪而不是采用相对速度，这一方法可以通过在文本窗口输入以下命令实现：define/models/dpm/tracking/track-inabsolute-frame。需要注意的是，计算结果会与符合参考系下壁面的位置有很大关系。颗粒的跟踪是在哪个参考坐标系下进展的，颗粒的入射速度就要在哪个参考坐标系下定义的。默认情况下，颗粒速度是基于当地坐标系定义的，如果你激活了track-in-absolut

29、e-frame方法如前所述，颗粒速度就基于绝对坐标系来定义。3DPM 模型的傻瓜用法所谓的傻瓜用法，就是不用考虑细节，甚至不必知道模型设置面板中每一项的意义所在，而只给出相应参数的设定来进展求解。我们不提倡这样的做法，但这也确实是能让新手尽快上路的好方法，当然，有可能计算的结果不准确，但对于简单的流场来讲，应该还可以承受，对于稍复杂的情况，即便是老手，也不敢保证一次建模、一次计算就能得到满意的结果，所以，慢慢调试吧!对于稳态的工况，为了确保计算结果的收敛，可以暂时先不在流场中添加离散相，而仅仅进展连续相的迭代，一直迭代到连续相收敛再参加离散相。当然，也可在计算得到收敛趋势时参加离散相。本节只讨

30、论DPM 模型面板的设定。下面说明傻瓜用法的操作步骤：1. 通过DefineModelsDiscrete Phase 来打开DPM 模型的控制面板，2. 选中interaction with Continuous Phase；3. 将Number of Continuous Phase Iterations per DPM Iteration 置为20；4. 选中Specify Length Scale，将Length Scale 置为，注意LengthScale 后面的单位是m；5. 粗略估计颗粒的行程，然后用该行程除以Length Scale，得到的值就是Max. Number Of St

31、eps 要输入的值。实际上，Length Scale 与Max.Number Of Steps 的乘积即为跟踪颗粒轨迹的最大长度，如果你想观察颗粒在整个流场中的流动，那么这个乘积的值就要大于颗粒的轨迹长度，所以此时可以适当地扩大Max. Number Of Steps 的值。6. 点击面板下方的injections，弹出Injections 面板，再点击Create，弹出Set Injection Properties 面板，在此面板中设定颗粒的属性。7. 在Point Properties 下输入颗粒的各种参数；8. 在Turbulent Dispersion 下激活Stochastic Tracking 选项，将Number ofTries 改成10。至此，DPM 模型的根本设定就全部完毕了。接下来的任务就是针对自己模型的特点，有针对性的到帮助文件中去寻找解决问题的方法。_我们先看看燃烧中的组分输运和有反响流动该如何处理。这是燃烧问题中很重要的一局部，前人开展了很多模型来处理不同的具体问题：a) 通用有限速度模型该方法基于组分质量分数的输运方程，采用你所定义的化学反响机制，对化学反响进展模拟。反响速度在这种方法中是以源项的形式出现在组分输运方程中的，计算反