In this paper, the Pseudo shock structure in a convergent-long divergent duct is investigated using large eddy simulation on the basis of Smagorinsky-Lilly, Wall-Adapting Local Eddy-Viscosity and Algebraic Wall-Modeled LES subgrid models. The first objective of the study is to apply different subgrid models to predict the structure of Lambda form shocks system, while the ultimate aim is to obtain further control of the shock behavior. To achieve these goals, the dynamic grid adaption and hybrid initialization techniques are applied under the 3D investigation to reduce numerical errors and computational costs. The results are compared to the existing experimental data and it is found that the WMLES subgrid model results in more accurate predictions when compared to the other subgrid models. Subsequently, the influences of the divergent section length with the constant ratio of the outlet to throat area and, the effects of discontinuity of the wall temperature on the flow physics are investigated. The results indicate that the structure of compressible flow in the duct is affected by varying these parameters. This is then further discussed to provide a deeper physical understanding of the mechanism of Pseudo shock motion.