Melt casting of energetic materials is investigated, and a numerical model is formulated for the analysis of the coupled fluid flow, heat transfer, and stress fields involved in this phase-change process. The numerical model is based on a conservative multi block control volume method. The SIMPLE algorithm is employed along with an enthalpy method approach to model the solidification process. Results from the model are verified against analytical solutions, experimental results, and published numerical results for simplified cases. In the melt casting of RDX-binder mixtures, the very high viscosity of the melt limits the influence of melt convection. The impacts of different cooling conditions on the velocity, temperature and stress distributions, as well as on the solidification time, are discussed. The present model can be used to improve the quality of cast explosives, by optimizing and controlling the processing conditions.