| Abstract: | Finite element solutions of the Euler and Navier-Stokes equations are presented, using a simple dissipation model. The discretization is based on the weak-Galerkin weighted residual method and equal interpolation functions for all the unknowns are permitted. The nonlinearity is iterated upon using a Newton method and at each iteration the linear algebraic system is solved by a direct solver with all unknowns fully coupled. Results are presented for two-dimensional transonic inviscid flows and two- and three-dimensional incompressible viscous flows. Convergence of the algorithm is shown to be quadratic, reaching machine accuracy in very few iterations. The inviscid results demonstrate the existence of nonunique numerical solutions to the steady Euler equations. |
