A D3Q19 Hybrid Recursive Regularized Pressure based Lattice Boltzmann Method (HRR-P LBM) is assessed for the simulation of complex transonic flows. Mass and momentum conservation equations are resolved through a classical LBM solver coupled with a finite volume resolution of entropy equation for a complete compressible solver preserving stability, accuracy and computational costs. An efficient treatment for wall and open boundaries is coupled with a grid refinement technique and extended to the HRR-P LBM in the scope of compressible aerodynamics. A Vreman subgrid turbulence model and an improved coupling of immersed boundary method with turbulence wall model on Cartesian grid accounts for unresolved scales by Large-Eddy Simulation (LES). The validity of the present method for transonic applications is investigated through various test cases with increasing complexity starting from an inviscid flow over a 10% bump and ending with a turbulent flow over a ONERA M6 three-dimensional wing.