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

《应用两流体模型的偏滤器极限热流预计方法研究》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、应用两流体模型的偏滤器极限热流预计方法研究土海滨(武夷山人学)摘要：临界热流预测准确性对水冷偏滤器设计十分重耍。针对两流体模型在过冷沸腾计算中无法对膜态沸腾阶段进行冇效计算的问题，采用含汽呆分布函数的方法来描述气泡流、搅样流和弥散流等不同阶段两相流动机制的变化，从而提高不同流动机制下两相界面传输和壁面换热的计算准确度，以实现模型对过冷膜态沸腾阶段的冇效计算。改进前后模型的计算结果农明，改进諦的两流体沸腾模型无法模拟出热负荷达到临界热流时的壁面温度突升现彖，而修」E厉的模型有效模拟出了这一现彖；通过与Ronald等人的试验数据对比，计算得到的临界热流值误差约为40%。关键词:水冷偏滤

2、器；临界热流；数值模拟；沸腾换热；两流体模熨中图分类号：TL62+6文献标识码：A1引言偏滤器是托卡马克聚变装置中排除杂质和输出反应能的关键部件，冃丽大多采用水冷，其临界热流(CriticalHeatFlux,CHF)是影响聚变装置安全性的重要限制性参数。对于偏滤器CHF的预测国内外学者通过理论和实验研究，提出了许多半机理模型和经验关系式⑴习。这些模型和关系式对其试验范围内的工况和换热结构适用性强，而对丁•新型换热结构的偏滤器临界热流预测准确度不足，数值模拟CHF预测方法是新型偏滤器冷却结构设计优化的有效手段⑷。高热流管流中临界热流发生在过冷沸腾状态下，对过冷沸腾进行数值模拟目前主

3、要冇两种方法，分别为采用两流体沸腾模型⑸和欧拉均相流沸腾模型⑹进行数值计算，其中两流休沸腾模型可通过其自定义子模型完整地描述过冷沸腾状况下的各种复杂机制，因此被广泛地应用丁•各种规则截而管道内沸腾两相流的计算与分析⑺。而EckhardKrepper等人⑷和Y.Bournonville等人⑻在对不同的过冷沸腾工况进行数值模拟时，均发现现冇模型对于过冷核态沸腾阶段的计算比较准确，但无法对过冷沸腾膜态阶段进行有效计算。针对该问题，木文采用含汽量分布函数的方法〔9］来描述气泡流、搅拌流和弥散流等不同阶段两相流动机制的变化，从而提高不同流动机制下两相界而传输和壁而换热的计算准确度，以实现模型

4、对过冷膜态沸腾阶段的有效计算。并拟采用Ronald等人［⑼的低热流试验数据对改进方法的准确性进行验证，以期为新型水冷偏滤器设计提供冇效的临界热流预测方法。2过冷沸腾两流体模型及修正2.1相间热质传递根据RPI模型［习，通过固体壁面传递给流体的热通量可以分为3部分：液相对流传热心，由于液体冲击气泡脱离处壁面导致的激冷传热@Q和蒸发传热心，即:其中:qc=hc(Tw-TlXl-Ab)Qe=^dlNwpvifgf/6其中短是单相对流换热系数，W/(m4g(pi-pJy］k)；心和”分别为壁面温度和液体温度,Ab=min(1,K竺加)是壁面上气泡所占面积比；俭是液相导热系数，W/(m・k)

5、；4T是时间周期，s；入i=kJpiCpi是扩散系数，m2/s；如是由Tolubinsky关系式［⑴(5)决定的气泡脱离直径，m；Nw=C1Q8TW-Tsaty力用兀，湍流耗散力鵝，壁面润滑力丘1。两相间总的作用力耳为以上各力z和:是气化核心密度，几是气体密度，kg/m这样，近壁单元中存在的液体蒸发传质速率由蒸发热流决定：mE=，fg十CplATsub在液体主流区，从壁而脱离的气泡被冷凝，其传质、传热速率可以表示为:rht=max严(皿-门)化o)HlvQiv=miHlv(7)式屮，幅为由Ranz-Marshall关系式决定的相间传热系数,4山=6a/db,其屮心为气泡直径，a为

6、空泡份额。2.2相间动量传递过冷沸腾汽:液两相间的动量传输通常表示为界面力的形式，包括曳力岛，升；ifg是产生蒸汽的汽化潜热，J/kg；f是气泡脱离频率，/=^=其屮各力的具体表达参考文献［ll］o2.3模型修正在壁而沸腾屮，临界热流现彖发生的特征是换热系数的急剧下降和舉而温度的迅速上升。在临界热流条件卜：壁面将不再被液体润湿，而被气体占据，此时能量直接由壁面传入气体。这样，壁面传热机制发生改变，式1将不再适用。另外，气体占据壁面使壁面附近的气相由离散的变为连续的，壁面附近的流动机制发生改变，这导致基于气泡是离散的这一前提提出的对各个物理量进行计算各种经验关系式不再适用。此时，需对

7、两流体沸腾模型进行修正。(1)壁面热流分配模型修止由于临界热流发生后，壁面人部分被气体占据，使用下述函数來对壁面热流分布进行修正〔典Qw=(Oc+@Q+0E)f(如)+(1-f(m))@v(9)其中Qc、0Q和蚣分别为液相对流热流、激冷热流和蒸发热流。Qv是由涨面传入气体的热流，该值的计算与液体对流热流心类似，其传热系数如通过壁面函数计算。函数门⑷)决定了壁面的热流分布，其表达式〔⑺为：其中气体体积分数⑰和临界气体体积分数Q",］和⑰,2共同决定了壁面热流分布：当时，