A finite element algorithm for turbulent flow processing a k-epsilon model

The first results in the development of a three-dimensional Navier-Stokes (NS) finite element code using the splitting up technique for simulating industrial flows are presented. The use of the Ronat algorithm in the splitting up of operators is discussed and applied to the NS equations. The treatment of the two main difficulties due to turbulence modeling for a three-dimensional finite-element code, the reformation of the matrix due to the change of the viscosity field at each time step, and the use of a wall function to avoid the need for an excessively thin grid along the wall due to the boundary layer is discussed. The use of the k-epsilon eddy viscosity model for turbulence modeling is addressed.