碳酚醛-铝板中二维X射线热击波数值模拟

文章利用平面应变正交各向异性材料率相关性的弹塑性本构模型和各向同性材料理想弹塑性本构模型,采用有限元方法编写程序,对X射线辐照碳酚醛-铝板时诱导的热击波进行数值模拟,讨论了热击波的传播规律。结果表明,在1 keV的 X射线辐照下,X射线穿透极浅,材料表面发生气化现象,应力峰值较大;在3 keV的X射线辐照下,X射线穿透较深,热击波呈现双波结构,材料的层裂更明显。     

等离子体鞘层加速模拟质子辐照连续谱 数值仿真研究

质子辐照是导致空间飞行器热控涂层性能衰退的重要原因。目前地面多采用单一能量的粒子替代空间能量连续分布的粒子来开展质子辐照模拟试验。文章提出了在一个脉冲宽度内获得连续能量质子谱的方法,即脉冲偏压等离子体鞘层加速方法。文章利用细胞粒子（particle-in-cell,PIC）模型对在3种脉冲偏压三角波作用下等离子体鞘层加速质子以获得连续能量质子谱的动力学过程进行了数值仿真研究,分析了3种脉冲三角波形对鞘层扩展、离子加速及能量分布的影响。结果表明,电势空间分布和离子运动状态紧密联系,鞘层内离子密度变化受到脉冲波形和离子热扩散运动的综合影响,通过调整脉冲偏压三角波形,能够获得不同的能量-剂量分布,从而为下一步工作打下基础。     

高能固体推进剂燃烧转爆轰的数值模拟

利用一维两相反应流模型,建立了NEPE高能推进剂在颗粒床中燃烧转爆轰的控制方程和辅助方程,用M ac-Corm ack差分格式进行数值求解,并与实验值进行了比较。结果表明,数值预测与试验结果有较好的一致性,在DDT的各个阶段,颗粒床都存在不同程度的动态压缩,压缩波和燃烧波的相互作用是NEPE推进剂燃烧向爆轰转变的内在原因。     

火星探测器表面材料催化特性数值模拟研究

针对火星探测器高超声速进入过程中的表面材料催化作用及其对气动热环境影响这一问题,建立了变壁面温度的火星大气表面材料催化作用模型,并基于火星大气物理化学模型和求解三维热化学非平衡N-S方程的数值方法,对典型火星探测器防热大底进行了数值模拟,获得了不同催化特性下的高超声速非平衡流场和气动热数据,分析了表面材料催化特性对气动热环境影响的规律性。研究结果表明：表面催化特性对壁面附近组分分布影响很大,催化反应进程主要受O原子浓度限制;有限催化热流随催化效率增大而增大,完全催化峰值热流比催化效率为1的有限催化峰值热流高25%~64%;表面温度随催化特性的变化规律与热流变化规律类似。有限催化模型能根据表面材料的催化特性精细化预测表面热流和温度,为防热设计提供更精确合理的参考标准。     

50 kW级高功率霍尔推力器放电通道数值模拟研究

高功率霍尔推力器兼顾了比冲高、推力大、寿命长等特点。为了提高设计效率并考察热负荷问题,以50 kW级霍尔推力器为对象,采用单元粒子法(PIC)/蒙特卡罗碰撞模型(MCC)/直接模拟蒙特卡罗碰撞模型(DSMC)混合算法,建立二维对称计算模型。基于准电中性假设、中性原子考虑为背景气体,计算得到标准工况下(功率50 kW,流量86.4 mg/s)推力为2.2 N,比冲为2 598 s,与同类推力器实验结果对比,误差分别为5.18%和3.35%。在此基础上,考察了多种工况下(工作电压400~600 V,工质流量69.12~103.68 mg/s)放电通道内离子数密度、离子轴向运动速度、电子温度分布等参数。结果表明:增大工作电压会增强粒子间相互作用及离子加速喷出效果,流量调节影响电子温度和离子数密度分布;从推力器性能方面来看,增大工作电压,推力比冲随之增大,流量增大、推力增大,推力器的热损失功率占比达到15.94%。研究结果为高功率霍尔推力器的设计和实验提供了一定的参考依据。     

SRAM的SEU效应及加固技术三维数值模拟

利用TCAD软件,构建了特征尺寸为0.18μm的三维器件模型,采用器件与电路联合仿真的方法,对重离子在静态随机存储器(SRAM)中引起的单粒子翻转(SEU)效应进行了模拟,分析了SRAM单粒子翻转的机理。仿真了各种阻值的反馈电阻对SRAM抗SEU的效果,确定了SRAM单元抗SEU反馈电阻的阻值。仿真结果表明,搭建的仿真平台可为加固型SRAM电路的研制提供仿真平台和设计依据。     

异形构件化学气相沉积SiC涂层的数值模拟

根据化学气相沉积(CVD)工艺制备SiC陶瓷涂层的工艺特点和典型异形构件的结构特点,建立了异形构件表面化学气相沉积SiC涂层的数学模型。利用该模型,对CVD法在典型异形表面制备SiC涂层进行了数值计算和分析。计算结果显示,带有斜面的构件对CVD SiC沉积过程有显著影响,在反应器大小允许的情况下,构件放置时,斜面与水平面的夹角越小越好,并尽可能将构件长的一面与水平面平行,这样有利于沉积的涂层均匀。此外,对于带有台阶的构件来说,正放的构件表面浓度大于倒放的构件,而浓度梯度则小于倒放的构件。因此,在实际应用中,应尽量使台阶部分放在气流的下游。上述研究结果对CVD工艺制备SiC涂层的优化具有一定的指导意义。     

