The physics of wall shear stress and boundary layer flow from jetting cavitation bubbles on unheated and heated surfaces was investigated numerically using a homogeneous mixture model based on a fully compressible Navier–Stokes solver. A spatial–temporal shear stress map and velocity field on unheated surfaces under several standoffs were considered first. Complex boundary layer flow, including separation vortex caused by an adverse pressure gradient and the generation of a thermal boundary layer because of skin friction with the shear flow, was discussed in detail. The tendency of different maximum shear stress for each standoff section was discussed with liquid film thickness and jet velocity. Also, we suggested maximum shear stress correlation equation for minimum film thickness and wall impact pressure as dimensionless form. As an expansion of previous studies, we examined bubble collapse–induced wall shear stress on a heated surface under various surface temperatures. The peak shear stress decreased with the lower Prandtl number flow. In the collapse phase, surface cooling by downward microjet and heat transfer with ring vortex were also captured. Consequently, we suggested correlation equations for primary and secondary peak shear stresses as function of Prandtl number and confirmed that the boundary layer flow including heat transfer can be well simulated.