爆炸冲击波强间断问题的高阶伪弧长算法研究

来源 :爆炸与冲击 | 被引量 : 0次 | 上传用户:morningwind2009
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
为了提高对冲击波强间断处的分辨率,通过引入弧长参数,使网格自适应地朝着间断处移动,并结合高精度WENO数值格式,进而达到了对大梯度物理量的高分辨率捕捉。针对网格移动造成的非均匀和非正交现象,通过坐标变换,使得计算过程在均匀正交的计算空间中进行。通过和有限体积下的数值结果对比,结合数值误差分析,可以看到高阶伪弧长数值算法不仅保证了高精度而且对间断的捕捉更加明显,在间断附近解的整体光滑性较好,网格的自适应移动使得解的奇异性得到了削弱,因此可以削弱高阶格式容易引起数值振荡这个缺点。最后采用高阶伪弧长算法计算了化学反应流问题,结果表明高阶伪弧长算法有着较快的收敛率,对于解决爆炸与冲击强间断问题有着较为明显的优势。“,”In this paper, in order to improve the resolution of capturing discontinuities, we introduce the pseudo arc-length parameter to make the mesh move to the discontinuities adaptively. By combining the high-precision WENO scheme with the pseudo arc-length algorithm, the advantages of both schemes can be shown, on the one hand, the solution has a higher convergence rate, on the other hand, it has a higher resolution for the region solution with larger physical variation . Because the traditional high-order scheme is based on Cartesian grid, and the grid in the pseudo arc-length numerical calculation is deformed. In view of the non-uniform grid and non-orthogonal deformation grid caused by the grid moving, the original deformed physical space is mapped to the uniform orthogonal arc-length calculation space by introducing coordinate transformation, and then the classical higher order scheme is used to solve the governing equations in the computational coordinate system. Through the comparison of some numerical examples and the analysis of numerical errors, it can be found that the pseudo arc-length algorithm is better than the finite volume method with fixed mesh. The high-order pseudo arc-length algorithm has a very high resolution to capture discontinuities, and the density of the grid near the discontinuities is very high. The adaptive grid movement weakens the singularity of the governing equation near the discontinuity, so the whole solution is smooth and the numerical oscillation is not obvious. This shows that the pseudo arc-length algorithm can overcome the shortcomings of high-order schemes which easily cause numerical oscillations. Finally, the chemical reaction flow problem is calculated. The results show that the high-order pseudo arc-length numerical algorithm with less mesh number has faster convergence rate and higher discontinuous resolution. Therefore, the high-order pseudo arc-length algorithm has obvious advantages in dealing with the strong discontinuity problem of explosion and shock.
其他文献
针对爆炸加载下金属柱壳膨胀断裂破片软回收的研究需求,本文通过理论分析和初步的数值模拟设计了由低密度聚氨酯泡沫与水介质为主体的回收装置.与传统单一材料为主的回收装置相比,该回收装置既能在破片高速阶段将低阻抗聚氨酯泡沫对破片的冲击压力减小到约为水对破片冲击压力的1/3,又使破片速度全程持续地较大幅度衰减,还能在破片低速阶段又能充分利用水介质密度大的优势,减小以聚氨酯泡沫单一材料为主的回收装置尺寸.依托该装置开展了炸药加载下304不锈钢柱壳膨胀断裂回收实验.通过测量回收池外壁速度、检查实验后的回收池外观,发现回
针对复合圆柱壳在炸药爆轰作用下的动力学响应及在此过程中伴随的失稳问题,研究了其制造工艺中可能出现的缺陷以及圆柱壳中铜线螺旋角和直径对复合圆柱壳稳定性产生的影响.采用SPH-FEM耦合算法,建立了复合圆柱壳二维细节模型,并提出了一种基于圆柱壳内壁粒子速度历史的失稳判据,计算了在不同参数条件下复合圆柱壳的失稳时间及对应的压缩率,对影响复合圆柱壳稳定性的因素进行了评估.分析结果表明,在复合圆柱壳制备过程中存在的折返层缺陷和铜线直径对复合圆柱壳的稳定性有较大影响,而螺旋角度对其稳定性影响不大.
探究了菌渣的水热液化转换成生物油燃料的过程.结果 表明,抗生素菌渣在260℃、保留时间是135 min时,获得最大的生物油产率(28.01%).通过6种不同的催化剂进行催化,加入催化剂后,生物油产率最大的是Na2CO3 (36.06%)和NaOH(36.31%).碱催化的生物油的含氮化合物的质量分数在41.16% ~ 49.74%之间,而酸催化产生的生物油含氮化合物的量在57.62%~ 59.32%之间.通过调节催化剂Na2CO3、NaOH的添加量发现,在投加量为8%时,生物油含氮量均最低,Na2CO3和
针对弹性楔形体砰击载荷及结构响应问题,采用理论解析计算方法以及任意拉格朗日-欧拉(arbitrary Lagrangian-Eulerian,ALE)流固耦合数值计算方法开展研究.分别分析不同边界条件、入水速度、板厚以及斜升角对弹性结构所受砰击力及结构响应的影响及变化规律,并探讨砰击载荷时间历程和物面分布特性.结果 表明,通过增大斜升角可有效降低弹性楔形体砰击载荷和结构响应,斜升角自10°增大至30°,无量纲砰击力峰值减小至6.9%,结构变形极值减小至6.5%.可通过水弹性因数RF=CBtanβ√EI/(
为探讨轻型陶瓷复合装甲抗侵彻过程中陶瓷的碎裂行为,采用12.7 mm穿燃弹对不同背板厚度及陶瓷厚度下陶瓷/金属复合装甲进行弹道冲击试验.通过观测回收的靶体陶瓷宏观破坏特征,分析不同厚度组合与陶瓷主要破坏特征之间的关系;并通过对陶瓷碎块的多级筛分称重,分析不同厚度组合下陶瓷面板的碎块尺度分布规律.结果 表明,陶瓷锥是陶瓷面板的主要破坏形态,其宏观裂纹主要有:径向裂纹、环向裂纹和锥形裂纹.陶瓷锥内可细分为由高压缩应力引起的粉末状较小陶瓷碎块组成的陶瓷粉碎锥和由应力波造成的较大片状陶瓷碎块组成的陶瓷破碎锥.冲击
为了研究中心点火管火焰在药床中的传播规律,设计了可视化模拟试验平台,开展了不同点火药量、不同装药结构的中心点传火试验。采用高速图像采集系统记录了中心点火管火焰在药床中的传播过程,采用瞬态压力记录仪记录膛内压力的时空变化。结果表明,点火药量为20 g时,出火时间为0.6 ms;点火药量为30 g时,出火时间为1.5 ms;杆状装药床的传火时间平均为2.2 ms,粒状装药床的传火时间平均为3.4 ms,而杆粒混装药床的传火时间为3.1 ms。可见,点火药量对药床出火时间影响显著,较大的点火药量导致药床出火时间
高熵合金,以其独特的合金设计和优异的综合性能,成为当下材料研究的热点。本文利用高真空电弧熔炼法成功制备出了CoCrFeNiAlx(x=0, 0.6, 1)系高熵合金,并通过分离式霍普金森压杆对其进行一系列不同应变速率下的动态压缩试验。通过X射线、扫描电镜和透射电镜分析,深入探索了该合金系的晶体结构、微观组织和变形特征。最后,利用修正后的Johnson-Cook (J-C)本构模型,获得了该体系高熵合金的动态本构关系。“,”High-entropy alloys (HEAs), due to their u
火灾与爆炸通常相伴发生,对工程结构安全造成了严重威胁。为研究高温下钢管混凝土柱抗爆性能,采用ABAQUS有限元软件建立了ISO 834标准火灾作用下钢管混凝土柱抗爆模型。在验证有限元模型可靠性基础上,首先分析了标准火灾作用下钢管混凝土柱抗爆工作机理;其次重点研究了受火时间、材料强度、含钢率以及爆炸当量对构件在标准火灾下抗爆性能的影响。研究结果表明:火灾作用下两端固结的钢管混凝土柱受爆炸荷载时,柱两端首先发生剪切破坏,随后整体发生受弯破坏;随着受火时间增加,钢管耗能占比降低,混凝土塑性变形逐渐成为主要耗能机
为了研究当量比对汽油燃料两相旋转爆轰发动机工作特性的影响,开展了以高总温空气为氧化剂的气液两相旋转爆轰实验研究。旋转爆轰发动机环形燃烧室外径、内径和长度分别为202、166和155 mm。汽油和高温空气采用高压雾化喷嘴与环缝对撞喷注的方式混合,以此提高推进剂的掺混效果与活性,采用预爆轰管作为点火装置。实验通过改变汽油质量流量改变推进剂当量比,并基于燃烧室内测得的高频动态压力和平均静压,对气液两相旋转爆轰波的传播模态和传播特性以及发动机的工作特性进行了详细分析。实验结果表明:在当量比为0.79~1.25时,
利用分离式霍普金森压杆(SHPB),对粗砂岩、中等粒径砂岩和细砂岩进行了应变率为69~83 s–1的动态单轴抗压实验,研究了粒径尺寸效应对砂岩动力特性的影响。通过三维数字图像相关(3D-DIC)技术分析高速摄像图像,获得了砂岩的实时应变场,据此分析了动态荷载下3种粒径砂岩的动力变形特性和裂纹开展行为。结果表明,砂岩弹性应变储能可逆释放的临界应变率随着粒径的减小而升高,动态压缩强度随着粒径减小而增大,动态强度应变率敏感度则与强度规律相反。相较于静态条件下,中等粒径砂岩和细砂岩的动态弹性模量增长了2~3倍,粗