改进型WENO格式的一维爆轰波数值模拟

基于3种改进型加权本质无振荡（WENO）格式,以5阶的WENO-Z格式为参考,采用与经典WENO格式不同形式的非线性权重计算方法,对稳定、中度不稳定和高度不稳定状态下的1维爆轰波系统进行了数值模拟,克服了经典的WENO格式在计算临界点时精度降低的问题.数值实验表明,对于3种状态下的爆轰系统,基于Lagrange多项式重构的WENO-Zη格式计算稳定性好,模拟效果与WENO-Z格式较为一致,比较适合于爆轰波的数值模拟.对于稳定和中度不稳定的爆轰系统,WENO-NS格式和WENO-P格式在数值模拟中会出现小的数值振荡,但计算稳定;而对于高度不稳定爆轰系统,WENO-NS格式计算不稳定,其改进形式WENO-P格式计算稳定性较好.      

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

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

PETN炸药爆轰参数的数值模拟

用吉布斯自由能最小原理,通过解化学平衡方程组,求解PETN炸药爆轰产物系统的平衡组分,计算结果与用BKW和LJD方法计算的结果相近.用自编的程序从碳的石墨相、金刚石相、类石墨液相和类金刚石液相4种相态中确定出炸药爆轰产物中游离碳更可能存在的相态,并用此相态计算碳的Gibbs自由能,以WCA状态方程作为爆轰气相产物的物态方程,对PETN炸药爆轰参数作了预言,爆轰CJ点的爆速、爆压和爆温的计算结果与实验值吻合得很好.     

TVD格式在三维非对称流数值模拟中的应用

数值格式为有限体积的隐式ＴＶＤ格式，具有高分辨率、快速稳定收敛、适用于复杂型面流场的特点。该文对典型的锥形喷管和抛物线扩张段喷管内的三维流场进行了数值模拟，得出喷管内物理参数的分布和倾斜入流对喷管侧向力及其力矩影响的变化曲线，计算结果和实验数据吻合良好。该方法可以推广到任意曲面型面喷管内流场的数值计算。为进一步进行曲面型微推偏喷管的设计奠定基础。     

预制破片战斗部爆轰产物泄露数值模拟

预制破片战斗部具有破片质量和形状近似100%可控的优点,但其爆轰产物的过早泄露会影响破片的加速过程.为了分析爆轰产物泄露对破片威力的影响,利用数值模拟方法对某预制破片战斗部的爆轰泄露过程进行了研究,并分析了炸药采用不同仿真算法所得结果与试验的吻合情况.结果表明:端面中心起爆时,爆轰泄露主要发生在周向,对破片飞散形态影响较大;流固耦合算法计算得到的破片速度和飞散角与试验吻合最好,最小平均相对误差分别为1.12%和0.56%;综合计算鲁棒性、计算精度、计算时间等,流固耦合算法较适用于爆轰产物泄露的研究.      

含铝炸药爆轰数值模拟研究

对含铝炸药25nm和50nm圆试验进行了二维爆轰数值模拟,标定了含铝炸药JWL状态方和反应速率方程参数,为检验含铝炸药爆轰模型及其相关参数的合理性,对含铝炸药驱动金属平板试验进行了二维爆轰数值模拟,得到了与实验值相符合的结果,分析了铝粉在炸药爆轰中的反应情况和约束条件对含铝炸药爆轰性能的影响。结果表明,铝粉主要在爆轰后期与炸药爆轰产物反应释放能量。     

炸药爆轰合成超微金刚石的数值模拟

在相变模型基础上对炸药爆轰过程中超微金刚石的合成进行数值模拟，用简单反应速率函数估算炸药爆轰反应区的热力学参数，用动力学有限元程序DYNA2D模拟爆轰产物的膨胀过程，按均匀成核和长大两步模拟爆轰过程中游离碳相变成金刚石的过程，同时考虑金刚石在爆轰产物膨胀过程中的石墨化。计算了不同装药条件和保护条件下金刚石的成核率，长大线速度，石墨化率，金刚石得率及粒径分布等信息，对计算结果进行分析和讨论。     

旋转磁场电磁搅拌数值模拟

对旋转磁场搅拌凝固过程中金属液的流动进行了数值模拟,结果表明,金属液所受Ｌｏｒｅｎｔｚ力切向分量及径向分量沿半径方向逐渐增大,但数值大小相差悬殊;切向流动速度沿半径方向呈抛物线规律变化,最大值发生在０．６ｒ０处,这种三维流动将有利于温度、成分均匀化的作用。     

石油射孔弹爆轰波相互作用数值模拟分析

采用动力学有限元方法,以某型号的射孔弹为例,对两发射孔弹的顺序起爆、爆轰波的相互作用效应进行数值模拟,得到相互干扰对聚能射流的影响规律,从而为提高射孔效率提供一定技术支持.     

矿井底板突水的数值模拟

本文根据岩石力学和渗流力学理论，用广义的运动稳定性原理，建立了矿井底板突水的数学模型和突水发生的判别准则，并进行了数值模拟。     

凝聚相炸药爆轰波冲击空穴塌陷过程的高精度数值模拟

利用高精度数值格式对凝聚相炸药爆轰波冲击空穴塌陷过程进行了研究. 采用RGFM和Level-set相结合的物质界面处理方法,克服了爆炸流场中高密度比、高压强比的物质界面容易引起数值振荡这一问题. 采用5阶WENO高精度有限差分格式进行空间离散,采用3阶TVD Runger-Kutta法进行时间离散,开发了高精度动态并行程序,并运用该程序对凝聚相炸药爆轰波冲击空穴后的塌陷过程进行了数值模拟,给出了PBX-9404炸药在不同入射压力下的起爆距离及冲击起爆过程中的压力历史曲线图,并分别与实验结果对比后进行了验证. 对凝聚相炸药中的空穴塌陷过程进行了高精度数值计算,给出了不同时刻的密度、压力和粒子速度云图,研究结果表明由于空穴的存在,爆轰波的传播过程中会在空穴处形成射流现象以及密度和能量的汇聚,进一步产生高温高压的热点.      

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

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

