水下近自由面爆炸过程中的多相流运动特性研究_余俊.pdf

水下近自由面爆炸过程中的多相流运动特性研究余俊1，2，王海坤1，2，刘国振1，郝轶1，张得扬1（1.中国船舶科学研究中心，江苏 无锡 214082；2.深海技术科学太湖实验室，江苏 无锡 214082）摘要：水下近自由面爆炸现象涉及到爆轰气体产物、水以及空气等多种流体之间的相互作用过程，同时还伴随有明显的汽液两相之间的相变转换。为了理解和掌握近水面爆炸过程中的多流体非线性耦合特性，引入基于相变转换的多相可压缩流体计算模型进行分析。采用一维激波管问题对计算模型的精确性和收敛性进行测试验证，在此基础上利用二维轴对称模型模拟近水面爆炸过程，获得冲击波传播、空化产生与溃灭、爆炸气泡运动等过程。爆炸气泡运动的计算结果与试验结果具有一致性。对于由水下爆炸引起的近自由面空化现象，数值模拟清晰展示了空化域的产生、发展及其溃灭的完整过程。空化域由开始的单连通域演化为涡环形态，并向外扩展直至完全溃灭。将蒸汽体积分数0.5作为空化域的识别判据时，计算获得的空化域的演化过程与试验结果较为一致。数值模拟结果还进一步揭示了空化域内部的一些重要的物理特性，本文的研究结果能够为水面舰艇结构的抗爆抗冲击防护研究提供技术支撑。关键词：水下爆炸；爆炸气泡；空泡；可压缩流体；相变转换；空泡溃灭中图分类号：O359 O383文献标识码：Adoi:10.3969/j.issn.1007-7294.2023.02.001Motion characteristics of multiphase flow inunderwater explosion near free surfaceYU Jun1,2,WANG Hai-kun1,2,LIU Guo-zhen1,HAO Yi1,ZHANG De-yang1(1.China Ship Scientific Research Center,Wuxi 214082,China;2.Taihu Laboratory ofDeep-sea Technological Science,Wuxi 214082,China)Abstract:Underwater explosion(UNDEX)near free surface involves the interaction between detonation gasproducts,water,air and other fluids,and also contains obvious phase transition.This interaction is a typicalmultiphase compressible fluid movement problem with phase effect,which is very nonlinear and difficult tostudy.In order to understand and master the nonlinear interaction characteristics of multi-component in underwater explosion near free surface,a compressible multi-phase fluid model based on phase transition wasintroduced.The control system equations of the model were solved by a simple fractional step approach,which was divided into hyperbolic step and relaxation step.A MUSCL-Hancock method with HLLC approximate Riemann solver was applied to discretize the homogenous hyperbolic equations in hyperbolic step,andthe nonlinear system was solved by Newton-Raphson iterative method in the following relaxation step.Theaccuracy and convergence of the numerical model were tested and verified by using the 1D shock tube example.On this basis,underwater explosion near free surface was simulated by using the 2D axisymmetric model,and the process of shock wave propagation,cavitation generation and collapse,and explosion bubble motion第27卷第2期船舶力学Vol.27 No.22023年2月Journal of Ship MechanicsFeb.2023文章编号：1007-7294（2023）02-0163-10收稿日期：2022-08-21作者简介：余俊（1984-），男，硕士，高级工程师，通讯作者，E-mail:；王海坤（1980-），男，博士，研究员；刘国振（1976-），男，硕士，高级工程师；郝轶（1983-），男，硕士，高级工程师；张得扬（1983-）男，硕士。were clearly captured.For the explosion bubble movement,the numerical results of bubble interface are ingood agreement with the experimental results.For the cavitation near free surface caused by underwater explosion,numerical simulation clearly shows the generation,development and collapse of the cavitation.Thecavitation domain evolves from the initial single connected domain to the vortex ring form and expands outward until it completely collapses.When the vapor volume fraction is greater than 0.0005(the identificationcriterion of cavitation domain),the numerical results of cavitation evolvement are consistent with the experimental results.The numerical results further reveal some important physical properties in the cavitation domain,which can provide some technical support for the research of anti-explosion and anti-impact protection of surface warship structures.Key words:underwater explosion;explosion bubble;cavitation;compressible fluid;phase transition model;cavitation collapse0 引言随着现代海战武器的高速发展，实战条件下舰船遭受的威胁越来越大，对舰船水下冲击与防护性能的研究受到了高度的重视。水面舰艇除了会遭受船底爆炸冲击损伤之外，近水面爆炸也会产生重要的威胁1。近水面爆炸现象的本质特征属于多相流体（气-液-汽）相互作用及其与结构的耦合运动过程，其非线性耦合效应主要体现在两个重要方面：首先是爆轰气体产物、水、空气等多组分流体之间的对流运动,其可归纳为多相可压缩流体界面捕捉问题；其次是自由面的存在使得流场中会出现稀疏波的传播，水中会出现典型空化现象,可以归结为相变问题。要想准确理解和掌握近水面爆炸过程中的多流体非线性耦合效应，需要同时重点处理好多相可压缩流体界面捕捉和相变问题。目前国内外已经对近水面爆炸现象进行了大量的研究2-4。在实验研究方面，Marcus等5开展了水下近自由面爆炸试验研究，发现水下爆炸产生的片空化区域在闭合溃灭时会产生较大的压力，获得了近水面爆炸过程中包括空化效应的冲击波和气泡阶段的压力曲线;Brett等6开展了一系列水下圆柱壳近距离爆炸试验，对冲击波和空化联合加载下圆柱壳的变形进行了测量分析，发现空化载荷的幅值与冲击波幅值相当，而且空化溃灭载荷造成的壳体变形要超过整个冲击波阶段产生的总变形；Kleine等7利用纹影法拍摄到了近水面爆炸试验过程中液体空化的密度变化过程图象，较为清晰地展示了爆炸气泡膨胀以及空化域内部演化图像；Cui等8开展了一系列不同边界条件下近水面小药量药包爆炸试验，获得了爆炸气泡与自由面运动的完整过程。由于实验研究获得的数据有限，在稳定性和重复性上难以满足要求，难以获得流场的精细特征。数值模拟能够较为全面地展示流场的精细特征，使得其在近水面爆炸研究中具有一定的优势。这里着重介绍基于多相可压缩流体计算模拟的研究情况，对于基于声学单元和边界单元的数值方法不作讨论。自从Aanhold等9提出cut-off空化模型之后，由于其使用简单方便，被广泛应用在水下近自由面或近结构爆炸分析中。Wardlaw等10利用cut-off模型分析了充水圆筒中心药包爆炸过程中冲击波在圆筒内部产生的局部空化现象，获得了空化溃灭载荷的典型特征；Shukla11、Yu12等采用cut-off空化模型与五方程相结合的方法对带有空化效应的近水面爆炸多流体运动问题进行了研究；Wu等13采用cut-off模型与流体局部间断迦流金（LDG）法相结合分析了水下远场冲击波与自由面作用产生空化的过程。为了克服cut-off模型无法描述空化域内液相和汽相组分的问题，Liu14、Schmidt15和Xie16等相继提出了isentropic模型、Schmidt模型和Modified Schmidt模型。这些模型本质上都属于单流体空化模型，在流体压力降低到饱和态压力之前均视为纯流体，当降到饱和态压力之后，各种模型的具体差别体现在空化域内汽液混合流体状态方程的构造方法不同。除了上述单流体空化模型之外，双流体模型也在不断发展，主要以Chiapolino17-18、Saurel19、Pelanti20-21等提出的四方程、六方程等模型为主，该模型通过持续追踪空化质量分数的演化来捕捉空化的产生、164船舶力学第27卷第2期发展与溃灭过程。双流体模型则将初始流体视为液相和汽相的某种混合，在流场发生改变时（如冲击波传播和气泡运动等包含间断时）液-汽两相之间会发生对流运动，同时还伴随有质量和热量交换的相变过程，空化被视为流体区域内部蒸汽含量不断累积的过程。因此双流体模型描述的空化演化过程更符合物理本质和思维逻辑，但双流体模型相比单流体模型，双流体计算过程非常复杂繁琐，导致目前水下爆炸领域内绝大部分的空化研究都采用单流体模型。本文拟采用Chiapolino等18提出的基于多相可压缩流体四方程的双流体空化模型来重点研究近水面爆炸过程中的多相流运动。首先对计算模型的控制方程及其数值方法进行了简要介绍，然后采用一维激波管问题对计算模型中的相变转换以及激波捕捉功能进行了考核，最后针对近水面爆炸现象，将计算结果与试验结果进行对比，并重点分析空化域内部的典型物理特性。1 计算模型简介1.1 控制方程不考虑粘性和热传导效应的多相可压缩流体运动方程可表示为18|t+()u=0()ut+()u u+pI=0()Et+()()E+p u=0Ykt+()Yku=0(1)式中，、u、p和E分别代表混合流体的密度、速度矢量、压力和单位质量总能，E=e+0.5uu，这里e为单位质量内能；Yk代表混合流体中各相介质的质量分数。控制方程（1）中参数k代表了混合流体中各相成分，为了计算模型中表达简便起见，这里约定如下：k=1代表液相成分；k=2代表蒸汽相成分；k=3,，N代表其它不发生相变的气相或液相成分。本文中令k=3代表空气（不可冷凝），k=4代表爆轰气体产物。在流体计算模型中普遍采用的是SG（Stiffed-Gas）状态方程，该方程能够描述多种液相或者气相流体的运动特性。近年来NASG（Nobel-Able Stiffed-Gas）状态方程逐渐被采纳，该方程修正了SG方程对液相介质中分子排斥力效应描述的不足22。该方程形式为|pk()k,ek=()k-1()ek-qkk-bk-kpkTk()pk,k=()pk+pk()k-bk()k-1 Cv,kgk()pk,Tk=()kCv,k-qk Tk-Cv,kTklnTkk()pk+pkk-1+bkpk+qkck()pk,k=k2k()pk+pkk-bk(2)式中，k=1/k、Tk、gk、ck分别代表k相流体介质的比体积、温度、吉布斯（Gibbs）自由能与声速，参数k、pk、Ck、qk，qk、bk等可以通过流体热动力学试验中与温度相关的饱和态压力、各相比体积、各相的焓等数据来拟合得到22。在饱和平衡态下液相与汽相的Gibbs自由能相等，即g1=g2，化简可得饱和态下温度第2期余俊等：水下近自由面爆炸过程中的 165与压力的关系式：ln(psat+p2)=As+Bs+Es psatTsat+Cs ln(Tsat)+Ds ln()psat+p1(3)式中下标“sat”代表饱和态，参数As、Bs、Cs、Ds、Es分别为As=Cp,1-Cp,2+q2-q1Cp,2-Cv,2,Bs=q1-q2Cp,2-Cv,2,Cs=Cp,2-Cp,1Cp,2-Cv,2,Ds=Cp,1-Cv,1Cp,2-Cv,2,Es=b1-b2Cp,2-Cv,2(4)式中Cp,k表示第k相流体定压比热容。液态水、水蒸汽以及空气的NASG状态方程参数如表1所示，其中W为物质的摩尔质量。本文算例涉及到水中空化问题所采用的计算参数均采用表1中的参数22。表1 液态水、水蒸气和空气的NASG状态方程参数Tab.1 NASG coefficients for liquid,vapor and air名称Cp/(Jkg-1K-1)Cv/(Jkg-1K-1)p/Pa液态水428536101.197.028E8汽态水14019551.470空气10077191.40名称q/(Jkg-1)q/(Jkg-1K-1)b/(m3kg-1)W/(gmol-1)液态水-1 177 78806.61E-418汽态水2 077 61614 317018空气000291.2 相变转换模型控制方程（1）并未考虑到各相之间发生的相变过程，对于流体中发生液相与其对应汽相之间的物质与热量转换，可以参考化学反应过程中判断反应方向的判据，认为在达到平衡态之后液相及其蒸汽相之间的吉布斯自由能相等。同时结合控制方程（1）的假定，认为系统达到平衡态后满足如下关系式18：|g1=g2,T=Tk,p=pk=k=1NYkk,e=k=1NYkek(5)将公式（2）和（3）代入，得ppartial=Y*2/W2Y*2/W2+3NYk/Wkp*=psat()T*(6)=Y*11()T*,p*+Y*22()T*,p*+k=3NYkk()T*,p*(7)e=Y*1e1()T*,p*+Y*2e2()T*,p*+k=3NYkek()T*,p*(8)这里ppartial代表气体混合物中蒸汽相的分压，由公式可知其与蒸汽相的摩尔质量分数成正比18。非线性方程（6）（8）可以采用迭代方法进行求解，如Newton-Raphson方法。1.3 数值离散控制方程（1）与状态方程（2）以及空化相变方程（6）（8）共同构成了考虑相变效应的多相可压缩流体控制方程组，可以利用相变松弛法则进行分步求解23。第一步求解齐次双曲型方程组（1），需要结合状态方程（2），利用二阶MUSCL-Hancock方法以及HLLC近似黎曼求解器。通过该步获得中间变量（v,e,p,T,Yk）；第二步采用 Newton-Raphson 迭代方法求得相变转换平衡态时的状态量（p*,T*,Yk=1,2）。2 计算与讨论2.1 一维激波管测试与验证本算例主要用来与Chiapolino等相同计算工况18进行对比验证。其中整个流场的蒸汽相初始质166船舶力学第27卷第2期量分数设置为0.2，液相与气相的质量分数分别为0.1和0.7，计算域为0,1 m。流场中密度和压力的初始间断面位于x=0.5 m，左、右侧压力分别为2105Pa和105Pa，两侧密度分别为1.94 kg/m3、1.02 kg/m3。图1为t=1 ms时刻流场中密度、压力、速度、温度、液相与蒸汽相质量分数分布曲线。其中P-T和no P-T曲线分别代表是否考虑相变转换时的计算工况，initial代表初始条件，Chiapolino代表对比的数据（网格数为100）。同时为了考察计算模型的网格收敛性，分别采用100、200和400个网格单元均匀划分计算域。由密度和压力分布曲线可知，考虑液相和汽相之间的相变转换之后，向右传播的压缩波和向左传播的稀疏波的波速都稍有下降。由温度曲线可知，向右传播的压缩波波后温度在两种工况下都有所升高，但波后稳定区的温度在考虑相变转换时为346.3K,而不考虑相变时为362.6 K，差异较大。同理，向左传播的稀疏波波后稳定区的温度都有所下降，但是考虑相变转换后为344.7 K，而不考虑时则为329.4 K。这是由于在压缩波传播的波阵面以及稀疏波传播过程中液相与汽相发生了质量和热量交换，在达到平衡态的过程中温度存在变换。而汽-液两相的质量转换可以参考质量分数曲线得知，压缩波过后流体中液相发生了蒸发，使得汽相质量分数增加。该计算结果与通常理解的冲击波作用后由于压力升高容易引起液相发生冷凝的认识相矛盾，这里需要结合温度变化进行联合分析。由于考虑了相变转换，冲击波作用后波后温度由初始的337.5 K上升到346.3 K,虽然波后压力也由初始的105Pa上升到1.4105Pa，但是从饱和态曲线P-T（公式（3）来分析会发现温度升高引起的蒸发效应要大于压力升高带来的冷凝效应，最终效果就是压缩波后处于蒸发状态。同理，向左传播稀疏波的最终效果是波后处于冷凝状态。由图1可知，随着网格数的增加，计算结果的收敛性较好，本文的计算结果与Chiapolino的结果吻合较好。IIIIII图1 1 ms时刻流场各物理量分布曲线图Fig.1 Diagram of different variables in fluid at t=1 ms2.2 水下近自由面爆炸现象本节利用上面介绍的计算模型与Cui等8的试验观察结果进行对比。该试验在边长2 m的立方体水槽内开展，采用等效于5.2 g TNT的药包起爆，药包中心位于水下0.13 m处。药包采用等效爆轰模型近似处理，根据Geer-Hunter理论模型24以及数值经验，初始球形爆轰气体产物半径设置为R0=0.009 m，爆轰气体产物内部密度和压力分别为1606 kg/m3和109Pa，采用理想气体状态方程，=1.8，CV=695，W=290。空气、水和爆炸气泡内的初始蒸汽质量分数分别为10-9、10-8和10-6。空气、水和爆炸气泡的初始温度分别设置为295 K、295 K和1120 K,其中空气和爆轰气体产物处于过热状态，水处于饱和状态。空气、水和爆炸气泡内的初始蒸汽体积分数分别为 1.5610-9、1.3910-5和8.0710-7。为了节省计算资源，采用二维轴对称模型进行模拟，如图2所示。其中y轴为对称旋转轴，y=0为初始第2期余俊等：水下近自由面爆炸过程中的 167水平面位置，药包位于(0,-0.13)m处。计算域采用非均匀结构性网格划分，如图2所示（图中显示的网格密度为实际的1/2）。其中区域0,0.4-0.6,0.2 m2内采用均匀网格，dx=dy=0.001 m，区域外采用等比数列格式的递增步长网格，等比为1.02。整个计算域网格为6501200，计算域大小为0,7.34-7.5,1.1 m2。y轴边界设置为对称条件，其它如空气层上方、水域下方以及右侧边界均设置为无反射条件，CFL设置为0.5。图 3 显示了在考虑相变条件下 t=0.084ms、0.168 ms、0.334 ms、0.5 ms、0.666 ms、1.0ms、4.176 ms和6.176 ms时刻流场中的密度、压力和蒸汽体积分数分布云图，其中所有图片采用统一大小的观察窗，尺寸为-0.6,0.6-0.7,0.35 m2。由于图3中的各物理量的内部分布跨度非常大，未单独列出各云图的刻度线，压力云图中的黑色虚线表示对应时刻爆炸气泡界面的位置。0.084 ms时刻，初始冲击波刚到达自由面不久，其产生的反射稀疏波使得自由面下方附近出现蒸汽聚集，图中对应时刻蒸汽体积分数云图显示的局部高亮区域内数值范围为0.42.5。随着稀疏波不断向下传播，0.168 ms时刻稀疏波已经达到爆炸气泡上壁面，并反射压缩波，可以发现此时爆炸气泡表面附近的局部范围处于相对高压区。随着稀疏波在爆炸气泡表面反射成压缩波以及自由面不断反射稀疏波，旋转轴上的蒸汽含量不断降低（对应空泡溃灭过程），而旋转轴外侧的蒸汽含量不断增加（对应空化产生过程），如图3中的0.334 ms和0.500 ms所示。随着空泡的不断产生和溃灭，空泡中心区域逐渐变薄，并向外扩张，如图3中的0.666 ms和1.0 ms所示。在4.167 ms和6.167 ms时刻空泡基本上完全消失，但爆炸气泡一直处于向外膨胀过程，并不断抬高水面位置。在整个过程中空化由开始的单连通域不断演化为多连通域，显示出涡环形状，涡环逐渐向外扩展，并逐渐变薄直至最终完全消失。0.084 ms0.168 ms0.334 ms0.500 ms图2 近水面爆炸二维轴对称计算模型与网格划分示意图Fig.2 Schematic of 2D axisymmetric model and mesh size ofunderwater explosion near free surface密度压力蒸汽体积分数密度168船舶力学第27卷第2期0.666 ms1.000 ms4.167 ms6.167 ms图3 典型时刻流场的密度、压力和蒸汽体积分数分布云图Fig.3 Contour plots of the density,pressure and vapor volume fractions at typical times在上述整个时间段内，自由面上方附近的空气中蒸汽含量一直比较高，并且包含蒸汽的区域不断扩展。这是由于水中初始冲击波作用自由面后向空气透射压缩波、水中的液相和汽相同时向空气中扩散的结果，在扩散过程中仍然存在液-汽两相之间的相变转换。这是因为水面上方初始空气中是不包含液相水，而且水蒸汽含量极低（10-9），而图3中蒸汽体积分数云图中红色区域代表自由面上方，其内部的液相水的质量分数范围在0.60.8，水蒸汽质量分数范围在0.10.3，剩下的为空气。文献8中只记录了通过观察窗获得的部分流场图像，如图4中第一行所示。试验获得的图片的中间纵剖面尺寸与数值模拟中的-0.4,0.4-0.68,0.2 m2基本相当。在第一行的试验图片中左下角两位数值代表图片序号，右下角为记录的时刻（单位为ms），其中0.167 ms代表装药起爆时刻。由图4可知，小药量装药在近水面爆炸时在水面下方产生的空化域形态与“云空化”非常相似，水中的空泡体积非常小，形成云带状，其发展过程中能够观察到明显的涡环形状。为了与试验进行对比，需要对“云空化”区域内的蒸汽含量进行判断，选取蒸汽体积分数是一个比较合适的对象。这里假定当流体单元中的蒸汽体积分数大于0.5时就标记为空化单元，单元所在的区域即为空化域。有关空化域内蒸汽含量的判据目前仍然没有统一的标准，这里选取的0.5只是数值计算上的一个参考。将图 3中蒸汽体积分数云图转换为空化域纹影图，如图4第二行所示。其中红色区域为水下空化域，黄色区域为自由面上方的高含量蒸汽集中区域，绿色为空气，蓝色为水域，浅蓝色为爆炸气泡。经过对比可知，计算获得的水中空化域与试验观察到的“云空化”范围基本相当，两者都清晰地展示了空化域形成涡环，并逐渐向外扩张和变薄的过程，最终完全溃灭。图5显示了空化域内的最大、最小和平均温度变化以及空化域总体积时程曲线。由图 5可知空泡总体积的膨胀过程与收缩过程的时间基本相当，整个空化运动周期约压力蒸汽体积分数0 ms0.168 ms0.334 ms0.500 ms0.666 ms1.000 ms4.167 ms6.167 ms图4 典型时刻的试验与计算结果对比Fig.4 Comparison of experimental and numerical resultsat typical times第2期余俊等：水下近自由面爆炸过程中的 169为1.14 ms。空化域膨胀和收缩过程中内部最高温度变化范围为294.6309.7 K，最低温度变化范围为283.7294.6 K。膨胀到最大体积时空化域内部的最高和最低温度分别为309.7 K和283.7 K，但在空化域变化过程中其内部体积平均温度基本恒定为294.6 K，与室温相当。图6为流场中三个坐标点处压力时程曲线图，分别为水下0.05 m处距离爆心不同水平距离处的测点。由图可知，三个测点处的空化溃灭过程中最大压力分别为2.15 MPa、0.98 MPa和0.92 MPa。图7 为空化域内最大、最小和平均压力变化时程曲线。由图可知，空化域内最大压力在体积膨胀阶段的大部分时刻保持在0.018 MPa左右，而在收缩阶段则主要在 0.018 MPa 上方振荡，最大值达到0.029 MPa。空化域内最小压力在中间的绝大部分时间内变化较为平缓，在32004800 Pa范围内变动，只是在空化刚开始产生和最后溃灭阶段变化梯度较大。空化域内体积平均压力变化范围为900018 000 Pa，其变化趋势与最小压力变化趋势基本保持一致。由此可见，空化域内压力并不是维持在饱和态压力不变，其具有分布不均匀性、变化梯度大等特点。同时结合图5可知其内部温度的变化范围在283.7309.7 K，空化域内部发生较为明显的相变转换现象。图8为空化域内的蒸汽相体积分数的最大和平均值变化时程曲线，空化域内蒸汽相体积分数的最小值为前面设定参考值的0.5。由图8可知，空化域内的蒸汽相体积分数在0.5173，其峰值并未出现在空化体积膨胀到最大的时刻。尽管空化整个运动过程中蒸汽相最大体积分数并不低，但是空化域内部整体体积平均值变化区间在0.510，说明蒸汽的平均含量比较低，这与前面介绍的实验中观察到的“云空化”较为一致，同时也说明了本文计算模型具有一定的合理性和精确性。图9为空化域内密度的最大、最小和平均值变化时程曲线。由图9可知，空化域内的最大密度变化范围为10511061 kg/m3，体积平均密度变化范围为10391051 kg/m3，最小密度变化范围8601061 kg/m3，其峰值也未出现在空化体积膨胀到最大的时刻。结合图8可知，空化域内蒸汽相的含量总体上非常小，对于流体密度降低的程度有限，空化域内总体体积平均密度变化较小。图5 空化域内最大、最小和平均温度以及空化域体积的时程曲线Fig.5 Time history of the maximun,minimun,averagetemperatures and cavitation domain volumes incavitation domain0.00.51.01.52.00102030400102030400.40.81.2012012p(MPa)t(ms)(0.0,-0.05)m(0.2,-0.05)m(0.4,-0.05)m图6 三个测点处压力时程曲线图Fig.6 Time history of pressure at threemeasuring points图7 空化域内最大、最小和平均压力时程曲线Fig.7 Time history of the maximum,minimum andaverage pressures in cavitation domain170船舶力学第27卷第2期3 结论针对近水面爆炸过程中出现的爆轰气体产物、水以及空气等多流体复杂运动，以及各种流体内部可能出现的液相水和汽相水之间的相变转换等现象，本文采用了基于相变效应的多流体运动计算模型进行处理，获得了如下结论：（1）利用一维激波管问题对计算模型精确性和收敛性测试表明，本文的模型计算结果与Chiapolino的结果吻合较好，能够精确地捕捉冲击波与稀疏波的传播过程，展示了冲击波作用之后流体蒸发以及稀疏波作用之后流体冷凝的过程，最终的状态往往是两种过程的综合作用的结果；（2）近水面爆炸过程中爆炸气泡运动过程的计算结果与试验结果符合较好，说明本文建立的多流体耦合计算模型能够捕捉多种界面之间的复杂运动过程。（3）多流体计算模型中引入了相变转换，可获得近水面爆炸过程中水面下方的液体内部发生的空化现象；采用蒸汽体积分数0.5作为空化域识别的判据，数值模拟获得的空化域范围与试验结果较为一致，空化域都是由开始的单连通域演化为涡环形态，并向外扩展直至完全溃灭。（4）数值模拟结果发现，空化域在运动过程中其内部蒸汽的体积分数范围为0.5173，但内部平均值的范围仅为0.510，说明蒸汽的平均体积含量比较低，这与实验中观察到的“云空化”现象具有一致性。同时该现象也可以从空化域内部密度的平均体积含量来得到印证，其变化范围为10391051 kg/m3，说明空化域内部平均密度与常态下液相水的密度非常接近。另外，空化域内的压力变化范围为320029 000 Pa，温度变化范围为283.7309.7 K，说明空化域内部的压力和温度一直处于动态变化之中，并不是维持在某个恒定的饱和态不变。近水面爆炸过程中除了存在多流体之间的非线性耦合效应，还会与周围结构发生流固耦合作用，其过程更加复杂，将在后续计划中研究。参考文献：1 刘建湖.舰船非接触水下爆炸动力学的理论与应用D.无锡:中国船舶科学研究中心,2002.Liu Jianhu.Theory and its applications of ship dynamic responses to non-contact underwater explosionsD.Wuxi:ChinaShip Scientific Research Center,2002.(in Chinese)2 盛振新,刘建湖,张显丕,等.水下爆炸气泡射流载荷阵列测量技术探索J.爆炸与冲击,2021,41(3):52-59.Sheng Zhenxin,Liu Jianhu,Zhang Xianpi,et al.On an array-sensor technology for measuring bubble jet load generated byunderwater explosionJ.Explosion and Shock Waves,2021,41(3):52-59.(in Chinese)222图8 空化域内蒸汽体积分数2最大值、平均值变化曲线Fig.8 Time history of the maximum and average vaporvolume fractions2in cavitation domain图9 空化域内密度的最大值、最小值和平均值的时程曲线Fig.9 Time history of the maximum,minimum andaverage densities in cavitation domain第2期余俊等：水下近自由面爆炸过程中的 1713 余俊,刘国振,汪俊,等.水下爆炸气泡射流载荷计算的新方法研究J.计算力学学报,2021,38(1):120-125.Yu Jun,Liu Guozhen,Wang Jun,et al.Research on a new method of bubble jet load in underwater explosionJ.ChineseJournal of Computational Mechanics,2021,38(1):120-125.(in Chinese)4 金键,朱锡,侯海量,等.大型舰船在水下接触爆炸下的毁伤与防护研究综述J.爆炸与冲击,2020,40(11):11-36.Jin Jian,Zhu Xi,Hou Hailiang,et al.Review on the damage and protection of large naval warships subjected to underwatercontact explosionsJ.Explosion and Shock Waves,2020,40(11):11-36.(in Chinese)5 Marcus M H.The response of a cylindrical shell to bulk cavitation loadingR.NSWC TR 81-295,1983.6 Brett J M,Yiannakopoulos G,van der Schaaf P J.Time-resolved measurement of the deformation of the submerged cylinders subjected to loading from a nearby explosionJ.International Journal of Impact Engineering,2000,24(9):875-890.7 Kleine H,Tepper S,Takehara K,et al.Cavitation induced by low-speed underwater impactJ.Shock Waves,2009,19:895-900.8 Cui P,Zhang A M,Wang S P.Small-charge underwater explosion bubble experiments under various boundary conditionsJ.Physics of Fluids,2016,28:117103.9 Aanhold J E,Meijer G J,Lemmen P P M.Underwater shock response analysis on a floating vesselJ.Shock and Vibration,1998,5:53-59.10 Wardlaw A B,Luton J A.Fluid-structure interaction mechanisms for close-in explosionsJ.Shock and Vibration,2000,7:265-275.11 Shukla R K,Carlos P,Jonathan B F.An interface capturing method for the simulation of multi-phase compressible flowsJ.Journal of Computational Physics,2010,229(19):7411-7439.12 Yu J,Liu G Z,Wang J,Wang