Múltiples aplicaciones en la industria aeroespacial utilizan cargas pirotécnicas para su funcionamiento, y estas cargas son predominantemente en forma de gránulos. Una de las cargas más utilizadas es el nitrato de potasio y boro (BPN), y el presente estudio se centra en modelar matemáticamente la combustión granular, su recreación experimental y llevar a cabo un estudio comparativo sobre tres esquemas de upwind para su simulación numérica. En este artículo se presenta una versión personalizada de la formulación multifluido compresible de siete ecuaciones para modelar matemáticamente la combustión granular. Se utilizan tres esquemas de upwind, a saber, HLLC, AUSM-up y HLLC-AUSM, para la comparación numérica. Utilizando estos, se desarrolla un código axisimétrico para el estudio comparativo. Para replicar experimentalmente la combustión granular, se dispara BPN granular en una instalación de prueba de bomba cerrada, y la historia de presión experimental se utiliza para las comparaciones numéricas. El código desarrollado puede predecir adecuadamente la física de la combustión granular, y los tres esquemas son igualmente capaces de realizar predicciones numéricas.