基于查表插值算法的高压CO_(2)管道泄漏瞬态特性数值模拟  

Numerical study on the transient behaviour of high-pressure CO_(2) pipeline leakage based on lookup table interpolation algorithm

作　　者：鲁寨军[1,2,3] 王甲强 张梓轩 孙永龙 王根达 刘东润 LU Zhaijun;WANG Jia-qiang;ZHANG Zi-xuan;SUN Yong-long;WANG Gen-da;LIU Dong-run(Key Laboratory of Traffic Safety on Track,Central South University,Ministry of Education,Changsha 410075,China;Joint International Research Laboratory of Key Technology for Rail Traffic Safety,Changsha 410075,China;National&Local Joint Engineering Research Center of Safety Technology for Rail Vehicle,Changsha 410075,China)

机构地区：[1]中南大学轨道交通安全教育部重点实验室,长沙410075 [2]轨道交通安全关键技术国际合作联合实验室,长沙410075 [3]轨道交通列车安全保障技术国家地方联合工程研究中心,长沙410075

出　　处：《安全与环境学报》2023年第2期442-450,共9页Journal of Safety and Environment

摘　　要：为准确高效地模拟高压CO_(2)管道泄漏的瞬态特性,基于Fluent仿真平台,利用用户自定义真实气体模型(User Defined Real Gas Model,UDRGM)和用户自定义函数(User Defined Function,UDF),结合查表法和双线性插值法建立CO_(2)的真实气体模型,并将压力驱动的Lee模型通过用户自定义函数嵌入Fluent求解器来模拟CO_(2)的非平衡相变过程,建立了高压CO_(2)管道泄漏的非平衡相变数值模型。通过与Botros等的试验数据进行对比分析,验证了该模型的准确性。在此基础上,对比了上述模拟方法与编译S-W(Span-Wagner)状态方程模拟方法的精度和效率,最后使用本模型研究了不同初始压力对高压CO_(2)管道泄漏瞬态特性的影响。结果表明:两种模拟方法精度接近,最大相差为7.37%,但提出的模拟方法效率明显优于编译S-W状态方程的模拟方法,计算时间相较缩短约86.9%;初始压力为11.27 MPa的最大总出口质量流量比4.36 MPa的大7.24 kg/s,而对应的气相CO_(2)出口质量流量却低0.14 kg/s;初始压力为11.27 MPa和4.36 MPa对应的近场最大射流流速分别为155.12 m/s和175.50 m/s;管内初始压力越低,泄漏后到达亚稳态时过热度越大,管内及泄漏口附近相变程度越剧烈,气相出口质量流量越大,近场射流峰值速度越大。这项研究可为工业规模管道泄漏的模拟和管道泄漏的三维模拟提供更高效的方式。This paper established an efficient and accurate non-equilibrium phase change numerical model for calculating high-pressure CO2pipeline leakage based on ANSYS Fluent. In this model, the look-up table interpolation method for pure CO2was established using the User Defined Real Gas Model(UDRGM) and User Defined Functions(UDFs). Additionally, the pressure-driven Lee model was embedded in the Fluent solver through UDF to simulate the non-equilibrium phase change process of CO2. The accuracy of the model was verified by comparative analysis with the experimental data of Botros et al. Furthermore, the accuracy and efficiency of the simulation method used in this paper and the simulation method of compiling the S-W(Span-Wagner) equation of state were compared. Finally, the model was used to study the effect of different initial pressures on the transient characteristics of high-pressure CO2pipeline leakage. The calculation results show that the accuracy of the two simulation methods is close, and the maximum difference is 7.37%. However, the efficiency of the simulation method in this paper is significantly better than that of the S-W equation of state, and the calculation time is reduced by about 86.9%. When the initial pressure is 11.27 MPa and 4.36 MPa, the pressure drop rate Δp/Δt is 22 975 MPa/s and 5 975 MPa/s, the temperature drop rate ΔT/Δt is 17 325.68 K/s and 5 200 K/s, and the initial decompression wave speed is 592 m/s and 450 m/s. The maximum total outlet mass flow with an initial pressure of 11.27 MPa is 7.24 kg/s larger than that of 4.36 MPa, while the corresponding gas-phase CO2outlet mass flow is 0.14 kg/s lower. When the initial pressure is 11.27 MPa and 4.36 MPa, the maximum flow velocity in the near field is 155.12 m/s and 175.50 m/s. The results show that the lower the initial pressure in the pipe, the more severe the phase change in the pipe and near the leakage opening, the greater the gas outlet mass flow rate, and the greater the peak velocity of the near-field jet. This research can provide

关 键 词：安全工程 环境工程 二氧化碳 非平衡相变 计算流体力学 管道泄漏 气液两相流 数值模拟 

分 类 号：X937[环境科学与工程—安全科学]