DG-XFEM formulation for the unsteady incompressible Navier-Stokes equations

Abstract

This Thesis proposes a combined formulation of the Discontinuous Galerkin Method (DG) with solenoidal basis functions and the eXtended Finite Element Method (XFEM), in order to solve the incompressible Navier-Stokes equations for unsteady flows around a solid object, providing high orders of accuracy in space and time. This DG-XFEM formulation simplifies the meshing process using structured meshes that also do not need to be updated at every time step if the object moves, reducing the computational cost.

In the DG-XFEM formulation a fixed structured mesh is used and its elements are classified in three groups, which receive a different treatment. First, the elements inside the solid object are excluded in the calculations since it is treated as a void. Second, the elements belonging to the fluid are calculated as in the DG solenoidal formulation. Third, for the elements cut by the interface integration is modified using XFEM in order to take into account only the fluid region, considering curved integration cells to accurately compute integrals in high-order elements; straight afterwards it is solved again with the DG solenoidal formulation. In the DG solenoidal formulation incompressible flows are first solved for velocity and only part of the pressure's degrees of freedom (hybrid pressure), reducing the overall size of the system to be solved, while the rest of pressure degrees of freedom (interior pressure) is computed as a postprocessing. A numerical validation of the method is given with the simulation of the classical benchmark test of the flow past a cylinder, showing its good performance in several cases tested.

En aquesta Tesina es proposa una formulació combinada de dos mètodes numèrics: ”Discontinuous Galerkin Method” (DG) amb funcions de base solenoïdals i ”eXtended Finite Element Method”(XFEM), per tal de resoldre les equacions de Navier-Stokes per a fluxos no estacionaris al voltant d’un objecte sòlid, proporcionant alts ordres de precisió en l’espai i el temps. Aquesta formulació combinada DG-XFEM simplifica el procés de mallat, utilitzant malles estructurades que a més no necessiten ser actualitzades en cada pas de temps si l’objecte es mou, reduint per tant el cost computacional. En la formulació DG-XFEM s’empra una malla estructurada fixa i els seus elements es classifiquen en tres grups, els quals reben tractaments diferents. Primerament, els elements dins del sòlid són exclosos en els càlculs ja que es tracten com a buits. A continuació, els elements que pertanyen al fluid es calculen igual que en la formulació de DG solenoïdal. Per últim, en els elements tallats per la interfase es modifica la integració fent servir XFEM per tal de tenir en compte només la regió de fluid, considerant cèl·lules d’integració corbes per calcular de forma acurada les integrals en elements d’alt ordre; tot seguit es resol de nou emprant la formulació de DG solenoïdal. En la formulació de DG solenoïdal els fluxos incompressibles es resolen en primer lloc per a les velocitats i només una part dels graus de llibertat de la pressió (pressió híbrida), reduint la grandària global del sistema d’equacions a resoldre, mentre que la resta de graus de llibertat (pressió interior) es calculen com un post-procés. El mètode està validat amb la simulació del clàssic problema de test d’un flux de fluid al voltant d’un cilindre, mostrant el seu bon comportament en vàries simulacions.