An All-Mach Number HLLC-Based Scheme for Multi-Phase Flow with Surface Tension

Abstract

This paper presents an all-Mach method for two-phase inviscid flow in the presence of surface tension. A modified version of the Hartens–Lax–van Leer Contact (HLLC) solver is developed and combined for the first time with a widely used volume-of-fluid (VoF) method: the compressive interface capturing scheme for arbitrary meshes (CICSAM). This novel combination yields a scheme with both HLLC shock capturing as well as accurate liquid–gas interface tracking characteristics. It is achieved by reconstructing non-conservative (primitive) variables in a consistent manner to yield both robustness and accuracy. Liquid–gas interface curvature is computed via height functions and the convolution method. We emphasize the use of VoF in the interest of interface accuracy when modelling surface tension effects. The method is validated using a range of test-cases available in the literature. The results show flow features that are in sensible agreement with previous experimental and numerical work. In particular, the use of the HLLC-VoF combination leads to a sharp volume fraction and energy field with improved accuracy.

Application of Central-Weighted Essentially Non-Oscillatory Finite-Volume Interface-Capturing Schemes for Modeling Cavitation Induced by an Underwater Explosion

“…We performed a quantitative investigation on the medium grid for the obtained density and pressure profiles along the horizontal axis (x = 0) and compared them with published results that utilize sharpening techniques and anti-diffusive techniques in conjunction with high-order schemes [9,49]. The results, as shown in Figures 7 and 8, are consistent with previously published findings [9,49,50].…”

“…Underwater explosion near a free surface; pressure plot along horizontal axis x = 0 at instants t = (0.4; 0.8) ms. Comparison are made with the results of[49,50,52]…”

Sloshing Wing Dynamics - 2nd Year Project Overview

“…Due the state-of-the-art nature of this work, it resulted in several conference and journal publications, Ref. [16][17][18][19].…”

A higher-order accurate VOF interface curvature computation scheme for 3D non-orthogonal structured meshes

“…The proposed scheme was implemented in the Elemental ® multiphysics code [39][40][41][42][43]. The first step in demonstrating higher-order accuracy of VOF interface reconstruction is the VOF field initialisation.…”