2. 深海技术科学太湖实验室,江苏无锡 214082
2. Taihu Laboratory of Deepsea Technological Science, Wuxi 214082, China
船舶中央冷却系统仅中央冷却器使用海水冷却,具有腐蚀性的海水不直接接触和冷却各种散热设备,可有效提高各设备的可靠性和寿命,已成为各型船舶冷却系统设计的首要选择。作为全船冷却系统,在实际应用过程中,因淡水侧各支路流量分配不均导致设备不能得到有效冷却的现象经常出现[1]。因此,有必要对船舶中央冷却系统水力特征进行分析。
刘方等[2]以传统的管路流量计算为基础,提出多分枝管路流量分配迭代补偿计算方法。宋振国[3]以图论为抓手,建立了船舶管网水力计算的数学模型,并通过实船试验,验证了数学模型的正确性。贾志强等[4]将SIMPLE算法应用到水力管网进行动态模拟,实现了瞬态计算。陆蕾[5]分析了钻井平台冷却水系统达到流量平衡的关键技术,并使用标准矩阵求解技术等方法对冷却系统进行优化。多数研究在对中央冷却系统计算过程中进行简化,并且需要代入特定模型进行求解,为实际工程应用带来了不便。
本文对中央冷却系统进行介绍,从简单管路的特性出发,提出基于阻抗平衡法的管路特性及流量分配计算方法,对冷却水管路进行定量分析,并采用Flowmaster软件对计算方法进行验证,最后利用阻抗平衡法反向计算了设计流量状态下的各支管阻抗以及局部阻力系数,可在实际工程中进行应用。
1 船舶中央冷却系统介绍船舶中央冷却系统属于闭式冷却水系统的一种,利用海水冷却泵从海水门中抽取海水并送至中央冷却器,对中央冷却器的淡水回路进行冷却,海水升温后流出舷外,完成海水侧循环。经中央冷却器冷却后的低温淡水,经淡水冷却泵加压后沿各冷却支管流至主机、空调等设备的冷却器,带走舱内各设备产生的热量,随后经冷却总管汇集至中央冷却器,完成淡水侧循环。中央冷却水系统使用海水冷却的冷却器仅1套,材料要求低;需耐腐蚀的管路也较短;同时舱内管路不直接承压,提高了管路及其附件的安全性。船舶中央冷却系统可分为混流式中央冷却水系统、独立式中央冷却水系统、自流式中央冷却水系统[6]。图1为某PLG船中央冷却系统原理图[7]。
![]() |
图 1 某PLG船中央冷却系统原理图 Fig. 1 Schematic diagram of central cooling system of a PLG ship |
如图1所示,中央冷却系统海水侧管路形式单一,不存在流量分配不均问题。淡水侧由淡水冷却泵统一向各冷却用户输送冷却水,由于各冷却支管通径、长度以及所使用的管路附件不同,在阻抗平衡的原则下,管内实际流量与设计流量会有较大差别,必须进行支管阻力匹配计算和设计。中央冷却系统淡水侧为闭式循环,流体的位能与压力能、速度能互相转化,因此在计算过程中可不考虑各管路空间走向,视为平面系统进行计算。
2.1 物理模型以4台冷却设备为例,将中央冷却系统淡水侧计算模型简化,如图2所示。各冷却设备的流量分别为Q1~Q4,总流量为Q,分析可得:
![]() |
图 2 中央冷却系统淡水侧计算模型图 Fig. 2 Calculation model of fresh water side of central cooling system |
{Qb2b1=Qa1a2=Q1,Qb3b2=Qa2a3=Q1+Q2,Qb4b3=Qa3a4=Q1+Q2+Q3,Q=Q1+Q2+Q3+Q4。 | (1) |
冷却管路在船舱内多沿舷侧布置,
黏性流体管道计算的伯努利方程为:
z1+p1γ+u212g=z2+p2γ+u222g+Hw。 | (2) |
式中:
压头损失可分为沿程压力损失和局部阻力损失,对于简单管路计算如下:
Hw=λLd×u22g+Σζu22g=(λLd+Σζ)u22g。 | (3) |
式中:
将
Hw=8(λLd+∑ζ)π2d4gQ2。 | (4) |
定义管路阻抗[8]:
S=8(λLd+∑ζ)π2d4g。 | (5) |
在冷却系统实际应用中,由于管路附件的布置、管路腐蚀等原因,管路内流体基本处于完全湍流区,雷诺数对阻力影响很小,此时可视为局部阻力系数
λ=1(1.74+2lgd2Δ)2。 | (6) |
依据式(5)及式(6)可计算出各分支管路的阻抗,结合式(4),则有:
Hw=SQ2。 | (7) |
船舶冷却管路为简单管路经过串联和并联组成的管网。由不同管径或不同粗糙度的简单管路头尾相连组成串联管路,如图3所示,其各管段中流量相同,总压头损失等于串联各管段压头损失之和,总阻抗为串联各管段阻抗之和,即
![]() |
图 3 串联管路示意图 Fig. 3 Schematic diagram of tandem pipeline |
{Q=Q1+Q2+Q3,Hw=Hw1+Hw2+Hw3,S=S1+S2+S3。 | (8) |
由2个或2个以上简单管路入口与入口相连接、出口与出口相连接组成并联管路,如图4所示。可知,其总流量等于各分支管段流量之和,总压头损失等于各分支管段压头损失,总阻抗的平方根倒数等于各分支管段阻抗的平方根倒数之和,即
![]() |
图 4 并联管路示意图 Fig. 4 Schematic diagram of parallel pipeline |
{Q=Q1+Q2+Q3,Hw=Hw1=Hw2=Hw3,1√S=1√S1+1√S2+1√S3。 | (9) |
由式(9)可进一步推导出,并联各分支管段流量之比等于阻抗的平方根倒数之比,即
Q1:Q2:Q3=1√S1:1√S2:1√S3。 | (10) |
在冷却管路已设计完成的情况下,利用阻抗平衡法对淡水侧进行水力计算,参考图2,计算步骤如下:
步骤1 首先计算
步骤2 计算管段
Sb2b1a1a2=Sb2b1+Sb1a1+Sa1a2。 | (11) |
步骤3 将管段
1√S1,2=1√Sb2b1a1a2+1√Sb2a2。 | (12) |
经过上述步骤后,冷却系统计算模型可简化为图5所示。
![]() |
图 5 计算步骤分析图① Fig. 5 Analysis of calculation steps ① |
步骤4 参考图5,计算管段
Sb3b2a2a3=Sb3b2+S1,2+Sa2a3。 | (13) |
步骤5 将管段
1√S1,2,3=1√Sb3b2a2a3+1√Sb3a3。 | (14) |
经过上述步骤后,冷却系统计算模型可简化为图6所示。
![]() |
图 6 计算步骤分析图② Fig. 6 Analysis of calculation steps ② |
步骤6 参考图6,计算管段
Sb4b3a3a4=Sb4b3+S1,2,3+Sa3a4。 | (15) |
步骤7 再根据并联管路流量与阻抗分配关系,求解
{Q4:Q1,2,3=1√S4:1√Sb4b3a3a4,Q4+Q1,2,3=Q。 | (16) |
步骤8 参考图5,求解
{Q3:Q1,2=1√S3:1√S1,2,Q3+Q1,2=Q1,2,3。 | (17) |
步骤9 参考图2,求解
{Q2:Q1=1√S2:1√Sb2b1a1a2,Q2+Q1=Q1,2。 | (18) |
经过步骤1~步骤9,根据各管路阻抗,并结合简单管路串并联基本运算方法计算出各支路流量,利用式(7)可计算出各管段的压头损失。
2.3 算例分析以图2为算例进行演示,在进行淡水侧冷却管路水力及流量计算前,依据各设备供应商提供的散热量及系统冷却水温升数据,可计算出设备所需的冷却水流量,进一步根据管路流速要求,计算出各支管管径,各管段基本参数如表1所示。
![]() |
表 1 各管段基本参数 Tab.1 Basic parameters of each pipe section |
管路绝对粗糙度为0.2 mm,利用阻抗平衡法,计算出各管段流量、流速等数据如表2所示。
![]() |
表 2 计算出的各管段水力数据 Tab.2 Calculated hydraulic data of each pipe section |
计算结果表明,各管段计算出的实际流量与实际流速均与设计阶段有较大差距,可通过各调节支管阀门开度达到预期目的。
3 仿真验证为验证阻抗平衡法计算的正确性,采用Flowmaster软件对上述算例进行仿真。建立淡水侧冷却水系统的仿真模型,如图7所示。利用恒流源Flow元件模拟淡水冷却泵,以一恒定流量持续输出冷却水;为提高计算效率,建模过程中忽略各管段的弯头,各管段中设备、阀门及其他管路附件均简化为Discrete Loss元件,其阻力系数等于各管段中所有设备及附件阻力系数之和;选用Cylindrical Rigid元件模拟各管段,输入表1中各管段的规格参数,采用稳态计算模式,对淡水侧水力特性进行仿真。仿真数据如表3所示。
![]() |
图 7 Flowmaster仿真模型 Fig. 7 Flowmaster simulation model |
![]() |
表 3 Flowmaster仿真计算出的各管段数据 Tab.3 Data of each pipe section calculated by flowmaster |
分析可知,仿真流量与理论计算流量最大误差为1.23%,仿真流速与理论计算流速最大误差为1.31%,仿真得出的压头损失与理论计算压头损失最大误差为3.51%,验证了阻抗平衡法理论计算数据的正确性。
4 流量调节分析由计算结果可知,在各管路经过串并联组成管网的情况下,管段
管段
![]() |
图 8
增加管段 |
在流量调节过程中,为达到设计流量要求,可转变为已知各管段规格参数及流量,求解各管段局部阻力系数的问题,采用阻抗平衡法从反向进行计算,步骤如下:
S1+Sb2b1+Sa1a2S2=Q22Q21, | (19) |
1√S1,2=1√S2+1√S1+Sb2b1+Sa1a2, | (20) |
S1,2+Sb3b2+Sa2a3S3=Q23(Q1+Q2)2, | (21) |
1√S1,2,3=1√S3+1√S1,2+Sb3b2+Sa2a3, | (22) |
S1,2,3+Sb4b3+Sa3a4S4=Q24(Q1+Q2+Q3)2。 | (23) |
通过式(19)~式(23),无法直接求出各管段阻抗值,但可分析得出
利用上述方法,逐步求解出各管段阻抗以及局部阻力系数,如表4所示。
![]() |
表 4 各支管阻抗以及局部阻力系数 Tab.4 Impedance and local resistance coefficient of each branch |
分析管段
1)阻抗计算法可根据管路基本特征,准确计算出在初始状态下各支管的流量、流速、压头损失等参数,利用Flowmaster软件对算例进行了仿真,仿真流量与理论计算流量最大误差为1.23%,仿真流速与理论计算流速最大误差为1.31%,仿真计算得出的压头损失与理论计算压头损失最大误差为3.51,二者吻合良好,验证了阻抗平衡法理论计算数据的正确性。
2)管路流量可采用改变阀门开度进行调节,可利用阻抗平衡法反向计算出在设定流量下的各支管阻抗以及局部阻力系数,精准进行流量控制。
3)阻抗平衡法在计算过程中,充分考虑到每一管段的影响因素,并且无需求解复杂方程,可作为中央冷却系统水力的计算工具。
[1] |
李宝仁, 曹博, 仇艳凯, 等. 基于图论的船舶海水冷却管网水力特性[J]. 舰船科学技术, 2014, 36(2): 102-107. LI Bao-ren, CAO Bo, QIU Yan-kai, et al. Analysis on hydraulic characteristics of seawater cooling pipelines network of the ship based on graph theory[J]. Ship Science and Technology, 2014, 36(2): 102-107. |
[2] |
刘方, 吴鹏飞. 船舶淡水冷却系统流量分配的探析[J]. 船海工程, 2012, 41(5): 87-90,94. LIU F, WU P F. The flow distribution analysis of vessel's fresh water cooling system[J]. Ship & Ocean Engineering, 2012, 41(5): 87-90,94. DOI:10.3963/j.issn.1671-7953.2012.05.024 |
[3] |
宋振国. 船舶管网水力计算方法研究[D]. 大连: 大连海事大学, 2011.
|
[4] |
贾志强, 董葳. 船舶中央冷却系统水力管网动态模拟[J]. 船舶工程, 2012, 34(3): 43-46. JIA Z Q, DONG W. Dynamic simulation of water network in central cooling system of vessel[J]. Ship Engineering, 2012, 34(3): 43-46. DOI:10.3969/j.issn.1000-6982.2012.03.012 |
[5] |
陆蕾. 半潜式钻井平台冷却水系统设计及流体分析软件的应用[J]. 船舶工程, 2015, 37(5): 78-81. LU L. Design of semi-submersible drilling unit water cooling system and application of fluid analysis software[J]. Ship Engineering, 2015, 37(5): 78-81. |
[6] |
张博强. 基于变量调节的船舶冷却水系统节能优化研究[D]. 大连: 大连海事大学, 2017.
|
[7] |
中国造船工程学会. 船舶设计实用手册轮机分册[M]. 北京: 国防工业出版社, 2013.
|
[8] |
魏龙. 热工与流体力学基础[M]. 北京: 化学工业出版社, 2017.
|