钛合金管件接头热成形工艺参数的数值模拟

采用PAMSTAMP数值模拟技术对TC4管件接头热成形的过程进行模拟,分析了成形温度、零件结构对成形性能的影响,成形温度越高,零件内圆角越大,零件型面质量越好。通过模拟得出较优的工艺参数,改善了零件的成形效果。     

气动加热下金属蜂窝板热响应特性数值模拟

提出一种计算非稳态气动加热下金属蜂窝板内辐射导热耦合换热数值计算方法。将金 属蜂窝夹芯板作为三层平板结构,建立沿厚度方向一维传热控制方程,采用Monte
Carlo方法求解蜂窝腔内辐射换热,有限体积法求解蜂窝板内耦合换热,考虑蜂窝金属与空 气热物性随温度变化,对某结构的正六角形蜂窝板在模拟气动加热条件下两表面温度响应进 行了数值预测,数值预测结果与文献中的实验结果对比表明该方法能够准确模拟蜂窝夹心板 内传热过程;对采用Swann & Pittman当量导热系数经验公式进行蜂窝非稳态表面温度响应 计算可靠性进行了验证,结果表明该公式在计算非稳态气动加热下蜂窝表面温度时存在较大 误差,给出的热面温度过高,冷面温度过低,两表面温差最大相对误差高达140％。     

强激光对靶材烧蚀效应的数值模拟研究

强激光烧蚀是一个复杂的物理化学过程,包括质量迁移、相变、运动边界等诸多复杂因素。如果靶为复合材料,烧蚀机理更加复杂,而且很难找到共性规律,这对数值计算方法提出了较高的要求。文章介绍了强激光烧蚀的理论基础、烧蚀模型以及对强激光烧蚀进行数值模拟的无网格方法（光滑粒子动力学方法,即SPH方法）。采用了SPH方法对典型靶的激光烧蚀进行了数值模拟,并给出了单层铝靶和环氧树脂/铝双层靶在激光辐照下的烧蚀形貌与温度场。计算结果表明,在不同材料的交界面,烧蚀将沿径向发展。     

液体推进剂爆炸冲击波的数值模拟

针对现有液体推进剂爆炸实验危险性高、耗费大以及实验条件严格及重复性差等缺点及难点,利用Autodyn软件分别对一维和二维两种情况下液体推进剂爆炸产生的冲击波进行数值模拟,并与经验预测及相关实验研究进行比较分析。研究结果表明:在距离大于10 m时,利用软件模拟计算的结果与经验公式一致性较好,两者偏差不超过5%。数值模拟不仅可以准确捕捉冲击波峰值压力,而且还能够较准确地得到液体推进剂在空气中爆炸时冲击波的衰减过程,为液体推进剂爆炸事故的毁伤评估和安全防护奠定了基础。     

高加速度冲击下固体推进剂药柱轴向形变的数值模拟

研究了截面形状、长径比及模量对固体推进剂药柱受到高加速度冲击时发生形变的影响。数值模拟表明，当药柱为厚壁柱壳、星形及圆柱形时，药柱截面形状对发生小变形药柱的应变影响不显著，三种形状药柱的轴向形变、应力、应变值差别不大。而轴向形变与长径比成正比关系，并与模量成反比关系。     

第二边界条件下梯度功能材料一维温度场的数值模拟

采用数值模拟的方法，对梯度功能材料在第二边界条件下温度场进行了分析，并根据所建立的数值方法对新研制的C／C／Al2O3炭陶梯度功能材料和C／C材料进行了比较，对具有不同组分分布曲线的梯度功能材料进行了比较。结果表明，C／C／Al2O3炭陶梯度功能材料的绝热性能优于C／C材料，可通过改变C／C／Al2O3炭陶梯度功能材料的组分分布，优化其内部的温度场分布。     

大气环境电子束材料表面强化机理数值模拟

建立了大气环境电子束材料表面强化过程的多阶段数值模型,用于对电子束与材料表面层的相互作用过程和结构与相变过程进行研究。通过应用这一模型,使得对材料表面强化中电子束传输、温度演化和相变动力学等多参数研究成为可能,由此建立了电子束输入参数与结构相变结果之间的关系。     

非结构网格技术应用于固体火箭发动机内流场数值模拟

从二维欧拉方程出发，采用非结构网格技术对固体火箭发动机内流场进行数值模拟，研究了点火初期、燃烧过程中和燃烧结束时固体火箭发动机燃烧室不规则物理区域内形成的复杂流动，分析了流场结构及特性，并对在这一领域应用非结构网格技术存在的问题进行分析。     

固冲发动机沉积数值模拟与试验研究

分析了液滴撞壁模型的发展概况,基于Mundo模型,分析了粘附、反弹、飞溅3种状态的固冲发动机沉积模型,并对发动机内沉积问题进行了数值模拟与试验研究.结果表明,发动机内沉积主要集中在补燃室头部、进气道下游、冲压喷管收敛段;颗粒粒径分布对沉积位置影响很大;由试验压强-时间曲线换算得到的冲压喷管喉部沉积速率与数值模拟结果基本相符.