拉氏程序使用减敏模型模拟炸药绕爆问题

利用二维拉氏非结构网格爆轰流体程序对高能炸药的爆轰波绕爆问题进行数值模拟.拉氏程序对绕爆问题进行模拟时,会引起网格严重扭曲,导致网格大变形.使用JWL状态方程和三项式反应率,拉氏程序得到的绕爆数值结果与实验定性相同,较好地模拟了爆轰波绕爆过程中的漩涡现象.为了模拟一些钝感炸药在绕爆过程中的死区现象,在炸药反应率中加入减敏因子,模拟结果捕捉到死区,并与国外文章中的数值模拟结果定性相同.     

X射线衍射法(XRD)分析煅烧白云石的物相组成

采用X射线衍射技术(XRD)对白云石原矿和不同温度煅烧的白云石原矿进行了物相分析,白云石原矿中主相为CaMg(CO3)2,相对含量约为98.6%,杂质相为SiO2,相对含量约为1.4%。白云石原矿经过300和500℃,煅烧2h未见主相CaMg(CO3)2分解,在750℃煅烧2h,主相CaMg(CO3)2部分分解为碳酸钙和氧化镁,在1000℃煅烧2hCaMg(CO3)2全部分解为氧化镁和氧化钙。     

科学计算程序的验证与确认

科学计算程序的验证和确认（Verification＆Validation，简称V＆V）是数值模拟的一个重要组成部分，其目的是为了评估模拟程序和物理模型的可靠性，并量化数值模拟程序计算结果的置信度。这是解决全系统试验无法进行或者代价十分昂贵情况下，确保数值模拟置信度的一种系统方法。     

任意区域上的粒子均匀分布方法

给出了任意区域内粒子的Voronoi面积均匀划分的方法,对给定的区域边界点,先由边界点得到边界粒子,再利用分布函数结合Delaunay加点算法得到内部粒子,这种方法可以得到一般的区域上粒子的均匀分布。这里给出了一个流体力学中利用SPH方法模拟内爆问题的例子,算例表明,均匀分布可以得到更高的无网格方法数值模拟精度。     

老版本FORTRAN程序的结构化改造

许多程序是在早版本的Fortran77语言或Fortran60环境下设计的，也由于计算机计算存储能力和计算速度限制，程序不是结构化的。这给程序维护与改造造成了很大困难。在充分分析程序特点的基础上，尝试编制软件对老程序进行结构化改造，取得了很好的效果。     

氧的低温吸附模型与计算

受化学氧碘激光器(COIL)高效率小型化及新型压力恢复技术发展的推动,考虑COIL中气体氧的低温吸附问题。从质量守恒观点出发考察了气体流过沸石多孔床的压力与速度的关系,建立了吸附剂吸收的物质与气流中物质关系的吸附平衡方程。由吸附物较少的特点,得出在确定的压力梯度下速度为常数的近似。由吸附平衡方程与吸附速率方程,构成低温吸附模型方程,其中,亨利等温线方程作为辅助关系。对低温吸附模型方程作数值求解,得到了气流中氧质量浓度及吸收的氧质量浓度随时间变化的规律,以及这些质量浓度分布的长时间渐近特征。     

辐射流体力学Lagrange方程组一类人为解构造方法

针对一维Lagrange辐射流体力学方程组,基于物理量在Euler空间和Lagrange空间的微分关系,提出了一种人为解构造方法,并构造了一类一维Lagrange辐射流体力学方程组人为解.构造的人为解在整个计算区域光滑可微,质量方程无源项.将构造的人为解应用至二维辐射流体力学Lagrange程序中,从数值误差、收敛阶方面验证了程序的正确性,展示了人为解的可行性和适用性.     

针对流体大变形问题网格邻域可变技术

结合大量实际问题,研究基于邻域可变的拉氏局部重分新方法,允许网格邻域在计算中局部跟随流体可变,虽然可能会导致任意多边形网格,但比标准的拉氏+重分方法有更好的拉氏特点.网格邻域可变的局部重分策略是基于计算过程中处理网格相交的思想,即在网格相交处改变网格邻域关系,使网格不会出现相交,保持计算正常进行.通过几个数值例子论证方法的特点,并与标准的拉氏+重分方法进行比较,显示其优点.     

爆轰计算JWL状态方程参数的不确定度

针对爆轰计算JWL状态方程参数,开展不确定度量化研究。首先从JWL状态方程性质出发,结 合数学、物理分析确定JWL状态方程中主要参数R1 的概率分布范围,然后利用非嵌入多项式混沌展开方法 对PBX9404炸药驱动飞层外界面速度进行了不确定度量化分析。给出的分析方法与分析流程对爆轰计算及 其他工程数值模拟不确定度量化有一定参考作用。     

TVD自相容流体力学拉氏格式研究

流体力学数值模拟格式总体上可分为Eulerian(欧氏)、Lagrangian(拉氏)和ALE(Arbitrary LagrangianEulerian),TVD广泛应用于Eulerian格式。本文利用具有TVD保持性质的Runge-Kutta型时间离散方法,构造了流体力学Lagrangian(拉氏)自相容格式,应用von Neumann小扰动技术分析了该格式的稳定性,并进行了相应的数值模拟,较好地抑制了激波波后非物理振荡。