Deformation,fragmentation and acceleration of a controlled fragmentation charge casing

Two different finite element software, LS-DYNA and Impetus, have been evaluated to test their ability to predict the deformation, fragmentation and acceleration of a controlled fragmentation charge casing. The general-purpose program LS-DYNA was used with a multi-material ALE formulation and a mass-preserving erosion criterion coupled to a Johnson-Cook fracture criterion. In the Impetus simulations, a third order Lagrangian element formulation was used for the casing and a node-splitting element erosion treatment coupled to a Cockcroft-Latham failure criterion was used to describe casing fracture. The high-explosive gases were described by a discrete particle formalism.In order to acquire data to validate our computational tools and constitutive models, a series of experiments have been performed using a laboratory charge with an internal grooved casing. In the test series, the charge geometry was fixed except that the groove depth were varied from very shallow to very deep resulting in different deformation patterns, fracture modes and terminal velocities. Various diagnostic tools captured the different stages of the expansion and fragmentation of the casing. A high-speed framing camera depicted the deformation pattern before fragmentation and was used to determine the moment when the casing failed. Three different complementary techniques were used to follow the acceleration of the fragments; a Photon Doppler velocimetry to determine the initial acceleration of the casing, double exposed radiographs to estimate the fragment velocity after break-up and a high-speed video to determine the terminal velocity of the fragment after leaving the fireball. In addition, the fragments were soft recovered in a set of sawdust pit tests and their final shape and weight were measured. A SEM was used to characterise the fracture surfaces and to determine the modus of fracture (tensile or shear failure).Comparisons to experiments show that both software can predict the change in deformation behaviour when the groove depth increases, from tangential necking for shallow grooves to radial punching for deep groves. Both software could also reasonable well predict the acceleration of the fragments, though both overestimates the terminal velocity for the charge with the deepest grooves.     

Fragment spatial distribution of prismatic casing under internal explosive loading

Non-cylindrical casings filled with explosives have undergone rapid development in warhead design and explosion control. The fragment spatial distribution of prismatic casings is more complex than that of traditional cylindrical casings. In this study, numerical and experimental investigations into the fragment spatial distribution of a prismatic casing were conducted. A new numerical method, which adds the Lagrangian marker points to the Eulerian grid, was proposed to track the multi-material interfaces and material dynamic fractures. Physical quantity mappings between the Lagrangian marker points and Eulerian grid were achieved by their topological relationship. Thereafter, the fragment spatial distributions of the prismatic casing with different fragment sizes, fragment shapes, and casing geometries were obtained using the numerical method. Moreover, fragment spatial distribution experiments were conducted on the prismatic casing with different fragment sizes and shapes, and the experimental data were compared with the numerical results. The effects of the fragment and casing geometry on the fragment spatial distributions were determined by analyzing the numerical results and experimental data. Finally, a formula including the casing geometry parameters was fitted to predict the fragment spatial distribution of the prismatic casing under internal explosive loading.     

在分离式Hopkinson拉杆实验中确定临界实验条件的经验方法

在分离式Hopkinson拉杆实验中,常常需要获得试样在经历一个拉力脉冲拉断情况下的最低应变率,而在分离式Hopkinson拉杆装置中确定对应最低应变率的发射压力需要进行多次尝试性实验。为了更好地解决这一问题,提出一种简单的经验方法,可以通过一次拉断的实验结果计算出一个拉力脉冲拉断试样需要的最小发射压力及相应的应变率,进一步可以预测不同发射压力下的应变率。对4种铝合金试样进行拉伸实验,结果表明,利用提出的经验方法得到的应变率结果与实验结果吻合较好,可以有效地对发射压力和实验应变率进行预估。     

Simulation of natural fragmentation of rings cut from warheads

Natural fragmentation of warheads that detonates causes the casing of the warhead to split into various sized fragments through shear or radial fractures depending on the toughness,density,and grain size of the material.The best known formula for the prediction of the size distribution is the Mott formulae,which is further examined by Grady and Kipp by investigating more carefully the statistical most random way of portioning a given area into a number of entities.We examine the fragmentation behavior of radially expanding steel rings cut from a 25 mm warhead by using an in house smooth particle hydrodynamic(SPH) simulation code called REGULUS.Experimental results were compared with numerical results applying varying particle size and stochastic fracture strain.The numerically obtained number of fragments was consistent with experimental results.Increasing expansion velocity of the rings increases the number of fragments.Statistical variation of the material parameters influences the fragment characteristics,especially for low expansion velocities.A least square regression fit to the cumulative number of fragments by applying a generalized Mott distribution shows that the shape parameter is around 4 for the rings,which is in contrast to the Mott distribution with a shape parameter of 1/2.For initially polar distributed particles,we see signs of a bimodal cumulative fragment distribution.Adding statistical variation in material parameters of the fracture model causes the velocity numerical solutions to become less sensitive to changes in resolution for Cartesian distributed particles.     

