轴对称地维氢氧爆轰驱动激波管的模拟

对由前向爆轰产生驱动气体并具有局部锥形收缩截面的轴对称激波管利用两相常比热完全气体的欧位方程组和有限体积TVD格式在局部下在交的网格上进行了二维数值模拟。将消除两种介质界面处数值振荡的新方法成功地推广到二维情况。分析了该激波管力学特性。     

激波诱导变形液滴外流场数值模拟

基于Delaunay非结构网格生成方法和局部阵面推进重新生成网格的方法,应用边界边剖分与合并技术,对二维含自由表面流动中运动边界网格的生成方法进行了研究,改善了局部网格重新生成过程中的网格质量;对激波诱导变形水液滴在高速气流作用下的受力情况做了数值模拟.结果表明:该文在处理二维自由表面的动边界问题中采用的方法是可行的.     

变截面激波管中激波运动规律的初步数值研究

利用有限差分方法求解变截面激波管中存在激波的二维流体流动问题。着重研究初始流场为静止和超音速两种情形。分析讨论所得计算结果,并和已有实验结果进行了比较。     

激波-多涡列-声场相互作用及其POD分析

采用5阶WENO格式,通过数值求解二维非定常欧拉方程,模拟了激波与四涡组成的涡列的相互作用过程.通过研究马赫数为1.05的激波与强度为0.25的多涡列相互作用,进一步揭示了激波-涡列相互作用的动力学特性以及声波的产生过程.     

双曲型守恒律的一类局部化的高效差分格式

构造了一维非线性双曲型守恒律的一类局部化的高效全离散差分格式,并将该格式推广到一维守恒方程组及二维守恒方程(组).最后,给出了几个标准算例.数值计算结果表明此格式具有高精度高分辨激波、稀疏波和接触间断,且边界条件易于处理等优点.     

二维激波边界层干扰的数值分析

该文利用二维可压缩雷诺平均Navier-Stokes方程，采用k-ε湍流模式，运用半离散有限体积法，对捆绑火箭助推器和芯级之间的激波边界层干扰进行了二维数值模拟．数值结果表明，在激波边界层的干扰下，壁面边界层流动出现了分离和再附过程，压力和温度出现激烈变化．     

入射和反射激波与火焰相互作用的实验和数值显示

以等物质的量比的甲烷/空气为研究对象,利用YA-16多闪光高速照相机对火球在入射激波及其反射激波作用下的变形进行了实验研究.基于带化学反应的二维轴对称Navier-Stokes方程,采用带有Superbee限制函数的波传播算法,对该现象进行了数值模拟,并用计算光学对数值结果进行了处理,得出的计算阴影图和实验阴影照片比较吻合.结果表明,火球在入射激波的作用下,出现了Richtmyer-Meshkov不稳定、斜压效应和Helmholtz不稳定,导致火球不断变形,在局部区域还出现了湍流,加剧了燃烧.随后,反射激波再次作用火球,加速火焰变形的同时,还使火球在波后高能量的可燃气中剧烈燃烧,湍流强度也大幅度增加,同时由于壁面的限制,使得火球发展呈蘑菇云形状.而激波穿越火焰时,轴线上的激波强度也会得到加强.     

含强激波二维流动的高分辨率数值模拟

采用A.Harten的高分辨率总交差不增(TVD)格式,对激波反射问题和带有台阶的二维风洞中高超声速绕流问题进行了数值模拟。计算结果表明:TVD格式比传统的数值格式如Godunov,BBC,Muscl,及PPM等有较高的激波分辨率,且伪振荡小。     

柱状汇聚激波冲击重气柱不稳定性的数值研究

采用有限体积法结合网格自适应技术的VAS2D程序,数值研究了二维重气柱在柱状汇聚激波及其反射激波作用下的Richtmyer-Meshkov(RM)不稳定性发展,并比较了平面激波条件和不同汇聚角度条件下激波和界面的演化过程.气柱界面内气体为六氟化硫(SF6),环境气体为空气,柱状汇聚激波在汇聚中心形成反射激波,平面激波与右端固壁作用形成反射激波.结果表明,柱状汇聚激波及其反射激波作用下的气柱演化、界面位置变化、尺寸变化及环量变化与平面激波及其反射激波作用下的情形有着显著差别,同时,不同汇聚角度也具有一些影响,揭示了柱状汇聚激波及其反射激波作用下的界面演化机理.     

利用激波管进行水蒸汽凝结过程的理论和实验研究

水蒸汽凝结过程是一种含有潜热的两相非等熵过程.本文在一台激波管装置上,对水蒸汽的凝结问题进行了实验研究.本文采用经典成核理论,将特征线方法应用到水蒸汽凝结过程,进行了理论分析和数值计算.本文的实验结果与国外其他作者在激波管上所作的同类结果能吻合较好.     

Ghost Fluid方法在可压缩多介质流动中的应用

在对可压多介质流动进行数值模拟时,界面两边流体状态方程不同会给计算带来很大麻烦,该文采用Ghost Fluid方法(GFM)成功地解决了这一问题.利用间断有限元Galerkin方法和Level Set方法,对2种不同比热比介质激波管问题进行了模拟,其数值结果和精确解非常吻合;对气体/水激波管、空气/氦气泡相互作用、空气/R22相互作用等问题进行了模拟,其结果与相关文献的结果一致,表明利用该方法进行可压缩多介质流动数值模拟是可行的.     

冷凝激波管内水,乙醇凝结两相流的理论和实验研究

