多孔介质模型.docx

《多孔介质模型.docx》由会员分享，可在线阅读，更多相关《多孔介质模型.docx（38页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、分量来定义。 6. 定义辐射参数 FLUENT6.1全攻略 图8-26 Solid（固体）面板 如果使用DO模型计算辐射过程，可以在Participates in Radiation（是否参与辐射）选 项中确定固体区域是否参与辐射过程。 8.19 多孔介质条件 很多问题中包含多孔介质的计算，比如流场中包括过滤纸、分流器、多孔板和管道集 阵等边界时就需要使用多孔介质条件。在计算中可以定义某个区域或边界为多孔介质，并 通过参数输入定义通过多孔介质后流体的压力降。在热平衡假设下，也可以确定多孔介质 的热交换过程。 在薄的多孔介质面上可以用一维假设“多孔跳跃（porous jump）”定义速度和压强的

2、降 落特征。多孔跳跃模型用于面区域，而不是单元区域，在计算中应该尽量使用这个模型， 因为这个模型可以增强计算的稳定性和收敛性。 9 FLUENT6.1全攻略 8.19.1 多孔介质模型的假设和限制条件 多孔介质模型采用经验公式定义多孔介质上的流动阻力。从本质上说，多孔介质模型 就是在动量方程中增加了一个代表动量消耗的源项。因此，多孔介质模型需要满足下面的 限制条件： （1）因为多孔介质的体积在模型中没有体现，在缺省情况下，FLUENT在多孔介质内 部使用基于体积流量的名义速度来保证速度矢量在通过多孔介质时的连续性。如果希望更 精确地进行计算，也可以让FLUENT在多孔介质内部使用真实速度，详情

3、见8.19.7节。 （2）多孔介质对湍流的影响仅仅是近似。 （3）在移动坐标系中使用多孔介质模型时，应该使用相对坐标系，而不是绝对坐标系， 以保证获得正确的源项解。 8.19.2 多孔介质的动量方程 在动量方程中增加一个动量源项可以模拟多孔介质的作用。源项由两部分组成：一个 粘性损失项，即方程（8-45）右端第一项；和一个惯性损失项，即方程（8-45）右端第二 项。 = + = 3 3 1 SiDij vjC v v （8-45） ijmagj 12 j j 1 1 式中Si是第i个（x、y或z方向）动量方程中的源项，D和C是给定矩阵。负的源项 又被称为“汇”，动量汇对多孔介质单元动量梯度的贡

4、献，在单元上产生一个正比于流体速 度（或速度平方）的压力降。 在简单、均匀的多孔介质上，还可以使用下面的数学模型： = + 1 2 （8-46） Si vi C vmagvi 2 式中 代表多孔介质的渗透性，C2是惯性阻力因子，将D和C分别定义为由1/和 C2为对角单元的对角矩阵。 = 1 = 1 CC Si C vC vv 1 （8-47） 式中C0和 C1为用户自定义的经验常数。其中压力降是各向同性的，C0的单位为国际 FLUENT中还可以用速度的指数律作为源项的模型，即： () 00 i 10 单位制。 1. 多孔介质的Darcy定律 FLUENT6.1全攻略 在流过多孔介质的层流中，压

5、力降正比于速度，常数C2可以设为零。忽略对流加速和 扩散项，多孔介质就简化为Darcy定律： pvG = （8-48） 3 = px v n jx j 1 xj FLUENT在x、y、z三个坐标方向计算出的压力降为： = 3 = = jy j 1 yj py v n 3 = = jz j 1 zj pz v n （8-49） 式中的1/ij 就是方程（8-45）中的D，vj是 x、y、z三个坐标方向的速度分量，nx、 ny、 nz是多孔介质在 x、y、z 三个坐标方向的真实厚度。如果计算中使用的厚度值不 等于真实厚度值，则需要对1/ij 做出调整。 2. 多孔中的惯性损失 在流速很高时，方程（

6、8-45）中的常数C2可以对惯性损失做出修正。C2可以被看作 流动方向上单位长度的损失系数，这样一来就可以将压力降定义为动压头的函数。 如果计算的是多孔板或管道阵列，在一些情况下可以略去渗透项，而只保留惯性损失 = 3 1 1 pC ij vjvmag （8-50） 或者写成分量形式： 2 2 j 项，并产生下面简化形式的多孔介质方程： j = 11 FLUENT6.1全攻略 3 1 = pxC xj nx vjvmag = 22 22 j 1 3 1 pyC yj ny vjvmag （8-51） = 22 j 1 1 3 1 22 pzC zj nz vjvmag 22 j 1 再次说明，

