Numerical simulation of wall shear stress and boundary layer flow from jetting cavitation bubble on unheated and heated surfaces

• Published in: International journal of heat and mass transfer 2024-05, Vol.222, p.125189, Article 125189
• Main Authors: Park, Seong-Ho, Phan, Thanh-Hoang, Nguyen, Van-Tu, Duy, Trong-Nguyen, Nguyen, Quang-Thai, Park, Warn-Gyu
• Format: Article
• Language: English
• Subjects: Boundary layer flow; Heated surface; Jetting cavitation bubble; Wall shear stress
• ISSN: 0017-9310; 1879-2189
• DOI: 10.1016/j.ijheatmasstransfer.2024.125189

•Wall shear stress distribution on the unheated and heated surfaces during the bubble expansion and collapse were investigated using a fully compressible mixture model.•Complex boundary layer flow including adverse pressure gradient and heat transfer in the liquid film were simulated.•Effects of microjet velocity and liquid film thickness on peak shear stress were discussed.•Correlation equations for the peak shear stress on heated surface under different Prandtl number flow were suggested. 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.