馬為峰, 路 駿, 萬榮華, 韓勇軍

非對稱加熱條件下鍋爐反應(yīng)器蒸汽生成過程參數(shù)影響研究

馬為峰, 路 駿, 萬榮華, 韓勇軍

(中國船舶集團有限公司 第705研究所, 陜西 西安, 710077)

針對魚雷Li/SF6閉式循環(huán)動力系統(tǒng)鍋爐反應(yīng)器復(fù)雜物理化學(xué)過程中參數(shù)耦合、參數(shù)時變造成的系統(tǒng)穩(wěn)定域度小等問題, 文中基于傳熱學(xué)理論, 采用分段式處理方法, 建立了鍋爐反應(yīng)器螺旋換熱管流動換熱模型, 對傳熱過程進行了數(shù)值仿真, 探討了非對稱加熱條件下鍋爐反應(yīng)器入口工質(zhì)溫度、出口工質(zhì)溫度、出口工質(zhì)壓力、出口工質(zhì)流量等參數(shù)的變化對蒸汽生成過程的影響規(guī)律。研究結(jié)果可為系統(tǒng)性能仿真和調(diào)節(jié)控制優(yōu)化, 提升系統(tǒng)穩(wěn)定域度提供支撐, 為鍋爐反應(yīng)器的工程設(shè)計、優(yōu)化和試驗驗證提供理論參考, 對同類裝置的研究和設(shè)計也具有一定的參考價值。

魚雷; 閉式循環(huán)動力系統(tǒng); 非對稱加熱; 蒸汽生成; 鍋爐反應(yīng)器

0 引言

基于Li/SF6燃料的魚雷閉式循環(huán)動力采用Li金屬燃料穩(wěn)定燃燒與高效傳熱生成過熱蒸汽, 驅(qū)動渦輪機做功產(chǎn)生魚雷推力, 采用過熱蒸汽作為循環(huán)工質(zhì), 在實現(xiàn)大功率輸出的同時, 實現(xiàn)了高比能動力系統(tǒng)的閉式循環(huán), 系統(tǒng)無工質(zhì)排放, 性能不受航行深度影響[1]。美國率先研制成功此動力裝置, 并于20世紀(jì)90年代裝備于MK50魚雷。

鍋爐反應(yīng)器是Li/SF6閉式循環(huán)動力系統(tǒng)的重要組成部分, 是Li和SF6發(fā)生反應(yīng)的燃燒室, 也是蒸汽生成的熱交換器, 其材質(zhì)、直徑和管長等對工質(zhì)在管內(nèi)的流動換熱性能影響較大[2]。鍋爐反應(yīng)器的蒸汽生成涵蓋了水工質(zhì)的預(yù)熱、沸騰、氣化和過熱等過程, 涉及熔池兩相流燃燒、螺旋管外流固傳熱以及螺旋管內(nèi)流動換熱, 而螺旋管內(nèi)流動換熱又涉及液態(tài)金屬同螺旋管管壁的換熱, 金屬管壁的導(dǎo)熱, 螺旋管管壁與單相水、單相蒸汽、氣液兩相流的換熱等, 而鍋爐反應(yīng)器的螺旋非對稱加熱條件更加劇了這些復(fù)雜物理化學(xué)過程的參數(shù)耦合, 使得各參數(shù)隨時而變, 極大壓縮了系統(tǒng)的穩(wěn)定域度, 增加了系統(tǒng)的調(diào)控難度。

鄭邯勇[3]對Li/SF6鍋爐反應(yīng)器內(nèi)流體的流動結(jié)構(gòu)與溫度分布進行了研究, 提出了鍋爐反應(yīng)器工作過程的階段劃分方法, 并指出了有利于鍋爐反應(yīng)器穩(wěn)定反應(yīng)與傳熱的條件。王樹峰[4]基于系統(tǒng)效率的理論計算, 提出了鍋爐反應(yīng)器的快速啟動溫度、輸出過熱蒸汽的溫度和壓力、入口水工質(zhì)溫度參數(shù)值。白潔等[5]建立了水下航行器螺旋管蒸發(fā)器穩(wěn)態(tài)一維分布參數(shù)模型, 獲得了工質(zhì)沿流動方向的狀態(tài)參數(shù)。但在鍋爐反應(yīng)器設(shè)計與試驗過程中, 僅通過負(fù)反饋控制方法簡單調(diào)節(jié)SF6流量實現(xiàn)系統(tǒng)調(diào)節(jié)控制, 缺乏蒸汽生成過程的參數(shù)變化規(guī)律支撐, 制約了鍋爐反應(yīng)器的優(yōu)化設(shè)計。