Fracture mechanism of a cylindrical shell cut by circumferential detonation collision

The failure mechanism of a cylindrical shell cut into fragments by circumferential detonation collision was experimentally and numerically investigated. A self-designed detonation wave regulator was used to control the detonation and cut the shell. It was found that the self-designed regulator controlled the fragment shape. The macrostructure and micro-characteristics of fragments revealed that shear fracture was a prior mechanism, the shell fractured not only at the position of detonation collision, but the crack also penetrated the shell at the first contact position of the Chapmen-Jouguet (C-J) wave. The effects of groove number and outer layer thickness on the fracture behavior were tested by simulations. When the thickness of the outer layer was 5–18 mm, it has little effect on fragmentation of the shell, and shells all fractured at similar positions. The increase of the groove number reduced the fracture possibility of the first contact position of the C-J wave. When the groove number reached 7 with a 10 mm outer layer (1/4 model), the fracture only occurred at the position of detonation collision and the fragment width rebounded.     

稀薄气体动理论中介观尺度数值模拟的加速方法

针对滑移流、早期过渡流区域采用离散速度法(discrete velocity method, DVM)求解Boltzmann方程收敛速度极慢、计算资源消耗大的难题,提出全流场耦合介观/宏观方程加速方法。在介观层面基于有限差分的DVM求解Boltzmann方程,在宏观层面基于有限体积的压力耦合方程组半隐式方法求解矩方程,充分利用纳维-斯托克斯-傅里叶/R26矩方程在低克努森数下的快速收敛特性,对介观方程进行加速。基于高阶Hermite多项式重构分布函数,完成宏观方程与介观方程之间的数据传输。仿真结果表明:全流场耦合介观/宏观方程的加速方法在滑移流、早期过渡流区域具有显著加速性能,最大降低了95.28%的计算时长;但对中、大克努森数流域,加速性能大幅度下降。     

平流层长航时气球上升过程超冷现象影响因素分析

平流层长航时气球超冷现象是指内部浮升气体温度低于外界大气温度,超冷会引起浮力损失进而阻碍气球上升过程。建立了长航时气球辐射、对流等热模型和上升过程动力学模型,仿真分析了初始净浮力、蒙皮热物性参数、放飞时间和放飞日期等因素对超冷现象的影响规律。研究结果表明,上升过程内外温差随初始净浮力增大而增大且变化显著,可见光吸收率和红外吸收率增大时,内外温差值总体上呈减小趋势;放飞时刻和放飞日期对超冷现象影响较小,但放飞时间不同,气球上升至设计驻空高度的时间差别较大。该结论可为平流层浮空器总体方案设计和放飞试验提供有益参考。     

Accident Scenarios Involving Pebble Bed High Temperature Reactors

Proponents of high temperature gas cooled reactors argue that the reactor type is inherently safe and that severe accidents with core damage and radioactive releases cannot occur. The argument is primarily based on the safety features of the special form of the fuel. This paper examines some of the assumptions underlying the safety case for high temperature gas cooled reactors and highlights ways in which there could be fuel failure even during normal operations of the reactor; these failures serve to create a radioactive inventory that could be released under accident conditions. It then describes the severe accident scenarios that are the greatest challenge to high temperature gas cooled reactor safety: ingress of air or water into the core. Then, the paper offers an overview of what could be learned from the experiences with high temperature gas cooled reactors that have been built; their operating history indicates differences between actual operations and theoretical behavior. Finally, the paper describes some of the multiple priorities that often drive reactor design, and how safety is compromised in the process of optimizing other priorities.     

环状两相流空泡份额新析

对圆管内环状两相流的空泡份额进行了分析.首先推导出气相摩擦乘子的一般表达式,然后求得空泡份额与滑速比的解析式,并进行了分析讨论.与有关文献结果的比较表明:在一定的压力范围内文中所得的解是正确、可信的.     

Optimizing natural gas pipeline transmission with nonuniform elevation: A new initialization approach

