Discontinuous Galerkin (dG) methods are well-suited for the implicit Large Eddy Simulation (iLES) of compressible turbulent flows, where the unfiltered Navier–Stokes equations are solved, and the dissipation of the numerical scheme plays the role of a subgrid-scale model that dissipates the smallest eddies. While LES is feasible for low- and moderate-Reynolds-number flows, its use in high-Reynolds-number configurations remains computationally prohibitive due to the stringent near-wall resolution requirements. To alleviate this cost, the inner part of the boundary layer can be modeled. In this work, we adopt a wall-modeled iLES strategy in which the no-slip boundary condition is replaced by a slip wall coupled with an algebraic wall-shear-stress model based on the Reichardt law. The model is implemented within an explicit entropy-stable dG solver that advances conservative variables in time and uses an L_2 projection onto entropy variables to assemble the spatial discretization. The formulation balances robustness and efficiency by adopting a physical-space orthonormal basis, which yields an identity mass matrix that simplifies explicit updates, while enforcing the second law of thermodynamics at the discrete level via entropy projection at a moderate additional cost~\cite{Alberti.ea:2024}. The approach is assessed in the turbulent channel flow case, demonstrating encouraging agreement with available reference data. The investigation will be extended to an adverse pressure gradient flow problem proposed in the ROSAS project (https://www.rosas-project.eu/). Additionally, preliminary results using more sophisticated wall models that solve coupled ordinary differential equations in the wall normal direction will be shown.