Development of a 3D compressible flow solver for PFEM fluid simulation

Abstract

In this master thesis, a 2D/3D free surface flow solver based on the PFEM method is developed. While an existing 2D implicit incompressible solver developed in a previous work by Marco-Lucio Cerquaglia is available at the MTFC and MN2L lab of Uliège, an explicit solver is developed especially in the goal of simulating 3D flows efficiently. To do so, the incompressible Navier-Stokes equations are replaced by weakly-compressible Navier-Stokes equations using a Tait-Murnaghan equation of state. The Computational Geometry Algorithms Library (CGAL) is used to handle everything related to Delaunay triangulation and alpha-shapes algorithm in 2D and 3D. The simulations results are then compared to analytical solutions and to Cerquaglia’s solver. The solvers seems able to correctly represent flow not to far from incompressibility, and 3D simulations are working.