| Citation: | DUAN Chaowei, SONG Pu, HU Hongwei, FENG Haiyun. Simple Method for the Calculation of Bubble Pulsation Period in Underwater Explosion[J]. Chinese Journal of High Pressure Physics, 2022, 36(1): 015101. doi: 10.11858/gywlxb.20210782
|
近年来随着液化天然气(LNG，Liquefied Natural Gas)产业在全球迅速发展，天然气的液化技术和设备也在不断发展完善、日渐成熟。我国LNG领域内的相关研究起步较晚，许多技术远远落后发达国家水平，在天然气液化工艺及装置的生产等方面缺乏自主产权。因此，开展天然气液化工艺及装置的研究，对于实现液化装置的国产化、高效化有十分重要的意义[1-3]。
超声速旋流分离技术是一种新兴的天然气加工处理技术，被较为广泛地用于天然气脱水、脱重烃、脱酸等方面，近年来开始逐渐应用于天然气液化方面[4-7]。天然气超声速液化的原理是：高压天然气混合物在Laval喷管内达到一定的温度、压力条件，开始凝结成核，最终凝结成液滴，在后续工艺中进行进一步气液分离。与传统的天然气液化技术相比，具有结构工艺简单、支持无人操作(适用于海底天然气处理)、对水合物抑制剂依赖性小、投资和运行成本低等优势[8-9]。
为了探究入口复杂多变的压力条件对天然气超声速液化特性的影响，对甲烷-乙烷气体混合物的超声速凝结流动特性进行研究，在凝结成核与生长理论的基础上建立了适用于甲烷-乙烷双可凝气体混合物的凝结流动数学模型，重点研究了入口压力对天然气混合物在Laval喷管内主要流动与凝结参数的影响规律。
Laval喷管结构主要包括入口段、渐缩段、喉部及扩张段4部分[10-11]。各部分参数如表 1所示，L0为入口长度，r1为渐缩段入口半径，L1为渐缩段长度，rcr为喉部截面半径，L2为渐扩段长度，r2为渐扩段出口半径。为尽量减小流场涡流的影响，渐缩段采用双三次曲线设计，喉部采用一段平缓光滑的圆弧作为过渡曲线，渐扩段采用等膨胀率设计，膨胀率取为10 000 s-1。考虑到实验加工方便，保证曲面的精度，且能够更加直观地观察Laval喷管内部的流场分布情况，所设计Laval喷管截面采用矩形截面，三维结构如图 1所示。
| L0/mm | r1/mm | L1/mm | rcr/mm | L2/mm | r2/mm |
| 50.00 | 20.00 | 56.01 | 2.50 | 71.28 | 6.15 |
采用欧拉双流体模型开展数值计算，控制方程主要包括气相流动方程和液相流动方程。在无滑移假设及欧拉双流体模型的前提下分别建立气相及液相流动控制方程组，液滴数目守恒方程及液滴半径、数目、湿度关系式分别添加到对应源相方程中[12-14]。
气相流动控制方程组
|
∂ρv∂t+∂∂xj(ρvuj)=Sm |
(1) |
|
∂ρc2∂t+∂∂xj(ρc2uj)=Sm,c2 |
(2) |
|
∂∂t(ρvui)+∂∂xj(ρvujui)=−∂pv∂xi+∂∂xj[μ(∂uj∂xi+∂ui∂xj−23δij∂uj∂xj)]+∂∂xj(−ρv¯u′iu′j)+Su |
(3) |
|
∂∂t(ρvE)+∂∂xj(ρvujE+ujpv)=∂∂xj(keff∂T∂xj+uiτeff)+Sh |
(4) |
液相流动控制方程组
|
∂∂t(ρY)+∂∂xj(ρujY)=SY |
(5) |
|
∂ρN∂t+∂∂xj(ρNuj)=J |
(6) |
|
rd=3√3Y/(4πρLN) |
(7) |
式中:ui、uj为时均速度分量，m/s；ρv为气相密度，kg/m3；ρ为气液混合相密度，kg/m3；pv为时均压力，Pa；μ为黏度，kg/(m·s)；δij为Kronecker delta数；E为总能量，J/kg；keff为有效导热系数，W/(m·K)；τeff为有效应力张量，无量纲；Y为液相质量分数，无量纲；rd为液滴半径，m；drd/dt为液滴生长速率，m/s；N为液滴数目，kg-1。
成核模型采用文献[15-16]中提出的双组分气体自发凝结成核模型修正方法。液滴生长过程采用Gyarmathy液滴生长模型，模型中液滴与气体间的传热系数[17-18]为
|
kr=λvrd11+2√8π1.5Pvγ1+γKn |
(8) |
依据传热、传质过程，可推导得到液滴生长速率计算模型
|
drddt=λvρLhLV(Ts−T)(1−rcrd)rd(1+2√8π1.5Pvγ1+γKn) |
(9) |
式中:λv为气体导热系数，W/(m·K)；Pv为气体Prandtl数；γ为气体比热比；hLV为凝结潜热，J/kg；Ts为气体压力对应的饱和温度，K；Kn表示Kundsen数。由于双组分气体不存在压力对应的饱和温度Ts这一概念，将双组分相图中露点线类比于单组分中饱和曲线。
针对气体状态方程的选择，由于低温气体已偏离理想气体假设，本研究采用了NIST真实气体模型进行计算。
湍流发生时会导致流体之间相互交换动量、能量，也会造成浓度的改变。本研究建模时忽略相间速度的滑移，即液滴产生不影响湍流，因此只考虑气相的湍流方程。FLUENT中提供了以下几种湍流模型：S-A模型、标准k-ε模型、RNG k-ε模型、Realizable k-ε模型、k-ω模型以及雷诺应力模型。S-A模型主要应用流动分离区附近模拟，标准k-ε模型、RNG k-ε模型一般用于各向同性的均匀湍流，k-ω模型可用于带压梯度的流动模拟和跨声速激波模拟，雷诺应力模型主要用于龙卷风、燃烧室等强烈旋转流动的模拟。对Laval喷管内跨声速流动，采用k-ω模型可以获得较为理想的计算精度和计算速度，故本研究采用该模型进行数值计算。
气体在Laval喷管中的流动属于高速可压缩流动，采用密度基进行求解，流动控制方程组、湍流动能方程、湍流耗散率方程均采用二阶迎风格式进行离散。
根据双组分气体在Laval喷管内的高速可压缩的流动特性，入口和出口边界设置为压力入口边界和压力出口边界条件，对于气体在Laval喷管内的超声速流动，由于所有的流动参数都可从Laval喷管内部外推得到，故在出口处不进行相应设置，固体壁面边界设置为无滑移、无渗流、绝热边界条件。
在数学模型中，由于气相方程添加了源相方程，液相方程定义了标量以及引入的真实气体方程，这些仅靠在FLUENT自带的模型和材料物性无法满足要求，需要编写相应的用户自定义函数(UDF)。本研究编写UDF时，分别定义DEFINE AJUST、DEFINE SOURCE和DEFINE PROPERTY 3个宏函数。DEFINE AJUST宏用来定义过饱和度、过冷度、成核速率、液滴生长率、液滴半径、液滴质量以及液滴表面张力等参数，DEFINE SOURCE宏用来定义控制方程中的质量、动量和能量源相，DEFINE PROPERTY用来定义数值计算中用到的真实气体的热力学参数如黏度系数、导热系数等。
为验证所建立的双组分气体凝结数学模型及数值计算方法的准确性，采用本研究所设计的Laval喷管结构，在中国石油大学(华东)超声速气体凝结流动实验系统开展了水-乙醇双可凝组分气体凝结相变实验研究。实验条件为：Laval喷管入口压力0.586 MPa，入口温度288.05 K，气体湿度98.1%，水与乙醇摩尔体积比84:16，气体流量为323.78标方每小时，实验测得的Laval喷管沿程压力分布如图 2所示，可以看出，压力分布实验结果与数值计算结果吻合较好，说明本研究所建立的双组分气体超声速凝结流动特性数学模型及数值计算方法具有一定的准确性和可靠性。
保持入口温度及组成(甲烷体积分数90%、乙烷体积分数10%)不变，研究不同的入口压力对Laval喷管内部甲烷-乙烷双组分气体凝结过程中压力、温度、成核率、液滴生长率、液滴半径、液相质量分数的影响。在数值计算中设定的入口温度为270 K，设定入口压力分别为5.5、6.0和6.5 MPa。Laval喷管内双组分气体凝结参数的变化趋势及对比如图 3~图 8所示。
从压力与温度分布可以看出，气体进入Laval喷管后压力、温度不断降低，当达到一定过冷度时，气体发生凝结并释放潜热，但凝结突跃现象对压力造成的影响并不显著，压力在Laval喷管渐扩段减小到了一个比较稳定状态，温度在减小到最小值后又略微上升，这主要是由于液滴凝结释放潜热引起的。随着入口压力的增大，出口压力略微升高，温升位置有所提前，出口温度也越高，这是因为，随着入口压力的增大，液滴成核与凝结量也随之增多，释放的潜热也就越多。
从成核率分布可以看出，保持其他条件一致，当压力发生变化时，成核速率的变化趋势几乎相同，在刚进入Laval喷管的一段距离内为零，在某一位置处开始，成核率从零开始突跃一直增大到峰值后迅速减小至零。随着入口压力从5.5 MPa增大到6.5 MPa，成核的发生位置(Wilson点)不断向前移动，逐渐向Laval喷管喉部靠拢，且成核率的最值逐渐增大。当压力为5.5 MPa时，成核发生位置为x=0.120 6 m，且在x=0.147 5 m处达到极限成核，为4.044×1020 m-3·s-1；当压力为6 MPa时，成核发生位置较5.5 MPa时向前移动，为x=0.119 9 m，极限成核位置也随压力的增大而前移，在x=0.139 3 m处达到8.062×1020 m-3·s-1；当压力继续增大为6.5 MPa时，成核发生位置较6 MPa时更加靠近喉部，为x=0.118 2 m，极限成核位置在x=0.132 8 m处，为9.015×1020 m-3·s-1。
从液滴半径分布可以看出，随着入口压力的增大，Laval喷管内平均液滴半径越大，出口液滴半径也随之增大，当压力为5.5、6.0、6.5 MPa时，对应最大液滴半径尺寸分别为415.86、447.88和477.44 nm。由此可知，压力的升高有利于液滴的生长。
从液滴生长率分布可以看出，液滴生长率在气体刚进入Laval喷管时一直为零，当液滴开始发生成核凝结时液滴生长率开始突增，变化到最大值后又迅速减小，最终减小为零。综合图 6和图 7还可以看出，随着入口压力的升高，在成核开始时液滴生长率较大，液滴半径增长速度较快，但一段距离后液滴生长率下降更大，液滴半径增长速度也明显放缓。
从液相质量分数分布可以看出，伴随着混合气体的凝结成核，液相质量分数也不断增大，且随着入口压力的升高，Laval喷管出口处的湿度值随之增大，当压力为5.5 MPa时，湿度的最大值为3.989 2%，当压力增大到6.5 MPa时湿度最终增大到7.382 0%。
(1) 建立了三维双组分天然气混合物超声速凝结流动数学模型，对Laval喷管内双组分混合物凝结流动进行了数值模拟，得出沿Laval喷管轴向的参数分布，通过开展双可凝组分气体凝结相变实验，对比发现数值模拟与实验结果基本一致，说明了所建立数学模型及计算方法的正确性。
(2) 利用数值模型研究了入口参数对天然气混合物超声速液化特性的影响，结果表明，保持Laval喷管入口温度及组成不变，增大入口压力，混合气体成核位置前移，成核率、平均液滴半径、液相质量分数均随之增大，即入口压力越大，混合气体在Laval喷管内越易发生凝结。在实际生产中，可以通过调节入口压力来促进天然气的凝结，提高Laval喷管的液化效率。
| [1] |
COLE R H. Underwater explosion [M]. Princeton, NJ: Princeton University Press, 1948.
|
| [2] |
姚熊亮, 刘文韬, 张阿漫, 等. 水下爆炸气泡及其对结构毁伤研究综述 [J]. 中国舰船研究, 2016, 11(1): 36–45. doi: 10.3969/j.issn.1673-3185.2016.01.006
YAO X L, LIU W T, ZHANG A M, et al. Review of the research on underwater explosion bubbles and the corresponding structural damage [J]. Chinese Journal of Ship Research, 2016, 11(1): 36–45. doi: 10.3969/j.issn.1673-3185.2016.01.006
|
| [3] |
姚熊亮, 汪玉, 张阿漫. 水下爆炸气泡动力学[M]. 哈尔滨: 哈尔滨工程大学出版社, 2012: 22−76.
YAO X L, WANG Y, ZHANG A M. Underwater explosion bubble dynamics [M]. Harbin: Harbin Engineering University Press, 2012: 22−76.
|
| [4] |
张阿漫, 汪玉, 闻雪友, 等. 水下爆炸气泡动态特性研究综述 [J]. 船舶力学, 2009, 13(5): 828–836. doi: 10.3969/j.issn.1007-7294.2009.05.022
ZHANG A M, WANG Y, WEN X Y, et al. Review of the dynamics of the underwater explosion bubble [J]. Journal of Ship Mechanics, 2009, 13(5): 828–836. doi: 10.3969/j.issn.1007-7294.2009.05.022
|
| [5] |
黄毅, 张弩. 水下爆炸气泡脉动载荷影响因素分析 [J]. 中国舰船研究, 2013, 8(6): 33–39.
HUANG Y, ZHANG N. Analysis of the impact parameters of the bubble pulsation load in underwater explosion [J]. Chinese Journal of Ship Research, 2013, 8(6): 33–39.
|
| [6] |
张姝红, 周华, 权琳, 等. 水下爆炸气泡脉动试验研究 [J]. 爆破, 2014, 31(1): 106–109. doi: 10.3963/j.issn.1001-487X.2014.01.025
ZHANG S H, ZHOU H, QUAN L, et al. Experimental study for bubble impulsion in underwater explosion [J]. Blasting, 2014, 31(1): 106–109. doi: 10.3963/j.issn.1001-487X.2014.01.025
|
| [7] |
李玉民, 倪芝芳. 水中爆炸气泡脉动流场的数值计算 [J]. 爆炸与冲击, 1996, 16(4): 377–381.
LI Y M, NI Z F. Numerical calculation on flow field of the bubble pulse by explosion underwater [J]. Explosion and Shock Waves, 1996, 16(4): 377–381.
|
| [8] |
张志江, 徐更光, 王廷增. 炸药水中爆炸气泡脉动分析计算 [J]. 爆破, 2007(1): 17–20. doi: 10.3963/j.issn.1001-487X.2007.01.005
ZHANG Z J, XU G G, WANG T Z. Analysis and calculation on bubble-oscillation by underwater explosion [J]. Blasting, 2007(1): 17–20. doi: 10.3963/j.issn.1001-487X.2007.01.005
|
| [9] |
王振宇, 唐兆田, 王善, 等. 水下爆炸气泡脉动的数值计算 [J]. 哈尔滨工程大学学报, 2008, 29(1): 22–27. doi: 10.3969/j.issn.1006-7043.2008.01.005
WANG Z Y, TANG Z T, WANG S, et al. Numerical calculations of the bubble pulsation forces of underwater explosions [J]. Joumal of Harbin Engineering University, 2008, 29(1): 22–27. doi: 10.3969/j.issn.1006-7043.2008.01.005
|
| [10] |
马海洋, 顾文彬, 夏卫国, 等. 气泡脉动法计算炸药能量的探讨 [J]. 西部探矿工程, 2003, 15(9): 117–118, 128. doi: 10.3969/j.issn.1004-5716.2003.09.057
MA H Y, GU W B, XIA W G, et al. Discussion of method of calculating explosive energy using the bubble moving principle [J]. West-China Exploration Engineering, 2003, 15(9): 117–118, 128. doi: 10.3969/j.issn.1004-5716.2003.09.057
|
| [11] |
赵铮, 陶钢, 杜长星. 爆轰产物JWL状态方程应用研究 [J]. 高压物理学报, 2009, 23(4): 277–282. doi: 10.3969/j.issn.1000-5773.2009.04.007
ZHAO Z, TAO G, DU C X. Application research on JWL equation of state of detonation products [J]. Chinese Journal of High Pressure Physics, 2009, 23(4): 277–282. doi: 10.3969/j.issn.1000-5773.2009.04.007
|
| [12] |
徐维铮, 吴卫国. JWL状态方程及其等效多方状态方程在内爆炸计算中的应用分析 [J]. 中国舰船研究, 2019, 14(3): 83–91.
XU W Z, WU W G. Application analysis of JWL EOS and the equivalent polytropic EOS in internal explosion calculation [J]. Chinese Journal of Ship Research, 2019, 14(3): 83–91.
|
| [13] |
孙美玲. 常微分方程数值解法的Matlab计算与可视比较 [J]. 高教学刊, 2016(19): 60–61.
SUN M L. Matlab calculation and visual comparison of ordinary differential equation numerical solution method [J]. Journal of Higher Education, 2016(19): 60–61.
|
| [14] |
宋浦, 杨凯, 梁安定, 等. 国内外TNT炸药的JWL状态方程及其能量释放差异分析 [J]. 火炸药学报, 2013, 36(2): 42–45. doi: 10.3969/j.issn.1007-7812.2013.02.010
SONG P, YANG K, LIANG A D, et al. Difference analysis on JWL-EOS and energy release of different TNT charge [J]. Chinese Journal of Explosives & Propellants, 2013, 36(2): 42–45. doi: 10.3969/j.issn.1007-7812.2013.02.010
|
| [15] |
王建灵, 赵东奎, 郭炜, 等. 水下爆炸能量测试中炸药入水深度的确定 [J]. 火炸药学报, 2002, 25(2): 30–31, 44. doi: 10.3969/j.issn.1007-7812.2002.02.012
WANG J L, ZHAO D K, GUO W, et al. Determination of the reasonable depth of explosives in water to measure underwater explosive energy [J]. Chinese Journal of Explosives & Propellants, 2002, 25(2): 30–31, 44. doi: 10.3969/j.issn.1007-7812.2002.02.012
|
| [16] |
郅斌伟, 马硕, 张志江. 水底爆炸气泡脉动特性 [J]. 爆破, 2009, 26(3): 22–24, 72. doi: 10.3963/j.issn.1001-487X.2009.03.006
ZHI B W, MA S, ZHANG Z J. Characteristic of bubble pulsation by underwater explosion [J]. Blasting, 2009, 26(3): 22–24, 72. doi: 10.3963/j.issn.1001-487X.2009.03.006
|
| [17] |
周霖, 谢中元, 陈勇. 炸药水下爆炸气泡脉动周期工程计算方法 [J]. 兵工学报, 2009, 30(9): 1202–1205. doi: 10.3321/j.issn:1000-1093.2009.09.010
ZHOU L, XIE Z Y, CHEN Y. An engineering calculation method on the bubble pulse period of underwater explosion [J]. Acta Armamentarii, 2009, 30(9): 1202–1205. doi: 10.3321/j.issn:1000-1093.2009.09.010
|
| [18] |
郅斌伟, 张志江, 马硕, 等. 水底爆炸冲击波峰值压力数值仿真 [J]. 爆破, 2009, 26(1): 22–24, 28. doi: 10.3963/j.issn.1001-487X.2009.01.006
ZHI B W, ZHANG Z J, MA S, et al. Research on shock wave peak pressure of under water explosion [J]. Blasting, 2009, 26(1): 22–24, 28. doi: 10.3963/j.issn.1001-487X.2009.01.006
|
| [19] |
张远平, 池家春, 龚晏青, 等. 含铝炸药水下爆炸性能的实验研究 [J]. 高压物理学报, 2010, 24(4): 316–320. doi: 10.11858/gywlxb.2010.04.013
ZHANG Y P, CHI J C, GONG Y Q, et al. Experimental study on underwater explosion performance of aluminized explosive [J]. Chinese Journal of High Pressure Physics, 2010, 24(4): 316–320. doi: 10.11858/gywlxb.2010.04.013
|
| [20] |
安丰江, 吴成, 王宁飞. 水下爆炸能量耗散特性分析研究 [J]. 北京理工大学学报, 2011, 31(4): 379–382.
AN F J, WU C, WANG N F. A research on the energy dissipation of underwater explosion [J]. Transactions of Beijing Institute of Technology, 2011, 31(4): 379–382.
|
Figures(3) / Tables(4)
Supported by: Beijing Renhe Information Technology Co. Ltd support:
info@rhhz.net
BIAN Jiang, CAO Xuewen, YANG Wen, YU Hongxi, YIN Pengbo. Influence of Inlet Pressure on Supersonic Liquefaction of Natural Gas Mixtures[J]. Chinese Journal of High Pressure Physics, 2018, 32(3): 031101. doi: 10.11858/gywlxb.20170639
| L0/mm | r1/mm | L1/mm | rcr/mm | L2/mm | r2/mm |
| 50.00 | 20.00 | 56.01 | 2.50 | 71.28 | 6.15 |
下载:
DownLoad:
