High performance computing of stiff bubble collapse on CPU-GPU heterogeneous platform

Assessment of volume penalization and immersed boundary methods for compressible flows with various thermal boundary conditions

“…On the other hand, to properly model the wave reflection at the fluid-solid interface, Liu and Vasilyev [26] stated that the porosity parameter ϕ must be as small as possible. However, this leads to stability issues since the density penalization term is treated explicitly because of the use of a HLLC scheme for convective terms [37,38]. Liu & Vasilyev [26] conducted a stability analysis for subsonic flows assuming that u = O(η) within the penalized area.…”

“…The viscous flux is computed with a second-order centered scheme. More details about the present solver are given in [37,38]. The solution at time t n is first updated to an intermediate solution with the convective and viscous fluxes at first order in time.…”

Prediction of shock heating during ultrasound-induced bubble collapse using real-fluid equations of state

“…To account for these deficiencies, further CFD studies have progressively been incorporating compressibility effects in bubble collapse simulations. In the majority of them, the Stiffened Gas (SG) EoS [59] , [60] , [61] , [62] , [63] , [64] and Tait [65] , [66] , [67] , [68] , [69] , [70] , [71] EoSs have been employed. Although these two EoSs consider liquid compressibility in a straightforward manner, they are unable to describe the liquid thermodynamic behaviour appropriately as the first one uses an unphysical specific heat ratio while the latter does not account for density changes due to temperature variations as explained in Section 2 .…”