对冷凝激波管内的水、乙醇蒸汽凝结两相流进行了较深入的理论与实验研究，蒸汽在自发凝结时，伴有潜热释放，气液两相之间产生传质、传热及动量传递，从而存在耦合效应。本文从两相流的基本方程组出发，导出非定常凝结两相流的特征线及其相容方程组，从而对激波管内凝结两相流进行了理论分析及数值计算。在一凝结激波管试验台上分别对水、乙醇、苯３种工质的凝结两相流进行了实验测试，用光测系统对液滴进行光散射测量，并引入先进的     

战术火箭子母战斗部第一次抛撒分离多体干扰流场数值模拟

以三维Navier-Stokes方程为基础,用计算流体力学(CFD)数值模拟方法对俄9M55K型火箭弹子母战斗部的第一次抛撒分离过程--子弹筒从母弹侧向抛出进行了研究.数值模拟结果表明,当子弹筒离母弹很近时,子弹筒头部激波在母弹上的反射是母弹对子弹筒气动干扰的主要原因;当子弹筒穿越母弹头部激波时,母弹头部激波是引起子弹筒气动特性剧烈变化的主要原因;多个子弹筒相距较近时,子弹筒头部激波之间的彼此干扰也很重要.这些结论可为子母型武器的总体设计、分离参数选择提供气动依据.     

激波数值摸拟的一个简单有效的方法(英文)

本文全面讨论了用拟特征线法求解激波问题的方案、格式和具体算例.尤其是在本文中首次提出了采用分裂的拟特征线法计算多维激波问题的简便格式.其格式的简单、计算量的节省和程序实现的灵活多样是其它许多方法所没有的.最后以Mach 反射为例做了实际计算,所得结果与已有的成功结果相比也是满意的.     

振荡管内非定常流动的数值模拟

建立了振荡管内非定常流动与传热的数学模型 ,分别采用Godunov格式和控制容积法对气动方程及传热方程进行离散 ,得到了管内压力波形、气流速度及温度随时间的变化历程 .所得管内压力波在波形和相位上与实验结果吻合较好 ;排气初期的气流速度最大且温度最低 ;管内出现强反射激波 .提出的求解Riemann问题迭代初值可使迭代快速、稳定收敛     

长持时爆炸冲击波作用下钢梁动力响应分析

大规模爆炸时引起的长持时平面爆炸波会造成结构的整体破坏,这与接触爆炸或近场爆炸的破坏模式有显著差异,但相关研究不多。为研究长持时冲击波作用下钢梁的动力响应,首先基于耦合欧拉-拉格朗日方法建立了大型数值激波管模型,然后采用该数值激波管得到作用在钢梁上的反射超压时程,最后通过显示动力学方法研究了冲击波入射方向、持时、超压等参数对钢梁破坏模式的影响规律。结果表明：本文提出的数值激波管模型能较好的模拟长持时爆炸波的绕射现象,可得到可靠的反射超压和冲量时程;随着冲击波入射角度的增加,简支钢梁跨中位移增大;冲击波超压峰值相同时,随着持时的增加,简支梁跨中的最大位移有增大的趋势,而最大位移的增长速率则会变缓;简支钢梁在长持时爆炸冲击波作用下容易发射受弯破坏,而梁跨中翼缘板会出现局部屈曲现象。     

变截面管道对瓦斯爆炸特性影响的数值模拟

以气体爆炸理论为基础,利用高精度的加权本质无振荡(WENO)格式对变截面管道中的瓦斯爆炸进行了数值模拟,探讨了变截面管道对瓦斯爆炸火焰传播的影响规律,得到了变截面管道造成瓦斯爆炸强度增大的结论. 在此基础上,分析了障碍物、壁面和三波结构对瓦斯二次爆炸的影响. 结果表明,经壁面反射后的激波和三波点的碰撞都能够诱导二次爆炸的产生. 这些结论为瓦斯爆炸的预防、安全评估和防火防爆提供了重要的理论依据.     

用于降低冲击噪声的非线性消声器设计

将非线性振子和亥姆霍兹消声器进行组合,设计出一种具有反射和吸收声压能量两种降噪方式的新型非线性消声器模型,用来降低声管中的冲击噪声。首先基于两种不同的组合方式,设计出两种结构不同的非线性消声器模型a和b;然后从物态方程出发,建立了非线性消声器的数学描述方程;并通过数值仿真,研究了非线性消声器降低冲击噪声的消声性能。结果表明,恰当的设计模型参数,在较宽的频带范围内,非线性消声器模型a和b均能以吸收小部分和反射大部分能量的方式有效降低冲击噪声。两个消声器模型b配合工作时,能以吸收部分能量的方式有效降低冲击噪声。     

单循环脉冲爆轰发动机的数值模拟

对高温火团引发的甲烷-空气混合气体的爆轰过程及管内、外流场做了二维轴对称数值模拟,考虑了CH4-O2-N2的详细化学反应动力学机理,该机理包含了19个基元反应和14种组分。采用分裂格式处理带化学反应的Navier-Stokes方程,使用改进的波传播方法求解流场,使用基于隐式Gear算法的微分方程解法器求解化学反应过程。模拟结果显示了爆轰波在管内的成长、稳定传播、进入喷嘴后衰减为激波以及进入外流场的全过程。还对轴线上的压力分布和封闭端的推力等进行了讨论,为脉冲爆轰发动机的开发研制提供信息。     

梯度功能材料薄壳动态热挠度的Laplace变换有限元分析

采用Laplace变换有限元法分析了梯度功能材料薄壳结构的热冲击挠度。首先从控制方程出发,推导出梯度功能材料热冲击问题的泛函及变分原理;然后对梯度功能材料薄壳的热冲击问题推导出有限元方程;利用Laplace变换引入初值条件后推导出拉氏域上的有限元计算公式,经求解后再通过Laplace数值逆变换求得时间域上各给定时刻的解答。计算实例证明了该方法的有效性,其计算结果可供工程设计参考。