City Research Online

Solution of cavitating compressible flows using Discontinuous Galerkin discretisation

Papoutsakis, A. ORCID: 0000-0002-5449-5921, Koukouvinis, P. ORCID: 0000-0002-3945-3707 & Gavaises, M. ORCID: 0000-0003-0874-8534 (2020). Solution of cavitating compressible flows using Discontinuous Galerkin discretisation. Journal of Computational Physics, 410, article number 109377. doi: 10.1016/


A methodology for modelling cavitating flows using a high-order Adaptive Mesh Refinement (AMR) approach based on the Discontinuous Galerkin method (DG) is presented. The AMR implementation used features on-the-fly adaptive mesh refinement for unstructured hybrid meshes. The specific implementation has been developed for the resolution of complex multi-scale phenomena where high accuracy p-adaptive discretisations are combined with an h-adaptive data structure. This approach accommodates the fine spatial resolution for the interface discontinuities and the shock waves observed in compressible cavitating flows. The Tait equation of state is used for the modelling of the liquid phase while an isentropic path is assumed for the liquid/vapour mixture. Second order spatial and a third order non-oscillatory temporal discretisation are used for the integration of the mass and momentum conservation equations, in order to resolve the flow structures responsible for the formation of cavitation bubbles and the resulting compression waves. Assessment of the developed methodology is performed for the one-dimensional advancement of a compressible liquid-vapour interface and the symmetric collapse of a spherical vapour bubble.

Following, results obtained with the developed multi-scale modelling AMR approach has revealed a complex bubble collapse mechanism near a rigid wall, providing evidence of processes that have been unknown before due to reduced resolution and dissipative nature of past simulations. The impinging jet accompanying the collapse of a bubble near a wall, was found to induce vortical structures, which result to the formation of a secondary cavitation of a wall-attached bubble at the vicinity of the impingement jet shear layer. At the final stages of the initial bubble collapse, the impinging jet was found to penetrate the centre-line of the wall bubble inducing its partial collapse. This secondary collapse results to a rich spatial structure of shock waves, interacting with the secondary bubbles. Moreover, the calculated pressure level are found to be much higher than those reported from previous methodologies.

Publication Type: Article
Additional Information: © 2020. This manuscript version is made available under the CC-BY-NC-ND 4.0 license
Publisher Keywords: Cavitation, Adaptive Mesh Refinement, Barotropic model, Discontinuous Galerkin method
Subjects: T Technology > TJ Mechanical engineering and machinery
Departments: School of Science & Technology > Engineering
SWORD Depositor:
[thumbnail of JCP_200220_accepted.pdf]
Text - Accepted Version
Available under License Creative Commons Attribution Non-commercial No Derivatives.

Download (4MB) | Preview


Add to AnyAdd to TwitterAdd to FacebookAdd to LinkedinAdd to PinterestAdd to Email


Downloads per month over past year

View more statistics

Actions (login required)

Admin Login Admin Login