不同直径圆球诱导燃烧的振荡机制与频率特性

为了分析不同直径圆球诱导振荡燃烧的规律，并揭示圆球大小在振荡燃烧现象中所发挥的深层次作用，本文采用二维轴对称Euler方程和基元反应模型，对不同直径的圆球在H2/air预混气体中诱导振荡燃烧的现象开展数值模拟研究。研究发现，振荡频率并不是简单地随直径增大而逐渐从高频向低频连续过渡，而是存在两次突变，形成了超高频、高频以及低频三种振荡燃烧模态。在两种模态间过渡时，振荡达到稳定状态前，会存在一段双频耦合的振荡阶段。三种不同振荡燃烧模态的产生是受到了不同振荡机制的作用，而两种模态间过渡时的双频耦合现象则是两种机制相互竞争的结果。     

可压缩多介质流动问题的高精度数值模拟方法

在可压缩流动问题的数值模拟领域,激波的高分辨率计算已取得重要进展。但是包含物质界面的可压缩多介质流动的数值模拟还存在诸多数值挑战,主要表现为界面处数值耗散过大和非物理振荡等问题。界面处流体性质的不连续性是造成可压缩多介质流动问题物理建模与数值方法困难的主要原因。为了建立一套高效的可压缩多介质流动问题的高精度数值模拟方案,本文从数值框架的选择、非守恒方程相容离散、高精度有界格式构造、界面压缩、界面-单介质分区计算方法等多个维度出发,综述近几年我们在可压缩多介质流动问题高精度数值模拟方法方面的研究进展。通过上述多个维度的工作,我们建立了一套适用于可压缩多介质流动问题的低耗散、基本无数值振荡的高精度欧拉数值方法,并成功应用于可压缩多介质大变形流动和界面不稳定性诱导湍流混合等问题的数值模拟。相关数值方法研究成果已集成至武器物理等领域工程数值模拟软件中,为相关工程任务提供了重要技术支撑。     

基于XFEM-CZM耦合法的胶接接头裂纹扩展分析及强度预测

采用扩展有限元法（XFEM）和内聚力模型（CZM）相耦合的方法，分析胶接接头胶层内部裂纹扩展、界面脱粘分层现象。采用内聚力单元和内聚力接触描述胶层/板材界面，建立单/双搭接接头有限元模型。预测拉伸载荷下接头的强度性能并与已有试验数据对比分析，验证XFEM-CZM耦合法的可行性及内聚力单元和内聚力接触2种界面建模方法的有效性。模拟裂纹从胶层内部扩展至胶层/板材界面并引起界面脱粘分层的过程，分析其损伤失效机理。讨论初始裂纹长度和界面刚度、强度及应变能释放率对胶接接头强度性能的影响。结果表明:胶接接头强度随初始裂纹长度增加而降低，且在双搭接接头模型中表现更为明显；界面刚度、强度对胶接接头强度影响较大而应变能释放率的影响较小。      

EPDM绝热层组分对硫化胶中残留硬脂酸含量的影响

为降低EPDM绝热层中残留硬脂酸含量,分析了EPDM绝热层中硬脂酸与氧化锌的原材料性能和反应活性,当绝热层中无其他组分存在时,两者在橡胶高温硫化过程中可反应完全生成硬脂酸锌,此时硫化胶中硬脂酸基本不发生残留。通过液相色谱等手段测定绝热层硫化胶片中残留硬脂酸含量,研究了EPDM绝热层组分对硬脂酸与氧化锌反应的影响程度和影响机理,发现在一些常用的工艺助剂,如补强剂、阻燃剂和树脂等影响下,氧化锌与硬脂酸会发生可逆化学反应,导致绝热层硫化胶片中游离硬脂酸的含量增加,影响程度与工艺助剂本身特性和用量有关。     

铝钛异种合金电子束焊接接头组织与界面分析

采用真空电子束技术获得了7075铝合金和TC4钛合金异种接头。分别通过金相显微镜、扫描电镜和能谱分析仪等检测仪器观察了焊接接头的宏观形貌、微观组织特征和元素分布情况,同时分析了过渡层中金属间化合物的形成过程及其影响因素。试验结果表明:熔合区铝合金侧微观组织以柱状晶、等轴晶组织为主,还有少量的细晶组织,钛合金侧会析出呈针状或片层状的α相和α'相。由于Ti Al3的自由能最小,在铝焊缝和钛焊缝的界面生成的金属间化合物主要成分为TiAl3,元素扩散不充分导致少量的Ti Al、TiAl2等中间相生成,且随着热输入量的增加,金属间化合物层逐渐变厚,但化合物种类不变。试验结果为铝钛异种合金焊接界面研究提供了借鉴参考。     

圆弧激波缓聚点火及后续燃烧传播

