液化石油气船舶液货管网的动态水力数值计算.docVIP

——————————————————————————————————————————————— 液化石油气船舶液货管网的动态水力数值计算 第28卷　第4期 2002年11月 大连海事大学学报 JournalofDalianMaritimeUniversity Vol.28,No.4Nov.,2002 文章编号:100627736(2002)0420029204 徐帮林,胡俊权,吴宛青 1 2 2 Ξ (1.中远(集团)总公司,北京;21大连海事大学) 摘要:,将吉布斯定律运用于两相流中,得 出两相平衡态下压力、,建立了稳定性和,,为对液化石油气船舶液货系统仿真操作培训提供了理论依据. 关键词:液化石油气船舶;液货系统;管网有限元;数值计算;仿真 中图分类号:U674.133;O242.21 　　文献标识码:A 0　引　言 根据《STCW78Π95公约》和《中华人民共和国散装液体货船船员特殊培训、考试和发证办法》,LPG船员岗前培训的一项重要内容是船舶液货装卸操作训练.对LPG??舶液货管网进行动态仿真是完成计算机数值仿真培训的基础.本文结合当前最具代表性的全冷式LPG船,针对其液货管网的工作特点,分别建立了液货单相流和气-液两相流水力模型,结合管网有限元理论,最终较好地解决了全冷式LPG船舶液货系统的动态数值仿真这一复杂问题. 　　以流量形式按Lagrange单元方式构造插值函数为 Qi= 25πdi λ8ρi(li+li当)Δpi 3 2 Δpi(2) 式中,Qi为管段i的体积流量,mΠs;Δpi为管段的压降,MPa;hwi为管段i的损失水头,m;λi为管段的沿程阻力系数;li为管段i的长度,m;li当为管段i的局部阻力损失当量长度,m;di为管段i的直径,m.1.2　两相流体 管网中一维气液两相流动的控制方程为两个相的连续方程、动量守恒方程、能量守恒方程以及气体状态方程.按Lockhart和Martinelli方法可以得出两相流单段管元的元素方程 4fi(li+li)mi Δpi= di2ρL 2 [2] 1　液货管网水力模型 1.1　单相流体 , 对于液货管网,单相液体流动符合质量守衡 和能量守衡定律.根据连续性方程、伯努利方程和达西公式等 [1] (3) ,可以得出单段管元的元素方程 Δpi=ρghwi = υ =ρg(λ+ξii)di2g 4Qi 2 2 　　以流量形式按Lagrange单元方式构造插值函数为 (1) Qi= 2fi(li+li当)ρLΔp2 i(li+li当)λπdi 2 Δpi(4) Ξ收稿日期:20022072121 基金项目:交通部重点科技计划资助项目(9506020). 作者简介:徐帮林(19742),男,河南光山人,研究生,主要从事特种液货船舶的研究. 大连海事大学学报 　第28卷　30式中,fi为与摩擦因子有关的变量;ρL为液相密度,kgΠm;mi为管段i中的流体质量流率,kgΠs. 不管单相流体还是两相流体,都可以表示成 (5)Qi=kiΔpi 式中,ki为待求未知插值函数,随压头水差而变化. 1.3　液货状态的确定 由工程热力学中著名的吉布斯相律I=C-P+2(I表示系统中独立强度参数的数目或称自由度;C表示系统包含的组元数;P表示系统包含 [3] 的相数)可知,对于两相平衡共存的单元系统,只需要知道一个强度参数,其余强度参数,即只要压力、 温度3 内能、比焓中的任何一个参数被确定,就能确定其 他的所有参数. 液货处于何种状态,确定其临界点是关键.在水力计算中,通过各管元结点理论压力的求解,可 t 获得该点平衡态下的理论温度{Ti},如果该点实际温度{Ti}低于{Ti},则说明液货处于单相即液态;反之,则说明液货处于两相状态. t 22.LPG船舶液货系统,按管,可对其进行管元和节点的系统拓( 见图1). 图1　典型LPG船舶液货管网系统拓扑图 2.2　总体方程式和总体矩阵方程式 根据前人的研究成果,针对两相流,本文采用有限管元法.为讨论问题的方便,以三管元组成子