Shock initiation of explosives: the idealized condensed-phase model

Many current models of condensed-phase explosives employ reactionrate law models where the form of the rate has a power-law dependenceon pressure (i.e. proportional to pⁿ where n is an adjustableparameter). Here, shock-induced ignition is investigated usinga simple model of this form. In particular, the solutions arecontrasted with those from Arrhenius rate law models as studiedpreviously. A large n asymptotic analysis is first performed,which shows that in this limit the evolution begins with aninduction stage, followed by a sequence of pressure runaways,resulting in a forward propagating, decelerating, shocklesssupersonic reaction wave (a weak detonation). The theory predictssecondary shock and super-detonation formation once the weakdetonation reaches the Chapman–Jouguet speed. However,it is found that secondary shock formation does not occur untilthe weak detonation has reached a point close to the initiatingshock, whereas for Arrhenius rate laws the shock forms closerto the piston. Numerical simulations are then conducted forO(1) values of n, and it is shown that the idealized condensed-phasemodel can qualitatively describe a wide range of experimentallyobserved behaviours, from growth mainly at the shock, to smoothgrowth of a pressure pulse behind the shock, to cases wherea secondary shock and possibly a super-detonation form. Thenumerics are used to reveal the different evolutionary mechanismsfor each of these cases. However, the evolution is found tobe sensitive to n, with the whole range of behaviours coveredby varying n from about 3 to 5. The simulations also confirmthe predictions of the theory that pressure-dependent rate lawsare unable to describe homogeneous explosive scenarios wherea super-detonation forms very close to the point of initialrunaway.     

炸药爆炸作用下飞片的运动

刚性飞片在炸药爆炸作用下的一维抛掷问题,仅当爆炸气体多方指数等于三时可以求得解析解;一般情况下须利用计算机求出数值解本文利用了爆炸气体中反射冲击波的“弱”击波特性,使飞片运动和飞片后方流场之间相互耦合的复杂问题解耦而归结为求解常微分方程问题;然后用小参数摄动法求出多方指数接近于三的各种炸药驱动飞片问题的近似解析解,所得飞片终速和数值解符合很好;从而给出了用爆速和多方指数等两个炸药示性参量表出的估算飞片运动的良好近似公式.     

Bildung schwacher stationärer Stosswellen in relaxierenden Gasgemischen

Summary The formation of a fully dispersed shock wave in a binary mixture of relaxing gases is investigated. It is assumed that the gas is bounded by a piston at the left, and at timet=0 is in static thermodynamic equilibrium. The piston velocity changes from zero att=0 to a constant non-zero value att>0. A uniformly valid approximation for the resulting wave motion is found by the method of matched asymptotic expansions. It turns out that the formation of a steady shock wave depends on the temperature dependence of the vibrational specific heats. If the coefficient of the non-linear term in Burgers equation is greater than zero, a steady shock wave is formed. However, for 0 this is not the case. One finds that the cases 0 are not realized in nature. Numerical results for the shock formation in relaxing air indicate, how the accuracy of acoustic theory decreases due to non-linear effects as time proceeds.     

Self-similar convergence of a shock wave in a heat conducting gas

The problem of the convergence of a spherical shock wave (SW) to the centre, taking into account the thermal conductivity of the gas in front of the SW, is considered within the limits of a proposed approximate model of a heat conducting gas with an infinitely high thermal conductivity and a small temperature gradient, such that the heat flux is finite in a small region in front of the converging SW. In this model, there is a phase transition in the surface of the SW from a perfect gas to another gas with different constant specific heat and the heat outflow. The gas is polytropic and perfect behind the SW. Constraints are derived which are imposed on the self-similarity indices as a function of the adiabatic exponents on the two sides of the SW. In front of the SW, the temperature and density increase without limit. In the general case, a set of self-similar solutions with two self-similarity indices exists but, in the case of strong SW close to the limiting compression, there are two solutions, each of which is completely determined by the motion of the spherical piston causing the self-similar convergence of the SW.     

Stability of viscous contact wave for the radiative and reactive gas with free boundary

The free boundary problem about the stability of viscous contact wave for the radiative and reactive gas is established by a basic energy method under the small perturbation. The present pressure includes a fourth order term about the absolute temperature from radiation effect as well as the ideal polytropic part, which brings the main difficulty to prove the asymptotic stability of the viscous contact wave.     

