Multiphase growth during crystallization severely affects deliverable output of explosive materials. Appearance and incomplete transformation of metastable phases are a major source of polymorphic impurities. This article presents a methodical and molecular level understanding of the metastable phase transformation mechanism during crystallization of cyclic nitramine explosives, viz. RDX, HMX and CL-20. Instantaneous reverse precipitation yielded metastable γ-HMX and β-CL-20 which undergo solution mediated transformation to the respective thermodynamic forms, β-HMX and ε-CL-20, following 'Ostwald's rule of stages'. However, no metastable phase, anticipated as β-RDX, was evidenced during precipitation of RDX, which rather directly yielded the thermodynamically stable α-phase. The γ→β-HMX and β→ε-CL-20 transformations took 20 and 60 minutes respectively, whereas formation of α-RDX was instantaneous. Density functional calculations were employed to identify the possible transition state conformations and to obtain activation barriers for transformations at wB97XD/6-311++G(d,p)(IEFPCM)//B3LYP/6-311G(d,p) level of theory. The computed activation barriers and lattice energies responsible for transformation of RDX, HMX and CL-20 metastable phases to thermodynamic ones conspicuously supported the experimentally observed order of phase stability. This precise result facilitated an understanding of the occurrence of a relatively more sensitive and less dense β-CL-20 phase in TNT based melt-cast explosive compositions, a persistent and critical problem unanswered in the literature. The crystalline material recovered from such compositions revealed a mixture of β- and ε-CL-20. However, similar compositions of RDX and HMX never showed any metastable phase. The relatively long stability with the highest activation barrier is believed to restrict complete β→ε-CL-20 transformation during processing. Therefore a method is suggested to overcome this issue.