火箭推进

为了获得电爆管输出压力及其冲击波传播规律,针对电爆管爆炸燃气建立了一维流动模型,并采用自适应二代小波配点法进行数值计算。分析了电爆管冲击波沿一维等截面和局部变截面通道的传播过程及规律。结果表明:自适应小波配点法结合人工黏性技术可以准确计算电爆管强冲击波传播问题; 电爆管爆炸冲击波沿一维等截面传火通道传播满足一定相似律,建立的相似律表达式能够估算同一类电爆管冲击波在爆炸近区的传播规律; 传火通道结构形式和结构参数变化对传播规律影响显著,针对同一类型的电爆管和变截面传火通道建立的相似律表达式,在1

In order to obtain the output pressure of cartridge and the propagation law of its shock wave, a one-dimensional flow model of the explosion gas was established, and the adaptive second-generation wavelet collocation method was applied in the numerical simulation.The propagation law of the shock wave along the one-dimensional equal section channel and the local variable section channel were analyzed. The results show that the adaptive wavelet collocation method combined with the artificial viscosity can precisely calculate the strong shock wave dashing out from the cartridge into the passage.The blast shock wave of the cartridge along the one-dimensional equal-section passage meets a certain similitude law, and the proposed similitude law formulas can estimate the propagation characteristics of the shock wave near the explosion field.The structural type and parameter of the passage have a significant effect on the propagation law.The similitude law expressions for the same type of cartridges and variable-section passage are proposed, and the estimated error is less than 10% under the condition of 1

引言
1 自适应小波配点法
2 电爆管传火通道模型
3 计算结果和分析
4 结论

图1 二代小波变换示意图<br/>Fig.1 Forward and inverse transformation of second-generation wavelet

图1 二代小波变换示意图
Fig.1 Forward and inverse transformation of second-generation wavelet

图2 传火通道横截面积变化示意图<br/>Fig.2 Cross-sectional area for different passage types

图2 传火通道横截面积变化示意图
Fig.2 Cross-sectional area for different passage types

图3 一维等截面传火通道模型<br/>Fig.3 One-dimensional physical model of the cartridge with an equal section passage

图3 一维等截面传火通道模型
Fig.3 One-dimensional physical model of the cartridge with an equal section passage

表1 电爆管模型初始条件<br/>Tab.1 Initial conditions of the cartridge model

表1 电爆管模型初始条件
Tab.1 Initial conditions of the cartridge model

图4 冲击波管的黎曼解与小波数值解<br/>Fig.4 Riemann and wavelet numerical solutions of the shock tube

图4 冲击波管的黎曼解与小波数值解
Fig.4 Riemann and wavelet numerical solutions of the shock tube

图5 等截面传火通道中压力变化和分布<br/>Fig.5 Pressure variation and pressure distribution in the equal-sectional passage

图5 等截面传火通道中压力变化和分布
Fig.5 Pressure variation and pressure distribution in the equal-sectional passage

图6 等截面传火通道不同位置的静压峰值<br/>Fig.6 Static pressure peaks at different locations of equal-section passage

图6 等截面传火通道不同位置的静压峰值
Fig.6 Static pressure peaks at different locations of equal-section passage

表2 等截面通道峰值压力估算验证<br/>Tab.2 Verification of the estimated peak pressures for equal-section passage

表2 等截面通道峰值压力估算验证
Tab.2 Verification of the estimated peak pressures for equal-section passage

图7 不同结构传火通道的压力特性<br/>Fig.7 Pressure characteristics of different passages

图7 不同结构传火通道的压力特性
Fig.7 Pressure characteristics of different passages

表3 局部拉瓦尔结构传火通道峰值压力估算验证<br/>Tab.3 Verification of the estimated peak pressures for the passages with a Laval structure

表3 局部拉瓦尔结构传火通道峰值压力估算验证
Tab.3 Verification of the estimated peak pressures for the passages with a Laval structure

表4 局部收缩或扩张传火通道峰值压力估算验证<br/>Tab.4 Verification of the estimated peak pressures for the passages with a convergent or divergent structure

表4 局部收缩或扩张传火通道峰值压力估算验证
Tab.4 Verification of the estimated peak pressures for the passages with a convergent or divergent structure

[1] 袁洪滨, 胡攀, 王志永, 等.浅析国外电爆阀电爆机理研究进展[J].火箭推进, 2016, 42(2):86-90.
[2] PAUL B H, GONTHIER K A,OKHUYSEN B.Analysis of gas dynamic waves in explosively actuated valves[J].Journal of Propulsion & Power, 2015, 26(3):479-496.
[3] DIBBERN A, SIPES W, HAGOPIAN M, et al.Implications of dynamic pressure transducer mounting variations on measurements in pyrotechnic test apparatus[C]//45th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit.Reston,Virigina:AIAA, 2009.
[4] REGELE J D, VASILYEV O V.An adaptive wavelet-collocation method for shock computations[J].International Journal of Computational Fluid Dynamics, 2009, 23(7):503-518.
[5] 唐玲艳.双曲型方程数值解的小波方法研究[D].长沙:国防科学技术大学, 2007:43-84.
[6] VASILYEV O V, BOWMAN C.Second-generation wavelet collocation method for the solution of partial differential equations[J].Journal of Computational Physics, 2000, 165(2):660-693.
[7] 王一博, 杨慧珠.方位各向异性介质的多尺度有限差分法波场模拟[J].清华大学学报(自然科学版), 2006, 46(2):293-296.
[8] KASSOY D, REGELE J, VASILYEV O.Detonation initiation on the microsecond time scale:one and two dimensional results obtained from adaptive wavelet-collocation numerical methods[C]//45th AIAA Aerospace Sciences Meeting and Exhibit.Reston, Virigina:AIAA, 2007.
[9] REGELE J D,KASSOY D R, VASILYEV O V.Numerical modeling of acoustic timescale detonation initiation[C]//46th 46th AIAA Aerospace Sciences Meeting and Exhibit. Reno, Nevasa:AIAA, 2008.
[10] REGELE J D,KASSOY D R,VEZOLAINEN A, et al.Purely gasdynamic multidimentioanal indirect detonation initiation using localized acoustic timescale power deposition[C]// 51st AIAA Aerospace Sciences Meeting Including the New Horizons Forum & Aerospace Exposition. Grapevine, Texas:AIAA, 2013.
[11] SCHEIDER K, VASILYEV O V.Wavelet Methods In Computational Fluid Dynamics[J].Annu.Rev.Fluid Mech.,2010(42):473-503.
[12] REGELE J D, KASSOY D R, ASLANI M, et al.Evolution of detonation formation initiated by a spatially distributed, transient energy source[J].Journal of Fluid Mechanics, 2016, 802:305-332.
[13] 许志宇, 谭永华, 李小明.基于自适应小波配点法的水击数值计算[J].火箭推进, 2019,45(5):32-37.
[14] 许志宇, 谭永华, 李小明.二维激波的小波多尺度数值计算[C]//中国航天第三专业信息网第40届技术交流会暨第四届空天动力联合会议论文集.昆明:中国航天第三专业信息网, 2019.

备注

引言

1 自适应小波配点法

2 电爆管传火通道模型

3 计算结果和分析

4 结论

期刊信息