一个高效的离散Boltzmann爆轰模型

构建一个适用于爆轰过程模拟的离散Boltzmann模型.该模型由一个离散Boltzmann方程和一个唯象反应率方程构成;在物理建模上,它等效于一个传统Navier-Stokes模型外加一个关于热动非平衡行为的粗粒化模型.与传统流体模型相比,它能够提供更多的动力学和动理学信息.该模型采用16个离散速度,相比于使用33个离散速度的模型具有更高的运算效率,模型中引入了额外自由度,通过调节额外自由度的数目,可以模拟各种不同比热比的爆轰.采用爆轰问题中的一些经典算例对所建立的模型进行数值验证.结果表明：该模型不仅对传统流体模型所能模拟的爆轰问题有效,而且能够用于一些传统流体模型不能描述的非平衡过程,有利于对爆轰问题的深入研究.     

A parabolic linear evolution equation for cellular detonation instability

Using the combined limits of a large activation energy and a ratio of specific heats close to unity, a dispersion relation has recently been derived which governs the stability of a steady Chapman - Jouguet detonation wave to two-dimensional linear disturbances. The analysis considers instability evolution time scales that are long on the time scale of fluid particle passage through the main reaction layer. In the following, a simplified polynomial form of the dispersion relation is derived under an appropriate choice of a distinguished limit between an instability evolution time scale that is long on the time scale of particle passage through the induction zone and a transverse disturbance wavelength that is long compared to the hydrodynamic thickness of the induction zone. A third order in time, sixth order in space, parabolic linear evolution equation is derived which governs the initial dynamics of cellular detonation formation. The linear dispersion relation is shown to have the properties of a most unstable wavenumber, leading to a theoretical prediction of the initial detonation cell spacing and a wavenumber above which all disturbances decay, eliminating the growth of small-wavelength perturbations. The role played by the curvature of the detonation front in the dynamics of the cellular instability is also highlighted.     

一维爆轰波在扰动来流中传播的数值研究

爆轰波在静止气体或定常来流中的传播得到了广泛研究, 然而在扰动来流中的传播研究较少。这方面的研究不仅是爆轰传播机制的重要组成部分, 还可为爆轰发动机的应用提供参考。文章基于两步诱导-放热总包反应模型, 开展了一维爆轰波在正弦密度扰动来流中的传播数值模拟。通过对数值结果分析, 获得了放热反应控制参数与爆轰波内在不稳定性的关系, 并在此基础上研究了扰动波长和幅值对一维爆轰波动力学过程的影响。研究发现, 在波前施加连续扰动会诱导爆轰波表现出更复杂的动力学行为, 且影响过程与爆轰波的内在不稳定性相关。对于稳定爆轰波, 扰动只在特定波长范围内引起前导激波后的压力振荡。对于不稳定爆轰波, 扰动会进一步强化其内在不稳定性。扰动幅值越大, 对爆轰波动力学过程的影响越显著。      

Stability of the square-wave detonation in a model system

Using a set of model equations for reactive flow, we study the stability of a “square-wave” detonation, in which each particle of the fluid reacts instantaneously after an induction time which depends on how hard it was shocked. We obtain a differential-difference equation for the shock velocity, valid for small perturbations about the steady solution. This equation is of so-called “advanced” type, in which the velocity at a given time depends on both velocity and acceleration at an earlier time.     

On turbulent chemical explosions into dilute aluminum particle clouds

We use a hybrid two-phase numerical methodology to investigate the flow-field subsequent to the detonation of a spherical charge of TNT with an ambient distribution of a dilute cloud of aluminum particles. Rayleigh–Taylor instability ensues on the contact surface that separates the inner detonation products and the outer shock-compressed air due to interphase interaction, which grows in time and results in a mixing layer where the detonation products afterburn with the air. At early times, the ambient particles are completely engulfed into the detonation products, where they pick up heat and ignite, pick up momentum and disperse. Subsequently, as they disperse radially outwards, they interact with the temporally growing Rayleigh–Taylor structures, and the vortex rings around the hydrodynamic structures results in the clustering of the particles by also introducing local transverse dispersion. Then the particles leave the mixing layer and quench, yet preserve their hydrodynamic ‘footprint’ even until much later; due to this clustering, preferential heating and combustion of particles is observed. With a higher initial mass loading in the ambient cloud, larger clusters are observed due to stronger/larger hydrodynamic structures in the mixing layer – a direct consequence of more particles available to perturb the contact surface initially. With a larger particle size in the initial cloud, clustering is not observed, but when the initial cloud is wider, fewer and degenerate clusters are observed. We identify five different phases in the dispersion of the particles: (1) engulfment phase; (2) hydrodynamic instability-interaction phase; (3) first vortex-free dispersion phase; (4) reshock phase; and (5) second vortex-free dispersion phase. Finally, a theoretical Buoyancy-Drag model is used to predict the growth pattern of the ‘bubbles’ and is in agreement with the simulation results. Overall, this study has provided some useful insights on the post-detonation explosive dispersal of dilute aluminum particle clouds.     

