“…The code uses MOL, which is an efficient semi-discrete approach for the solution of time-dependent partial differential equations (PDEs), in conjunction with (i) a higher-order spatial discretization scheme which chooses biased-upwind or biased-downwind schemes in a zone of dependence manner; (ii) a parabolic algorithm for the computation of axial pressure gradient which does not require the solution of an elliptic equation for pressure; (iii) an elliptic grid generator using body-fitted coordinate system for application to complex geometries. The validity and the predictive ability of the code were tested by applying it to the simulation of laminar/turbulent, isothermal/non-isothermal incompressible flows and comparing its predictions with either [3,4]. In successive studies by the same group [5,6], the code was further developed by incorporation of; the solution of species equations using finite rate chemistry model together with a Total Variation Diminishing (TVD) flux limiter based discretization scheme for the computation of convective derivatives [5]; a radiation submodel to account for radiative heat transfer [6], for the simulation of transient reacting radiating flows.…”