Prediction of multiphase flows with sharp interfaces using anisotropic mesh optimisation

摘要

We propose an integrated, parallelised modelling approach to solve complex multiphase flow problems with sharp interfaces. This approach is based on a finite-element, double control-volume methodology, and employs highly-anisotropic mesh optimisation within a framework of high-order numerical methods and algorithms, which include adaptive time-stepping, metric advection, flux limiting, compressive advection of interfaces, multi-grid solvers and preconditioners. Each method is integral to increasing the fidelity of representing the underlying physics while maximising computational efficiency, and, only in combination, do these methods result in the accurate, reliable, and efficient simulation of complex multiphase flows and associated regime transitions. These methods are applied simultaneously for the first time in this paper, although some of the individual methods have been presented previously. We validate our numerical predictions against standard benchmark results from the literature and demonstrate capabilities of our modelling framework through the simulation of laminar and turbulent two-phase pipe flows. These complex interfacial flows involve the creation of bubbles and slugs, which involve multi-scale physics and arise due to a delicate interplay amongst inertia, viscous, gravitational, and capillary forces. We also comment on the potential use of our integrated approach to simulate large, industrial-scale multiphase pipe flow problems that feature complex topological transitions.