7、上式中的厚度即模型中定义的厚度。 8.19.3 多孔介质能量方程的处理 多孔介质对能量方程的影响体现在对对流项和时间导数项的修正上。在多孔介质对对 流项的计算中采用了有效对流函数，在时间导数项中则计入了固体区域对多孔介质的热惯 性效应： GG t+ 1 + + = + f E fs Esvf E f pkeff ThiJivS f ()( ()()h i （8-52） 式中Ef 为流体总能，Es为固体介质总能， 为介质的多孔性，keff 介质的有效热导 率，Shf 流体焓的源项。 FLUENT中使用的有效导热率keff 是流体导热率和固体导热率的体积平均值： keff = kf + (1)ks

8、 （8-53） 式中 为介质的多孔率，kf 为流体的热导率，ks为固体的热导率。kf 和 ks可以用自 定义的函数计算。各向异性的热导率也可以用自定义函数定义。 8.19.4 多孔介质模型对湍流的处理 在缺省情况下，FLUENT 在多孔介质计算中通过求解标准守恒型方程计算湍流变量。 在计算过程中，通常假设固体介质对湍流的生成和耗散没有影响。在多孔介质的渗透率很 大，因而介质的几何尺度对湍流涡结构没有影响时，这个假设是合理的。在另外一些情况 12 FLUENT6.1全攻略 中，可能不需要考虑流动中的湍流，即假定流动为层流。 如果计算中使用的湍流模型是k 、k 或 Spalart-Allmaras

9、 模型中的一种，则可 以通过将湍流粘度t设为零的方式忽略湍流的影响。如果湍流粘度t被设置为零，则在 计算过程中仍然会将湍流变量输运到介质的另一面，但是湍流对动量输运过程的影响则完 全被消除。在Fluid（流体）面板中将多孔介质区设为Laminar Zone（层流区）选项就可以 将湍流粘度t设为零。 8.19.5多孔介质对瞬态标量方程的影响 对于多孔介质的瞬态计算，多孔介质对时间导数项的影响体现在对所有变量的输运方 程和连续性方程求解中。在考虑多孔介质的影响时，时间导数项变为() t ，式中为 任意流动变量， 为多孔率。在瞬态流动计算中，多孔介质的影响是被自动加入计算过程 的，其多孔率在缺省设置

10、中等于1。 8.19.6 多孔介质计算中用户的输入参数 多孔介质计算中需要输入的项目如下： （1）定义多孔介质区域。 （2）定义多孔介质速度函数形式。 （3）定义流过多孔介质区的流体属性。 （4）设定多孔区的化学反应。 （5）设定粘性阻力系数。 （6）设定多孔介质的多孔率。 （7）在计算热交换的过程中选择多孔介质的材料。 （8）设定多孔介质固体部分的体热生成率。 （9）设定流动区域上的任意固定值的流动参数。 （10）需要的话，将多孔区流动设为层流，或取消湍流计算。 （11）定义旋转轴或区域的运动。 所有参数设置均在流体（Fluid）面板中完成，见图8-27。 1. 定义多孔区 选中Porous

11、 Zone（多孔介质区）选项就可以将流体区域设为多孔介质。 13 2. 定义多孔介质速度公式 FLUENT6.1全攻略 图8-27 多孔计算选项 在Solver（求解器）面板中有一个Porous Formulation（多孔公式）区可以确定在多孔 介质区域上使用名义速度或物理速度。缺省设置为名义速度。 3. 定义流过多孔介质的流体 在Material Name（材料名称）中选择所需的流体名称即可。可以用编辑功能改变流体 的参数设置。组元计算或多相流计算中的流体不在这里定义，而是在Species Model（组元 模型）面板中定义。 4. 在多孔区域上设置化学反应 在Fluid（流体）面板上选中

12、Reaction（反应）选项，再从Reaction Mechanism（反应 14 FLUENT6.1全攻略 机制）中选择合适的反应就可以在多孔介质区域的计算中加入化学反应。 如果化学反应中包括表面反应，则需要设定Surface to Volume Ratio（面体比）。面体比 是多孔介质单位体积上拥有的表面积，因此可以作为催化反应强度的度量。根据这个参数， FLUENT可以计算出体积单元上总的表面积。 5. 定义粘性和惯性阻力系数 粘性和惯性阻力系数的定义方式是相同的。在直角坐标系中定义阻力系数的办法是： 在二维问题中定义一个方向矢量，或在三维问题中定义两个方向矢量，然后再在每个方向 上定义

13、粘性和惯性阻力系数。在二维计算中的第二个方向，即没有被显式定义的那个方向， 是与被定义的方向矢量相垂直的方向。与此类似，在三维问题中的第三个方向为垂直于前 两个方向矢量构成平面的方向。在三维问题中，被定义的两个方向矢量应该是相互垂直的， 如果不垂直的话，FLUENT 会将第二个方向矢量中与第一个方向矢量平行的分量删除，强 制令二者保持垂直。因此第一个方向矢量必须准确定义。 用UDF也可以定义粘性和惯性阻力系数。在UDF被创建并调入FLUENT后，相关的 用户定义选项就会出现在下拉列表中。需要注意的是，用 UDF 定义的系数必须使用 DEFINE_PROPERTY宏。 如果计算的问题是轴对称旋转

14、流，可以为粘性和惯性阻力系数定义一个附加的方向分 量。这个方向应该与其他两个方向矢量相垂直。 在三维问题中，还允许使用圆锥（或圆柱）坐标系。需要提醒的是，多孔介质流中计 算粘性或惯性阻力系数时采用的是名义速度。 6. 定义阻力的步骤 定义阻力系数的步骤如下： （1）定义方向矢量。 1）如果采用直角坐标系，在二维计算中只要定义第一方向矢量（Direction-1 Vector） 即可，在三维计算中则还要定义第二方向矢量（Direction-2 Vector）。矢量的起点都 是（0,0）或（0,0,0），终点则用坐标（x,y）或（x,y,z）定义。 在无法确定多孔介质的方向矢量时，可以使用三维的平

15、面工具或二维的线工具定义 方向矢量： 在多孔介质边界上启动plane tool（面工具）或line tool（线工具）。 旋转面工具或线工具的轴，使之与多孔介质相吻合。 （从面工具更新）或Update From Line Tool（从线工具更新），FLUENT就会将工具的红箭头方向设为第 一方向矢量。在三维问题中，同时将绿色箭头方向设为第二方向矢量。 2）在三维问题中还可以使用圆锥坐标系完成设置，具体步骤如下： 选择Conical（圆锥）选项。 定义Cone Axis Vector（圆锥轴矢量）和Point on Cone Axis（圆锥轴上的点）。 从 Fluid（流体）面板中点击Updat

16、e From Plane Tool 15 FLUENT6.1全攻略 圆锥轴矢量的方向由起点为（0,0,0）终点为（x,y,z）的矢量确定。而圆锥轴 上的点则用于将阻力转换为直角坐标系中的值。 设置Cone Half Angle（半锥角）。如果使用的是圆柱坐标系，则将半锥角设为 零。在事先不知道圆锥轴的方向和锥面上一点坐标的时候，可以用面工具设 置圆锥轴矢量和锥面上一点的坐标。一种方法是： （a）在Snap to Zone（转换到区域）按钮旁边选择垂直于圆锥轴的边界区域。 （b）点击转换到区域按钮，FLUENT会自动将面工具设在边界上。同时设定 的有Cone Axis Vector（圆锥轴矢量）

17、、Point on Cone Axis（圆锥轴上一点）和 半锥角。 另一种方法是： （a）将面工具设在多孔区域上。 （b）转动、平移面工具的轴，使得红色箭头指向圆锥轴的方向，同时使面工 具的原点与圆锥轴的原点重合。 （c）在二者的轴和原点全部重合后，点击Update From Plane Tool（从面工具 更新）按钮，FLUENT就会自动设定圆锥轴矢量和圆锥上一点的坐标。 （2）在Viscous Resistance（粘性阻力）下，在每个方向上定义粘性阻力系数1/，在 Inertial Resistance（惯性阻力）下，在每个方向上定义惯性阻力系数C2。如果使用圆锥定 义方法，则方向一（D

18、irection-1）为圆锥轴方向，方向二（Direction-2）为垂直于锥面的方 向，对应于圆柱坐标系中的径向，方向三（Direction-3）为周向。 在三维问题中有三类系数，二维问题中则有两类： 1）在各向同性问题中，所有方向的阻力系数是相同的，比如海绵。在这类问题中， 必须明确地将各方向的阻力系数设为同一个值。 2）在某个方向的阻力系数与其他方向有所不同时，应该谨慎设置阻力参数。比如， 如果一个多孔介质区由圆柱形吸管组成，并且吸管的开口方向与流动方向平行，则 流体会比较容易地流过吸管，但是在其他两个方向则不易流过。如果有一个与流动 方向垂直的平板，则流体完全不能流过平板，而是沿平板向

19、其他两个方向流动。 3）三维问题中另一种可能情况是所有三个方向都是不同的。比如，如果多孔介质区 由不规则排列的物体组成，则流体在其中流过时在各个方向上的过程都是不同的。 因此需要在各个方向上设置不同的系数。 7. 粘性和惯性阻力系数的推导 推导粘性和惯性阻力系数的方法如下： （1）采用已知的压强损失，基于名义速度推导多孔介质参数。 使用多孔介质模型必须注意的是，FLUENT 中假设多孔介质单元是完全浸没在流体中 的，而用户定义的1/ij 和C2ij 都是基于这个假设得出的。假定真实容器中的压强降与速度 16 FLUENT6.1全攻略 的关系是已知的，而其中的多孔介质却是部分浸没在流体中的。通过

20、下面的练习可以掌握 一种计算C2的方法。 例题：一块多孔板有 25浸没在流体中，压强降为动压头的 0.5 倍，损失因子KL的 定义为： =2 225% 1 p KL v （8-54） 因此KL应该等于0.5，计算动压使用的速度是经过25浸没在流体中的多孔介质上的 速度。计算C2时应该注意： 1）在假定多孔板100浸没在流体中的前提下计算其速度。 2）损失系数必须被转化为单位长度上动压损失。 因此损失因子KL适用于下面的公式： =2 2100% 1 p KL v （8-55） 式中v25% = 4v100% ， 2 2 25% = = = v 4 0.5 8 v 8 2 100% 1 K L K

21、L （8-56） 根据第 2）条要求，还需要将KL转换为单位厚度多孔板上的损失系数。假定板的厚 度为1mm，则惯性损失因子为： C = KL 8 m 1 210 3 8000 厚度 （8-57） 对各向异性介质，上述过程需要在各坐标方向上重复进行。 （2）用Ergun方程计算充填床的多孔介质参数。 作为第二个例子，可以研究一下充填床问题。在湍流中，充填床的数学模型是用穿透 率和惯性损失系数来定义的。计算相关常数的一种办法是使用半经验公式 Ergun 方程，这 个方程适用的雷诺数范围很广，同时也使用于多种填充物： 17 FLUENT6.1全攻略 ()() 2 33 p 2 150 11.75 1

22、 L Dv D pp 2 =+v 2 （8-58） 在计算层流时，（8-58）式中右端第二项可以被去掉，Ergun 方程随之简化为 Blake-Kozeny方程： =v L D p 2 150 1 （8-59） () 2 3 p 在上述方程中 为粘度，Dp为粒子的平均直径，L为床的深度， 为空腔比率，定 义为空腔与填充床的体积比。对比方程（8-48）、（8-50）和（8-58）可知，各方向上的渗透 率和惯性损失系数为： 23 150 1 = Dp （8-60） () 23 ()2 3.5 1 = C （8-61） Dp Dp （3）用经验公式推导穿过多孔板的流动参数 作为第三个例子，下面用Va

23、n Winkle方程计算带方孔的多孔板上压强的损失。方程的 提出者认为该方程适用于呈三角形分布的等距方孔板的湍流计算，具体形式如下： m =CAf (2p)/(1(Af / Ap)2) （8-62） 式中m为通过板的质量流量，Af 为孔的总面积，Ap板的总面积（固体与孔的和）， D/t孔直径与板厚之比，C是随雷诺数和D/t变化的系数，其值可以通过查表获得。在 t/D 1.6，且Re 4000时，C近似等于 0.98，其中雷诺数是用孔的直径做特征长，孔 中流体的速度做特征速度求出的。 12 2 2 p1p /f1 2 A A = （8-63） 其中v是名义速度。与方程（8-50）对比后可以发现：

24、 () v 将方程m = vAp代入（8-62），并除以板厚x =t，可得： () x C t = （8-64） （4）用表格数据计算流过纤维垫的多孔介质参数 1/2 1 22 A A Cpf C t 18 FLUENT6.1全攻略 在计算流体流过纤维垫或过滤器这类问题时，除了Blake-Kozeny方程外，还可以用试 验数据进行计算，见下表： 表 8-6 纤维计算试验数据表格 固体材料的体积比 纤维的无量纲渗透率B 0.262 0.25 0.258 0.26 0.221 0.40 0.218 0.41 0.172 0.80 表中B = /a2， a为纤维直径。 的含义与（8-48）式中相同，

25、对于给定纤维直径 和体积比的多孔介质是容易求出的。 （5）用压强和速度的试验数据计算多孔介质系数 在已知多孔介质上的速度与压强降的试验数据时，可以通过插值求出多孔介质上的系 数。下面通过一个例子说明这一点。假设试验数据如表8-7所示。 采用表8-7的数据可以拟合出一条“速度压强降”曲线，其方程为： p = 0.28296v2 4.33539v （8-65） 其中p为压强降，v为速度。这个方程等价于方程（8-46），形式上是速度v的多项 式，对比两个方程可知： 1 0.28296=C2 （8-66） 2 表 8-7 速度、压强降的试验数据 速度（m/s） 压强降（Pa） 20.0 78.0 50

26、.0 487.0 80.0 1432.0 110.0 2964.0 在 =1.225kg/m3时，惯性阻力因子C2 = 0.462。同样地， 4.33539= （8-67） 在 =1.7894105 时，可知1/ = 242282。 19 8. 采用指数律模型 FLUENT6.1全攻略 如果采用指数律模型计算多孔介质的动量源项，则需要输入的系数只有C0和 C1。在 Fluid（流体）面板的Power Law Modal（指数律模型）下，输入C0和C1。 指数律可以与Darcy和惯性模型同时使用，C0必须使用国际单位制。 9. 定义多孔性 在Fluid（流体）面板中的Fluid Porosity

27、（流体多孔率）下设置多孔率，即可定义计算 中的多孔率参数。 定义多孔率的另一个方法是使用 UDF 函数。在创建了相关函数并将其载入 FLUENT 后就可以在计算中使用了。 10. 定义多孔介质材料 在Fluid（流体）面板中的Fluid Porosity（流体多孔率）下，选择Solid Material Name （固体材料名称），然后直接进行选择即可。如果固体材料的属性参数不符合计算要求，可 以对其进行编辑，比如可以用UDF函数编辑材料的各向异性导热率。 11. 定义源项 如果在计算中需要考虑多孔介质上的热量生成，可以开启Source Terms（源项）选项， 并设置一个非零的Energy（

28、能量）源项。求解器将把用户输入的源项值与多孔介质的体积 相乘获得总的热量生成量。 12. 定义固定值 如果有些变量的值不需要由计算得出，就可以选择Fixed Values（固定值）选项，并认 为设定这些参数。 13. 在多孔区域中取消湍流计算 在Fluid（流体）面板中，开启Laminar Zone（层流区）选项，就可以将湍流粘度设为 零，从而使相关区域中的流动保持层流状态。 14. 定义旋转轴和区域运动 其方法与标准流体区域上的设置相同，这里不再重复。 20 FLUENT6.1全攻略 8.19.7 用物理速度定义多孔介质 在多孔介质的计算中，FLUENT 在缺省情况下用名义速度计算多孔介质区

29、域的流动。 名义速度可以在总体上保证流动的连续性，但是对于流动细节的考虑则不如物理速度，因 此如果要更准确地计算多孔介质流动应该使用物理速度。名义速度与物理速度的关系如下： GG （8-68） 在Solver（求解器）面板中的Porous Formulation（多孔介质公式）中可以定义物理速 v名义 v物理 度。 采用物理速度公式，假定为通用的流体变量，则各向同性多孔介质的控制方程为： ()()() G vS t += + （8-69） 假定流动为单相流，则体积平均的质量和动量应满足的守恒型方程如下： tvG （8-70） ()()( )v v ( )+()=0 GGGGGG G 2 +=

30、+ + v C vvpB f t 2 （8-71） （8-71）式的最后一项是由多孔介质产生的粘性和惯性阻力项。 8.19.8 多孔介质流动的求解策略 在多数情况下，可以用标准的求解次序计算多孔介质流动，但是在计算过程中可以发 现计算的收敛速度大大放慢，特别是通过多孔介质区的压强降很大时更是如此。收敛速度 放慢的主要原因是多孔介质引起的压强降在动量方程中是作为一个源项出现的，而源项则 可以看作对计算增加了一个扰动，使原有格式的收敛性有所降低。 解决收敛性问题的方法是提高初始流场的准确性，初始流场的准确性可以加速收敛过 程。比如在多孔介质前后设置更接近真实压强值的压强条件，或者在计算中先用不带多

31、孔 介质模型的模型计算出初始流场，再把多孔介质模型设置到计算中，这样在一个更接近于 最终结果的初始流场中开始的计算能够更快地收敛。 在各向异性多孔介质的计算中，如果某个方向上的阻力系数远远大于其他两个方向， 也无需将阻力因子设置成很大的值，而最好是将它限制在其他两个方向的阻力系数的2、3 个量级以内。比如多孔介质的一个方向上的阻力系数可能是无穷大，但是在计算中不要将 阻力系数设置成无穷大，而是需要将它限制在主流方向阻力系数的1000倍以内。 21 FLUENT6.1全攻略 8.19.9 多孔介质的后处理 能够反映多孔介质流动特点的参数是速度和压强。这两个参数在FLUENT内建的后处 理模块中都

32、可以找到，并很容易进行图形显示等操作。 需要注意的是多孔介质的热力学变量的计算结果不能反映固体介质的性质。报告中显 示的密度、热容、热导率和焓，都是多孔介质区域中流体的属性，与固体介质无关。 如果多孔介质流动计算中存在表面化学反应，可以用 Reactions.（反应）项中的 Arrhenius Rate of Reaction-n（第n个反应的Arrhenius反应速率）来表示。 8.20 风扇边界条件 在已知风扇几何特征和流动特征的条件下，风扇的这些特征可以被参数化用于计算风 扇风扇对流场的影响。在风扇边界条件中可以输入一条确定风扇前后压力头与速度关系的 经验曲线。同时可以确定的量还包括风扇

33、旋转速度的径向与周向速度分量。风扇模型不是 绕风扇流动的精确模型，但是它可以计算流过风扇的流量。风扇可以与其他类型的源项共 同使用，也可以作为唯一的源项使用。在作为唯一源项时，系统流量的计算是在考虑系统 损失和前面提到的风扇曲线的前提下完成的。 8.20.1 风扇方程 1. 经过风扇的压强跃升模型 风扇的物理模型是一个无限薄的面，压强在经过这个面时出现跃升，而跃升的值是速 度的函数。 N = n 1 （8-72） pfnv 1 n 在函数关系取为多项式的情况下，压强跃升函数的形式为： = 式中p为压强跃升量， fn为多项式系数，v为速度。需要注意的是在这个关系式中， 速度v的值应该保持为正值，

34、因为这代表流体是向前流过风扇，所以才会产生压强跃升。 另一个确定压强跃升的方法是用垂直于风扇的质量平均速度作为计算压强跃升的速度，算 出风扇区域的压强跃升。 22 18.3.6 连接多相的UDFs到FLUENT 1 DEFINE_CAVITATION_RATE 一旦通译或编译了DEFINE_CAVITATION_RATE UDF，函数的名字将会在FLUENT的 User-Defined Function Hooks面板（如图18.3.6.1所示）中变地可视与可选。请注意，气蚀 率的UDFs仅适用于混合物的多相模型。 要连接UDF到FLUENT中，首先，在Multiphase Model面板（如

35、图18.3.6.2）中激活 Mixture（混合物）模型，然后在Interphase Mass Transfer下面选择Cavitation。 Define Models Multiphase 下一步，打开User-Defined Function Hooks面板， Define User-Defined Function Hooks 图18.3.6.1：User-Defined Function Hooks 面板 图18.3.6.2：Multiphase Model面板 在Cavitatin Mass Rate Function的下拉菜单中选择函数的名字（例如，user_cavitation

36、_rate）， 然后点击OK。 2 DEFINE_EXCHANGE_PROPERTY 一旦通译或编译了 DEFINE_EXCHANGE_PROPERTY UDF，函数的名字将会在 FLUENT的User-Defined Function面板（如图18.3.6.4所示）中变地可视与可选。 用户化的质量输运UDFs可以用于VOF，混合物和Eulerian多相模型。阻力系数UDFs 可以用于混合物和Eulerian模型，而热传递和升力系数UDFs仅适用于Eulerian模型。在连 接函数之前需要激活多相模型。 要连接一个交换特性的 UDF 到 FLUENT，首先需要在 Phases 面板中通过点击 Interactions.打开Phase Interaction面板（如图18.3.6.3所示）。 Define Phases 下一步，在Phase Interaction面板中点击合适的标签（如，Drag），然后从所需要的相 应交换特性（如，Drag Coefficient）的下拉菜单中选择user-defined。这样将打开User-Defined Functions面板。 Define Dynamic Mesh Parameters 图18.3.5.3：Phase Interac