模拟冲击波传播规律的自由拉氏法

 介绍了模拟可压缩流体运动规律的自由拉氏法，编制了一维程序，对激波管问题和飞片打靶问题进行了数值模拟，并将数值结果与传统拉氏方法的结果进行了比较。     

水下爆炸冲击波在含吸收层结构中的传播规律的数值模拟

 利用ANSYS/LS-DYNA软件,对水下爆炸冲击波在含有吸收层的多层结构中的传播规律进行了数值模拟。结果表明：冲击波在吸收层中的衰减是由泡沫材料自身的本构引起的;在含泡沫材料吸收层的结构中峰值压力仅为不含吸收层的14%,冲击波传播波形发生改变;在含有泡沫混凝土的多层结构中应力波峰值的下降幅度比含泡沫铝的大,但冲击波传播波形的变化不是很明显;泡沫混凝土中的总能量大于泡沫铝中的总能量,因此泡沫混凝土作为防护吸收层要优于泡沫铝。     

气液爆轰波及冲击波传播的数值研究

采用高分辨率ＴＶＤ格式、ＭａｃＣｏｒｍａｃｋ格式和二维轴对称气液两相方程组数值模拟了气液爆轰的发展过程及冲击波误减规律,计算结果与国内外研究结果符合良好。     

岩石中爆炸波传播的数值模拟

采用一维球对称流体弹塑性模型模拟了一系列ＴＮＴ填实爆炸引起的周围岩石动力学参量的变化,根据计算结果拟合了岩石中峰值应力、粒子速度、加速度及位移的近似公式,并与一些文献中提供的实测结果进行了比较,说明数值模拟可为爆破工程、安全评估提供参考。     

激光驱动冲击波在铝-金阻抗匹配靶中的传播稳定性

 针对神光Ⅱ第九路激光条件,利用1维JB程序和多台阶靶技术对冲击波在铝 金阻抗匹配靶中的传播稳定性进行了理论和实验研究,实验结果与理论结果基本吻合。结果表明,激光驱动冲击波在铝 金阻抗匹配靶中传播时,高阻抗待测材料金中的冲击波稳定传播最大距离急剧减小。因此,在进行铝-金阻抗匹配靶物理参数设计时,应保证高阻抗材料金台阶厚度满足冲击波传播稳定性,然后再按照阻抗匹配实验中两种材料的冲击波速度比来确定低阻抗标准材料铝的台阶厚度。根据神光Ⅱ第九路激光条件,铝-金阻抗匹配靶中铝基底厚度选取为30 μm左右较好,金台阶和铝台阶厚度应分别小于10 μm和17 μm。     

三维电磁扩散场数值模拟及磁化效应的影响

采用矢量有限元法实现了三维电磁扩散场数值模拟,并成功将其应用在大地电磁的正演研究中.为灵活精确地拟合起伏地形和地下不规则构造,采用由不规则四面体单元组成的非结构化网格,可根据模型设计的需要调整网格的大小.引入了基于二次场理论,将解析的一次场从总场中扣除,直接计算二次场,使得误差仅局限于相对较小的二次场,以提高总场计算精度.常规的节点有限元法不满足电性分界面上法向电场不连续和无源区单元内电流密度无散,违反麦克斯韦方程组.为克服节点有限元法的弊端,使用矢量有限元法求解基于二次电场的偏微分方程.另外,在算法设计中,考虑了磁导率参数的变化,可以模拟磁导率不均匀的模型.通过与COMMEMI模型已发表的结果对比,证明了本文算法的正确性和精确性.为突显非结构网格优势,计算了椭球异常体模型和任意地形模型的MT响应,并详细讨论了地形和磁化效应对三维数值模拟结果的影响.     

冲击波加载速率对微射流影响的数值模拟

材料表面在发生熔化前，微射流可能是微喷射的主要物理机制之一。曾鉴荣等在纯铅的实验中发现，当靶板中出现三波结构（即弹性先驱波、相变波和塑性波）时，测得峰值压力为22GPa时纯铅样品的微喷射量比峰值压力为20GPa的单次冲击加载喷射量几乎减少了1／2。Asay在铝平面样品的微喷射实验中，也发现随着冲击波加载速率的减小（上升沿宽度增加），喷射量大致按指数规律减小。对于自由面上缺陷平均尺度为5lain的样品，在冲击加载变到35ns波阵面宽度的加载条件时，喷射量约降低了2个数量级。     