In this paper, we develop an iterative piecewise linear approximation approach with a novel initialization method to solve natural gas pipeline transmission problems with the nonuniform network elevation. Previous approaches, such as energy minimization methods, cannot be applied to solve problems with the nonuniform network elevation because they exclude pressure range constraints, and thus provide solutions far from optimum. We propose a new initialization model that considers pressure range constraints and improves the optimality of the solutions and the computational efficiency. Furthermore, we extend the energy minimization methods and provide the necessary conditions under which the extended methods operate in networks with the nonuniform elevation. We test the performances of the methods with previously reported pipeline networks from the literature, with the open data set GasLib, and with our industrial collaborator. The initialization approach is shown to be more efficient than the method with fixed initial breakpoints. The newly proposed initialization approach generates solutions with a higher accuracy than the extended energy minimization methods, especially in large‐size networks. The proposed method has been applied to natural gas transmission planning by the China National Petroleum Corporation and has brought a direct profit increase of 330 million U.S. dollars in 2015‐2017.     

Effects of nano-sized aluminum on detonation characteristics and metal acceleration for RDX-based aluminized explosive

Nano-sized aluminum(Nano-Al)powders hold promise in enhancing the total energy of explosives and the metal acceleration ability at the same time.However,the near-detonation zone effects of reaction between Nano-Al with detonation products remain unclear.In this study,the overall reaction process of 170 nm Al with RDX explosive and its effect on detonation characteristics,detonation reaction zone,and the metal acceleration ability were comprehensively investigated through a variety of experiments such as the detonation velocity test,detonation pressure test,explosive/window interface velocity test and confined plate push test using high-resolution laser interferometry.Lithium fluoride(LiF),which has an inert behavior during the explosion,was used as a control to compare the contribution of the reaction of aluminum.A thermochemical approach that took into account the reactivity of aluminum and ensuing detonation products was adopted to calculate the additional energy release by afterburn.Combining the numerical simulations based on the calculated afterburn energy and experimental results,the param-eters in the detonation equation of state describing the Nano-Al reaction characteristics were calibrated.This study found that when the 170 nm Al content is from 0％to 15％,every 5％increase of aluminum resulted in about a 1.3％decrease in detonation velocity.Manganin pressure gauge measurement showed no significant enhancement in detonation pressure.The detonation reaction time and reaction zone length of RDX/Al/wax/80/15/5 explosive is 64 ns and 0.47 mm,which is respectively 14％and 8％higher than that of RDX/wax/95/5 explosive(57 ns and 0.39 mm).Explosive/window interface velocity curves show that 170 nm Al mainly reacted with the RDX detonation products after the detonation front.For the recording time of about 10 μs throughout the plate push test duration,the maximum plate velocity and plate acceleration time accelerated by RDX/Al/wax/80/15/5 explosive is 12％and 2.9 μs higher than that of RDX/LiF/wax/80/15/5,respectively,indicating that the aluminum reaction energy significantly increased the metal acceleration time and ability of the explosive.Numerical simulations with JWLM explosive equation of state show that when the detonation products expanded to 2 times the initial volume,over 80％of the aluminum had reacted,implying very high reactivity.These results are significant in attaining a clear understanding of the reaction mechanism of Nano-Al in the development of aluminized explosives.     

高焓高压空气加热器壁面传热过程数值计算

以一种高焓高压空气加热器为研究对象,对其冷却通道的流动和传热进行了三维数值仿真,冷却剂采用液态水,燃烧室和喷管分别采用不同的壁面材料,考虑壁面材料物性随温度的变化,并通过加热器热试验证了数值仿真结果的正确性。在此基础上,对比分析了气体辐射、冷却通道结构以及冷却水的流动方式对壁面换热的影响。结果表明:气体辐射对加热器燃烧室段壁面换热影响较大,对喷管壁面换热无明显的影响。在传热计算中,忽略气体辐射会引起较大的误差;冷却通道数和宽度存在最优组合,使得壁面换热量最大。冷却水的流动方式对燃气侧壁面温度影响较小,但对冷却液侧壁面温度的影响较大。     

高速艇自由变深水动力试验研究

通过模型试验,研究了高速艇从压载半潜状态自由变深到高速滑行状态的水动力现象.结果表明,当容积傅氏数大于1.25时,与舷外相通的水能完全自动排出.高速艇自由变深主要有以下5个原因:惯性作用、压载水逐渐排出导致的艇体上抬作用、艇艉由于负压形成的空穴对压载水的吸附作用、艇体动升力导致的艇体上抬作用、有利的艉倾作用.其中对具有较快加速运动的情形,惯性作用起重要地位,艇体动升力在排出流道水方面所发挥的作用相对较小,但在保持流道水排空方面发挥主导作用.     

轨道-线圈复合型电磁炮交叉作用研究