Nonself-similar flow with a shock wave reflected from the center of symmetry and new self-similar solutions with two reflected shocks

In some problems concerning cylindrically and spherically symmetric unsteady ideal (inviscid and nonheat-conducting) gas flows at the axis and center of symmetry (hereafter, at the center of symmetry), the gas density vanishes and the speed of sound becomes infinite starting at some time. This situation occurs in the problem of a shock wave reflecting from the center of symmetry. For an ideal gas with constant heat capacities and their ratio γ (adiabatic exponent), the solution of this problem near the reflection point is self-similar with a self-similarity exponent determined in the course of the solution construction. Assuming that γ on the reflected shock wave decreases, if this decrease exceeds a threshold value, the flow changes substantially. Assuming that the type of the solution remains unchanged for such γ, self-similarity is preserved if a piston starts expanding from the center of symmetry at the reflection time preceded by a finite-intensity reflected shock wave propagating at the speed of sound. To answer some questions arising in this formulation, specifically, to find the solution in the absence of the piston, the evolution of a close-to-self-similar solution calculated by the method of characteristics is traced. The required modification of the method of characteristics and the results obtained with it are described. The numerical results reveal a number of unexpected features. As a result, new self-similar solutions are constructed in which two (rather than one) shock waves reflect from the center of symmetry in the absence of the piston.     

DEPENDENCE OF QUALITATIVE BEHAVIOR OF THE NUMERICAL SOLUTIONS ON THE IGNITION TEMPERATURE FOR A COMBUSTION MODEL

We study the dependence of qualitative behavior of the numerical solutions （obtained by a projective and upwind finite difference scheme） on the ignition temperature for a combustion model problem with general initial condition. Convergence to weak solution is proved under the Courant-Friedrichs-Lewy condition. Some condition on the ignition temperature is given to guarantee the solution containing a strong detonation wave or a weak detonation wave. Finally, we give some numerical examples which show that a strong detonation wave can be transformed to a weak detonation wave under some well-chosen ignition temperature.     

Traveling waves for a boundary reaction–diffusion equation

We prove the existence of a traveling wave solution for a boundary reaction–diffusion equation when the reaction term is the combustion nonlinearity with ignition temperature. A key role in the proof is plaid by an explicit formula for traveling wave solutions of a free boundary problem obtained as singular limit for the reaction–diffusion equation (the so-called high energy activation energy limit). This explicit formula, which is interesting in itself, also allows us to get an estimate on the decay at infinity of the traveling wave (which turns out to be faster than the usual exponential decay).     

A REMARK ON HOFER-ZEHNDER SYMPLECTIC CAPACITY IN SYMPLECTIC MANIFOLDS M×R~(2n)

ＡＲＥＭＡＲＫＯＮＨＯＦＥＲＺＥＨＮＤＥＲＳＹＭＰＬＥＣＴＩＣＣＡＰＡＣＩＴＹＩＮＳＹＭＰＬＥＣＴＩＣＭＡＮＩＦＯＬＤＳＭ×Ｒ２ｎＭＡＲＥＮＹＩＡｂｓｔｒａｃｔＴｈｅａｕｔｈｏｒｓｔｕｄｉｅｓｔｈｅＨｏｆｅｒＺｅｈｎｄｅｒｃａｐａｃｉｔｙａ...     

Global supersonic conic shock wave for the steady supersonic flow past a cone: Polytropic gas

In this paper, we establish the global existence and stability of a steady conic shock wave for the symmetrically perturbed supersonic flow past an infinitely long conic body as long as the vertex angle is less than a critical value. The flow is assumed to be polytropic, isentropic and described by a steady potential equation. Based on the delicate asymptotic expansion of the background solution, one can verify that the boundary conditions on the shock and the conic surface satisfy the “dissipative” property. From this property, by use of the reflected characteristics method and the special form of the shock equation, we show that the conic shock attached at the vertex of the cone exists globally in the whole space when the speed of the supersonic coming flow is appropriately large. On the other hand, we remove the smallness restriction on the sharp vertex angle in order to establish the global existence of a shock or a global weak solution, moreover, our proof approach is different from that in [Shuxing Chen, Zhouping Xin, Huicheng Yin, Global shock wave for the supersonic flow past a perturbed cone, Comm. Math. Phys. 228 (2002) 47-84] and [Zhouping Xin, Huicheng Yin, Global multidimensional shock wave for the steady supersonic flow past a three-dimensional curved cone, Anal. Appl. 4 (2) (2006) 101-132].     

