数值模拟步骤

1、数值模拟1、CFD方法简介利用CFD方法，采用流体力学分析软件 Flue nt对三相分离器的流场进行了研究 和分析，为实验研究提供理论支持。CFD是英文Computational Fluid Dynamics（计算流体动力学）的缩写,是一门用数值计算方法求解流动主控方程以发现各种流动现象规律的学科。用CFD技术进行 数值求解的基本思想是：把原来在空间与时间坐标中连续的物理量的场，用一系列 有限个离散点上的值的集合来代替，通过一定的原则来建立离散点上变量值之间关 系的代数方程，求解代数方程以获得所求解变量的近似值。 其主要用途是对流态进行 数值仿真模拟计算，因此，CFD技术的用途十分广泛，可用于

2、传质、传热、动量传递 及燃烧等方面的研究。流体机械的研究中多用CFD方法对分离器进行仿真模拟，其基本应用步骤如下:1）利用Gimbit进行前处理a. 根据分离的形状、结构及尺寸建立几何模型;b. 对所建立的几何模型进行网格划分；2）利用Flue nt进行求解a. 确定计算模型及材料属性;b. 对研究模型设置边界条件;C.对前期设置进行初始化，选择监视器，进行迭代计算;3）利用Flue nt进行后续处理，实现计算结果可视化及动画处理。上述迭代求解后的结果是离散后的各网格节点上的数值，这样的结果不直观。因 此需要将求解结果的速度场、温度场或浓度场等用计算机表示出来，这也是CFD技术 应用的必要组成

3、部分。利用CFD方法进行仿真模拟可以对分离器的结构设计及参数选择作出指导，保证 设计的准确度，也可以为分离器样机的试验提供理论参考。由于CFD仿真模拟的广泛 使用及其重要性，国内外很多学者，如 Mark D Turrell 、M.Narasimha、师奇威等都 对其进行了研究，尤其是 A.F. Nowakowski及Daniel J.SUASNABAR等人对CFD技术在旋流器模拟方面的应用做了详细的介绍，这些工作对CFD技术的发展起到了积极的 促进作用。2、控制方程流体流动要受物理定律的支配，基本的守恒定律包括：质量守恒定律，动量守恒 定律、能量守恒定律。如果流动处于湍流状态，系统还应要遵守附

4、加的湍流输运方程。 1、基本假设天然气的主要成分是是甲烷，故采用甲烷的替代天然气的性质; 在工作状态下，流动不随时间变化，流动为稳态； 水在管内的流动可以简化成二维流动；不考虑温度的影响，服从绝热流动基本方程。(1) 鉴于我国各主力油田采出液含水已达到 80%以上，故以水代替采出液进行分 析计算；(2)(3)(4)2、基本控制方程(1)连续性方程在直角坐标系下的质量守恒方程又称连续性方程：W 0(1-1)Xi(2)动量方程:(UiUj)X XjXj(1-2)2式(1-1，1-2)中U为流体速度，下标i表示方向，上标表示脉动值，上标“-” 代表对时间的平均，卩为动力粘度，P为流体微元体上的压力。