研究了轨道-线圈复合型电磁炮的驱动力、轨道-线圈交叉作用及其对系统性能的影响,提出了改进型的四轨复合炮.理论分析与数值计算表明:在普通的双轨复合炮中,交叉作用会使弹丸与炮管间存在较大的横向压力,导致弹丸运动不稳定及产生摩擦阻力,而四轨复合炮基本消除了交叉作用的不利影响,系统性能显著优于双轨复合炮.     

大气压力对火灾烟气影响的数值模拟

对环境大气压力不同的室内火灾分别进行了数值模拟,得到了两种环境条件下的火源热释放速率、温度、烟粒子浓度、CO浓度的变化数据,并对模拟结果进行了详细的分析与比较。研究结果表明：大气压力较低的地区,发生火灾时,燃烧不充分,热释放速率的两个峰值相距较近,稳定燃烧阶段短;烟气温度比常压地区低;烟粒子浓度变化不明显;CO浓度明显大于常压地区,因此烟气的毒性是低压地区威胁室内人员生命安全的重要因素。     

水压场引起海底地震波的特征

将舰船水压场激励的海底地震波简化为液固多层半无限空间中集中力引起的地震波动问题,基于波数积分方法,通过FFP数值积分求解了海底表面声压、位移和加速度幅值,得到了海底表面声压、位移和加速度的频率特性曲线。结果表明,水压场激励的海底地震波频率很低,属于超低频波动信号;水压场引起海底表面的声压、竖直加速度和水平加速度在频域内基本呈现单峰型特征。     

相对静止气氛下硼颗粒着火过程

针对相对静止气氛下单颗粒硼的着火过程展开了系统研究,考虑硼颗粒周围径向气相流动以及硼颗粒与周围环境间的传热和传质过程,建立了一维硼颗粒着火模型。分析了硼颗粒在实现着火和未能实现着火两种典型情形下的颗粒相、气相组分参数以及颗粒外表面上Stefan流的变化规律,并对其成因展开了分析。研究表明,在实现着火和未能实现着火两种典型情形下,硼颗粒外表面的Stefan流都会经历先由周围空间流向颗粒表面,而后变为由颗粒表面流向周围空间的过程;考虑颗粒外表面处气相Stefan流作用后,硼颗粒的着火延迟时间减少,且压力越大,影响越大。     

模型三组元喷嘴雾化SMD变化规律

为了充分研究新型三组元喷嘴的内混腔雾化特性 ,设计了模型喷嘴进行多工况试验。试验发现 ,维持气体或液体压降不变 ,改变另一种介质的压降 ,会得到不同的雾化平均直径随气体压降或液体压降 ,以及气液比ALR的变化曲线 ,这些曲线有着明显的规律。大量试验表明 ,虽然平均直径SMD随气体压降或液体压降的变化曲线明显不同 ,但其随气液比的变化规律却趋于相同     

船用锅炉排烟余热利用的研究

为了充分利用锅炉排烟余热,可将船用锅炉通风系统中配备一台燃气涡轮,在高负荷时进入烟囱前的锅炉排烟(即从经济器管束排出的烟气)流经该涡轮时所作的功可带动通风系统的压气机,但在低负荷时压气机需由燃气涡轮和辅助汽轮机联合驱动.经对该通风系统的压气机、燃气涡轮和汽轮机进行有关热力计算,结果表明:在锅炉通风系统中采用燃气涡轮来利用进入烟囱前的锅炉排烟余热,可以显著降低锅炉烟囱排烟的温度和节省蒸汽的消耗量.     

Experimental and numerical study of the fragmentation of expanding warhead casings by using different numerical codes and solution techniques

There has been increasing interest in numerical simulations of fragmentation of expanding warheads in 3D. Accordingly there is a pressure on developers of leading commercial codes, such as LS-DYNA, AUTODYN and IMPETUS Afea, to implement the reliable fracture models and the efficient solution techniques. The applicability of the Johnson–Cook strength and fracture model is evaluated by comparing the fracture behaviour of an expanding steel casing of a warhead with experiments. The numerical codes and different numerical solution techniques, such as Eulerian, Lagrangian, Smooth particle hydrodynamics (SPH), and the corpuscular models recently implemented in IMPETUS Afea are compared. For the same solution techniques and material models we find that the codes give similar results. The SPH technique and the corpuscular technique are superior to the Eulerian technique and the Lagrangian technique (with erosion) when it is applied to materials that have fluid like behaviour such as the explosive and the tracer. The Eulerian technique gives much larger calculation time and both the Lagrangian and Eulerian techniques seem to give less agreement with our measurements. To more correctly simulate the fracture behaviours of the expanding steel casing, we applied that ductility decreases with strain rate. The phenomena may be explained by the realization of adiabatic shear bands. An implemented node splitting algorithm in IMPETUS Afea seems very promising.