Skip Navigation


IMA Journal of Applied Mathematics Advance Access originally published online on January 8, 2008
IMA Journal of Applied Mathematics 2008 73(2):361-373; doi:10.1093/imamat/hxm065
This Article
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
73/2/361    most recent
hxm065v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Sharpe, G. J.
Right arrow Articles by Short, M.
Right arrow Search for Related Content
Social Bookmarking
Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© The Author 2008. Published by Oxford University Press on behalf of the Institute of Mathematics and its Applications. All rights reserved.

Shock initiation of explosives: the idealized condensed-phase model

G. J. Sharpe{dagger}

School of Mechanical Engineering, University of Leeds, Leeds LS2 9JT, UK

V. Gorchkov

Theoretical and Applied Mechanics, University of Illinois Urbana, IL 61801, USA

M. Short

Los Alamos National Laboratory, New Mexico, NM 87545, USA

{dagger} Email: g.j.sharpe{at}leeds.ac.uk

Received on July 6, 2006; Accepted on May 18, 2007

Many current models of condensed-phase explosives employ reaction rate law models where the form of the rate has a power-law dependence on pressure (i.e. proportional to pn where n is an adjustable parameter). Here, shock-induced ignition is investigated using a simple model of this form. In particular, the solutions are contrasted with those from Arrhenius rate law models as studied previously. A large n asymptotic analysis is first performed, which shows that in this limit the evolution begins with an induction stage, followed by a sequence of pressure runaways, resulting in a forward propagating, decelerating, shockless supersonic reaction wave (a weak detonation). The theory predicts secondary shock and super-detonation formation once the weak detonation reaches the Chapman–Jouguet speed. However, it is found that secondary shock formation does not occur until the weak detonation has reached a point close to the initiating shock, whereas for Arrhenius rate laws the shock forms closer to the piston. Numerical simulations are then conducted for O(1) values of n, and it is shown that the idealized condensed-phase model can qualitatively describe a wide range of experimentally observed behaviours, from growth mainly at the shock, to smooth growth of a pressure pulse behind the shock, to cases where a secondary shock and possibly a super-detonation form. The numerics are used to reveal the different evolutionary mechanisms for each of these cases. However, the evolution is found to be sensitive to n, with the whole range of behaviours covered by varying n from about 3 to 5. The simulations also confirm the predictions of the theory that pressure-dependent rate laws are unable to describe homogeneous explosive scenarios where a super-detonation forms very close to the point of initial runaway.

Keywords: detonation; explosives; numerical simulation; shock waves.


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?




Disclaimer: Please note that abstracts for content published before 1996 were created through digital scanning and may therefore not exactly replicate the text of the original print issues. All efforts have been made to ensure accuracy, but the Publisher will not be held responsible for any remaining inaccuracies. If you require any further clarification, please contact our Customer Services Department.