The main objectives of this work are the formulation, implementation and applications of a numerical algorithm to simulate turbulent, incompressible, isothermal flows. The transient three-dimensional flow is analyzed using an explicit Taylor-Galerkin scheme and the finite element method with hexahedrical eight-node element. Turbulence is simulated using Large Eddy Simulation. For sub-grid scales two models were implemented, the classical Smagorinsky's model and the dynamic eddy viscosity model. For the second filtration, which is necessary in the dynamic model, a new method was developed based on independent finite elements that involve each node in the original mesh. The implemented scheme is efficient and good results with low additional computational cost were obtained. Results for two classical problems, the driven cavity and the backward facing step are presented. Comments about the model applicability for flows with high Reynolds numbers are also presented.
Turbulence; large Eddy simulation; dynamic model; finite element method; computational fluid dynamics