文中針對鍋爐反應(yīng)器螺旋換熱管進行的流動換熱過程研究, 探索非對稱加熱條件下蒸汽生成過程的參數(shù)影響規(guī)律, 為鍋爐反應(yīng)器的工程設(shè)計、優(yōu)化和試驗驗證提供理論依據(jù)。

1 鍋爐反應(yīng)器工作原理

鍋爐反應(yīng)器采用雙圈螺旋管形式, 將2根不銹鋼鋼管并排纏繞成內(nèi)外圈螺旋管, 將外緣與端蓋焊接在一起, 構(gòu)成環(huán)形空間作為鍋爐反應(yīng)器腔室, 內(nèi)部裝載金屬Li和啟動劑。鍋爐反應(yīng)器結(jié)構(gòu)如圖1所示。

啟動時, 啟動劑燃燒產(chǎn)生熱量加熱熔融固態(tài)金屬Li, 產(chǎn)生的液態(tài)Li與噴入的SF6發(fā)生化學(xué)反應(yīng), 以釋放的反應(yīng)熱加熱螺旋管內(nèi)的工質(zhì)水, 使工質(zhì)預(yù)熱、汽化、過熱, 生成高溫高壓蒸汽, 進入渦輪機推動渦輪做功。

圖1 鍋爐反應(yīng)器結(jié)構(gòu)

2 計算方法和流程

作為鍋爐反應(yīng)器設(shè)計的已知條件首要的是出口工質(zhì)壓力、出口工質(zhì)溫度和出口工質(zhì)流量, 其次是入口工質(zhì)溫度和入口工質(zhì)壓力。因此, 計算以出口工質(zhì)壓力、出口工質(zhì)溫度、出口工質(zhì)流量和入口工質(zhì)溫度為已知條件, 反推鍋爐反應(yīng)器結(jié)構(gòu)參數(shù)和入口工質(zhì)壓力。計算流程如圖2所示。

蒸汽生成過程分為過冷段、飽和段和過熱段三部分。過冷段入口為具有一定溫度的水, 出口是干度為0的飽和水; 飽和段入口是干度為0的飽和水, 出口是干度為1的飽和蒸汽; 過熱段入口是干度為1的飽和蒸汽, 出口為具有一定過熱度的蒸汽。根據(jù)鍋爐反應(yīng)器出口工質(zhì)壓力、出口工質(zhì)溫度和出口工質(zhì)流量, 假定過熱段壓降和過熱段管壁平均溫度, 計算過熱段平均溫度和平均壓力下的工質(zhì)普朗特數(shù)、運動粘度、導(dǎo)熱系數(shù)和對數(shù)溫差等物性參數(shù), 獲得壓力損失和管壁平均溫度計算值, 并與假設(shè)值比較, 如果誤差在允許范圍內(nèi), 計算過熱段管長、高度等結(jié)構(gòu)參數(shù)。然后繼續(xù)進行飽和段和過冷段的計算, 否則, 選擇假設(shè)值重新計算, 飽和段和過冷段的計算與過熱段的計算類似。

3 計算模型

假設(shè)鍋爐反應(yīng)器換熱管與Li熔池相鄰的一半用于換熱, 另一半不參與傳熱, 且Li熔池內(nèi)溫度均勻一致, 到換熱管的換熱系數(shù)為常值, 建立如下計算模型。

1) 計算段加熱功率

2) 計算段單位流通面積的工質(zhì)質(zhì)量流量

3) 計算段加速壓力降

4) 計算段摩擦阻力壓力降[6-7]

5) 計算段換熱面積

6) 計算段平均干度

7) 計算段傳熱系數(shù)[8-11]

