Se informan simulaciones de grandes vórtices de llamas coaxiales de O líquido/H gaseoso en condiciones de presión subcrítica en este artículo. Estas simulaciones reproducen los casos experimentales Mascotte A1, A10 y A30, que operan a 1, 10 y 30 bar, respectivamente, y para los cuales se disponen de mediciones de temperatura y visualizaciones experimentales. El objetivo principal de este trabajo es evaluar la precisión del modelo de equilibrio homogéneo de múltiples fluidos (HEM) descrito en Pelletier et al. (Computers & Fluids, 2020) para aplicaciones de motores de cohetes. De particular interés es la comparación con las mediciones de temperatura experimentales de Grisch et al. (Aerospace science and technology, 2003). Para ello, se realizan simulaciones numéricas con cuidado, con el fin de asegurar una adecuada convergencia estadística y estimar la influencia de la resolución de la malla para cada caso. A pesar de las suposiciones crudas-sin tensión superficial y sin modelo de atomización, por ejemplo-que se hacen con el HEM utilizado en este trabajo, se encuentra que los resultados están en razonable acuerdo con las mediciones para el caso A10, incluso con la malla más gruesa. Para el caso A30, se requiere una resolución de malla fina para capturar la zona de recirculación de baja intensidad aguas abajo del chorro interno necesaria para reproducir la forma del perfil experimental. Finalmente, las simulaciones del caso A1, con el número de Weber más bajo, muestran grandes desviaciones con respecto a las mediciones experimentales. Se espera que esto se deba a una deficiencia del modelo para reproducir adecuadamente el flujo disperso de dos fases.