基于OpenFOAM的高温熔融锡液水力碎化数值模拟研究.pdf

International Journal of Fluid Dynamics 流体动力学流体动力学,2023,11(2),45-56 Published Online June 2023 in Hans.https:/www.hanspub.org/journal/ijfd https:/doi.org/10.12677/ijfd.2023.112005 文章引用文章引用:纪国剑,单梦琪,孟兴成,王政伟.基于 OpenFOAM 的高温熔融锡液水力碎化数值模拟研究J.流体动力学,2023,11(2):45-56.DOI:10.12677/ijfd.2023.112005 基于基于OpenFOAM的高温熔融锡液水力碎化数值的高温熔融锡液水力碎化数值模拟研究模拟研究 纪国剑纪国剑1,2*，单梦琪单梦琪1，孟兴成，孟兴成1，王政伟，王政伟1 1常州大学能源学院，江苏 常州 2江苏省绿色过程装备重点实验室，江苏 常州 收稿日期：2023年3月21日；录用日期：2023年3月28日；发布日期：2023年6月12日 摘摘 要要 为揭示高温熔融锡液与水相互作用过程中的碎化机理，基于开源软件包为揭示高温熔融锡液与水相互作用过程中的碎化机理，基于开源软件包OpenFOAM-v1806下的求解器下的求解器icoReactingMultiphaseInterFoam模拟熔融锡液入水后的水力碎化过程，分析瑞利模拟熔融锡液入水后的水力碎化过程，分析瑞利泰勒不稳定性、开泰勒不稳定性、开尔文尔文亥姆霍兹不稳定性及界面剪切力对熔融锡液柱碎化过程作用过程，探究水力碎化对液柱直径及液亥姆霍兹不稳定性及界面剪切力对熔融锡液柱碎化过程作用过程，探究水力碎化对液柱直径及液柱入水初速度对碎化剧烈程度的影响。研究表明液柱碎化剧烈程度、液柱单位时间内入水深度与液柱直柱入水初速度对碎化剧烈程度的影响。研究表明液柱碎化剧烈程度、液柱单位时间内入水深度与液柱直径成反比，液柱入水后碎化剧烈程度、锡凝固速度与液柱入水时初速度成正比。径成反比，液柱入水后碎化剧烈程度、锡凝固速度与液柱入水时初速度成正比。关键词关键词 高温熔融金属，水力碎化，计算流体力学，高温熔融金属，水力碎化，计算流体力学，OpenFoam模拟模拟 Numerical Simulation of Hydraulic Fragmentation of High-Temperature Molten Tin Based on OpenFOAM Guojian Ji1,2*,Mengqi Shan1,Xingcheng Meng1,Zhengwei Wang1 1School of Energy,Changzhou University,Changzhou Jiangsu 2Jiangsu Key Laboratory of Green Process Equipment,Changzhou University,Changzhou Jiangsu Received:Mar.21st,2023;accepted:Mar.28th,2023;published:Jun.12th,2023 Abstract In order to reveal the mechanism of fragmentation during the interaction between high-temperature *第一作者。纪国剑 等 DOI:10.12677/ijfd.2023.112005 46 流体动力学 molten tin and water,the icoReactingMultiphaseInterFoam solver within OpenFOAM-v1806 was used to simulate the hydraulic fragmentation process of molten tin entering water.The effects of Rayleigh-Taylor instability,Kelvin-Helmholtz instability and interfacial shear force on the frag-mentation process of molten tin were analyzed.And the sensitivity of the degree of fragmentation to the diameter of the liquid column and the initial velocity of the liquid column entering the wa-ter was explored.Research shows that the degree of fragmentation of the liquid column and un-derwater penetration in unit time are inversely proportional to the diameter of the liquid column.Meanwhile,the degree of fragmentation and the formation rate of solid tin are proportional to the initial speed of the liquid column entering water.Keywords High-Temperature Molten Metal,Hydraulic Fragmentation,Computational Fluid Dynamics,Numerical Simulation Based on OpenFOAM Copyright 2023 by author(s)and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License(CC BY 4.0).http:/creativecommons.org/licenses/by/4.0/1.引言引言 在工业生产中水作为最常见的冷却剂，与高温熔融金属接触反应过程中常伴随剧烈的传热传质过程、熔融金属碎化过程以及多种形式的能量转化过程。在高温熔融金属遇水碎化过程数值模拟研究中，周源 1认为水力学碎化过程发生在压力峰值的两端。袁明豪 2基于熔融液滴水力碎化模拟研究发现液滴变形与液滴表面的水力不稳定性是液滴碎化的主要机理。杨志 3利用 MPS 数值模拟程序将熔滴碎化区域分为凝固变形区、破碎快速增长区、破碎缓慢增长区，碎化程度随温度及 Weber 数的增长呈先增后减的趋势。钟明君 4根据 METRIC 程序对熔融物水力学碎化行为及其界面特性进行数值计算，研究表明在高 Weber数下液柱碎化主要受侧面不稳定波和剥离效应主导。沈正祥等 5采用有限体积法与水平集(Level-set)方法捕捉低速射流碎化分解过程中到界面的不稳定发展、变形、流场内部速度场及压力场分布等；研究获悉当入射速度较低时，流场内轴向速度随入射深度增加而变大，射流表面碎化范围内径向速度变化比较剧烈，对应的压力场为负压，射流的碎化长度与雷诺数(Re/We0.5)呈指数型关系。为揭示熔融金属与水相互作用机理，前人多利用计算流体动力学相关软件，如 ANSYS CFD、ANSYS CFX、MPS、MC3D 等来进行数值模拟研究。其中，OpenFOAM 软件运用较少，但其作为一种开源软件，可编译多种符合实际应用的求解器，用于计算分析不可压缩或可压缩流体的相变换热过程、燃烧过程、浮力驱动过程、固体应力分析等。本研究选取金属锡为研究材料，基于开源软件包 OpenFOAM-v1806下的 icoReactingMultiphaseInterFoam 求解器模拟熔融锡液柱入水后的水力碎化过程，分析瑞利泰勒(Rayleigh-Taylor,R-T)不稳定性、开尔文亥姆霍兹(Kelvin-Helmholtz,K-H)不稳定性及界面剪切力对熔融锡液柱碎化过程作用结果，探究水力碎化对液柱直径及液柱入水初速度对碎化剧烈程度的影响。2.数值模拟数值模拟方法方法 2.1.VOF 模型模型 VOF 模型对计算域内不互溶的多种流体求解同一动量方程，网格单元内各流体的体积分数由和为 1的条件来约束。i为 0 或 1 时，分别表示网格内无该流体或均为该流体，介于之间则表示存在其他流Open AccessOpen Access纪国剑 等 DOI:10.12677/ijfd.2023.112005 47 流体动力学 体。11jii=(1)因多相流系统内存在相变，以液相向汽相运输为例，相体积运输方程如下：()111liqliqliqliqliqvapmt+=u (2)密度、粘度、比热容等相关物理量与网格内介质的体积分数相关，物性表达公式如下：1jiii=(3)1jiii=(4)1jiii=(5)1jpipiiCC=(6)式中，为体积分数；t 为时间，s；u 为流体矢量流速，m/s；为密度，kg/m3；为动力粘度，kg/(ms)；为导热系数，W/(mK)；Cp为定压比热容，J/(kgK)；j 为计算域内的相数；i 为计算域内每一相对应的标号；下标 liq、vap 分别指液相、汽相。在 setFieldsDict 中定义初始相域，所有单元均需被定义且同一域中组分不可叠加，要求满足下式：11jiiY=(7)式中，Y 为组分分数。连续性方程是质量守恒定律在流体力学中的具体应用，使用前提是将流体看作连续介质，应用连续介质模型，认为速度和密度均为空间坐标和时间的连续可微函数 6。其中，不可压流体的连续性方程为：()0t+=u (8)Navier-Stokes 方程是牛顿第二定律的运用，核心是动量守恒。因此，Navier-Stokes 方程最基本的假设是把流体看作连续介质，从宏观方位分析时，可以忽略流体的随机分子运动，以描述流体的多个宏观物理量 7。结合连续性方程，考虑重力、表面张力两种体积力的守恒型不可压缩动量方程如下：()()()Tpt+=+uuuuugf (9)式中，p 为正应力，Pa；g 为重力加速度矢量，m/s2；f 为通过连续表面应力模型(Continuum Surface Force,CSF)计算所得的表面张力，N/m。由于相变过程中伴随着能量的传递，需要建立并求解能量方程。能量方程基于热力学第一定律，在OpenFOAM 的多相流求解器中，不可压缩流体的能量守恒方程表达形式如下：()ppTmLTTtCC+=u (10)纪国剑 等 DOI:10.12677/ijfd.2023.112005 48 流体动力学 式中，T 为温度，K；L 为单位质量蒸汽冷凝所释放的相变潜热，J/kg；m 为单位体积蒸发率，kg/(m3s)。2.2.相变模型相变模型 求解器中编译有两种相变模型Lee 模型和 KineticGasEvaporation model 模型，可针对每一类型的相变过程引入特定相变模型。本次研究涉及金属凝固过程，选取 Lee 相变驱动模型。1)由 Lee 模型驱动的相间传质方程 当 C 0 时为熔化/蒸发过程：activateactivateactivate for TTmCTTT=(11)当 C 0 时为凝固/冷凝过程：activateactivateactivate for TTmCTTT=0 时为蒸发过程，当 C 0 时为冷凝过程。2.3.湍流模型湍流模型 在 icoReactingMultiphaseInterFoam 的 turbulenceProperties 中，可以设置用于求解的流动类型。这里有三种不同的模型以供选择，包括 laminar、LES、RAS。LES 和 RAS 是由多个子类模型组成的。其中，LES 分为“等容的 LES 湍流模型”和“各向异性的 LES 湍流模型”10。RAS 分为两类，一类是不可压缩的，一类是可压缩的。本次研究指定 RAS 中的-模型，流动方程见下式。1)湍流动能方程()()DDPDt=+(16)式中，为湍流动能，m2/s2；D为湍流动能的有效扩散系数；P为湍流动能产生率，m2/s3；为湍流动能耗散率，m2/s3。2)湍流动能的耗散率方程 纪国剑 等 DOI:10.12677/ijfd.2023.112005 49 流体动力学 ()()213223CDDPCCDt=+u (17)式中，D为湍流动能耗散率的有效扩散系数；C1、C2和 C3是模型系数。3)湍流粘度方程 2vC=(18)式中，为湍流粘度，m2/s；C为湍流粘度的模型系数。4)OpenFOAM 实现 OpenFOAM 软件对上述方程做了进一步变化，以实现湍流模型在各类求解器中的调用；修正后的湍流动能方程和湍流动能的耗散率方程如下：()()()2fvOptions23DGSSt+=+uu (19)()()()21132fvOptions23DCGCCCSSt+=+uu(20)式中，G 是由于雷诺应力张量的各向异性部分而产生的湍流动能产生率，m2/s3；在 turbulenceProperties字典中，模型系数值分别为 C=0.09、C1=1.44、C2=1.92、C3=0.0、=1.0、=1.3。3.建模建模 3.1.数值计算模型数值计算模型 Figure 1.Numerical model for hydraulic fragmentation of molten metal 图图 1.熔融金属水力碎化数值计算模型 为研究熔融锡液柱水力碎化过程的界面行为及相关物理现象，做出以下假设：熔融锡液柱周围冷却剂不发生相变，引入熔融锡液凝固相变模型，即忽略多相流中水的蒸发冷凝效应，只考虑水力碎化作用。纪国剑 等 DOI:10.12677/ijfd.2023.112005 50 流体动力学 高温熔融锡液柱水力碎化模拟的物理模型及网格划分如图 1 所示。建立二维矩形计算域，计算域的长、宽分别为 700 mm 140 mm，进口处为 20 mm。计算模型的两侧及底部黑色边界类型定义为无滑移、零速度梯度的壁面(wall)，用于约束流体及固体域；顶部两端的灰色边界定义为压力出口(pressure-outlet)，当迭代期间出现回流时会收敛速度更好；顶部的红色边界定义为熔融锡液柱的速度入口(velocity-inlet)，锡液作为不可压缩流体，速度入口可进一步控制流量。整个模拟不考虑空气的夹带效应，流动模型为层流。熔融锡液柱经顶部进口以一定速度流入水域内，如图 1 所示红色域为高温熔融锡液，蓝色域为水；在网格划分时，对锡液柱及液柱附近的网格进行局部加密。在常温常压环境下，锡液初始温度为 1050、初始速度为 3 m/s；计算域内水的初始温度为 27。边界条件根据模拟情况设置如下：进口处压力为calculated 边界条件，压力由其他边界条件及参数来定义；温度与速度为固定值，因此均为 fixedValue 边界条件；相体积分数为 inletOutlet 边界条件，该边界条件由 mixed 边界条件衍生而来，当流体流出时指定为 zeroGradient 条件，当流体流入时则是 fixedValue 条件。出口处速度为 pressureInletOutletVelocity 边界，流体流出时指定 zeroGradient 条件，流入时则根据计算块(patch)法线方向上的通量分配速度，边界条件详细信息见表 1。Table 1.Set up a series of boundary condition 表表 1.边界条件设置 进口 出口 壁面 p calculated zeroGradient zeroGradient T fixedValue zeroGradient zeroGradient U fixedValue pressureInletOutletVelocity noslip alpha inletOutlet zeroGradient zeroGradient 3.2.网格相关性验证网格相关性验证 为准确捕获界面，同时节约计算时间成本，设置三种网格以进行网格相关性研究，这三种四面体网格质量均为 1。监测点设置在水平中心 0.07 m、纵向高度 0.6 m 处；如图 2 所示，比较监测点处温度随时间的变化曲线，确定后续计算分析所用网格单元格总数为 276,202 个。Figure 2.Variation curves of temperature at measuring points according to different grid numbers 图图 2.不同网格数下监测点处温度变化曲线 纪国剑 等 DOI:10.12677/ijfd.2023.112005 51 流体动力学 4.仿真结果与分析仿真结果与分析 4.1.水力碎化过程分析水力碎化过程分析 Figure 3.Simulation diagrams of molten tin fragmentation at diameter is 20 mm and velocity is 3 m/s 图图 3.液柱直径为 20 mm、初速度为 3 m/s 熔融锡液柱碎化相图 如图 3 所示为液柱直径为 20 mm、初速度为 3 m/s 熔融金属锡液碎化相图。在 0.015 s 时，由于顶部相对高压区域及与金属锡和水密度差相关的 R-T 不稳定性的综合作用，金属锡液柱前端发生形变、横向扩张形成伞状；伞状顶部两侧金属细丝向外翻转，受边界层剪切效应发生断裂。由于熔融锡液与水刚开始直接接触，锡液总量小，此时碎化效应并不明显。从 0.015 s 开始，金属锡液柱持续以 3 m/s 的速度流入水中，液柱前端横向扩张的伞状现象因熔融金属锡液柱入水及液柱碎化造成的扰动不再显著，受与速度差相关的 K-H 不稳定波动影响而造成的碎化现象更加突出。在 0.015 s0.039 s 内，熔融锡液柱侧面出现明显的水力不稳定波动，随着该不稳定波持续生长，从 0.03 s 开始可以明显地观察到液柱两侧出现多个不对称的倒钩状金属细丝。金属细丝在 0.033 s 时发生明显的断裂现象，剥离液柱表面。Figure 4.Disturbance near the wire on the metal side 图图 4.液柱侧面金属细丝附近的扰流 纪国剑 等 DOI:10.12677/ijfd.2023.112005 52 流体动力学 图 4 为 0.039 s 时熔融锡液附近的扰流速度矢量图，因 K-H 不稳定性波动形成的倒钩状金属细丝附近往往伴随着流体涡旋流动。该扰流对金属碎化的影响表现在扰流能促使液柱表面金属细丝拉长而加剧细丝断裂、碎化。断裂首先发生在倒钩状细丝尖端，该部分剥离液柱主体的同时，由于熔融锡液内部粘性力的存在，剩余部分将继续发展，再次经历拉长、断裂、碎化的过程，最终使液柱发生整体断裂。熔融锡液遇水后的碎化过程为高密度比(液柱密度与水的密度之比)模拟，对比如图 5 所示陈竞覃等 11高密度比下水银射流碎化过程拍摄图、原理示意图以及如图 6 所示 Iwasawa 等 12低熔点合金液柱射流碎化实验拍摄图。由于模拟不考虑空气夹带和蒸发相变，观察图 5 中液柱入水后侧面的不稳定性波动情况及图 6 中液柱前端的横向扩张和液柱侧面的剥离现象，与数值模拟结果基本一致。Figure 5.Photo and schematic diagram of jet breakup behavior 图图 5.射流碎化过程拍摄图及原理示意图 Figure 6.Photo of jet breakup behavior 图图 6.射流碎化实验拍摄图 Figure 7.Simulation diagrams of molten tins solidification at diameter is 20 mm and velocity is 3 m/s 图图 7.液柱直径为 20 mm、初速度为 3 m/s 时熔融锡液柱凝固相图 纪国剑 等 DOI:10.12677/ijfd.2023.112005 53 流体动力学 如图 7 所示，熔融锡液柱与水直接接触，液柱顶部及两侧率先发生凝固。0.039 s 时对比凝固相图与熔融锡液柱碎化相图，明显发现脱离液柱的锡液滴在失去液柱内部锡液的热量供给后，凝固速度加快，凝固相体积分数较大。图 8 为固态锡体积分数随时间的变化曲线，固态锡生成速率随反应的深入而提升，这是由于锡液碎化剧烈，形成众多小尺寸金属颗粒，这些颗粒储热少，极容易被周围冷却水冷却凝固。Figure 8.Curve for theaveragealphas of solid tin per unit time 图图 8.固态锡体积分数均值随时间的变化曲线 4.2.液柱直径对碎化过程的影响液柱直径对碎化过程的影响 Figure 9.Fragmentation processes of molten tin at diameter is 10 mm 图图 9.液柱直径为 10 mm 时熔融锡液柱碎化过程 Figure 10.Fragmentation processes of molten tin at diameter is 30 mm 图图 10.液柱直径为 30 mm 时熔融锡液柱碎化过程 纪国剑 等 DOI:10.12677/ijfd.2023.112005 54 流体动力学 为研究金属液柱直径对水力碎化的影响，除了进口尺寸 20 mm 外，另设置了 10 mm 和 30 mm 两种进口尺寸，分别进行同一数值工况的水力碎化模拟研究，如图 9、图 10 所示。R-T 不稳定性波动与相对速度无关故而保持不变，因此相图重点关注由 K-H 不稳定性波动造成的液柱侧面的水力碎化现象。当熔融锡液直径为 10 mm 时，0.0135 s 时液柱已经下落至水深 0.19 m 处；而熔融液柱直径为 20 mm 时，下落至水深 0.19 m 处需要耗时 0.039 s。熔融锡液柱在水内下落过程中，由于液柱直径较小，内部粘性力弱，液柱侧面受到 K-H 不稳定性波动的影响，多处出现断裂现象并发生剧烈的碎化。当熔融锡液直径为 30 mm时，由于进口处水压的存在，直径为 30 mm 的液柱下落高度进口处所受阻力较大，熔融锡液柱下落高度耗时 0.056 s 仅下落了 0.13 m，低于相同时间内直径为 10 mm 和 20 mm 时熔融锡液柱下落高度。反应持续至 0.054 s 时，熔融锡液在水下形成直径较为稳定的液柱，液柱两侧受 K-H 不稳定性波动的影响，两侧出现明显的弯钩状金属细丝并断裂，液柱碎化现象较之前显著提升。因此，液柱入水初速度相同时，随着液柱直径的增大，液柱入水用于克服入水阻力的动能消耗越大，在相同时间内液柱入水深度越小，如图 11 所示。同时，相同入水时间内，液柱化程度与液柱的直径成反比，液柱直径越小，熔融锡液柱碎化得越剧烈。Figure 11.Curve for height of the top of liquid column entering water with time 图图 11.液柱顶部入水高度随时间的变化曲线 4.3.液柱入水初速度对碎化过程的影响液柱入水初速度对碎化过程的影响 Figure 12.Fragmentation processes of molten tin at velocity is 5 m/s 图图 12.初速度为 5 m/s 时熔融锡液柱碎化过程 纪国剑 等 DOI:10.12677/ijfd.2023.112005 55 流体动力学 Figure 13.Fragmentation processes of molten tin at velocity is 10 m/s 图图 13.初速度为 10 m/s 时熔融锡液柱碎化过程 为研究金属液柱下落速度对水力碎化的影响，除了初速度 3 m/s 外，另设置了 5 m/s 和 10 m/s 两种进口速度，分别进行同一数值工况的水力碎化模拟研究，如图 12、图 13 所示。对比分析初速度为 3 m/s、5 m/s、10 m/s 时熔融金属液柱碎化过程，当熔融锡液柱与水直接接触反应至 0.02 s 时，熔融锡液柱入水初速度增大意味着熔融液柱入水时的动能越大，液柱克服顶部高压及水的阻力作用向下流动过程中，初速度大的液柱顶部碎化更剧烈。同时，由于单位时间内流入水内的锡液总量增加，在液柱直径相同的情况下，热量传递速率加快，水力扰动变大，加速了熔融锡液柱的碎化。图 14 为固体锡初速度 10 m/s 时，同一时间段内固态锡生成量高于其他两种初速度下的固态锡生成量，即初速度为 10 m/s 的锡凝固速度较快。因此，熔融锡液柱入水后的碎化程度与液柱入水时的初速度成正比，液柱入水时初速度越大，水力波动越剧烈，熔融锡液柱顶端及侧面的碎化程度得越高，锡凝固速度越快。Figure 14.Curves for the averagealphas of solid tin per unit time at different velocities 图图 14.不同入水速度下固态锡体积分数均值随时间的变化曲线 5.结论结论 熔融锡液水力碎化模拟研究得到以下结论：1)水力碎化现象方面：熔融锡液柱碎化是三种水力因素共同作用的结果。液柱顶部受与密度相关的瑞利泰勒(Rayleigh-Taylor,R-T)不稳定性影响呈伞状向两端横向发展，液柱侧面受与相对速度相关的开纪国剑 等 DOI:10.12677/ijfd.2023.112005 56 流体动力学 尔文亥姆霍兹(Kelvin-Helmholtz,K-H)不稳定性影响生成多个不对称的倒钩状金属细丝，受界面剪切力的影响液柱伞状顶部横向发展的两个尖端及液柱侧面倒钩状金属细丝被拉长并断裂。2)水力碎化剧烈程度与液柱直径、液柱入水初速度的敏感性分析：液柱入水初速度相同时，随着熔融锡液柱直径的增大，相同时间内液柱入水深度越小；液柱碎化剧烈程度与液柱直径成反比。液柱入水直径相同时，熔融锡液柱入水后的碎化程度、锡凝固速度与液柱入水时初速度成正比。基金项目基金项目 国家重点研发计划项目(2017YFC0805101)；江苏省绿色过程装备重点实验室开放课题基金资助项目(GPE201901)；2022 年江苏省研究生科研与实践创新计划项目(SJCX22_1434)。参考文献参考文献 1 周源.蒸汽爆炸中熔融金属液滴热碎化机理及模型研究D:博士学位论文.上海:上海交通大学,2014:3-131.2 袁明豪.蒸汽爆炸相关现象机理的数值模拟研究D:博士学位论文.上海:上海交通大学,2008:1-106.3 杨志.钠冷快堆堆芯熔融材料-液态钠相互作用破碎特性及机理研究D:博士学位论文.哈尔滨:哈尔滨工程大学,2020:19-110.4 钟明君.熔融物与冷却剂相互作用过程及熔融物水力学碎化行为的数值模拟研究D:博士学位论文.上海:上海交通大学,2016:1-121.5 沈正祥,李金柱,吕中杰,等.水动力作用下低速射流碎化的数值模拟J.计算力学学报,2013,30(2):308-312.6 商晓波.湿式离合器接合过程理论分析及试验研究D:硕士学位论文.杭州:浙江大学,2019:8-64.7 毛菊珍.欧拉模型中多烟雾自由融合实时绘制算法研究D:硕士学位论文.秦皇岛:燕山大学,2016:8-16.8 Zientara,M.,Jakubczyk,D.,Litniewski,M.,et al.(2013)Transport of Mass at the Nanoscale during Evaporation of Droplets:The Hertz-Knudsen Equation at the Nanoscale.The Journal of Physical Chemistry C,117,1146-1150.https:/doi.org/10.1021/jp3091478 9 金明善,李文佐,宫宝安,等.Clausius-Clapeyron 方程在解题中的应用J.大学化学,2011,26(1):3.10 Roohi,E.,Zahiri,A.P.and Passandideh-Fard,M.(2013)Numerical Simulation of Cavitation around a Two-Dimen-sional Hydrofoil Using VOF Method and LES Turbulence Model.Applied Mathematical Modelling,37,6469-6488.https:/doi.org/10.1016/j.apm.2012.09.002 11 陈竞覃,周源,汪杨乐,等.气腔夹带作用下射流碎化过程实验研究J.核动力工程,2019,40(S2):26-30.12 Iwasawa,Y.,Abe,Y.,Kaneko,A.,et al.(2012)Jet Breakup Behavior with Surface Solidification.2012 20th Interna-tional Conference on Nuclear Engineering and the ASME 2012 Power Conference,Anaheim,July 30-August 3,2012,175-184.https:/doi.org/10.1115/ICONE20-POWER2012-54412