偏滤器极限热流数值预计方法对比研究.docx
《偏滤器极限热流数值预计方法对比研究.docx》由会员分享,可在线阅读,更多相关《偏滤器极限热流数值预计方法对比研究.docx(11页珍藏版)》请在三一文库上搜索。
1、偏滤器极限热流数值预计方法对比研究夫青山(武夷山大学)摘要:临界热流的数值模拟是新型偏滤器冷却结构设计优化的有效手段,而目前尚无成熟有效的水冷偏滤器临界热流数值预测方法,针对此问题,本文对比研究了均相流沸腾模型和两流体沸腾模型对水冷偏滤器临界热流进行数值模拟的适用性;通过两种模型含气量分布、壁面温度、壁面平均换热系数、临界热流值等计算结果的对比,发现针对本文的研究对象和研究参数范围,均相流沸腾模型对于临界热流值的计算普适性较好,对于不同工况的临界热流计算误差均不超过20%;两流体沸腾模型对于高压力工况(Pin=4.2Mpa)临界热流计算误差小于5%,但对于低压力工况(Pin=0.2Mpa)计算
2、误差达到40%;在两相体积分数和壁面换热系数分布上两流体模型计算结果比均相流模型更合理。关键词:水冷偏滤器;临界热流;数值模拟;沸腾换热中图分类号:TL62+6 文献标识码:A1 引言偏滤器是托卡马克聚变装置中排除杂质和输出反应能的关键部件,目前大多采用水冷,其临界热流(Critical Heat Flux,CHF)是影响聚变装置安全性的重要限制性参数。对于偏滤器CHF的预测国内外学者通过理论和实验研究,提出了许多半机理模型和经验关系式1-3。这些模型和关系式对其试验范围内的工况和换热结构适用性强,而对于新型换热结构的偏滤器临界热流预测准确度不足,数值模拟CHF预测方法是新型偏滤器冷却结构设计
3、优化的有效手段4。高热流管流中临界热流发生在过冷沸腾状态下,对过冷沸腾的数值模拟目前主要有两种方法,分别为采用欧拉均相流沸腾模型5和两流体沸腾模型6进行数值计算。欧拉均相流沸腾模型采用平均化的方法对两相进行处理,处理方法较为简单,适用性较好,对于过冷核态沸腾阶段和过冷膜态沸腾阶段均可以进行计算,但无法清楚表现出两相间的作用。郭雷、张树生等人7在使用该模型对竖直矩形细通道内沸腾换热进行计算时发现,模型可以模拟出典型临界热流现象,但准确性不足。对于两流体沸腾模型,该模型分别对两相进行求解,并通过相间作用将两相联系起来,对两相间的作用描述较为清楚,可以完整描述过冷沸腾的各种流动和换热机制8,其计算准
4、确性取决于对流动换热机制描述的正确性。Eckhard Krepper等人4和Y. Bournonville等人9在对不同的过冷沸腾工况进行数值模拟时,发现其对于过冷核态沸腾阶段的计算比较准确,但无法有效模拟出临界热流发生时的壁面温度突升现象。其原因主要在于临界热流发生前后管内的流动和换热机制发生变化,模型中缺乏对不同流动换热机制的处理方法。为了解决这一问题,Ioilev A,和Samigulin M等人10提出了采用分布函数来对不同流动机制进行描述的方法,但其对于偏滤器临界热流进行计算的可行性尚无具体研究。为获得水冷偏滤器临界热流数值预测方法,本文将对比研究均相流沸腾模型和两流体沸腾模型对水冷
5、偏滤器临界热流进行数值模拟的适用性;通过两种模型含气量分布、壁面温度、壁面平均换热系数、临界热流值等计算结果的对比探索偏滤器临界热流数值计算方法,为工程设计提供技术支持。2 数值模拟方法2.1 计算模型为便于与试验数据进行对比,本文对Ronald 11等人的水冷偏滤器试验模型进行计算,偏滤器模型如图1所示,计算工况如表1所示,其中工况1是Ronald 12等人的试验工况,其压力较低,为了更全面比较两种模型的适用性,令取高压工况进行计算,即工况2。计算中,采用速度进口和压力出口,并假设出口边界两相处于同一个压力场。加热壁给定热通量。壁面对液相水设定为无滑动边界条件。增加100mm长的进出口段,免
6、除进出口效应的影响。网格划分如图2所示,为消除计算结果对网格的依赖性,进行网格独立性试验,结果表明网格数加密至100万时,数值解稳定。图1 偏滤器实验模型图 图2 CFD数值计算网格划分表1 计算工况计算工况出口压力P/MPa进口流量v/ (ms-1)热流密度q/(MWm-2)进口水温t/10.2071.180.21.2 2624.541102002.2 均相流沸腾模型均相流模型将流动的气液两相看作一种具有特殊物性均匀介质,其物性参数取两相平均值。相变的实现通过在控制方程加入质量源项和能量源项实现。对于质量传递源项,使用Lee模型12来表示两相间的质量传递: mv=vllTl-TsatTsat
7、,TlTsat 0 ,TlTsat (2)其中,ml和mv分别为单位时间单位体积内气液两相间质量传递,kg/s m3;l和v分别为液相和气相的体积分数;Tl、Tv和 Tsat分别为液相温度、气相温度和当地饱和温度,K;l和v分别为液相和气相的时间松弛因子,为保证计算的准确性,采用Hertz -Knudsen公式13来计算其值。根据分子动力学理论,交界面上的表面蒸发冷凝质量流量可以表示为:F=M2RTsat(P*-Psat) (3)其中F为表面质量流量kg/s m2;P*为气体压力Pa;Psat为饱和压力Pa;Tsat为当地压力下水饱和温度K;R为气体常数;是与气体物理性质有关的调节因子;M是气
8、体摩尔质量,对于水是18g/mol。根据Clapeyron-Clausius公式,将压力和温度关系带入4式,可得表面蒸发质量流量和(Tl-Tsat)间的关系为 13:F=M2RTsatL(gll-g)Tl-TsatTsat (4)将该式与4式对比得:=6dM2RTsatLll-g (5)其中d为气泡直径m;L是当地压力下水的气化潜热,J/kg;l和g分别是是液态水和蒸汽密度,kg/m3。对于能量源项Sh,采用与质量相关的能量传递方程: e=Lm (6)其中,L为当地气化潜热,其值为当地压力下饱和水蒸汽和饱和液态水之间的焓差;m为气液转化质量源项,其形式为公式1、2。可见均相流沸腾模型对两相的处
9、理较为简单,对于不同的流动换热机制采取相同的处理方式,且气液转化效率只与当地过冷度和时间松弛因子有关,因此其适用性较好,对于过冷沸腾的不同阶段(如核态沸腾阶段、膜态沸腾阶段)均可以进行计算,且计算准确性主要取决于时间松弛因子,这也造成了其无法准确描述两相间的作用,计算准确性有限。2.3 两流体沸腾模型两流体模型假定流场中的气液两相均为相互穿透的连续介质;两相在时空上共存并满足各自的控制微分方程;气液两相间动量、能量以及质量的传输通过相间传输项耦合,主要包括相间热质传递模型和相间动量传递模型。2.3.1 相间热质传递根据RPI模型6,通过固体壁面传递给流体的热通量qW可以分为3部分:液相对流传热
10、qC,由于液体冲击气泡脱离处壁面导致的激冷传热qQ和蒸发传热qE;但在临界热流条件下,壁面将不再被液体润湿,而被气体占据,此时能量直接由壁面传入气体。这样,壁面传热机制发生改变,RPI模型将不再适用。此时壁面热流还有一部分传入气体,故使用下述函数来对壁面热流分布进行描述10:qw=qC+qQ+qEfl+1-flqV (7)其中:qC=hc(Tw-Tl)(1-Ab) (8)qQ=2kllT(Tw-Tl) (9)qE=dw3Nwvifgf/6 (10)qV=hV(Tw-TV)Ab (11)其中hc和hV分别是液相和气相对流换热系数,W/(m2k);Tw和Tl分别为壁面温度和液体温度,K;Ab=mi
11、n(1,KNwDw24)是壁面上气泡所占面积比;kl是液相导热系数,W/(mk);T是时间周期,s;l=kllCpl是扩散系数,m2/s;dw是由Tolubinsky关系式4决定的气泡脱离直径,m;Nw=C1.085(Tw-Tsat)210是气化核心密度,v是气体密度,kg/m3;ifg是产生蒸汽的汽化潜热,J/kg;f是气泡脱离频率,f=1T=4g(l-v)3lDw。函数fl决定了壁面的热流分布,其表达式10为:fv=1-fl=120 vv,2 (12)其中气体体积分数v和临界气体体积分数v,1和v,2共同决定了壁面热流分布,一般情况下,认为v,1=0.9,v,2=0.95。这样,近壁单元中
12、存在的液体蒸发传质速率由蒸发热流决定:mE=qEifg+Cp,lTsub (13)在液体主流区,从壁面脱离的气泡被冷凝,其传质、传热速率可以表示为:ml=max(hlvTsat-TlAlvHlv,0) (14)Qlv=mlHlv (15)式中,hlv为由Ranz-Marshall关系式决定的相间传热系数,Alv=6/db,其中db为气泡直径,为空泡份额。2.3.2 相间动量传递过冷沸腾汽液两相间的动量传输通常表示为界面力的形式,包括曳力FD,升力Flift,湍流耗散力Ftd,壁面润滑力Fwl。两相间总的作用力Flv为以上各力之和:Flv=FD+Flift+Ftd+Fwl (16)其中各力的具体
13、表达参考文献4。2.3.3 流动机制转变的描述Ioilev A和Samigulin M等人10指出,在临界热流状况下,管内流动机制发生变化,由气泡流转为弥散流,液体由连续相转为离散相,气体由离散相变为连续相。随着流动机制的转变,其界面面积,动量传递项(曳力,升力,湍流耗散力等),热量传递项等都发生变化。为了计算不同流动机制下的界面传输项(动量交换、能量交换、界面面积等),采用分布函数fv来描述。以表示界面传输项10,则:=1-fvbubbly+fvdroplet (17)其中fv与式12表达式相同,但v,1=0.3,v,2=0.7。bubbly和droplet分为为气泡流和弥散流条件下计算得到
14、的界面传输项。当液体被定义为主相后,在计算流动机制转化过程中,这是保持不变的。但是在计算bubbly时液体被当作主相,而计算droplet时气体被当作主相。因此,两流体沸腾模型对两相间的作用描述较为清楚,可以完整地描述过冷沸腾中的各种流动和换热机制。但该模型十分复杂,需要对各种流动换热机制了解十分清楚,其计算准确性也取决于对流动换热机制描述的正确性。3计算结果及讨论3.1两种模型含气量计算结果对比图3为两种模型在工况1下沿程和不同截面含气量计算结果的对比,其中图3a和3b为均相流模型计算结果,图3c和3d为两流体模型计算结果。由图可以看出,对于沿程含气量分布的计算,两种模型计算结果大致相似,都
15、较为不均匀,含气量变化趋势一致;对于周向含气量分布的计算结果差距较大,均相流模型的计算结果较为不均匀,且冷却通道四周都有气体分布,而两流体模型计算结果较为均匀,气体主要分布在冷却通道上表面附近,且仅分布在离壁面很近的区域内。其原因在于:(1)均相流模型没有考虑相间作用力和气相对于主流湍流流动的影响,只考虑了表面张力,气体的分布主要取决于当地液体温度,因此其周向分布随即性较大,分布较不均匀;(2)两流体模型对于两相间的作用考虑较为全面,气相受到液相多种作用力的影响,分布较为均匀,且由于重力的作用,气体主要集中在冷却管道上表面,冷却管道下半部分没有气体分布。a. 均相流模型轴向含气量计算结果b.均
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 滤器 极限 热流 数值 预计 方法 对比 研究
链接地址:https://www.31doc.com/p-14859591.html