Generation of zonal magnetic fields by drift waves in a current carrying nonuniform magnetoplasma

It is shown that zonal magnetic fields (ZMFs) can be nonlinearly excited by incoherent drift waves (DWs) in a current carrying nonuniform magnetoplasma. The dynamics of incoherent DWs in the presence of ZMFs is governed by a wave-kinetic equation. The governing equation for ZMFs in the presence of nonlinear advection force of the DWs is obtained from the parallel component of the electron momentum equation and the Faraday law. Standard techniques are used to derive a nonlinear dispersion relation, which depicts instability via which ZMFs are excited in plasmas. ZMFs may inhibit the turbulent cross-field particle and energy transport in a nonuniform magnetoplasma.     

一维非定常炸轰的特征线方法数值模拟

本文用特征线方法计算了一维非定常炸轰的发展和传播问题。为模拟炸药的化学反应,采用了一种冲击波引爆的唯象模型。以PBX-9404高级炸药为例给出了计算的POP图与实验结果的比较以及不同时刻波后流场的分布。     

The power of distinguished limits: A comment on ‘On the two-dimensional instability of square-wave detonations’

The stability of both Chapman-Jouguet and overdriven square-wave detonations is investigated in the limit of a large detonation Mach number and a ratio of specific heats close to one.     

Numerical simulations of pulsating detonations: I. Nonlinear stability of steady detonations

Very-long-time numerical simulations of an idealized pulsating detonation with one irreversible reaction having an Arrhenius form are performed using a hierarchical adaptive second-order Godunov-type scheme. The initial data are given by the steady solution and the truncation error produces the perturbation to trigger the instability. The detonation is allowed to run for thousands of half-reaction times of the underlying steady wave to ensure that the final amplitudes and periods of the nonlinear oscillations are achieved. Thorough resolution studies are performed for various representative regimes of the instability. It is shown that to obtain quantitatively good solutions over 50 numerical grid points in the half-reaction length of the steady detonation are required, while to obtain a converged solution over 100 points are required, even near the stability boundary. This is much higher resolution than has generally been used in previous papers in either one or two dimensions. Resolutions of less than approximately 20 points per half-reaction length give very poor predictions of the periods and amplitudes near the stability boundary or entirely spurious solutions for more unstable detonations. The evolution of the converged solutions as the activation energy increases, and the detonation becomes more unstable, is also investigated.     

Statistical analysis of detonation wave structure

Hydrocarbon fueled detonations are imaged in a narrow channel with simultaneous schlieren and broadband chemiluminescence at 5 MHz. Mixtures of stoichiometric methane and oxygen are diluted with various levels of nitrogen and argon to alter the detonation stability. Ethane is added in controlled amounts to methane, oxygen, nitrogen mixtures to simulate the effects of high-order hydrocarbons present in natural gas. Sixteen unique mixtures are characterized by performing statistical analysis on data extracted from the images. The leading shock front of the schlieren images is detected and the normal velocity is calculated at all points along the front. Probability distribution functions of the lead shock speed are generated for all cases and the moments of distribution are computed. A strong correlation is found between mixture instability parameters and the variance and skewness of the probability distribution; mixtures with greater instability have larger skewness and variance. This suggests a quantitative alternative to soot foil analysis for experimentally characterizing the extent of detonation instability. The schlieren and chemiluminescence images are used to define an effective chemical length scale as the distance between the shock front and maximum intensity location along the chemiluminescence front. Joint probability distribution functions of shock speed and chemical length scale enable statistical characterization of coupling between the leading shock and following reaction zone. For more stable, argon dilute mixtures, it is found that the joint distributions follow the trend of the quasi-steady reaction zone. For unstable, nitrogen diluted mixtures, the distribution only follows the quasi-steady solution during high-speed portions of the front. The addition of ethane is shown to have a stabilizing effect on the detonation, consistent with computed instability parameters.     

On the characteristics of two-dimensional double cellular detonations with two successive reactions model