計算段換熱管到工質(zhì)的換熱系數(shù)

計算段工質(zhì)平均溫度下的雷諾數(shù)

計算段曲度

飽和段修正系數(shù)

8) 計算段管壁平均溫度校對值

9) 計算段高度

4 計算結(jié)果和參數(shù)影響規(guī)律

4.1 入口工質(zhì)溫度影響分析

維持其他參數(shù)不變, 僅改變鍋爐反應(yīng)器入口工質(zhì)溫度, 過冷段參數(shù)變化曲線如圖3所示。

計算結(jié)果表明:

1) 鍋爐反應(yīng)器入口工質(zhì)溫度的增加相當(dāng)于工質(zhì)帶入了一部分熱量, 使得過冷段加熱功率減小, 換熱面積和長度減小, 由于過冷段長度變小, 過冷段工質(zhì)壓降變小, 入口工質(zhì)壓力降低;

圖3 過冷段參數(shù)隨入口工質(zhì)溫度變化曲線

2) 在維持出口工質(zhì)參數(shù)不變的情況下, 隨著鍋爐反應(yīng)器入口工質(zhì)溫度的增加, 過熱段和飽和段的換熱面積、長度和高度等參數(shù)和入口工質(zhì)參數(shù)不變;

3) 過冷段加熱功率減小使得鍋爐反應(yīng)器總加熱功率減小, 由過冷段組成的內(nèi)圈高度減小。

因此在構(gòu)建基于鍋爐反應(yīng)器的閉式循環(huán)動力時, 提高鍋爐反應(yīng)器入口工質(zhì)溫度, 不僅有利于提高系統(tǒng)效率, 還可減少鍋爐反應(yīng)器長度, 降低泵出口工質(zhì)壓力要求。

4.2 出口工質(zhì)溫度影響分析

維持鍋爐反應(yīng)器其他參數(shù)不變, 在穩(wěn)態(tài)值[–18%,18%]范圍內(nèi)改變出口工質(zhì)溫度, 計算結(jié)果如圖4～圖6所示。

計算結(jié)果表明, 隨著出口工質(zhì)溫度增加:

1) 過熱段加熱功率、換熱面積、長度和工質(zhì)壓降等參數(shù)大幅增加, 過熱段入口工質(zhì)壓力增加,為保證入口的干飽和蒸汽狀態(tài), 過熱段入口工質(zhì)溫度小幅增加;

圖4 過熱段參數(shù)隨出口工質(zhì)溫度變化曲線

圖5 飽和段參數(shù)隨出口工質(zhì)溫度變化曲線

圖6 過冷段參數(shù)隨出口工質(zhì)溫度變化曲線

2) 飽和段加熱功率、換熱面積、長度和工質(zhì)壓降等參數(shù)小幅減少, 入口工質(zhì)壓力增加, 為維持飽和段入口的飽和水狀態(tài), 入口工質(zhì)溫度小幅增加;

3) 過冷段加熱功率、換熱面積、長度和工質(zhì)壓降等參數(shù)小幅增加, 過冷段入口工質(zhì)壓力大幅增加。

綜合分析, 為滿足鍋爐反應(yīng)器出口工質(zhì)溫度的增加, 過熱段所需加熱功率大幅增加, 飽和段所需的加熱功率稍有降低, 過冷段所需加熱功率小幅增加, 飽和段降低的加熱功率幾乎與過冷段增加的加熱持平。由飽和段和過熱段組成的外圈高度大幅增加, 而由過冷段組成的內(nèi)圈高度稍有增加, 因此需適當(dāng)增加內(nèi)圈高度, 匹配外圈高度。

4.3 出口工質(zhì)壓力影響分析

維持鍋爐反應(yīng)器的其他參數(shù)不變, 在穩(wěn)態(tài)值[–15%, 15%]范圍內(nèi)改變出口工質(zhì)壓力, 計算結(jié)果如圖7～圖9所示。

圖7 過熱段參數(shù)隨出口工質(zhì)壓力變化曲線

圖8 飽和段參數(shù)隨出口工質(zhì)壓力變化曲線

圖9 過冷段參數(shù)隨出口工質(zhì)壓力變化曲線

