“…The Navier-Stokes equations for steady flowfields are solved via an implicit algorithm with second order spatial accuracy, higher order term relaxation with a relaxation factor of 0.25 is used to accelerate convergence. Flux vectors are computed according to the Advection Upstream Splitting Method (AUSM), which accurately captures shock discontinuities and has been used in a wide range of problems concerning hypersonic aerodynamics [23,24]. Combustion within the nozzle is represented by the Evans and Schexnayder supersonic hydrogen flame model which considers 12 species and 25 chemical reactions as presented in Table 3 [25].…”