高分辨率化学-生物地层学数值旋回模拟——以广西中-上泥盆统为例

应用数学分析方法并结合化学-生物地层学手段,对广西巴漆、谷闭、堵金等3个中-上泥盆统剖面进行了系统的化学元素分析测试数据的数值旋回模拟。分析结果表明:能反映一定区域环境意义的相关元素的组合是客观稳定的;数值分析可在地层剖面上揭示出稳定、清晰的波动旋回;这些旋回可以将牙形石谱系带进行更精细的带内划分,并且可在不同的剖面间进行牙形石谱系带内旋回与旋回之间的对比,提高了牙形石谱系带内地层划分对比的精度和分辨率。     

基于ENO格式的一维溃坝水流数值模拟

将ENO格式和Runge-Kutta时间离散的思想应用于一维Saint-Venant方程组的求解,数值模拟溃坝洪水,得出了水位和流速的沿程分布.经与理论解比较可知,数值解在间断渡附近没有出现数值振荡,水位和流速大小与理论解均吻合较好,表明采用ENO格式所建立的高分辨率模型能够很好地模拟溃坝波的演进过程.     

损伤炸药的冲击起爆数值模拟

通过对炸药化学反应速率方程的分析,建立了基于KIM弹粘塑性球壳塌缩热点模型原理的三项式整体化学反应速率方程模型. 运用遗传算法确定了反应速率方程相关参数,通过与Forest-Fire反应速率模型数值模拟结果对比验证所建模型的合理性. 将所建反应速率方程模型嵌入有限元程序对PBX炸药起爆过程进行数值模拟,数值模拟结果与试验结果吻合较好,可以描述分析受冲击加载造成孔隙率、颗粒尺寸等变化的损伤炸药的冲击起爆过程.     

六种微物理方案对湖北一次强对流过程的数值模拟

利用中尺度数值模式WRF(weather researchforecasting model),选定六种微物理方案(Lin、WSM6、WDM6、Morrison、Milbrandt-Yau、NSSL),对发生在2007年4月15日湖北地区的一次强对流天气进行模拟,对比分析了各方案所模拟出的降水差异及成因。结果显示:六种方案模拟出的降水区域均偏西,其中Milbrandt-Yau方案降水偏少,其余五种方案则偏多;不同方案所模拟的云中微物理结构不同,云中冰相过程的发展很大程度上影响了降水强度,降水主要来源于大的雹霰粒子下落融化。雹霰粒子半径不同是引起降水差异的主要原因,大粒子含量高,则会产生较强的降水;反之,则降水较少。     

惰性气体泄爆实验与数值模拟

对惰性气体在3种不同压力下的泄爆进行了实验,并分别获得了流场图像、实验中气源出口截面的气体静压与周围环境静压之比分别为1．2、2．0和2．5。同时应用改造后的高阶精度MUSCL差分格式,对泄爆流场进行了数值模拟,将计算结果和实验结果对比后,发现两者吻合。根据射流动力学理论,分析了产生这种波系结构和流场压力变化的原因。当压力比为1．2时,波系为驻波形式;随着压力比的增加,射流波系中出现激波结构;当压力比达到2．5时,流动变为高度欠膨胀,波系中的激波在强度和直径上得到了加强。结果对进一步研究可燃性介质的泄爆过程,提供了可供参照的图像和数据。     

一种基于WENO格式的一维溃坝水流数值模拟研究

采用一维Saint-Venant方程组,应用WENO格式和Runge-Kutta时间离散的思想,进行溃坝水流的数值模拟,得出了水位和流速的沿程分布,并与理论解比较,发现数值解在间断波附近没有出现数值振荡,水位和流速数值解与理论解吻合较好,表明WENO格式是一种进行溃坝水流模拟的非常理想的差分格式.     

结合高分辨率多光谱影像和LiDAR数据提取城区建筑

结合IKONOS影像和LiDAR数据,应用面向对象分类分析方法试验提取高分辨率多光谱卫星影像城区建筑物.处理流程包括如下步骤:(1)影像融合增强;(2)影像分割;(3)影像对象分类;(4)建筑物对象几何形状规则化处理.试验结果表明,面向对象分类分析是一种适于利用高空间分辨率遥感影像数据进行城区建筑物制图的有效方法.该方...     

改进型脉冲耦合神经网络高分辨率SAR图像分割

针对高分辨率SAR (合成孔径雷达)图像噪声强, 目标分割难度大的特点, 提出一种改进的脉冲耦合神经网络(pulse coupled neural network, PCNN)模型的SAR图像分割算法。首先根据SAR图像中相干斑噪声的特点, 采用复小波进行去噪。然后, 在传统PCNN模型的基础上, 对神经元的输入信号, 尤其是链接系数和 阈值的非线性衰减子因子进行了改进和简化, 同时对链接强度系数β进行理论上的近似推导, 并减少人工设置的参数。最后, 通过最佳阈值对其结果进行二值化处理得到感兴趣的目标图像。实验结果表明, 改进后的算法运行效率提高, 自适应性增强。与传统算法相比, 区域一致性提高0.013, 区域的对比度提高0.015, 效果优于传统的PCNN算法, 为高分辨率SAR图像分割提供了一种新策略。