結(jié)果表明, 為滿足鍋爐反應(yīng)器出口工質(zhì)壓力的增加:

1) 過熱段、飽和段和過冷段的入口工質(zhì)壓力相應(yīng)增加, 但在溫度不變的情況下蒸汽焓值有所降低, 導(dǎo)致總的加熱功率有所降低;

2) 由于過熱段、飽和段和過冷段的入口工質(zhì)壓力增加, 為達到相應(yīng)壓力下的工質(zhì)飽和, 工質(zhì)飽和溫度相應(yīng)增加, 因此鍋爐反應(yīng)器過熱段入口工質(zhì)溫度和飽和段入口工質(zhì)溫度增加;

3) 飽和段的換熱面積和長度等參數(shù)減小, 但過熱段的換熱面積和長度等參數(shù)稍有增加, 過冷段的換熱面積和長度等參數(shù)增加, 由過冷段組成的內(nèi)圈高度增加, 由飽和段和過熱段組成的外圈高度減小。

4.4 工質(zhì)流量影響分析

維持鍋爐反應(yīng)器其他參數(shù)不變, 在穩(wěn)態(tài)值[–20%, 20%]范圍內(nèi)改變鍋爐反應(yīng)器出口工質(zhì)流量, 計算結(jié)果如圖10～圖12所示。

圖10 過熱段參數(shù)隨出口工質(zhì)流量變化曲線

圖11 飽和段參數(shù)隨出口工質(zhì)流量變化曲線

圖12 過冷段參數(shù)隨出口工質(zhì)流量變化曲線

結(jié)果表明, 隨著鍋爐反應(yīng)器工質(zhì)流量的增加:

1) 鍋爐反應(yīng)器總的加熱功率增加, 由過冷段組成的內(nèi)圈高度增加, 由飽和段和過熱段組成的外圈高度增加;

2) 過熱段、飽和段和過冷段的加熱功率增加, 引起過熱段、飽和段和過冷段的換熱面積、長度等參數(shù)增加。

5 結(jié)束語

文中針對鍋爐反應(yīng)器螺旋換熱管進行的流動換熱過程研究, 探索了非對稱加熱條件下出入口參數(shù)變化對蒸汽生成過程的參數(shù)影響規(guī)律, 研究結(jié)果表明, 鍋爐反應(yīng)器入口工質(zhì)溫度的改變只影響過冷段參數(shù), 出口工質(zhì)溫度的改變對過熱段參數(shù)影響較大, 需要大幅改變過熱段加熱功率; 出口工質(zhì)壓力的改變對過熱段、飽和段和過冷段的參數(shù)影響都較大, 工質(zhì)流量的改變直接決定了加熱功率和加熱面積, 進而影響到壓力的變化。通過研究, 掌握了復(fù)雜物理化學(xué)過程中結(jié)構(gòu)參數(shù)和輸入輸出參數(shù)之間的耦合關(guān)系, 可為鍋爐反應(yīng)器的調(diào)節(jié)控制參數(shù)選取和穩(wěn)定調(diào)節(jié)提供支撐, 為提升系統(tǒng)穩(wěn)定域度提供參考。

[1] 查志武. 魚雷熱動力技術(shù)[M]. 北京: 國防工業(yè)出版社, 2006.

[2] 張小艷, 姜芳芳. 螺旋管換熱技術(shù)的研究現(xiàn)狀綜述[J]. 制冷與空調(diào), 2014, 28(1): 75-80.

Zhang Xiao-yan, Jiang Fang-fang. Present Situation of the Technical Research on Spiral Tube[J]. Refrigeration and Air Conditioning, 2014, 28(1): 75-80.

[3] 鄭邯勇. Li/SF6鍋爐反應(yīng)器內(nèi)流體的流動結(jié)構(gòu)與溫度分布[J]. 艦船科學(xué)技術(shù), 2001, 23(3): 44-47.

[4] 王樹峰. 鍋爐反應(yīng)器的參數(shù)選擇[J]. 艦船科學(xué)技術(shù), 2003, 25(z1): 40-42.

