A finite element analysis of the exact nonlinear formulation of a lifting surface in steady incompressible flow, with the evaluation of the correct wake geometry

The problem of steady incompressible flow for lifting surfaces is considered. An integral equation is solved relating the values of the potential discontinuity on the lifting surface and its wake to the values of the normal derivative of the potential which are known from the boundary conditions. The lifting surface and the wake are divided into small quadrilateral surface elements. The values of the potential discontinuity and the normal derivative of the potential are assumed to be constant within each lifting surface element and equal to their values at the centroids of the lifting surface elements. This yields a set of linear algebraic equations. An iteration procedure is used to obtain the wake geometry: the velocities at the corner points of the wake elements are calculated and the wake streamlines are aligned to be parallel to the velocity vector. The procedure is repeated until convergence is attained.