排序方式: 共有38条查询结果,搜索用时 203 毫秒
11.
利用松弛近似,将非线性的凝聚炸药爆轰控制方程转化为线性的松弛方程组,并采用五阶WENO格式和五阶线性多步显隐格式对线性松弛方程组进行空间方向和时间方向的离散,由此建立具有高精度和高分辨率性质的计算凝聚炸药爆轰的松弛方法。建立的松弛方法可以避免求解Riemann问题及计算非线性通量的Jacobi矩阵,同时无需分裂处理反应源项。通过对凝聚炸药的平面一维定常爆轰波结构及球面一维聚心、散心爆轰起爆和传播过程的数值模拟,验证了所建立的松弛方法能够很好地计算凝聚炸药爆轰问题。 相似文献
12.
13.
为了对金属约束条件下的定常非理想爆轰进行理论研究,对未反应炸药和爆轰产物采用JWL形
式状态方程,对金属采用p(,T)形式状态方程。采用过爆轰前导冲击波的流线偏转角在影响域内沿高度线
性变化的假设,并且由未反应炸药和金属的冲击波极曲线的交点确定爆轰冲击边缘角,则可从未反应炸药斜
冲击波极曲线关系式求出爆轰前导冲击波阵面的形状。采用爆轰前导冲击波阵面之后的流线是直线的假设,
则爆轰流动控制方程由偏微分方程变为沿流线的常微分方程,沿着所有流线求解便给出爆轰化学反应区内
声速线与化学反应结束线的位置。理论分析同时给出约束金属折射冲击波后面的流体的流动状态。理论结
果给出的爆轰化学反应区结构特征和约束金属内的流动状态特征与高精度数值模拟的结果符合良好,说明
本文中给出的理论方法具有良好的合理性和适用性。 相似文献
14.
15.
16.
首先用改进冲击波极曲线理论分析惰性介质对钝感炸药爆轰的约束作用类型。改进冲击波极曲线基于爆轰ZND模型建立在前导冲击波上,并且未反应炸药采用JWL状态方程,惰性介质采用p(ρ,T)形式状态方程。理论考察声速小于炸药CJ爆速且压缩性不同的6种典型惰性介质约束情况。然后用带三项式Lee-Tarver化学反应率的二维Lagrange流体力学方法数值模拟考察约束相互作用。数值考察约束介质的影响因素有:压缩性、厚度、典型双层介质组合约束。从数值结果看出,由介质压缩性的不同给出的约束作用方式共7种:其中6种出现在介质声速小于炸药CJ爆速条件下,可运用冲击波极曲线理论;另外一种出现在介质声速大于炸药CJ爆速条件下,不能使用冲击波极曲线理论。同时,介质厚度、双层介质组合方式也能够影响爆轰前导冲击波阵面形状以及爆轰化学反应流动状态。 相似文献
17.
多介质流体动力学过程的数值模拟往往涉及混合物状态方程的计算. 做图法和Newton 法是混合物状态方程计算常采用的方法, 前者虽直观精度却差, 后者计算效率高却只具有局部收敛性, 当解与其初始猜测值相差较远时Newton法不一定能够获得收敛解. 为此, 本文给出一种具有大范围收敛性的嵌入算法(imbedding method)求解混合物状态方程, 其基本思想是通过引入嵌入参数, 将待解的混合物状态方程和易解的混合物状态方程线性组合, 构成嵌入方程组, 当嵌入参数从0连续地变化到1 时, 嵌入方程组的解由易解的混合物状态方程的解连续地变化为待解的混合物状态方程的解. 嵌入方程组可由Newton法迭代求解, 也可转化为以嵌入参数为自变量的常微分方程组, 从而易于由成熟的计算方法如梯形法等进行求解. 进一步利用热力学基本关系, Maxwell形式的微分方程描述了压力和温度随嵌入参数的演化速率与应变速率和组分质量分数演化速率的关系. 对铅锡混合物热力学量的计算表明了本文算法的有效性. 相似文献
18.
19.
再入飞行器湍流尾迹流场研究 总被引:1,自引:0,他引:1
再入飞行器湍流尾迹流场状况,直接关系到飞行器的雷达散射特性。对再入飞行器湍流尾迹等离子体场理论模型,试图通过湍流模式理论来表达,即使用κ-ε-g模型方程来封闭平均化的全Navier-Stokes方程,从而准确获得流动平均场和脉动场信息。使用的N-S平均方程由质量加权平均过程产生,湍流模型方也经过可压缩性修正。真实气体效应重点考察空气处于局部热化学平衡状态。流动控制方程运用一个二阶TVD格式的有限体积法求解,以一典型小钝锥体零攻角再入飞行为例,计算了在两种高程(H-40km和H=30km)条件下的高超声速湍流尾迹流场。获得的尾迹流场参数与流动物理状况符合,并且湍流脉动参数与已有相应的实验结果定性一致,初步证实该方法合理。 相似文献
20.