Abstract
This work describes a Finite Element Newton Method for the solution of the stationary
Navier–Stokes equations for two-dimensional incompressible flows. We start from the weak
variational formulation of the problem and adopt an unequal order interpolation
\$P^1\$–\$P^2\$ for pressure and velocity. Rather general boundary conditions are considered. The Newton
method for the nonlinear system of coupled equations is written in a particularly
transparent incremental form and the Jacobian linear system is solved by means of a dire
ct algorithm (MUMPS). The results of some numerical tests are provided to demonstrate the correctness and capability of the method.
Users
Please
log in to take part in the discussion (add own reviews or comments).