采用激波管实验和准一维数值模拟的方法，对预混可燃气体中圆弧汇聚激波的自点火现象及后续燃烧波的传播特性进行研究。其中圆弧汇聚激波由平面运动激波通过精确设计的弧形过渡管段转变得到。研究表明:收缩段中圆弧汇聚激波波后的非均匀梯度环境由激波在平直段、弧形过渡段和扇形收缩段中传播所分别诱导的3个梯度区共同构成。随着圆弧汇聚激波的不断增强，圆弧激波后某处首先形成一个无激波的温和反应区。该反应区逆流锋面的初期运动速度远超Chapman-Jouguet（CJ）爆轰波速，而反应产物区流动则呈现出一定弱爆轰波特征。进一步分析发现，该反应锋面本质上是一种"自发反应波"（spontaneous reaction wave），而非常规意义上的动力学波，其速度与汇聚激波波后气流点火时间梯度的倒数吻合。而后，反应区的扩张速度很快降至CJ爆轰波速以下，伴随反应锋面附近激波的产生以及激波-火焰复合结构的形成。激波-火焰结构最终加速演变为反向传播的爆轰波。在一定的条件下，由于入射激波转变过程和汇聚所构造的特定点火环境，自发反应波可再次赶超爆轰波，成为新的燃烧波前；而当自发反应波速度再次低于CJ爆轰波速时，它将再次转变为爆轰波；在此过程中，原先的爆轰波阵面蜕变为反应产物中传播的激波。     

热/力耦合作用下基于应力分析的冰破坏准则

根据飞机热除冰的物理过程,考虑外部空气动力和蒙皮表面加热的作用,建立了NACA 0012翼型前缘冰层应力计算模型。采用有限元方法和平面三角形单元对控制方程组进行了求解,获得了外部空气动力和蒙皮表面加热对冰层黏附界面应力的影响规律。研究表明:蒙皮表面不加热时,来流速度影响了黏附界面应力的强度,来流攻角影响了黏附界面应力的分布,冰-蒙皮间黏附界面切应力最大值随来流速度呈近似线性增大趋势,但外部空气动力很难造成冰层破坏。蒙皮表面加热时,冰-蒙皮间黏附界面的耦合应力和冰层内部的主应力随着热流密度的增大而增大,很容易超过剪切强度,这是造成冰破坏的关键因素。耦合冰-蒙皮剪切强度随界面温度的变化关系,初步建立了基于应力分析和热/力耦合作用的冰破坏判断准则。外部空气动力产生的界面应力和蒙皮表面加热产生的界面热应力之和，必须大于与蒙皮表面温度相关的剪切强度，则冰层发生破坏，破坏位置是耦合应力超过剪切强度的区域。      

基于十二烷基苯磺酸固化季戊四醇丙烯醛树脂的粘结体系流变特征及反应动力学研究

为实现季戊四醇丙烯醛树脂(PEAR)/十二烷基苯磺酸(DBSA)体系在浇注PBX炸药中的应用以及获得该体系在工程应用中的工艺温度参数,采用粘度实验研究了体系的粘度特性,采用动态差示扫描量热法(DSC),通过模拟n级反应动力学模型、Kissinger微分法、Ozawa积分法以及Crane方程研究了体系的固化反应动力学。结果表明,50℃以上PEAR粘度几乎不受转速影响,PEAR与DBSA质量比大于25∶1,可保证浇注过程的顺利进行。PEAR/DBSA体系的凝胶化温度为345.92 K,固化温度为383.83 K,后处理温度为411.46 K。PEAR/DBSA体系固化反应为放热反应,反应的表观活化能为74.84 kJ/mol,指前因子为2.54×109min~(-1),反应级数为1.02,反应热为190.66 J/g。     

热障涂层的界面形貌对TGO层生长行为的作用机制

通过对陶瓷层/粘接层界面应力和热生长氧化物(thermally grown oxides,TGO)层的生长规律分析,研究涂层界面凸起对TGO生长行为的作用机制。结果表明:TGO在陶瓷层/粘接层界面凸起区域的生长速率高于其他区域;由ANSYS有限元软件应力分析可知,随着界面粗糙度的增加(10μm增至20μm),涂层的界面应力增加(185 MPa增到406 MPa);TGO层的厚度增加(1.6μm增至9.3μm)同样会使界面应力增加(142 MPa增大到574 MPa);此外,在高温氧化过程中,陶瓷层/粘接层界面凸起区域主要表现为拉应力,凹陷区域为压应力;拉应力能够促进TGO层的快速生长,压应力则会抑制TGO的生长速率;在保证涂层有效结合强度的前提下,降低粘接层的表面粗糙度,能够有利于降低涂层的界面应力以及减缓TGO的生长速率,从而提高热障涂层高温稳定性。