环境温度对高压储氢罐泄漏扩散影响的数值模拟

基于FLUENT软件的物种传输与反应模块建立了高压储氢罐泄漏扩散的模型,提出了研究高压储氢罐泄漏扩散的数值模拟方法.考虑到不同地区和季度温差较大,对不同环境温度下高压储氢罐发生泄漏扩散进行数值模拟.模拟结果表明:随着环境温度的升高,氢气泄漏后扩散的速度逐渐增大,但从数值上看相差较小;在不同环境温度下,泄漏口处射流速度基本相同.     

压力恢复系统超声速扩散段三维流场数值模拟

 采用数值模拟的手段，对压力恢复系统超声速扩散段三维流场进行研究。数值模拟使用LU分解法和NND差分格式求解全Navier-Stokes 方程，并加入了湍流模型。对得到的流场结构进行了分析，为下一步工作打下了基础。     

颗粒湍流扩散的数值模拟

1引言颗粒湍流扩散是气固两相流中的一个重要的研究课题,近年来,国外有大量文献报道这方面的工作。因为颗粒的湍流扩散的研究的一个基本问题就是湍流的描述,它一方面可以来检验湍流描述的正确性;另一方面,该问题是气固两相流理论中,特别是颗粒群轨道模型中一个难点。描述湍流大体上可分为两种方法。一是直接数值模拟方法（DNS）,J．B．Mclaughlin对管流内颗粒沉降的研究就基于直接数值模拟山,K．D．squires和J．K．Eaton对颗粒扩散的研究也基于DNS［2］,DNS一般只适用于低Reynolds数。另一方面,用随机方法来模拟湍流,该…     

Cu(111)三维表面岛对表面原子扩散影响的分子动力学研究

采用嵌入原子方法的原子间相互作用势，利用分子动力学方法计算了同质外延生长中不同层数的三维Cu(111)表面岛上表面原子扩散激活能，分析了三维表面岛的层数对表面原子交换扩散和跳跃扩散势垒的影响. 研究结果表明，二维Enrilich-Schwoebel(ES)势垒小于三维ES势垒，且三维ES势垒不随表面岛层数的增加而显著变化. 对于侧向表面为(100)的表面岛，表面原子沿〈011〉方向上的扩散行为，随表面岛层数增加而逐渐变化；在表面岛层数达到3层时，扩散路径上的势垒变化趋于稳定，表面原子扩散以下坡扩散为主. 对于侧面取向为(111)的表面岛，当表面岛层数大于3层后，开始呈现上坡扩散的可能. 关键词： 表面原子 扩散 分子动力学模拟     

核爆聚变电站发电时的冲击波超压传播规律

 核爆聚变电站在发电时会产生冲击波,其超压传播规律是核爆聚变电站洞壁设计的重要基础。现有公式在描述冲击波峰值超压时存在一定的不足（如相对误差大等）。通过理论分析得到新的冲击波峰值超压公式、正压作用时间τ₊和冲击波传播系数ω,最终得到冲击波超压的传播规律。通过对比分析表明,新公式比现有公式更加适合求解核爆聚变电站发电时所产生的冲击波峰值超压。     

煤矿巷道瓦斯爆炸冲击波与高温气流的关系

 为了给瓦斯爆炸后煤尘二次爆炸的深入研究提供理论依据,应用计算流体力学方法,对煤矿巷道内瓦斯爆炸的瞬态流场进行了数值模拟,得到了冲击波与高温气流流动的时空关系,并借助实验对数值方法进行了验证。研究表明：在瓦斯爆炸后的一定时间内,近场区域和远场的部分区域极有可能引发煤尘二次爆炸。给出了可能发生煤尘二次爆炸的区域随瓦斯区长度的函数关系式,以及远场中峰值温度和峰值超压到达时间的间隔随轴向距离和瓦斯区长度的分布特性。     