Wang Shu-feng. The Choice of the Boiler’s Parameter[J]. Ship Science and Technology, 2003, 25(z1): 40-42.

[5] 白杰, 黨建軍, 李代金, 等. 螺旋管蒸發(fā)器工質(zhì)流動換熱建模及仿真[J]. 魚雷技術(shù), 2015, 23(4): 301-304.

Bai Jie, Dang Jian-Jun, Li Dai-Jin, et al. Modeling and Simulation of Flow and Heat Transfer of Working Medium in Spiral Tube Evaporator[J]. Torpedo Technology, 2015, 23(4): 301-304.

[6] Jun Lu. One-dimensional Model of a Closed-Loop Underwater Propulsion System[J/OL]. AIAA Modelling and Simulation Technologies Conference. [2015-06-18]. https: //arc.aiaa.org/doi/abs/10.2514/6.2015-3362.

[7] Li D, Luo K, Dang J J. A Moving Boundary Model for Two-phase Flow Heat Exchanger Incorporated with Relative Velocities between Boundaries and Fluid[J]. Inter- national Journal of Heat and Mass Transfer, 2016, 95: 35-44.

[8] 鄔田華. 工程傳熱學(xué)[M]. 武漢: 華中科技大學(xué)出版社, 2011.

[9] 陶文權(quán). 傳熱學(xué)基礎(chǔ)[M]. 北京: 電力工業(yè)出版社, 1981.

[10] 景思睿. 流體力學(xué)[M]. 西安: 西安交通大學(xué), 2001.

[11] 袁媛, 彭敏俊, 夏庚磊, 等.螺旋管式直流蒸汽發(fā)生器熱工水力分析模型[J]. 原子能科學(xué)技術(shù), 2014, 48(z1): 251-256.

Yuan Yuan, Peng Min-jun, Xia Geng-lei, et al. Thermal-hydraulic Model for Helical-coil Once-through Steam Generator[J]. Atomic Energy Science and Technology, 2014, 48(z1): 251-256

[12] 韓勇軍, 楊赪石, 彭博, 等. 基于閉式循環(huán)動力系統(tǒng)的魚雷殼體冷凝器一維建模與仿真[J]. 魚雷技術(shù), 2015, 23(4): 291-295.

Han Yong-jun, Yang Cheng-shi, Peng Bo, et al. One-Di- mensional Modelling and Simulation of Torpedo Shell Condenser in Closed-loop Thermal Propulsion System[J]. Torpedo Technology, 2015, 23(4): 291-295.

Study on Parameter Effect of the Boiler Reactor Steam Generation Process under Asymmetric Heating

MA Wei-feng, LU Jun, WAN Rong-hua, HAN Yong-jun

(The 705 Research Institute, China State Shipbuilding Corporation Limited, Xi’an 710077, China)

To solve the small stability domain problems caused by the coupling of parameters and the time-varying parameters in the complex physicochemical process of a boiler reactor used in a torpedo Li/SF6closed-cycle power system,based on heat transfer theory, the heat transfer model of a helical exchange tube in a boiler reactor was established using segment processing. A numerical simulation of the heat transfer was conducted, and the influence of inlet-working fluid temperature, outlet-working fluid temperature, outlet-working fluid pressure, and outlet-working fluid flow on the boiler reactor steam generation process under asymmetric heating conditions was explored. The results can provide support for the system performance simulation, control optimization, and stability domain improvement, and provide a theoretical reference for regulation control for the design, optimization, and experimental verification of boiler reactors, and have certain reference values for the research and design of similar devices.

torpedo; closed-cycle power system; asymmetric heating; steam generation; boiler reactor

TJ630.32

A

2096-3920(2021)03-0326-07

10.11993/j.issn.2096-3920.2021.03.012

馬為峰, 路駿, 萬榮華, 等. 非對稱加熱條件下鍋爐反應(yīng)器蒸汽生成過程參數(shù)影響研究[J]. 水下無人系統(tǒng)學(xué)報, 2021, 29(3): 326-332.

2020-02-02;

2020-02-28.

馬為峰(1977-), 男, 碩士, 高工, 主要研究方向為水下能源動力技術(shù).

(責(zé)任編輯: 陳 曦)