From Mesh to Meshlessa Generalized Meshless Formulation Based on Riemann Solvers for Computational Fluid Dynamics

Resumen

From mesh to meshless: A generalized meshless formulation based on Riemann solvers for Computational Fluid Dynamics This thesis deals with the development of high accuracy meshless methods for the simulation of compressible and incompressible flows. Meshless methods were conceived to overcome the constraints that mesh topology impose on traditional mesh-based numerical methods. Despite the fact that meshless methods have achieved a relative success in some particular applications, the truth is that mesh-based methods are still the preferred choice to compute flows that demand high-accuracy. Instead of assuming that meshless and mesh-based methods are groups of methods that follow independent development paths, in this thesis it is proposed to increase the accuracy of meshless methods by taking guidance of some successful techniques adopted in the mesh-based community. The starting point for the development is inspired by the SPH-ALE scheme proposed by Vila. Especially, the flexibility of the ALE framework and the introduction of Riemann solvers are essential elements adopted. High accuracy is obtained by using the Moving Least Squares (MLS) technique. MLS serves multiple tasks in the implemented scheme: high order reconstruction of Riemann states, more accurate viscous flux evaluation and the replacement of the limited kernel approximation by MLS approximation with polynomial degree consistency by design. The stabilization of the scheme for compressible flows with discontinuities is based on a posteriori stabilization technique (MOOD) that introduces a great improvement compared with the traditional a priori flux limiters. The MLSPH-ALE scheme is the first proposed meshless formulation that uses high order consistent MLS approximation in a versatile ALE framework. In addition, the procedure to obtain the semi-discrete formulation keeps track of a boundary term, which eases the implementation of the boundary conditions. Another important contribution is related with the general concept of the MLSPHALE formulation. The MLSPH-ALE scheme is proved to be a global meshless formulation that under some particular settings provides the same semi-discrete equations that other meshless formulations published. The MLSPH-ALE scheme has been tested for the computation of turbulent flows. The low dissipation inherent to the Riemann solver is compatible with the implicit LES turbulent model. The proposed formulation is able to capture the energy cascade in the subsonic regime where traditional SPH formulations are reported to fail.