Numerical workflow for scale-resolving computations of space launcher afterbody flows with and without jets

A comprehensive framework for high fidelity computations of two-species compressible turbulent flows

“…In order to contribute to the development of advanced numerical methods for bi-species supersonic mixing layers, two bi-species ZDES [22] computations are performed in this study on a well-known supersonic coaxial jets experiment run by Cutler et al [27,1] and reproduced numerically by several authors in the past with RANS [1,18] or RANS/LES formulations [17,23,19,21]. To meet the numerical challenges described above, the computations are run with an original hybrid numerical scheme [28] based on Ducros's sensor [29] associated to a low dissipation version of the AUSM scheme [30,31]. This new scheme is able to switch automatically from a low-dissipation formulation in presence of vortical structures to a more robust formulation in presence of shock waves (see Sec.…”

“…In order to test different numerical strategies, we have adapted three mono-species schemes to bi-species flow: the AUSM+ [30], a robust upwind scheme, an adapted low-dissipative AUSM+(P) [31] (called here AUSM+(P) L-D) and a hybrid shock-capturing scheme combining the two previous AUSM type schemes to perform high-fidelity, scale-resolving simulations of compressible flows [28].…”

“…Based on the two numerical schemes presented above, a hybrid scheme has been recently proposed in [28].…”

“…A secondorder Minmod MUSCL reconstruction [52] is employed with the AUSM+ scheme for the internal flows and a third-order Koren-type MUSCL reconstruction [53] is used with the hybrid scheme to treat the mixing and the stimulated boundary layers. See [28] for more details on these MUSCL procedures as the adaptation to bi-species flows only consists in the straightforward extension of these procedures to six equations instead of five.…”

“…Air Argon For the ZDES computations, the convective fluxes are computed thanks to the presented hybrid scheme [28] and a Koren MUSCL reconstruction. The timestep for ZDES is set to ∆t CF D = 5.10 −8 s thus enabling to ensure that CF L = |u+c|∆t CF D ∆n < 10, with ∆ n the radius of the inscribed sphere inside a cell, everywhere in the inner mixing layer.…”