The random projection method for a model problem of combustion with stiff chemical reactions

In this paper we extend the random projection method, recently proposed by the author and S. Jin [J. Comput. Phys. 163 (2000) 216] for under resolved numerical simulations of a qualitative model problem for combustion with stiff chemical reactions:

In this problem, the reaction time is small, making the problem numerically stiff. A classic spurious numerical phenomenon – the incorrect shock speed – occurs when the reaction time scale is not properly resolved numerically. The random projection method is introduced recently to handle this kind of numerical difficulty. The key idea in this method is to randomize the ignition temperature in a suitable domain. Several numerical experiments demonstrate the reliability and robustness of this method.     

Rapid cylindrically and spherically symmetric strong compression of a perfect gas

The problem of the rapid cylindrically and spherically symmetric strong compression of a perfect (non-viscous and non-heat-conducting) gas is solved. The term “rapid” denotes that the compression time is much less than the run time of a sound wave across the initial cylindrical or spherical volume, while the term “strong” in this case means the simultaneous attainment of as large a density and temperature as desired. By definition, rapid compression must begin in a strong shock wave, which propagates to the axis or centre of symmetry. When the shock wave approaches the centre of symmetry this flow is described by the self-similar Guderley equation with an unbounded rise in temperature, pressure and velocity and a finite increase in the density at the centre of symmetry both behind the arriving and behind the reflected shock waves. To obtain as high an increase in the density as desired one must add on a centred compression wave with focus at the centre of symmetry to the overtaking shock wave at the instant it arrives at the centre of symmetry C⁻-characteristic. Outside a small neighbourhood of the focus one can calculate, by the method of characteristics, the centred wave and the trajectory of the piston which produces it. As for any centred wave, this calculation must be carried out from the centre of symmetry. Since some of the parameters at the focus (certainly the pressure, temperature and velocity of the gas) are unbounded, it is necessary to preface the calculation by the method of characteristics by constructing an analytic solution which holds in a small neighbourhood of the centre of symmetry. Below, after constructing the required solution, the centred waves corresponding to it and the trajectories of the piston producing them are calculated.     

Nonlinear Spatial Evolution of Small Disturbances from Boundary Conditions in Flat Inclined-Channel Flow

The asymptotic behavior of small disturbances as they evolve spatially from boundary conditions in a flat inclined channel is determined. These small disturbances develop into traveling waves called roll waves, first discussed by Dressler in 1949. Roll waves exist if the Froude number F exceeds 2, which consist of a periodic pattern of bores, or discontinuities. After confirming the instability condition for   F > 2  for the linearized equations in the boundary value case, the nonlinear boundary value problem for the weakly unstable region of F slightly larger than 2 is studied. Multiple scales and the Fredholm alternative theorem are applied to determine the evolution of the solution in space. It is found that the solution is dominated by the evolution of the disturbance along one characteristic. The shock conditions governing the asymptotic solution are determined and these conditions are used to determine the approximate shape of the resulting traveling wave from the solution. Both asymptotic and numerical results for periodic disturbances are presented.     

On the acceleration of shock waves and concentration of energy

By a series of simple examples related to exact solutions of problems in gas dynamics and magnetohydrodynamics, possible mechanisms of acceleration of shock waves and concentration of energy are elucidated. The acceleration of a shock wave is investigated in the problem of motion of a plane piston at a constant velocity in the case when the initial density of the medium drops in the presence of constant counterpressure. It is shown that in this situation a “blow-up” regime is induced by a shock wave going to infinity in finite time even for limited work of the piston. A simple spherically symmetric solution with a converging shock wave is constructed and shown to lead to the concentration of energy. A general method for solving one-dimensional non-self-similar problems related to matching the equilibrium state to a motion with homogeneous deformation on a shock wave is discussed; this method leads to a solution in quadratures.     