5、C模型，其原理是用紊动能 而流体的有效粘性系数即是流体分子由于RNG k-c模型在工程流场计算中具有诸多优势，故对结构内部的数值模拟 采用RNG- c模型进行。RN(k- c模型是改进的标准k- k和紊动能耗散率c来表示流体湍流粘性系数， 粘性系数和湍流粘性系数之和。 k方程和c方程分别为：(k)t(kUi)XiXjkkUeff GkXj(1-3)()(tUi)Xi*C1XjGkC2(1-4)其中:effk2C ，C 0.0854, k1.39Ci*Ci(1lG 1.42, C21.68(1-5)V k1 ILUiEMJEij寺 0 4.377, 0012RNG- £模型通过修正湍流

6、粘度，考虑了平均流动中的旋转和旋流流动情况， 在 6方程中增加了一项，从而反映了主流的时均应变率 Ej，这样，RNG k-£模型中产 生项不仅与流动有关，而且在同一问题中也还是空间坐标的函数。故 RNG k-£模型 可以更好的处理高应变率及流线弯曲程度较大的流动。3、网格划分借助gambit软件对直管段划分结构网格，对非直管段划分非结构网格。四种结 构的网格划分如图所示：结构二结构一结构三结构四5、模型的选择及边界条件的确定应用多相流中的欧拉模型；控制方程采用 RNG k弋 模型；入口边界条件设置为 速度入口，出口边界条件设置为自由出流；由于入口两相中含夜体积分数极小，湍动

7、 粘度影响远大于重力，固忽略重力的影响；湍流强度设置为10%水力直径为圆管内流 道截面直径。6模拟结果分析6.1结构一的模拟1）入口流速3m/s,液相体积分数分别为：0. 001%、0.003%、0.005% 三种情况气相浓度场分布如图所示：1普SMItIMIhh pr 21-1-'JFVt fs Jn tfim «ubw_三种情况速度场分布如图所示:hnHUN 九-i-il. 'Ms一 95一莓fi-1 -£ -=_£&ii±i±-T iKBBlsie-IZA 二_*三聲童ae叵 ftn-t-=>_-1-=-&#

8、163;«& nbfH-h 三工一二壬 二=#=E-亠三二一？凰 三I昌t 一邑世音匕一宫m fit >-<?>-> =-L=-p 1- -ZKH9-r H >l_Bi=_-i-KO=-£ _-nvl"-BEXr=aMB-w 套克_艾二E一二-w il-kr更三=_二三玉 slpwywF 鱼匡£ 左一T运兰二二二 一¥H云二三-里丐* _島 £_4I=T 三fi n-li-«£-3-IK 二一_睾_y ,A > I t - 1 7 J .1 j c u .二二-ggiCH

9、FUftwTO j-；fl 腳h fc*TWi n：J£.s wSitIZ-tZ-zLiit三T =- h¥*1iX.WI1二 i r Tr £1±n1imE.E r»i1>1s1>1 Bf . 迁J.MV.!S：jiLi £亠i.斗k£TI-I irT_Up - p"3 h：J iRB-iV441：I.Ms*'门K_，MM父 rLjCr i%7:r.-e-'JLntvDC.sziiH-LPtT IJ>£«.pkii. Wtw ii'iMlI IltR-

10、 '.屈-lU匸护- 三*齐 VW冬13 二一匚盂* -冒*二1= -二宇-工二*二 三 Tv-Fimf 巨上-mhiq三出=“1_三苣玉 可二亍兰畫邑si*-Iit=-Kw-T豈4云星壬E=1=9 -=TS-I -E-i-Tr-三vu-fml-厘m詣鼻一一Hli- 耳工T_-=<3三1二二爭二=一 一 兰=工昱工兰亠01= eiwhi mmu?-=£=-4L三T呂更=二£l_i品 二量三I三m*一一兀 一*4 三T>v=兰一亠.0L in.'lt-1!:Iti-wr IMi<< jRfl-rd- <11 13 :2）固定液相浓

11、度为0.003%,变入口速度分别为：1m/s、3m/s、5m/s。气相浓度场 分布如图所示：=hM *1 HlhM» 1|lw.flr.JAMFl <V "袅 鼻 1 *, %哄 <1*7-壬三三二三三二二#'i:山 lwT ±21 iirb.r i三种情况速度场分布如图所示r - ir-1 u 因 團 型 一 曜二i_c_-Q 三百星iihLJ旬 iE昏垦|一量辰饪 £A±>f-l-I-x WK-L$-2_ix=T&1-V =-l-IM 一 豊7-#亀专_£= ItwI-rB-II* thah

12、63;v -a*=t 壬IF R一一 V _ 一 w一学一-二- _ = Tn>- -T> Til-F旦蔭二见 7- m一B誉二：=益三品 VCBM- vt-s-二亘二二旦 Y一胃£富匸一言丘;汙二这一寸亘吕 一一置工一耳二-二- K一尼-二一丘 |寸| 二 £住白亞上_一一 IBtwt垂X二耳戏-号討芒益 i -_1*累艮書孟 -iB&AEIVtt-_9_ ll«tnl l-ls- FIj-ualpgEE盲-1iBfBtf-PmHwtjnrwr-tsuHl-_!豊理ts-l-ls -n-su£ Fg三駅三n雪ti-三臣二nV-I c

13、=一書三_g= B5H-EF =f=iw 菱-丹曽PMmJI -ww二二 M T E-匚吐A"二5龍-L-舍=-in 亍 - - -丄昌= F-zrT LT 舍二 hlkam三円丘一區ti三 卍二*二耳用_ e ;兰士 _!_ -=f仝 r?N-二 国he匾虎_ 三1!亍_ y/Hh-H- -sc-_-_- Bssva-_ -:三：三云 -z3T?-_ 三九旦盘三k -qMFUaJI 全 帮卜三一旦上 二一 三卽旦亍 2=-.? -_ - L-三二哼二 £ki氏"旦可鱼环MLdh-JSLf-一 FLH = lj-*p二 nL丽11- =一一丘尼 斎迅三童匚斗 -a

14、*-=3-£ r三壬尹JN -?-3-f 产二亍-窑-dll 73-7FA 一岸-1=- £-r =* = -c'-£- 巨三尋=_ 一：一一 JrN 一一 tt-=m=-二鼻-T-3-± -= -i - ELkl=rmff -=3-< kl-L-lfr-6.2结构二入口流速ru上Ni 115 Ml. 7i rwpt'3m/s,液相体积分数分别为：0. 001% 0.003%、0.005%。三种情况气相浓度场分布如图所示：PF1K-F HRil/u-rWrl； jjrfL_lEhtf3|A iLlHia r0ir 广 Al"

15、; FbtfnTf 釘血 hM-iU HAM*三种情况速度场分布如图所示waI >；tf-i'i-r>I i M-h-I -tn I Ft.-*H- I I工鵲藝陆;II Uab. I> J Mi"? ： Tnrf1: "n 'H1 ：>i ；工!_IF Vf)4_ri" t1 q讣.L訂专-thKaripnwt- ix>Itah-oc Jdl IFiAWr 1 j I 注MMW1 订 4 iEhE By i-CKtyrvTtJK- iptwf- Illi 4*'fi榊*>*»村 riUPfr

16、心 ijdi irim 無*阿“护r*-« 201 fcPlUVHT ( J -|5i ifcMd.2）固定液相浓度为0.003%,变入口速度分别为：1m/s、3m/s、5m/s。气相浓度场 分布如图所示：ME"二眄二二二 二：”Bal'd:*Znxn ri W-p TH：- rl-n 'ta J-.ZI'lF.巧 r -'. n4|iv-门4r切 IJ 三种情况速度场分布如图所示I ：-J：IIII Ln IM'仁Z'P 1 H：iT：bra zi'a -fcwT.JM-k> HI：!_- -tr .鼻hr 1

17、1 _町s- m -!-I：»3Nd" 11 :2li .3 JV-JII :坯D fej-11 汕 lill- 耐y I"!齐ipwtt- >j in* *'FJKTi JiMi irim:pbnv-i-vWvct' Gjta* Wa ciy Hta|fflvJ-i i|taH- Li iig M.31ItVjfM 4恤巒 *. fl自F f 刖ipg i-rHEiy 匚EmlDj '/rb lEfHi-''1 miUn HiHTFU口! 机hJV £Wi1 ww»6.3结构三1）入口流速3m/s

18、,液相体积分数分别为：0. 001%、0.003%、0.005%。三种情况 气相浓度场分布如图所示：IPi i世-»l-d谢冠-存 K二匸二-二匚；gJMtTSfc_y37&fgjr:TLrrWcik- ,.' _ -= r_Z_a二-二 : 'vHrFci.li- 則¥比 AhFrr Iw »i- IJC厂,3淀尸.31 b- 1三种情况速度场分布如图所示- I,51” H1 -f- /筈-u！II:41-li11|siHN-.h；-i- -,比1 >!.- -,电！ I a B»r ； .* 三丄壬 _I ' _

19、*丄二-*:-Z". L J二:圭NAgrv 丄 wn -；L- - hB>17rS".L-» - I , I 一一=E .J J. 士 - r= - >- -二 ； -宀 二二 七于宀iu寸二；二r="F.r7：尹王：-一宀Un n,l1町UV桥4*|：4宀iHhiiiZi |*HHH i-a- 5j.h宁"宀一_ 亠 _. . _；:S- C匸二二九疗8 -二"-二一一.- -二 5eH 升 .尸亍;,r"," 我 k 1= ；. 贰、=Lnr»g2riii<w< = 宀 r.

20、 =2詳护步站:： + =： =：1- ! . - - - J - r 'rTn, - - -_ .;匚FT三匸 .-. F占冇七丸/*十*f氧；. .f 三 = ”=.,=："：_rsg-.子"fXE 吋" m- .- .='h 畑八.- : F-4- rl-*nn| ihi2）固定液相浓度为0.003%,变入口速度分别为：1m/s、3m/s、5m/s。气相浓度场分布如图所示:二二二二二二二上匕伫=三种情况速度场分布如图所示:二I-,量:-: E - I1. ” n:tk-X-u-x"BkrF-H-E-n-u；- Mm 耆壬云曲心-：于

21、 -4- 4 * 玄-卄；_=严 r ：.-一-=. I. " I 1 P二-號免?=宾牛._ . 2_： -心d'SM “:-L-_- .-：花蒔；5：£再：口：匸 .绘 T匕5»yJ T; -r J- - -.-.,严.4 j!.v=N»T2r屮 * 25.-'.r *卢”. i-r眄；_11亠“- , 片- -L L-|V=_W-LL :- r 匕 L 一- mb- IH 二丁；.二T- ” ”' -，: ±1'- J.:nSi«'宀5“ '三'： 护=py_ ' - " .V F U ： *-1 - - I * > I I -. I''> ." L，. ,«-«-_!-.- . -1'jc*I -r _一一'. -I.4 J F_I. - r Hp- ,r 厂 p PJ-M.-J U _尸| !>, k| ILP r _T_ruu-L-；逻卡善主芒7TV=；6.4结构四1）入口流速3m/s,液相体积分数分别为：0. 001%、0.003%、0.005%。三种情况气相浓度场分布如图所示：.7 .hI IIHruCr.&i1曲