炸药爆轰的连续介质本构模型和数值计算方法

假定炸药和爆轰产物处于局部热力学平衡状态, 即它们的压力和温度相同, 利用热力学基本关系建立炸药爆轰过程的连续介质本构模型的一般理论框架. 在此框架下, 炸药爆轰本构模型由一组常微分方程构成, 包括炸药和爆轰产物的状态方程、简单混合法则、化学反应速率方程和能量守恒方程, 易于由成熟的计算方法如梯形法等进行求解. 一组广义Maxwell型非线性固体本构形式的微分方程描述了压力和温度随时间的演化速率与应变率和化学反应速率的关系, 借助简单混合物理论, 其中的系数由炸药和爆轰产物的材料参数确定. 未反应的炸药和爆轰产物采用JWL状态方程, 化学反应率方程采用Lee-Tarver点火-燃烧二项式模型, 模拟PBX-9404炸药的一维冲击波起爆过程和爆轰波传播过程. 计算结果表明了本文给出的本构模型和相应计算方法的有效性. 关键词： 炸药爆轰 本构模型 化学反应率方程 数值模拟     

松弛方法在计算凝聚炸药爆轰问题中的应用

利用松弛近似,将非线性的凝聚炸药爆轰控制方程转化为线性的松弛方程组,并采用五阶WENO格式和五阶线性多步显隐格式对线性松弛方程组进行空间方向和时间方向的离散,由此建立具有高精度和高分辨率性质的计算凝聚炸药爆轰的松弛方法。建立的松弛方法可以避免求解Riemann问题及计算非线性通量的Jacobi矩阵,同时无需分裂处理反应源项。通过对凝聚炸药的平面一维定常爆轰波结构及球面一维聚心、散心爆轰起爆和传播过程的数值模拟,验证了所建立的松弛方法能够很好地计算凝聚炸药爆轰问题。     

磺胺甲基异噁唑稀土络合物的研究

用磺胺甲基异噁唑钠(NaSMZ)与氯化稀土在水溶液中相互作用,制得了15个磺胺甲基异噁唑稀土络合物,组成为Ln(SMZ)_3·4H_2O。研究了它们的红外光谱、核磁共振谱、热谱、X射线粉末衍射等物理化学性质,并对化合物的结构进行了讨论。     

基于特征理论的二维无粘Lagrange流体力学有限体积法

提出一个求解二维无粘Lagrange流体力学方程的中心型有限体积方法.采用特征理论求解网格节点处的速度及压力,并利用这些物理量更新节点位置及计算网格界面通量.方法适用于结构网格与非结构网格.典型数值实验的结果表明,格式具有较好的收敛性、对称性、能量守恒性及鲁棒性,且能自然地求解多物质流动问题.     

基于特征理论的二维可压缩流动的二阶拉氏算法

提出一种求解二维拉氏可压缩流体力学方程的中心型二阶精度有限体积方法.利用特征理论构造网格节点处的局部近似演化算子,算子用来求解网格节点处的速度及压力,利用这些物理量更新节点位置及计算网格界面通量.通过结合一定的重构方案,该方法达到时、空二阶精度,并且形式简单、计算量小,适用于结构网格与非结构网格.典型数值实验表明,本文格式具有良好的收敛性、对称性及鲁棒性,且能自然地求解多物质流动问题.     

碱性溶液中Ag3O4的电化学性质及反应机理

首次研究高价银氧化物Ag3O4(可看作为由Ag(Ⅲ)和Ag(Ⅱ)*组成)在碱性水溶液中的电化学性质和反应机理.循环伏安和XRD测试表明,Ag3O4在碱性溶液中的电化学还原过程比较复杂:在较慢的放电条件下,Ag3O4中的Ag(Ⅲ)按照Ag(Ⅲ)→Ag(Ⅱ)→Ag(Ⅰ)→Ag反应途径逐步还原为单质银;在较快的放电条件下,Ag(Ⅲ)可以直接被还原为Ag(Ⅰ),即Ag(Ⅲ)→Ag(Ⅰ)→Ag.而Ag(Ⅱ)*可直接还原成金属Ag.Ag3O4的理论放电容量可以达到553.1 mAh/g,比通常锌-氧化银电池所用AgO的电容量高出27.8%.在119C放电倍率下,Ag3O4的放电容量依然达到理论容量的83%.显示了作为新型化学电源材料的应用前景.     

保持对称性的二维Lagrange流体力学有限差分方法

广泛应用的二维直角坐标系下的Wilkins有限差分格式在计算一维柱面问题时,通过等角度划分周向网格能够获得严格的对称性,非等角度划分周向网格会产生较严重的不对称性.通过分析Wilkins有限差分格式在处理非等角度划分周向网格的一维柱面问题时破坏对称性的原因,指出周向网格的非等角度划分产生了周向压力分量,从而产生了周向加速度分量和周向运动速度,以此为基础提出一种对该有限差分格式进行修正的方法,将节点处的周向压力分量做算术平均运算,以消除周向压力分量,只剩径向压力分量起作用.因而该修正方法在以任意角度划分周向网格的条件下都能够保持严格的对称性.通过几个典型算例验证该结论,对对称流动,修正方法与原始方法所获得的结果一致,对非对称流动,二者有微小差异.     

惯性-视觉复合传感器同步内触发设计

通过同步内触发惯性 视觉复合传感器,实现对远距离、高速动态目标的测量。将高精度惯性传感单元ADIS16355与普通单目摄像机同步触发采集,通过任意位置变化、光轴角度旋转,达到双目视觉效果,有效解决了三维空间运动目标定位问题。由激光按键控制图像起始帧的标注,同时产生高电平信号,与摄像机视频信号分离出的奇/偶场同步脉冲信号一起通过与门电路输出,作为单片机标记惯性数据的中断信号,实现同步测量,可有效消除错帧、丢帧等现象,精度可达微秒级,为惯性传感器与普通单目摄像机同步采集提供了一种有效的解决方案。     

铁-铍介质界面上激波折射现象的数值研究

用数值方法研究铁-铍介质界面上的激波折射现象.运用激波极曲线理论分析不同强度的激波从正规折射过渡到非正规折射的临界角变化.运用一个具有二阶精度和波传播性质的激波捕捉法,数值求解激波折射运动的流体力学方程组.对正规折射,数值结果与激波极曲线理论一致;对非正规折射,不同强度的激波大都存在前驱的折射激波,并且入射激波的强度不同、入射角度不同,激波折射的图像也不同.     

混合物状态方程的计算

多介质流体动力学过程的数值模拟往往涉及混合物状态方程的计算. 做图法和Newton 法是混合物状态方程计算常采用的方法, 前者虽直观精度却差, 后者计算效率高却只具有局部收敛性, 当解与其初始猜测值相差较远时Newton法不一定能够获得收敛解. 为此, 本文给出一种具有大范围收敛性的嵌入算法(imbedding method)求解混合物状态方程, 其基本思想是通过引入嵌入参数, 将待解的混合物状态方程和易解的混合物状态方程线性组合, 构成嵌入方程组, 当嵌入参数从0连续地变化到1 时, 嵌入方程组的解由易解的混合物状态方程的解连续地变化为待解的混合物状态方程的解. 嵌入方程组可由Newton法迭代求解, 也可转化为以嵌入参数为自变量的常微分方程组, 从而易于由成熟的计算方法如梯形法等进行求解. 进一步利用热力学基本关系, Maxwell形式的微分方程描述了压力和温度随嵌入参数的演化速率与应变速率和组分质量分数演化速率的关系. 对铅锡混合物热力学量的计算表明了本文算法的有效性.