Simulación por Elementos Finitos de las Ecuaciones de Navier-Stokes en 3D

Abstract

We describe the computational and numerical aspects of a FEM code solving the Navier-Stokes equations in 3Ddomains. We use a discretization method belonging to the actual trends to solve incompressible flow as a limit of a compressible one with equal order interpolation functions for velocity and pressure variables and solve the equations appeared by an explicit pseudo transient integration that allow to simulate problems of medium complexity using minimal computational resources. The algorithm include automatic local time stepping to improve the convergence rate specially in refined meshes. Finally we plot results to validate the code for the lid driven square cavity at Re from 1 to 1000 and present a second result over the multiply connected square cavity at Re from 0 to 250.