Double cellular detonations were numerically investigated using two-dimensional Euler equations with two successive chemical reactions, whose reaction lengths differ one order of magnitude. Simulated soot track images showed the double cellular structure with two cell widths that differ one order of magnitude, as well as previous experiments and numerical simulations. We successfully divided the double cellular detonation with two successive exothermic reactions into two detonations, primary and secondary detonations, with a single exothermic reaction, based on p–v relation of Rayleigh line and Hugoniot curves with the addition of the hypothetical condition of intermediate initial state. The ratio of cell widths of primary and secondary detonations showed good agreement with that caused by the first and second reactions of double cellular detonation, and there was no interaction between two successive chemical reactions. The linear stability analysis of planar detonation and the soot rack images of double, primary and secondary detonations showed that instabilities of primary and secondary detonations are dominant to that of double cellular detonation with two successive reactions. We confirmed the validity of division of two successive reactions to clarify the detonation instability and its cellular structure.     

Linear and nonlinear excitations in complex plasmas with nonadiabatic dust charge fluctuation and dust size distribution

Both linear and nonlinear excitation in dusty plasmas have been investigated including the nonadiabatic dust charge fluctuation and Gaussian size distribution dust particles. A linear dispersion relation and a Korteweg-de Vries-Burgers equation governing the dust acoustic shock waves are obtained. The relevance of the instability of wave and the wave evolution to the dust size distribution and nonadiabatic dust charge fluctuation is illustrated both analytically and numerically. The numerical results show that the Gaussian size distribution of dust particles and the nonadiabatic dust charge fluctuation have strong common influence on the propagation of both linear and nonlinear excitations.     

Clustering and combustion of dilute aluminum particle clouds in a post-detonation flow field

A hybrid two-phase numerical methodology is used to investigate the flow-field subsequent to the detonation of a spherical charge of TNT with an ambient distribution of a dilute cloud of aluminum particles. The interaction of the particle cloud with the contact surface results in Rayleigh–Taylor instability, which grows in time and gives rise to a mixing layer where the detonation products mix with the air and afterburn. At early times, the ambient particles get engulfed into the detonation products and ignite. Subsequently, they catch up with the Rayleigh–Taylor structures, and the vortex rings around the hydrodynamic structures cause transverse dispersion that results in the clustering of particles. Then, the particles leave the mixing layer and quench, yet preserve their hydrodynamic foot print. Preferential heating and combustion of particles occurs due to clustering. A higher initial mass loading in the ambient cloud results in larger clusters due to stronger/larger vortex rings around the hydrodynamic structures. A larger particle size results in the formation of fewer and degenerate clusters when the initial width of the cloud is larger. A theoretical model is used to predict the bubble amplitudes, and are in good accordance with the simulation results. Overall, this study has provided some useful insights on the explosive dispersal of dilute aluminum particle clouds and the gas dynamics of the flow field in the mixing layer.     

Dust acoustic waves in collisional uniform dense magnetoplasma

In order to study the characteristics of dust acoustic waves in a uniform dense dusty magnetoplasma system, a nonlinear dynamical equation is deduced using the quantum hydrodynamic model to account for dust–neutral collisions. The linear dispersion relation indicates that the scale lengths of the system are revised by the quantum parameter, and that the wave motion decays gradually leading the system to a stable state eventually. The variations of the dispersion frequency with the dust concentration, collision frequency, and magnetic field strength are discussed. For the coherent nonlinear dust acoustic waves, new analytic solutions are obtained, and it is found that big shock waves and wide explosive waves may be easily produced in the background of high dusty density, strong magnetic field, and weak collision. The relevance of the obtained results is referred to dense dusty astrophysical circumstances.     

Formation of transverse waves in oblique detonations

The structure of oblique detonation waves stabilized on a hypersonic wedge in mixtures characterized by a large activation energy is investigated via steady method of characteristics (MoC) calculations and unsteady computational flowfield simulations. The steady MoC solutions show that, after the transition from shock-induced combustion to an overdriven oblique detonation, the shock and reaction complex exhibit a spatial oscillation. The degree of overdrive required to suppress this oscillation was found to be nearly equal to the overdrive required to force a one-dimensional piston-driven detonation to be stable, demonstrating the equivalence of two-dimensional steady oblique detonations and one-dimensional unsteady detonations. Full unsteady computational simulations of the flowfield using an adaptive refinement scheme showed that these spatial oscillations are transient in nature, evolving in time into transverse waves on the leading shock front. The formation of left-running transverse waves (facing upstream) precedes the formation of right-running transverse waves (facing downstream). Both sets of waves are convected downstream away from the wedge in the supersonic flow behind the leading oblique front, such that the mechanism of instability must continuously generate new transverse waves from an initially uniform flow. Together, these waves define a cellular structure that is qualitatively similar to a normal propagating detonation.     

Hydrodynamic waves in regions with smooth loss of convexity of isentropes: general phenomenological theory

