This paper presents zonal detached-eddy simulation of a round underexpanded sonic jet exhausting from a realistic aircraft afterbody and controlled by four radial injections. The comparison with a former experiment proves the capability of zonal detached-eddy simulation to reproduce mean flow induced by the complex interaction between the compressible waves and the turbulent jet shear layer. However, the very near-field jet structure indicates that the shear-layer transition to a fully turbulent state is delayed. Therefore, a modification of the characteristic length scale used in the modeling is proposed which is shown to improve the prediction of the shear-layer growth. The control injectors are aimed to enhance jet mixing. Experiment indicated the efficiency of this additional device but also expressed the need for more data to get a deeper and more complete insight into the whole flowfield. Both spatial and temporal information are provided by zonal detached-eddy simulation which is then a good candidate. This method is shown to correctly simulate streamwise vortices generated by secondary injections, their action in jet distortion, as well as their subsequent decay. Finally, computational results allow us to evaluate the action of the hypermixer on turbulent activity and its consequences on jet dilution. Nomenclature a = critical speed of the primary flow D = nozzle throat diameter d = distance to the closest wall k = turbulent kinetic energy k = nondimensionalized turbulent kinetic energy integrated in jet cross sections Y; Z M = Mach number N = normalized vorticity vector p = pressure hp i jp i0 i = nondimensionalized stagnation pressure averaged in the jet volume Pr = Prandtl number Q = mass flow rate S ij = components of the strain tensor T = temperature U = longitudinal velocity u = velocity vector V c = core volume V j = jet volume X, Y, Z = Cartesian coordinates linked to the scale model X, R, = cylindrical coordinates x, y, z = local coordinates linked to the grid cells = mesh characteristic length max = based on the maximum cell size vol = based on the cell volume ! = based on the vorticity x, y, z = cell sizes in the local coordinate system = small parameter = small vector = molecular dynamic viscosity t = turbulent dynamic viscosity t = turbulent kinematic viscosity = density d ij = components of the deviatoric part of the Reynolds-averaged Navier-Stokes/subgridscale stress tensor R xr , RM xr , RR xr = total, modeled, and resolved Reynolds shear stress in the plane (X;R) X = normalized streamwise vorticity ! xD=2a ! = vorticity vector r u Subscripts i = isentropic stagnation value 0 = ambient flow value 1 = primary flow value 2 = ventilation flow value 3 = control flow value IntroductionA ERONAUTICAL applications and more generally many engineering issues are concerned with high Mach and Reynolds-number flows featuring attached and detached boundary layers in internal as well as external parts. Experiment has proved to be essential in the past while computational tools were still not mature, and ...