A novel penalized wall function method for simulations of wall-bounded compressible turbulent flows is proposed. The new approach based on the Reynolds-averaged Navier–Stokes (RANS) equations to model outer region boundary layer, while inner part approximated by equilibrium model. differential formulation match external and solutions reformulated in a form generalized characteristic-based volum...