气液爆轰波及冲击波传播的数值研究

采用高分辨率TVD格式、MacCormack格式和二维轴对称气液两相方程组数值模拟了气液爆轰的发展过程及冲击波衰减规律,计算结果与国内外研究结果符合良好。     

爆炸冲击波作用下预制孔靶板塑性变形规律的研究

利用非线性动力学软件AUTODYN,对爆炸冲击波作用下无孔及预制孔靶板的塑性变形进行了数值模拟,获得了不同孔数、孔径情况下靶板中心点挠度的变化规律;设计了模型实验,其结果与数值结果符合较好;并以爆炸冲击波作用下无孔靶板中心点挠度计算公式和数值计算结果为基础,通过拟合给出了预制孔靶板中心点挠度与孔密度、孔径之间的函数关系,该关系可为有孔平板目标的毁伤评估及破片、冲击波对平板目标联合毁伤研究提供参考。     

冲击波垂向加载下PZT 95/5铁电陶瓷的去极化和放电过程分析

电极化后的PZT 95/5铁电陶瓷能够在冲击波作用下快速去极化并释放束缚电荷,形成高功率的瞬态输出电能。对于垂直于极化方向的冲击波加载情况,通过将去极化过程中的铁电陶瓷等效为电流源、电容和电导的并联电路,综合考虑冲击波压力对波速和去极化相变过程的影响,以及冲击波前、后铁电陶瓷的介电常数和电导率变化,建立了描述冲击波垂向加载下PZT 95/5铁电陶瓷去极化和放电过程的模型,解析获得了铁电陶瓷的放电电流表述。在此模型基础上,开展了短路和电阻负载条件下PZT 95/5铁电陶瓷在冲击放电过程中的输出电流特征分析,并与相关实验结果进行了对比。结果表明:模型能较好地模拟实验观测的铁电陶瓷PZT 95/5的冲击放电过程,以及冲击波压力、负载电阻等对冲击放电输出电流的影响规律。     

基于SIMPLE求解对流-扩散方程的一种新方法

针对SIMPLE系列算法,通过分别求解控制容积界面和节点两个位置的对流-扩散方程来提高数值计算精度,提出了一种流动与传热数值模拟的新方法。采用新的数值方法对顶盖驱动流与正方腔自然对流进行了数值模拟。数值模拟结果表明,在相同网格划分时,新的数值方法相对迎风格式、乘方格式、QUICK格式SIMLE算法的计算精度高;而在计算精度基本相同时,新方法有较高的计算效率。     

机匣仿生顶室对压气机间隙流动影响研究

动叶叶顶产生的间隙泄漏流对压气机的效率和稳定性产生明显影响。为了减小泄漏流对流道内部流场结构的影响,本文基于蜻蜓翅翼的翅室和褶皱结构,在跨声速压气机叶栅的机匣上布置仿生顶室结构。采用数值模拟方法,研究不同高度的仿生顶室结构对间隙泄漏量和流动损失的影响。结果表明:流体流经仿生顶室时,顶室内部存在明显的驻留涡,抑制泄漏涡沿节距方向的发展,减小间隙的泄漏流量和流动损失;在间隙高度为2 mm情况下,仿生顶室对泄漏量和流动损失的控制效果最为明显,最佳方案能够使总压损失系数下降1.2%,间隙泄漏量降低1%。     

马赫数对等离子体碰撞冲击波结构的影响——全动理学模拟研究

利用全动理学粒子(PIC)模拟研究不同马赫(Mach)数条件下等离子体碰撞冲击波的结构。研究发现：在低Mach数条件下,冲击波波阵面位置物理量的分布较为平缓(努森数较小),等离子体的粘性和热流可由经典输运理论描述,这种情况下数值求解双流体方程组得到的冲击波结构与PIC模拟一致。随着冲击波Mach数的增加,冲击波波阵面位置物理量的分布变陡,努森数增加,动理学效应对等离子体输运的影响变得显著。高Mach数条件下,动理学效应主要体现在“先驱离子”对离子粘性和热流的增强以及电子非局域输运对电子热流的影响。通过影响等离子体输运行为,动理学效应可以显著影响冲击波的结构特征。