应用两流体模型的偏滤器极限热流预计方法研究.docVIP

应用两流体模型的偏滤器极限热流预计方法研究 王海滨 (武夷山大学) 摘要：临界热流预测准确性对水冷偏滤器设计I?分重要。针对两流体模型在过冷沸腾计算中无法对膜态沸 腾阶段进行有效让算的问题,采用含汽最分布函数的方法來描述气泡流、搅拌流和弥散流等不同阶段两相 流动机制的变化，从而提高不同流动机制下两相界面传输和壁面换热的计算准确度，以实现模型对过冷膜 态沸腾阶段的有效计算。改进前后模型的计算结果表明，改进前的两流体沸腾模型无法模拟出热负荷达到 临界热流时的壁面温度突升现象，而修正后的模型有效模拟出了这一现象；通过与Ronald等人的试验数据 对比，计算得到的临界热流值误差约为40%。 关键词：水冷偏滤器：临界热流：数值模拟：沸腾换热：两流体模熨 屮图分类号：TL62+6文献标识码：A 1引言 偏滤器是托卡马克聚变装置小排除杂质和输出反应能的关键部件，H前大多 釆用水冷,其临界热流(Critical Heat Flux, CHF)是影响聚变装置安全性的重 要限制性参数。对于偏滤器CHF的预测国内外学者通过理论和实验研究，提出 了许多半机理模型和经验关系式〔口］。这些模型和关系式对其试验范I韦I内的工况 和换热结构适用性强，血对于新型换热结构的偏滤器临界热流预测准确度不足， 数值模拟CHF预测方法是新型偏滤器冷却结构设计优化的有效手段⑷。高热流 管流屮临界热流发生在过冷沸腾状态下，对过冷沸腾进行数值模拟n前主要有两 种方法，分别为采用两流体沸腾模型⑸和欧拉均相流沸腾模型⑹迸?行数值计算， 其屮两流体沸腾模型可通过其口定义子模型完整地描述过冷沸腾状况下的各种 复杂机制，因此被广泛地应用于各种规则截而管道内沸腾两相流的计算与分析⑺。 而EckhardKrepper等人同和Y. Boumonvillc等人凶在对不同的过冷沸腾工况进行 数值模拟吋，均发现现有模型对于过冷核态沸腾阶段的计算比较准确，但无法对 过冷沸腾膜态阶段进行有效计算。针对该问题，木文采用含汽量分布函数的方法 切来描述气泡流、搅拌流和弥散流等不同阶段两相流动机制的变化，从而提高不 同流动机制下两相界而传输和壁而换热的计算准确度，以实现模型对过冷膜态沸 腾阶段的有效计算。并拟采用Ronald等人的低热流试验数据对改进方法的准 确性进行验证，以期为新型水冷偏滤器设计提供有效的临界热流预测方法。 2过冷沸腾两流体模型及修匸 2.1相间热质传递 根据RPi模型⑸，通过固体壁而传递给流体的热通量q”可以分为3部分：液相对流传热0C，由于液体冲击气泡脱离处壁面导致的激冷传热0Q和蒸发传热厲，即：其中:= Qc 根据RPi模型⑸，通过固体壁而传递给流体的热通量q”可以分为3部分：液相对流传热0C，由于液体冲击气泡脱离处壁面导致的激冷传热0Q和蒸发传热厲， 即：其中: = Qc + 0Q + Qe ⑴ = —G(1 —%) ⑵ 力惴(心—G⑶@e = nd^Nwpvifgf/6 (4) 其屮很是单相对流换热系数，W/(m2k)；心和7｝分别为壁面温度和液体温度，K； % = min (1,K竺匹)是壁面上气泡所占面积比；俭是液相导热系数，W/(m -k)； 4T是时间周期，s；入i = kjpi5是扩散系数，m/s；血是由Tolubinsky关系式［川 决定的气泡脱离直径，m； Nw = C1QQ\TW Tsat)21Q是气化核心密度，几是气体 密度，kg/m3； ifg是产生蒸汽的汽化潜热，J/kg； /?是气泡脱离频率，/ = ^ = 4g(pi-pM寸 3pQw 这样，近壁单元屮存在的液体蒸发传质速率由蒸发热流决定: Qe ifg + Cp,lA7sub (5) 在液体主流区，从壁面脱离的气泡被冷凝，其传质、传热速率可以表示为: mi = max (塑如匚处,o)Hlv (6) Qtv = (7) 式小，尬为由Ranz-Marshall关系式决定的相间传热系数,4山=6a/db,其屮为 气泡直径，ex为空泡份额。 2?2相间动量传递 过冷沸腾汽液两相间的动量传输通常表示为界而力的形式，包括曳力忌，升 力用化，湍流耗散力嘉，壁而润滑力為。两相间总的作用力耳为以上各力z和: Av =爲 + + Ed + Avl 其屮各力的具体表达参考文献［ll］o 2.3模型修正 在壁面沸腾屮，临界热流现象发生的特征是换热系数的急剧下降和壁而温度 的迅速上升。在临界热流条件下，壁面将不再被液体润湿，而被气体占据，此时 能量直接由壁面传入气体。这样，壁面传热机制发生改变，式1将不再适用。另 外，气体占据壁面使壁面附近的气相由离散的变为连续的，壁面附近的流动机制 发生改变，这导致基于气泡是离散的这一前提提出的对各个物理量进行计算各种 经验关系式不再适用。此时，需对两流体沸腾模型进行修正。 (1)壁面热