The depressurization of a vessel containing saturated or subcooled liquid may occur in a variety of industrial processes and often poses a potentially hazardous situation. A ID plane numerical model was developed for estimating the thermodynamic and the dynamic state of the boiling liquid during a boiling liquid expanding vapor explosion (BLEVE) event. Based on the choice of the initial nucleation sites density, the model predicts, simultaneously, the bubble growth processes in the liquid at the superheat-limit state, the front velocity of the expanding liquid, and the shock wave pressure formed by the liquid expansion through the air. Conditions of shock formation were found to be normally associated with high initial temperatures that can bring the liquid to its superheat-limit state during the initial depressurization. Furthermore, the high initial temperature also induces a generation of higher vapor pressures that forces a rapid mixture expansion. Model predictions of the shock wave strengths, in terms of TNT equivalence, were compared against those obtained by simple energy models. As expected, the simple energy models over predicts the shock wave strength. However, the simple model which accounts for the expansion irreversibility, produces results which are closer to current model predictions. (c) 2007 Elsevier Ltd. All rights reserved.