A general phenomenological theory of hydrodynamic waves in regions with smooth loss of convexity of isentropes is developed, based on the fact that for most media these regions in the p-V plane are anomalously small. Accordingly the waves are usually weak and can be described in the manner analogous to that for weak shock waves of compression. The corresponding generalized Burgers equation is derived and analyzed. The exact solution of the equation for steady shock waves of rarefaction is obtained and discussed.     

Pulsating instability of detonations with a two-step chain-branching reaction model: theory and numerics

The nonlinear dynamics of Chapman–Jouguet pulsating detonations are studied both numerically and asymptotically for a two-step reaction model having separate induction and main heat release layers. For a sufficiently long main heat release layer, relative to the length of the induction zone, stable one-dimensional detonations are shown to be possible. As the extent of the main reaction layer is decreased, the detonation becomes unstable, illustrating a range of dynamical states including limit-cycle oscillations, period-doubled and four-period solutions. Keeping all other parameters fixed, it is also shown that detonations may be stabilized by increasing the reaction order in the main heat release layer. A comparison of these numerical results with a recently derived nonlinear evolution equation, obtained in the asymptotic limit of a long main reaction zone, is also conducted. In particular, the numerical solutions support the finding from the analytical analysis that a bifurcation boundary between stable and unstable detonations may be found when the ratio of the length of the main heat release layer to that of the induction zone layer is O(1/ε), where ε (?1) is the inverse activation energy in the induction zone.     

线性偏振激光在相对论等离子体中的调制不稳定性

 从相对论等离体中电磁波的非线性色散方程出发,利用Karpman方法获得了线性偏振波模所满足的非线性控制方程,在非线性色散方程和非线性控制方程的基础上对线性偏振激光在相对论等离体中传播的调制不稳定性进行分析,给出了调制不稳定的时间增长率与扰动态波数之间的函数关系。     

The effects of mixture preburning on detonation wave propagation

Pressure gain combustion in the form of continuous detonations can provide a significant increase in the efficiency of a variety of propulsion and energy conversion devices. In this regard, rotating detonation engines (RDEs) that utilize an azimuthally-moving detonation wave in annular systems are increasingly seen as a viable approach to realizing pressure gain combustion. However, practical RDEs that employ non-premixed fuel and oxidizer injection need to minimize losses through a number of mechanisms, including turbulence-induced shock-front variations, incomplete fuel-air mixing, and premature deflagration. In this study, a canonical stratified detonation configuration is used to understand the impact of preburning on detonation efficiency. It was found that heat release ahead of the detonation wave leads to weaker shock fronts, delayed combustion of partially-oxidized fuel-air mixture, and non-compact heat release. Furthermore, large variations in wave speeds were observed, which is consistent with wave behavior in full-scale RDEs. Peak pressures in the compression region or near triple points were considerably lower than the theoretically-predicted values for ideal detonations. Analysis of the detonation structure indicates that this deflagration process is parasitic in nature, reducing the detonation efficiency but also leading to heat release far behind the wave that cannot directly strengthen the shock wave. This parasitic combustion leads to commensal combustion (heat release far downstream of the wave), indicating that it is the root cause of combustion efficiency losses.     

Detonability limits in thin annular channels

In this paper, detonability limits in two-dimensional annular channels are investigated. Since the channel heights are small in comparison to the tube diameter, curvature effects can be neglected and the annular channels can be considered to be essentially two-dimensional. Mixtures that are highly diluted with argon are used since previous investigations seem to indicate that detonations in such mixtures are “stable” in that cellular instabilities play minor roles on the propagation of the detonation. For stable detonations where the ZND structure is valid, boundary layer effects can be modeled as a flow divergence term in the conservation of mass equation following the pioneering work of Fay [J.A. Fay, Phys. Fluids 2(3) (1959) 283–289]. Expansion due to flow divergence in the reaction zone results in a velocity deficit. There exists a maximum deficit when an eigenvalue detonation velocity can no longer be found, which can be taken as the onset of the detonability limits. Experimentally, it was found that unlike “unstable” detonations, the detonability limits for “stable” detonations are well-defined. No unstable near-limit phenomena (e.g., galloping detonations) was observed. Good agreement is found between the theoretical predictions and the experimentally obtained velocity deficits and limits in the two channel heights of 2.2 and 6.9 mm for hydrogen–oxygen and acetylene–oxygen mixtures diluted with over 50% argon. It may be concluded that at least for these special mixtures where the detonation is “stable,” the failure mechanism is due to flow divergence caused by the negative displacement thickness of the boundary layer behind the leading shock front of the detonation wave.