Numerical simulation for aluminum/air two-dimensional viscous two-phase detonation
-
摘要: 为了深入研究爆轰波形成和传播的机理与特性,建立了管内铝粉/空气二维黏性两相爆轰过程的数学模型,采用守恒元与求解元方法进行数值计算,并对其物理参数的分布进行了分析。结果表明:管内燃烧转爆轰的初期,压力沿径向变化明显,与壁面碰撞有明显的反射波;燃烧转爆轰的中后期,压力沿轴向变化明显,但径向效应仍不能忽视,碰撞形成的反射波对最终稳定爆轰波的形成影响较大。研究结果同时表明:不仅铝粉颗粒初始半径对爆轰波的形成与传播有一定的影响;而且气体的黏性作用在研究爆轰管内近壁面处流场时不容忽视。研究结果有利于进一步揭示铝粉燃烧转爆轰的机理。Abstract: A two-dimensional mathematical model of viscous aluminum/air two-phase detonation is established. The flow field inside the detonation tube was calculated by the method of conservation element and solution element. And the distribution of physical parameters is analyzed. Numerical results show that in the early stages of the deflagration-to-detonation transition process in the tube, pressure change obviously along the radius and significant reflection wave is collided off the wall. In the later period, pressure along the axis changes obviously, but the radial effect still can't be ignored. The reflection wave from the collision is important on the formation of steady detonation wave. Numerical results also show that the initial radius of the aluminum powder particles have certain influence on the formation and propagation of detonation wave. And certain influence of gas viscosity on flow field near the wall in the detonation tube was exerted. The results of the study can be utilized to reveal the mechanism of the deflagration-to-detonation transitim.
-
Key words:
- mechanics of explosion /
- two phase flow /
- CE/SE method /
- aluminum dust /
- detonation
-
铝粉含能量高, 并且在空气中易产生爆轰, 一方面在军事上可将其作为固体推进剂的重要成分, 另一方面在生产过程中易发生粉尘爆轰的危害也不容忽视。因此对铝粉/空气两相爆轰这个复杂而又不稳定过程的研究备受关注。E.L.Dreizin[1]、S.Goroshin等[2]对不同氧浓度热气体中的铝粉颗粒的燃烧进行了大量的实验研究。刘晓利等[3]、陈志华等[4]、李小东等[5]在铝粉粒径、铝粉浓度和点火延迟时间对爆轰参数的影响方面进行了实验研究。在理论研究方面, B.Veyssiere等[6]和A.V.Fedorov等[7]用两相流体力学模型研究铝粉爆轰的定常问题。韦伟等[8-9]和洪滔等[10-11]对铝粉爆轰的一维和二维非定常问题进行了研究。
在前人研究的基础上, 本文中考虑气固两相爆轰中气体的黏性作用, 并以此建立管内非定常二维气固两相黏性爆轰的理论模型, 应用CE/SE方法数值求解气固两相爆轰模型。分析铝粉/空气两相爆轰波管内形成和传播的过程, 研究铝粉颗粒半径对爆轰波形成和强度的影响, 讨论考虑气体黏性作用的重要性。
1. 爆轰两相控制方程
为了简化所研究的管内爆轰气/固两相非定常二维轴对称流场, 作以下假定:(1)铝粉颗粒为初始直径相同的球形, 在气体中均匀分布可视其为拟流体, 但忽略其颗粒间的相互作用力, 以及与壁面的作用; (2)考虑气相的黏性作用, 忽略颗粒相的黏性作用; (3)燃烧产物氧化铝粉末视作气相产物, 但不考虑其对压力的影响; (4)化学反应释放的能量都被气相吸收。基于以上假设, 可得到含组分变化的两相爆轰控制方程组:
\frac{\partial \boldsymbol{U}}{\partial t}+\frac{\partial \boldsymbol{F}}{\partial x}+\frac{\partial \boldsymbol{G}}{\partial y}=\frac{\partial \boldsymbol{F}_{\mathrm{v}}}{\partial x}+\frac{\partial \boldsymbol{G}_{\mathrm{v}}}{\partial y}+\boldsymbol{R}-\frac{\boldsymbol{G}-\boldsymbol{G}_{\mathrm{v}}}{y} (1) \boldsymbol{U}=\left[\begin{array}{c} \varphi_{\mathrm{g}} \rho_{\mathrm{g}} \\ \varphi_{\mathrm{s}} \rho_{\mathrm{s}} \\ \varphi_{\mathrm{g}} \rho_{\mathrm{g}} u_{\mathrm{g}} \\ \varphi_{\mathrm{g}} \rho_{\mathrm{g}} v_{\mathrm{g}} \\ \varphi_{\mathrm{s}} \rho_{\mathrm{s}} u_{\mathrm{s}} \\ \varphi_{\mathrm{s}} \rho_{\mathrm{s}} v_{\mathrm{s}} \\ \varphi_{\mathrm{g}} \rho_{\mathrm{g}} E_{\mathrm{g}} \\ \varphi_{\mathrm{s}} \rho_{\mathrm{s}} E_{\mathrm{s}} \\ \varphi_{\mathrm{g}} \rho_{\mathrm{g}} Y_{0} \\ \varphi_{\mathrm{g}} \rho_{\mathrm{g}} Y_{\mathrm{Al}} \end{array}\right], \boldsymbol{F}=\left[\begin{array}{c} \varphi_{\mathrm{g}} \rho_{\mathrm{g}} u_{\mathrm{g}} \\ \varphi_{\mathrm{s}} \rho_{\mathrm{s}} u_{\mathrm{s}} \\ \varphi_{\mathrm{g}} \rho_{\mathrm{g}} u_{\mathrm{g}}^{2} \\ \varphi_{\mathrm{g}} \rho_{\mathrm{g}} u_{\mathrm{g}} v_{\mathrm{g}} \\ \varphi_{\mathrm{s}} \rho_{\mathrm{s}} u_{\mathrm{s}}^{2} \\ \varphi_{\mathrm{s}} \rho_{\mathrm{s}} u_{\mathrm{s}} v_{\mathrm{s}} \\ \varphi_{\mathrm{g}} u_{\mathrm{g}}\left(\rho_{\mathrm{g}} E_{\mathrm{g}}+p\right) \\ \varphi_{\mathrm{s}} u_{\mathrm{s}}\left(\rho_{\mathrm{s}} E_{\mathrm{s}}+p\right) \\ \varphi_{\mathrm{g}} \rho_{\mathrm{g}} u_{\mathrm{g}} Y_{0} \\ \varphi_{\mathrm{g}} \rho_{\mathrm{g}} u_{\mathrm{g}} Y_{\mathrm{Al}} \end{array}\right], \boldsymbol{G}=\left[\begin{array}{c} \varphi_{\mathrm{g}} \rho_{\mathrm{g}} v_{\mathrm{g}} \\ \varphi_{\mathrm{s}} \rho_{\mathrm{s}} v_{\mathrm{s}} \\ \varphi_{\mathrm{g}} \rho_{\mathrm{g}} u_{\mathrm{g}} v_{\mathrm{g}} \\ \varphi_{\mathrm{g}} \rho_{\mathrm{g}} v_{\mathrm{g}}^{2} \\ \varphi_{\mathrm{s}} \rho_{\mathrm{s}} u_{\mathrm{s}} v_{\mathrm{s}} \\ \varphi_{\mathrm{s}} \rho_{\mathrm{s}} v_{\mathrm{s}}^{2} \\ \varphi_{\mathrm{g}} v_{\mathrm{g}}\left(\rho_{g} E_{\mathrm{g}}+p\right) \\ \varphi_{\mathrm{s}} v_{\mathrm{s}}\left(\rho_{s} E_{\mathrm{s}}+p\right) \\ \varphi_{\mathrm{g}} \rho_{\mathrm{g}} v_{\mathrm{g}} Y_{0} \\ \varphi_{\mathrm{g}} \rho_{\mathrm{g}} v_{\mathrm{g}} Y_{\mathrm{Al}} \end{array}\right], \boldsymbol{F}_{\mathrm{v}}=\left[\begin{array}{c} 0 \\ 0 \\ \varphi_{\mathrm{g}} \tau_{\mathrm{g}, x x} \\ \varphi_{\mathrm{g}} \tau_{\mathrm{g}, x y} \\ \varphi_{\mathrm{s}} \tau_{\mathrm{s}, x x} \\ \varphi_{\mathrm{s}} \tau_{\mathrm{s}, x y} \\ \varphi_{\mathrm{g}}\left(u_{\mathrm{g}} \tau_{\mathrm{g}, x x}+v_{\mathrm{g}} \tau_{\mathrm{g}, x y}+r_{\mathrm{c}} \frac{\partial T_{\mathrm{g}}}{\partial x}\right) \\ \varphi_{\mathrm{s}}\left(u_{\mathrm{s}} \tau_{\mathrm{s}, x x}+v_{\mathrm{s}} \tau_{\mathrm{s}, x y}+r_{\mathrm{c}} \frac{\partial T_{\mathrm{s}}}{\partial x}\right) \\ 0 \\ 0 \end{array}\right] \boldsymbol{G}_{\mathrm{v}}=\left[\begin{array}{c} 0 \\ 0 \\ \varphi_{\mathrm{g}} \tau_{\mathrm{g}, y x} \\ \varphi_{\mathrm{g}} \tau_{\mathrm{g}, y y} \\ \varphi_{\mathrm{s}} \tau_{\mathrm{s}, y x} \\ \varphi_{\mathrm{s}} \tau_{\mathrm{s}, y y} \\ \varphi_{\mathrm{g}}\left(u_{\mathrm{g}} \tau_{\mathrm{g}, y x}+v_{\mathrm{g}} \tau_{\mathrm{g}, y y}+r_{\mathrm{c}} \frac{\partial T_{\mathrm{g}}}{\partial y}\right) \\ \varphi_{\mathrm{s}}\left(u_{\mathrm{s}} \tau_{\mathrm{s}, x x}+v_{\mathrm{s}} \tau_{\mathrm{s}, y y}+r_{\mathrm{c}} \frac{\partial T_{\mathrm{s}}}{\partial y}\right) \\ 0 \\ 0 \end{array}\right] \boldsymbol{R}=\left[\begin{array}{c} I_{\mathrm{d}} \\ -I_{\mathrm{d}} \\ I_{\mathrm{d}} u_{\mathrm{g}}-F_{\mathrm{d}, x} \\ I_{\mathrm{d}} v_{\mathrm{g}}-F_{\mathrm{d}, y}+\varphi_{\mathrm{g}}\left(p-\tau_{\mathrm{g}, \theta \theta}\right) / y \\ -I_{\mathrm{d}} u_{\mathrm{s}}+F_{\mathrm{d}, x} \\ -I_{\mathrm{d}} v_{\mathrm{s}}+F_{\mathrm{d}, y}+\varphi_{\mathrm{s}}\left(p-\tau_{\mathrm{s}, \theta \theta}\right) / y \\ -Q_{\mathrm{d}}+Q_{\mathrm{c}}-\left(F_{\mathrm{d}, x} u_{\mathrm{g}}+F_{\mathrm{d}, y} v_{\mathrm{g}}\right)+I_{\mathrm{d}}\left(E_{\mathrm{g}}+p / \rho_{\mathrm{g}}\right) \\ Q_{\mathrm{d}}+\left(F_{\mathrm{d}, x} u_{\mathrm{s}}+F_{\mathrm{d}, y} v_{\mathrm{s}}\right)-I_{\mathrm{d}}\left(E_{\mathrm{s}}+p / \rho_{\mathrm{s}}\right) \\ -\frac{8}{9} I_{\mathrm{d}} \\ \frac{17}{9} I_{\mathrm{d}} \end{array}\right] 式中:
\tau_{\mathrm{g}, x x}=\mu_{\mathrm{g}}\left(2 \frac{\partial u_{\mathrm{g}}}{\partial x}-\frac{2}{3} \nabla u_{\mathrm{g}}\right), \tau_{\mathrm{s}, x x}=\mu_{\mathrm{s}}\left(2 \frac{\partial u_{\mathrm{s}}}{\partial x}-\frac{2}{3} \nabla u_{\mathrm{s}}\right), \tau_{\mathrm{g}, y y}=\mu_{\mathrm{g}}\left(2 \frac{\partial u_{\mathrm{g}}}{\partial x}-\frac{2}{3} \nabla u_{\mathrm{g}}\right) \tau_{\mathrm{s}, y y}=\mu_{\mathrm{s}}\left(2 \frac{\partial u_{\mathrm{s}}}{\partial x}-\frac{2}{3} \nabla u_{\mathrm{s}}\right) \tau_{\mathrm{g}, x y}=\tau_{g, y x}=2 \mu_{g}\left(\frac{\partial u_{g}}{\partial y}+\frac{\partial v_{g}}{\partial x}\right), \tau_{\mathrm{s}, x y}=\tau_{\mathrm{s}, y x}=2 \mu_{\mathrm{s}}\left(\frac{\partial u_{\mathrm{s}}}{\partial y}+\frac{\partial v_{\mathrm{s}}}{\partial x}\right) \tau_{\mathrm{g}, \theta \theta}=\mu_{\mathrm{g}}\left(2 \frac{v_{\mathrm{g}}}{y}-\frac{2}{3} \nabla u_{\mathrm{g}}\right), \tau_{\mathrm{s}, \theta \theta}=\mu_{\mathrm{s}}\left(2 \frac{v_{\mathrm{s}}}{y}-\frac{2}{3} \nabla u_{\mathrm{s}}\right), \nabla u_{\mathrm{g}}=\frac{\partial u_{\mathrm{g}}}{\partial x}+\frac{\partial\left(y v_{\mathrm{g}}\right)}{y \partial y}, \nabla u_{\mathrm{s}}=\frac{\partial u_{\mathrm{s}}}{\partial x}+\frac{\partial\left(y v_{\mathrm{s}}\right)}{y \partial y} 下标g、s表示气相和固相; φ、ρ、u、v、μ、T、E分别表示体积分数比、密度、轴向速度、径向速度、黏性系数、温度和内能, 其中φg+φs=1;p为爆轰管内的压力; YO、YAl分别为气相部分中O2和Al2O3的体积分数; Id是单位体积内铝粉颗粒因燃烧引起的质量变化率; Fd是气相与固体颗粒相间的相互作用力; Qd是气相与固体颗粒间的对流传热量; Qc是铝粉燃烧释放的热量; τij(i, j=x, y, θ), 当i=j时表示正应力, 当i≠j时表示剪切应力, x、y分别表示爆轰管的轴向和径向。
2. 数值计算方法
2.1 CE/SE方法
通过之前的研究, 发现二维CE/SE方法在扑捉爆轰波这类强间断波的能力很强[8-9]。本文中尝试采用该方法处理黏性的问题。
对于守恒方程组(1), 令
, 根据Chang[14]对时间和空间的统一处理思想, 令x1=x, x2=y, x3=t作为三维欧几里得空间E3上的3个坐标, 利用高斯定理对其积分可得[14]:
\oint_{S(V)} \boldsymbol{h} \mathrm{d} s-\oint_{S(V)} \boldsymbol{h}_{\mathrm{V}} \mathrm{d} s=\int_{V} \boldsymbol{I} \mathrm{d} V (2) 式中:S(V)是V在E3上任意时空区域的边界, h=(F, G, U)是时间-空间通量流密度矢量。对于任意的(x, y, t)∈SE(i, j, n), U(x, y, t)、F(x, y, t)、G(x, y, t)用离散量U*(x, y, t; i, j, n)、F*(x, y, t; i, j, n)、G*(x, y, t; i, j, n)来替代。利用泰勒级数展开式见参考文献[14]。由此方程(12)在守恒元CE(Q)上用离散量表示为[14]:
\oint_{S(\mathrm{CE}(Q))} \boldsymbol{h}^{*} \mathrm{~d} s-\oint_{S(\mathrm{CE}(Q))} \boldsymbol{h}_{\mathrm{V}} \mathrm{d} s=\int_{\mathrm{CE}(Q)} \boldsymbol{I} \mathrm{d} v (3) 由上式表明穿过
所有边界面的h*的积分通量是守恒的, 因此利用守恒元和求解元思想可将
进一步推导可得到数值计算格式, 具体参见文献[14]。
2.2 源项处理
由于化学反应的特征时间远小于对流的特征时间, 因此源项为刚性。本文中采用分裂方法求解, 即先不考虑源项的影响, 用CE/SE方法求解(U)jn。再将(U)jn作为初值, 应用四阶Runge-Kutta法求解常微分方程组
\frac{\mathrm{d} \boldsymbol{U}}{\mathrm{d} t}=\boldsymbol{R}-\frac{\boldsymbol{G}-\boldsymbol{G}_{\mathrm{v}}}{y} 2.3 初始条件与边界条件
计算初始条件为:根据实验设备的实际尺寸, 爆轰管直径取0.06 m, 长取1 m。因流场轴对称, 因此只取流场的一半为计算区域。初始时爆轰管内充满空气和固体颗粒, 铝粉和空气的初始质量按其化学当量比取值。初始压力取一个标准大气压, 空气和固体颗粒的温度取常温300 K。铝粉颗粒半径取0.5 μm。爆轰管起爆端轴线处取10×15个网格模拟初始点火, 点火条件为常温常压的15倍, 压力取1.5 MPa, 温度取4 500 K。
计算边界条件为:计算中心轴上采用轴向对称边界条件, 壁面上采用无滑移边界条件[6]。
3. 结果分析
3.1 管内压力分布
通过计算得到爆轰管内不同时刻压力空间分布云图, 如图 1所示。图 1(a)为铝粉在空气中被点燃后产生的高温高压燃烧产物压缩未燃混合物, 形成的燃烧波在爆轰管内沿轴向与径向传播, 压力峰值基本不变(p=1.1 MPa)。图 1(b)显示燃烧波传播过程中高压区域不断扩大, 当沿径向运动的燃烧波与下壁面碰撞时, 形成一个局部爆轰点和向上壁面运动的反射波, 其压力峰值上升到2.2 MPa。局部爆轰点周围的铝粉以及反射波扫过区域的未反应物被迅速点燃, 铝粉燃烧释放热量速率加快, 燃烧波在传播过程中能量不断增加。图 1(c)显示爆轰点高压范围不断扩大, 并向上壁面运动, 与上壁面碰撞反射, 逐渐形成一道明显的爆轰波, 传播至距爆轰管封闭端0.19 m处。图 1(d)中爆轰波基本达到稳定状态, 此处距左边封闭端0.546 m, 压力峰值达到2.8 MPa。
通过以上对压力云图的分析可见, 铝粉尘在管内燃烧转爆轰的初期, 压力效应以径向为主, 而随着爆轰波逐渐形成, 压力轴向效应已起主导作用, 而径向效应明显减弱, 这与实际的铝粉尘燃烧转爆轰情况相符。铝粉在空气中被点燃后形成的燃烧波上下壁面的反射作用较管内液体燃料的爆轰过程[13]明显。该作用加快了铝粉燃烧速率, 对铝粉/空气爆轰过程中局部爆轰点的产生和最终稳定爆轰波形成有较大的关系。
3.2 铝粉颗粒半径变化对流场影响
图 2为不同时刻爆轰管内铝粉颗粒半径空间分布曲面图。图中铝粉颗粒半径d为0.5 μm, 表示铝粉尚未开始燃烧, 等于0表示燃烧完全。图 2(a)为16 μs时刻铝粉颗粒燃烧情况。此时刚开始点火不久, 燃烧波后的铝粉没有燃烧完全, d > 0。因此爆轰能量不够, 稳定的爆轰波尚未形成, 见图 1(b)。图 2(b)为44 μs时刻铝粉颗粒燃烧情况。此时爆轰波扫过后的铝粉在高温高压环境下迅速燃烧完毕, 颗粒半径由0.5 μm变化到0, 已经形成爆轰波, 见图 1(d)。
通过计算表明, 不同的颗粒半径对稳定爆轰波形成有一定的影响。比较颗粒半径分别为0.4、0.5、0.6 μm的初始起爆条件, 如图 3所示。颗粒半径越小, 越早达到稳定爆轰, 爆轰波峰值较低; 颗粒半径越大, 越迟达到稳定爆轰, 爆轰波峰值较高。但颗粒半径达到一定大小时对爆轰波的影响减弱。
3.3 黏度对速度及压力的影响
通过比较分析爆轰波阵面附近速度云图, 从图 4可见, 爆轰波传播的速度因气体黏性的作用而明显下降, 由1 300 m/s下降至1 100 m/s。受黏性力影响, 轴线处的轴向速度和径向速度的峰值速度更为平稳。其中轴线处轴向速度是指在爆轰管的轴线处气相的轴向速度, 轴线处径向速度是指在爆轰管的轴线处气相的径向速度。
受气体黏性的影响, 管壁处轴向速度u整体明显下降, 如图 5(a)~(b)所示; 而径向速度v整体明显上升, 如图 5(c)~(d)所示。其中管壁处轴向速度是指接近爆轰管管壁处气相的轴向速度, 管壁处径向速度是指接近爆轰管管壁处气相的径向速度。可见速度受黏性的影响在爆轰管壁面处更为显著, 如图 4(b)所示。
比较同一时刻不考虑气体黏性作用的爆轰管内压力分布云图(图 6)与考虑气体黏性作用的爆轰管内压力分布云图(图 1(d))可见。同样是在t=44 μs时刻, 不考虑气体黏性作用时爆轰波已传至距封闭端0.66 m处, 压力峰值可达3.4 MPa, 考虑气体黏性作用时爆轰波仅传至距封闭端0.54 m处, 压力峰值下降为2.8 MPa。证实考虑气体黏性作用, 使得爆轰波在传播过程中阻力增加, 传播速度下降, 而且峰值压力也受到影响。可见由于黏性的存在, 使得一部分的能量损失, 因此爆轰波强度有所减弱。
4. 结论
通过对铝粉/空气二维黏性两相管内爆轰过程分析, 建立了合理的数学模型, 不仅考虑了气相的组分, 更考虑了气体黏性因素。采用守恒元与求解元方法对其进行数值计算, 对流场各物理量进行了较为详细的研究。结果表明:
(1) 二维黏性CE/SE方法在捕捉铝粉/空气两相爆轰波方面的能力较强。
(2) 管内燃烧转爆轰初期, 压力沿径向变化明显, 与壁面碰撞有明显的反射波; 燃烧转爆轰中后期, 压力沿轴向变化明显, 但径向效应仍不能忽视, 碰撞形成的反射波对最终稳定爆轰波的形成影响较大。
(3) 初始铝粉颗粒半径对爆轰形成与传播有一定的影响。
(4) 气体的黏性作用使得爆轰波峰值速度更平稳, 对近爆轰管内壁面处速度的影响较轴线处明显。受黏性作用影响, 管壁处径向速度显著增加, 而轴向速度明显下降。壁面处呈现明显的黏性效应。气体黏性作用使得管内爆轰波传播速度变慢, 而且峰值压力下降。
-
-
[1] Dreizin E L. Experimental study of aluminum particle flame evolution in normal and micro-gravity[J]. Combustion and Flame, 1999, 116(3): 323-333. http://www.sciencedirect.com/science/article/pii/S0010218097003313 [2] Goroshin S, Bidabadi M, Lee J H S. Quenching distance of laminar flame in aluminum dust cIouds[J]. Combustion and Flame, 1996, 105: 147-160. doi: 10.1016/0010-2180(95)00183-2 [3] 刘晓利, 李鸿志, 郭建国, 等.铝粉-空气混合物燃烧转爆轰(DDT)过程的实验研究[J].爆炸与冲击, 1995, 15(3): 217-228. http://www.cqvip.com/Main/Detail.aspx?id=1693987Liu Xiao-li, Li Hong-zhi, Guo Jian-guo, et al. Deflagration to detonatiion transition(DDT)in aluminum dust-air mixture[J]. Explosion and Shock Waves, 1995, 15(3): 217-228. http://www.cqvip.com/Main/Detail.aspx?id=1693987 [4] 陈志华, 范宝春, 李鸿志.燃烧管内悬浮铝粉燃烧爆炸过程的研究[J].高压物理学报, 2006, 20(2): 157-162. http://www.cnki.com.cn/Article/CJFDTotal-GYWL200602007.htmChen Zhi-hua, Fan Bao-chun, Li Hong-zhi. Investigations on combustion and explosion of Suspended aluminum particles in a large combustion tube[J]. Journal of High Pressure Physics, 2006, 20(2): 157-162. http://www.cnki.com.cn/Article/CJFDTotal-GYWL200602007.htm [5] 李小东, 刘庆明, 白春华, 等.铝粉-空气混合物的燃烧转爆轰过程[J].火炸药学报, 2009, 32(6): 58-61.Li Xiao-dong, Liu Qing-ming, Bai Chun-hua, et al. Deflagration to detonation transition process in aluminum dust-air mixture[J]. Chinese Journal of Explosives & Propellants, 2009, 32(6): 58-61. [6] Veyssiere B, Khasainov B A. Stead, Plane, double-front detonations in gaseous detonable mixtures containing a suspension of aluminum particles[J]. Combustion and Flame, 1991, 85: 241-253. doi: 10.2514/5.9781600865886.0284.0299 [7] Fedorov A V, Khmel T A, Fomin V M. Non-equilibrium model of steady detonations in aluminum particles-oxygen suspensions[J]. Shock Waves, 1999, 9(5): 3l3-318. doi: 10.1007/s001930050191 [8] 韦伟, 翁春生.基于CE/SE方法的铝粉尘爆轰一维两相数值计算[J].南京师范大学学报, 2012, 12(2): 53-56. http://d.wanfangdata.com.cn/Periodical/ddxb201204022Wei Wei, Weng Chun-sheng. One-dimension two-phase-flow numerical simulation of aluminum-dust detonation based on CE/SE method[J]. Journal of Nanjing Normal University, 2012, 12(2): 53-56. http://d.wanfangdata.com.cn/Periodical/ddxb201204022 [9] 韦伟, 翁春生.基于CE/SE方法的铝粉尘爆轰二维两相数值计算[J].弹道学报, 2012, 24(4): 99-102. http://d.wanfangdata.com.cn/Periodical/ddxb201204022Wei Wei, Weng Chun-sheng. Two-dimension two-phase-flow numerical simulation of aluminum-dust detonation based on CE/SE method[J]. Journal of Ballistics, 2012, 24(4): 99-102. http://d.wanfangdata.com.cn/Periodical/ddxb201204022 [10] 洪涛, 秦承森.爆轰波管中铝粉尘爆轰的数值模拟[J].爆炸与冲击, 2004, 24(3): 193-200. http://www.cnki.com.cn/article/cjfdtotal-bzcj200403001.htmHong tao, Qin Cheng-sen. Numerical simulation of dust detonation of aluminum powder in explosive tubes[J]. Explosion and Shock Waves, 2004, 24(3): 193-200. http://www.cnki.com.cn/article/cjfdtotal-bzcj200403001.htm [11] 洪涛, 秦承森.悬浮铝粉尘爆轰波参数[J].含能材料, 2004, 12(3): 129-133. http://d.wanfangdata.com.cn/Periodical/hncl200403001Hong tao, Qin Cheng-sen. Parameters of detonation in suspended aluminum dust[J]. Chinese Journal of Energetic Materials, 2004, 12(3): 129-133. http://d.wanfangdata.com.cn/Periodical/hncl200403001 [12] 白桥栋, 翁春生.二维CE/SE方法在内弹道湍流两相流动中的应用[J].兵工学报, 2009, 30(6): 682-687. http://d.wanfangdata.com.cn/Periodical/bgxb200906006Bai Qiao-dong, Weng Chun-sheng. The application of two dimensional space-time conservation element and solution element method(CE/SE)to turbulent two-phase flow in interior ballistics[J]. Acta Armamentarii, 2009, 30(6): 682-687. http://d.wanfangdata.com.cn/Periodical/bgxb200906006 [13] Ma Dan-hua, Weng Chun-sheng. Application of two-dimensional viscous CE/SE method in calculation of two-phase detonation[J]. Journal of China Ordnance, 2010, 6(1): 5-9. http://www.cnki.com.cn/Article/CJFDTotal-BAXY201001003.htm [14] 翁春生, 王浩.计算内弹道[M].北京: 国防工业出版社, 2006: 276-290. 期刊类型引用(7)
1. 苏晓杰,续晗,翁春生,倪晓冬,肖博文,郑权. 粉末爆轰发动机的粉末燃料供应特性研究. 推进技术. 2022(10): 418-426 . 百度学术
2. 刘国秋,孙宇航,高明,桂柯. 铝粉燃料火箭发动机的爆轰特性研究. 科技风. 2021(16): 175-176 . 百度学术
3. 刘龙,夏智勋,黄利亚,?马立坤,那旭东. 镁颗粒-空气混合物一维稳态爆震波特性数值模拟. 物理学报. 2019(24): 216-226 . 百度学术
4. 刘龙,夏智勋,黄利亚. 气相氛围中悬浮粉末燃料爆震燃烧研究进展. 宇航学报. 2018(03): 239-248 . 百度学术
5. 陈曦,陈先锋,张洪铭,刘晅亚,张英,牛奕,胡东涛. 惰化剂粒径对铝粉火焰传播特性影响的实验研究. 爆炸与冲击. 2017(04): 759-765 . 本站查看
6. 韦伟,翁春生. 以铝粉为燃料的脉冲爆轰发动机数值研究. 固体火箭技术. 2017(01): 37-40+59 . 百度学术
7. 韦伟,翁春生. 基于三维两相CE/SE方法的点火位置对固体燃料PDE的影响研究. 弹道学报. 2016(03): 65-70 . 百度学术
其他类型引用(4)
-