驱动管中柱状装药爆轰过程的数值模拟

数值模拟了爆炸驱动管中柱状装药内爆轰波的传播过程.计算采用欧拉型有限体积方法,炸药及爆轰产物均采用JWL状态方程,空气采用理想气体状态方程,采用"点火-生长"模型计算化学反应速率.计算得到了驱动管内波系结构的发展过程,爆速与经验公式符合得较好.计算表明,驱动管侧壁的压力峰值在800MPa以上,而在管底中心处,由于激波的汇聚,压力峰值高达12.4GPa.     

基于格结构假设空间证据推理的双向传播计算模型

在专家系统中,基于证据理论的推理网络一般为格结构形式,将证据理论的现存形式直接应用于具有复杂结构的推理过程存在很大的困难。本文基于格结构假设空间,采用不确定知识的可信度表示方法,给出了初始基本概率赋值函数的一种计算方法,并在扩充的Shenoy-Shafer结点处理模型的基础上,提出了证据推理的双向传播计算模型。     

基于新型非结构有限体积解耦算法的激波诱导燃烧数值模拟

把基于结构网格有限差分方法建立的化学非平衡流动模拟的新型解耦算法推广到基于非结构网格的有限体积法,从而可用于模拟复杂几何构型下的化学反应流动.对H2/ Air预混气体中激波诱导振荡燃烧的Lehr实验进行了数值模拟,计算得到的振荡频率与实验结果符合很好,表明计算方法具有时间和空间二阶精度.通过对不同几何外形发射体的计算发...     

热射流起爆过程的数值模拟

采用改进的化学非平衡流解耦方法处理Euler反应流方程,并以H2/air 9组分21方程模型对竖直喷注的热射流起爆过程进行了二维数值模拟。对流项采用五阶WENO离散格式,时间推进采用二阶Runge-Kutta方法。详细分析了热射流以不同速度、入射位置、入射宽度和入射倾角喷射时在爆震管内形成的流场,总结了射流参数状态影响起爆的一般规律;解释了激波、火焰的相互作用对"热点"的形成以及转变为爆震的影响,特别是激波反射对"热点"形成的促进作用。结果表明,爆震波通常由火焰面和固壁附近狭长未燃区域中的"热点"产生。"热点"向爆震波发展的过程是过驱爆震阶段,存在三波结构。为了实现快速起爆,应当增大射流入射速度、贴近侧壁并以适当入射角度喷射热射流。     

一种新型的计算化学非平衡流动的解耦方法

给出一种计算化学非平衡流动的解耦方法。利用该方法从薄层近似N S方程出发 ,采用NND2M差分格式 ,数值模拟了H2 /O2 燃烧的超音速绕流场 ;计算表明这种方法可以有效地解决刚性问题 ,而且与现有的隐式算法相比 ,具有公式推导简便 ,占用内存少 ,计算效率高等优点 ,特别适合于多组分和多反应的情况     

伪谱法模拟各向异性/非均匀介质中弹性波的传播

数值模拟一直是研究弹性波在复杂介质中传播问题的有效手段之一。各向异性问题和非均匀问题日益引起关注。伪谱法是近年来发展较快、建立在弹性波动方程基础上的一种数值模拟方法。它的主要特征是将波动方程中的空间微分变换成频域中的乘法运算 ,从而只需要较少的空间格点就可以得到较高的计算精度。介绍了伪谱法及其稳定性条件 ;编制了一个适合于小型工作站或者PC机的模拟地震波在二维各向异性均匀 /非均匀介质中传播的数值模拟程序 ,模拟了各向异性 /非均匀介质中弹性波传播理论 ,分析了裂纹以及各向异性界面对波传播的影响     

装药结构对破片初速影响特性

为了满足低附带战斗部的毁伤效果,将炸药与破片混合形成新型混合装药模式。在混合装药模式中,为了研究不同装药结构对破片初速的影响特性,设计了3种不同装药结构并进行数值模拟和试验;对比发现,数值模拟与试验结果规律吻合,破片的初速与装药结构有一定的关系。根据试验数据拟合了3种不同装药结构中破片最大速度曲线对比图,通过分析得出:不同装药结构中,炸药与破片分层越多,爆轰瞬间,相邻炸药所产生爆轰效果会发生叠加,爆轰波对破片的驱动效果会增强,则破片的初速也就越大。     

飞船再入舱三维化学非平衡流数值模拟

对飞船再入舱高超声速化学非平衡三维流场进行数值模拟。控制方程为含化学反应源项的全Ｎａｖｉｅｒ－ Ｓｔｏｋｅｓ 方程组, 化学模型为高温空气的七组元模型, 其组元成分为Ｎ２ 、Ｏ２、ＮＯ、Ｎ、Ｏ、ＮＯ＋ 、ｅ－ 。差分格式采用张涵信院士提出的ＮＮＤ激波捕捉格式并用全耦合方法处理流动方程及化学反应方程, 对化学反应源项采用全隐式处理, 运用时间预处理技术以加快收敛。同时, 也对完全气体模型进行计算, 以分析真实气体效应。与实验数据对比, 说明本文结果是可信的。     

有限水深Kelvin源格林函数及其导数的快速计算

应用线性船舶水动力学理论,将有限水深Kelvin移动兴波源格林函数分解成3部分:简单Rankine源集合、局部扰动项和波函数项。对各部分在亚临界和超临界航速时的快速计算分别作了讨论,其中对计算费时的局部扰动项采用了复数中的最陡下降积分法,加快了计算速度,进而引出各部分偏导数的计算。最后,比较了格林函数的差分值与偏导数值,计算了薄船兴波阻力与表面兴波,讨论了不同面元网格划分的收敛性,证实了亚临界、超临界航速兴波分别以横波、散波为主。     

水下结构声辐射FEM/BEM简化计算方法研究

为提高水下结构声辐射FEM/BEM方法的计算速度,首先以节点位移作为中间变量对湿表面处结构和流体的运动方程进行解耦;然后,采用单元的外法向矢量来表示相关系数,并简化了积分奇异性的处理;最后,以脉动球和水下加筋圆柱壳为算例进行了计算,计算结果、理论结果与试验结果吻合较好,由此验证了方法的有效性。