一维可压粘性气体粘性激波的渐近稳定性

本文讨论一维粘性热传导多方气体粘性激波的渐近稳定性,如果初始扰动以及δ=|u+-u-|适当的小,则解在最大模的意义下趋于粘性激波.     

Stability of viscous contact wave for compressible navier-stokes system of general gas with free boundary

In this paper, we study the large time behavior of solutions to the nonisentropic Navier-Stokes equations of general gas, where polytropic gas is included as a special case, with a free boundary. First we construct a viscous contact wave which approximates to the contact discontinuity, which is a basic wave pattern of compressible Euler equation, in finite time as the heat conductivity tends to zero. Then we prove the viscous contact wave is asymptotic stable if the initial perturbations and the strength of the contact wave are small. This generalizes our previous result [6] which is only for polytropic gas.     

Numerical solutions of a model for the propagation of a surface-catalysed flame in a tube

Numerical simulations of a surface-catalysed flame in a tubeare performed, corresponding to an experiment where a premixedfuel is fed into a tube whose inner surface is coated with acatalyst. In these experiments, subsequent to ignition, a reactionwave can be seen as a red-hot region which propagates back alongthe tube towards the inlet, and is due to low temperature combustionoccurring only on the inner surface of the tube where the catalystis present. The solutions of a mathematical model for this behaviourshow that initial-value problems do indeed result in such steadilypropagating waves. The numerically obtained wave speeds andsteady solution are compared to a previous large Damköhlernumber (D_a) asymptotic analysis using a simple reaction ratemodel, and agreement is very good even for moderately largevalues of D_a. However, for such Damköhler numbers, thewave speeds are found to be much larger than observed experimentally.Indeed, the simulations show that O(1) values of D_a are requiredto obtain the lower experimental wave speeds. Nevertheless,the wave speeds as a function of flow rate through the tubedo not agree well with the preliminary experimental resultsfor any choice of the parameters. A more realistic, Arrheniusreaction rate model is then considered. The Arrhenius modelpredicts a rapid change in temperature at the wave front, inmuch better agreement with the experiments than for the simplerreaction model.     

有限变形弹性杆中三种非线性弥散波

在一维弹性细杆拉压、扭转和弯曲波的经典线性理论基础上,分别计入有限变形和弥散效应,借助Hamilton变分原理,由统一的方法导出了3种非线性弥散波的演化方程．对3种演化方程进行了定性分析．结果表明,这些方程在相平面上存在同宿轨道或异宿轨道,分别相应于孤波解或冲击波解．根据齐次平衡原理,用Jacobi椭圆函数展开对这些演化方程进行了求解,在一定的条件下它们均可能存在孤立波解或冲击波解,这与方程的定性分析完全一致．     

Focusing of shock waves in a highly viscous fluid

The method of combined asymptotic expansions is used to solve the problem of the focusing of a shock wave (in a weakly compressible medium of high viscosity. Asymptotic forms of the solution are constructed in a number of spatial zones. The focusing zone is described by its asymptotic form obtained by combining it with the solution corresponding to viscous geometrical acoustics. The reflection of a shock wave formed as a result of velocity jump near one of the foci of the ellipsoid of revolution is discussed as an example. Analytical relationships descrbing the focusing zone around the second focus are obtained. It is shown that at the focus itself the wave profile has an antisymmetric form, and the compression wave is followed by a rarefaction wave of the same form.     

Vanishing pressure limits of Riemann solutions to the isentropic relativistic Euler system for Chaplygin gas

The Riemann solutions to the isentropic relativistic Euler system for Chaplygin gas with a small parameter are considered. Unlike the polytropic or barotropic gas cases, we find that firstly, as the parameter decreases to a certain critical number, the two-shock solution converges to a delta shock wave solution of the same system. Moreover, as the parameter goes to zero, that is, the pressure vanishes, the solution is nothing but the delta shock wave solution to the zero-pressure relativistic Euler system. Meanwhile, the two-rarefaction wave solution tends to the vacuum solution to the zero-pressure relativistic system, and the solution containing one rarefaction wave and one shock wave tends to the contact discontinuity solution to the zero-pressure relativistic system as pressure vanishes.