大功率汽轮机湿蒸汽级两相凝结流动的数值模拟与分析.pdf

20 1 0年第2期 圣 言 芘粉。 II D o N G FA N G T U R B IN E 1 7 巫志华吴其林赵世全 (东方汽轮机有限公司 ，四川 德阳 ，618000) 摘 要：采用真实气体平衡态计算模式和非平衡态自发凝结两相流动计算模式两种不同的计算方法对某大功率汽轮机低 压末两级内三维湿蒸汽两相凝结流动进行 了模拟和分析。计算表明非平衡效应引起湿蒸汽级组内流量、焓降在各级之间 的分配 、叶栅出 口气流角、各级反动度都有发生变化 。 关键词：大功率汽轮机；平衡态；非平衡态；湿蒸汽级；两相凝结流动 N um e ri c a l S i m ul ati o n a nd A na l ysi s of T w o-p ha se C o nd e nsatio n F l ow s i n th e W e t S te a m S ta g e s of a H i g h P o w e r T u rb i ne W u Zhi hua， W u Q i l i n， Z hao Shi quan (D ongfang T urb i ne C o L td D eyang Si c huan 6 1 8000) A b stra c t：W i th the equ i l i b ri u m c o n den sati o n m od el an d no n eq ui l i brium c o nd en sati on m o d el ，th e nu m eri c a l si m u l ati on an d anal y si s of tw o-p hase w et steam flow s i n the l ast tw o stages o f a hi gh p ow er tu rb i n e w as p erfo rm ed T h e resul ts sho w ed that the no n,-equ i l i b ri u m eff ec t l ead ed to th e v ari ati o n o f the m ass fl o w rate o f the turb i ne stag es,， th e en th al p y di stri b uti o n， th e ou tl et fl o w ang l e o f c asc ad e an d the reac ti on deg ree o f th e turb i n e stag e K e y w o rd s ： h i g h p ow er steam tu rb i n e， eq ui l i brium ， no n-eq u i l i brium ， w et steam stage ， tw o-p ha se c o nd ensati on flow 0 引言 湿蒸汽两相凝 结流动广泛存在于火 电和核电 汽轮机 中。湿蒸汽两相凝结流动 的出现一方面降 低 了流动 的效率；另一方面水滴对叶片冲蚀将 导 致结构安全 问题。虽然湿汽 问题对汽轮机的设计 运行 非常重要 ，但三维湿蒸汽 的凝结流动 由于存 在相变而导致其精确的实验测量和数值模拟都 比 较 困难 ，在实际设计过程中，常采用 由Baumann于 1921提 出的经验公式 来考虑湿 汽损失 。这一经 验公式并没有坚 实的理论支持而可能无法满足现 代大功率汽轮机 设计 中以高效率和高可靠性 为 目 标的要求 。在现代大功率火 电汽轮机中，湿蒸汽 级的做功能力能 占整机 的10以上 ，而对压水堆 核 电汽轮 机，进 口蒸汽 则为饱和湿蒸汽。鉴于 当 前汽轮机仍在电力工业 中占据着主导的地位 ，对 湿蒸汽凝结流动进行广泛深入 的研究是很有必要 作者 简介：巫志 华 (1978一 ) ，男，博士，2006年毕业于西安交通大学能源与动力工程学 院，主要从事 汽轮机 的设计研 发工作 学兔兔 w w w .x u e t u t u .c o m 1 8 f索 言 芘 粉拨 I D ON G FA N G T UR B IN E 的，国外许 多透平制造企业联合大学、研 究所对 湿蒸汽凝结流动进行 了大量的研究 ，并取得 了不 少成果 。 上述研究中，对湿蒸汽凝结流动的研究工作 主要集中在喷管和单列叶栅上。从公开发表的文 献来看，对汽轮机级组环境下的三维两相凝结流 动试验研究和数值模拟都很少 ，文献 5中对一台 小功率汽轮机末三级 内的流动进行 了三 维数值分 析 ，比较 了有凝 结与无凝结两种流动计算结果， 但该文无凝结的计算模型中，工质性质 的计算并 没有采用真实气体性质模型，而在大功率汽轮机 设计和计算中，对湿蒸汽凝 结流动通常采用平衡 态计 算模型来进行 ，为进一步分析有凝 结与无凝 结两种流动之间的差异 ，本 文就平衡态模式 以及 非平衡凝结模式两种不 同的计算方法对某大功率 汽轮机低压末两级的湿蒸汽 两相凝结流动进行了 数值模拟，比较 了两种模型计算结果的差异。 数理模型和计算方法 1 1 平衡态计算的数理模型 本文 中采用 的平衡态计算方法是把汽液混合 物当作单相流体处理 ，这样流动 的控制方程就是 通 常的Navi erStokes (NS)方程 。在 直角坐标系 下，这一方程组的守恒形式可以写成： 一0 U +一O E c +堡 +丝： 堡 +堡 +亟f11 at 8x 8v a z bc 8v 8Z E p u p “ 十P p u v p u w pH u p v p uv l p +P p p H v G o= 0 ， c G = I T g p w puw p v w p w P p H w 20 1 0年第2期 其中，P 表示混合物密度，“ ，v，W 为速度分量， P 是静压， ， 为应力张量，E =p e+O ( + 2十 2)表 示总能。为 了使得方程组封 闭，需要补充状态方 程来描述压力 。在完全气体计算模式 中，可以采 用完全气体的状态方程(P =( 一 1)pe)，其中选取 适 当的气体常数和 比热比，如对湿蒸汽 ，通常把 气体 常数取为461 5J (kg K )，比热比取 1 135， 这种方法得到的工质性质与实际湿蒸汽性质间存 在差异 ，在本文 中的平衡态模型计算 中对压力的 描述采用如下表达式： P = p(e， P) (3) 即对平衡态气体的描述不再采用类似完全气 体状态方程而直接从水和蒸汽性质表 中求得 。采 用 比能P和密度 两个独立的变量，水和蒸汽性质 通过IAPWSIF97得到。 1 2 非平衡凝结流动的数理模型 采用文献 7中建立了完全欧拉 欧拉模型来 描述凝结流动 。该模型假设汽液两相 间无滑移， 湿蒸汽两相流动可 以采用两相分开处理 ，其中混 合物 中的蒸汽 为汽相 ，液相 为大量分散在汽相中 凝结水滴 ，相间的作用通过相应 的源项来表示。 气相的控制方程可用通常的可压缩Ns方程加表示 相变引起 的源项来描述 。在直角坐标系下可 以写 成 ： au + + + 警= +警 +警 + c4 + +三 + = + +二 + l d、 p p ,u p P Q ： l p,v I =l p E l 0 l I E =l J I g F = p ,uw p gH u 0 y x q y G = p p p p = l pgV2+Pla =I p p S g = p P p ,H w P m P u m V P W m - pro (hfg-ht) (5) 其中，Pg，P ， ， ，h ， 龟 分别表示气相 密度 、水蒸气混合物 的密度、蒸汽相的总能、应 o I 0 = 学兔兔 w w w .x u e t u t u .c o m 学兔兔 w w w .x u e t u t u .c o m 2 0 l 圣 右 汽 扮 拨 f DO N GFA N G TU R B INE 道 内流动的变化使得动 叶进 口攻角也发生 了明显 的变化。 Z (m )04 1 图2a 次末级中叶展 处压力分布 Z (m) 图2b 末级 中叶展处压力分布 图3给 出了两种计算模型下透平级中各列 叶栅 出口的湿度沿叶 葛 分布 曲线。从 图上可以发现两 种计算模型下 ，次末级 的湿度分布有更为明显的 差别：从次末级静叶栅 出口的湿度分布来看 ，采 用非平衡凝结计算模 型的结果表 明，在叶栅顶部 基本上还没有 出现凝 结；非平衡模型计算得 到的 次末级动叶出 口湿度分布趋势也与平衡态模 型的 有 明显差异 。但末级的静动叶出 口湿度分布趋势 则基本一致，这主要是 由于末级进 口蒸汽 己经带 有一定的湿度所致 。我们知道 ，非平衡凝结过程 中包含了成核和水滴成长两个过程 ，如果非平衡 凝结流动过程 中仅仅包含水滴成长过程 ，则非平 20 1 0年第2期 衡效应的影 响会 比较小，这也意味着如果进 13 是 带有一定湿度的蒸汽 ，则非平衡凝结的影响会减 弱。 图3 叶栅出口湿度沿叶高分布 图4给 出了两级各列叶栅出口气流角沿叶高分 布情况 。从 图上可 以发现两种计算模型下 ，静叶 栅出口气流角相差最大可达1度 以上，这对通流设 计的影响是很明显的。而从图5的级反动度分布曲 线上也充分反映了 由于非平衡凝结引起气流角的 变化进而导致反动度 的变化 。从图5的反动度分布 来看 ，相对于末级而言，非平衡凝结对次末级的 影响更为显著 ，这与前面 的压力分布 以及湿度分 布是一致的。从非平衡凝结 的计算结果来看 ，次 末级的反动度 明显减小，而末级的反动度分布虽 然变化不大，但根部反动度也有所减小，这在设 计 中需要引起重视。 “(Degree) 图4a 次末级静叶栅出口气流角分布 学兔兔 w w w .x u e t u t u .c o m 2O1 0年 第2期 d (Degree) 图4b 末级静叶栅出口气流角分布 B (Degree) 图4c次末级动叶栅出口相对气流角沿叶高分布 图4d 末级动叶栅出口相对气流角沿叶高分布 圣言 乞 扮揍 D O N G FA N G TU R B IN E 图5a 次末级反动度分布 2 1 3 结论 图5b 末级反动度分布 采用真实气体平衡态计算模 式和非平衡态 自 发凝结两相流动模式两种不同的计算模型，对某 大功率低压末两级 内湿蒸汽两相凝结流动进行 了 模拟，分析计算结果可得出如下结论： (1)非平衡效应使透平级的通流能力增大 ，级 的焓降分配也发生了改变 。 (2)多级环境下的非平衡凝结流动比单一喷管 或者单叶栅通道 内的流动要复杂得多，通过喷 管计算或试验得出的结论不能简单应用于多级环 境下的流动。 (下转第27页) 学兔兔 w w w .x u e t u t u .c o m 201 0年第2期 4 结论 本文研究 了透平进 口温度 、压力和环境温度 对有机 朗肯循环输 出功的影响，在给定余热条件 下 ，采 用模拟 退火 算法对 循环 的参 数进行 了优 化 ，得到 了最优运行状态所对应 的参数和功率 ， 得到以下结论： (1) 在本文所使用的有机工质 朗肯循环中， 循环净输出功 随着透平进 口压力 的提高 ，先是增 加 ，而后减少 ，随着透平进 口温度从相应压力下 的饱和温度增加到最高温度 ，循环净输 出功不 断 减小 。环境温度对有机 朗肯循环影响较大，大概 有30 的功率随着环境温度升高而减少 。 (2 ) 在给 定的余热条 件下 ，通过优化 ，得 出：采 用R236fa为工 质 ，循 环 的输 出净 功率 最 大 。 参考文献 1】 H D M adh aw a， H etti arac h c hi ， M i h ajl o G ol ubovi c ， Wi l l i am M Worek ， Y asu yuki Ikegam i O p ti m um desi gn c riteria fo r an O rgani c R an ki n e 圣 言 浅扮接 D O N G FA N G T U R B IN E 2 7 c yc l e usi n g l o w tem p erature geoth erm al heat sourc esJ Energy： 2007： 1 9 2】 U l l i D resc her， D i eter B ruggem ann Fl ui d sel ec ti on fo r the O rgani c R anki ne C yc l e (O R C ) i n bi om ass pow er and heat pl antsJ A ppfi ed T h e rm al e n g i n e erin g ， 2 0 0 7 ， 2 7 ： 2 2 3 2 2 8 3Y C hen， P Lundqvi st， A Johansson， P H atel 1 A c o m parati v e stu dy of th e c arbo n di o xi d e transc ri ti c al pow er c y c l e c o m pared w i th an organ i c ranki n e c y c l e w i th P 123 as w orki ng fl u i d i n w aste heat rec overyJ A ppl i ed them al engineefi ng， 2006， 26： 2142 2147 4刘志刚， 刘咸定， 赵冠春 工质热物理性质计算 程 序 的 编制 及应 用fM 】 北京 ： 科 学 出版 社 ， 19 9 2 5张家荣， 赵廷元 工程 常用物质的热物理性质 手册M 1 北京： 新时代 出版社， 1987 66 康立 山， 谢 云， 尤 矢勇， 罗祖华 非数值并行 算 法(第一册) 模拟退火算法fM 北京： 科学出 版社 1994 -_ _-_ _一 (上接第21页) (3)多级环境下，非平衡效应对叶栅 出口气流 角有 明显影响，使级反动度 发生变化 。在本文考 察的算例中，非平衡效应导致根部反动度减小。 参考文献 【 1】M oore M J et a1 A erotherm odyrnam i c s of Low Pressure Steam T u rb i n es an d C o n den sers U SA ： H em i sp h ere P ub sN ng C orp orati on ，198 7 ： 87 122 【 2】 B akhtar F， Ebrahi m i M ， W ebb R A O n the Perf orm anc e of a C asc ade of T urb i n e R otor T i p Sefion B l ad i n g i n N u c l eating Steam P art 1：S urf ac e P ressure D i stri b utio n s P ro c In sm Mec h E n grs，P art C ，199 5 ， 20 9： 115 124 【 3J B ol m D ， K erp i c c i H ， R en J， Surken N H om ogen