爆炸驱动液体分层现象和液滴尺寸模拟

爆炸驱动液体介质外界面的分散和破碎是气溶胶云团形成的重要过程。采用基于维数分裂的欧拉程序和Youngs混合界面处理方法,对中心药爆炸驱动甘油和水介质流场的液体分层现象进行了数值模拟。结合试验结果推断提出了液滴形成过程的三种并存机制:外层射流破碎、内层R-T失稳和中间液层"空化"破碎,分别建立了不同液层破碎液滴的尺寸模拟方法。对比给出抛撒甘油和水装置初级液滴的尺寸分布及最外层理论射流量。     

欧拉计算程序实现非均质炸药冲击起爆的数值模拟

爆轰的数值模拟研究多采用Lagrangian方法，这是因为Lagrangian网格随质点运动，介质之间的界面用滑移面来描述，故能精确地计算介质的界面。但在遇到网格扭曲和网格变形时，计算往往无法继续进行，这时就要用到网格重分等技术，这无疑增加了人为干预和计算误差。Eulerian方法能计算大变形问题，但不能精确地描述物质的界面，对于模拟平面撞击实验，爆轰波在单种介质中传播，对物质界面几乎没有要求。所以，Eulerian方法在模拟此类问题较Lagrangian方法具有明显的优势。     

耦合AMR的GEL算法

数值模拟中，欧拉方法能计算大变形流场，但不能精确地区分物质界面，拉氏方法的单元边界即为物质边界，因此可以精确区分不同的物质，但当计算单元变形较大时计算精度变差甚至无法进行。如果在流场内不同区域采用不同的计算方法，在计算区域交界处进行合理的数据交换，则既能计算大变形流场又能在流场内保持清晰的物质界面。     

冲击压力及加载速率对沟槽微射流的影响

通过二维弹塑性流体力学欧拉程序MEPH,系统研究冲击压力及加载速率对沟槽微射流喷射量的影响.计算结果显示,微射流的喷射量随冲击压力的增大而增大,但喷射量对压力的变化并不敏感;计算结果和实验结果均表明,加载速率对沟槽微射流的喷射量有较大影响,随着加载速率的减小,微粒子喷射量大致按指数规律减少,所得数值计算结果得到实验验证.     

冲击作用下非均质炸药热点形成的离散元方法

用离散元方法模拟以HMX为基的塑料粘结炸药与HMX颗粒炸药在冲击载荷作用下的细观响应过程.结果表明对于HMX颗粒炸药,炸药颗粒界面处的温度及压力远高于颗粒内部;对于HMX为基的塑料粘结炸药,粘结剂处的温度及压力亦高于颗粒内部;两种冲击响应过程的对比分析表明,粘结剂有效地降低了炸药颗粒边界处的温度及压力.     

冲击加载下孔洞形成微射流的最大侵彻深度

应用二维弹塑性流体力学欧拉程序MEPH,研究冲击载荷作用下孔洞形成的微射流.主要关注孔洞微射流的形成过程,重点分析冲击压力及孔洞半径等因素对于孔洞微射流最大侵彻深度的影响.分析孔洞塌缩形成射流的空间速度分布特征,并应用非定常射流理论,对不同冲击压力及孔洞半径条件下,微射流的最大侵彻深度进行估算,估算结果与数值模拟结果符合较好.     

计算域可变的CEL方法

为了更好地兼顾拉氏方法和欧拉方法各自的特长,提出一种可将拉氏介质映射到欧拉计算域的耦合欧拉-拉格朗日(CEL)方法。通过这种映射,将欧拉-拉格朗日重叠区域的接触面协调问题转换为欧拉区域内的多介质计算问题,简化了CEL方法的构造过程。通过与侵彻实验和结构对爆炸冲击波响应实验的比较,验证了新算法,计算结果与实验数据符合较好。     

含氦泡金属铝层裂响应的数值分析

因自辐照效应的影响,一些材料内部会产生大量的氦泡,关注这些氦泡对材料力学性能的影响是目前损伤破坏研究中的重要问题之一。结合相关文献的实验结果,采用耦合材料初始损伤、孔洞尺寸及惯性影响的损伤模型,对该问题进行了数值分析。结果显示:氦泡的内压及材料变形中温度的变化对损伤发展的影响很小;材料的初始损伤越大,材料内部应力减小得越快,损伤增长得越慢;因惯性的影响,初始氦泡越大,损伤增长相对较慢。因此,分析含氦泡材料的层裂损伤问题需要重点关注材料初始氦泡大小、初始损伤以及损伤演化过程中惯性的影响。     

柱面内爆驱动金属界面不稳定性的数值模拟研究

对柱面爆轰驱动内壁刻有正弦扰动的金属钢壳与内部硅橡胶界面产生不稳定性问题进行数值模拟,计算结果与实验结果定性符合.与不考虑金属强度情况对比分析认为,未熔化状态下金属强度对不稳定性具有较强抑制作用,在某些加载条件下会使扰动增长率随扰动模数增加而减小.之后,对强度因素影响下内爆压缩驱动金属不稳定性问题的扰动发展规律进行了总结.在聚心反射波到达壳体之前,造成初始界面反转的RM不稳定性起主导作用,随着扰动模数增加扰动由呈近似线性发展到基本不发展变化,基本不变化后的扰动振幅也随模数增加而减小.聚心反射波作用到壳体内界面后,减速RT不稳定性作用明显增强,与强度等因素共同作用造成扰动呈明显非线性发展.无论是前期RM不稳定性主导阶段还是之后以减速RT不稳定性为主的扰动发展阶段,强度因素均能造成未熔化状态下金属不稳定性截止波长存在.