“…The inviscid flux at each cell face is evaluated using limited linear reconstruction [43] and Riemann-solver based flux functions [44,45], while the viscous flux is evaluated utilizing a hybrid average gradient-diamond path method [46]. A standard, explicit, two-stage, second-order-accurate, RungeKutta, time-marching scheme is used to integrate forward in time the non-linear, coupled-system of ordinary differential equations resulting from the finite-volume spatial discretization procedure.…”