This paper develops methods for interface-capturing in multi-phase flows. The main novelties of these methods are: (a) multi-component modelling that embeds interface structures into the continuity equation; (b) a new family of triangle/tetrahedron finite elements, in particular the P1DG − P2 (linear discontinuous between elements velocity and quadratic continuous pressure); (c) an interface-capturing scheme based on compressive control volume advection methods and high-order finite element interpolation methods; (d) a time stepping method that allows use of relatively large time step sizes; and (e) application of anisotropic mesh adaptivity to focus the numerical resolution around the interfaces and other areas of important dynamics. This modelling approach is applied to a series of pure advection problems with interfaces as well as to the simulation of the standard computational fluid dynamics benchmark test cases of a collapsing water column under gravitational forces (in two and three dimensions) and sloshing water in a tank. Two more test cases are undertaken in order to demonstrate the many-material and compressibility modelling capabilities of the approach. Numerical simulations are performed on coarse unstructured meshes to demonstrate the potential of the methods described here to capture complex dynamics in multi-phase flows.
|Number of pages||27|
|Journal||International Journal for Numerical Methods in Fluids|
|Early online date||27 Aug 2015|
|Publication status||Published - 10 Feb 2016|
We would like to thank the EPSRC MEMPHIS Multiphase Programme grant, the EPSRC Computational Modelling for Advanced Nuclear Power Plants project and the EU FP7 projects THINS and GoFastR for helping to fund this work.
- Anisotropic mesh adaptivity
- Compositional modelling
- Compressive advection
- Discontinuous Galerkin