A mathematical model that predicts hydroxylamine nitrate (HAN) (NH2OH center dot HNO3) stability is applied to aqueous solutions containing HAN, nitric acid and plutonium that are used in plutonium purification processes. The model estimates the stability of these solutions with respect to the rapid, hazardous, autocatalytic reaction of HAN with nitric acid that generates heat and gas. It also accounts for reaction kinetics, temperature changes, gas generation rates, solution volumes and flow rates, and distribution of plutonium and nitric acid between aqueous and organic phases. The model is applied to three typical process vessels used in solvent extraction purification of plutonium a countercurrent aqueous/organic plutonium stripping column, an oxidation column used for HAN and hydrazine destruction, and a plutonium rework tank. Both normal and off-normal process scenarios are modeled. Two of the off-normal scenarios lead to the rapid autocatalytic reaction of HAN with nitric acid where heat and gas are generated and that could lead to damage of the process equipment and/or release of hazardous plutonium solution from the vessel. In these two cases, stationary aqueous solutions containing HAN, Pu(III), and nitric acid were allowed to slowly react until conditions for the autocatalytic reaction were reached. (c) 2015 Elsevier Ltd. All rights reserved.