ABSTRACT High-order accurate Large-eddy simulations (LES) have been carried out for compressible, turbulent flow through two bell-shaped convergent-divergent nozzles- one with a short divergent portion and the other with a longer one. An explicit filtering approach based on approximate deconvolution method is used for the LES. The nozzles have isothermal walls and circular cross-sections. The incoming flow is a turbulent, subsonic fully developed pipe flow at centreline Mach number of around 0.4 and friction Reynolds number 216. An adaptive filter is used to capture the shocks appearing in the divergent portion of the nozzles. The complicated behaviour of the mean flow, Reynolds stresses, production and pressure-strain terms of the axial Reynolds stress budget in the shock train region is shown. The broadband nature of the shock train oscillations having multiple dominant frequencies in these nozzles is observed in the present LES.