燃烧系统的离散Boltzmann建模与模拟研究进展

燃烧系统的诸多模拟依托于流体建模, 离散Boltzmann方法(discrete Boltzmann method, DBM) 是近年来发展起来的一种新的流体介观建模方法. 本文简要评述DBM发展的两个方向——Navier-Stokes等偏微分方程的数值逼近解法和复杂系统的微介观动理学建模. 主要介绍在燃烧系统模拟方面DBM已有的工作、新近的思路、与传统流体建模的异同以及近期的研究成果. 本文重点传递的信息为: 作为复杂系统微介观动理学建模出现的DBM在模拟过程中同时给出“流动”及其相伴随的、关系最密切的那部分“热动”非平衡效应; 它为燃烧等复杂系统中各类非平衡行为的描述、非平衡信息的提取、非平衡程度的度量提供了一种简洁、有效的方法; 它所提供的热动非平衡测量量有两类: 一类是直接比较分布函数和平衡态分布函数的动理学矩关系得到的, 一类是来自于Chapman-Enskog多尺度分析给出的热传导和黏性项. 基于第二类DBM, 可以实现(燃烧等)一大类复杂流体系统的多尺度物理建模.     

爆轰波在炸药-金属界面上的折射分析

针对爆轰波在炸药-金属界面上折射时由实验获得的金属折射冲击波压力与经典爆轰波极曲线理论预测的压力存在显著差异这一问题, 本文展开了进一步的理论和数值模拟分析研究. 首先通过分析指出经典爆轰波极曲线理论的缺陷, 并对爆轰波极曲线理论进行了改进, 改进爆轰波极曲线理论给出了炸药爆轰波折射类型以及折射冲击作用点处的压力值. 然后发展了一个基于次特征理论来数值求解爆轰反应流动控制方程的二阶中心型Lagrange方法, 并数值模拟了一个典型的炸药爆轰波折射实验. 改进爆轰波极曲线理论和数值模拟分析结果表明, 爆轰波折射类型有三种:反射冲击波的正规折射、带Mach反射的非正规折射、无反射波的正规折射, 并且金属折射冲击波压力值随入射角增大而单调减小.     

PBX炸药粘弹性统计微裂纹本构模型的改进

基于描述PBX炸药本构关系的粘弹性统计微裂纹模型(Visco-SCRAM),通过修正其微裂纹方向矢量的取值范围、初始缺陷尺寸和体量本构关系,得到模型的改进形式。将改进的模型嵌入有限元程序中,开展PBX炸药JO-9159的平板撞击实验数值模拟研究,讨论了微裂纹数密度参数对计算结果的影响。结合实测曲线,对模型改进前、后的计算结果进行对比和分析。结果表明,改进的模型能够更合理地描述低速冲击下PBX炸药的损伤和破坏。     

学习于敏院士的科学精神

 2015年1月9日,习近平主席亲自授予于敏院士国家最高科学技术奖,表彰他在国防科研上作出杰出成就,为国家安全和科技进步作出的卓越贡献。     

一种结合Roe格式的高精度半拉氏方法

针对守恒形式的欧拉方程组, 构造了一种结合Roe格式的守恒型有限体 积形式的半拉氏方法. 通过发展一种基于Roe特征速度的拉氏质点回溯方法, 由此来计算 半点的流量并作为边界通量的近似, 使得这种半拉氏方法在时空离散上达到二阶精度, 并 且保证了守恒性. 其中回溯点处物理量采用本质无振荡格式(ENO)方法进行插值重 构得到, 不需要增加人工黏性且避免了有限体积多矩半拉氏方法中限制器选择的问题, 又能够达到时空的高阶精度. 方法简便, 易于实现, 兼具拉氏方法和欧拉方法的优点. 一维和二维数值算例表明, 此方法对激波和接触间断都取得了满意的模拟效果, 可用于可压缩复杂流动问题的计算.     

波后慢放能过程对钝感炸药爆轰过程的影响

采用唯象反应率模型与慢放能函数相结合方法,研究钝感炸药波后缓慢放能过程对爆轰过程的影响,计算结果与实验对比分析表明,慢放能反应率模型可更好地描述钝感炸药的爆轰物理过程.     

带随机初值的随机偏微分方程整体解的存在性

研究了一类带随机初值并且由分数次Brownian运动驱动的随机偏微分方程.借助于Kolmogorov准则,建立了整体Lipschitz条件下此类随机偏微分方程的一个解.同时证明了局部Lipschitz条件下整体解的存在性.     

强磁场中氢原子的能级结构

发展了一套简单高效的非微扰理论方法研究强磁场中的原子能级结构.作为例子,给出了磁场强度从0到1000个原子单位,氢原子基态和低激发态的结合能以及四极矩等重要原子能级结构参数.结果表明:相对于其他的高精度计算方法,该方法不仅能计算出高精度的能级位置而且可方便给出精确的电子波函数.此外,该法还具有很强的普适性,可直接推广到原子与任意方向的交叉电磁场相互作用的研究. 关键词： 强磁场 CWDVR谱方法 氢原子能级结构 四极矩     

Z箍缩丝阵负载区域的似稳态磁场分布

 将负载区域的电流（丝阵电流、阴极板电流、阳极板电流和回流柱电流）离散成电流线或电流面等电流微元,根据毕奥-萨伐尔定律,计算所有电流微元在指定场点的磁场,再通过叠加给出该点的总磁场。研究结果发现:在丝阵外围区域,仅由丝阵电流所产生的磁场偏离无限长直导线磁场公式的值,但全部电流所产生的总磁场与公式给出的值很接近。同时,研究了不同负载结构参数下的磁场分布,结果表明：增加丝根数有助于减小单根丝表面的局部磁场,改善丝阵外围磁场分布的均匀性。     

高速气流中激光加热平板数值模拟与分析

 采用流固耦合方法,数值模拟了高速流场中激光作用下来流速度对平板温度分布的影响。结果表明:无激光辐照时,高速气流中平板有较高的气动生热平衡温度,且平板-气流之间的换热系数随来流速度增大而增大;在平板前沿换热系数增长最快,沿平板长度方向增速趋于平缓。分析了激光辐照时高速气流中激光加热平板的温度分布情况,考察了来流速度不同时,气动生热、散热和激光辐照对平板温度的影响,给出了激光辐照后的温升情况和温度分布,分析了在不同速度来流下,对流散热、摩擦生热和激光加热之间的竞争关系,结果表明,平板温度具体分布主要是加热过程竞争的结果。