
Citation: | ZHU Wenrui, WU Xingxing, LIU Jianhu, WANG Jun, ZHAO Yanjie, LI Tianran. Failure of Square Plate under the Influence of Boundary Conditions Subjected to Shock Loading[J]. Chinese Journal of High Pressure Physics, 2021, 35(1): 014202. doi: 10.11858/gywlxb.20200565 |
不同流体间的扰动界面在冲击波作用下会发生不稳定性增长,称为Richtmyer-Meshkov(RM)不稳定性[1-2]。当RM不稳定性发展进入非线性阶段后,扰动界面会发展为“蘑菇”状的尖钉/气泡结构;之后由于界面两侧切向速度不匹配而产生剪切效应,导致界面上又发生Kelvin-Helmholtz(KH)不稳定性[3],并诱发更多小尺度涡。KH不稳定性的出现会破坏“蘑菇”状界面结构的对称性,使之破碎,并加速湍流混合。流体界面不稳定性和湍流混合涉及许多重要领域,如武器物理、惯性约束聚变、超音速燃烧、天体物理等,因此具有重大的研究意义。
较早之前,对于RM界面不稳定性和湍流混合,研究人员关注的重点是湍流混合区宽度的增长规律及影响因素[4-16]。随着实验诊断技术、高精度数值算法及并行计算技术的发展,研究重点逐渐深入到湍流混合的内在机制和统计规律[17-23]以及以此为基础的湍流混合建模[24]。然而,由于湍流混合的多尺度特性,实验中所能获取的流场信息依然十分有限,无法给出足够多的湍流混合区瞬态数据和精细复杂的混合区流场结构;而数值模拟基本上能够给出所需要的有关湍流混合区流场的所有信息,特别是借助超大规模的高性能计算,可以在不同尺度上研究湍流混合的物理机制及其演化的统计特性。Cabot等[25]就曾在IBM的BlueGene/L超级计算机上利用6万多个核实现了Rayleigh-Taylor(RT)不稳定性及湍流混合的近300亿网格(3 0723)的超大规模直接数值模拟,研究了湍流混合区增长的标度率问题和雷诺数效应。
AWE开展的“反尖端”RM不稳定性激波管实验中初始的SF6/空气右界面呈“∨”形,整体上类似于反的士兵臂章,称为“反尖端”界面。本研究采用自研的可压缩多介质黏性流动和湍流大涡模拟程序MVFT[6, 21],对“反尖端”界面不稳定性及湍流混合进行大规模并行数值模拟,在高分辨率的网格下细致地研究“反尖端”界面不稳定性的发展过程和规律,包括波系作用、湍流混合区结构的演化、湍流混合区统计量的变化等。
可压缩多介质黏性流动和湍流大涡模拟程序MVFT所求解的流场控制方程是经过Favre滤波后的可压缩多介质黏性流动Navier-Stokes方程组
{∂ˉρ∂t+∂ˉρ˜uj∂xj=0∂ˉρ˜ui∂t+∂(ˉρ˜ui˜uj+ˉpδij)∂xj=∂σij∂xj−∂τij∂xj∂ˉρ˜E∂t+∂(ˉρ˜uj˜E+ˉp˜uj)∂xj=∂(σij−τij)˜ui∂xj−∂(qlj+QTj)∂xj∂˜Y(s)∂t+˜uj∂˜Y(s)∂xj=∂∂xj(˜D∂˜Y(s)∂xj)−∂QYj∂xj |
(1) |
式中:下标i、j分别代表x、y、z 3个方向,遵循张量运算法则;
σij=μl[∂˜ui∂xj+∂˜uj∂xi−23δij(∂˜uk∂xk)] |
(2) |
式中:μl为流体黏性。气体介质的状态方程采用理想气体状态方程。
具体计算过程为:首先采用算子分裂技术将(1)式描述的物理过程分解为3个子过程进行计算,即整个通量分解为无黏通量、黏性通量和热通量3部分。无黏通量的计算采用多介质高精度PPM(Piecewise Parabolic Method)方法,两步Lagrange-Remapping型的PPM方法分4个步骤进行:(1)物理量分段抛物插值,(2)近似Riemann问题求解,(3)Lagrange方程组推进求解,(4)将物理量映射到静止的欧拉网格上。然后在无黏通量的基础上,采用二阶空间中心差分方法和两步Runge-Kutta时间推进方法求解黏性通量、热通量及标量输运通量。
本研究所采用的模型是英国AWE开展的“反尖端”RM不稳定性激波管实验[27],计算模型如图1所示,由细金属丝制成的网状支架支撑微米量级的硝化纤维薄膜将实验气体隔开,形成预设的初始界面。激波管内截面尺寸为20 cm×10 cm,“反尖端”区气体为SF6,其宽度为15 cm,两侧气体是空气,“反尖端”界面平衡位置距离激波管尾端20 cm。另外,在初始界面上设置振幅小于0.01 cm、波长为0.5~5 cm的小尺度随机扰动,以模拟支架的扰动效应。初始入射冲击波马赫数为1.26。计算域:(x, y, z)∈[0, 80 cm]×[–10 cm, 10 cm]×[0, 10 cm]。实验得到的结果十分有限,只是不同时刻混合区的演化图像及壁面气泡和尖钉位置(距激波管末端的距离),不足以深入分析界面不稳定性及湍流混合发展规律。为此,本研究借助大规模数值模拟,利用1 024个CPU,采用粗(网格大小1.00 mm)、中(网格大小0.50 mm)、细(网格大小0.25 mm)3种不同分辨率的网格,详细研究“反尖端”界面不稳定性及湍流混合的演化特性。细网格下的网格总量为1.024×109。气体初始参数如表1所示,其中γ为气体比热比。
Gas | ρ/(kg·m–3) | p/MPa | γ | μl/(Pa·s) | Diffusion coefficient/(m2·s–1) |
SF6 | 5.97 | 0.1 | 1.09 | 1.474 6×10–5 | 0.97×10-5 |
Air | 1.18 | 0.1 | 1.40 | 1.852 6×10–5 | 2.04×10–5 |
冲击波从低阻抗介质(如空气)向高阻抗介质(如SF6)方向加载物质界面时,会产生一个透射冲击波和一个反射冲击波;冲击波从高阻抗介质向低阻抗介质方向加载物质界面,会产生一个透射冲击波和一个反射稀疏波。图2给出了以密度显示的“反尖端”计算模型轴线上的一维近似波谱,可以看出复杂的波与界面以及波与波之间的相互作用,冲击波加载使流场密度增大,稀疏波卸载使流场密度减小。首先,入射冲击波从左侧加载空气/SF6左界面,产生一个左行反射冲击波和一个右行透射冲击波;该右行透射冲击波加载SF6/空气右界面,在SF6气体中产生一个左行反射稀疏波,向右边空气中透射一个右行冲击波;右边空气中的右行冲击波会从激波管尾端反射回来再次与SF6/空气右界面作用,向SF6中产生一个左行透射冲击波(会在左界面上发生分解,向SF6中反射一个右行稀疏波),并向右边空气中产生一个右行反射冲击波,该右行反射冲击波会反复从激波管尾端反射回来加载SF6/空气右界面,在界面上发生分解;前面第一次在SF6/空气右界面上分解向SF6中产生的左行反射稀疏波和空气/SF6左界面作用后向SF6中产生一个反射压缩波,该压缩波和SF6/空气右界面作用分解产生一个左行反射稀疏波,之后在SF6左右界面之间反复分解成压缩波和稀疏波。如图2所示,在SF6气体区域的胞格结构显示了冲击波、压缩波、稀疏波、界面之间的复杂相互作用,使得界面不稳定性发展更复杂,而波和界面的每一次作用都会促使混合区快速发展[22]。
图3给出了“反尖端”界面及流场演化的实验图像和数值模拟结果比较,左列是实验结果,右3列从左至右依次是网格为1.00、0.50、0.25 mm时数值模拟的沿展向平均的密度场图像,清晰地显示了波和界面的相互作用及波系在流场中的相互作用和演化过程。实验和数值模拟二者之间,包括界面的形状和位置及波阵面的形状和位置均吻合很好。当冲击波与SF6区域左右界面作用后,界面上的小尺度随机扰动开始增长;而且在冲击波与SF6区域的右界面作用时,由于冲击波是从高阻抗介质向低阻抗介质方向加载界面,所以右界面的大尺度“反尖端”界面发生反相,从中心轴处生长出大尺度尖钉,靠近激波管上下壁处形成大尺度气泡结构,并且尖钉和气泡逐渐长大。对比3种不同网格分辨率数值模拟的展向平均密度场可以看出,界面大尺度结构和波系的演化差别很小,只是在中后期界面小尺度结构之间有较小的差别,这是由于展向平均抹去了一些细节结构造成的。图4、图5和图6分别给出了2.0、3.0、4.0 ms时刻不同网格分辨率下以SF6体积分数(YSF6)显示的湍流混合区三维图像。图中显示了复杂的三维湍流混合区演化,包括小尺度气泡/尖钉的长大、变形、融合等;不同分辨率网格计算所得的湍流混合区图像差别很大,其中高分辨率网格捕捉到了更精细的湍流混合区结构。图7给出了大尺度壁面气泡和中心尖钉位置随时间变化的曲线,该位置是相对于激波管尾端的距离,不同颜色曲线表示不同网格分辨率下的数值计算结果。可以看出,在不同网格分辨率下计算所得的大尺度壁面气泡和中心尖钉的位置差别较小,这是因为大尺度壁面气泡和尖钉的位置由湍流混合区的包络线确定,而湍流混合区的包络线主要由大尺度结构决定,湍流混合区大尺度结构的平均特征在不同网格分辨率下的差别很小;虽然差别较小,但是细网格下的模拟结果与实验结果吻合得更好。
图8和图9显示了1.0、2.0、3.0、4.0 ms时中心轴上流场密度和SF6体积分数在不同网格分辨率下的分布。图中清晰地显示了冲击波的传播及其在界面上的分解,以及界面的发展演化,比如1.0 ms时刻的密度分布显示,冲击波加载使流场密度增大,SF6气体中的右行透射冲击波与SF6/空气界面作用后向SF6气体中反射的左行稀疏波的卸载效应使流场密度减小。在前期(1.0~2.0 ms),不同网格分辨率计算的流场密度和SF6体积分数的分布差别很小,说明混合区流场湍流发展不充分,小尺度脉动还未完全发展起来;后期(3.0~4.0 ms)不同网格分辨率的计算结果差别增大,表明在冲击波的多次复杂加载作用下混合区流场已经进入完全湍流阶段,更多小尺度脉动出现,细网格对这些小尺度含能涡结构有更高的分辨率,而粗网格将大量的小尺度含能涡结构耗散掉。
湍动能和拟涡能可以用来表征流场湍流脉动和湍涡的强度和发展。图10和图11分别为1.0、2.0、3.0、4.0 ms时刻无量纲化的流场湍动能K和拟涡能Ω在不同网格分辨率下沿冲击波运动方向的分布,无量纲化速度参数取初始入射冲击波后的流场速度133.6 m/s,无量纲化长度参数取计算模型宽度0.1 m。可以看出:随着混合区的发展,小尺度脉动逐渐发展起来;随着网格分辨率的增大,对这些小尺度脉动的捕捉能力也逐渐提高,捕捉到了更强的湍流脉动和湍涡。所以更小尺寸的网格对流场有更高的分辨率,可以捕捉到更精细的混合区湍流结构,这对于研究RM不稳定性及湍流混合演化的动力学行为非常重要。
利用可压缩多介质黏性流动和湍流大涡模拟程序MVFT,对“反尖端”界面不稳定性及其诱发的湍流混合进行了大规模并行三维数值模拟分析,揭示了界面不稳定性及湍流混合的复杂发展过程和规律。冲击波和空气/SF6/空气界面作用并发生分解而产生冲击波、稀疏波、压缩波,这些次生波在SF6气体中运动并相互作用,而且多次加载物质界面以及波和界面的每一次作用都会加速湍流混合区的发展和物质混合。高/低阻抗构型的“反尖端”界面受冲击加载后发生反相而发展为大尺度的壁面气泡和中心轴尖钉结构,该大尺度气泡和尖钉结构基本确定了湍流混合区的平均几何特征和包络范围而不依赖计算网格。采用高分辨率的计算网格能够捕捉到更精细的混合区小尺度湍涡结构及更强的湍流脉动,显示出湍流混合区的复杂结构和特征。
[1] |
NURICK G N, SHAVE G C. The deformation and tearing of thin square plates subjected to impulsive loads—an experimental study [J]. International Journal of Impact Engineering, 1996, 18(1): 99–116. doi: 10.1016/0734-743X(95)00018-2
|
[2] |
RAMAJEYATHILAGAM K, VENDHAN C P. Deformation and rupture of thin rectangular plates subjected to underwater shock [J]. International Journal of Impact Engineering, 2004, 30(6): 699–719. doi: 10.1016/j.ijimpeng.2003.01.001
|
[3] |
NURICK G N, RADFORD A M. Deformation and tearing of clamped circular plates subjected to localised central blastloads [C]//REDDY B D. Recent developments in computational and applied mechanics. A Volume in Honour of John B Martin, 1997: 276–301.
|
[4] |
RAJENDRAN R, NARASIMHAN K. Damage prediction of clamped circular plates subjected to contact underwater explosion [J]. International Journal of Impact Engineering, 2001, 25(4): 373–386. doi: 10.1016/S0734-743X(00)00051-8
|
[5] |
JACOB N, YUEN S C K, NURICK G N, et al. Scaling aspects of quadrangular plates subjected to localised blast loads—experiments and predictions [J]. International Journal of Impact Engineering, 2004, 30(8/9): 1179–1208.
|
[6] |
LONGÈRE P, GEFFROY-GRÈRE A-G, LEBLÉ B, et al. Ship structure steel plate failure under near-filed air-blast loading: numerical simulations vs experiment [J]. International Journal of Impact Engineering, 2013, 62: 88–98.
|
[7] |
李典, 郑羽, 陈长海, 等. 空爆载荷下舰船典型结构损伤研究进展 [J]. 船舶力学, 2020, 24(4): 543–557. doi: 10.3969/j.issn.1007-7294.2020.04.015
LI D, ZHENG Y, CHEN C H, et al. Review on damage of typical ship protective structures under explosion load [J]. Journal of Ship Mechanics, 2020, 24(4): 543–557. doi: 10.3969/j.issn.1007-7294.2020.04.015
|
[8] |
JACOB N, NURICK G N, LANGDON G S. The effect of stand-off distance on the failure of fully clamped circular mild steel plates subjected to blast loads [J]. Engineering Structure, 2007, 29(10): 2723–2736. doi: 10.1016/j.engstruct.2007.01.021
|
[9] |
TEELING-SMITH R G, NURICK G N. The deformation and tearing of thin circular plates subjected to impulsive loads [J]. International Journal of Impact Engineering, 1991, 11(1): 77–91. doi: 10.1016/0734-743X(91)90032-B
|
[10] |
SHEN W Q, JONES N. Dynamic response and failure of fully clamped circular plates under impulsive loading [J]. International Journal of Impact Engineering, 1993, 13(2): 259–278. doi: 10.1016/0734-743X(93)90096-P
|
[11] |
BONORCHIS D, NURICK G N. The influence of boundary conditions on the loading of rectangular plates subjected to localised blast loading—importance in numerical simulations [J]. International Journal of Impact Engineering, 2009, 36(1): 40–52. doi: 10.1016/j.ijimpeng.2008.03.003
|
[12] |
NURICK G N, GELMAN M E, MARSHELL N S. Tearing of blast loaded plates with clamped boundary conditions [J]. International Journal of Impact Engineering, 1996, 18(7/8): 803–827.
|
[13] |
MENKES S B, OPAT H J. Broken beams—tearing and shear failures in explosively loaded clamped beams [J]. Experimental Mechanics, 1973, 13(11): 480–486. doi: 10.1007/BF02322734
|
[14] |
BAO Y B, WIERZBICKI T. On fracture locus in the equivalent strain and stress triaxiality space [J]. International Journal of Mechanical Sciences, 2004, 46(1): 81–98. doi: 10.1016/j.ijmecsci.2004.02.006
|
[15] |
TENG X, WIERZBICKI T. Evaluation of six fracture models in high velocity perforation [J]. Engineering Fracture Mechanics, 2006, 73(12): 1653–1678. doi: 10.1016/j.engfracmech.2006.01.009
|
[16] |
BAO Y B, WIERZBICKI T. On the cut-off value of negative triaxiality for fracture [J]. Engineering Fracture Mechanics, 2005, 72(7): 1049–1069. doi: 10.1016/j.engfracmech.2004.07.011
|
[1] | YUAN Jianan, LI Jianfu, WANG Xiaoli. Research Progress of High Energy Density Nitrogen[J]. Chinese Journal of High Pressure Physics, 2024, 38(4): 040102. doi: 10.11858/gywlxb.20230797 |
[2] | LI Feng, ZHANG Chenyu, WANG Yue, WANG Bo, ZHANG Mengyu, JING Yadong. Optimization Design of Gas-Liquid Conveying Pipe Structure for 20 L Spherical Explosion Experimental Device[J]. Chinese Journal of High Pressure Physics, 2023, 37(4): 045301. doi: 10.11858/gywlxb.20230651 |
[3] | LIAO Fang, LI Shiqiang, WU Guiying. Topological Optimization and Dynamic Response of Periodic Porous Sandwich Structure under Impact Load[J]. Chinese Journal of High Pressure Physics, 2022, 36(5): 054201. doi: 10.11858/gywlxb.20220560 |
[4] | ZHAO Xiaoming, YANG Yumin, JIANG Nan, CAI Zhongwei, OUYANG Song. Optimization of Charging Structure of Surrounding Holes in Smooth Blasting of Deep Diversion Tunnel[J]. Chinese Journal of High Pressure Physics, 2022, 36(4): 045301. doi: 10.11858/gywlxb.20220503 |
[5] | LI Teng, ZHANG Chenfan, DENG Qingtian, LI Xinbo, WEN Jinpeng. Optimized Design and Energy Absorption of TPU Material Based on Hierarchical Structure[J]. Chinese Journal of High Pressure Physics, 2022, 36(6): 064104. doi: 10.11858/gywlxb.20220542 |
[6] | YAO Chenhui, YANG Gang, ZHANG Zhe, LI Anqi. Reliability Optimization Design of Anti-Penetration Perforated Armor[J]. Chinese Journal of High Pressure Physics, 2022, 36(4): 045102. doi: 10.11858/gywlxb.20220507 |
[7] | LIU Feiming, LEI Jianyin, QIAO Li, LIU Zhifang. Crashworthiness Optimization of Horsetail-Bionic Thin-Walled Structures[J]. Chinese Journal of High Pressure Physics, 2022, 36(5): 054205. doi: 10.11858/gywlxb.20220516 |
[8] | LI Jianping, LIU Siqi. Structure Optimization of Square Explosive Dispersion Device[J]. Chinese Journal of High Pressure Physics, 2020, 34(2): 023301. doi: 10.11858/gywlxb.20190835 |
[9] | QIAO Yu, CHEN Li, ZHANG Qingming, JU Yuanyuan, LU Yangyu. Dynamic Response of the Polyetherimide Sabot of a Long Rod Projectile under Bore Pressure Load[J]. Chinese Journal of High Pressure Physics, 2020, 34(6): 065106. doi: 10.11858/gywlxb.20200558 |
[10] | XIAO Hongcheng, YIN Dongmei, LIN Qinghua, LI Baoming. Structural Optimum Design and Pre-stress Simulation of Fiber Housing for Railgun[J]. Chinese Journal of High Pressure Physics, 2018, 32(5): 055107. doi: 10.11858/gywlxb.20170594 |
[11] | LIU Hongjie, WANG Weili, MIAO Run, WU Shiyong. Optimum Design of Annular Double Done Shaped Charge Structure[J]. Chinese Journal of High Pressure Physics, 2018, 32(6): 065105. doi: 10.11858/gywlxb.20180539 |
[12] | FU Bi-Wei, AI Yu, XI Yan-Qing. Flow Field Analysis and Structure Optimization of the Nozzle Tail Flow[J]. Chinese Journal of High Pressure Physics, 2016, 30(4): 328-334. doi: 10.11858/gywlxb.2016.04.010 |
[13] | REN Guo-Wu, TANG Tie-Gang, GUO Zhao-Liang, LI Qing-Zhong. Statistical Study on Dynamic Fragmentation of Multiple Metal Rings[J]. Chinese Journal of High Pressure Physics, 2015, 29(2): 117-122. doi: 10.11858/gywlxb.2015.02.005 |
[14] | CHEN Chuang, LI Wei-Bing, WANG Xiao-Ming, LI Wen-Bin, WU Wei. Optimization Design of Precursor K-Charge Structure of Tandem Warhead[J]. Chinese Journal of High Pressure Physics, 2014, 28(1): 73-78. doi: 10.11858/gywlxb.2014.01.012 |
[15] | CHEN Jun, LONG Yao, CHEN Dong-Quan. Study on the Thermal Conductivity of HMX Explosive: Method Development and Theoretical Calculation[J]. Chinese Journal of High Pressure Physics, 2013, 27(2): 199-204. doi: 10.11858/gywlxb.2013.02.005 |
[16] | GAO Xiang, CHEN Peng-Wan, LIU Jian-Jun, XIONG Guang-Yuan. Experimental Study on Phase Transformation of TiO2 Induced by High Energy Milling[J]. Chinese Journal of High Pressure Physics, 2012, 26(4): 421-425. doi: 10.11858/gywlxb.2012.04.010 |
[17] | SUN Zhan-Feng, XU Hui, LI Qing-Zhong, ZHANG Chong-Yu. Further Study on JWL Equation of State of Detonation Product for Insensitive High Explosive[J]. Chinese Journal of High Pressure Physics, 2010, 24(1): 55-60 . doi: 10.11858/gywlxb.2010.01.010 |
[18] | GAO Ling-Ling, MA Yan-Mei, LIU Dan, HAO Jian, JIN Yun-Xia, WANG Feng, WANG Qiu-Shi, ZOU Guang-Tian, CUI Qi-Liang. Raman Spectra Characterization of Cycloheptane under High Pressure[J]. Chinese Journal of High Pressure Physics, 2008, 22(2): 192-196 . doi: 10.11858/gywlxb.2008.02.013 |
[19] | HU Dong, YUAN Chang-Ying, LI Ping, CHENG Xin-Lu, LIU Jin-Chao, SUN Zhu-Mei. Spectroscopic Studies of Epoxypropane Ignition[J]. Chinese Journal of High Pressure Physics, 2003, 17(3): 169-172 . doi: 10.11858/gywlxb.2003.03.002 |
[20] | LI Peng-Yuan, GENG Man, TIE Jun, HUANG Pei-Ji, CUI Xi-Rong, YAN Bing, WANG Jing-Quan, JIANG Yuan-Qi. A Study of the Heat Resistance of N+ Implanted Synthetic Diamond[J]. Chinese Journal of High Pressure Physics, 1997, 11(4): 303-306 . doi: 10.11858/gywlxb.1997.04.012 |
Gas | ρ/(kg·m–3) | p/MPa | γ | μl/(Pa·s) | Diffusion coefficient/(m2·s–1) |
SF6 | 5.97 | 0.1 | 1.09 | 1.474 6×10–5 | 0.97×10-5 |
Air | 1.18 | 0.1 | 1.40 | 1.852 6×10–5 | 2.04×10–5 |